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

AAATTTATCAAGTGCTTAAATAATTAAATCTGTGCTAAAAACCAGGTAAGGATCAGTAGGTCAGCACTGCCGCCTGGACTGAGATTTCGTTCGATACACTCCCTGTCGAACTGCCGGAGATAATCGAGATCGGCGGGGGTTCGAATGCCCCCTTTTTGCAATAATGTTTGCGCCTCGCGCTGTAGCCAGCGCAGGCCCCCCTCGCCACCGCGCGATGCAACGTTGGTATCGCCGTTGATCGCCATCAGTAGGAGCAAGGTATCGAGCAATGCCAGTTCA  >  NC_000913/646608‑646886
                                                                                                                                           |                                                                                                                                            
aaaTTTATCAAGTGCTTAAATAATTAAATCTGTGCTAAAAACCAGGTAAGGATCAGTAGGTCAGCACTGCCGCCTGGACTGAGATTTCGTTCGATACACTCCCTGTCGAACTGCCGGAGATAATCGAGATCGGCGGGGGt                                                                                                                                              <  2:789069/140‑1 (MQ=255)
 aaTTTATCAAGTGCTTAAATAATTAAATCTGTGCTAAAAACCAGGTAAGGATCAGTAGGTCAGCACTGCCGCCTGGACTGAGATTTCGTTCGATACACTCCCTGTCGAACTGCCGGAGATAATCGAGATCGGCGGGGGtt                                                                                                                                             <  1:835901/140‑1 (MQ=255)
      aTCAAGTGCTTAAATAATTAAATCTGTGCTAAAAACCAGGTAAGGATCAGTAGGTCAGCACTGCCGCCTGGACTGAGATTTCGTTCGATACACTCCCTGTCGAACTGCCGGAGATAATCGAGATCGGCGGGGGTTCGatt                                                                                                                                        >  1:668219/1‑138 (MQ=255)
           gTGCTTAAATAATTAAATCTGTGCTAAAAACCAGGTAAGGATCAGTAGGTCAGCACTGCCGCCTGGACTGAGATTTCGTTCGATACACTCCCTGTCGAACTGCCGGAGATAATCGAGATCGGCGGGGGGTGCAGTGcccc                                                                                                                                  >  2:515921/1‑140 (MQ=255)
              cTTAAATAATTAAATCTGTGCTAAAAACCAGGTAAGGATCAGTAGGTCAGCACTGCCGCCTGGACTGAGATTTCGTTCGATACACTCCCTGTCGAACTGCCGGAGATAATCGAGATCGGCGGGGGGTTGAATTtccccct                                                                                                                               >  1:1504952/1‑138 (MQ=255)
              cTTAAATAATTAAATCTGTGCTAAAAACCAGGTAAGGATCAGTAGGTCAGCACTGCCGCCTGGACTGAGATTTCGTTCGATACACTCCCTGTCGAACTGCCGGAGATAATCGAGATCGGCGGGGGGTCGGATGcccccct                                                                                                                               >  2:1261629/1‑138 (MQ=255)
              cTTAAATAATTAAATCTGTGCTAAAAACCAGGTAAGGATCAGTAGGTCAGCACTGCCGCCTGGACTGAGATTTCGTTCGATACACTCCCTGTCGAACTGCCGGAGATAATCGAGATCGGCGGGGGGTCGAATGcccccct                                                                                                                               >  1:394824/1‑138 (MQ=255)
                      atTAAATCTGTGCTAAAAACCAGGTAAGGATCAGTAGGTCAGCACTGCCGCCTGGACTGAGATTTCGTTCGATACACTCCCTGTCGAACTGCCGGAGATAATCGAGATCGGCGGGGGGTCGAATGCCCCCTTTTTTCaat                                                                                                                       >  1:798874/1‑140 (MQ=255)
                      atTAAATCTGTGCTAAAAACCAGGTAAGGATCAGTAGGTCAGCACTGCCGCCTGGACTGAGATTTCGTTCGATACACTCCCTGTCGAACTGCCGGAGATAATCGAGATCGGCGGGGGGTCGAATGCCCCCCTTTTTCaaa                                                                                                                       >  1:1430972/1‑139 (MQ=255)
                      atTAAATCTGTGCTAAAAACCAGGTAAGGATCAGTAGGTCAGCACTGCCGCCTGGACTGAGATTTCGTTCGATACACTCCCTGTCGAACTGCCGGAGATAATCGAGATCGGCGGGGGGGTGAAAACCCCCcttttttaaa                                                                                                                       >  1:494748/1‑135 (MQ=255)
                       ttAAATCTGTGCTAAAAACCAGGTAAGGATCAGTAGGTCAGCACTGCCGCCTGGACTGAGATTTCGTTCGATACACTCCCTGTCGAACTGCCGGAGATAATCGAGATCGGCGGGGGGTCGAATGCCCCCCTTTTTCaata                                                                                                                      >  1:1142759/1‑140 (MQ=255)
                        tAAATCTGTGCTAAAAACCAGGTAAGGATCAGTAGGTCAGCACTGCCGCCTGGACTGAGATTTCGTTCGATACACTCCCTGTCGAACTGCCGGAGATAATCGAGATCGGCGGGGGGTTGGATGGCCCCCCTTTTCaataa                                                                                                                     >  2:1266845/1‑140 (MQ=255)
                        tAAATCTGTGCTAAAAACCAGGTAAGGATCAGTAGGTCAGCACTGCCGCCTGGACTGAGATTTCGTTCGATACACTCCCTGTCGAACTGCCGGAGATAATCGAGATCGGCGGGGGGGCGAATGCCCCCCTTTTTCaataa                                                                                                                     >  1:671651/1‑140 (MQ=255)
                               gtgCTAAAAACCAGGTAAGGACCAGTGGGTCAGCACTGCCGCCTGGACTGAGATTTCGTTCGATACACTCCCTGTCGAACTGCCGGAGATAATCGAGATCGGCGGGGGTTCGAATGCCCCCTTTTTGCAATAATGTTTgc                                                                                                              <  2:1240299/140‑1 (MQ=255)
                                tgCTAAAAACCAGGTAAGGATCAGTAGGTCAGCACTGCCGCCTGGACTGAGATTTCGTTCGATACACTCCCTGTCGAACTGCCGGAGATAATCGAGATCGGCGGGGGGTCGGATAGCCCCCTTTTTTCATAATGTTTtcg                                                                                                             >  2:919434/1‑140 (MQ=255)
                                tgCTAAAAACCAGGTAAGGATCAGTAGGTCAGCACTGCCGCCTGGACTGAGATTTCGTTCGATACACTCCCTGTCGAACTGCCGGAGATAATCGAGATCGGCGGGGGGTCGGATGGCCCCCTTTTTCAATAATGttttcc                                                                                                             >  2:878299/1‑137 (MQ=255)
                                tgCTAAAAACCAGGTAAGGATCAGTAGGTCAGCACTGCCGCCTGGACTGAGATTTCGTTCGATACACTCCCTGTCGAACTGCCGGAGATAATCGAGATCGGCGGGGGGTCGGATGCCCCCCTTTTTCAATAAAGTTTtcg                                                                                                             >  1:1431599/1‑140 (MQ=255)
                                 gCTAAAAACCAGGTAAGGATCAGTAGGTCAGCACTGCCGCCTGGACTGAGATTTCGTTCGATACACTCCCTGTCGAACTGCCGGAGATAATCGAGATCGGCGGGGGGGGGGGTGGCCCCCTTTTTTAATAATGTTTtcgc                                                                                                            >  2:903098/1‑140 (MQ=255)
                                 gCTAAAAACCAGGTAAGGATCAGTAGGTCAGCACTGCCGCCTGGACTGAGATTTCGTTCGATACACTCCCTGTCGAACTGCCGGAGATAATCGAGATCGGCGGGGGGGCGGATTTCCCCCTTTTTCTATAATAtttttgc                                                                                                            >  2:422180/1‑136 (MQ=255)
                                    aaaaaCCAGGTAAGGATCAGTAGGTCAGCACTGCCGCCTGGACTGAGATTTCGTTCGATACACTCCCTGTCGAACTGCCGGAGATAATCGAGATCGGCGGGGGTTCGAATGCCCCCTTTTTGCAATAATGTTTGCGCCTc                                                                                                         <  1:210518/140‑1 (MQ=255)
                                         ccAGGTAAGGATCAGTAGGTCAGCACTGCCGCCTGGACTGAGATTTCGTTCGATACACTCCCTGTCGAACTGCCGGAGATAATCGAGATCGGCGGGGGTTCGAATGCCCCCTTTTTGCAATAATGTTTGCGCCTCGCGCt                                                                                                    <  1:721757/140‑1 (MQ=255)
                                                 ggATCAGTAGGTCAGCACTGCCGCCTGGACTGAGATTTCGTTCGATACACTCCCTGTCGAACTGCCGGAGATAATCGAGATCGGCGGGGGCTGGAGAGGCCCCCTTTTGTAAGAAAGTTGTCGCGCCGCGCTGTGGCCac                                                                                            >  2:272569/1‑139 (MQ=255)
                                                      aGTAGGTCAGCACTGCCGCCTGGACTGAGATTTCGTTCGATACACTCCCTGTCGAACTGCCGGAGATAATCGAGATCGGCGGGGGGTCGAATGCCCCCCTTTTTCAATAATGTTTTCGCCTCGCGCCGTAGCCAGCGCAg                                                                                       >  1:851478/1‑140 (MQ=255)
                                                       gTAGGTCAGCACTGCCGCCTGGACTGAGATTTCGTTCGATACACTCCCTGTCGAACTGCCGGAGATAATCGAGATCGGCGGGGGGTGGAATGCCCCCCTTTTTCAATAAAGTTTTTGCCTCCCCCCGGAGCCCGCGCCgg                                                                                      >  2:349237/1‑140 (MQ=255)
                                                       gTAGGTCAGCACTGCCGCCTGGACTGAGATTTCGTTCGATACACTCCCTGTCGAACTGCCGGAGATAATCGAGATCGGCGGGGGGTCGGATGCCCCCCTTTTTTAATAATGTTTTCGCCTCCCCCCGGAGCCAGCGCCgg                                                                                      >  2:1169858/1‑140 (MQ=255)
                                                       gTAGGTCAGCACTGCCGCCTGGACTGAGATTTCGTTCGATACACTCCCTGTCGAACTGCCGGAGATAATCGAGATCGGCGGGGGGTCGGAAGCCCCCCTTTTTCAAAAAAATTTTCGCCCCCCCCCGTCGCCACCGCCgg                                                                                      >  2:1027443/1‑140 (MQ=255)
                                                         aGGTCAGCACTGCCGCCTGGACTGAGATTTCGTTCGATACACTCCCTGTCGAACTGCCGGAGATAATCGAGATCGGCGGGGGTTCGAATGCCCCCTTTTTGCAATAATTTTTGCGCCTCGCGCTGTAGCCAGCGCAGGcc                                                                                    <  1:1342552/140‑1 (MQ=255)
                                                             cAGCACTGCCGCCTGGACTGAGATTTCGTTCGATACACTCCCTGTCGAACTGCCGGAGATAATCGAGATCGGCGGGGGGTTGGATGCCCCCCTTTTTTAATAATGTTTTCCCCTCGCGCCGTAGCCCGCGCCGGcccccc                                                                                >  1:694299/1‑140 (MQ=255)
                                                             cAGCACTGCCGCCTGGACTGAGATTTCGTTCGATACACTCCCTGTCGAACTGCCGGAGATAATCGAGATCGGCGGGGGGTCGAATGCCCCCTTTTTTCAATAAAGTTTTCGCCTCGCGCTGTAGCCAGCGCAGGcccccc                                                                                >  1:676275/1‑140 (MQ=255)
                                                               gCACTGCCGCCTGGACTGAGATTTCGTTCGATACACTCCCTGTCGAACTGCCGGAGATAATCGAGATCGGCGGGGGGTCGGATGCCCCCCTTTTTCAATAATGTTTTCGCCCCGCGCTGTAGCCAGCGCAGGcccccccc                                                                              >  1:1504008/1‑138 (MQ=255)
                                                                cACTGCCGCCTGGACTGAGATTTCGTTCGATACACTCCCTGTCGAACTGCCGGAGATAATCGAGATCGGCGGGGGGTCGAATGCCCCCCTTTTTCAATAAAATTTTCGCCCCCCGCCGTAGCCAGCGCCGGccccccccc                                                                             >  2:1178864/1‑137 (MQ=255)
                                                                  cTGCCGCCTGGACTGAGATTTCGTTCGATACACTCCCTGTCGAACTGCCGGAGATAATCGAGATCGGCGGGGGGTGGAGTACCCCCCCTTTTTAATAAAGTTTTCGCACCGCGCTGTGCCCCCCGCCAGccccccccccc                                                                           >  2:865664/1‑135 (MQ=255)
                                                                     ccgccTGGACTGAGATTTCGTTCGATACACTCCCTGTCGAACTGCCGGAGATAATCGAGATCGGCGGGGGGTCGGATGCCCCCCTTTTTCAATAATGTTTGCGCCCCGCGCCGTAGCCCGCGCAGGGCCCCCCCCCCAcc                                                                        >  1:704993/1‑140 (MQ=255)
                                                                        ccTGGACTGAGATTTCGTTCGATACACTCCCTGTCGAACTGCCGGAGATAATCGAGATCGGCGGGGGTTCGAATGCCCCCTTTTTGCAATAATGTTTGCGCCTCGCGCTGTAGCCAGCGCAGGCCCCCCTCGCCACcgcg                                                                     <  1:29023/140‑1 (MQ=255)
                                                                            gACTGAGATTTCGTTCGATACACTCCCTGTCGAACTGCCGGAGATAATCGAGATCGGCGGGGGTTCGAATGCCCCCTTTTTGCAATAATGTTTGCGCCTCGCGCTGTAGCCAGCGCAGGCCCCCCTCGCCACCGCGCGAt                                                                 <  1:422180/140‑1 (MQ=255)
                                                                             aCTGAGATTTCGTTCGATACACTCCCTGTCGAACTGCCGGAGATAATCGAGATCGGCGGGGGTTCGAATGCCCCCTTTTTGCAATAATGTTTGCGCCTCGCGCTGTAGCCAGCGCAGGCCCCCCTCGCCACCGCGCGATg                                                                <  1:256502/140‑1 (MQ=255)
                                                                              cTGAGATTTCGTTCGATACACTCCCTGTCGAACTGCCGGAGATAATCGAGATCGGCGGGGGTTCGAATGCCCCCTTTTTGCAATAATGTTTGCGCCTCGCGCTGTAGCCAGCGCAGGCCCCCCTCGCCACCGCGCGATGc                                                               <  1:515921/140‑1 (MQ=255)
                                                                              cTGAGATTTCGTTCGATACACTCCCTGTCGAACTGCCGGAGATAATCGAGATCGGCGGGGGTTCGAATGCCCCCTTTTTGCAATAATGTTTGCGCCTCGCGCTGTAGCCAGCGCAGGCCCCCCTCGCCACCGCGCGATGc                                                               <  1:1474386/140‑1 (MQ=255)
                                                                                 agaTTTCGTTCGATACACTCCCTGTCGAACTGCCGGAGATAATCGAGATCGGCGGGGGGTCGGATGTCCCCCTTTTTCAAAAAAGTTTTCGTCTCGGGCCGTTGCCGGCGGCCGCCCCCCCCGCCCCCGCGCGGTggtac                                                            >  2:1359310/1‑136 (MQ=255)
                                                                                           cgATACACTCCCTGTCGAACTGCCGGAGATAATCGAGATCGGCGGGGGGGCGGATGCCCCCCTTTTTCAATAAATTTTTCGCCCCCCGCCGTAGCCCGCGCCGGCCCCCCCCCCCCCCCCCCGGCGCAACGGTGGGATCg                                                  >  2:1111721/1‑140 (MQ=255)
                                                                                                  cTCCCTGTCGAACTGCCGGAGATAATCGAGATCGGCGGGGGTTCGAATGCCCCCTTTTTGCAATAATGTTTGCGCCTCGCGCTGTAGCCAGCGCAGGCCCCCCTCGCCACCGCGCGATGCAACGTTGGTATCGCCGTTGa                                           <  1:1096263/140‑1 (MQ=255)
                                                                                                   tCCCTGTCGAACTGCCGGAGATAATCGAGATCGGCGGGGGTTCGAATGCCCCCTTTTTGCAATAATGTTTGCGCCTCGCGCTGTAGCCAGCGCAGGCCCCCCTCGCCACCGCGCGATGCAACGTTGGTATCGCCGTTGAt                                          <  2:1335388/140‑1 (MQ=255)
                                                                                                    cccTGTCGAACTGCCGGAGATAATCGAGATCGGCGGGGGTTCGAATGCCCCCTTTTTGCAATAATGTTTGCGCCTCGCGCTGTAGCCAGCGCAGGCCCCCCTCGCCACCGCGCGATGCAACGTTGGTATCGCCGTTGATc                                         <  1:355555/140‑1 (MQ=255)
                                                                                                     ccTGTCGAACTGCCGGAGATAATCGAGATCGGCGGGGGGTCGGATGCCCCCCTTTTTCTATAATGTTTTTGCCTCCCCCCGTAGCCAGCGCCCGCCCCCCCCCCCCCCCCCCCGTGCAACGGTGGTATCGCCGCTGATCg                                        >  2:1304763/1‑140 (MQ=255)
                                                                                                        gTCGAACTGCCGGAGATAATCGAGATCGGCGGGGGTTCGAATGCCCCCTTTTTGCAATAATGTTTGCGCCTCGCGCTGTAGCCAGCGCAGGCCCCCCTCGCCACCGCGCGATGCAACGTTGGTATCGCCGTTGATCGCCa                                     <  1:406742/140‑1 (MQ=255)
                                                                                                        gTCGAACTGCCGGAGATAATCGAGATCGGCGGGGGTTCGAATGCCCCCTTTTTGCAATAATGTTTGCGCCTCGCGCTGTAGCCAGCGCAGGCCCCCCTCGCCACCGCGCGATGCAACGTTGGTATCGCCGTTGATCGCCa                                     <  1:1485030/140‑1 (MQ=255)
                                                                                                        gTCGAACTGCCGGAGATAATCGAGATCGGCGGGGGTTCGAATGCCCCCTTTTTGCAATAATGTTTGCGCCTCGCGCTGTAGCCAGCGCAGGCCCCCCTCGCCACCGCGCGATGCAACGTTGGTATCGCCGTTGATCGCCa                                     <  1:1169858/140‑1 (MQ=255)
                                                                                                        gTCGAACTGCCGGAGATAATCGAGATCGGCGGGGGGTGGGATGTCCCCCTTTTTCAATAATTTTTTCGCCTCGCGCTGTAGCCAGCGCAGGGCCCCCCCCCCCCCCCGCGGGGCAACGTTGGGGTCGCCCTTGTACGccc                                     >  2:238392/1‑139 (MQ=255)
                                                                                                        gTCGAACTGCCGGAGATAATCGAGATCGGCGGGGGGTCGAATGCCCCCCCTTTTTCATAAAGTTTTCGGCCCCCGCTGTGGACCACGCCGGGCCCCCCCCCCCCCCCCCGCAGCCACGGTCGGGTCGCCGGTGAGCGccc                                     >  2:1034277/1‑139 (MQ=255)
                                                                                                         tCGAACTGCCGGAGATAATCGAGATCGGCGGGGGTTCGAATGCCCCCTTTTTGCAATAATGTTTGCGCCTCGCGCTGTAGCCAGCGCAGGCCCCCCTCGCCACCGCGCGATGCAACGTTGGTATCGCCGTTGATCGCCAt                                    <  1:2746/140‑1 (MQ=255)
                                                                                                             aCTGCCGGAGATAATCGAGATCGGCGGGGGTTCGAATGCCCCCTTTTTGCAATAATGTTTGCGCCTCGCGCTGTAGCCAGCGCAGGCCCCCCTCGCCACCGCGCGATGCAACGTTGGTATCGCCGTTGATCGCCATCAGt                                <  2:500105/140‑1 (MQ=255)
                                                                                                             aCTGCCGGAGATAATCGAGATCGGCGGGGGGTCGAAAGCCCCCCTTTTTCCATAATGTTTTCGCCCCCCGCCGGTGCCCGCGCCGGCCCCCCCCCCCCCCCCGCGCCGCAACGTTGGGATCGCCGCTGTTCGCCCTCCGt                                >  2:1053267/1‑140 (MQ=255)
                                                                                                             aCTGCCGGAGATAATCGAGATCGGCGGGGGGGTGGATGGCCCCCTTTTTCAATAATATTTTGGCCCCCCCCCGTAGCCCGCGCCGGCCCCCCCCCCCCCCCCCCGGTGCAACGGTGGGGTCGCCGGTGTTCGccccccgg                                >  2:1170928/1‑134 (MQ=255)
                                                                                                                gCCGGAGATAATCGAGATCGGCGGGGGTTCGAATGCCCCCTTTTTGCAATAATGTTTGCGCCTCGCGCTGTAGCCAGCGCAGGCCCCCCTCGCCACCGCGCGATGCAACGTTGGTATCGCCGTTGATCGCCATCAGTAgg                             <  2:851478/140‑1 (MQ=255)
                                                                                                                    gagaTAATCGAGATCGGCGGGGGTTCGAATGCCCCCTTTTTGCAATAATGTTTGCGCCTCGCGCTGTAGCCAGCGCAGGCCCCCCTCGCCACCGCGCGATGCAACGTTGGTATCGCCGTTGATCGCCATCAGTAGGAGCa                         <  1:1266845/140‑1 (MQ=255)
                                                                                                                    gagaTAATCGAGATCGGCGGGGGTTCGAATGCCCCCTTTTTGCAATAATGTTTGCGCCTCGCGCTGTAGCCAGCGCAGGCCCCCCTCGCCACCGCGCGATGCAACGTTGGTATCGCCGTTGATCGCCATCAGTAGGAGCa                         <  1:349237/140‑1 (MQ=255)
                                                                                                                      gaTAATCGAGATCGGCGGGGGTTCGAATGCCCCCTTTTTGCAATAATGTTTGCGCCTCGCGCTGTAGCCAGCGCAGGCCCCCCTCGCCACCGCGCGATGCAACGTTGGTATCGCCGTTGATCGCCATCAGTAGGAGCAAg                       <  1:528965/140‑1 (MQ=255)
                                                                                                                          aTCGAGATCGGCGGGGGTTCGAATGCCCCCTTTTTGCAATAATGTTTGCGCCTCGCGCTGTAGCCAGCGCAGGCCCCCCTCGCCACCGCGCGATGCAACGTTGGTATCGCCGTTGATCGCCATCAGTAGGAGCAAGGTAt                   <  2:1286062/140‑1 (MQ=255)
                                                                                                                          aTCGAGATCGGCGGGGGTTCGAATGCCCCCTTTTTGCAATAATGTTTGCGCCTCGCGCTGTAGCCAGCGCAGGCCCCCCTCGCCACCGCGCGATGCAACGTTGGTATCGCCGTTGATCGCCATCAGTAGGAGCAAGGTAt                   <  1:1508432/140‑1 (MQ=255)
                                                                                                                               gaTCGGCGGGGGTTCGAATGCCCCCTTTTTGCAATAATGTTTGCGCCTCGCGCTGTAGCCAGCGCAGGCCCCCCTCGCCACCGCGCGATGCAACGTTGGTATCGCCGTTGATCGCCATCAGTAGGAGCAAGGTATCGAGc              <  1:435755/140‑1 (MQ=255)
                                                                                                                                    gcggGGGTTCGAATGCCCCCTTTTTGCAATAATGTTTGCGCCTCGCGCTGTAGCCAGCGCAGGCCCCCCTCGCCACCGCGCGATGCAACGTTGGTATCGCCGTTGATCGCCATCAGTAGGAGCAAGGTATCGAGCAATGc         <  1:520133/140‑1 (MQ=255)
                                                                                                                                    gcggGGGTTCGAATGCCCCCTTTTTGCAATAATGTTTGCGCCTCGCGCTGTAGCCAGCGCAGGCCCCCCTCGCCACCGCGCGATGCAACGTTGGTATCGCCGTTGATCGCCATCAGTAGGAGCAAGGTATCGAGCAATGc         <  1:1225891/140‑1 (MQ=255)
                                                                                                                                    gcggGGGTTCGAATGCCCCCTTTTTGCAATAATGTTTGCGCCTCGCGCTGTAGCCAGCGCAGGCCCCCCTCGCCACCGCGCGATGCAACGTTGGTATCGCCGTTGATCGCCATCAGTAGGAGCAAGGTATCGAGCAATGc         <  2:121773/140‑1 (MQ=255)
                                                                                                                                    gcggGGGTTCGAATGCCCCCTTTTTGCAATAATGTTTGCGCCTCGCGCTGTAGCCAGCGCAGGCCCCCCTCGCCACCGCGCGATGCAACGTTGGTATCGCCGTTGATCGCCATCAGTAGGAGCAAGGTATCGAGCAATGc         <  1:1304763/140‑1 (MQ=255)
                                                                                                                                    gcggGGGTTCGAATGCCCCCTTTTTGCAATAATGTTTGCGCCTCGCGCTGTAGCCAGCGCAGGCCCCCCTCGCCACCGCGCGATGCAACGTTGGTATCGCCGTTGATCGCCATCAGTAGGAGCAAGGTATCGAGCAATGc         <  2:1504008/140‑1 (MQ=255)
                                                                                                                                    gcggGGGTTCGAATGCCCCCTTTTTGCAATAATGTTTGCGCCTCGCGCTGTAGCCAGCGCAGGCCCCCCTCGCCACCGCGCGATGCAACGTTGGTATCGCCGTTGATCGCCATCAGTAGGAGCAAGGTATCGAGCAATGc         <  1:674102/140‑1 (MQ=255)
                                                                                                                                    gcggGGGTTCGAATGCCCCCTTTTTGCAATAATGTTTGCGCCTCGCGCTGTAGCCAGCGCAGGCCCCCCTCGCCACCGCGCGATGCAACGTTGGTATCGCCGTTGATCGCCATCAGTAGGAGCAAGGTATCGAGCAATGc         <  1:878299/140‑1 (MQ=255)
                                                                                                                                     cggGTGATCGAATGCTCCCTTTTTGCAATTATGGTTGCGCCTCCCGCTGTAGCCAGCGCAGGCCTCCCTCGTCACCGCGCGTTGCAACGTTGGTATCGCCGTTGATCGCCATCAGTAGGAGCAAGGTATCGAGCAATGcc        <  2:668219/140‑1 (MQ=255)
                                                                                                                                          gTTCGAATGCCCCCTTTTTGCAATAATGTTTGCGCCTCGCGCTGTAGCCAGCGCAGGCCCCCCTCGCCACCGCGCGATGCAACGTTGGTATCGCCGTTGATCGCCATCAGTAGGAGCAAGGTATCGAGCAATGCCAGTTc   <  2:394824/140‑1 (MQ=255)
                                                                                                                                           ttCGAATGCCCCCTTTTTGCAATAATGTTTGCGCCTCGCGCTGTAGCCAGCGCAGGCCCCCCTCGCCACCGCGCGATGCAACGTTGGTATCGCCGTTGATCGCCATCAGTAGGAGCAAGGTATCGAGCAATGCCAGTTCa  <  1:599695/140‑1 (MQ=255)
                                                                                                                                           ttCGAATGCCCCCTTTTTGCAATAATGTTTGCGCCTCGCGCTGTAGCCAGCGCAGGCCCCCCTCGCCACCGCGCGATGCAACGTTGGTATCGCCGTTGATCGCCATCAGTAGGAGCAAGGTATCGAGCAATGCCAGTTCa  <  2:1431599/140‑1 (MQ=255)
                                                                                                                                           |                                                                                                                                            
AAATTTATCAAGTGCTTAAATAATTAAATCTGTGCTAAAAACCAGGTAAGGATCAGTAGGTCAGCACTGCCGCCTGGACTGAGATTTCGTTCGATACACTCCCTGTCGAACTGCCGGAGATAATCGAGATCGGCGGGGGTTCGAATGCCCCCTTTTTGCAATAATGTTTGCGCCTCGCGCTGTAGCCAGCGCAGGCCCCCCTCGCCACCGCGCGATGCAACGTTGGTATCGCCGTTGATCGCCATCAGTAGGAGCAAGGTATCGAGCAATGCCAGTTCA  >  NC_000913/646608‑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: