Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_000913646,7470TG40.8% 118.3 / 11.9 72T255P (ACC→CCC) citG2‑(5''‑triphosphoribosyl)‑3'‑dephosphocoenzyme‑A synthase
Reads supporting (aligned to +/- strand):  ref base T (9/33);  new base G (29/0);  total (39/33)
Fisher's exact test for biased strand distribution p-value = 2.51e-12
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

ATTTCCAAATTTATCAAGTGCTTAAATAATTAAATCTGTGCTAAAAACCAGGTAAGGATCAGTAGGTCAGCACTGCCGCCTGGACTGAGATTTCGTTCGATACACTCCCTGTCGAACTGCCGGAGATAATCGAGATCGGC‑GGGGGTTCG‑AAT‑G‑CCCCCTTTTTGCAATAATGTTTGCGCCTCGCGCTGTAGCCAGCG‑CAGGCCCCCCTCGCCACCGCGCGATGCAACGTTGGTATCGCCGTTGATCGCCATCAGTAGGAGCAAGGTATCGAGCAATGCCAGTTCAGGATCTAACC  >  NC_000913/646602‑646896
                                                                                                                                                  |                                                                                                                                                         
atTTCCAAATTTATCAAGTGCTTAAATAATTAAATCTGTGCTAAAAACCAGGTAAGGATCAGTAGGTCAGCACTGCCGCCTGGACTGAGATTTCGTTCGATACACTCCCTGTCGAACTGCCGGAGATAATCGAGATCGGC‑GGGGGTTCG‑a                                                                                                                                                      <  2:260860/150‑1 (MQ=255)
    ccAAATTTATCAAGTGCTTAAATAATTAAATCTGTGCTAAAAACCAGGTAAGGATCAGTAGGTCAGCACTGCCGCCTGGACTGAGATTTCGTTCGATACACTCCCTGTCGAACTGCCGGAGATAATCGAGATCGGC‑GGGGGTTCG‑AAT‑G‑c                                                                                                                                                <  2:812030/150‑1 (MQ=255)
      aaaTTTATCAAGTGCTTAAATAATTAAATCTGTGCTAAAAACCAGGTAAGGATCAGTAGGTCAGCACTGCCGCCTGGACTGAGATTTCGTTCGATACACTCCCTGTCGAACTGCCGGAGATAATCGAGATCGGC‑GGGGGGTCG‑AAA‑G‑cccc                                                                                                                                             >  2:234315/1‑151 (MQ=255)
       aaTTTATCAAGTGCTTAAATAATTAAATCTGTGCTAAAAACCAGGTAAGGATCAGTAGGTCAGCACTGCCGCCTGGACTGAGATTTCGTTCGATACACTCCCTGTCGAACTGCCGGAGATAATCGAGATCGGC‑GGGGGGTCG‑GAA‑G‑acccc                                                                                                                                            >  1:1278118/1‑151 (MQ=255)
         tttATCAAGTGCTTAAATAATTAAATCTGTGCTAAAAACCAGGTAAGGATCAGTAGGTCAGCACTGCCGCCTGGACTGAGATTTCGTTCGATACACTCCCTGTCGAACTGCCGGAGATAATCGAGATCGGC‑GGGGGGTTG‑AAA‑G‑cccccct                                                                                                                                          >  2:774658/1‑149 (MQ=255)
          ttATCAAGTGCTTAAATAATTAAATCTGTGCTAAAAACCAGGTAAGGATCAGTAGGTCAGCACTGCCGCCTGGACTGAGATTTCGTTCGATACACTCCCTGTCGAACTGCCGGAGATAATCGAGATCGGC‑GGGGGGTCG‑AAT‑G‑acccc                                                                                                                                            >  1:8545/1‑148 (MQ=255)
               aaGTGCTTAAATAATTAAATCTGTGCTAAAAACCAGGTAAGGATCAGTAGGTCAGCACTGCCGCCTGGACTGAGATTTCGTTCGATACACTCCCTGTCGAACTGCCGGAGATAATCGAGATCGGC‑GGGGGTTCG‑AAT‑G‑CCCCCTTTTTGCa                                                                                                                                    <  2:595871/151‑1 (MQ=255)
                  tGCTTAAATAATTAAATCTGTGCTAAAAACCAGGTAAGGATCAGTAGGTCAGCACTGCCGCCTGGACTGAGATTTCGTTCGATACACTCCCTGTCGAACTGCCGGAGATAATCGAGATCGGC‑GGGGGGTCG‑AAT‑G‑CCCCCctttttct                                                                                                                                    >  1:1048737/1‑145 (MQ=255)
                    cTTAAATAATTAAATCTGTGCTAAAAACCAGGTAAGGATCAGTAGGTCAGCACTGCCGCCTGGACTGAGATTTCGTTCGATACACTCCCTGTCGAACTGCCGGAGATAATCGAGATCGGC‑GGGGGGTCG‑GAT‑G‑GCCCCCTTTTTtaataa                                                                                                                                >  1:1268465/1‑150 (MQ=255)
                     ttAAATAATTAAATCTGTGCTAAAAACCAGGTAAGGATCAGTAGGTCAGCACTGCCGCCTGGACTGAGATTTCGTTCGATACACTCCCTGTCGAACTGCCGGAGATAATCGAGATCGGC‑GGGGGTTTG‑AAT‑G‑CCCCCTTTTTGCaataa                                                                                                                                >  2:510558/1‑149 (MQ=255)
                     ttAAATAATTAAATCTGTGCTAAAAACCAGGTAAGGATCAGTAGGTCAGCACTGCCGCCTGGACTGAGATTTCGTTCGATACACTCCCTGTCGAACTGCCGGAGATAATCGAGATCGGC‑GGGGGTGCT‑AAA‑A‑CCCCCtctttttaa                                                                                                                                   >  2:665437/1‑142 (MQ=255)
                     ttAAATAATTAAATCTGTGCTAAAAACCAGGTAAGGATCAGTAGGTCAGCACTGCCGCCTGGACTGAGATTTCGTTCGATACACTCCCTGTCGAACTGCCGGAGATAATCGAGATCGGC‑GGGGGGTTG‑GAA‑G‑CCCCCCTTTTTCtataa                                                                                                                                >  2:717870/1‑149 (MQ=255)
                      tAAATAATTAAATCTGTGCTAAAAACCAGGTAAGGATCAGTAGGTCAGCACTGCCGCCTGGACTGAGATTTCGTTCGATACACTCCCTGTCGAACTGCCGGAGATAATCGAGATCGGC‑GGGGGTTCG‑AAT‑G‑CCCCCTTTTTGCaataat                                                                                                                               <  1:905094/149‑1 (MQ=255)
                        aataatTAAATCTGTGCTAAAAACCAGGTAAGGATCAGTAGGTCAGCACTGCCGCCTGGACTGAGATTTCGTTCGATACACTCCCTGTCGAACTGCCGGAGATAATCGAGATCGGC‑GGGGGTTCG‑AAT‑G‑CCCCCTTTTTGCAATAATGttt                                                                                                                           <  2:759483/151‑1 (MQ=255)
                         ataatTAAATCTGTGCTAAAAACCAGGTAAGGATCAGTAGGTCAGCACTGCCGCCTGGACTGAGATTTCGTTCGATACACTCCCTGTCGAACTGCCGGAGATAATCGAGATCGGC‑GGGGGTTCG‑AAT‑G‑CCCCCTTTTTGCAATAATGTTTg                                                                                                                          <  2:1025974/151‑1 (MQ=255)
                               aaaTCTGTGCTAAAAACCAGGTAAGGATCAGTAGGTCAGCACTGCCGCCTGGACTGAGATTTCGTTCGATACTCTCCCTGTCGAACTGCCGGAGATAATCGAGATCGGC‑GGGGGTTCG‑AAT‑G‑CCCCCTTTTTGCAATAATGTTTgcgctt                                                                                                                     <  1:488796/150‑3 (MQ=255)
                                aaTCTGTGCTAAAAACCAGGTAAGGATCAGTAGGTCAGCACTGCCGCCTGGACTGAGATTTCGTTCGATACACTCCCTGTCGAACTGCCGGAGATAATCGAGATCGGCGGGGGGTTGG‑AAG‑A‑CCCCCTTTTT‑CTATAATTTTTTTGCCTc                                                                                                                    >  2:538575/1‑150 (MQ=255)
                                  tCTGTGCTAAAAACCAGGTAAGGATCAGTAGGTCAGCACTGCCGCCTGGACTGAGATTTCGTTCGATACACTCCCTGTCGAACTGCCGGAGATGATCGAGATCGGC‑GGGGGGTGG‑GAT‑G‑GCCCCTTTTTTCTATAATTTTTTCGcc                                                                                                                      >  1:106037/1‑146 (MQ=255)
                                  tCTGTGCTAAAAACCAGGTAAGGATCAGTAGGTCAGCACTGCCGCCTGGACTGAGATTTCGTTCGATACACTCCCTGTCGAACTGCCGGAGATAATCGAGATCGGC‑GGGGGGTCG‑GAT‑G‑CCCCCCTTTTTCAATAATGtttttgcccccc                                                                                                                  >  2:157386/1‑141 (MQ=255)
                                    tgtgCTAAAAACCAGGTAAGGATCAGTAGGTCAGCACTGCCGCCTGGACTGAGATTTCGTTCGATACACTCCCTGTCGAACTGCCGGAGATAATCGAGATCGGC‑GGGGGGTCG‑AAT‑A‑ACCCCCTTTTTCAAAAATGttttttcctccc                                                                                                                  >  2:586480/1‑139 (MQ=255)
                                     gtgCTAAAAACCAGGTAAGGATCAGTAGGTCAGCACTGCCGCCTGGACTGAGATTTCGTTCGATACACTCCCTGTCGAACTGCCGGAGATAATCGAGATCGGC‑GGGGGTTCG‑AAT‑G‑CCCCCTTTTTGCAATAATGTTTGCGCCTcgcg                                                                                                                 <  1:1078229/148‑1 (MQ=255)
                                     gtgCTAAAAACCAGGTAAGGATCAGTAGGTCAGCACTGCCGCCTGGACTGAGATTTCGTTCGATACACTCCCTGTCGAACTGCCGGAGATAATCGAGATCGGC‑GGGGGTTCG‑AAT‑G‑CCCCCTTTTTGCAATAATGTTTGCGCCTcgc                                                                                                                  <  1:1279321/147‑1 (MQ=255)
                                      tgCTAAAAACCAGGTAAGGATCAGTAGGTCAGCACTGCCGCCTGGACTGAGATTTCGTTCGATACACTCCCTGTCGAACTGCCGGAGATAATCGAGATCGGC‑GGGGGTTTG‑ACA‑G‑CCCCCCCTTTTCTAAAATGtttttgccccccg                                                                                                                 >  2:1017474/1‑137 (MQ=255)
                                            aaaCCAGGTAAGGATCAGTAGGTCAGCACTGCCGCCTGGACTGAGATTTCGTTCGATACACTCCCTGTCGAACTGCCGGAGATAATCGAGATCGGC‑GGGGGGTTG‑GAT‑G‑CCCCCCTTTTTTAATAATTTTTTTGCCTccccc                                                                                                                >  1:898975/1‑138 (MQ=255)
                                              aCCAGGTAAGGATCAGTAGGTCAGCACTGCCGCCTGGACTGAGATTTCGTTCGATACACTCCCTGTCGAACTGCCGGAGATAATCGAGATCGGC‑GGGGGTTCG‑AAT‑G‑CCCCCTTTTTGCAATAATGTTTGCGCCTCGCGCTGTAGCCAgcg                                                                                                     <  1:601650/151‑1 (MQ=255)
                                                       ggATCAGTAGGTCAGCACTGCCGCCTGGACTGAGATTTCGTTCGATACACTCCCTGTCGAACTGCCGGAGATAATCGAGATCGGC‑GGGGGTTCG‑AAT‑G‑CCCCCTTTTTGCAATAATGTTTGCGCCTCGCGCTGTAGCCAGCG‑CAGGccccc                                                                                           <  2:64709/151‑1 (MQ=255)
                                                         aTCAGTAGGTCAGCACTGCCGCCTGGACTGAGATTTCGTTCGATACACTCCCTGTCGAACTGCCGGAGATAATCGAGATCGGCGGGGGGTTGG‑AAG‑C‑CCCCCTTTTTTCATAAATTTTTTGGCC‑CCCGCTTTAGCCCGCC‑CCGGgcccccc                                                                                         >  1:1011105/1‑150 (MQ=255)
                                                         aTCAGTAGGTCAGCACTGCCGCCTGGACTGAGATTTCGTTCGATACACTCCCTGTCGAACTGCCGGAGATAATCGAGATCGGC‑GGGGGGTCG‑AAA‑G‑CCCCCCTTTTTTTATAAATTTTTTGCCTCCCGCCGTTGCCCGCG‑CCGGccccccc                                                                                         >  2:1275710/1‑150 (MQ=255)
                                                           cAGTAGGTCAGCACTGCCGCCTGGACTGAGATTTCGTTCGATACACTCCCTGTCGAACTGCCGGAGATAATCGAGATCGGC‑GGGGGTTCG‑AAT‑G‑CCCCCTTTTTGCAATAATGTTTGCGCCTCGCGCTGTAGCCAGCG‑CAGGCCCCCCTCg                                                                                       <  2:1128564/151‑1 (MQ=255)
                                                            aGTAGGTCAGCACTGCCGCCTGGACTGAGATTTCGTTCGATACACTCCCTGTCGAACTGCCGGAGATAATCGAGATCGGC‑GGGGGGTCGAAAT‑G‑CCCCCTTTTT‑CAATAAAGTTTTCGCCCCGCGCTCTAGACCGCG‑CAGGcccccc                                                                                          >  2:613590/1‑147 (MQ=255)
                                                              tAGGTCAGCACTGCCGCCTGGACTGAGATTTCGTTCGATACACTCCCTGTCGAACTGCCGGAGATAATCGAGATCGGC‑GGGGGCTTG‑AAG‑GCCCCCCTTTTTGTATTAATGTTTTC‑CCCCGCCCTGTTGGCTGCG‑CAGGgccccccccccc                                                                                    >  2:780429/1‑145 (MQ=255)
                                                                   cAGCACTGCCGCCTGGACTGAGATTTCGTTCGATACACTCCCTGTCGAACTGCCGGAGATAATCGAGATCGGC‑GGGGGTTCG‑AAT‑G‑CCCCCTTTTTGCAATAATGTTTGCGCCTCGCGCTGTAGCCAGCG‑CAGGCCCCCCTCGCCACcgcg                                                                               <  2:19118/151‑1 (MQ=255)
                                                                          gccgccTGGACTGAGATTTCGTTCGATACACTCCCTGTCGAACTGCCGGAGATAATCGAGATCGGC‑GGGGGTTCG‑AAT‑G‑CCCCCTTTTTGCAATAATGTTTGCGCCTCGCGCTGTAGCCAGCG‑CAGGCCCCCCTCGCCACCGCGCGATg                                                                          <  2:173672/149‑1 (MQ=255)
                                                                            cgccTGGACTGAGATTTCGTTCGATACACTCCCTGTCGAACTGCCGGAGATAATCGAGATCGGC‑GGGGGGTTG‑AAT‑G‑CCCCCCTTTTTTAATAATTTTTTCGCCCCCCGCCGTAGGCCGCG‑CCAGCCCCCCCCCCCCCCCCCCGCTGCaa                                                                       >  1:1066959/1‑150 (MQ=255)
                                                                            cgccTGGACTGAGATTTCGTTCGATACACTCCCTGTCGAACTGCCGGAGATAATCGAGATCGGC‑GGGGGGTCG‑AAA‑G‑CCCCCCTTTTTTAATAATGTTTTCGCCCCCCGCTGTAGCCAGCG‑CAGGCCCCCCCCCCCCCCCCGCGATGCaa                                                                       >  2:565881/1‑150 (MQ=255)
                                                                             gccTGGACTGAGATTTCGTTCGATACACTCCCTGTCGAACTGCCGGAGATAATCGAGATCGGC‑GGTGGGTCG‑ACT‑G‑CCCCCCTTTTTCAATAATATTTTTTCCCGTCGCTGTGACCAGCG‑CCGGACCCCCCCCCCAccccccggagcc                                                                        >  1:694050/1‑138 (MQ=255)
                                                                                      gagaTTTCGTTCGATACACTCCCTGTCGAACTGCCGGAGATAATCGAGATCGGC‑GGGGGTTCG‑AAT‑G‑CCCCCTTTTTGCAATAATGTTTGCGCCTCGCGCTGTAGCCAGCG‑CAGGCCCCCCTCGCCACCGCGCGATGCAACGTTGGTATCg                                                            <  1:1192767/151‑1 (MQ=255)
                                                                                      gagaTTTCGTTCGATACACTCCCTGTCGAACTGCCGGAGATAATCGAGATCGGC‑GGGGGGTCG‑AAT‑G‑CCCCCCTTTTTTAATAATGTTTTCGCCCCGCCCTGGAGCCAGCG‑CAGGGCCCCCCCGCCCCCCCGCCGTGCAACGGTGGGAt                                                              >  2:436215/1‑149 (MQ=255)
                                                                                       agaTTTCGTTCGATACACTCCCTGTCGAACTGCCGGAGATAATCGAGATCGGC‑GGGGGGTCC‑CAC‑G‑GCCCCCTTTTGCAATAATTATTTCGGCTCCCGCGGCGGGCAGCG‑CAGGCCCCCCCCACCACCCCGCGCTGCAACGtttggat                                                              >  1:357202/1‑143 (MQ=255)
                                                                                        gaTTTCGTTCGATACACTCCCTGTCGAACTGCCGGAGATAATCGAGATCGGC‑GGGGGTTCG‑AAT‑G‑CCCCCTTTTTGCAATAATGTTTGCGCCTCGCGCTGTAGCCAGCG‑CAGGCCCCCCTCGCCACCGCGCGATGCAACGTTGGTATCTcc                                                          <  2:235657/151‑1 (MQ=255)
                                                                                          tttCGTTCGATACACTCCCTGTCGAACTGCCGGAGATAATCGAGATCGGC‑GGGGGGTTG‑AAA‑G‑CCCCCTTTTTGCAATAAATTTTTTGCCTCCCGCCGTAGCCAGGG‑CACGCCCCCCCCCCCACCCCCCGGTGCAAAGTTGGTATTGcc                                                          >  2:1171430/1‑149 (MQ=255)
                                                                                                tcgATACACTCCCTGTCGAACTGCCGGAGATAATCGAGATCGGC‑GGGGGTTCG‑AAT‑G‑CCCCCTTTTTGCAATAATGTTTGCGCCTCGCGCTGTAGCCAGCG‑CAGGCCCCCCTCGCCACCGCGCGATGCAACGTTGGTATCGCCGTTGATc                                                   <  2:938418/150‑1 (MQ=255)
                                                                                                tcgATACACTCCCTGTCGAACTGCCGGAGATAATCGAGATCGGC‑GGGGGTTCG‑AAT‑G‑CCCCCTTTTTGCAATAATGTTTGCGCCTCGCGCTGTAGCCAGCG‑CAGGCCCCCCTCGCCACCGCGCGATGCAACGTTGGTATCGCCGTTGATc                                                   <  1:4144/150‑1 (MQ=255)
                                                                                                 cgATACACTCCCTGTCGAACTGCCGGAGATAATCGAGATCGGC‑GGGGGGTCG‑AAT‑G‑CCCCCCTTTTTCAATAATATTTTCGCCCCCCGCCGTCGCCAGCG‑CAGGCCCCCCCCCCCCCCCCCCGGTGCAACGTTGGGATCGCCGTTGAt                                                    >  2:1208179/1‑148 (MQ=255)
                                                                                                     acacTCCCTGTCGAACTGCCGGAGATAATCGAGATCGGC‑GGGGGGTCG‑TAT‑G‑CCCCCCTTTTTTCATAACTTTTTTGCCCCCCCCTGTAGCCAGCGCCCGGCCCCCCTC‑CCCCCCCCCGCTGCGAAGTGGGGGTCGCCGCTGATCGcc                                                >  1:750220/1‑148 (MQ=255)
                                                                                                      cacTCCCTGTCGAACTGCCGGAGATAATCGAGATCGGC‑GGGGGTTCG‑AAT‑G‑CCCCCTTTTTGCAATAATGTTTGCGCCTCGCGCTGTAGCCAGCG‑CAGGCCCCCCTCGCCACCGCGCGATGCAACGTTGGTATCGCCGTTGATCGCCATCa                                            <  2:694984/151‑1 (MQ=255)
                                                                                                         tCCCTGTCGAACTGCCGGAGATAATCGAGATCGGC‑GGGGGGTTG‑ACA‑G‑CCCCCCTTTTTCAATAAAATTTTTGCCCCCCGCCGTAGCCAGCG‑CCGGGCCCCCCCCCCCCCCCGCGGAGCACCGTTTGGATCGCCGGTGGTCGCCCTc                                             >  1:16498/1‑147 (MQ=255)
                                                                                                         tCCCTGTCGAACTGCCGGAGATAATCGAGATCGGC‑GGGGGGTTG‑AAA‑G‑CCCCCCTTTTTCAATAAAGTTTTCGCCTCCCGCCGGAGCCAGCG‑CCGGGCCCCCCCCCCCCCCCCCGGTGCCACGGTGGGATCGCCGGTGATCGCCCTc                                             >  2:658478/1‑147 (MQ=255)
                                                                                                          cccTGTCGAACTGCCGGAGATAATCGAGATCGGC‑GGGGGGTCG‑GAT‑G‑CCCCCCTTTTTCAATAATATTTTTGCCCCCCGCTGTAGCCAGCG‑CCGGCCCCCCCCCCCCCCCCCCCATGCAACGTTGGGATTGCCGTTGGTCGCCCTCAGt                                          >  1:798353/1‑149 (MQ=255)
                                                                                                           ccTGTCGAACTGCCGGAGATAATCGAGATCGGC‑GGGGGTTCG‑AAT‑G‑CCCCCTTTTTGCAATAATGTTTGCGCCTCGCGCTGTAGCCAGCG‑CAGGCCCCCCTCGCCACCGCGCGATGCAACGTTGGTATCGCCGTTGATCGCCATCAGTAgg                                       <  2:1270077/151‑1 (MQ=255)
                                                                                                              gTCGAACTGCCGGAGATAATCGAGATCGGC‑GGGGGTTCG‑AAT‑G‑CCCCCTTTTTGCAATAATGTTTGCGCCTCGCGCTGTAGCCAGCG‑CAGGCCCCCCTCGCCACCGCGCGATGCAACGTTGGTATCGCCGTTGATCGCCATCAGTAGGAg                                     <  2:264718/150‑1 (MQ=255)
                                                                                                              gTCGAACTGCCGGAGATAATCGAGATCGGC‑GGGGGTTCG‑AAT‑G‑CCCCCTTTTTGCAATAATGTTTGCGCCTCGCGCTGTAGCCAGCG‑CAGGCCCCCCTCGCCACCGCGCGATGCAACGTTGGTATCGCCGTTGATCGCCATCAGTAGGAg                                     <  2:696232/150‑1 (MQ=255)
                                                                                                              gTCGAACTGCCGGAGATAATCGAGATCGGC‑GGGGGTTCG‑AAT‑G‑CCCCCTTTTTGCAATAATGTTTGCGCCTCGCGCTGTAGCCAGCG‑CAGGCCCCCCTCGCCACCGCGCGATGCAACGTTGGTATCGCCGTTGATCGCCATCAGTAGGAg                                     <  2:468992/150‑1 (MQ=255)
                                                                                                              gTCGAACTGCCGGAGATAATCGAGATCGGC‑GGGGGGTCG‑GATAG‑CCCCCTTTTT‑CAAAAAAATTTTCGCCTCCCGCCGTGGCCCGCG‑CCGGGCCCCCCCCCCACCCCCCCATGCAACGCCGGGATCGCCCTTGGTCGCCCCCAGgaggag                                     >  1:755691/1‑150 (MQ=255)
                                                                                                                 gAACTGCCGGAGATAATCGAGATCGGC‑GGGGGTTCG‑GAT‑G‑CCCCCCTTTTTCAATAATGTTTTCGCCCCGCGCCGTAGCCAGCG‑CAGGCCCCCCCCCCCCCCCCGCGCTGCCACGGTGGGATCGCCGTTGGTCGCCCTCAGTAGGAg                                     >  2:269387/1‑147 (MQ=255)
                                                                                                                  aaCTGCCGGAGATAATCGAGATCGGC‑GGGGGGTCG‑AAA‑G‑CCCCCCTTTTTCAATAATGATTTCGCCCCCCGCTGTAGGCAGCG‑CCGGCCCCCCCCCCCCCCCCCCGCGGCCACGGTTGGATCGCCGTTGTTCGCCATCaggagggggaa                                  >  2:434381/1‑144 (MQ=255)
                                                                                                                  aaCTGCCGGAGATAATCGAGATCGGC‑GGGGGGTCC‑GAT‑G‑CCCCCCTTTTTTAAAAAATTTTTCCCCCCCCCCTGTAGCCAGCG‑CCGGCCCCCCCCCCCACCCCCCGGAGCAAAGTTTGTATCGCCCTTTGTCGCCATCCGTAGGAGGAAgg                                >  2:1210050/1‑151 (MQ=255)
                                                                                                                   aCTGCCGGAGATAATCGAGATCGGC‑GGGGGGTCG‑GAT‑G‑GCCCCCTTTTTTAATAATATTTGCGCCTCGCGCGGGTGGCCGGG‑CCGGGCCCCCCCCCCACCCCCCGACGGAAAGGTGGGATCGCCCCTGGTCGCCATCCGTAGGAgggagg                                >  1:783870/1‑145 (MQ=255)
                                                                                                                      gCCGGAGATAATCGAGATCGGC‑GGGGGGTCG‑AAT‑T‑CCCCCCTTTTGTAATAATGTTTTCGCCTCCCGCTGTAGCCAGCG‑CAGGACCCCCCCGCCCCCCCGCCATGCCACGTTGGGATCGCCGTTGATTGCCCTCAATAGGGGCAAgg                                >  2:493691/1‑147 (MQ=255)
                                                                                                                       ccGGAGATAATCGAGATCGGC‑GGGGGTTCG‑AAT‑G‑CCCCCTTTTTGCAATAATGTTTGCGCCTCGCGCTGTAGCCAGCG‑CAGGCCCCCCTCGCCACCGCGCGATGCAACGTTGGTATCGCCGTTGATCGCCATCAGTAGGAGCAAGGTATCg                           <  1:335496/151‑1 (MQ=255)
                                                                                                                        cGGAGATAATCGAGATCGGC‑GGGGGTTCG‑AAT‑G‑CCCCCTTTTTGCAATAATGTTTGCGCCTCGCGCTGTAGCCAGCG‑CAGGCCCCCCTCGCCACCGCGCGATGCAACGTTGGTATCGCCGTTGATCGCCATCAGTAGGAGCAAGGTATCGa                          <  2:527347/151‑1 (MQ=255)
                                                                                                                         ggAGATAATCGAGATCGGC‑GGGGGTTCG‑AAT‑G‑CCCCCTTTTTGCAATAATGTTTGCGCCTCGCGCTGTAGCCAGCG‑CAGGCCCCCCTCGCCACCGCGCGATGCAACGTTGGTATCGCCGTTGATCGCCATCAGTAGGAGCAAGGTATCGAg                         <  2:495148/151‑1 (MQ=255)
                                                                                                                          gagaTAATCGAGATCGGC‑GGGGGTTCG‑AAT‑G‑CCCCCTTTTTGCAATAATGTTTGCGCCTCGCGCTGTAGCCAGCG‑CAGGCCCCCCTCGCCACCGCGCGATGCAACGTTGGTATCGCCGTTGATCGCCATCAGTAGGAGCAAGGTATCGAGc                        <  1:739211/151‑1 (MQ=255)
                                                                                                                           agaTAATCGAGATCGGC‑GGGGGTTCG‑AAT‑G‑CCCCCTTTTTGCAATAATGTTTGCGCCTCGCGCTGTAGCCAGCG‑CAGGCCCCCCTCGCCACCGCGCGATGCAACGTTGGTATCGCCGTTGATCGCCATCAGTAGGAGCAAGGTATCGAGCa                       <  1:117939/151‑1 (MQ=255)
                                                                                                                                      aTCGGC‑GGGGGTTCG‑AAT‑G‑CCCCCTTTTTGCAATAATGTTTGCGCCTCGCGCTGTAGCCAGCG‑CAGGCCCCCCTCGCCACCGCGCGATGCAACGTTGGTATCGCCGTTGATCGCCATCAGTAGGAGCAAGGTATCGAGCAATGCCAGTTCa            <  1:538560/151‑1 (MQ=255)
                                                                                                                                       tCGGC‑GGGGGTTCG‑AAT‑G‑CCCCCTTTATGCAATAATGTTTGCGCCTCGCGCTGTAGCCAGCG‑CAGGCCGCCCTCGCCACCGCGCGATGCAACGTTGGTATCGCCGTTGATCGCCATCAGTAGGAGCAAGGTATCGGGCAATGCCAGTTCa            <  1:971158/150‑1 (MQ=255)
                                                                                                                                        cggc‑ggGGGTTCG‑AAT‑G‑CCCCCTTTTTGCAATAATGTTTGCGCCTCGCGCTGTAGCCAGCG‑CAGGCCCCCCTCGCCACCGCGCGATGCAACGTTGGTATCGCCGTTGATCGCCATCAGTAGGAGCAAGGTATCGAGCAATGCCAGTTCAg           <  1:164710/150‑1 (MQ=255)
                                                                                                                                           c‑ggGGGTTCG‑AAT‑G‑CCCCCTTTTTGCAATAATGTTTGCGCCTCGCGCTGTAGCCAGCG‑CAGGCCCCCCTCGCCACCGCGCGATGCAACGTTGGTATCGCCGTTGATCGCCATCAGTAGGAGCAAGGTATCGAGCAATGCCAGTTCAGGATc       <  2:258255/151‑1 (MQ=255)
                                                                                                                                           c‑ggGGGTTCG‑AAT‑G‑CCCCCTTTTTGCAATAATGTTTGCGCCTCGCGCTGTAGCCAGCG‑CAGGCCCCCCTCGCCACCGCGCGATGCAACGTTGGTATCGCCGTTGATCGCCATCAGTAGGAGCAAGGTATCGAGCAATGCCAGTTCAGGATc       >  1:1217454/1‑151 (MQ=255)
                                                                                                                                              ggggTTCG‑AAT‑G‑CCCCCTTTTTGCAATAATGTTTGCGCCTCGCGCTGTAGCCAGCG‑CAGGCCCCCCTCGCCACCGCGCGATGCAACGTTGGTATCGCCGTTGATCGCCATCAGTAGGAGCAAGGTATCGAGCAATGCCAGTTCAGGATCt      >  1:1205001/1‑150 (MQ=255)
                                                                                                                                              ggggTTCG‑AAT‑G‑CCCCCTTTTTGCAATAATGTTTGCGCCTCGCGCTGTAGCCAGCG‑CAGGCCCCCCTCGCCACCGCGCGATGCAACGTTGGTATCGCCGTTGATCGCCATCAGTAGGAGCAAGGTATCGAGCAATGCCAGTTCAGGATCTa     <  1:605403/151‑1 (MQ=255)
                                                                                                                                                 gTTCG‑AAT‑G‑CCCCCTTTTTGCAATAATGTTTGCGCCTCGCGCTGTAGCCAGCG‑CAGGCCCCCCTCGCCACCGCGCGATGCAACGTTGGTATCGCCGTTGATCGCCATCAGTAGGAGCAAGGTATCGAGCAATGCCAGTTCAGGATCTAAcc  >  2:444597/1‑151 (MQ=255)
                                                                                                                                                  |                                                                                                                                                         
ATTTCCAAATTTATCAAGTGCTTAAATAATTAAATCTGTGCTAAAAACCAGGTAAGGATCAGTAGGTCAGCACTGCCGCCTGGACTGAGATTTCGTTCGATACACTCCCTGTCGAACTGCCGGAGATAATCGAGATCGGC‑GGGGGTTCG‑AAT‑G‑CCCCCTTTTTGCAATAATGTTTGCGCCTCGCGCTGTAGCCAGCG‑CAGGCCCCCCTCGCCACCGCGCGATGCAACGTTGGTATCGCCGTTGATCGCCATCAGTAGGAGCAAGGTATCGAGCAATGCCAGTTCAGGATCTAACC  >  NC_000913/646602‑646896

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