Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009134,133,2320AG42.9% 50.3 / 15.5 42E206G (GAA→GGA) metF5,10‑methylenetetrahydrofolate reductase
Reads supporting (aligned to +/- strand):  ref base A (4/20);  new base G (18/0);  total (22/20)
Fisher's exact test for biased strand distribution p-value = 2.12e-08
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

CGGAAGCAAAAAGCGCTCAGGCGGATTTGCTTAATCTGAAACGCAAAGTGGATGCCGGAGCCAACCGCGCGATTACTCAGTTCTTCTTCGATGTCGAAAGCTACCTGCGTTTTCGTGACCGCTGTGTATCGGCGGGCATTG‑ATGTGGAAATTA‑TTCCGGGAATTTTGCCGGT‑ATCTAACTTTAAACAGGCGAAGAAATTTGCCGATATGACCAACGTGCGTATTCCGGCGTGGATGGCGCAAATGTTCGACGGTCTGGATGATGATGCCGAAACCCGCAAACTGGTTGGCGCGAA  >  NC_000913/4133085‑4133379
                                                                                                                                                    |                                                                                                                                                     
cGGAAGCAAAAAGCGCTCAGGCGGATTTGCTTAATCTGAAACGCAAAGTGGATGCCGGAGCCAACCGCGCGATTACTCAGTTCTTCTTCGATGTCGAAAGCTACCTGCGTTTTCGTGACCGCTGTGTATCGGCGGGCATTG‑ATGTGGgaa                                                                                                                                                     >  1:553060/1‑150 (MQ=255)
  gAAGCAAAAAGCGCTCAGGCGGATTTGCTTAATCTGAAACGCAAAGTGGATGCCGGAGCCAACCGCGCGATTACTCAGTTCTTCTTCGATGTCGAAAGCTACCTGCGTTTTCGTGACCGCTGTGTATCGGCGGGCATTG‑ATGTGGAAatta                                                                                                                                                  <  1:698770/151‑1 (MQ=255)
       aaaaaGCGCTCAGGCGGATTTGCTTAATCTGAAACGCAAAGTGGATGCCGGAGCCAACCGCGCGATTACTCAGTTCTTCTTCGATGTCGAAAGCTACCTGCGTTTTCGTGACCGCTGTGTATCGGCGGGCCGTG‑GTTTGGGAaata‑tt                                                                                                                                               >  1:387625/1‑148 (MQ=255)
        aaaaGCGCTCAGGCGGATTTGCTTAATCTGAAACGCAAAGTGGATGCCGGAGCCAACCGCGCGATTACTCAGTTCTTCTTCGATGTCGAAAGCTACCTGCGTTTTCGTGACCGCTGTGTATCGGCGGTCCGTG‑ATGTGGGAatta‑tt                                                                                                                                               >  1:686607/1‑147 (MQ=255)
         aaaGCGCTCAGGCGGATTTGTTTAATCTGAAACGCAAAGGGGATGCCGGAGCCAACCGCGCGATTACTCAGTTCTTCTTCGATGTCCAAAGCTACCTGCGTTTTCGTGACCGCTGTGTATCGGCGGGCATTG‑ATGTGGAAATTA‑TTCCggg                                                                                                                                          <  1:609816/151‑1 (MQ=255)
            gcgcTCAGGCGGATTTGCTTAATCTGAAACGCAAAGTGGATGCCGGAGCCAACCGCGCGATTACTCAGTTCTTCTTCGATGTCGAAAGCTACCTGCGTTTTCGTGACCGCTGTGTATCGGCGGGCATTT‑GTGTGGGAATTATTTCCggg                                                                                                                                          >  2:768098/1‑149 (MQ=255)
               cTCAGGCGGATTTGCTTAATCTGAAACGCAAAGTGGATGCCGGAGCCAACCGCGCGATTACTCAGTTCTTCTTCGATGTCGAAAGCTACCTGCGTTTTCGTGACCGCTGTGTATCGGCGGGGATTT‑ATGTTGGAATTT‑TTTCTGGAAtttt                                                                                                                                    >  2:43258/1‑151 (MQ=255)
                    gCGGATTTGCTTAATCTGAAACGCAAAGTGGATGCCGGAGCCAACCGCGCGATTACTCAGTTCTTCTTCGATGTCGAAAGCTACCTGCGTTTTCGTGACCGCTGTGTATCGGCGGGCATTG‑ATGTGGAAATTA‑TTCCGGGAATTTTGCCgg                                                                                                                               <  2:546319/151‑1 (MQ=255)
                                    tGAAACGCAAAGTGGATGCCGGAGCCAACCGCGCGATTACTCAGTTCTTCTTCGATGTCGAAAGCTACCTGCGTTTTCGTGACCGCTGTGTATCGGCGGGCCTTT‑ATGTTGGAAATA‑TTTCGGGGATTTTTTCGGG‑ATCTAACTTTtaa                                                                                                                >  2:27449/1‑149 (MQ=255)
                                     gAAACGCAAAGTGGATGCCGGAGCCAACCGCGCGATTACTCAGTTCTTCTTCGATGTCGAAAGCTACCTGCGTTTTCGTGACCGCTGTGTATCGGCGGGCATTG‑ATGTGGAAATTA‑TTCCGGGAATTTTGCCGGT‑ATCTAACTTTAAAc                                                                                                               <  2:144587/149‑1 (MQ=255)
                                     gAAACGCAAAGTGGATGCCGGAGCCAACCGCGCGATTACTCAGTTCTTCTTCGATGTCGAAAGCTACCTGCGTTTTCGTGACCGCTGTGTATCGGCGGGCATTG‑ATGTGGAAATTA‑TTCCGGGAATTTTGCCGGT‑ATCTAACTTTAAACa                                                                                                              <  1:389560/150‑1 (MQ=255)
                                      aaaCGCAAAGTGGATGCCGGAGCCAACCGCGCGATTACTCAGTTCTTCTTCGATGTCGAAAGCTACCTGCGTTTTCGTGACCGCTGTGTATCGGCGGGCATTG‑ATGTGGAAATTA‑TTCCGGGAATTTTGCCGGT‑ATCTAACTTTAAACAgg                                                                                                            <  2:99325/151‑1 (MQ=255)
                                           cAAAGTGGATGCCGGAGCCAACCGCGCGATTACTCAGTTCTTCTTCGATGTCGAAAGCTACCTGCGTTTTCGTGACCGCTGTGTATCGGCGGGCATTG‑ATTTGGGAATTATTTTCGGGAGTTTT‑CCGGG‑ATCTTACTTTtaa                                                                                                                >  2:308598/1‑142 (MQ=255)
                                                 ggATGCCGGAGCCAACCGCGCGATTACTCAGTTCTTCTTCGATGTCGAAAGCTACCTGCGTTTTCGTGACCGCTGTGTATCTGCTGGCATTG‑ATGTGGAAATTA‑TTCCGGGAATTTTGCCGGT‑ATCTAACTTTAAACAGGCGAAGAAAtt                                                                                                  <  2:461391/150‑1 (MQ=255)
                                                 ggATGCCGGAGCCAACCGCGCGATTACTCAGTTCTTCTTCGATGTCGAAAGCTACCTGCGTTTTCGTGACCGCTGTGTATCGGCGGGCATTG‑ATGTGGAAATTA‑TTCCGGGAATTTTGCCGGT‑ATCTAACTTTAAACAGGCGAAGAAAttt                                                                                                 <  1:137954/151‑1 (MQ=255)
                                                              aaCCGCGCGATTACTCAGTTCTTCTTCGATGTCGAAAGCTACCTGCGTTTTCGTGACCGCTGTGTATCGGCGGGCCTTT‑AAGTGGGAATTT‑TTTCCGGAAGTTTGCCCGC‑GTCCTAATTTATAAAGGCGGAGAAAATTGGCggaatta                                                                                       >  1:50897/1‑142 (MQ=255)
                                                              aaCCGCGCGATTACTCAGTTCTTCTTCGATGTCGAAAGCTACCTGCGTTTTCGTGACCGCTGTGTATCGGCGGGCATAGTAGTTGGGATTTT‑TTCTGGGAATTTT‑CCGCG‑GTCTTACCTTTAACAAGCGGATAAATTTgacca                                                                                            >  2:217619/1‑139 (MQ=255)
                                                                 cgcgcgATTACTCAGTTCTTCTTCGATGTCGAAAGCTACCTGCGTTTTCGTGACCGCTGTGTATCGGCGGGGATTG‑ATGTGGGAATTA‑TTTCGGGGATTTTTCCGGG‑ATCTTACTTTTAAAAGGCGAAGAAATTTGCCGATATGACCaa                                                                                   >  1:116941/1‑149 (MQ=255)
                                                                  gcgcgATTACTCAGTTCTTCTTCGATGTCGAAAGCTACCTGCGTTTTCGTGACCGCTGTGTATCGGCGGGCGTTT‑GTGTGGGAAATT‑TTTCGGGGATTTTTTCGGG‑ATCTAACTTTAAAAAGGGGAAGAAATTTGGCGATATggacaa                                                                                   >  2:388172/1‑143 (MQ=255)
                                                                   cgcgATTACTCAGTTCTTCTTCGATGTCGAAAGCTACCTGCGTTTTCGTGACCGCTGTGTATCGGCGGGCCTTG‑GTGTGGGAATTA‑TTCCGGGAATTTTTTCGGG‑ATCTTAATTTTAACAAGGGGAGAAAATTGCCGATATGaac                                                                                     >  1:619049/1‑143 (MQ=255)
                                                                      gATTACTCAGTTCTTCTTCGATGTCGAAAGCTACCTGCGTTTTCGTGACCGCTGTGTATCGGCGGGCGTTG‑TTGTGGGAATTA‑TTTCGGGGGTTTTTTCGGC‑GTCTAAATTTTAAAAGGGGAAGAAATTTTCCGGGATGAACAACGAGCg                                                                             >  2:366818/1‑150 (MQ=255)
                                                                      gATTACTCAGTTCTTCTTCGATGTCGAAAGCTACCTGCGTTTTCGTGACCGCTGTGTATCGGCGGGCATTG‑ATGTGGAAATTA‑TTCCGGGAATTTTGCCGGT‑ATCTAACTTTAAACAGGCGAAGAAATTTGCCGATATGACCAACGTGCGt                                                                            <  2:29022/151‑1 (MQ=255)
                                                                             cAGTTCTTCTTCGATGTCGAAAGCTACCTGCGTTTTCGTGACCGCTGTGTATCGGCGGGCATTG‑ATGTGGGAATTA‑TTTCGGGGAATTTTTCGGG‑ATGTAACTTTAAAAAGGAGAAGAAATTTGCCGATATGatccacacgtgtatt                                                                         >  2:136917/1‑133 (MQ=255)
                                                                             cAGTTCTTCTTCGATGTCGAAAGCTACCTGCGTTTTCGTGACCGCTGTGTATCGGCGGGCATTG‑ATGTGGAAATTA‑TTCCGGGAATTTTGCCGGT‑ATCTAACTTTAAACAGGCGAAGAAATTTGCCGATATGACCAACGTGCGTATTcc                                                                       <  1:711140/149‑1 (MQ=255)
                                                                              aGTTCTTCTTCGATGTCGAAAGCTACCTGCGTTTTCGTGACCGCTGTGTATCGGCGGGGCTTG‑ATGTGGGAATTA‑TTTTGGGGAATTTTCCGGT‑ATCTTACTTTTAAAATGCGAAGAAATTTGCCGGTATTACCAACCTGCGtttt                                                                         >  2:725255/1‑146 (MQ=255)
                                                                                    tcttcGATGTCGAAAGCTACCTGCGTTTTCGTGACCGCTGTGTATCGGCGGGCATTG‑ATGTGGAAATTA‑TTCCGGGAATTTTGCCGGT‑ATCTAACTTTAAACAGGCGAAGAAATTTGCCGATATGACCAACGTGCGTATTCCGGCGTgg                                                                <  1:397377/149‑1 (MQ=255)
                                                                                        cGATGTCGAAAGCTACCTGCGTTTTCGTGACCGCTGTGTATCGGCGGGCATTG‑ATGTGGAAATTA‑TTCCGGGAATTTTGCCGGT‑ATCTAACTTTAAACAGGCGAAGAAATTTGCCGATATGACCAACGTGCGTATTCCGGCGTGGATGgc                                                           <  1:639829/150‑1 (MQ=255)
                                                                                          aTGTCGAAAGCTACCTGCGTTTTCGTGACCGCTGTGTATCGGCGGGCATTG‑ATGTGGAAATTA‑TTCCGGGAATTTTGCCGGT‑ATCTAACTTTAAACAGGCGAAGAAATTTGCCGATATGACCAACGTGCGTATTCCGGCGTGGATGGCGCa                                                        <  2:696324/151‑1 (MQ=255)
                                                                                               gAAAGCTACCTGCGTTTTCGTGACCGCTGTGTATCGGCGGGGGATG‑ATGAGGAGATTA‑TTCCGGGAGTTTTTGCGGG‑CTCTACCTTTAAACAGGCGGAGAAAATTGTCGCAATTACAAACGTGCGGTTTTCGGCGCgggcggagccaa                                                      >  2:503582/1‑138 (MQ=255)
                                                                                               gAAAGCTACCTGCGTTTTCGTGACCGCTGTGTATCGGCGGGCATTG‑ATGTGGGAATTA‑TTTCGGGGATTTTTCCGGG‑ATCTAACTTTTAACAGGCGGAGAAATTTGCCGATATGACCAACGTGCGTATTCCGGCGGGGGTGGCGCaaa                                                      >  1:158223/1‑148 (MQ=255)
                                                                                                   gCTACCTGCGTTTTCGTGACCGCTGTGTATCGGCGGGCATTG‑ATGTGGAAATTA‑TTCCGGGAATTTTGCCGGT‑ATCTAACTTTAAACAGGCGAAGAAATTTGCCGATATGACCAACGTGCGTATTCCGGCGTGGATGGCGCAAATGTTCGa                                               <  1:752911/151‑1 (MQ=255)
                                                                                                    cTACCTGCGTTTTCGTGACCGCTGTGTATCGGCGGGCATTG‑TTGTGGAAAATT‑TTTCGGGGATTTTTTCGGG‑ATCTAAATTTTAAAAGGGGAAGAAAATTGTCGGTATGACCAACGTGCGTATTTTGGCGGGGGTGGGGCAAAAGTTTGa                                               >  2:375138/1‑150 (MQ=255)
                                                                                                      aCCTGCGTTTTCGTGACCGCTGTGTATCGGCGGGCATTG‑ATGTGGAAATTA‑TTCCGGGAATTTTGCCGGT‑ATCTAACTTTAAACAGGCGAAGAAATTTGCCGATATGACCAACGTGCGTATTCCGGCGTGGATGGCGCAAATGTTCGACg                                             <  2:352346/150‑1 (MQ=255)
                                                                                                            gTTTTCGTGACCGCTGTGTATCGGCGGGGATTG‑ATGTGGGAAATA‑TTTCCGGGATTTTTCCGGTAATTTAATTTTAAAAAGG‑GAAGAAAATTTGCGATATTACCCAAGGGGGGATTCCCGGGGGGGTGGGGGAAAAGTTCGAAGGTCTgg                                       >  2:515956/1‑150 (MQ=255)
                                                                                                            gTTTTCGTGACCGCTGTGTATCGGCGGGCCTTG‑TTTTGGGAAATT‑TTTCGGGGATTTTGCCGGT‑ATCTAACTTTTAACAGGCGAAGAAAATTTCCGCTATGAACAAAGTGCGGATTTCGGCGTGGGTGGGGCAAATGTTCGGCGGGCTgg                                       >  1:96965/1‑150 (MQ=255)
                                                                                                                tCGTGACCGCTGTGTATCGGCGGGCAATG‑TTGTGGGAATTA‑TTTCGGGGATTTTGCCGGG‑ATCTAACTTTAAACAGGCGAAGAAATTTGCCGGTATGACCAAAGTGCGTATTTCGGCGGGGGTGGGGCAAATGTTCGACGGTCTg                                        >  1:582668/1‑145 (MQ=255)
                                                                                                                       cGCTGTGTATCGGCGGGCATTG‑ATGTGGAAATTA‑TTCCGAGAATTTTGCCGGT‑ATCTAACTTTAAACAGGCGAAGAAATTTGCCGATATGACCAACGTGCGTATTCCGGCGTGGATGGCGCAAATGTTCGACGGTCTGGATGATGATGCCg                           <  2:789896/151‑1 (MQ=255)
                                                                                                                           gtgtATCGGCGGGCATTG‑ATGTGGAAATTA‑TTCCGGGAATTTTGCCGGT‑ATCTAACTTTAAACAGGCGAAGAAATTTGCCGATATGACCAACGTGCGTATTCCGGCGTGGATGGCGCAAATGTTCGACGGTCTGGATGATGATGCCGaa                         <  2:642590/149‑1 (MQ=255)
                                                                                                                                    cggGCATTG‑ATGTGGAAATTA‑TTCCGGGAATTTTGCCGGT‑ATCTAACTTTAAACAGGCGAAGAAATTTGCCGATATGACCAACGTGCGTATTCCGGCGTGGATGGCGCAAATGTTCGACGGTCTGGATGATGATGCCGAAACCCGCAAACt              <  2:7946/151‑1 (MQ=255)
                                                                                                                                      ggCATTG‑ATGTGGAAATTA‑TTCCGGGAATTTTGCCGGT‑ATCTAACTTTAAACAGGCGAAGAAATTTGCCGATATGACCAACGTGCGTATTCCGGCGTGGATGGCGCAAATGTTCGACGGTCTGGATGATGATGCCGAAACCCGCAAACTgg            >  1:378316/1‑151 (MQ=255)
                                                                                                                                              aTGTGGAAATTA‑TTCCGGGAATTTTGCCGGT‑ATCTAACTTTAAACAGGCGAAGAAATTTGCCGATATGACCAACGTGCGTATTCCGGCGTGGATGGCGCAAATGTTCGACGGTCTGGATGATGATGCCGAAACCCGCAAACTGGTTGgcgc     <  2:30869/151‑1 (MQ=255)
                                                                                                                                                 tgGAAATTA‑TTCCGGGAATTTTGCCGGT‑ATCTAACTTTAAACAGGCGAAGAAATTTGCCGATATGACCAACGTGCGTATTCCGGCGTGGATGGCGCAAATGTTCGACGGTCTGGATGATGATGCCGAAACCCGCAAACTGGTTGGCGCGa   >  2:586851/1‑150 (MQ=255)
                                                                                                                                                  ggAAATTA‑TTCCGGGAATTTTGCCGGT‑ATCTAACTTTAAACAGGCGAAGAAATTTGCCGATATGACCAACGTGCGTATTCCGGCGTGGATGGCGCAAATGGTCGACGGTCTGGATGATGATGCCGAAACCCGCAAACTGGTTGGCGCGaa  <  2:463409/150‑1 (MQ=255)
                                                                                                                                                    |                                                                                                                                                     
CGGAAGCAAAAAGCGCTCAGGCGGATTTGCTTAATCTGAAACGCAAAGTGGATGCCGGAGCCAACCGCGCGATTACTCAGTTCTTCTTCGATGTCGAAAGCTACCTGCGTTTTCGTGACCGCTGTGTATCGGCGGGCATTG‑ATGTGGAAATTA‑TTCCGGGAATTTTGCCGGT‑ATCTAACTTTAAACAGGCGAAGAAATTTGCCGATATGACCAACGTGCGTATTCCGGCGTGGATGGCGCAAATGTTCGACGGTCTGGATGATGATGCCGAAACCCGCAAACTGGTTGGCGCGAA  >  NC_000913/4133085‑4133379

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