Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_000913646,7470TG27.0% 61.3 / 16.3 37T255P (ACC→CCC) citGtriphosphoribosyl‑dephospho‑CoA synthase
Reads supporting (aligned to +/- strand):  ref base T (2/25);  new base G (10/0);  total (12/25)
Fisher's exact test for biased strand distribution p-value = 1.89e-07
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.80e-01
Rejected as consensus: Frequency below/above cutoff threshold.

TTATCAAGTGCTTAAATAATTAAATCTGTGCTAAAAACCAGGTAAGGATCAGTAGGTCAGCACTGCCGCCTGGACTGAGATTTCGTTCGATACACTCCCTGTCGAACTGCCGGAGATAATCGAGATCGGCGGGGGTTCGAATGCCCCCTTTTTGCAATAATGTTTGCGCCTCGCGCTGTAGCCAGCGCAGGCCCCCCTCGCCACCGCGCGATGCAACGTTGGTATCGCCGTTGATCGCCATCAGTAGGAGCAAGGTATCGAGCAATGCCAGTTCA  >  NC_000913/646612‑646886
                                                                                                                                       |                                                                                                                                            
ttATCAAGTGCTTAAATAATTAAATCTGTGCTAAAAACCAGGTAAGGATCAGTAGGTCAGCACTGCCGCCTGGACTGAGATTTCGTTCGATACACTCCCTGTCGAACTGCCGGAGATAATCGAGATCGGCGGGGGTTCGa                                                                                                                                          <  2:50603/140‑1 (MQ=255)
        tGCTTAAATAATTAAATCTGTGCTAAAAACCAGGTAAGGATCAGTAGGTCAGCACTGCCGCCTGGACTGAGATTTCGTTCGATACACTCCCTGTCGAACTGCCGGAGATAATCGAGATCGGCGGGGGGTCGAATGccccc                                                                                                                                 >  1:731686/1‑140 (MQ=255)
          cTTAAATAATTAAATCTGTGCTAAAAACCAGGTAAGGATCAGTAGGTCAGCACTGCCGCCTGGACTGAGATTTCGTTCGATACACTCCCTGTCGAACTGCCGGAGATAATCGAGATCGGCGGGGGtt                                                                                                                                             <  1:205076/127‑1 (MQ=255)
          cTTAAATAATTAAATCTGTGCTAAAAACCAGGTAAGGATCAGTAGGTCAGCACTGCCGCCTGGACTGAGATTTCGTTCGATACACTCCCTGTCGAACTGCCGGAGATAATCGAGATCGGCGGGGGTTCGAATGCCCCCtt                                                                                                                               >  1:305066/1‑140 (MQ=255)
          cTTAAATAATTAAATCTGTGCTAAAAACCAGGTAAGGATCAGTAGGTCAGCACTGCCGCCTGGACTGAGATTTCGTTCGATACACTCCCTGTCGAACTGCCGGAGATAATCGAGATCGGCGGGGGGTCGAATGcccccct                                                                                                                               >  1:292169/1‑138 (MQ=255)
             aaaTAATTAAATCTGTGCTAAAAACCAGGTAAGGATCAGTAGGTCAGCACTGCCGCCTGGACTGAGATTTCGTTCGATACACTCCCTGTCGAACTGCCGGAGATAATCGAGATCGGCGGGGGGTCGAATGCCCCCctttt                                                                                                                            >  1:78573/1‑140 (MQ=255)
                 aatTAAATCTGTGCTAAAAACCAGGTAAGGATCAGTAGGTCAGCACTGCCGCCTGGACTGAGATTTCGTTCGATACACTCCCTGTCGAACTGCCGGAGATAATCGAGATCGGCGGGGGTTCGAATGCCCCCTTTTTGCaa                                                                                                                        <  2:588866/140‑1 (MQ=255)
                    tAAATCTGTGCTAAAAACCAGGTAAGGATCAGTAGGTCAGCACTGCCGCCTGGACTGAGATTTCGTTCGATACACTCCCTGTCGAACTGCCGGAGATAATCGAGATCGGCGGGGGGTCGGATGCCCCCCTTTTGCaataa                                                                                                                     >  2:232415/1‑140 (MQ=255)
                      aaTCTGTGCTAAAAACCAGGTAAGGATCAGTAGGTCAGCACTGCCGCCTGGACTGAGATTTCGTTCGATACACTCCCTGTCGAACTGCCGGAGATAATCGAGATCGGCGGGGGTTTGAATGCCCCCCTTTTGTcataatt                                                                                                                   >  1:60254/1‑139 (MQ=255)
                           gtgCTAAAAACCAGGTAAGGATCAGTAGGTCAGCACTGCCGCCTGGACTGAGATTTCGTTCGATACACTCCCTGTCGAACTGCCGGAGATAATCGAGATCGGCGGGGGGTTGGATGGCCCCCTTTTTtaataaatttttt                                                                                                              >  2:345228/1‑133 (MQ=255)
                             gCTAAAAACCAGGTAAGGATCAGTAGGTCAGCACTGCCGCCTGGACTGAGATTTCGTTCGATACACTCCCTGTCGAACTGCCGGAGATAATCGAGATCGGCGGGGGGTCGAATGCCCCCCTTTTTCTATAAAGTTTtcgc                                                                                                            >  2:67014/1‑140 (MQ=255)
                                         gTAAGGATCAGTAGGTCAGCACTGCCGCCTGGACTGAGATTTCGTTCGATACACTCCCTGTCGAACTGCCGGAGATAATCGAGATCGGCGGGGGGTCGGATGCCCCCCTTTTTTTATAATTTTTGTGCCTCCCGCCGttg                                                                                                >  2:674164/1‑138 (MQ=255)
                                             ggATCAGTAGGTCAGCACTGCCGCCTGGACTGAGATTTCGTTCGATACACTCCCTGTCGAACTGCCGGAGATAATCGAGATCGGCGGGGGGTCTAATGCCCCCCTTTTTCAATAATTTTTTCGCCCCCCGCTGGAGCCAg                                                                                            >  2:623125/1‑140 (MQ=255)
                                                 cAGTAGGTCAGCACTGCCGCCTGGACTGAGATTTCGTTCGATACACTCCCTGTCGAACTGCCGGAGATAATCGAGATCGGCGGGGGTTCGAATGCCCCCTTTTTGCAATAATGTTTGCGCCTCGCGCTGTAGCCAGCGCa                                                                                        <  2:687673/140‑1 (MQ=255)
                                                   gTAGGTCAGCACTGCCGCCTGGACTGAGATTTCGTTCGATACACTCCCTGTCGAACTGCCGGAGATAATCGAGATCGGCGGGGGGGCGGATGCCCCCCTTTTTCAATAATTTTTTTGCCCCCCGCTGTAGGCCGCGCCgg                                                                                      >  2:311104/1‑140 (MQ=255)
                                                     aGGTCAGCACTGCCGCCTGGACTGAGATTTCGTTCGATACACTCCCTGTCGAACTGCCGGAGATAATCGAGATCGGCGGGGGTTCGAATGCCCCCTTTTTGCAATAATGTTTGCGCCTCGCGCTGTAGCCAGCGCAGGcc                                                                                    <  1:437146/140‑1 (MQ=255)
                                                     aGGTCAGCACTGCCGCCTGGACTGAGATTTCGTTCGATACACTCCCTGTCGAACTGCCGGAGATAATCGAGATCGGCGGGGGTTCGAATGCCCCCTTTTTGCAATAATGTTTGCGCCTCGCGCTGTAGCCAGCGCAGGcc                                                                                    <  1:345228/140‑1 (MQ=255)
                                                         cAGCACTGCCGCCTGGACTGAGATTTCGTTCGATACACTCCCTGTCGAACTGCCGGAGATAATCGAGATCGGCGGGGGGTCGGATGCCCCCCTTTTTTAATAATGTTTTCGCCCCCCCCCCTAGCCCGCGCCGGcccccc                                                                                >  2:96209/1‑140 (MQ=255)
                                                                   gccTGGACTGAGATTTCGTTCGATACACTCCCTGTCGAACTGCCGGAGATAATCGAGATCGGCGGGGGTTCGAATGCCCCCTTTTTGCAATAATGTTTGCGCCTCGCGCTGTAGCCAGCGCAGGCCCCCCTCGCCACcgc                                                                      <  1:201156/140‑1 (MQ=255)
                                                                    ccTGGACTGAGATTTCGTTCGATACACTCCCTGTCGAACTGCCGGAGATAATCGAGATCGGCGGGGGTTCGAATGCCCCCTTTTTGCAATAATGTTTGCGCCTCGCGCTGTAGCCAGCGCAGGCCCCCCTCGCCACcgcg                                                                     <  1:623125/140‑1 (MQ=255)
                                                                    ccTGGACTGAGATTTCGTTCGATACACTCCCTGTCGAACTGCCGGAGATAATCGAGATCGGCGGGGGTTCGAATGCCCCCTTTTTGCAATAATGTTTGCGCCTCGCGCTGTAGCCAGCGCAGGCCCCCCTCGCCACcgcg                                                                     <  2:305066/140‑1 (MQ=255)
                                                                    ccTGGACTGAGATTTCGTTCGATACACTCCCTGTCGAACTGCCGGAGATAATCGAGATCGGCGGGGGTTCGAATGCCCCCTTTTTGCAATAATGTTTGCGCCTCGCGCTGTAGCCAGCGCAGGCCCCCCTCGCCACcgcg                                                                     <  1:674164/140‑1 (MQ=255)
                                                                          cTGAGATTTCGTTCGATACACTCCCTGTCGAACTGCCGGAGATAATCGAGATCGGCGGGGGTTCGAATGCCCCCTTTTTGCAATAATGTTTGCGCCTCGCGCTGTAGCCAGCGCAGGCCCCCCTCGCCACCGCGCGATGc                                                               <  2:682173/140‑1 (MQ=255)
                                                                            gagaTTTCGTTCGATACACTCCCTGTCGAACTGCCGGAGATAATCGAGATCGGCGGGGGGGCGGATACGCCCCTTTTTGTATAATCTTTGCGCCTCCCGCCGGGGCCCCCGCAGGGCCCCCCCCCCCCcccgcgcatcga                                                             >  2:201156/1‑134 (MQ=255)
                                                                                              cTCCCTGTCGAACTGCCGGAGATAATCGAGATCGGCGGGGGTTCGAATGCCCCCTTTTTGCAATAATGTTTGCGCCTCGCGCTGTAGCCAGCGCAGGCCCCCCTCGCCACCGCGCGATGCAACGTTGGTATCGCCGTTGa                                           <  1:711740/140‑1 (MQ=255)
                                                                                                    gTCGAACTGCCGGAGATAATCGAGATCGGCGGTGGGTCGCATTTCCCCTTTTTTTAATAATTTTTTTCCCTCCCCCCCTCGGCAGCGCAGGCCCCCCCCCCCCCCCCCCGCTCCAACCGTGGTATACCCGCTGATCggcc                                     >  2:375159/1‑137 (MQ=255)
                                                                                                    gTCGAACTGCCGGAGATAATCGAGATCGGCGGGGGGTTGGATGCCCCCCTTTTTCAATAATTTTTTCGCCCCGCGCCGTAGCCCGCGCCGGCCCCCCCCCCCCCCCCCCGCCGCAACGGTGGGATCGCCCTTGGTccccc                                     >  2:701597/1‑139 (MQ=255)
                                                                                                    gTCGAACTGCCGGAGATAATCGAGATCGGCGGGGGGTTGGATGCCCCCCTTTTTCAATAAATTTTTCGCCCCCCCCCGTAGCCAGCGCCGGCCCCCCCCCCCCCCCCGCGCTGCAACCGTGGGATTGCCCTTGATCGccc                                     >  2:512933/1‑139 (MQ=255)
                                                                                                    gTCGAACTGCCGGAGATAATCGAGATCGGCGGGGGGTCGAATGCCCCCCTTTTTTAATAATGTTTTTGCCTCGCGCTGTAGCCAGCGCCGGCCCCCCCCCCCCCCCCGCGGTGCAACGTTGGGATTGCCGTTGGTCGccc                                     >  2:190760/1‑139 (MQ=255)
                                                                                                    gTCGAACTGCCGGAGATAATCGAGATCGGCGGGGGGTCGAATGCCCCCCTTTTTTAATAAATTTTTTGCCCCCCGCCGTAGGCAGCGCCGGCCCCCCCCCCCCCCCCGCGCTGCAACGCTGGTATCGCCGTTGATCGccc                                     >  2:776472/1‑139 (MQ=255)
                                                                                                    gTCGAACTGCCGGAGATAATCGAGATCGGCGGGGGGTCGAATGCCCCCCTTTTTCAATAATTTTTTTGCCCCCCGCCGTAGCCCGCGCCGGGCCCCCCCCCCCCCCCGCGGTGCAACGCTGGTATCGCCCCTGGTCGccc                                     >  2:573190/1‑139 (MQ=255)
                                                                                                    gTCGAACTGCCGGAGATAATCGAGATCGGCGGGGGGTCGAATGCCCCCCTTTTTCAATAAATTTTTCGCCCCGCGCCGTTGCCCGCGCCGGCCCCCCCCCCCCCCCCCCGGCGCCACCTTGGGATTGCCGTTGGTCGccc                                     >  2:715408/1‑139 (MQ=255)
                                                                                                         aCTGCCGGAGATAATCGAGATCGGCGGGGGTTCGAATGCCCCCTTTTTGCAATAATGTTTGCGCCTCGCGCTGTAGCCAGCGCAGGCCCCCCTCGCCACCGCGCGATGCAACGTTGGTATCGCCGTTGATCGCCATCAGt                                <  1:48149/140‑1 (MQ=255)
                                                                                                         aCTGCCGGAGATAATCGAGATCGGCGGGGGTTCGAATGCCCCCTTTTTGCAATAATGTTTGCGCCTCGCGCTGTAGCCAGCGCAGGCCCCCCTCGCCACCGCGCGATGCAACGTTGGTATCGCCGTTGATCGCCATCAGt                                <  2:339365/140‑1 (MQ=255)
                                                                                                           tGCCGGAGATAATCGAGATCGGCGGGGGGTTGGATGCCCCCCTTTTTTAATAATTTTTTCGCCCCCCGCTGTAGCCCGCGCAGGCCCCCCCCCCCCCCCCCCGCTGCAACGTTGGTATCGCCGTTGATCGCCCTCAGGAg                              >  2:258985/1‑140 (MQ=255)
                                                                                                           tGCCGGAGATAATCGAGATCGGCGGGGGGTCGGATGCCCCCCTTTTTCAATAATGTTTTTGCCCCGCGCCGTAGCCAGCGCCGGCCCCCCCCCCCCCCCCCCCGTGCAACCTTGGTGTCGCCGCTGGTCGCCCTccggaa                              >  2:166905/1‑135 (MQ=255)
                                                                                                             ccGGAGATAATCGAGATCGGCGGGGGTTCGAATGCCCCCTTTTTGCAATAATGTTTGCGCCTCGCGCTGTAGCCAGCGCAGGCCCCCCTCGCCACCGCGCGATGCAACGTTGGTATCGCCGTTGATCGCCATCAGTAGGa                            <  2:462218/140‑1 (MQ=255)
                                                                                                             ccGGAGATAATCGAGATCGGCGGGGGTTCGAATGCCCCCTTTTTGCAATAATGTTTGCGCCTCGCGCTGTAGCCAGCGCAGGCCCCCCTCGCCACCGCGCGATGCAACGTTGGTATCGCCGTTGATCGCCATCAGTAGGa                            <  1:776472/140‑1 (MQ=255)
                                                                                                             ccGGAGATAATCGAGATCGGCGGGGGGTTGGATGGCCCCCTTTTTCAATAATTTTTTCGCCTCGCGCTGGAGCCAGCGCAGGCCCCCCCCCCCCCCCCGCGCGGCAACGTTGGTATCGCCGTTGGTCGCCCTCaggaggg                            >  2:709123/1‑139 (MQ=255)
                                                                                                                gagaTAATCGAGATCGGCGGGGGTTCGAATGCCCCCTTTTTGCAATAATGTTTGCGCCTCGCGCTGTAGCCAGCGCAGGCCCCCCTCGCCACCGCGCGATGCAACGTTGGTATCGCCGTTGATCGCCATCAGTAGGAGCa                         <  1:239654/140‑1 (MQ=255)
                                                                                                                gagaTAATCGAGATCGGCGGGGGTTCGAATGCCCCCTTTTTGCAATAATGTTTGCGCCTCGCGCTGTAGCCAGCGCAGGCCCCCCTCGCCACCGCGCGATGCAACGTTGGTATCGCCGTTGATCGCCATCAGTAGGAGCa                         <  2:265407/140‑1 (MQ=255)
                                                                                                                gagaTAATCGAGATCGGCGGGGGTTCGAATGCCCCCTTTTTGCAATAATGTTTGCGCCTCGCGCTGTAGCCAGCGCAGGCCCCCCTCGCCACCGCGCGATGCAACGTTGGTATCGCCGTTGATCGCCATCAGTAGGAGCa                         <  2:25036/140‑1 (MQ=255)
                                                                                                                gagaTAATCGAGATCGGCGGGGGTTCGAATGCCCCCTTTTTGCAATAATGTTTGCGCCTCGCGCTGTAGCCAGCGCAGGCCCCCCTCGCCACCGCGCGATGCAACGTTGGTATCGCCGTTGATCGCCATCAGTAGGAGCa                         <  1:701597/140‑1 (MQ=255)
                                                                                                                gagaTAATCGAGATCGGCGGGGGTTCGAATGCCCCCTTTTTGCAATAATGTTTGCGCCTCGCGCTGTAGCCAGCGCAGGCCCCCCTCGCCACCGCGCGATGCAACGTTGGTATCGCCGTTGATCGCCATCAGTAGGAGCa                         <  1:166905/140‑1 (MQ=255)
                                                                                                                  gaTAATCGAGATCGGCGGGGGTTCGAATGCCCCCTTTTTGCAATAATGTTTGCGCCTCGCGCTGTAGCCAGCGCAGGCCCCCCTCGCCACCGCGCGATGCAACGTTGGTATCGCCGTTGATCGCCATCAGTAGGAGCAAg                       <  1:67014/140‑1 (MQ=255)
                                                                                                                      aTCGAGATCGGCGGGGGTTCGAATGCCCCCTTTTTGCAATAATGTTTGCGCCTCGCGCTGTAGCCAGCGCAGGCCCCCCTCGCCACCGCGCGATGCAACGTTGGTATCGCCGTTGATCGCCATCAGTAGGAGCAAGGTAt                   <  1:258985/140‑1 (MQ=255)
                                                                                                                                gcggGGGTTCGAATGCCCCCTTTTTGCAATAATGTTTGCGCCTCGCGCTGTAGCCAGCGCAGGCCCCCCTCGCCACCGCGCGATGCAACGTTGGTATCGCCGTTGATCGCCATCAGTAGGAGCAAGGTATCGAGCAATGc         <  1:232415/140‑1 (MQ=255)
                                                                                                                                      gTTCGAATGCCCCCTTTTTGCAATAATGTTTGCGCCTCGCGCTGTAGCCAGCGCAGGCCCCCCTCGCCACCGCGCGATGCAACGTTGGTATCGCCGTTGATCGCCATCAGTAGGAGCAAGGTATCGAGCAATGCCAGTTc   <  1:573190/140‑1 (MQ=255)
                                                                                                                                      gTTCGAATGCCCCCTTTTTGCAATAATGTTTGCGCCTCGCGCTGTAGCCAGCGCAGGCCCCCCTCGCCACCGCGCGATGCAACGTTGGTATCGCCGTTGATCGCCATCAGTAGGAGCAAGGTATCGAGCAATGCCAGTTc   <  1:715408/140‑1 (MQ=255)
                                                                                                                                       ttCGAATGCCCCCTTTTTGCAATAATGTTTGCGCCTCGCGCTGTAGCCAGCGCAGGCCCCCCTCGCCACCGCGCGATGCAACGTTGGTATCGCCGTTGATCGCCATCAGTAGGAGCAAGGTATCGAGCAATGCCAGTTCa  <  2:60254/140‑1 (MQ=255)
                                                                                                                                       ttCGAATGCCCCCTTTTTGCAATAATGTTTGCGCCTCGCGCTGTAGCCAGCGCAGGCCCCCCTCGCCACCGCGCGATGCAACGTTGGTATCGCCGTTGATCGCCATCAGTAGGAGCAAGGTATCGAGCAATGCCAGTTCa  <  1:417767/140‑1 (MQ=255)
                                                                                                                                       |                                                                                                                                            
TTATCAAGTGCTTAAATAATTAAATCTGTGCTAAAAACCAGGTAAGGATCAGTAGGTCAGCACTGCCGCCTGGACTGAGATTTCGTTCGATACACTCCCTGTCGAACTGCCGGAGATAATCGAGATCGGCGGGGGTTCGAATGCCCCCTTTTTGCAATAATGTTTGCGCCTCGCGCTGTAGCCAGCGCAGGCCCCCCTCGCCACCGCGCGATGCAACGTTGGTATCGCCGTTGATCGCCATCAGTAGGAGCAAGGTATCGAGCAATGCCAGTTCA  >  NC_000913/646612‑646886

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