Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_00091332,6530TC47.4% 41.7 / 13.9 38S613P (TCC→CCC) carBcarbamoyl‑phosphate synthase large subunit
Reads supporting (aligned to +/- strand):  ref base T (17/3);  new base C (0/18);  total (17/21)
Fisher's exact test for biased strand distribution p-value = 4.02e-08
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

GCTTGGCGGCGGCCCGAACCGTATCGGTCAGGGTATCGAATTCGACTACTGTTGCGTACACGCCTCGCTGGCGCTGCGCGAAGACGGTTACGAAACCATTATGGTTAACTGTAACCC‑GGAAA‑CCGTCTCCACCGACTACGACACTTCCGACCGCCTCTACTTCGAGCCGGTAACTCTGGAAGATGTGCTGGAAATCGTGCGTATCGAGAAGCCGAAAGGCGTTATCGTCCAGTACGGCGGTCAGACCCCGCTGAAACTGGCGCGCGCGCTGGAAGCTGCTGGCGTACCGGTTATC  >  NC_000913/32508‑32802
                                                                                                                                                   |                                                                                                                                                     
gCTTGGCGGCGGCCCGAACCGTATCGGTCAGGGTATCGAATTCGACTACTGTTGCGTACACGCCTCGCTGGCGCTGCGCGAAGACGGTTACGAAACCATTATGGTTAACTGTAACCC‑GGAAA‑CCGTCTCCACCGACTACGACACTTCCGa                                                                                                                                                   <  2:607104/150‑1 (MQ=255)
           gcCCGAACCGTATCGGTCAGGGTATCGAATTCGACTACTGTTGCGTACACGCCTCGCTGGCGCTGCGCGAAGACGGTTACGAAACCATTATGGTTAACTGTAACCC‑GGAAA‑CCGTCTCCACCGACTACGACACTTCCGACCGCCTCTACtt                                                                                                                                       >  1:440248/1‑151 (MQ=255)
               gAACCGTATCGGTCAGGGTATCGAATTCGACTACTGTTGCGTACACGCCTCGCTGGCGCTGCGCGAAGACGGTTACGAAACCATTATGGTTAACTGTAACCC‑GGAAA‑CCGTCTCCACCGACTACGACACTTCCGACCGCCTCTACTTCg                                                                                                                                     >  2:322267/1‑149 (MQ=255)
                        cGGTCAGGGTATCGAATTCGACTACTGTTGCGTACACGCCTCGCTGGCGCTGCGCGAAGACGGTTACGAAACCATTATGGTTAACTGTAACCC‑GGAAA‑CCGTCTCCACCGACTACGACACTTCCGACCGCCTCTACTTCGAGCCGGTAAct                                                                                                                          >  1:383649/1‑151 (MQ=255)
                         gttCAGGTTATAAAATTCGACTACTGTGGCGTACACGCCTCGCTGGCGCTGCGCGAAGACGGTTACAAACCCATTATGGTTAACTGTACCCC‑GGAAC‑CCGTCCCCCCCGACTACGACACTTCCGCCCGCCTCTACTTCGAGCCGGTAAc                                                                                                                           <  2:579051/147‑1 (MQ=255)
                         ggccAGGGTATCGAATTCGACTCCTTTTGCGTACCCGCCTCGCTGGCCCTGCGCGAAGCCGTTTACAAAACCTTTAGGTTTACTTGTACCCC‑GAAAC‑CCGTCTCCCCCGACTCCGACCCTCCGGACCGCCTCTACTTCGAGCCGGTAAct                                                                                                                          <  2:214036/147‑1 (MQ=255)
                          gTCAGGGTATCGAATTCGACTACTGTTGCGTACACGCCTCGCTGGCGCTGCGCGAAGACGGTTACGAAACCATTATGGTTAACTGTAACCC‑GGAAA‑CCGTCTCCACCGACTACGACACTTCCGACCGCCTCTACTTCGAGCCGGTAActct                                                                                                                        >  1:69653/1‑151 (MQ=255)
                              gggTATGAAATTCGACTACTTTTGCGTCCCCGCCTCGCTGGCGCTGCGCGAAGACGGTTACAAAACCATTATGTTAACCTGTACCCC‑GAAAC‑CCGTCCCCCCCGCCTACAACACTCCCGACCGCCTCTACTTCGAGCCGGTAActct                                                                                                                        <  1:67673/147‑1 (MQ=255)
                              gggTATCGATTTCGCTTTCTTTTGGGTTCCCGCCTGGCCGGCGCCGGCCGAGGGGGGTTCCCAACCCTTATTGGTAAACTGACACCA‑GAAAC‑CCGCCCCCACCGACTACCGCATTCCCGCCCGCCTCTACTTCGAGCCGGTAACTCTGGa                                                                                                                     <  2:315623/150‑1 (MQ=255)
                                      atttCGACTACTGTTGCTTCCCCGCCTCGCTGGCGCTGCGCAAAGACGTTTACGAACCCATTATGTTTAACTGTACCCC‑GAAAC‑CCTTCTCCACCGCCTACGACACTCCCGACCGCCTCTACTTCGAGCCGGTAACTCTGGAAGAtgtg                                                                                                              <  2:477882/147‑1 (MQ=255)
                                       aTTCGACTACTGTTGCGTACACGCCTCGCTGGCGCTGCGCGAAGACGGTTACGAAACCATTATGGTTAACTGTAACCC‑GGAAA‑CCGTCTCCACCGACTACGACACTTCCGACCGCCTCTACTTCGAGCCGGTAACTCTGGAAGATGTGCTg                                                                                                           >  1:320137/1‑151 (MQ=255)
                                          cGACTACTGTTCCGTACCCGCCTCGCTGGCGCTGCGCGAGGCCGTTTACAAACCCTTTATGTTTAACTGTACCCC‑GGAAC‑CGGTCTCCCCCGACTACGACCCTTCCGACCGCCTCTACTTCGAGCCGGTAACTCTGGAAGATGTGCTGGa                                                                                                         <  2:381006/150‑1 (MQ=255)
                                                 tGTTGCGTACACGCCTCGCTGGCGCTGCGCGAAGACGGTTACGAAACCATTATGGTTAACTGTAACCC‑GGAAA‑CCGTCTCCACCGACTACGACACTTCCGACCGCCTCTACTTCGAGCCGGTAACTCTGGAAGATGTGCTGGAAATCGTGc                                                                                                 >  2:160225/1‑151 (MQ=255)
                                                     gCGTACACGCCGCGCGGGCGCTGCGCGAAGGCGGTTACGAACCCATTATGGTAAACTGTAACCC‑GAAAC‑CCGTCTCCCCCGACTAAGACACTCCCGACCGCCTCTACTTCGAGCCGGTAACTCTGGAAGATGTGCTGGAAATCGTGCGt                                                                                               <  1:167964/149‑1 (MQ=255)
                                                       gTACCCCCCCCGCGGGCGCTGCGCAAAGACGTTTACAAACCCTTTATGTTTAACTGTACCCC‑GAAAC‑CCGTCCCCCCCGACTACGACACTCCCGACCGCCTCTACTTCGAGCCGGTAACTCTGGAAGATGTGCTGGAAATCGTGCGTATCg                                                                                           <  2:320170/151‑1 (MQ=255)
                                                          cccGCCTCGCGGGCGCTGCGCGAAGCGGTTTACAAACCCATTATGTTTAACTGT‑ACCC‑GGAAACCCGTCCCCCCCGACTACGACCCTCCGGACCGCCTCTACTTCGAGCCGGTAACTCTGGAAGATGTGCTGGAAATCGTGCGTAt                                                                                             <  1:542587/144‑1 (MQ=255)
                                                          cccGCCCCGGTGGCGCTGCGCGAAGGGGGTTACAAAACCATTTTGTTTAACTGT‑ACCC‑GGAAACCCGTCTCCCCCGACTACGACCCTCCCGACCGCCTCTACTTCGAGCCGGTAACTCTGGAAGATGTGCTGGAAATCGTGCGTATCg                                                                                           <  1:372515/146‑1 (MQ=255)
                                                           acGCCTCGCTGGCGCTGCGCGAAGACGGTTACGAAACCATTATGGTTAACTGTAACCC‑GGAAA‑CCGTCTCCACCGACTACGACACTTCCGACCGCCTCTACTTCGAGCCGGTAACTCTGGAAGATGTGCTGGAAATCGTGCGTATCgaga                                                                                        >  1:365767/1‑150 (MQ=255)
                                                             gCCTCGCTGGCGCTGCGCGAAGACGGTTACGAAACCATTATGGTTAACTGTAACCC‑GGAAA‑CCGTCTCCACCGACTACGACACTTCCGACCGCCTCTACTTCGAGCCGGTAACTCTGGAAGATGTGCTGGAAATCGTGCGTATCGAGAAGc                                                                                     >  2:409825/1‑151 (MQ=255)
                                                                     ggCGCGGCGGGAAGA‑GGTTGAAAAACCCTTTGGGTTAAATGAAACCC‑GGAAACCCGTCCCCCCCGATAACGACACTCCCAACCGCCTCTACTTCGAGCCGGTAACTCTGGAAGATGTGCTGGAAATCGTGCGTATCGAGAAGCCGaaa                                                                                <  1:52986/148‑1 (MQ=255)
                                                                              cgAAGACGGTTACGAAACCATTATGGTTAACTGTAACCC‑GGAAA‑CCGTCTCCACCGACTACGACACTTCCGACCGCCTCTACTTCGAGCCGGTAACTCTGGAAGATGTGCTGGAAATCGTGCGTATCGAGAAGCCGAAAGGCGTTATCGTc                                                                    >  2:192826/1‑151 (MQ=255)
                                                                                aaGACGGTTACGAAACCATTATGTTTACCTGTACCCC‑GAAAA‑CCGTCCCCCCCGCCTACGACACTCCCGACCGCCTCTACTTCGAGCCGGTAACTCTGGAAGATGTGCTGGAAATCGTGCGTATCGAGAAGCCGAAAGGCGTTATCg                                                                      <  2:157079/147‑1 (MQ=255)
                                                                                       ttACAAAACCTTTATGGTTAACTGTACCCC‑GAAAC‑CCTTCCCCCCCGACTACGACCCTCCGGACCGCCTCTACTTCGAGCCGGTAACTCTGGAAGATGTGCTGGAAATCGTGCGTATCGAGAAGCCGAAAGGCGTTATc                                                                       <  1:173341/139‑1 (MQ=255)
                                                                                        tACGAAACCATTATGGTTAACTGTAACCC‑GGAAA‑CCGTCTCCACCGACTACGACACTTCCGACCGCCTCTACTTCGAGCCGGTAACTCTGGAAGATGTGCTGGAAATCGTGCGTATCGAGAAGCCGAAAGGCGTTATCGTCCAGTAcggcg                                                          >  1:319192/1‑151 (MQ=255)
                                                                                            aacccATTAGGTTAAACTGT‑ACCC‑GGAAACCCGTCCCCACCGACCACGACACTCCCGACCGCCTCTACTTCGAGCCGGTAACTCTGGAAGATGTGCTGGAAATCGTGCGTATCGAGAAGCCGAAAGGCGTTATCGTCCAGTACGGCGGTCa                                                      <  1:670761/148‑1 (MQ=255)
                                                                                            aaaCCATTATGTTTACCTGTACCCC‑GGAAC‑CCGTCCCCCCCGACTACGACCCTCCCGACCGCCTCTACTTCGAGCCGGTAACTCTGGAAGATGTGCTGGAAATCGTGCGTATCGAGAAGCCGAAAGGCGTTATCGTCCAGTAcggc                                                           <  1:192805/146‑1 (MQ=255)
                                                                                            aaaCCATTATGGTTAACTGTAACCC‑GGAAA‑CCGTCTCCACCGACTACGACACTTCCGACCGCCTCTACTTCGAGCCGGTAACTCTGGAAGATGTGCTGGAAATCGTGCGTATCGAGAAGCCGAAAGGCGTTATCGTCCAGTACGGCGGTc                                                       >  1:617495/1‑150 (MQ=255)
                                                                                              aCCATTATGGTTAACTGTAACCC‑GTAAA‑CCGTCTCCACCGACTACGACACTTCCGACCGCCTCTACTTCGAGCCGGTAACTCTGGAAGATGTGCTGGAAATCGTGCGTATCGAGAAGCCGAAAGGCGTTATCGTCCAGTAcggcgg                                                         >  1:192408/1‑146 (MQ=255)
                                                                                                 aTTATGGTTAACTGTAACCC‑GGAAA‑CCGTCTCCACCGACTACGACACTTCCGACCGCCTCTACTTCGAGCCGGTAACTCTGGAAGATGTGCTGGAAATCGTGCGTATCGAGAAGCCGAAAGGCGTTATCGTCCAGTACGGCGGTCAGAccc                                                 >  1:653560/1‑151 (MQ=255)
                                                                                                     tGGTTAACTGTAACCC‑GGAAA‑CCGTCTCCACCGACTACGACACTTCCGACCGCCTCTACTTCGAGCCGGTAACTCTGGAAGATGTGCTGGAAATCGTGCGTATCGAGAAGCCGAAAGGCGTTATCGTCCAGTACGGCGGTCAGACCCCGCt                                             >  2:532027/1‑151 (MQ=255)
                                                                                                       gTTAACTGTAACCC‑GGAAA‑CCGTCTCCACCGACTACGACACTTCCGACCGCCTCTACTTCGAGCCGGTAACTCTGGAAGATGTGCTGGAAATCGTGCGTATCGAGAAGCCGAAAGGCGTTATCGTCCAGTACGGCGGTCAGACCCCGCt                                             >  1:93852/1‑149 (MQ=255)
                                                                                                             ctgaaCCCGGGAAC‑CCGCCCCCCCCGACAACGCCCCTCCCGACCGCCTCTACTTCGAGCCGGTAACTCTGGAAGATGTGCTGGAAATCGTGCGTATCGAGAAGCCGAAAGGCGTTATCGTCCAGTACGGCGGTCAGACCCCGCTGaa                                          <  1:6394/144‑1 (MQ=255)
                                                                                                                            cgttcccccccGCCTACGCCCCTCCCGCCCGCCTCTACTTCGAGCCGGTAACTCTGGAAGATGTGCTGGAAATCGTGCGTATCGAGAAGCCGAAAGGCGTTATCGTCCAGTACGGCGGTCAGACCCCGCTGAAACTGGCGCGCGCGCt                           <  1:271360/139‑1 (MQ=255)
                                                                                                                                cccccccGACTACGACCCTCCCGACCGCCTCTACTTCGAGCCGGTAACTCTGGAAGATGTGCTGGAAATCGTGCGTATCGAGAAGCCGAAAGGCGTTATCGTCCAGTACGGCGGTCAGACCCCGCTGAAACTGGCGCGCGCGCTGGAAgct                    <  2:447304/146‑1 (MQ=255)
                                                                                                                                  cccccGACTACGACCCTCCCGACCGCCTCTACTTCGAGCCGGTAACTCTGGAAGATGTGCTGGAAATCGTGCGTATCGAGAAGCCGAAAGGCGTTATCGTCCAGTACGGCGGTCAGACCCCGCTGAAACTGGCGCGCGCGCTGGAAg                      <  1:221002/144‑1 (MQ=255)
                                                                                                                                    aCCGACTACGACACTTCCGACCGCCTCTACTTCGAGCCGGTAACTCTGGAAGATGTGCTGGAAATCGTGCGTATCGAGAAGCCGAAAGGCGTTATCGTCCAGTACGGCGGTCAGACCCCGCTGAAACTGGCGCGCGCGCTGGAAgctgct                 >  2:629563/1‑150 (MQ=255)
                                                                                                                                     cggACAACGCCCCTCCCGACCGCCTCTACTTCGAGCCGGTAACTCTGGAAGATGTGCTGGAAATCGTGCGTATCGAGAAGCCGAAAGGCGTTATCGTCCAGTACGGCGGTCAGACCCCGCTGAAACTGGCGCGCGCGCTGGAAGctgctg                <  1:578397/148‑1 (MQ=255)
                                                                                                                                      cgacaccgccctttCCGCCCGCCTCTTCTTCGAGCCGGTAACTCTGGAAGATGTGCTGGAAATCGTGCGTATCGAGAAGCCGAAAGGCGTTATCGTCCAGTACGGCGGTCAGACCCCGCTGAAACTGGCGCGCGCGCTGGAAGCTGCTGGc              <  2:520439/139‑1 (MQ=255)
                                                                                                                                                acTTCCGACCGCCTCTACTTCGAGCCGGTAACTCTGGAAGATGTGCTGGAAATCGTGCGTATCGAGAAGCCGAAAGGCGTTATCGTCCAGTACGGCGGTCAGACCCCGCTGAAACTGGCGCGCGCGCTGGGAGCTGCTGGCGTACCGGt      >  1:117989/1‑149 (MQ=255)
                                                                                                                                                  ttCCGACCGCCTCTACTTCGAGCCGGTAACTCTGGAAGATGTGCTGGAAATCGTGCGTATCGAGAAGCCGAAAGGCGTTATCGTCCAGTACGGCGGTCAGACCCCGCTGAAACTGGCGCGCGCGCTGGAAGCTGCTGGCGTACCGGTTATc  >  2:568454/1‑151 (MQ=255)
                                                                                                                                                   |                                                                                                                                                     
GCTTGGCGGCGGCCCGAACCGTATCGGTCAGGGTATCGAATTCGACTACTGTTGCGTACACGCCTCGCTGGCGCTGCGCGAAGACGGTTACGAAACCATTATGGTTAACTGTAACCC‑GGAAA‑CCGTCTCCACCGACTACGACACTTCCGACCGCCTCTACTTCGAGCCGGTAACTCTGGAAGATGTGCTGGAAATCGTGCGTATCGAGAAGCCGAAAGGCGTTATCGTCCAGTACGGCGGTCAGACCCCGCTGAAACTGGCGCGCGCGCTGGAAGCTGCTGGCGTACCGGTTATC  >  NC_000913/32508‑32802

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