Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_000913646,7650GT42.4% 43.1 / 22.6 33Q249K (CAA→AAA) citG2‑(5''‑triphosphoribosyl)‑3'‑dephosphocoenzyme‑A synthase
Reads supporting (aligned to +/- strand):  ref base G (3/16);  new base T (14/0);  total (17/16)
Fisher's exact test for biased strand distribution p-value = 9.77e-07
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

TAATTAAATCTGTGCTAAAAACCAGGTAAGGATCAGTAGGTCAGCACTGCCGCCTGGACTGAGATTTCGTTCGATACACTCCCTGTCGAACTGCCGGAGATAATCGAGATCGGCGGGGGTTCGAATGCCCCCTTTTTGCAATAATGTTTGCGCCTCGCGCTGTAGCCAGCGCAGGCCCCCCTCGCCACCGCGCGATGCAACGTTGGTATCGCCGTTGATCGCCATCAGTAGGAGCAAGGTATCGAGCAATGCCAGTTCAGGATCTAACCCCTGATCCAGCAGAGTG  >  NC_000913/646628‑646913
                                                                                                                                         |                                                                                                                                                    
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:225636/151‑1 (MQ=255)
                aaaaaCCAGGTAAGGATCAGTAGGTCAGCACTGCCGCCTGGACTGAGATTTCGTTCGATACACTCCCTGTCGAACTGCCGGAGATAATCGAGATCGGCGGGGGTTCGAATGCCCCCTTTTTGCAATAATGTTTGCGCCTCGCGCTGTAGcc                                                                                                                         <  2:1688/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)
                                                                                                                            aTGCCCCCTTGTTGCAATAATGGTTGCGCCTCGCGCTGTAGCCAGCGCAGGCCCCCCTCGCCACCGCGCGATGCAACGTTGGTATCGCCGTTGATCGCCATCAGTAGGAGCAAGGTATCGAGCAATGCCAGTTCAGGATCTATCCCCTGAt             <  2:8177/151‑1 (MQ=255)
                                                                                                                              gCCCCCTTTTTGCAATAATGTTTGCGCCTCGCGCTGTAGCCAGCGCAGGCCCCCCTCGCCACCGCGCGATGCAACGTTGGTATCGCCGTTGATCGCCATCAGTAGGAGCAAGGTATCGAGCAATGCCAGTTCAGGATCTAACCCCTGATcc           >  2:300231/1‑151 (MQ=255)
                                                                                                                                     ttttGCAATAATGTTTGCGCCTCGCGCTGTAGCCAGCGCAGGCCCCCCTCGCCACCGCGCGATGCAACGTTGGTATCGCCGTTGATCGCCATCAGTAGGAGCAAGGTATCGAGCAATGCCAGTTCAGGATCTAACCCCTGATCCAGCaga     <  1:402313/150‑1 (MQ=255)
                                                                                                                                        tGCAATAATGTTTGCGCCTCGCGCTGTAGCCAGCGCAGGCCCCCCTCGCCACCGCGCGATGCAACGTTGGTATCGCCGTTGATCGCCATCAGTAGGAGCAAGGTATCGAGCAATGCCAGTTCAGGATCTAACCCCTGATCCAGCAGAGTg  <  2:231629/150‑1 (MQ=255)
                                                                                                                                         |                                                                                                                                                    
TAATTAAATCTGTGCTAAAAACCAGGTAAGGATCAGTAGGTCAGCACTGCCGCCTGGACTGAGATTTCGTTCGATACACTCCCTGTCGAACTGCCGGAGATAATCGAGATCGGCGGGGGTTCGAATGCCCCCTTTTTGCAATAATGTTTGCGCCTCGCGCTGTAGCCAGCGCAGGCCCCCCTCGCCACCGCGCGATGCAACGTTGGTATCGCCGTTGATCGCCATCAGTAGGAGCAAGGTATCGAGCAATGCCAGTTCAGGATCTAACCCCTGATCCAGCAGAGTG  >  NC_000913/646628‑646913

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