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

AAGTGCTTAAATAATTAAATCTGTGCTAAAAACCAGGTAAGGATCAGTAGGTCAGCACTGCCGCCTGGACTGAGATTTCGTTCGATACACTCCCTGTCGAACTGCCGGAGATAATCGAGATCGGCGGGGGTTCGAATGCCCCCTTTTTGCAATAATGTTTGCGCCTCGCGCTGTAGCCAGCGCAGGCCCCCCTCGCCACCGCGCGATGCAACGTTGGTATCGCCGTTGATCGCCATCAGTAGGAGCAAGGTATCGAGC  >  NC_000913/646617‑646874
                                                                                                                                   |                                                                                                                                
aaGTGCTTAAATAATTAAATCTGTGCTAAAAACCAGGTAAGGATCAGTAGGTCAGCACTGCCGCCTGGACTGAGATTTCGTTCGATACACTCCCTGTCGAACTGCCGGAGATAATCGAGATCGGCGGGGGTTCGAATGcc                                                                                                                         <  1:923670/140‑1 (MQ=255)
     cTTAAATAATTAAATCTGTGCTAAAAACCAGGTAAGGATCAGTAGGTCAGCACTGCCGCCTGGACTGAGATTTCGTTCGATACACTCCCTGTCGAACTGCCGGAGATAATCGAGATCGGCGGGGGGTGGGATTGCCCCCt                                                                                                                     >  2:595066/1‑140 (MQ=255)
      agacagaatTAAATCTGTGCTAAAAACCAGGTAAGGATCAGTAGGTCAGCACTGCCGCCTGGACTGAGATTTCGTTCGATACACTCCCTGTCGAACTGCCGGAGATAATCGAGATCGGCGGGGGTTCGAATGCCCCCttt                                                                                                                   <  2:911091/134‑1 (MQ=255)
       gacagaatTAAATCTGTGCTAAAAACCAGGTAAGGATCAGTAGGTCAGCACTGCCGCCTGGACTGAGATTTCGTTCGATACACTCCCTGTCGAACTGCCGGAGATAATCGAGATCGGCGGGGGTTCGAATGCCCCCtttt                                                                                                                  <  1:1138071/135‑1 (MQ=255)
               tAAATCTGTGCTAAAAACCAGGTAAGGATCAGTAGGTCAGCACTGCCGCCTGGACTGAGATTTCGTTCGATACACTCCCTGTCGAACTGCCGGAGATAATCGAGATCGGCGGGGGGTCGGATGCCCCCCTTTTTCtataa                                                                                                          >  2:464841/1‑140 (MQ=255)
                  aTCTGTGCTAAAAACCAGGTAAGGATCAGTAGGTCAGCACTGCCGCCTGGACTGAGATTTCGTTCGATACACTCCCTGTCGAACTGCCGGAGATAATCGAGATCGGCGGGGGTTCGAATGCCCCCTTTTTGCAATAATGt                                                                                                       <  1:19574/140‑1 (MQ=255)
                       tgCTAAAAACCAGGTAAGGATCAGTAGGTCAGCACTGCCGCCTGGACTGAGATTTCGTTCGATACACTCCCTGTCGAACTGCCGGAGATAATCGAGATCGGCGGGGGTTCGAATGCCCCCTTTTTGCAATAATGTTTgcg                                                                                                  <  1:744793/140‑1 (MQ=255)
                       tgCTAAAAACCAGGTAAGGATCAGTAGGTCAGCACTGCCGCCTGGACTGAGATTTCGTTCGATACACTCCCTGTCGAACTGCCGGAGATAATCGAGATCGGCGGGGGGTTGGATGCCCCCCTTTTTTTATAATGtttttg                                                                                                  >  2:604898/1‑137 (MQ=255)
                       tgCTAAAAACCAGGTAAGGATCAGTAGGTCAGCACTGCCGCCTGGACTGAGATTTCGTTCGATACACTCCCTGTCGAACTGCCGGAGATAATCGAGATCGGCGGGGGGTCGGATGGCCCCCTTTTTTAATAAAGTTTtcg                                                                                                  >  2:748016/1‑140 (MQ=255)
                       tgCTAAAAACCAGGTAAGGATCAGTAGGTCAGCACTGCCGCCTGGACTGAGATTTCGTTCGATACACTCCCTGTCGAACTGCCGGAGATAATCGAGATCGGCGGGGGGTCGGATGCCCCCCTTTTTCCATAATGTCTgtt                                                                                                  >  2:1138071/1‑138 (MQ=255)
                       tgCTAAAAACCAGGTAAGGATCAGTAGGTCAGCACTGCCGCCTGGACTGAGATTTCGTTCGATACACTCCCTGTCGAACTGCCGGAGATAATCGAGATCGGCGGGGGGTCGGAAGCCCCCCTTTTTTAATAATAtttttg                                                                                                  >  1:356978/1‑137 (MQ=255)
                       tgCTAAAAACCAGGTAAGGATCAGTAGGTCAGCACTGCCGCCTGGACTGAGATTTCGTTCGATACACTCCCTGTCGAACTGCCGGAGATAATCGAGATCGGCGGGGGGTCGAATGCCCCCCTTTTTCAATAATGTTTtcg                                                                                                  >  1:713644/1‑140 (MQ=255)
                        gCTAAAAACCAGGTAAGGATCAGTAGGTCAGCACTGCCGCCTGGACTGAGATTTCGTTCGATACACTCCCTGTCGAACTGCCGGAGATAATCGAGATCGGCGGGGGGTTGGATGCCCCCCTTTTTCAATAATGttttccc                                                                                                 >  1:56626/1‑136 (MQ=255)
                         cTAAAAACCAGGTAAGGATCAGTAGGTCAGCACTGCCGCCTGGACTGAGATTTCGTTCGATACACTCCCTGTCGAACTGCCGGAGATAATCGAGATCGGCGGGGGGTCGAATGCCCCCCTTTTTCAATAATGTTTGCGcc                                                                                                >  1:64569/1‑140 (MQ=255)
                          tAAAAACCAGGTAAGGATCAGTAGGTCAGCACTGCCGCCTGGACTGAGATTTCGTTCGATACACTCCCTGTCGAACTGCCGGAGATAATCGAGATCGGCGGGGGTTCGAATGCCCCCTTTTTGCAATAATGTTTGCGCCt                                                                                               <  2:1253967/140‑1 (MQ=255)
                           aaaaaCCAGGTAAGGATCAGTAGGTCAGCACTGCCGCCTGGACTGAGATTTCGTTCGATACACTCCCTGTCGAACTGCCGGAGATAATCGAGATCGGCGGGGGTTCGAATGCCCCCTTTTTGCAATAATGTTTGCGCCTc                                                                                              <  1:604898/140‑1 (MQ=255)
                                 cAGGTAAGGATCAGTAGGTCAGCACTGCCGCCTGGACTGAGATTTCGTTCGATACACTCCCTGTCGAACTGCCGGAGATAATCGAGATCGGCGGGGGGTCGAATGCCCCCCTTTTTTAATAATTTTTTGGCCTCCGGCTg                                                                                        >  2:744793/1‑140 (MQ=255)
                                              gTAGGTCAGCACTGCCGCCTGGACTGAGATTTCGTTCGATACACTCCCTGTCGAACTGCCGGAGATAATCGAGATCGGCGGGGGTTCGAATGCCCCCTTTTTGCAATAATGTTTGCGCCTCGCGCTGTAGCCAGCGCAgg                                                                           <  2:1116419/140‑1 (MQ=255)
                                              gTAGGTCAGCACTGCCGCCTGGACTGAGATTTCGTTCGATACACTCCCTGTCGAACTGCCGGAGATAATCGAGATCGGCGGGGGTTCGAATGCCCCCTTTTTGCAATAATGTTTGCGCCTCGCGCTGTAGCCAGCGCAgg                                                                           <  1:346971/140‑1 (MQ=255)
                                              gTAGGTCAGCACTGCCGCCTGGACTGAGATTTCGTTCGATACACTCCCTGTCGAACTGCCGGAGATAATCGAGATCGGCGGGGGTTCGAATGCCCCCTTTTTGCAATAATGTTTGCGCCTCGCGCTGTAGCCAGCGCAgg                                                                           <  1:748016/140‑1 (MQ=255)
                                                aGGTCAGCACTGCCGCCTGGACTGAGATTTCGTTCGATACACTCCCTGTCGAACTGCCGGAGATAATCGAGATCGGCGGGGGTTCGAATGCCCCCTTTTTGCAATAATGTTTGCGCCTCGCGCTGTAGCCAGCGCAGGcc                                                                         <  1:768772/140‑1 (MQ=255)
                                                    cAGCACTGCCGCCTGGACTGAGATTTCGTTCGATACACTCCCTGTCGAACTGCCGGAGATAATCGAGATCGGCGGGGGGTCGAATACCCCCCTTTTTTAATAATATTTTCGCCTCCCGCTGTAGCCCGCCCACGcccccc                                                                     >  2:1193200/1‑140 (MQ=255)
                                                       cACTGCCGCCTGGACTGAGATTTCGTTCGATACACTCCCTGTCGAACTGCCGGAGATAATCGAGATCGGCGGGGGGGTGGATGTCCCCCTTTTTCCATAATTTTTTTGGCTCGCGCCGGAGGCAGCGCCGGccccccccc                                                                  >  2:102706/1‑137 (MQ=255)
                                                         cTGCCGCCTGGACTGAGATTTCGTTCGATACACTCCCTGTCGAACTGCCGGAGATAATCGAGATCGGCGGGGGGTCGAATGCCCCCCTTTTTCAATAATGTTTTCGCCTCGCGCTGTAGCCAGCGCCGGccccccccccc                                                                >  1:1032805/1‑135 (MQ=255)
                                                            ccgccTGGACTGAGATTTCGTTCGATACACTCCCTGTCGAACTGCCGGAGATAATCGAGATCGGCGGGGGGTCGGATTCCCCCCTTTTTTAATAATGTTTTCGCCTCCCGCTGTAGCCAGCGCAGGGCCCCCCCGCcccc                                                             >  2:469533/1‑140 (MQ=255)
                                                            ccgccTGGACTGAGATTTCGTTCGATACACTCCCTGTCGAACTGCCGGAGATAATCGAGATCGGCGGGGGGTCGAAAGCCCCCCTTTTTCAATAATGTTTTCGCCTCGCGCTGTAGCCAGCGCAGGCCCCCCCCCCCAcc                                                             >  1:494058/1‑140 (MQ=255)
                                                               ccTGGACTGAGATTTCGTTCGATACACTCCCTGTCGAACTGCCGGAGATAATCGAGATCGGCGGGGGTTCGAATGCCCCCTTTTTGCAATAATGTTTGCGCCTCGCGCTGTAGCCAGCGCAGGCCCCCCTCGCCACcgcg                                                          <  2:35850/140‑1 (MQ=255)
                                                                                ttcgATACACTCCCTGTCGAACTGCCGGAGATAATCGAGATCGGCGGGGGTTCGAATGCCCCCTTTTTGCAATAATGTTTGCGCCTCGCGCTGTAGCCAGCGCAGGCCCCCCTCGCCACCGCGCGATGCAACGTTGGTAt                                        <  1:595066/140‑1 (MQ=255)
                                                                                ttcgATACACTCCCTGTCGAACTGCCGGAGATAATCGAGATCGGCGGGGGGTGGGATGCCCCCCTTTTTCAATAAAATTTTTGCCCCCCGCCGTAGCCAGCGCAGGCCCCCCCCCCCCCCCCGCGATGCAACGTTGGTAt                                        >  2:87443/1‑140 (MQ=255)
                                                                                     tACACTCCCTGTCGAACTGCCGGAGATAATCGAGATCGGCGGGGGGTTGAATGCCCCCCTTTTTCAATAATTTTTTTGCCCCCCGCCGTAGGCAGCGCCGGCCCCCCCCCCCCCCCCGCGGTGCAACGTTGGGTTCGccc                                   >  2:332022/1‑139 (MQ=255)
                                                                                         cTCCCTGTCGAACTGCCGGAGATAATCGAGATCGGCGGGGGTTCGAATGCCCCCTTTTTGCAATAATGTTTGCGCCTCGCGCTGTAGCCAGCGCAGGCCCCCCTCGCCACCGCGCGATGCAACGTTGGTATCGCCGTTGa                               <  2:920318/140‑1 (MQ=255)
                                                                                         cTCCCTGTCGAACTGCCGGAGATAATCGAGATCGGCGGGGGTTCGAATGCCCCCTTTTTGCAATAATGTTTGCGCCTCGCGCTGTAGCCAGCGCAGGCCCCCCTCGCCACCGCGCGATGCAACGTTGGTATCGCCGTTGa                               <  2:356978/140‑1 (MQ=255)
                                                                                              tGTCGAACTGCCGGAGATAATCGAGATCGGCGGGGGTTCGAATGCCCCCTTTTTGCAATAATGTTTGCGCCTCGCGCTGTAGCCAGCGCAGGCCCCCCTCGCCACCGCGCGATGCAACGTTGGTATCGCCGTTGATCGcc                          <  1:1043437/140‑1 (MQ=255)
                                                                                               gTCGAACTGCCGGAGATAATCGAGATCGGCGGGGGTTCGAATGCCCCCTTTTTGCAATAATGTTTGCGCCTCGCGCTGTAGCCAGCGCAGGCCCCCCTCGCCACCGCGCGATGCAACGTTGGTATCGCCGTTGATCGCCa                         <  2:91463/140‑1 (MQ=255)
                                                                                               gTCGAACTGCCGGAGATAATCGAGATCGGCGGGGGTTCGAATGCCCCCTTTTTGCAATAATGTTTGCGCCTCGCGCTGTAGCCAGCGCAGGCCCCCCTCGCCACCGCGCGATGCAACGTTGGTATCGCCGTTGATCGCCa                         <  1:332022/140‑1 (MQ=255)
                                                                                               gTCGAACTGCCGGAGATAATCGAGATCGGCGGGGGGTTGGATGCCCCCCTTTTTTTATAATGTTTTCGCCCCGCCCCGTAGCCAGCGCAGGGCCCCCCCCCCCCCCCCCGGTGCCACGTTGGTGTTGCCGTTGGTCGccc                         >  2:137861/1‑139 (MQ=255)
                                                                                               gTCGAACTGCCGGAGATAATCGAGATCGGCGGGGGGTTGAATGCCCCCCTTTTTCAATAATGTTTTCGCCCCGCGCCGTAGCCAGCGCAGGCCCCCCCCCCCCCCCCGCGAGGCAACGCTGGTATCGCCGGTGGTccccc                         >  2:395084/1‑139 (MQ=255)
                                                                                               gTCGAACTGCCGGAGATAATCGAGATCGGCGGGGGGTTGAATGCCCCCCTTTTTCAATAAAGTTTTTGCCTCCCCCTGTAGCCAGCGCCGGCCCCCCCCCCCCCCCCCCGCGGGAACCGTTGGATCGCCGCTGATCGccc                         >  2:669802/1‑139 (MQ=255)
                                                                                               gTCGAACTGCCGGAGATAATCGAGATCGGCGGGGGGTCGCATGCCCCCCTTTTTCAATAATGTTTTCGCCCCGCGCCGTAGCCAGCGCCGGCCCCCCCCCCCCCCCCCCGGTGCAACGGTGGTATCGCCCGTGGTCGccc                         >  2:1043437/1‑139 (MQ=255)
                                                                                                tCGAACTGCCGGAGATAATCGAGATCGGCGGGGGTTCGAATGCCCCCTTTTTGCAATAATGTTTGCGCCTCGCGCTGTAGCCAGCGCAGGCCCCCCTCGCCACCGCGCGATGCAACGTTGGTATCGCCGTTGATCGCCAt                        <  2:1197711/140‑1 (MQ=255)
                                                                                                tCGAACTGCCGGAGATAATCGAGATCGGCGGGGGTTCGAATGCCCCCTTTTTGCAATAATGTTTGCGCCTCGCGCTGTAGCCAGCGCAGGCCCCCCTCGCCACCGCGCGATGCAACGTTGGTATCGCCGTTGATCGCCAt                        <  1:102706/140‑1 (MQ=255)
                                                                                                gacagcTGCCGGAGATAATCGAGATCGGCGGGGGTTCGAATGCCCCCTTTTTGCAATAATGTTTGCGCCTCGCGCTGTAGCCAGCGCAGGCCCCCCTCGCCACCGCGCGATGCAACGTTGGTATCGCCGTTGATCGCCAt                        <  1:542083/135‑1 (MQ=255)
                                                                                                   aaCTGCCGGAGATAATCGAGATCGGCGGGGGTTCGAATGCCCCCTTTTTGCAATAATGTTTGCGCCTCGCGCTGTAGCCAGCGCAGGCCCCCCTCGCCACCGCGCGATGCAACGTTGGTATCGCCGTTGATCGCCATCAg                     <  2:781973/140‑1 (MQ=255)
                                                                                                    aCTGCCGGAGATAATCGAGATCGGCGGGGGTTCGAATGCCCCCTTTTTGCAATAATGTTTGCGCCTCGCGCTGTAGCCAGCGCAGGCCCCCCTCGCCACCGCGCGATGCAACGTTGGTATCGCCGTTGATCGCCATCAGt                    <  2:64569/140‑1 (MQ=255)
                                                                                                    aCTGCCGGAGATAATCGAGATCGGCGGGGGTTCGAATGCCCCCTTTTTGCAATAATGTTTGCGCCTCGCGCTGTAGCCAGCGCAGGCCCCCCTCGCCACCGCGCGATGCAACGTTGGTATCGCCGTTGATCGCCATCAGt                    <  1:981790/140‑1 (MQ=255)
                                                                                                    aCTGCCGGAGATAATCGAGATCGGCGGGGGTTCGAATGCCCCCTTTTTGCAATAATGTTTGCGCCTCGCGCTGTAGCCAGCGCAGGCCCCCCTCGCCACCGCGCGATGCAACGTTGGTATCGCCGTTGATCGCCATCAGt                    <  1:137861/140‑1 (MQ=255)
                                                                                                        ccGGAGATAATCGAGATCGGCGGGGGGTTGGATGGCCCCCTTTTTTAATAATATTTTCGCCCCCCGCTGTAGCCCGCGCAGGCCCCCCCCCCCCCCCCCCGCGGCCAAGGTGGTATCGCCCGTGATCGCCACCAGTAggg                >  2:128681/1‑139 (MQ=255)
                                                                                                        ccGGAGATAATCGAGATCGGCGGGGGGTCGGATGGCCCCCTTTTTTAATAATTTTTTCGCCCCGCGCTGTAGCCAGCGCAGGCCCCCCCCCCCCCCCCGCGATGCCACGTTGGTATCGCCGGTGATTGCCATCaggaggg                >  2:1174097/1‑139 (MQ=255)
                                                                                                           gagaTAATCGAGATCGGCGGGGGTTCGAATGCCCCCTTTTTGCAATAATGTTTGCGCCTCGCGCTGTAGCCAGCGCAGGCCCCCCTCGCCACCGCGCGATGCAACGTTGGTATCGCCGTTGATCGCCATCAGTAGGAGCa             <  1:1043758/140‑1 (MQ=255)
                                                                                                           gagaTAATCGAGATCGGCGGGGGTTCGAATGCCCCCTTTTTGCAATAATGTTTGCGCCTCGCGCTGTAGCCAGCGCAGGCCCCCCTCACCACCGCGCGATGCAACGTTGGTATCGCCGTTGATCGCCATCAGTAGGAGCa             <  1:930051/140‑1 (MQ=255)
                                                                                                                aaTCGAGATCGGCGGGGGTTCGAATGCCCCCTTTTTGCAATAATGTTTGCGCCTCGCGCTGTAGCCAGCGCAGGCCCCCCTCGCCACCGCGCGATGCAACGTTGGTATCGCCGTTGATCGCCATCAGTAGGAGCAAGGTa        <  1:270975/140‑1 (MQ=255)
                                                                                                                aaTCGAGATCGGCGGGGGGTCGGATGCCCCCCTTTTTCTATAATGTTTTCGCCTCCCGCCGTAGCCAGCGCCGGCCCCCCCCCCCCCCCCGCCGTGCCACGGGGTATCGCCCGTGATCGCCCTCAGTAGGAGGaaaggaa       >  2:137710/1‑135 (MQ=255)
                                                                                                                aaTCGAGATCGGCGGGGGGTCGGATGCCCCCCTTTTTCAATAATTTTTTCGCCTCGCGCCGTAGCCAGCGCAGGCCCCCCCCCCCCCCCCCCGCTGCAACGTTGGGATCGCCCCTGATCGCCATCACGAGGGGCactgtc        >  2:542083/1‑135 (MQ=255)
                                                                                                                 aTCGAGATCGGCGGGGGTTCGAATGCCCCCTTTTTGCAATAATGTTTGCGCCTCGCGCTGTAGCCAGCGCAGGCCCCCCTCGCCACCGCGCGATGCAACGTTGGTATCGCCGTTGATCGCCATCAGTAGGAGCAAGGTAt       <  2:713644/140‑1 (MQ=255)
                                                                                                                      gaTCGGCGGGGGTTCGAATGCCCCCTTTTTGCAATAATGTTTGCGCCTCGCGCTGTAGCCAGCGCAGGCCCCCCTCGCCACCGCGCGATGCAACGTTGGTATCGCCGTTGATCGCCATCAGTAGGAGCAAGGTATCGAGc  <  1:496642/140‑1 (MQ=255)
                                                                                                                      gaTCGGCGGGGGTTCGAATGCCCCCTTTTTGCAATAATGTTTGCGCCTCGCGCTGTAGCCAGCGCAGGCCCCCCTCGCCACCGCGCGATGCAACGTTGGTATCGCCGTTGATCGCCATCAGTAGGAGCAAGGTATCGAGc  <  1:869556/140‑1 (MQ=255)
                                                                                                                                   |                                                                                                                                
AAGTGCTTAAATAATTAAATCTGTGCTAAAAACCAGGTAAGGATCAGTAGGTCAGCACTGCCGCCTGGACTGAGATTTCGTTCGATACACTCCCTGTCGAACTGCCGGAGATAATCGAGATCGGCGGGGGTTCGAATGCCCCCTTTTTGCAATAATGTTTGCGCCTCGCGCTGTAGCCAGCGCAGGCCCCCCTCGCCACCGCGCGATGCAACGTTGGTATCGCCGTTGATCGCCATCAGTAGGAGCAAGGTATCGAGC  >  NC_000913/646617‑646874

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: