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

ATTTCCAAATTTATCAAGTGCTTAAATAATTAAATCTGTGCTAAAAACCAGGTAAGGATCAGTAGGTCAGCACTGCCGCCTGGACTGAGATTTCGTTCGATACACTCCCTGTCGAACTGCCGGAGATAATCGAGATCGGCGGGGGTTC‑GAATGCCCCCTTTTTGCAATAATGTTTGCGCCTCGCGCTGTAGCCAGCG‑CAGGCCCCCCTCGCCACCGCGCGATGCAACGTTGGTATCGCCGTTGATCGCCATCAGTAGGAGCAAGGTATCGAGCAATGCCAGTTCAGGATCTAACCC  >  NC_000913/646602‑646897
                                                                                                                                                 |                                                                                                                                                        
atTTCCAAATTTATCAAGTGCTTAAATAATTAAATCTGTGCTAAAAACCAGGTAAGGATCAGTAGGTCAGCACTGCCGCCTGGACTGAGATTTCGTTCGATACACTCCCTGTCGAACTGCCGGAGATAATCGAGATCGGCGGGGGGTC‑gt                                                                                                                                                     >  2:46325/1‑149 (MQ=255)
   tCCAAATTTATCAAGTGCTTAAATAATTAAATCTGTGCTAAAAACCAGGTAAGGATCAGTAGGTCAGCACTGCCGCCTGGACTGAGATTTCGTTCGATACACTCCCTGTCGAACTGCCGGAGATAATCGAGATCGGCGGGGGTTC‑GAATg                                                                                                                                                  <  2:1108362/150‑1 (MQ=255)
    ccAAATTTATCAAGTGCTTAAATAATTAAATCTGTGCTAAAAACCAGGTAAGGATCAGTAGGTCAGCACTGCCGCCTGGACTGAGATTTCGTTCGATACACTCCCTGTCGAACTGCCGGAGATAATCGAGATCGGCGGGGGTTC‑GAATGc                                                                                                                                                 <  2:328170/150‑1 (MQ=255)
      aaaTTTATCAAGTGCTTAAATAATTAAATCTGTGCTAAAAACCAGGTAAGGATCAGTAGGTCAGCACTGCCGCCTGGACTGAGATTTCGTTCGATACACTCCCTGTCGAACTGCCGGAGATAATCGAGATCGGCGGGGGTTC‑GAATGcccc                                                                                                                                              <  2:365624/151‑1 (MQ=255)
       aaTTTATCAAGTGCTTAAATAATTAAATCTGTGCTAAAAACCAGGTAAGGATCAGTAGGTCAGCACTGCCGCCTGGACTGAGATTTCGTTCGATACACTCCCTGTCGAACTGCCGGAGATAATCGAGATCGGCGGGGGGTC‑GGATGccccc                                                                                                                                             >  2:1292483/1‑151 (MQ=255)
       aaTTTATCAAGTGCTTAAATAATTAAATCTGTGCTAAAAACCAGGTAAGGATCAGTAGGTCAGCACTGCCGCCTGGACTGAGATTTCGTTCGATACACTCCCTGTCGAACTGCCGGAGATAATCGAGATCGGCGGCGGGTccacattgccc                                                                                                                                              >  2:184889/1‑141 (MQ=255)
          ttATCAAGTGCTTAAATAATTAAATCTGTGCTAAAAACCAGGTAAGGATCAGTAGGTCAGCACTGCCGCCTGGACTGAGATTTCGTTCGATACACTCCCTGTCGAACTGCCGGAGATAATCGAGATCGGCGGGGGGTC‑GGATGCCCCCctt                                                                                                                                          >  1:631240/1‑151 (MQ=255)
          ttATCAAGTGCTTAAATAATCAAATCTCTGCTAAAAACCAGGTAATGATCAGTAGGTCAGCACTGCCGCATGGACGGAGATTTCGTTCGATATACTCCCTGTCGAACTGCCGGAGATAATCGAGATCGGCGGGGGTTC‑GAAATACCCCtt                                                                                                                                           >  1:631285/1‑150 (MQ=255)
           tATCAAGTGCTTAAATAATTAAATCTGTGCTAAAAACCAGGTAAGGATCAGTAGGTCAGCACTGCCGCCTGGACTGAGATTTCGTTCGATACACTCCCTGTCGAACTGCCGGAGATAATCGAGATCGGCGGGGGGTT‑GAATGCCCCCcttt                                                                                                                                         >  2:776920/1‑151 (MQ=255)
            aTCAAGTGCTTAAATAATTAAATCTGTGCTAAAAACCAGGTAAGGATCAGTAGGTCAGCACTGCCGCCTGGACTGAGATTTCGTTCGATACACTCCCTGTCGAACTGCCGGAGATAATCGAGATCGGCGGGGGTTC‑GAATGCCCCCttttt                                                                                                                                        <  2:1258468/151‑1 (MQ=255)
                aGTGCTTAAATAATTAAATCTGTGCTAAAAACCAGGTAAGGATCAGTAGGTCAGCACTGCCGCCTGGACTGAGATTTCGTTCGATACACTCCCTGTCGAACTGCCGGAGATAATCGAGATCGGCGGGGGTTC‑GAATGCCCCCTTTTTGCaa                                                                                                                                    <  1:811295/151‑1 (MQ=255)
                      tAAATAATTAAATCTGTGCTAAAAACCAGGTAAGGATCAGTAGGTCAGCACTGCCGCCTGGACTGAGATTTCGTTCGATACACTCCCTGTCGAACTGCCGGAGATAATCGAGATCGGCGGGGGTTC‑GAATGCCCCCTTTTTGCAATAATg                                                                                                                               <  1:614921/150‑1 (MQ=255)
                        aataatTAAATCTGTGCTAAAAACCAGGTAAGGATCAGTAGGTCAGCACTGCCGCCTGGACTGAGATTTCGTTCGATACACTCCCTGTCGAACTGCCGGAGATAATCGAGATCGGCGGGGGTTC‑GAATGCCCCCTTTTTGCAATAATGttt                                                                                                                            <  2:99779/151‑1 (MQ=255)
                         ataatTAAATCTGTGCTAAAAACCAGGTAAGGATCAGTAGGTCAGCACTGCCGCCTGGACTGAGATTTCGTTCGATACACTCCCTGTCGAACTGCCGGAGATAATCGAGATCGGCGGGGGTTC‑GAATGCCCCCTTTTTGCAATAATGTTTg                                                                                                                           <  1:581432/151‑1 (MQ=255)
                         ataatTAAATCTGTGCTAAAAACCAGGTAAGGATCAGTAGGTCAGCACTGCCGCCTGGACTGAGATTTCGTTCGATACACTCCCTGTCGAACTGCCGGAGATAATCGAGATCGGCGGGGGGTC‑GCATTCCCCCCTTTTTtaataaaatttt                                                                                                                           >  1:425093/1‑145 (MQ=255)
                               aaaTCTGTGCTAAAAACCAGGTAAGGATCAGTAGGTCAGCACTGCCGCCTGGACTGAGATTTCGTTCGATACACTCCCTGTCGAACTGCCGGAGATAATCGAGATCGGCGGGGGTTC‑GAATGCCCCCTTTTTGCAATAATGTTTGCGCCTc                                                                                                                     <  2:601633/151‑1 (MQ=255)
                                 aTCTGTGCTAAAAACCAGGTAAGGATCAGTAGGTCAGCACTGCCGCCTGGACTGAGATTTCGTTCGATACACTCCCTGTCGAACTGCCGGAGATAATCGAGATCGGCGGGGGGTC‑GAATGCCCCCCTTTTTTAATAAAGTTTTCGCCTc                                                                                                                     >  1:726150/1‑149 (MQ=255)
                                   cTGTGCTAAAAACCAGGTAAGGATCAGTAGGTCAGCACTGCCGCCTGGACTGAGATTTCGTTCGATACACTCCCTGTCGAACTGCCGGAGATAATCGAGATCGGCGGGGGGTC‑GAATGCCCCCCTTTTTCTATAATTTTTTCGcccc                                                                                                                     >  2:812625/1‑145 (MQ=255)
                                     gtgCTAAAAACCAGGTAAGGATCAGTAGGTCAGCACTGCCGCCTGGACTGAGATTTCGTTCGATACACTCCCTGTCGAACTGCCGGAGATAATCGAGATCGGCGGGGGTTC‑GAATGCCCCCTTTTTGCAATAATGTTTGCGCCTCgcgc                                                                                                                 <  1:1062172/149‑1 (MQ=255)
                                     gtgCTAAAAACCAGGTAAGGATCAGTAGGTCAGCACTGCCGCCTGGACTGAGATTTCGTTCGATACACTCCCTGTCGAACTGCCGGAGATAATCGAGATCGGCGGGGGTTC‑GAATGCCCCCTTTTTGCAATAATGTTTGCGCCTCGCGCt                                                                                                                <  1:583776/150‑1 (MQ=255)
                                          aaaaaCCAGGTAAGGATCAGTAGGTCAGCACTGCCGCCTGGACTGAGATTTCGTTCGATACACTCCCTGTCGAACTGCCGGAGATAATCGAGATCGGCGGGGGTTC‑GAATGCCCCCTTTTTGCAATAATGTTTGCGCCTCGCGCTGTAGc                                                                                                           <  2:1111156/150‑1 (MQ=255)
                                          aaaaaCCAGGTAAGGATCAGTAGGTCAGCACTGCCGCCTGGACTGAGATTTCGTTCGATACACTCCCTGTCGAACTGCCGGAGATAATCGAGATCGGCGGGGGGTC‑GAATTCCCCCCTTTTTCAATAATGTTTGCGCCTCGCGCCGTAGcc                                                                                                          >  2:156288/1‑151 (MQ=255)
                                          aaaaaCCAGGTAAGGATCAGTAGGTCAGCACTGCCGCCTGGACTGAGATTTCGTTCGATACACTCCCTGTCGAACTGCCGGAGAT‑ATCGAGATCGGCGGGGGGTC‑GGATGCCCCCCTTTTTCAATAATGTTTTCGGCTCGCCCCGTAGcc                                                                                                          >  2:252603/1‑150 (MQ=255)
                                            aaaCCAGGTAAGGATCAGTAGGTCAGCACTGCCGCCTGGACTGAGATTTCGTTCGATACACTCCCTGTCGAACTGCCGGAGATAATCGAGATCGGCGGGGGGTC‑GAAAGCCCCCTTTTTGCTATAAAGTTTGCGCCTCCCGCCGTAGcc                                                                                                          >  1:378532/1‑149 (MQ=255)
                                              aCCAGGTAAGGATCAGTAGGTCAGCACTGCCGCCTGGACTGAGATTTCGTTCGATACACTCCCTGTCGAACTGCCGGAGATAATCGAGATCGGCGGGGGTTC‑GAATGCCCCCTTTTTGCAATAATGTTTGCGCCTCGCGCTGTAGCCAgcg                                                                                                      <  1:103226/151‑1 (MQ=255)
                                              aCCAGGTAAGGATCAGTAGGTCAGCACTGCCGCCTGGACTGAGATTTCGTTCGATACACTCCCTGTCGAACTGCCGGAGATAATCGAGATCGGCGGGGGTTC‑GAATGCCCCCTTTTTGCAATAATGTTTGCGCCTCGCGCTGTAGCCAgcg                                                                                                      <  1:981912/151‑1 (MQ=255)
                                               ccAGGTAAGGATCAGTAGGTCAGCACTGCCGCCTGGACTGAGATTTCGTTCGATACACTCCCTGTCGAACTGCCGGAGATAATCGAGATCGGCGGGGGGTC‑GAATGCCCCCCTTTTTTAATAATGTTTTCGCCTCGCCCCGTTGcc                                                                                                          >  1:714048/1‑146 (MQ=255)
                                                cAGGTAAGGATCAGTAGGTCAGCACTGCCGCCTGGACTGAGATTTCGTTCGATACACTCCCTGTCGAACTGCCGGAGATAATCGAGATCGGCGGGGGTTC‑GAATGCCCCCTTTTTGCAATAATGTTTGCGCCTCGCGCTGTAGCCAGCG‑Ca                                                                                                   <  2:1118860/151‑1 (MQ=255)
                                                cAGGTAAGGATCAGTAGGTCAGCACTGCCGCCTGGACTGAGATTTCGTTCGATACACTCCCTGTCGAACTGCCGGAGATAATCGAGATCGGCGGGGGGTT‑GAATGCCCCCCTTTTTCAAAAATGTTTTTGCCCCGCGCTGTAGcccccgc                                                                                                     >  2:677123/1‑145 (MQ=255)
                                                  ggTAAGGATCAGTAGGTCAGCACTGCCGCCTGGACTGAGATTTCGTTCGATACACTCCCTGTCGAACTGCCGGAGATAATCGAGATCGGCGGGGGTTC‑GAATGCCCCCTTTTTGCAATAATGTTTGCGCCTCGCGCTGTAGCCAGCG‑CAgg                                                                                                 <  1:954729/151‑1 (MQ=255)
                                                  ggTAAGGATCAGTAGGTCAGCACTGCCGCCTGGACTGAGATTTCGTTCGATACACTCCCTGTCGAACTGCCGGAGATAATCGAGATCGGCGGGGGTTC‑GAATGCCCCCTTTTTGCAATAATGTTTGCGCCTCGCGCTGTAGCCAGCG‑CAgg                                                                                                 <  1:627939/151‑1 (MQ=255)
                                                  ggTAAGGATCAGTAGGTCAGCACTGCCGCCTGGACTGAGATTTCGTTCGATACACTCCCTGTCGAACTGCCGGAGATAATCGAGATCGGCGGGGGGTC‑GAATGCCCCCCTTTTTTAATAATGTTTTTGCCCCGCGCCGTAGCCAGCG‑CAgg                                                                                                 >  2:1114426/1‑151 (MQ=255)
                                                   gTAAGGATCAGTAGGTCAGCACTGCCGCCTGGACTGAGATTTCGTTCGATACACTCCCTGTCGAACTGCCGGAGATAATCGAGATCGGCGGGGGGTG‑GAATGGCCCCCTTTTTTAAAAAAATTTTTGCCCCGCGCGGTAGCCAgcg‑c                                                                                                    >  1:506035/1‑147 (MQ=255)
                                                   gTAAGGATCAGTAGGTCAGCACTGCCGCCTGGACTGAGATTTCGTTCGATACACTCCCTGTCGAACTGCCGGAGATAATCGAGATCGGCGGGGGGTC‑GAATGGCCCCCTTTTTCAATAATGTTTTTGCCCCCCGCTGGAGCCAGCG‑CCgg                                                                                                 >  1:873428/1‑150 (MQ=255)
                                                   gTAAGGATCAGTAGGTCAGCACTGCCGCCTGGACTGAGATTTCGTTCGATACACTCCCTGTCGAACTGCCGGAGATAATCGAGATCGGCGGGGGGCT‑GGATTCCCCCCTCTTGCGATAAAATGTGCACGTcgcgtttccca                                                                                                           >  1:118750/1‑134 (MQ=255)
                                                      aGGATCAGTAGGTCAGCACTGCCGCCTGGACTGAGATTTCGTTCGATACACTCCCTGTCGAACTGCCGGAGATAATCGAGATCGGCGGGGGGTC‑GAATGCCCCCTTTTTTTAATAATGTTTTCGCCTCGCCCCGTAGCCCCCC‑CCGGcccc                                                                                             >  2:1143916/1‑151 (MQ=255)
                                                           cAGTAGGTCAGCACTGCCGCCTGGACTGAGATTTCGTTCGATACACTCCCTGTCGAACTGCCGGAGATAATCGAGATCGGCGGGGGTTC‑GAATGCCCCCTTTTTGCAATAATGTTTGCGCCTCGCGCTGTAGCCAGCG‑CAGGCCCCCCTc                                                                                         <  1:26910/150‑1 (MQ=255)
                                                            aGTAGGTCAGCACTGCCGCCTGGACTGAGATTTCGTTCGATACACTCCCTGTCGAACTGCCGGAGATAATCGAGATCGGCGGGGGTTC‑GAATGCCCCCTTTTTGCAATAATGTTTGCGCCTCGCGCTGTAGCCAGCG‑CAGGCCCCCCTCGc                                                                                       <  2:896410/151‑1 (MQ=255)
                                                            aGTAGGTCAGCACTGCCGCCTGGACTGAGATTTCGTTCGATACACTCCCTGTCGAACTGCCGGAGATAATCGAGATCGGCGGGGGGTCGGAATCCCCCCTTTTT‑CAATAATGTTTTCGCCCCGCCCTGTAGCCACCGCCCGGccccccccc                                                                                        >  2:1078385/1‑148 (MQ=255)
                                                               aGGTCAGCACTGCCGCCTGGACTGAGATTTCGTTCGATACACTCCCTGTCGAACTGCCGGAGATAATCGAGATCGGCGGGGGTTC‑GAATGCCCCCTTTTTGCAATAATGTTTGCGCCTCGCGCTGTAGCCAGCG‑CAGGCCCCCCTCGCCAc                                                                                    <  1:505005/151‑1 (MQ=255)
                                                                ggTCAGCACTGCCGCCTGGACTGAGATTTCGTTCGATACACTCCCTGTCGAACTGCCGGAGATAATCGAGATCGGCGGGGGTTC‑GAATGCCCCCTTTTTGCAATAATGTTTGCGCCTCGCGCTGTAGCCAGCG‑CAGGCCCCCCTCGCCAcc                                                                                   <  2:1170041/151‑1 (MQ=255)
                                                                  tCAGCACTGCCGCCTGGACTGAGATTTCGTTCGATACACTCCCTGTCGAACTGCCGGAGATAATCGAGATCGGCGGGGGGTC‑GAATGGCCCCTTTTTTCAATAAAGTTTTCGCCTCCCGCTGTAGCCAGCG‑CAGGGCCCCCCCCCCAcccc                                                                                 >  2:562481/1‑149 (MQ=255)
                                                                   cAGCACTGCCGCCTGGACTGAGATTTCGTTCGATACACTCCCTGTCGAACTGCCGGAGATAATCGAGATCGGCGGGGGTTC‑GAATGCCCCCTTTTTGCAATAATGTTTGCGCCTCGCGCTGTAGCCAGCG‑CAGGCCCCCCTCGCCACcgc                                                                                 <  2:150321/150‑1 (MQ=255)
                                                                    aGCACTGCCGCCTGGACTGAGATTTCGTTCGATACACTCCCTGTCGAACTGCCGGAGATAATCGAGATCGGCGGGGGGTC‑GAATGCCCCCCTTTTTTAATAATTTTTTTGCCTCGCGCGGTAGCCAGCG‑CAGGCCCCCCCCCCCCCcccgc                                                                               >  2:874624/1‑151 (MQ=255)
                                                                     gCACTGCCGCCTGGACTGAGATTTCGTTCGATACACTCCCTGTCGAACTGCCGGAGATAATCGAGATCGGCGGGGGTTC‑GAATGCCCCCTTTTTGCAATAATGTTTGCGCCTCGCGCTGTAGCCAGCG‑CAGTCCCCCCTCGCCACcgcgc                                                                               <  1:979946/150‑1 (MQ=255)
                                                                       aCTGCCGCCTGGACTGAGATTTCGTTCGATACACTCCCTGTCGAACTGCCGGAGATAATCGAGATCGGCGGGGGTTC‑GAATGCCCCCTTTTTGCAATAATGTTTGCGCCTCGCGCTGTAGCCAGCG‑CAGGCCCCCCTCGCCACcgcgc                                                                               <  1:1042709/148‑1 (MQ=255)
                                                                         tGCCGCCTGGACTGAGATTTCGTTCGATACACTCCCTGTCGAACTGCCGGAGATAATCGAGATCGGCGGGGGGTT‑GAATGCCCCCCTTTTTCCATAATTTTTTCGCCCCCCCCCGTAGCCCGCC‑CAGGGCCCCCCCCCCACcgcgccggt                                                                           >  2:275057/1‑146 (MQ=255)
                                                                          gccgccTGGACTGAGATTTCGTTCGATACACTCCCTGTCGAACTGCCGGAGATAATCGAGATCGGCGGGGGTTC‑GAATGCCCCCTTTTTGCAATAATGTTTGCGCCTCGCGCTGTAGCCAGCG‑CAGGCCCCCCTCGCCACCGCGCGATg                                                                           <  1:282920/149‑1 (MQ=255)
                                                                              ccTGGACTGAGATTTCGTTCGATACACTCCCTGTCGAACTGCCGGAGATAATCGAGATCGGCGGGGGTTC‑GAATGCCCCCTTTTTGCAATAATGTTTGCGCCTCGCGCTGTAGCCAGCG‑CAGGCCCCCCTCGCCACCGCGCGATGCAAc                                                                       <  1:1250556/149‑1 (MQ=255)
                                                                                 ggACTGAGATTTCGTTCGATACACTCCCTGTCGAACTGCCGGAGATAATCGAGATCGGCGGGGGTTC‑GAATGCCCCCTTTTTGCAATAATGTTTGCGCCTCGCGCTGTAGCCAGCG‑CAGGCCCCCCTCGCCACCGCGCGATGCAACGTTg                                                                   <  2:402242/150‑1 (MQ=255)
                                                                                  gACTGAGATTTCGTTCGATACACTCCCTGTCGAACTGCCGGAGATAATCGAGATCGGCGGGGGTTC‑GAATGCCCCCTTTTTGCAATAATGTTTGCGCCTCGCGCTGTAGCCAGCG‑CAGGCCCCCCTCGCCACCGCGCGATGCAACGTTGGt                                                                 <  1:447550/151‑1 (MQ=255)
                                                                                           ttcgttcgATACACTCCCTGTCGAACTGCCGGAGATAATCGAGATCGGCGGGGGTTC‑GAATGCCCCCTTTTTGCAATAATGTTTGCGCCTCGCGCTGTAGCCAGCG‑CAGGCCCCCCTCGCCACCGCGCGATGCAACGTTGGTATCGCCGt                                                         <  2:1083735/150‑1 (MQ=255)
                                                                                              gttcgATACACTCCCTGTCGAACTGCCGGAGATAATCGAGATCGGCGGGGGGTC‑GAATGCCCCCCTTTTTCAATAAATTTTTCGCCCCCCGCCGTAGCCCGCG‑CCGGCCCCCCCCCCCCCCCCCCCGCGCCACGTTGGTATTGCCGCTGAt                                                     >  1:917028/1‑151 (MQ=255)
                                                                                               ttcgATACACTCCCTGTCGAACTGCCGGAGATAATCGAGATCGGCGGGGGGTT‑GAATGCCCCCCTTTTTCAATAATATTTTCGCCTCGCGCTTTAGCCAGCG‑CAGGCCCCCCCCCCCACCCCGCGATGCAACCTTGGTATCGCCGGTGATc                                                    >  2:177783/1‑151 (MQ=255)
                                                                                                tcgATACACTCCCTGTCGAACTGCCGGAGATAATCGAGATCGGCGGGGGGTT‑GAAAGCCCCCCTTTTTCAATAATGTTTTCGCCCCGCGCTGTAGCCCGCG‑CAGGCCCCCCCCCCCCCCCCCCGCTGCAACGGTGGGATCGCCCtttgt                                                     >  2:592091/1‑146 (MQ=255)
                                                                                                   aTACACTCCCTGTCGAACTGCCGGAGATAATCGAGATCGGCGGGGGTTC‑GAATGCCCCCTTTTTGCAATAATGTTTGCGCCTCGCGCTGTAGCCAGCG‑CAGGCCCCCCTCGCCACCGCGCGATGCAACGTTGGTATCGCCGTTGATCGCCa                                                <  1:421091/151‑1 (MQ=255)
                                                                                                    tACTCTCCCTGTCGAACTGCCGGAGATAATCGAGATCGGCGGGGGTTC‑GAATGCCCCCTTTTTGCAATAATGTTTGCGCCGCGCGCTGTAGCCAGCG‑GAGGCTCCCCTCGCCCCCGCGCGATGCAACGTTGGTATCGCTGTTGATCGccc                                                <  1:467024/150‑2 (MQ=255)
                                                                                                     acacTCCCTGTCGAACTGCCGGAGATAATCGAGATCGGCGGGGGGTC‑GCATGTCCCCCTTTTTCCATAATTTTTGCGCCCCGCGCTGTAGGCCGCG‑CCGGCCCCCCCCCCCCCCCCCCGATGCAACCTTGGGATCGCCCCTGGTCGccccc                                              >  1:1089772/1‑148 (MQ=255)
                                                                                                           ccTGTCGAACTGCCGGAGATAATCGAGATCGGCGGGGGTTC‑GAATGCCCCCTTTTTGCAATAATGTTTGCGCCTCGCGCTGTAGCCAGCG‑CAGGCCCCCCTCGCCACCGCGCGATGCAACGTTGGTATCGCCGTTGATCGCCATCAGTAgg                                        <  1:1152366/151‑1 (MQ=255)
                                                                                                             tGTCGAACTGCCGGAGATAATCGAGATCGGCGGGGGGTC‑GAATGCCCCCCTTTTTTAATAAAGTTTTTGCCCCGCCCTGTAGCCAGCG‑CCGGCCCCCCCCCCCACCCCGCGACGCAACGCTGGGGTTGCCGTTGATCGCCATCAGCAgggg                                      >  2:268025/1‑149 (MQ=255)
                                                                                                              gTCGAACTGCCGGAGATAATCGAGATCGGCGGGGGTTC‑GAATGCCCCCTTTTTGCAATAATGTTTGCGCCTCGCGCTGTAGCCAGCG‑CAGGCCCCCCTCGCCACCGCGCGATGCAACGTTGGTATCGCCGTTGATCGCCATCAGTAGGa                                       <  2:856598/149‑1 (MQ=255)
                                                                                                               tCGAACTGCCGGAGATAATCGAGATCGGCGGGGGGTC‑GAACGCCCCCCTTTTTCAATAATATTTGCGCCCCCCGCTGTAGCCAGCG‑CCGGCCCCCCCCCCCCCCCCGCGGTGCAAAGATGGTATCGCCGCTGATCGCCATCCGgaggag                                      >  2:772776/1‑149 (MQ=255)
                                                                                                                 gAACTGCCGGAGATAATCGAGATCGGCGGGGGGTC‑GAATGCCCCCCTTTTTCAATAATGTTTTCGCCCCCCGCTGTAGCCAGCG‑CAGGCCCCCCCCCCCCCCCCCCGATGCAAAGTTGGTATCGCCCGTGATCGCCCTCAGTAGGAg                                      >  2:892397/1‑147 (MQ=255)
                                                                                                                   aCTGCCGGAGATAATCGAGATCGGCGGGGGTTC‑GAATGCCCCCTTTTTGCAATAATGTTTGCGCCTCGCGCTGTAGCCAGCG‑CAGGCCCCCCTCGCCACCGCGCGATGCAACGTTGGTATCGCCGTTGATCGCCATCAGTAGGAGCAAGGt                                <  1:1240493/151‑1 (MQ=255)
                                                                                                                   aCTGCCGGAGATAATCGAGATCGGCGGGGGTTC‑GAATGCCCCCTTTTTGCAATAATGTTTGCGCCTCGCGCTGTAGCCAGCG‑CAGGCCCCCCTCGCCACCGCGCGATGCAACGTTGGTATCGCCGTTGATCGCCATCAGTAGGAGCAAGGt                                <  1:1151994/151‑1 (MQ=255)
                                                                                                                     tGCCGGAGATAATCGAGATCGGCGGGGGGTC‑GAAAGCCCCCCTTTTTCAATAAATTTTTCGCCTCCCGCTGGAGCCCGCG‑CCGGCCCCCCCCCCCCCCCCCCGATGCAAAAGTGGTATCGCCGCTGATCGCCATCCGTAGGAGGAAgg                                 >  1:981928/1‑148 (MQ=255)
                                                                                                                           agaTAATCGAGATCGGCGGGGGTTC‑GAATGCCCCCTTTTTGCAATAATGTTTGCGCCTCGCGCTGTAGCCAGCG‑CAGGCCCCCCTCGCCACCGCGCGATGCAACGTTGGTATCGCCGTTGATCGCCATCAGTAGGAGCAAGGTATCGAGCa                        <  1:689381/151‑1 (MQ=255)
                                                                                                                           agaTAATCGAGATCGGCGGGGGGTC‑GAATGCCCCCTTTTTTCAATAATGTTTTCGCCCCGCGCCGTAGCCCGCG‑CCGGCCCCCCCCGCCCCCCCCCGATGCAAAGGTGGGATCGCCGTTGATTGCCATCAGTAGGAGCAAGGGATCGAg                          >  2:575665/1‑149 (MQ=255)
                                                                                                                                aTCGAGATCGGCGGGGGTTC‑GAATGCCCCCTTTTTGCAATAATGTTTGCGCCTCGCGCTGTAGCCAGCG‑CAGGCCCCCCTCGCCACCGCGCGATGCAACGTTGGTATCGCCGTTGATCGCCATCAGTAGGAGCAAGGTATCGAGCAATGcc                   <  2:552694/151‑1 (MQ=255)
                                                                                                                                  cGAGATCGGCGGGGGTTC‑GAATGCCCCCTTTTTGCAATAATGTTTGCGCCTCGCGCTGTAGCCAGCG‑CAGGCCCCCCTCGCCACCGCGCGATGCAACGTTGGTATCGCCGTTGATCGCCATCAGTAGGAGCAAGGTATCGAGCAATGCCAg                 <  2:40766/151‑1 (MQ=255)
                                                                                                                                  cGAGATCGGCGGGGGTTC‑GAATGCCCCCTTTTTGCAATAATGTTTGCGCCTCGCGCTGTAGCCAGCG‑CAGGCCCCCCTCGCCACCGCGCGATGCAACGTTGGTATCGCCGTTGATCGCCATCAGTAGGAGCAAGGTATCGAGCAATGCCAg                 <  1:1048523/151‑1 (MQ=255)
                                                                                                                                     gaTCGGCGGGGGTTC‑GAATGCCCCCTTTTTGCAATAATGTTTGCGCCTCGCGCTGTAGCCAGCG‑CAGGCCCCCCTCGCCACCGCGCGATGCAACGTTGGTATCGCCGTTGATCGCCATCAGTAGGAGCAAGGTATCGAGCAATGCCAGTTc              <  2:748024/151‑1 (MQ=255)
                                                                                                                                        cggcggGGGTTC‑GAATGCCCCCTTTTTGCAATAATGTTTGCGCCTCGCGCTGTAGCCAGCG‑CAGGCCCCCCTCGCCACCGCGCGATGCAACGTTGGTATCGCCGTTGATCGCCATCAGTAGGAGCAAGGTATCGAGCAATGCCAGTTCAg            <  1:1228035/150‑1 (MQ=255)
                                                                                                                                          gcggGGGTTC‑GAATGCCCCCTTTTTGCAATAATGTTTGCGCCTCGCGCTGTAGCCAGCG‑CAGGCCCCCCTCGCCACCGCGCGATGCAACGTTGGTATCGCCGTTGATCGCCATCAGTAGGAGCAAGGTATCGAGCAATGCCAGTTCAGGAt         >  1:479059/1‑151 (MQ=255)
                                                                                                                                                 ttC‑GAATGCCCCCTTTTTGCAATAATGTTTGCGCCTCGCGCTGTAGCCAGCG‑CAGGCCCCCCTCGCCACCGCGCGATGCAACGTTGGTATCGCCGTTGATCGCCATCAGTAGGAGCAAGGTATCGAGCAATGCCAGTTCAGGATCTAAccc  >  2:943250/1‑151 (MQ=255)
                                                                                                                                                 |                                                                                                                                                        
ATTTCCAAATTTATCAAGTGCTTAAATAATTAAATCTGTGCTAAAAACCAGGTAAGGATCAGTAGGTCAGCACTGCCGCCTGGACTGAGATTTCGTTCGATACACTCCCTGTCGAACTGCCGGAGATAATCGAGATCGGCGGGGGTTC‑GAATGCCCCCTTTTTGCAATAATGTTTGCGCCTCGCGCTGTAGCCAGCG‑CAGGCCCCCCTCGCCACCGCGCGATGCAACGTTGGTATCGCCGTTGATCGCCATCAGTAGGAGCAAGGTATCGAGCAATGCCAGTTCAGGATCTAACCC  >  NC_000913/646602‑646897

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