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

AATATTTCCAAATTTATCAAGTGCTTAAATAATTAAATCTGTGCTAAAAACCAGGTAAGGATCAGTAGGTCAGCACTGCCGCCTGGACTGAGATTTCGTTCGATACACTCCCTGTCGAACTGCCGGAGATAATCGAGATCGGCGGGGGTTCGAATGCCCCCTTTTTGCAATAATGTTTGCGCCTCGCGCTGTAGCCAGCGCAGGCCCCCCTCGCCACCGCGCGATGCAACGTTGGTATCGCCGTTGATCGCCATCAGTAGGAGCAAGGTATCGAGCAATGCCAGTTCAGGATCT  >  NC_000913/646599‑646892
                                                                                                                                                    |                                                                                                                                                 
aaTATTTCCAAATTTATCAAGTGCTTAAATAATTAAATCTGTGCTAAAAACCAGGTAAGGATCAGTAGGTCAGCACTGCCGCCTGGACTGAGATTTCGTTCGATACACTCCCTGTCGAACTGCCGGAGATAATCGAGATCGGCGGGGGtt                                                                                                                                                  <  1:234817/150‑1 (MQ=255)
        cAAATTTATCAAGTGCTTAAATAATTAAATCTGTGCTAAAAACCAGGTAAGGATCAGTAGGTCAGCACTGCCGCCTGGACTGAGATTTCGTTCGATACACTCCCTGTCGAACTGCCGGAGATAATCGAGATCGGCGGGGGGTCGAATGccc                                                                                                                                         >  1:62125/1‑151 (MQ=255)
               aTCAAGTGCTTAAATAATTAAATCTGTGCTAAAAACCAGGTAAGGATCAGTAGGTCAGCACTGCCGCCTGGACTGAGATTTCGTTCGATACACTCCCTGTCGAACTGCCGGAGATAATCGAGATCGGCGGTGGGTCGAAATCCCCCctttt                                                                                                                                  >  1:122158/1‑151 (MQ=255)
               aTCAAGTGCTTAAATAATTAAATCTGTGCTAAAAACCAGGTAAGGATCAGTAGGTCAGCACTGCCGCCTGGACTGAGATTTCGTTCGATACACTCCCTGTCGAACTGCCGGAGATAATCGAGATCGGCGGGGGTTCGAATGCCCCCttttt                                                                                                                                  <  1:450559/151‑1 (MQ=255)
               aTCAAGTGCTTAAATAATTAAATCTGTGCTAAAAACCAGGTAAGGATCAGTAGGTCAGCACTGCCGCCTGGACTGAGATTTCGTTCGATACACTCCCTGTCGAACTGCCGGAGATAATCGAGATCGGCGGGGGGTTGGATGCCCCCctttt                                                                                                                                  >  2:149832/1‑151 (MQ=255)
                  aaGTGCTTAAATAATTAAATCTGTGCTAAAAACCAGGTAAGGATCAGTAGGTCAGCACTGCCGCCTGGACTGAGATTTCGTTCGATACACTCCCTGTCGAACTGCCGGAGATAATCGAGATCGGCGGGGGGTCGAATGCCCCCctttt                                                                                                                                  >  1:613891/1‑148 (MQ=255)
                     tGCTTAAATAATTAAATCTGTGCTAAAAACCAGGTAAGGATCAGTAGGTCAGCACTGCCGCCTGGACTGAGATTTCGTTCGATACACTCCCTGTCGAACTGCCGGAGATAATCGAGATCGGCGGGGGTTTTAAGTCCCCCtcttttct                                                                                                                               >  1:699228/1‑145 (MQ=255)
                             taatTAAATCTGTGCTAAAAACCAGGTAAGGATCAGTAGGTCAGCACTGCCGCCTGGACTGAGATTTCGTTCGATACACTCCCTGTCGAACTGCCGGAGATAATCGAGATCGGCGGGGGGTCGAATGCCCCCCTTTTTCAATAAAGttt                                                                                                                      >  1:689014/1‑149 (MQ=255)
                              aatTAAATCTGTGCTAAAAACCAGGTAAGGATCAGTAGGTCAGCACTGCCGCCTGGACTGAGATTTCGTTCGATACACTCCCTGTCGAACTGCCGGAGATAATCGAGATCGGCGGGGGTTCGAATGCCCCCTTTTTGCAATAATGTTTgcg                                                                                                                   <  2:94070/151‑1 (MQ=255)
                              aatTAAATCTGTGCTAAAAACCAGGTAAGGATCAGTAGGTCAGCACTGCCGCCTGGACTGAGATTTCGTTCGATACACTCCCTGTCGAACTGCCGGAGATAATCGAGATCGGCGGGGGTTCGAATGCCCCCTTTTTGCAATAATGTTTgc                                                                                                                    <  2:677772/150‑1 (MQ=255)
                                    aTCTGTGCTAAAAACCAGGTAAGGATCAGTAGGTCAGCACTGCCGCCTGGACTGAGATTTCGTTCGATACACTCCCTGTCGAACTGCCGGAGATAATCGAGATCGGCGGGGGGTCGCATGCCCCCCTTTTTCAATAAAGTTTTTGCCTc                                                                                                               >  1:754288/1‑149 (MQ=255)
                                     tCTGTGCTAAAAACCAGGTAAGGATCAGTAGGTCAGCACTGCCGCCTGGACTGAGATTTCGTTCGATACACTCCCTGTCGAACTGCCGGAGATAATCGAGATCGGCGGGGGGTCGAATGCCCCCCTTTTTCAATAATGTTTTCGcccc                                                                                                               >  1:157199/1‑146 (MQ=255)
                                         tgCTAAAAACCAGGTAAGGATCAGTAGGTCAGCACTGCCGCCTGGACTGAGATTTCGTTCGATACACTCCCTGTCGAACTGCCGGAGATAATCGAGATCGGCGGGGGGTCGAATGCCCCCCTTTTTCAATAAATTTTTCGcccc                                                                                                               >  2:444377/1‑142 (MQ=255)
                                          gCTAAAAACCAGGTAAGGATCAGTAGGTCAGCACTGCCGCCTGGACTGAGATTTCGTTCGATACACTCCCTGTCGAACTGCCGGAGATAATCGAGATCGGCGGGGGTTCGAATGCCCCCTTTTTGCAATAATGTTTGCGCCTCGCGCTGTa                                                                                                       <  2:166031/151‑1 (MQ=255)
                                             aaaaaCCAGGTAAGGATCAGTAGGTCAGCACTGCCGCCTGGACTGAGATTTCGTTCGATACACTCCCTGTCGAACTGCCGGAGATAATCGAGATCGGCGGGGGTTCGAATGCCCCCTTTTTGCAATAATGTTTGCGCCTCGCGCTGTAg                                                                                                      <  1:69972/149‑1 (MQ=255)
                                             aaaaaCCAGGTAAGGATCAGTAGGTCAGCACTGCCGCCTGGACTGAGATTTCGTTCGATACACTCCCTGTCGAACTGCCGGAGATAATCGAGATCGGCGGGGGTTCGAATGCCCCCTTTTTGCAATAATGTTTGCGCCTCGCGCTGTAGcc                                                                                                    <  2:1688/151‑1 (MQ=255)
                                             aaaaaCCAGGTAAGGATCAGTAGGTCAGCACTGCCGCCTGGACTGAGATTTCGTTCGATACACTCCCTGTCGAACTGCCGGAGATAATCGAGATCGGCGGGGGTTCGAATGCCCCCTTTTTGCAATAATGTTTGCGCCTCGCGCTGTAGcc                                                                                                    <  2:225636/151‑1 (MQ=255)
                                                  ccAGGTAAGGATCAGTAGGTCAGCACTGCCGCCTGGACTGAGATTTCGTTCGATACACTCCCTGTCGAACTGCCGGAGATAATCGAGATCGGCGGGGGGTCGAATGCCCCCCTTTTGCAATAATGTTTTCGCCCCCCGCTGTAGCCCgc                                                                                                 >  1:202382/1‑149 (MQ=255)
                                                         aGGATCAGTAGGTCAGCACTGCCGCCTGGACTGAGATTTCGTTCGATACACTCCCTGTCGAACTGCCGGAGATAATCGAGATCGGCGGGGGGTTGAATGGCCCCCTTTTTCAATAATGTTTTTGCCTCGCGCCGTAGCCAGCGCAGGcccc                                                                                        >  2:95893/1‑151 (MQ=255)
                                                               aGTAGGTCAGCACTGCCGCCTGGACTGAGATTTCGTTCGATACACTCCCTGTCGAACTGCCGGAGATAATCGAGATCGGCGGGGGTTCGAATGCCCCCTTTTTGCAATAATGTTTGCGCCTCGCGCTGTAGCCAGCGCAGGCCCCCCTCGc                                                                                  <  1:490702/151‑1 (MQ=255)
                                                                  aGGTCAGCACTGCCGCCTGGACTGAGATTTCGTTCGATACACTCCCTGTCGAACTGCCGGAGATAATCGAGATCGGCGGGGGGTCGGATGCCCCCCTTTTTCAAAAAAATTTTGGCCCCCCGCCGTTGCCAACGCCGGccccccccccc                                                                                 >  2:736904/1‑144 (MQ=255)
                                                                    gTCAGCACTGCCGCCTGGACTGAGATTTCGTTCGATACACTCCCTGTCGAACTGCCGGAGATAATCGAGATCGGCGGGGGGTCGAATGCCCCCCTTTTTCAATAATGTTTTCGCCTCGCGCTGTAGCCAGCGCAGGccccccccccccccc                                                                             >  1:582163/1‑142 (MQ=255)
                                                                          aCTGCCGCCTGGACTGAGATTTCGTTCGATACACTCCCTGTCGAACTGCCGGAGATAATCGAGATCGGCGGGGGTTCGAATGCCCCCTTTTTGCAATAATGTTTGCGCCTCGCGCTGTAGCCAGCGCAGGCCCCCCTCGCCACCGCGCGAt                                                                       <  1:474162/151‑1 (MQ=255)
                                                                            tGCCGCCTGGACTGAGATTTCGTTCGATACACTCCCTGTCGAACTGCCGGAGATAATCGAGATCGGCGGGGGGTCGAATGCCCCCCTTTTTTAATAATGTTTTCGCCCCGCGCCGTAGCCAGCGCAGGGCCCCCCCCCCCCCCCGCGATg                                                                      >  1:78504/1‑150 (MQ=255)
                                                                             gccgccTGGACTGAGATTTCGTTCGATACACTCCCTGTCGAACTGCCGGAGATAATCGAGATCGGCGGGGGTTCGAATGCCCCCTTTTTGCAATAATGTTTGCGCCTCGCGCTGTAGCCAGCGCAGGCCCCCCTCGCCACCGCGCGATg                                                                      <  1:494675/149‑1 (MQ=255)
                                                                                 ccTGGACTGAGATTTCGTTCGATACACTCCCTGTCGAACTGCCGGAGATAATCGAGATCGGCGGGGGTTCGAATGCCCCCTTTTTGCAATAATGTTTGCGCCTCGCGCTGTAGCCAGCGCAGGCCCCCCTCGCCACCGCGCGATGCAACGt                                                                <  1:730164/151‑1 (MQ=255)
                                                                                      aCTGAGATTTCGTTCGATACACTCCCTGTCGAACTGCCGGAGATAATCGAGATCGGCGGGGGGTCGAATGCCCCCCTTTTTCAATAAATTTTTCGCCTCTCGCCGTAATCCGCGCGGGCCCCCCTcccc                                                                                 >  1:141296/1‑129 (MQ=255)
                                                                                              ttcgttcgATACACTCCCTGTCGAACTGCCGGAGATAATCGAGATCGGCGGGGGTTCGAATGCCCCCTTTTTGCAATAATGTTTGCGCCTCGCGCTGTAGCCAGCGCAGGCCCCCCTCGCCACCGCGCGATGCAACGTTGGTATCGCCGt                                                    <  2:287432/150‑1 (MQ=255)
                                                                                              ttcgttcgATACACTCCCTGTCGAACTGCCGGAGATAATCGAGATCGGCGGGGGGTCGAATGCCCCCCTTTTTCAATAATGTTTTCGCCTCGCCCCGTAGCCAGCGCAGGCCCCCCCCCCCCCCCCCCGATGCAACGGTGGTATTGCCggt                                                   >  1:333629/1‑149 (MQ=255)
                                                                                                   tcgATACACTCCCTGTCGAACTGCCGGAGATAATCGAGATCGGCGGGGGTTCGAATGCCCCCTTTTTGCAATAATGTTTGCGCCTCGCGCTGTAGCCAGCGCAGGCCCCCCTCGCCACCGCGCGATGCAACGTTGGTATCGCCGTTGATc                                               <  1:669158/150‑1 (MQ=255)
                                                                                                     gATACACTCCCTGTCGAACTGCCGGAGATAATCGAGATCGGCGGGGGGTCGAATGCCCCCCTTTTTCAAAAATTTTTTTGCCCCGCCCCGTAGACAGCGCAGGCCCCCCCCCCCCCCCCGCGATGCAACGTTGGGATAGCCGGTGATCGc                                             >  2:345344/1‑150 (MQ=255)
                                                                                                         cacTCCCTGTCGAACTGCCGGAGATAATCGAGATCGGCGGGGGGTCGAATGGCCCCCTTTTTTAATAAAATTTTCGCCCCGCGCCGTAGCCAGCGCAGGCCCCCCCCCCCCCCCCGCGATGCCACGTTGGGATCGCCCTTGGTCGcc                                            >  1:655838/1‑147 (MQ=255)
                                                                                                                   cGAACTGCCGGAGATAATCGAGATCGGCGGGGGTTCGAATGCCCCCTTTTTGCAATAATGTTTGCGCCTCGCGCTGTAGCCAGCGCAGGCCCCCCTCGCCACCGCGCGATGCAACGTTGGTATCGCCGTTGATCGCCATCAGTAGGAGCaa                              <  2:328654/151‑1 (MQ=255)
                                                                                                                         gCCGGAGATAATCGAGATCGGCGGGGGGTCGAATTCCCCCCTTTTTCAATAAAGTTTTTGCCCCCCGCCGTAGCCAGCGCAGGCCCCCCCCCCCCCCCCCCGCTGCAACGTTGGTATCGCCGTTGTTCGCCCTCACGTGGAGGAAggg                           >  2:3258/1‑147 (MQ=255)
                                                                                                                           cGGAGATAATCGAGATCGGCGGGGGGTCGAATGGCCCCCTTTTTCAATAAATTTTTCGCCCCGCGCTGTAGCCCGCGCAGGCCCCCCCCCCCCCCCCCCGCTGGAAAGTTGGTATTGCCGTTGGTCGCCCCCAGTAGGAGGAAGGGAttg                       >  2:620818/1‑148 (MQ=255)
                                                                                                                                               gggggTTCGAATGCCCCCTTTTTGCAATAATGTTTGCGCCTCGCGCTGTAGCCAGCGCAGGCCCCCCTCGCCACCGCGCGATGCAACGTTGGTATCGCCGTTGATCGCCATCAGTAGGAGCAAGGTATCGAGCAATGCCAGTTCAGGATCt  >  2:733205/1‑151 (MQ=255)
                                                                                                                                                    |                                                                                                                                                 
AATATTTCCAAATTTATCAAGTGCTTAAATAATTAAATCTGTGCTAAAAACCAGGTAAGGATCAGTAGGTCAGCACTGCCGCCTGGACTGAGATTTCGTTCGATACACTCCCTGTCGAACTGCCGGAGATAATCGAGATCGGCGGGGGTTCGAATGCCCCCTTTTTGCAATAATGTTTGCGCCTCGCGCTGTAGCCAGCGCAGGCCCCCCTCGCCACCGCGCGATGCAACGTTGGTATCGCCGTTGATCGCCATCAGTAGGAGCAAGGTATCGAGCAATGCCAGTTCAGGATCT  >  NC_000913/646599‑646892

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