Predicted mutation
evidence seq id position mutation freq annotation gene description
RA W3110S.gb 188,973 A→G 9.5% H178H (CAT→CAC map ← methionine aminopeptidase

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*W3110S.gb188,9730AG 9.5% 105.6 / 2.8 42H178H (CAT→CACmapmethionine aminopeptidase
Reads supporting (aligned to +/- strand):  ref base A (22/16);  new base G (2/2);  total (24/18)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

GAACGTCATCCCAGGTTTCAGTACGACGTTGGTTTCACGGGAGTCATAGTGCAGCACCTGCGGTTCTTCATGGAAGCCGCGACCAATACCGTGTCCGCAATATTCACGAACGACGGAGAAGCCTTCTGCTTCGACAAATT  >  W3110S.gb/188904‑189043
                                                                     |                                                                      
gAACGTCATCCCAGGTTTCAGTACGACGTTGATTTCACGGGAGTCATAGTGCAGCACCTGCGGTTCTTCAt                                                                       >  1:410522/1‑71 (MQ=255)
 aaCGTCATCCCAGGTTTCAGTACGACGTTGGTTTCACGGGAGTCATAGTGCAGCACCTGCGGTTCTTCAt                                                                       <  1:1278238/70‑1 (MQ=255)
   cGTCATCCCAGGTTTCAGTACGACGTTGGTTTCACGGGAGTCATAGTGCAGCACCTGCGGTTCTTCATgg                                                                     >  1:2498770/1‑70 (MQ=255)
       aTCCCAGGTTTCAGTACGACGTTGGTTTCACGGGAGTCATAGTGCAGCACCTGCGGTTCTTCATGGAAGcc                                                                >  1:2677885/1‑71 (MQ=255)
        tCCCAGGTTTCAGTACGACGTTGGTTTCACGGGAGTCATAGTGCAGCACCTGCGGTTCTTCATGGAAGCcg                                                               >  1:4377711/1‑71 (MQ=255)
          ccAGGTTTCAGTACGACGTTGGTTTCACGGGAGTCATAGTGCAGCACCTGCGGTTCTTCATGGAAGCcgc                                                              >  1:4304527/1‑70 (MQ=255)
          ccAGGTTTCAGTACGACGTTGGTTTCACGGGAGTCATAGTGCAGCACCTGCGGTTCTTCATGGAAGCcgc                                                              >  1:2862042/1‑70 (MQ=255)
            aGGTTTCAGTACGACGTTGGTTTCACGGGAGTCATAGTGCAGCACCTGCGGTTCTTCATGGAAGCCGCGAc                                                           <  1:2723697/71‑1 (MQ=255)
             ggTTTCAGTACGACGTTGGTTTCACGGGAGTCATAGTGCAGCACCTGCGGTTCTTCATGGAAGCCGCGAc                                                           >  1:3942518/1‑70 (MQ=255)
              gTTTCAGTACGACGTTGGTTTCACGGGAGTCATAGTGCAGCACCTGCGGTTCTTCATGGAAGCCGCGACCa                                                         <  1:670026/71‑1 (MQ=255)
                 tCAGTACGACGTTGGTTTCACGGGAGTCATAGTGCAGCACCTGCGGTTCTTCa                                                                        >  1:1612814/1‑53 (MQ=255)
                 tCAGTACGACGTTGGTTTCACGGGAGTCATAGTGCAGCACCTGCGGTTCTTCATGGAAGCCGCGACCAAt                                                       <  1:1105876/70‑1 (MQ=255)
                   agcacgacgTTGGTTTCACGGGAGTCATAGTGCAGCACCTGCGGTTCTTCGTGGAAGCCGCGACCAATAc                                                     <  1:4151466/67‑1 (MQ=255)
                   agcacgacgTTGGTTTCACGGGAGTCATAGTGCAGCACCTGCGGTTCTTCGTGGAAGCCGCGACCAATAc                                                     <  1:96383/67‑1 (MQ=255)
                   aGTACGACGTTGGTTTCACGGGAGTCATAGTGCAGCACCTGCGGTTCTTCATTGAAGCCGCGACCAATAc                                                     <  1:2021254/70‑1 (MQ=255)
                   aGTACGACGTTGGTTTCACGGGAGTCATAGTGCAGCACCTGCGGTTCTTCATGGAAGCCGCGACCAATAcc                                                    <  1:3169305/71‑1 (MQ=255)
                   aGTACGACGTTGGTTTCACGGGAGTCATAGTGCAGCACCTGCGGTTCTTCATGGAAGCCGCGACCAATAcc                                                    <  1:883407/71‑1 (MQ=255)
                    gTACGACGTTGGTTTCACGGGAGTCATAGTGCAGCACCTGCGGTTCTTCATGGAAGCCGCGACCAATAcc                                                    <  1:4052533/70‑1 (MQ=255)
                      acgacgTTGGTTTCACGGGAGTCATAGTGCAGCACCTGCGGTTCTTCATGGAAGCCGCGACCAATACCgtg                                                 >  1:2509009/1‑71 (MQ=255)
                      acgacgTTGGTTTCACGGGAGTCATAGTGCAGCACCTGCGGTTCTTCATGGAAGCCGCGACCAATACCgtg                                                 >  1:1528832/1‑71 (MQ=255)
                       cgacgTTGGTTTCACGGGAGTCATAGTGCAGCACCTGCGGTTCTTCATGGa                                                                    >  1:3289502/1‑51 (MQ=255)
                             tGGTTTCACGGGAGTCATAGTGCAGCACCTGCGGTTCTTCATGGAAGCCGCGACCaa                                                        >  1:3853497/1‑57 (MQ=255)
                             tGGTTTCACGGGAGTCATAGTGCAGCACCTGCGGTTCTTCATGGAAGCCGCGACCaa                                                        >  1:4030788/1‑57 (MQ=255)
                             tGGTTTCACGGGAGTCATAGTGCAGCACCTGCGGTTCTTCATGGAAGCCGCGACCAATACCGTGTCCGCaa                                          >  1:492095/1‑71 (MQ=255)
                              ggTTTCACGGGAGTCATAGTGCAGCACCTGCGGTTCTTCATGGAAg                                                                  >  1:1738988/1‑46 (MQ=255)
                               gTTTCACGGGAGTCATAGTGCAGCACCTGCGGTTCTTCATGGAAGCCGCGAcc                                                          <  1:1324647/53‑1 (MQ=255)
                                  tCACGGGAGTCATAGTGCAGCACCTGCGGTTCTTCATGGAAGCCGCGACCAATACCGTGTCCGCAATAtt                                      >  1:197062/1‑70 (MQ=255)
                                        gAGTCATAGTGCAGCACCTGCGGTTCTTCATGGAAGCCGCGACCAATACCGTGTCCGCAATATTCacgaac                               <  1:1253754/71‑1 (MQ=255)
                                        gAGTCATAGTGCAGCACCTGCGGTTCTTCATGGAAGCCGCGACCAATACCGTGTCCGCAATATTCacgaac                               <  1:813085/71‑1 (MQ=255)
                                             aTAGTGCAGCACCTGCGGTTCTTCATGGAAGCCGCGACCAATACCGTGTCCTCAATATTCACGAaccacg                           <  1:3240072/70‑1 (MQ=255)
                                             aTAGTGCAGCACCTGCGGTTCTTCATGGAAGCCGCGACCAATACCGTGTCCGCAATATTCACGAacgacg                           >  1:2880832/1‑70 (MQ=255)
                                              tAGTGCAGCACCTGCGGTTCTTCATGGAAGCCGCGACCAATACCGTGTCCGCAATATTCacgaacga                             <  1:3696076/67‑1 (MQ=255)
                                                gTGCAGCACCTGCGGTTCTTCATGGAAGCCGCGACCAATAcc                                                    >  1:3448052/1‑42 (MQ=255)
                                                  gcagcaCCTGCGGTTCTTCATGGAAGCCGCGACCAATACCGTGTCCGCAat                                         >  1:433324/1‑51 (MQ=255)
                                                        ccTGCGGTTCTTCATGGAAGCCGCGACCAATACCGTGTCCGCAATATTCACGAACGACGGAGAAGCCTTCt               <  1:3752883/71‑1 (MQ=255)
                                                          tGCGGTTCTTCATGGAAGCCGCGACCAATACCGTGTCCGCAATATTCACGAACGACGgaga                       >  1:2717223/1‑61 (MQ=255)
                                                            cGGTTCTTCATGGAAGCCGCGACCAATACCGTGTCCGCAATATTCACGAACGACGGAGAAGCCTTCTGCt            <  1:2149343/70‑1 (MQ=255)
                                                             ggTTCTTCATGGAAGCCGCGACCAATACCGTGTCCGCAATATTCACGAACGACGGAGAAGCCTTCTGCTTc          <  1:1521052/71‑1 (MQ=255)
                                                              gTTCTTCGTGGAAGCGGCGACCAATACCGTGTCCACAatat                                       >  1:1338174/1‑41 (MQ=25)
                                                              gTTCTTCGTGGAAGCCGCGACCAATACCGTGTCCACAATAtt                                      >  1:2272180/1‑42 (MQ=38)
                                                                  ttcATGGAAGCCGCGACCAATACCGTGTCCGCAATATTCACGAACGACGGAGAAGc                    >  1:1182695/1‑56 (MQ=255)
                                                                   tcATGGAAGCCGCGACCAATACCGTGTCCGCAata                                        >  1:4431492/1‑35 (MQ=255)
                                                                    cATGGAAGCCGCGACCAATACCGTGTCCGCAATATTCac                                   >  1:75329/1‑39 (MQ=255)
                                                                     aTGGAAGCCGCGACCAATACCGTGTCCGCAATATTCACGAACGACGGAGAAGCCTTCTGCTTCGACAAAtt  >  1:171557/1‑71 (MQ=255)
                                                                     aTGGAAGCCGCGACCAATACCGTGTCCGCAATATTCACGAACGACGGAGAAGCCTTCTGCTTCGACAAAt   <  1:3962038/70‑1 (MQ=255)
                                                                     |                                                                      
GAACGTCATCCCAGGTTTCAGTACGACGTTGGTTTCACGGGAGTCATAGTGCAGCACCTGCGGTTCTTCATGGAAGCCGCGACCAATACCGTGTCCGCAATATTCACGAACGACGGAGAAGCCTTCTGCTTCGACAAATT  >  W3110S.gb/188904‑189043

Alignment Legend
Aligned base mismatch/match (shaded by quality score): ATCG/ATCG < 3 ≤ ATCG/ATCG < 21 ≤ ATCG/ATCG < 32 ≤ ATCG/ATCG < 36 ≤ ATCG/ATCG
Unaligned base: atcg    Masked matching base: atcg    Alignment gap:     Deleted base: