Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_000913646,7470TG34.6% 33.2 / 13.2 26T255P (ACC→CCC) citGtriphosphoribosyl‑dephospho‑CoA synthase
Reads supporting (aligned to +/- strand):  ref base T (1/16);  new base G (9/0);  total (10/16)
Fisher's exact test for biased strand distribution p-value = 3.20e-06
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00
Rejected as consensus: Frequency below/above cutoff threshold.

GTGCTTAAATAATTAAATCTGTGCTAAAAACCAGGTAAGGATCAGTAGGTCAGCACTGCCGCCTGGACTGAGATTTCGTTCGATACACTCCCTGTCGAACTGCCGGAGATAATCGAGATCGGCGGGGGTTCGAATGCCCCCTTTTTGCAATAATGTTTGCGCCTCGCGCTGTAGCCAGCGCAGGCCCCCCTCGCCACCGCGCGATGCAACGTTGGTATCGCCGTTGATCGCCATCAGTAGGAGCAAGGTATCGAGCAATGCC  >  NC_000913/646619‑646880
                                                                                                                                |                                                                                                                                      
gTGCTTAAATAATTAAATCTGTGCTAAAAACCAGGTAAGGATCAGTAGGTCAGCACTGCCGCCTGGACTGAGATTTCGTTCGATACACTCCCTGTCGAACTGCCGGAGATAATCGAGATCGGCGGGGGTTCGAATGcccc                                                                                                                            <  1:502804/140‑1 (MQ=255)
   cTTAAATAATTAAATCTGTGCTAAAAACCAGGTAAGGATCAGTAGGTCAGCACTGCCGCCTGGACTGAGATTTCGTTCGATACACTCCCTGTCGAACTGCCGGAGATAATCGAGATCGGCGGGGGTTCg                                                                                                                                     >  2:317577/1‑129 (MQ=255)
   cTTAAATAATTAAATCTGTGCTAAAAACCAGGTAAGGATCAGTAGGTCAGCACTGCCGCCTGGACTGAGATTTCGTTCGATACACTCCCTGTCGAACTGCCGGAGATAATCGAGATCGGCGGGGGTTCg                                                                                                                                     <  1:317577/129‑1 (MQ=255)
   cTTAAATAATTAAATCTGTGCTAAAAACCAGGTAAGGATCAGTAGGTCAGCACTGCCGCCTGGACTGAGATTTCGTTCGATACACTCCCTGTCGAACTGCCGGAGATAATCGAGATCGGCGGGGGGTCGGATGcccccct                                                                                                                         >  1:626658/1‑138 (MQ=255)
           atTAAATCTGTGCTAAAAACCAGGTAAGGATCAGTAGGTCAGCACTGCCGCCTGGACTGAGATTTCGTTCGATACACTCCCTGTCGAACTGCCGGAGATAATCGAGATCGGCGGGGGGTCGGATGCCCCCCTTTTTCaaa                                                                                                                 >  1:333919/1‑139 (MQ=255)
             tAAATCTGTGCTAAAAACCAGGTAAGGATCAGTAGGTCAGCACTGCCGCCTTGACTGAGATTTCGTTCGATACACTCCCTGTCGAACTGCCGGAGATAATCGAGATCGGCGGGGGTTCGAATGCCCCCTTTTTGCaataa                                                                                                               <  1:653435/140‑1 (MQ=255)
                     tgCTAAAAACCAGGTAAGGATCAGTAGGTCAGCACTGCCGCCTGGACTGAGATTTCGTTCGATACACTCCCTGTCGAACTGCCGGAGATAATCGAGATCGGCGGGGGGTCGAATGCCCCCCTTTTTCAATAATGttttcc                                                                                                       >  1:617181/1‑137 (MQ=255)
                     tgCTAAAAACCAGGTAAGGATCAGTAGGTCAGCACTGCCGCCTGGACTGAGATTTCGTTCGATACACTCCCTGTCGAACTGCCGGAGATAATCGAGATCGGCGGGGGGTCGAATGCCCCCCTTTTTCAATAATGTTTtcg                                                                                                       >  2:18453/1‑140 (MQ=255)
                      gCTAAAAACCAGGTAAGGATCAGTAGGTCAGCACTGCCGCCTGGACTGAGATTTCGTTCGATACACTCCCTGTCGAACTGCCGGAGATAATCGAGATCGGCGGGGGGTCGGATGCCCCCCTTTTTCAATAAAGtttttgc                                                                                                      >  2:690903/1‑136 (MQ=255)
                                     aGGATCAGTAGGTCAGCACTGCCGCCTGGACTGAGATTTCGTTCGATACACTCCCTGTCGAACTGCCGGAGATAATCGAGATCGGCGGGGGTTCGAATGCCCCCTTTTTGCAATAATGTTTGCGCCTCGCGCTGTAGCCa                                                                                       <  2:78276/140‑1 (MQ=255)
                                        aTCAGTAGGTCAGCACTGCCGCCTGGACTGAGATTTCGTTCGATACACTCCCTGTCGAACTGCCGGAGATAATCGAGATCGGCGGGGGTTCGAATGCCCCCTTTTTGCAATAATGTTTGCGCCTCGCGCTGTAGCCAgcg                                                                                    <  1:97246/140‑1 (MQ=255)
                                         tCAGTAGGTCAGCACTGCCGCCTGGACTGAGATTTCGTTCGATACACTCCCTGTCGAACTGCCGGAGATAATCGAGATCGGCGGGGGGGCCGATTTCCCCTTTTTTTAAAAAAATTTTCGCCCCCCGCCGTAGGCAgcgc                                                                                   >  1:751436/1‑140 (MQ=255)
                                          aggTAGGTCAGCACTGCCGCCTGGACTGAGATTTCGTTCGATACACTCCCTGTCGAACTGCCGGAGATAATCGAGATCGGCGGGGGTTCGAATGCCCCCTTTTTGCAATAATGTTTGCGCCTCGCGCTGTAGCCAGCGCa                                                                                  <  1:512235/138‑1 (MQ=255)
                                           aGTAGGTCAGCACTGCCGCCTGGACTGAGATTTCGTTCGATACACTCCCTGTCGAACTGCCGGAGATAATCGAGATCGGCGGGGGGTCGAAAGCCCCCCTTTTTCAATAATTTTTTCGCCCCCCGCCGTAGCCAGCGCAg                                                                                 >  1:563721/1‑140 (MQ=255)
                                            gTAGGTCAGCACTGCCGCCTGGACTGAGATTTCGTTCGATACACTCCCTGTCGAACTGCCGGAGATAATCGAGATCGGCGGGGGGTCGGATGCCCCCCTTTTTCAATAATATTTTCGCCTCCCGCTGTAGCCAGCGCAgg                                                                                >  1:639141/1‑140 (MQ=255)
                                                gTCAGCACTGCCGCCTGGACTGAGATTTCGTTCGATACACTCCCTGTCGAACTGCCGGAGATAATCGAGATCGGCGGGGGGTCGGATGGCCCCCTTTTTCAATAATTTTTTCGCCCCCCGCTGTAGCCAGCGCAGGcccc                                                                            >  1:485728/1‑140 (MQ=255)
                                                     cACTGCCGCCTGGACTGAGATTTCGTTCGATACACTCCCTGTCGAACTGCCGGAGATAATCGAGATCGGCGGGGGGTCGAATGCCCCCCTTTTTCAATAAATTTTTCGCCCCCCGCTGTAGCCCGCGCAGGccccccccc                                                                       >  1:30109/1‑137 (MQ=255)
                                                       cTGCCGCCTGGACTGAGATTTCGTTCGATACACTCCCTGTCGAACTGCCGGAGATAATCGAGATCGGCGGGGGGTTGCATGCCCCCCTTTTTGAATAATTTTTTCGCCCCGCGCTGTAGCCAGCGCACGccccccccccc                                                                     >  2:136185/1‑135 (MQ=255)
                                                       cTGCCGCCTGGACTGAGATTTCGTTCGATACACTCCCTGTCGAACTGCCGGAGATAATCGAGATCGGCGGGGGGGCGAATGGCCCCCTTTTTTTATTATGTTTTAGGCCCCCCCTGTTGGCAGCGGCGGCCCCCCCCGcc                                                                     >  2:512235/1‑140 (MQ=255)
                                                                ggACTGAGATTTCGTTCGATACACTCCCTGTCGAACTGCCGGAGATAATCGAGATCGGCGGGGGTTCGAATGCCCCCTTTTTGCAATAATGTTTGCGCCTCGCGCTGTAGCCAGCGCAGGCCCCCCTCGCCACCGCGCGa                                                            <  2:766757/140‑1 (MQ=255)
                                                                   cTGAGATTTCGTTCGATACACTCCCTGTCGAACTGCCGGAGATAATCGAGATCGGCGGGGGTTCGAATGCCCCCTTTTTGCAATAATGTTTGCGCCTCGCGCTGTAGCCAGCGCAGGCCCCCCTCGCCACCGCGCGATGc                                                         <  2:161485/140‑1 (MQ=255)
                                                                                        tCCCTGTCGAACTGCCGGAGATAATCGAGATCGGCGGGGGGTTGCAAGCCCCCCTTTTTCAATAATATTTTCGCCCCCCCCCGTAGCCAACGCAGGGCCCCCCCGCCACCCCGCGATGCCACCTTGGGATTGCCGTTGAt                                    >  2:55712/1‑140 (MQ=255)
                                                                                            tGTCGAACTGCCGGAGATAATCGAGATCGGCGGGGGTTCGAATGCCCCCTTTTTGCAATAATGTTTGCGCCTCGCGCTGTAGCCAGCGCAGGCCCCCCTCGCCACCGCGCGATGCAACGTTGGTATCGCCGTTGATCGcc                                <  1:257818/140‑1 (MQ=255)
                                                                                             gTCGAACTGCCGGAGATAATCGAGATCGGCGGGGGGTCGGATGCCCCCCTTTTTTAATAAAATTTTCGCCCCCCCCTGTAGCCAGCGCAGGCCCCCCCCCCCCCCCCGCGATGCCACGGTTGTATCGCCGCTGATTGccc                               >  2:257818/1‑139 (MQ=255)
                                                                                              tCGAACTGCCGGAGATAATCGAGATCGGCGGGGGTTCGAATGCCCCCTTTTTGCAATAATGTTTGCGCCTCGCGCTGTAGCCAGCGCAGGCCCCCCTCGCCACCGCGCGATGCAACGTTGGTATCGCCGTTGATCGCCAt                              <  1:136185/140‑1 (MQ=255)
                                                                                                  aCTGCCGGAGATAATCGAGATCGGCGGGGGTTCGAATGCCCCCTTTTTGCAATAATGTTTGCGCCTCGCGCTGTAGCCAGCGCAGGCCCCCCTCGCCACCGCGCGATGCAACGTTGGTATCGCCGTTGATCGCCATCAGt                          <  1:589898/140‑1 (MQ=255)
                                                                                                      ccGGAGATAATCGAGATCGGCGGGGGTTCGAATGCCCCCTTTTTGCAATAATGTTTGCGCCTCGCGCTGTAGCCAGCGCAGGCCCCCCTCGCCACCGCGCGATGCAACGTTGGTATCGCCGTTGATCGCCATCAGTAGGa                      <  1:95881/140‑1 (MQ=255)
                                                                                                         gagaTAATCGAGATCGGCGGGGGTTCGAATGCCCCCTTTTTGCAATAATGTTTGCGCCTCGCGCTGTAGCCAGCGCAGGCCCCCCTCGCCACCGCGCGATGCAACGTTGGTATCGCCGTTGATCGCCATCAGTAGGAGCa                   <  2:626658/140‑1 (MQ=255)
                                                                                                         gagaTAATCGAGATCGGCGGGGGTTCGAATGCCCCCTTTTTGCAATAATGTTTGCGCCTCGCGCTGTAGCCAGCGCAGGCCCCCCTCGCCACCGCGCGATGCAACGTTGGTATCGCCGTTGATCGCCATCAGTAGGAGCa                   <  1:18453/140‑1 (MQ=255)
                                                                                                         gagaTAATCGAGATCGGCGGGGGTTCGAATGCCCCCTTTTTGCAATAATGTTTGCGCCTCGCGCTGTAGCCAGCGCAGGCCACCCTCGCCACCGCGCGATGCAACGTTGGTATCGCCGTTGATCGCCATCAGTAGGAGCa                   <  1:676835/140‑1 (MQ=255)
                                                                                                              aaTCGAGATCGGCGGGGGTTCGAATGCCCCCTTTTTGCAATAATGTTTGCGCCTCGCGCTGTAGCCAGCGCAGGCCCCCCTCGCCACCGCGCGATGCAACGTTGGTATCGCCGTTGATCGCCATCAGTAGGAGCAAGGTa              <  2:333919/140‑1 (MQ=255)
                                                                                                                   agaTCGGCGGGGGTTCGAATGCCCCCTTTTTGCAATAATGTTTGCGCCTCGCGCTGTAGCCAGCGCAGGCCCCCCTCGCCACCGCGCGATGCAACGTTGGTATCGCCGTTGATCGCCAGCAGTAGGAGCAAGGTATCGAg         <  2:751436/140‑1 (MQ=255)
                                                                                                                     aTCGGCGGGGGGTCGAATGCCCCCCTTTTTTAATAAAGTTTTTGCCTCGCGCTGTAGCCAGCGCAGGCCCCCCCCCCCCCCCCCCGGTGCAACGGTGGTATCGCCGTTGGTCGCCCTCAGTAGGAGCAAGGGATCGAGCa       >  2:397304/1‑140 (MQ=255)
                                                                                                                       cggcggGGGGTCGGATGCCCCCCTTTTTCCAAAAATTTTTCGCCCCCCGCTGTAGCCAGCGCAGGCCCCCCCCCCCACCCCCCGGTGCAACGGTGGTGTTGCCGGTTGTCGCCCCCCGTAGGAGCAAGGGATCGAGCaaa     >  2:296472/1‑139 (MQ=255)
                                                                                                                          cggGGGGTCGAAGGCCCCCCTTTTTCAATAATGTTTTCGCCTCGCGCTGTAGCCAGCGCAGGCCCCCCCCCCCCCCCCGCGGTGCAACGTTGGTATCGCCGTTGATCGCCCTCAGTAGGAGCAAGGGATCGAGCAAAGcc  >  2:234784/1‑140 (MQ=255)
                                                                                                                                |                                                                                                                                      
GTGCTTAAATAATTAAATCTGTGCTAAAAACCAGGTAAGGATCAGTAGGTCAGCACTGCCGCCTGGACTGAGATTTCGTTCGATACACTCCCTGTCGAACTGCCGGAGATAATCGAGATCGGCGGGGGTTCGAATGCCCCCTTTTTGCAATAATGTTTGCGCCTCGCGCTGTAGCCAGCGCAGGCCCCCCTCGCCACCGCGCGATGCAACGTTGGTATCGCCGTTGATCGCCATCAGTAGGAGCAAGGTATCGAGCAATGCC  >  NC_000913/646619‑646880

Alignment Legend
Aligned base mismatch/match (shaded by quality score): ATCG/ATCG < 3 ≤ ATCG/ATCG < 15 ≤ ATCG/ATCG < 27 ≤ ATCG/ATCG < 33 ≤ ATCG/ATCG < 40 ≤ ATCG/ATCG
Unaligned base: atcg    Masked matching base: atcg    Alignment gap:     Deleted base: