Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_000913646,7470TG45.2% 52.7 / 12.1 42T255P (ACC→CCC) citG2‑(5''‑triphosphoribosyl)‑3'‑dephosphocoenzyme‑A synthase
Reads supporting (aligned to +/- strand):  ref base T (2/21);  new base G (19/0);  total (21/21)
Fisher's exact test for biased strand distribution p-value = 9.40e-10
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

ATATTTCCAAATTTATCAAGTGCTTAAATAATTAAATCTGTGCTAAAAACCAGGTAAGGATCAGTAGGTCAGCACTGCCGCCTGGACTGAGATTTCGTTCGATACACTCCCTGTCGAACTGCCGGAGATAATCGAGATCGGCGGGGGTTCG‑AATGCCCCCTTTTTGCAATAATGTTTGCGCCTCGCGCTGTAGCCAGCGCAGGCCCCCCTCGCCACCGCGCGATGCAACGTTGGTATCGCCGTTGATCGCCATCAGTAGGAGCAAGGTATCGAGCAATGCCAGTTCAGGATC  >  NC_000913/646600‑646891
                                                                                                                                                   |                                                                                                                                                 
atatTTCCAAATTTATCAAGTGCTTAAATAATTAAATCTGTGCTAAAAACCAGGTAAGGATCAGTAGGTCAGCACTGCCGCCTGGACTGAGATTTCGTTCGATACACTCCCTGTCGAACTGCCGGAGATAATCGAGATCGGCGGGGGGTCg                                                                                                                                                >  1:249937/1‑151 (MQ=255)
           tttATCAAGTGCTTAAATAATTAAATCTGTGCTAAAAACCAGGTAAGGATCAGTAGGTCAGCACTGCCGCCTGGACTGAGATTTCGTTCGATACACTCCCTGTCGAACTGCCGGAGATAATCGAGATCGGCGGGGGGTTG‑GATGgcccc                                                                                                                                      >  2:440362/1‑149 (MQ=255)
             tATCAAGTGCTTAAATAATTAAATCTGTGCTAAAAACCAGGTAAGGATCAGTAGGTCAGCACTGCCGCCTGGACTGAGATTTCGTTCGATACACTCCCTGTCGAACTGCCGGAGATAATCGAGATCGGCGGGGGGGCG‑GATGGCCCCcttt                                                                                                                                  >  2:616436/1‑151 (MQ=255)
               tCAAGTGCTTAAATAATTAAATCTGTGCTAAAAACCAGGTAAGGATCAGTAGGTCAGCACTGCCGCCTGGACTGAGATTTCGTTCGATACACTCCCTGTCGAACTGCCGGAGATAATCGAGATCGGCGGGGGGTCG‑AATGCCCCCtttttt                                                                                                                                >  1:558941/1‑150 (MQ=255)
                    tGCTTAAATAATTAAATCTGTGCTAAAAACCAGGTAAGGATCAGTAGGTCAGCACTGCCGCCTGGACTGAGATTTCGTTCGATACACTCCCTGTCGAACTGCCGGAGATAATCGAGATCGGCGGGGGGTTG‑AATGACCCCTTTTTTCaa                                                                                                                             >  1:178651/1‑149 (MQ=255)
                      cTTAAATAATTAAATCTGTGCTAAAAACCAGGTAAGGATCAGTAGGTCAGCACTGCCGCCTGGACTGAGATTTCGTTCGATACACTCCCTGTCGAACTGCCGGAGATAATCGAGATCGGCGGGGGGTCG‑GATGGCCCCCTTTTGCaataa                                                                                                                          >  1:180526/1‑150 (MQ=255)
                         aaaTAATTAAATCTGTGCTAAAAACCAGGTAAGGATCAGTAGGTCAGCACTGCCGCCTGGACTGAGATTTCGTTCGATACACTCCCTGTCGAACTGCCGGAGATAATCGAGATCGGCGGGGGTTTG‑AATTCCCCCCTTTTTCaaaaa                                                                                                                          >  1:434523/1‑147 (MQ=255)
                          aataatTAAATCTGTGCTAAAAACCAGGTAAGGATCAGTAGGTCAGCACTGCCGCCTGGACTGAGATTTCGTTCGATACACTCCCTGTCGAACTGCCGGAGATAATCGAGATCGGCGGGGGTTCGAAATCCCCCCTTTTGGCattaagttt                                                                                                                      >  1:179726/1‑147 (MQ=255)
                             aatTAAATCTGTGCTAAAAACCAGGTAAGGATCAGTTGGTCAGCACTGCCGCCTGGACTGAGATTTCGTTCGATACACTCCCTGTCGAACTGCCGGAGATAATCGAGATCGGCGGGGGTTCG‑AATGCCCCCTTTTTGCAATAATGTTTgc                                                                                                                   <  2:514913/150‑1 (MQ=255)
                             aatTAAATCTGTGCTAAAAACCAGGTAAGGATCAGTAGGTCAGCACTGCCGCCTGGACTGAGATTTCGTTCGATACACTCCCTGTCGAACTGCCGGAGATAATCGAGATCGGCGGGGGTTCG‑AATGCCCCCTTTTTGCAATAATGTTTgcg                                                                                                                  <  2:238102/151‑1 (MQ=255)
                                tAAATCTGTGCTAAAAACCAGGTAAGGATCAGTAGGTCAGCACTGCCGCCTGGACTGAGATTTCGTTCGATACACTCCCTGTCGAACTGCCGGAGATAATCGAGATCGGCGGGGGGTCG‑GATGCCCCCCTTTTTTTATAAAGttttttcc                                                                                                                >  1:701322/1‑145 (MQ=255)
                                 aaaTCTGTGCTAAAAACCAGGTAAGGATCAGTAGGTCAGCACTGCCGCCTGGACTGAGATTTCGTTCGATACACTCCCTGTCGAACTGCCGGAGATAATCGAGATCGGCGGGGGTTCG‑AATGCCCCCTTTTTGCAATAATGTTTGCGCCt                                                                                                               <  1:84853/150‑1 (MQ=255)
                                  aaTCTGTGCTAAAAACCAGGTAAGGATCAGTAGGTCAGCACTGCCGCCTGGACTGAGATTTCGTTCGATACACTCCCTGTCGAACTGCCGGAGATAATCGAGATCGGCGGGGGTTCG‑AATGCCCCCTTTTTGCAATAATGTTTGCGCCTcg                                                                                                             <  1:182831/151‑1 (MQ=255)
                                  aaTCTGTGCTAAAAACCAGGTAAGGATCAGTAGGTCAGCACTGCCGCCTGGACTGAGAGTTCGTTCGATACACTCCCTGTCGAACTGCCGGAGAGAATCGAGATCGGCGGGGGTTCG‑AATGCCCCCTTTTTGCAATAATGTTTGCGCCTcg                                                                                                             <  2:340331/151‑1 (MQ=255)
                                   aTCTGTGCTAAAAACCAGGTAAGGATCAGTAGGTCAGCACTGCCGCCTGGACTGAGATTTCGTTCGATACACTCCCTGTCGAACTGCCGGAGATAATCGAGATCGGCGGGGGGGCC‑AATGCCCCCCTTTTTCAATAATggcttattccca                                                                                                             >  1:64953/1‑139 (MQ=255)
                                     cTGTGCTAAAAACCAGGTAAGGATCAGTAGGTCAGCACTGCCGCCTGGACTGAGATTTCGTTCGATACACTCCCTGTCGAACTGCCGGAGATAATCGAGATCGGCGGGGGGTCG‑AATGCCCCCCTTTTGCaataaatttt                                                                                                                     >  2:255153/1‑135 (MQ=255)
                                            aaaaaCCAGGTAAGGATCAGTAGGTCAGCACTGCCGCCTGGACTGAGATTTCGTTCGATACACTCCCTGTCGAACTGCCGGAGATAATCGAGATCGGCGGGGGTTCG‑AATGCCCCCTTTTTGCAATAATGTTTGCGCCTCGCGCTGTAGcc                                                                                                   <  1:676757/151‑1 (MQ=255)
                                                   aGGTAAGGATCAGTAGGTCAGCACTGCCGCCTGGACTGAGATTTCGTTCGATACACTCCCTGTCGAACTGCCGGAGATAATCGAGATCGGCGGGGGGTCG‑AATTCCCCCCTTTTTTCATAATGtttttgcccccc                                                                                                            >  1:668230/1‑126 (MQ=255)
                                                    ggTAAGGATCAGTAGGTCAGCACTGCCGCCTGGACTGAGATTTCGTTCGATACACTCCCTGTCGAACTGCCGGAGATAATCGAGATCGGCGGGGGGTGG‑AATGCCCCCCTTTTTCCATTATTTTTTCGGCCCCCGCTTTTGCCAgc                                                                                                >  1:534866/1‑146 (MQ=255)
                                                       aaGGATCAGTAGGTCAGCACTGCCGCCTGGACTGAGATTTCGTTCGATACACTCCCTGTCGAACTGCCGGAGATAATCGAGATCGGCGGGGGGTCG‑CATGCCCCCCTTTTGCAATAATGTTTTCGCCTCGCGCGGTAGCCAGCGCCGGgcc                                                                                        >  2:304877/1‑151 (MQ=255)
                                                            tCAGTAGGTCAGCACTGCCGCCTGGACTGAGATTTCGTTCGATACACTCCCTGTCGAACTGCCGGAGATAATCGAGATCGGCGGGGGGTTG‑GATGCCCCCTTTTTTCAATAATGTTTTCGCCTCGCGCTGTAGCCCGCGCAGGcccccccc                                                                                   >  1:68921/1‑149 (MQ=255)
                                                            tCAGTAGGTCAGCACTGCCGCCTGGACTGAGATTTCGTTCGATACACTCCCTGTCGAACTGCCGGAGATAATCGAGATCGGCGGGGGGTCG‑TAACCCCCCCTTTTTTGATAATGTTTATGGCTCGCGCCGGAGCCAGCGCAGGcccccccc                                                                                   >  1:648979/1‑149 (MQ=255)
                                                                  ggTCAGCACTGCCGCCTGGACTGAGATTTCGTTCGATACACTCCCTGTCGAACTGCCGGAGATAATCGAGATCGGCGGGGGTTCG‑AATGCCCCCTTTTTGCAATAATGTTTGCGCCTCGCGCTGTAGCCAGCGCAGGCCCCCCTCGCCAcc                                                                             <  2:56337/151‑1 (MQ=255)
                                                                     cAGCACTGCCGCCTGGACTGAGATTTCGTTCGATACACTCCCTGTCGAACTGCCGGAGATAATCGAGATCGGCGGGGGTTCG‑AATGCCCCCTTTTTGCAATAATGTTTGCGCCTCGCGCTGTAGCCAGCGCAGGCCCCCCTCGCCACcgc                                                                           <  2:436042/150‑1 (MQ=255)
                                                                     cAGCACTGCCGCCTGGACTGAGATTTCGTTCGATACACTCCCTGTCGAACTGCCGGAGATAATCGAGATCGGCGGGGGGTCG‑GAGGCCCCCTTTTTTCAATAATGTTTTCGCCTCGCCGTGCAGACAGCGGCGGccccccccccccccccc                                                                          >  2:472191/1‑140 (MQ=255)
                                                                      aGCACTGCCGCCTGGACTGAGATTTCGTTCGATACACTCCCTGTCGAACTGCCGGAGATAATCGAGATCGGCGGGGGTTCG‑AATGCCCCCTTTTTGCAATAATGTTTGCGCCTCGCGCTGTAGCCAGCGCAGGCCCCCCTCGCCACcgcgc                                                                         <  1:89962/151‑1 (MQ=255)
                                                                      aGCACTGCCGCCTGGACTGAGATTTCGTTCGATACACTCCCTGTCGAACTGCCGGAGATAATCGAGATCGGCGGGGGGTTG‑AATGGCCCCCTTTTTTAATAATTTTTTTGCCTCCCGCCGTAGCCAGCGCCGGCCCCCCCCGCCAcccccc                                                                         >  1:417388/1‑147 (MQ=255)
                                                                            gccgccTGGACTGAGATTTCGTTCGATACACTCCCTGTCGAACTGCCGGAGATAATCGAGATCGGCGGGGGTTCG‑AATGCCCCCTTTTTGCAATAATGTTTGCGCCTCGCGCTGTAGCCAGCGCAGGCCCCCCTCGCCACCGCGCGATg                                                                     <  2:434296/149‑1 (MQ=255)
                                                                                         agaTTTCGTTCGATACACTCCCTGTCGAACTGCCGGAGATAATCGAGATCGGCGGGGGTTCG‑AATGCCCCCTTTTTGCAATAATGTTTGCGCCTCGCGCTGTAGCCAGCGCAGGCCCCCCTCGCCACCGCGCGATGCAACGTTGGTATCGc                                                      <  1:506359/151‑1 (MQ=255)
                                                                                              tcgttcgATACACTCCCTGTCGAACTGCCGGAGATAATCGAGATCGGCGGGGGTTCG‑AATGCCCCCTTTTTGCAATAATGTTTGCGCCTCGCGCTGTAGCCAGCGCAGGCCCCCCTCGCCACCGCGCGATGCAACGTTGGTATCGCCGtt                                                  <  2:34244/150‑1 (MQ=255)
                                                                                               cgttcgATACACTCCCTGTCGAACTGCCGGAGATAATCGAGATCGGCGGGGGTTCG‑AATGCCCCCTTTTTGCAATAATGTTTGCGCCTCGCGCTGTAGCCAGCGCAGGCCCCCCTCGCCACCGCGCGATGCAACGTTGGTATCGCCGTTGa                                                <  2:188311/151‑1 (MQ=255)
                                                                                                gttcgATACACTCCCTGTCGAACTGCCGGAGATAATCGAGATCGGCGGGGGGTCG‑AATGCCCCCTTTTTGCAAAAATGTTTTTGCCTCGCGCTCTAGCCCGCGCAGGCCCCCCCCCCCACCGCGCGATACAAAGGTGGTATCGCCGGTGAt                                               >  1:632343/1‑151 (MQ=255)
                                                                                                   cgATACACTCCCTGTCGAACTGCCGGAGATAATCGAGATCGGCGGGGGTTCG‑AATGCCCCCTTTTTGCAATAATGTTTGCGCCTCGCGCTGTAGCCAGCGCAGGCCCCCCTCGCCACCGCGCGATGCAACGTTGGTATCGCCGTTGATCGc                                            <  1:415154/151‑1 (MQ=255)
                                                                                                    gATACACTCCCTGTCGAACTGCCGGAGATAATCGAGATCGGCGGGGGGTCG‑GAAGCCCCCCTTTTTCCATAATGTTTTCGCCTCCCGCTGTAGCCCGCGCAGGCCCCCCCCCCCCCCCCCCGATGCAACGGTGGTATCGCCGTTGATCGcc                                           >  1:427985/1‑151 (MQ=255)
                                                                                                        cacTCCCTGTCGAACTGCCGGAGATAATCGAGATCGGCGGGGGTTCG‑AATGCCCCCTTTTTGCAATAATGTTTGCGCCTCGCGCTGTAGCCAGCGCAGGCCCCCCTCGCCACCGCGCGATGCAACGTTGGTATCGCCGTTGATCGCCATCa                                       <  1:378807/151‑1 (MQ=255)
                                                                                                                  cGAACTGCCGGAGATAATCGAGATCGGCGGGGGGTCG‑GATGCCCCCTTTTTTCTATAATGTTTTCGCCCCCCGCTGGAGCCAGCGCAGGCCCCCCCCGCCCCCCCGCGATGCAACGTTGGTATCGCCGGTGATCGCCCTCCGTAgg                                  >  1:35924/1‑146 (MQ=255)
                                                                                                                     aCTGCCGGAGATAATCGAGATCGGCGGGGGTTCG‑AATGCCCCCTTTTTGCAATAATGTTTGCGCCTCGCGCTGTAGCCAGCGCAGGCCCCCCTCGCCACCGCGCGATGCAACGTTGGTATCGCCGTTGATCGCCATCAGTAGGAGCAAGGt                          <  2:52677/151‑1 (MQ=255)
                                                                                                                        gCCGGAGATAATCGAGATCGGCGGGGGTTCG‑AATGCCCCCTTTTTGCAATAATGTTTGCGCCTCGCGCTGTAGCCAGCGCAGGCCCCCCTCGCCACCGCGCGATGCAACGTTGGTATCGCCGTTGATCGCCATCAGTAGGAGCAAGGTa                         <  2:318184/149‑1 (MQ=255)
                                                                                                                        gCCGGAGATAATCGAGATCGGCGGGGGTTCG‑AATGCCCCCTTTTTGCAATAATGTTTGCGCCTCGCGCTGTAGCCAGCGCAGGCCCCCCTCGCCACCGCGCGATGCAACGTTGGTATCGCCGTTGATCGCCATCAGTAGGAGCAAGGTATc                       <  1:255197/151‑1 (MQ=255)
                                                                                                                         ccGGAGATAATCGAGATCGGCGGGGGTTCG‑AATGCCCCCTTTTTGCAATAATGTTTGCGCCTCGCGCTGTAGCCAGCGCAGGCCCCCCTCGCCACCGCGCGATGCAACGTTGGTATCGCCGTTGATCGCCATCAGTAGGAGCAAGGTATCg                      <  2:679274/151‑1 (MQ=255)
                                                                                                                          cGGAGATAATCGAGATCGGCGGGGGTTCG‑AATGCCCCCTTTTTGCAATAATGTTTGCGCCTCGCGCTGTAGCCAGCGCAGGCCCCCCTCGCCACCGCGCGATGCAACGTTGGTATCGCCGTTGATCGCCATCAGTAGGAGCAAGGTATCGa                     <  2:508912/151‑1 (MQ=255)
                                                                                                                                             cggGGGTTCG‑AATGCCCCCTTTTTGCAATAATGTTTGCGCCTCGCGCTGTAGCCAGCGCAGGCCCCCCTCGCCACCGCGCGATGCAACGTTGGTATCGCCGTTGATCGCCATCAGTAGGAGCAAGGTATCGAGCAATGCCAGTTCAGGATc  <  1:370168/151‑1 (MQ=255)
                                                                                                                                                   |                                                                                                                                                 
ATATTTCCAAATTTATCAAGTGCTTAAATAATTAAATCTGTGCTAAAAACCAGGTAAGGATCAGTAGGTCAGCACTGCCGCCTGGACTGAGATTTCGTTCGATACACTCCCTGTCGAACTGCCGGAGATAATCGAGATCGGCGGGGGTTCG‑AATGCCCCCTTTTTGCAATAATGTTTGCGCCTCGCGCTGTAGCCAGCGCAGGCCCCCCTCGCCACCGCGCGATGCAACGTTGGTATCGCCGTTGATCGCCATCAGTAGGAGCAAGGTATCGAGCAATGCCAGTTCAGGATC  >  NC_000913/646600‑646891

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