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

GCTTGGCGGCGGCCCGAACCGTATCGGTCAGGGTATCGAATTCGACTACTGTTGCGTACACGCCTCGCTGGCGCTGCGCGAAGACGGTTACGAAACCATTATGGTTAACTGTAACCCGGAAACCGTCTCCACCGACTACGACACTTCCGACCGCCTCTACTTCGAGCCGGTAACTCTGGAAGATGTGCTGGAAATCGTGCGTATCGAGAAGCCGAAAGGCGTTATCGTCCAGTACGGCGGTCAGACCCCGCTGAAACTGGCGCGCGCGCTGGAAGCTGCTGGCGTACCGGTTATC  >  NC_000913/32508‑32802
                                                                                                                                                   |                                                                                                                                                     
gCTTGGCGGCGGCCCGAACCGTATCGGTCAGGGTATCGAATTCGACTACTGTTGCGTACACGCCTCGCTGGCGCTGCGCGAAGACGGTTACGAAACCATTATGGTTAACTGTAACCCGGAAACCGTCTCCACCGACTACGACACTTCCGa                                                                                                                                                   <  2:607104/150‑1 (MQ=255)
           gcCCGAACCGTATCGGTCAGGGTATCGAATTCGACTACTGTTGCGTACACGCCTCGCTGGCGCTGCGCGAAGACGGTTACGAAACCATTATGGTTAACTGTAACCCGGAAACCGTCTCCACCGACTACGACACTTCCGACCGCCTCTACtt                                                                                                                                       >  1:440248/1‑151 (MQ=255)
               gAACCGTATCGGTCAGGGTATCGAATTCGACTACTGTTGCGTACACGCCTCGCTGGCGCTGCGCGAAGACGGTTACGAAACCATTATGGTTAACTGTAACCCGGAAACCGTCTCCACCGACTACGACACTTCCGACCGCCTCTACTTCg                                                                                                                                     >  2:322267/1‑149 (MQ=255)
                        cGGTCAGGGTATCGAATTCGACTACTGTTGCGTACACGCCTCGCTGGCGCTGCGCGAAGACGGTTACGAAACCATTATGGTTAACTGTAACCCGGAAACCGTCTCCACCGACTACGACACTTCCGACCGCCTCTACTTCGAGCCGGTAAct                                                                                                                          >  1:383649/1‑151 (MQ=255)
                         gttCAGGTTATAAAATTCGACTACTGTGGCGTACACGCCTCGCTGGCGCTGCGCGAAGACGGTTACAAACCCATTATGGTTAACTGTACCCCGGAACCCGTCCCCCCCGACTACGACACTTCCGCCCGCCTCTACTTCGAGCCGGTAAc                                                                                                                           <  2:579051/147‑1 (MQ=255)
                         ggccAGGGTATCGAATTCGACTCCTTTTGCGTACCCGCCTCGCTGGCCCTGCGCGAAGCCGTTTACAAAACCTTTAGGTTTACTTGTACCCCGAAACCCGTCTCCCCCGACTCCGACCCTCCGGACCGCCTCTACTTCGAGCCGGTAAct                                                                                                                          <  2:214036/147‑1 (MQ=255)
                          gTCAGGGTATCGAATTCGACTACTGTTGCGTACACGCCTCGCTGGCGCTGCGCGAAGACGGTTACGAAACCATTATGGTTAACTGTAACCCGGAAACCGTCTCCACCGACTACGACACTTCCGACCGCCTCTACTTCGAGCCGGTAActct                                                                                                                        >  1:69653/1‑151 (MQ=255)
                              gggTATGAAATTCGACTACTTTTGCGTCCCCGCCTCGCTGGCGCTGCGCGAAGACGGTTACAAAACCATTATGTTAACCTGTACCCCGAAACCCGTCCCCCCCGCCTACAACACTCCCGACCGCCTCTACTTCGAGCCGGTAActct                                                                                                                        <  1:67673/147‑1 (MQ=255)
                              gggTATCGATTTCGCTTTCTTTTGGGTTCCCGCCTGGCCGGCGCCGGCCGAGGGGGGTTCCCAACCCTTATTGGTAAACTGACACCAGAAACCCGCCCCCACCGACTACCGCATTCCCGCCCGCCTCTACTTCGAGCCGGTAACTCTGGa                                                                                                                     <  2:315623/150‑1 (MQ=255)
                                      atttCGACTACTGTTGCTTCCCCGCCTCGCTGGCGCTGCGCAAAGACGTTTACGAACCCATTATGTTTAACTGTACCCCGAAACCCTTCTCCACCGCCTACGACACTCCCGACCGCCTCTACTTCGAGCCGGTAACTCTGGAAGAtgtg                                                                                                              <  2:477882/147‑1 (MQ=255)
                                       aTTCGACTACTGTTGCGTACACGCCTCGCTGGCGCTGCGCGAAGACGGTTACGAAACCATTATGGTTAACTGTAACCCGGAAACCGTCTCCACCGACTACGACACTTCCGACCGCCTCTACTTCGAGCCGGTAACTCTGGAAGATGTGCTg                                                                                                           >  1:320137/1‑151 (MQ=255)
                                          cGACTACTGTTCCGTACCCGCCTCGCTGGCGCTGCGCGAGGCCGTTTACAAACCCTTTATGTTTAACTGTACCCCGGAACCGGTCTCCCCCGACTACGACCCTTCCGACCGCCTCTACTTCGAGCCGGTAACTCTGGAAGATGTGCTGGa                                                                                                         <  2:381006/150‑1 (MQ=255)
                                                 tGTTGCGTACACGCCTCGCTGGCGCTGCGCGAAGACGGTTACGAAACCATTATGGTTAACTGTAACCCGGAAACCGTCTCCACCGACTACGACACTTCCGACCGCCTCTACTTCGAGCCGGTAACTCTGGAAGATGTGCTGGAAATCGTGc                                                                                                 >  2:160225/1‑151 (MQ=255)
                                                     gCGTACACGCCGCGCGGGCGCTGCGCGAAGGCGGTTACGAACCCATTATGGTAAACTGTAACCCGAAACCCGTCTCCCCCGACTAAGACACTCCCGACCGCCTCTACTTCGAGCCGGTAACTCTGGAAGATGTGCTGGAAATCGTGCGt                                                                                               <  1:167964/149‑1 (MQ=255)
                                                       gTACCCCCCCCGCGGGCGCTGCGCAAAGACGTTTACAAACCCTTTATGTTTAACTGTACCCCGAAACCCGTCCCCCCCGACTACGACACTCCCGACCGCCTCTACTTCGAGCCGGTAACTCTGGAAGATGTGCTGGAAATCGTGCGTATCg                                                                                           <  2:320170/151‑1 (MQ=255)
                                                          cccGCCTCGCGGGCGCTGCGCGAAGCGGTTTACAAACCCATTATGTTTAACTGTACCCGGAAACCCGTCCCCCCCGACTACGACCCTCCGGACCGCCTCTACTTCGAGCCGGTAACTCTGGAAGATGTGCTGGAAATCGTGCGTAt                                                                                             <  1:542587/144‑1 (MQ=255)
                                                          cccGCCCCGGTGGCGCTGCGCGAAGGGGGTTACAAAACCATTTTGTTTAACTGTACCCGGAAACCCGTCTCCCCCGACTACGACCCTCCCGACCGCCTCTACTTCGAGCCGGTAACTCTGGAAGATGTGCTGGAAATCGTGCGTATCg                                                                                           <  1:372515/146‑1 (MQ=255)
                                                           acGCCTCGCTGGCGCTGCGCGAAGACGGTTACGAAACCATTATGGTTAACTGTAACCCGGAAACCGTCTCCACCGACTACGACACTTCCGACCGCCTCTACTTCGAGCCGGTAACTCTGGAAGATGTGCTGGAAATCGTGCGTATCgaga                                                                                        >  1:365767/1‑150 (MQ=255)
                                                             gCCTCGCTGGCGCTGCGCGAAGACGGTTACGAAACCATTATGGTTAACTGTAACCCGGAAACCGTCTCCACCGACTACGACACTTCCGACCGCCTCTACTTCGAGCCGGTAACTCTGGAAGATGTGCTGGAAATCGTGCGTATCGAGAAGc                                                                                     >  2:409825/1‑151 (MQ=255)
                                                                     ggCGCGGCGGGAAGAGGTTGAAAAACCCTTTGGGTTAAATGAAACCCGGAAACCCGTCCCCCCCGATAACGACACTCCCAACCGCCTCTACTTCGAGCCGGTAACTCTGGAAGATGTGCTGGAAATCGTGCGTATCGAGAAGCCGaaa                                                                                <  1:52986/148‑1 (MQ=255)
                                                                              cgAAGACGGTTACGAAACCATTATGGTTAACTGTAACCCGGAAACCGTCTCCACCGACTACGACACTTCCGACCGCCTCTACTTCGAGCCGGTAACTCTGGAAGATGTGCTGGAAATCGTGCGTATCGAGAAGCCGAAAGGCGTTATCGTc                                                                    >  2:192826/1‑151 (MQ=255)
                                                                                aaGACGGTTACGAAACCATTATGTTTACCTGTACCCCGAAAACCGTCCCCCCCGCCTACGACACTCCCGACCGCCTCTACTTCGAGCCGGTAACTCTGGAAGATGTGCTGGAAATCGTGCGTATCGAGAAGCCGAAAGGCGTTATCg                                                                      <  2:157079/147‑1 (MQ=255)
                                                                                       ttACAAAACCTTTATGGTTAACTGTACCCCGAAACCCTTCCCCCCCGACTACGACCCTCCGGACCGCCTCTACTTCGAGCCGGTAACTCTGGAAGATGTGCTGGAAATCGTGCGTATCGAGAAGCCGAAAGGCGTTATc                                                                       <  1:173341/139‑1 (MQ=255)
                                                                                        tACGAAACCATTATGGTTAACTGTAACCCGGAAACCGTCTCCACCGACTACGACACTTCCGACCGCCTCTACTTCGAGCCGGTAACTCTGGAAGATGTGCTGGAAATCGTGCGTATCGAGAAGCCGAAAGGCGTTATCGTCCAGTAcggcg                                                          >  1:319192/1‑151 (MQ=255)
                                                                                            aacccATTAGGTTAAACTGTACCCGGAAACCCGTCCCCACCGACCACGACACTCCCGACCGCCTCTACTTCGAGCCGGTAACTCTGGAAGATGTGCTGGAAATCGTGCGTATCGAGAAGCCGAAAGGCGTTATCGTCCAGTACGGCGGTCa                                                      <  1:670761/148‑1 (MQ=255)
                                                                                            aaaCCATTATGTTTACCTGTACCCCGGAACCCGTCCCCCCCGACTACGACCCTCCCGACCGCCTCTACTTCGAGCCGGTAACTCTGGAAGATGTGCTGGAAATCGTGCGTATCGAGAAGCCGAAAGGCGTTATCGTCCAGTAcggc                                                           <  1:192805/146‑1 (MQ=255)
                                                                                            aaaCCATTATGGTTAACTGTAACCCGGAAACCGTCTCCACCGACTACGACACTTCCGACCGCCTCTACTTCGAGCCGGTAACTCTGGAAGATGTGCTGGAAATCGTGCGTATCGAGAAGCCGAAAGGCGTTATCGTCCAGTACGGCGGTc                                                       >  1:617495/1‑150 (MQ=255)
                                                                                              aCCATTATGGTTAACTGTAACCCGTAAACCGTCTCCACCGACTACGACACTTCCGACCGCCTCTACTTCGAGCCGGTAACTCTGGAAGATGTGCTGGAAATCGTGCGTATCGAGAAGCCGAAAGGCGTTATCGTCCAGTAcggcgg                                                         >  1:192408/1‑146 (MQ=255)
                                                                                                 aTTATGGTTAACTGTAACCCGGAAACCGTCTCCACCGACTACGACACTTCCGACCGCCTCTACTTCGAGCCGGTAACTCTGGAAGATGTGCTGGAAATCGTGCGTATCGAGAAGCCGAAAGGCGTTATCGTCCAGTACGGCGGTCAGAccc                                                 >  1:653560/1‑151 (MQ=255)
                                                                                                     tGGTTAACTGTAACCCGGAAACCGTCTCCACCGACTACGACACTTCCGACCGCCTCTACTTCGAGCCGGTAACTCTGGAAGATGTGCTGGAAATCGTGCGTATCGAGAAGCCGAAAGGCGTTATCGTCCAGTACGGCGGTCAGACCCCGCt                                             >  2:532027/1‑151 (MQ=255)
                                                                                                       gTTAACTGTAACCCGGAAACCGTCTCCACCGACTACGACACTTCCGACCGCCTCTACTTCGAGCCGGTAACTCTGGAAGATGTGCTGGAAATCGTGCGTATCGAGAAGCCGAAAGGCGTTATCGTCCAGTACGGCGGTCAGACCCCGCt                                             >  1:93852/1‑149 (MQ=255)
                                                                                                             ctgaaCCCGGGAACCCGCCCCCCCCGACAACGCCCCTCCCGACCGCCTCTACTTCGAGCCGGTAACTCTGGAAGATGTGCTGGAAATCGTGCGTATCGAGAAGCCGAAAGGCGTTATCGTCCAGTACGGCGGTCAGACCCCGCTGaa                                          <  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)
                                                                                                                                                   |                                                                                                                                                     
GCTTGGCGGCGGCCCGAACCGTATCGGTCAGGGTATCGAATTCGACTACTGTTGCGTACACGCCTCGCTGGCGCTGCGCGAAGACGGTTACGAAACCATTATGGTTAACTGTAACCCGGAAACCGTCTCCACCGACTACGACACTTCCGACCGCCTCTACTTCGAGCCGGTAACTCTGGAAGATGTGCTGGAAATCGTGCGTATCGAGAAGCCGAAAGGCGTTATCGTCCAGTACGGCGGTCAGACCCCGCTGAAACTGGCGCGCGCGCTGGAAGCTGCTGGCGTACCGGTTATC  >  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: