Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 1,551,282 G→T 8.6% intergenic (+41/+85) groL → / ← pykA chaperonin GroEL/pyruvate kinase II

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0029471,551,2820GT 8.6% 222.9 / 9.7 80intergenic (+41/+85)groL/pykAchaperonin GroEL/pyruvate kinase II
Reads supporting (aligned to +/- strand):  ref base G (38/35);  new base T (4/3);  total (42/38)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 5.18e-01

TGATGATCACCACCGAAGCCATGGTTGCCGACCTGCCGGAAGACAAGCCAGCTGCTGGCATGCCTGATATGGGTGGCATGGGTGGCATGGGCGGCATGATGTAAGCCAGCCTTACCCCCCTGCACCATCAAAAAGCCCCGGTTCGCGAGAACCGGGGCTTTTTTAT‑‑TGCTGTCGGGGCCTCTTCGCAGCACAAGGCTGCGCTTGAGCCTTGTGCTGTGAAACGAGTGGATCAACCGACCAGCGGATCACCGACGTGCAGGATCTTCATGCCATTGGTGCCACCGATGGTGTGG  >  NC_002947/1551138‑1551430
                                                                                                                                                 |                                                                                                                                                      
tgatgtTCACCTCCGTAGCCATGGTTGCCGACCTGCCGGAAGACAAGCCAGCTGCTGGCATGCCTGATATGGGTGGCATGGGTGTCATGGGCGGCATGATGTAAGCCAGCCTTACCCCCCTGCACCATCAAAAAGCCCCGGTTctcgc                                                                                                                                                     <  3:269537/148‑2 (MQ=255)
tgatgaTCACCACCGAAGCCATGGTTGCCGACCTGCCGGAAGACAAGCCAGCTGCTGGCATGCCTGATATGGGTGGCATGGGTGGCATGGGCGGCATGATGTAAGCCAGCCTTACCCCCCTGCACCATCAAAAAGCCCCGGTTcgc                                                                                                                                                       <  1:18528/146‑1 (MQ=255)
tgatgaTCACCACCGAAGCCATGGTTGCCGACCTGCCGGAAGACAAGCCAGCTGCTGGCATGCCTGATATGGGTGGCATGGGTGGCATGGGCGGCATGATGTAAGCCAGCCTTACCCCCCTGCACCATCAAAAAGCCCCGGTTcgc                                                                                                                                                       >  2:18528/1‑146 (MQ=255)
tgatgaTCACCACCGAAGCCATGGTTGACGACCTGCCGGAAGACAAGCCAGCTGCTGGCATGCCTGATATGGGTGGCATGGGTGGCATGGGCGGCATGATGTAAGCCAGCCTTACCCCCCTGCACCATCAAAAAGCCACGCTTagc                                                                                                                                                       >  8:170605/1‑146 (MQ=255)
      tCACCACCGAAGCCATGGTTGCCGACCTGCCGGAAGACAAGCCAGCTGCTGGCATGCCTGATATGGGTGGCATGGGTGGCATGGGCGGCATGATGTAAGCCAGCCTTACCCCCCTGCACCATCAAAAAGCCCCGGTTCGCGAGAACCgg                                                                                                                                              <  1:204417/149‑1 (MQ=255)
      tCACCACCGAAGCCATGGTTGCCGACCTGCCGGAAGACAAGCCAGCTGCTGGCATGCCTGATATGGGTGGCATGGGTGGCATGGGCGGCATGATGTAAGCCAGCCTTACCCCCCTGCACCATCAAAAAGCCCCGGTTCGCGAGAACCgg                                                                                                                                              <  6:58812/149‑1 (MQ=255)
        accgccGAAGCCATGGTTGCCGACCTGCCGGAAGACAAGCCAGCTGCTGGCATGCCTGATATGGGTGGCATGGGTGGCATGGGCGGCATGATGTAAGCCAGCCTTACCCCCCTGCACCATCAAAAAGCCCCGGTTCTCGCGAACCgggg                                                                                                                                            >  4:265008/1‑149 (MQ=255)
           accGAAGCCATGGTTACCGACCTGCCGGAAGACAAGCCAGCTGCTGGCATGCCTGATATGGGTGGCATGGGTGGCATGGGCGGCATGATGTAAGCCAGCCTTACCCCCCTGCACCATCAAAAAGCCCCGGTTCGCGAGAACCGGGGCtt                                                                                                                                         <  6:178295/149‑1 (MQ=255)
            ccGAAGCCATGGTTGCCGACCTGCCGGAAGACAAGCCAGCTGCTGGCATGCCTGATATGGGTGGCATGGGTGGCATGGGCGGCATGATGTAAGCCAGCCTTACCCCCCTGCACCATCAAAAAGCCCCGGTTCGCGAGAACCGGGGCttt                                                                                                                                        >  6:272145/1‑149 (MQ=255)
             cGAAGCCATGGTTGCCGACCTGCCGGAAGACAAGCCAGCTGCTGGCATGCCTGATATGGGTGGCATGGGTGGCATGGGCGGCATGATGTAAGCCAGCCTTACCCCCCTGCACCATCAAAAAGCCCCGGTTCTCGCGAACCGGGGCtttt                                                                                                                                       <  8:15437/149‑1 (MQ=255)
                aGCCATGGTTGCCGACCTGCCGGAAGACAAGCCAGCTGCTGGCATGCCTGATATGGGTGGCATGGGTGGCATGGGCGGCATGATGTAAGCCAGCCTTACCCCCCTGCACCATCAAAAAGCCCCGGTTCTCGCGAACCGGGGCtttttga                                                                                                                                    <  3:265008/149‑3 (MQ=255)
                    aTGGTTGCCGACCTGCCGGAAGACAAGCCAGCTGCTGGCATGCCTGATATGGGTGGCATGGGTGGCATGGGCGGCATGATGTAAGCCAGCCTTACCCCCCTGCACCATCAAAAAGCCCACGGTTCGCGAGAACCGGGGCTTTTTTAT‑‑Tg                                                                                                                               <  1:222938/149‑1 (MQ=255)
                    aTGGTTGCCGACCTGCCGGAAGACAAGCCAGCTGCTGGCATGCCTGATATGGGTGGCATGGGTGGCATGGGCGGCATGATGTAAGCCAGCCTTACCCCCCTGCACCATCAAAAAGCCCAGGTTCGCGAGAACCGGGGCTTTTTTAT‑‑TGc                                                                                                                              <  7:95977/149‑1 (MQ=255)
                          gCCGACCTGCCGGAAGACAAGCCAGCTGCTGGCATGCCTGATATGGGTGGCATGGGTGGCATGGGCGGCATGATGTAAGCCAGCCTTACCCCCCTGCACCATCAAAAAGCCCCGGTTCGCGAGAACCGGGGCTTTTTTAT‑‑TGCTGTCgg                                                                                                                        >  2:198508/1‑149 (MQ=255)
                          gCCGACCTGCCGGAAGACAAGCCAGCTGCTGGCATGCCTGATATGGGTGGCATGGGTGGCATGGGCGGCATGATGTAAGCCAGCCTTACCCCCCTGCACCATCAAAAAGCCCCGGTTCGCGAGAACCGGGGCTTTTTTAT‑‑TGCTGTCg                                                                                                                         <  8:46591/148‑1 (MQ=255)
                           ccGACCTGCCGGAAGACAAGCCAGCTGCTGGCATGCCTGATATGGGTGGCATGGGTGGCATGGGCGGCATGATGTAAGCCAGCCTTACCCCCCTGCACCATCAAAAAGCCCCGGTTCGCGAGAACCGGGGCTTTTTTAT‑‑TGCTGTCggg                                                                                                                       >  3:171035/1‑149 (MQ=255)
                             gACCTGCCGGAAGACAAGCCAGCTGCTGGCATGCCTGATATGGGTGGCATGGGTGGCATGGGCGGCATGATGTAAGCCAGCCTTACCCCCCTGCACCATCAAAAAGCCCCGGTTCGCGAGAACCGGGGGTTTTTTAT‑‑TGCTGTCGGGGc                                                                                                                     <  3:282596/149‑1 (MQ=255)
                             gACCTGCCGGAAGACAAGCCAGCTGCTGGCATGCCTGATATGGGTGGCATGGGTGGCATGGGCGGCATGATGTAAGCCAGCCTTACCCCCCTGCACCATCAAAAAGCCCCGGTTCGCGAGAACCGGGGCtttt                                                                                                                                       <  2:248033/133‑1 (MQ=255)
                             gACCTGCCGGAAGACAAGCCAGCTGCTGGCATGCCTGATATGGGTGGCATGGGTGGCATGGGCGGCATGAGGTAAGCCAGCCTTACCCCCCTGCACCATCAAAGAGCCCCGGTTCGCGAGAACCGGGGCTTTTTTAT‑‑TGCTGTCGGGGc                                                                                                                     <  5:332047/149‑1 (MQ=255)
                             gACCTGCCGGAAGACAAGCCAGCTGCTGGCATGCCTGATAAGGGTGGCATGGGTGGCATGGGCGGCATGATGTAAGCCAGCCTTACCCCCCTGCACCATCAAAAAGCCCCGGTTCGCGAGAACCGGGGCtttt                                                                                                                                       >  1:248033/1‑133 (MQ=255)
                                cTGCCGGAAGACAAGCCAGCTGCTGGCATGCCTGATATGGGTGGCATGGGTGGCATGGGCGGCATGATGTAAGCCAGCCTTACCCCCCTGCACCATCAAAAAGCCCCGGTTCGCGAGAACCGGGGCTTTTTTAT‑‑TGCTGTCGGGGCctc                                                                                                                  >  1:149149/1‑149 (MQ=255)
                                   ccGGAAGACAAGCCAGCTGCTGGCATGCCTGATATGGGTGGTATGGGTGGCATGGGCGGCATGATGTAAGCCAGCCTTACCCCCCTGCACCATCAAAAAGCCCCGGTTCGCGAGAACCGGGGCTTTTTTAT‑‑TGCTGTCGGGGCCTCtt                                                                                                                <  7:181067/148‑1 (MQ=255)
                                     ggAAGACAAGCCAGCTGCTGGCATGCCTGATATGGGTGGCATGGGTGGCATGGGCGGCATGATGTAAGCCAGCCTTACCCCCCTGCACCATCAAAAAGCCCCGGTTCGCGAGAACCGGGGCTTTTTTAT‑‑TGCTGTCGGGGCCTCTTCg                                                                                                              >  3:77835/1‑148 (MQ=255)
                                      gAAGACAAGCCAGCTGCTGGCATGCCTGATATGGGTGGCATGGGTGGCATGGGCGGCATGATGTAAGCCAGCCTTACCCCCCTGCACCATCAAAAAGCCCCGGTTCGCGAGAACCGGGGCtt                                                                                                                                         <  3:253737/122‑1 (MQ=255)
                                      gAAGACAAGCCAGCTGCTGGCATGCCTGATATGGGTGGCATGGGTGGCATGGGCGGCATGATGTAAGCCAGCCTTACCCCCCTGCACCATCAAAAAGCCCCGGTTCGCGAGAACCGGGGCtt                                                                                                                                         >  4:253737/1‑122 (MQ=255)
                                         gACAAGCCAGCTGCTGGCATGCCTGATATGGGTGGCATGGGTGGCAAGGGCGGCATGATGTAAGCCAGCCTTACCCCCCTGCACCATCAAAAAGCCCCGGTTCGCGAGAACCGGGGCTTTTTTAT‑‑TGCTGTCGGGGCCTCTTCgcagca                                                                                                         >  5:292673/1‑149 (MQ=255)
                                            aaGCCAGCTGCTGGCATGCCTGATATGGGTGGCATGGGTGGCATGGGCGGCATGATGTAAGCCAGCCTTACCCCCCTGCACCATCAAAAAGCCCCGGTTCGCGAGAACCGGGGCTTTTTTAT‑‑TGCTGTCGGGGCCTCTTCGCAGcac                                                                                                        >  3:216461/1‑147 (MQ=255)
                                            aaGCCAGCTGCTGGCATGCCTGATATGGGTGGCATGGGTGGCATGGGCGGCATGATGTAAGCCAGCCTTACCCCCCTGCACCATCAAAAAGCCCCGGTTCGCGAGAACCGGGGCTTTTTTAT‑‑TGCTGTCGGGGCCTCTTCGCAGcac                                                                                                        >  5:188242/1‑147 (MQ=255)
                                            aaGCCAGCTGCTGGCATGCCTGATATGGGTGGCATGGGTGGCATGGGCGGCATGATGTAAGCCAGCCTTACCCCCCTGCACCATCAAAAAGCCCCGGTTCGCGAGAACCGGGGCTTTTTTAT‑‑TGCTGTCGGGGCCTCTTCGCAGcac                                                                                                        >  1:267873/1‑147 (MQ=255)
                                              gCCAGCTGCTGGCATGCCTGATATGGGTGGCATGGGTGGCATGGGCGGCATGATGTAAGCCAGCCTTACCCCCCTGCACCATCAAAAAGCCCCGGTTCGCGAGAACCGGGGCTTTTTTAT‑‑TGCTGTCGGGGCCTCTTCGCAGCACAAgg                                                                                                    >  5:223438/1‑149 (MQ=255)
                                                cAGCTGCTGGCATGCCTGATATGGGTGGCATGGGTGGCATGGGCGGCATGATGTAAGCCAGCCTTACCCCCCTGCACCATCAAAAAGCCCCGGTTCGCGAGAACCGGGGCTTTTTTAT‑‑TGCTGTCGGGGCCTCTTCGCAGCACAAGGCt                                                                                                  >  2:128692/1‑149 (MQ=255)
                                                 aGCTGCTGGCATGCCTGATATGGGTGGCATGGGTGGCAAGGGCGGCATGATGTAAGCCAGCCTTACCCCCCTGCACCATCAAAAAGCCCCGGTTCGCGAGAACCGGGGCTTTTTTAT‑‑TGCTGTCGGGGCCTCTTCGCAGCACAAGGCTg                                                                                                 <  6:292673/149‑1 (MQ=255)
                                                    tgctgGCATGCCTGATATGGGTGGCATGGGTGGCATGGGCGGCATGATGTAAGCCAGCCTTACCCCCCTGCACCATCAAAAAGCCCCGGTTCGCGAGAACCGGGGCTTTTTTAT‑‑TGCTGTCGGGGCCTCTTCGCAGCACAAGGCTgcgc                                                                                              >  3:241311/1‑149 (MQ=255)
                                                      ctgGCATGCCTGATATGGGTGGCATGGGTGGCATGGGCGGCATGATGTAAGCCAGCCTTACCCCCCTGCACCATCAAAAAGCCCCGGTTCGCGAGAACCGGGGCTTTTTTAT‑‑TGCTGTCGGGGCCTCTTCGCAGCACAAGGCTgcgc                                                                                              <  7:170605/147‑1 (MQ=255)
                                                            tGCCTGATATGGGTGGCATGGGTGGCATGGGCGGCATGATGTAAGCCAGCCTTACCCCCCCGCACCATCAAAAAGCCCCGGTTCGCGAGAACCGGGGCTTTTTTAT‑‑CGCTGTCGGGGCCTCTTCGCAGCACAAGGCTCCGCTTGAGCCt                                                                                      >  8:83630/1‑149 (MQ=255)
                                                              ccTGATATGGGTGGCATGGGTGGCATGGGCGGCATGATGTAAGCCAGCCTTACCCCCCTGCACCATCAAAAAGCCCCGGTTCGCGAGAACCGGGGCTTTTTTAT‑‑TGCTGTCGGGGCCTCTTCGCAGCACAAgg                                                                                                    <  4:81956/133‑1 (MQ=255)
                                                              ccTGATATGGGTGGCATGGGTGGCATGGGCGGCATGATGTAAGCCAGCCTTACCCCCCTGCACCATCAAAAAGCCCCGGTTCGCGAGAACCGGGGCTTTTTTAT‑‑TGCTGTCGGGGCCTCTTCGCAGCACAAGGCTGCGCTTGAGCCTtg                                                                                    <  8:48571/149‑1 (MQ=255)
                                                                tGATATGGGTGGCATGGGTGGCATGGGCGGCATGATGTAAGCCAGCCTTACCCCCCTGCACCATCAAAAAGCCCCGGTTCGCGAGAACCGGGGCTTTTTTAT‑‑TGCTGTCGGGGCCTCTTCGCAGCACAAGGCTGCGCTTGAGCCTtgtg                                                                                  >  2:13200/1‑149 (MQ=255)
                                                                  atatGGGTGGCATGGGTGGCATGGGCGGCATGATGTAAGCCAGCCTTACCCCCCTGCACCATCAAAAAGCCCCGGTTCGCGAGAACCGGGGCTTTTTTAT‑‑TGCTGTCGGGGCCTCTTCGCAGCACAAGGCTGCGCTTGAGCCTtgtg                                                                                  <  4:241311/147‑1 (MQ=255)
                                                                   tatGGGTGGCATGGGTGGCATGGGCGGCATGATGTAAGCCAGCCTTACCCCCCTGCACCATAAAAAAGCCCCCGTTCGCGAGAACCCGGGCTTTTTTAT‑‑TTCTGTCGGGGCCTATTAGCCGCACAAGGCTg                                                                                                 >  2:32209/1‑131 (MQ=255)
                                                                    atgggtggcatgggtggcATGGGCGGAATGATGTAAGCCAGCCTTACCCCCCTGCACCATCAAAAAGCCCCGGTTCGGGAGAACCGGGGGTTTTTTAT‑‑TGCTGTCGGGGCCTCTTCGCAGCACTAGGCTGGGCTTGAGGGTTGTGCtgt                                                                              >  1:195992/1‑149 (MQ=255)
                                                                         tggcatgggtggcatgggCGGCATGATGTAAGCCAGCCTTACCCCCCTGCACCATCAAAAAGCCCCGGTTCGCGAGAACCGGGGCTTTTTTAT‑‑TGCTGTCGGGGCCTCTTCGCAGCACAAGGCTGCGCTTGAGCCTTGTGCTGTGAAAc                                                                         <  4:331965/149‑1 (MQ=255)
                                                                         tggcatgggtggcatgggCGGCATGATGTAAGCCAGCCTTACCCCCCTGCACCATCAAAAAGCCCCGGTTCGCGAGAACCGGGGCTTTTTTAT‑‑TGCTGTCGGGGCCTCTTCGCAGCACAAGGCTGCGCTTGAGCCTTGTGCTGTGAAAc                                                                         >  5:299218/1‑149 (MQ=255)
                                                                             atgggtggcatgggCGGCATGATGTAAGCCAGCCTTACCCCCCTGCACCATCAAAAAGCCCCGGTTCGCGAGAACCGGGGCTTTTTTAT‑‑TGCTGTCGGGGCCTCTTCGCAGCACAAGGCTGCGCTTGAGCCTTGTGCTGTGAAACGAGt                                                                     >  5:317613/1‑149 (MQ=255)
                                                                             atgggtggcatgggCGGCATGATGTAAGCCAGCCTTACCCCCCTGCACCATCAAAAAGCCCCGGTTCGCGAGAACCGGGGCTTTTTTAT‑‑TGCTGTCGGGGCCTCTTCGCAGCACAAGGCTGCGCTTGAGCCTTGTGCTGTGAAACGAGt                                                                     >  3:106706/1‑149 (MQ=255)
                                                                              tgggtggcatgggCGGCATGATGTAAGCCAGCCTTACCCCCCTGCACCATCAAAAAGCCCCGGTTCGCGAGAACCGGGGCTTTTTTAT‑‑TGCTGTCGGGGCCTCTTCGCAGCACAAGGCTGCGCTTGAGCCTTGTGCTGTGAAACGAGTg                                                                    >  6:112946/1‑149 (MQ=255)
                                                                                ggtggcatgggCGGCATGATGTAAGCCAGCCTTACCCCCCTGCACCATCAAAAAGCCCCGGTTCGCGAGAACCGGGGCTTTTTTAT‑‑TGCTGTCGGGGCCTCTTCGCAGCACAAGGCTGCGCTTGAGCCTTGTGCTGTGAAACGAGTGGa                                                                  <  1:198508/149‑1 (MQ=255)
                                                                                    gcatgggCGGCATGATGTAAGCCAGCCTTACCCCCCTGCACCATCAAAAAGCCCCGGTTCGCGAGAACCGGGGCTTTTTTAT‑‑TGCTGTCGGGGCCTCTTCGCAGCACAAGGCTGCGCTTGAGCCTTGTGCTGTGAAACGAGTGGATCa                                                               <  7:92014/148‑1 (MQ=255)
                                                                                      atgggCGGCATGATGTAAGCCAGCCTTACCCCCCTGCACCATCAAAAAGCCCCGGTTCGCGAGAACCGGGGCTTTTTTAT‑‑TGCTGTCGGGGCCTCTTCGCAGCACAAGGCTGCGCTTGAGCCTTGTGCTGTGAAACGAGTGGATCaa                                                              <  8:200357/147‑1 (MQ=255)
                                                                                      atgggCGGCATGATGTAAGCCAGCCTTACCCCCCTGCACCATCAAAAAGCCCCGGTTCGCGAGAACCGGGGCTTTTTTAT‑‑TGCTGTCGGGGCCTCTTCGCAGCACAAGGCTGCGCTTGAGCCTTGTGCTGTGAAACGAGTGGATCaa                                                              >  7:200357/1‑147 (MQ=255)
                                                                                      atgggCGGCATGATGTAAGCCAGCCTTACCCCCCTGCACCATCAAAAAGCCCCGGTTCGCGAGAACCGGGGCTTTTTTAT‑‑TGCTGTCGGGGCCTCTTCGCAGCACAAGGCTGCGCTTGAGCCTTGTGCTGTGAAACGAGTGGATCAAcc                                                            >  8:328885/1‑149 (MQ=255)
                                                                                          gcggcATGATGTAAGCCAGCCTTACCCCCCTGCACCATCAAAAAGCCCCGGTTCGCGAGAACCGGGGCTTTTTTAT‑‑TGCTGTCGGGGCCTCTTCGCAGCACAAGGCTGCGCTTGAGCCTTGTGCTGTGAAACGAGTGGATCAACCGAcc                                                        <  1:13200/149‑1 (MQ=255)
                                                                                             gcATGATGTAAGCCAGCCTTACCCCCCTGCACCATCAAAAAGCCCCGGTTCGCGAGAACCGGGGCTTTTTTAT‑‑TGCTGTCGGGGCCTCTTCGCAGCACAAGGCTGCGCTTGAGCCTTGTGCTGTGAAACGAGTGGATCAACCGACCAGc                                                     <  4:77835/149‑1 (MQ=255)
                                                                                             gcATGATCAAAGCCAGCCTTACCCCCCTGCACCATCAAAAACCCCCGGTTCGCGACAACCGGGGCTTTTTTAT‑‑TGCTGTCGGGGCCTCTTCGCAGCACAAGGCAGCGCTTGAGCCTTGTGCTGTGAAACGAGTGGATCAACCGACCAGc                                                     <  2:195992/149‑1 (MQ=255)
                                                                                               atgatgTAAGCCAGCCTTACCCCCCTGCACCATCAAAAAGCCCCGGTTCGCGAGAACCGGGGCTTTTTTAT‑‑TGCTGTCGGGGCCTCTTCGCAGCACAAGGCTGCGCTTGAGCCTTGTGCTGTGAAACGAGTGGATCAACCGACCAGCg                                                    >  8:133885/1‑148 (MQ=255)
                                                                                                tgatgTAAGCCAGCCTTACCCCCCTGCACCATCAAAAAGCCCCGGTTCGCGAGAACCGGGGCTTTTTTAT‑‑TGCTGTCGGGGCCTCTTCGCAGCACAAGGCTGCGCTTGAGCCTTGTGCTGTGAAACGAGTGGATCAACCGACCAGCGGa                                                  <  6:188242/149‑1 (MQ=255)
                                                                                                 gatgTAAGCCAGCCTTACCCCCCTGCACCATCAAAAAGCCCCGGTTCGCGAGAACCGGGGCTTTTTTAT‑‑TGCTGTCGGGGCCTCTTCGCAGCACAAGGCTGCGCTTGAGCCTTGTGCTGTGAAACGAGTGGATCAACCGACCAGCGGAt                                                 >  1:149731/1‑149 (MQ=255)
                                                                                                   tgTAAGCCAGCCTTACCCCCCTGCACCATCAAAAAGCCCCGGTTCGCGAGAACCGGGGCTTTTTTAT‑‑TGCTGTCGGGGCCTCTTCGCAGCACAAGGCTGCGCTTGAGCCTTGTGCTGTGAAACGAGTGGATCAACCGACCAGCGGATCa                                               <  5:108030/149‑1 (MQ=255)
                                                                                                     tAAGCCAGCCTTACCCCCCTGCACCATCAAAAAGCCCCGGTTCGCGAGAACCGGGGCTTTTTTAT‑‑TGCTGTCGGGGCCTCTTCGCAGCACAAGGCTGCGCTTGAGCCTTGTGCTGTGAAACGAGTGGATCAACCGACCAGCGGATCAcc                                             >  1:74386/1‑149 (MQ=255)
                                                                                                     tAAGCCAGCCTTACCCCCCTGCACCATCAAAAAGCCCCGGTTCGCGAGAACCGGGGCTTTTTTAT‑‑TGCTGTCGGGGCCTCTTCGCAGCACAAGGCTGCGCTTGAGCCTTGTGCTGTGAAACGAGTGGATCAACCGACCAGCGGATCAc                                              <  7:328885/148‑1 (MQ=255)
                                                                                                      aaGCCAGCCTTACCCCCCTGCACCATCAAAAAGCCCCGGTTCGCGAGAACCGGGGCTTTTTTAT‑‑TGCTGTCGGGGCCTCTTCGCAGCACAAGGCTGCGCTTGAGCCTTGTGCTGTGAAACGAGTGGATCAACCGACCAGCGGATCAcc                                             <  4:106706/148‑1 (MQ=255)
                                                                                                      aaGCCAGCCTTACCCCCCTGCACCATCAAAAAGCCCCGGTTCGCGAGAACCGGGGCTTTTTTAT‑‑TGCTGTCGGGGCCTCTTCGCAGCACAAGGCTGCGCTTGAGCCTTGTGCTGTGAAACGAGTGGATCAACCGACCAGCGGATCACCg                                            >  3:167258/1‑149 (MQ=255)
                                                                                                        gccagccTTACCCCCCTGCACCATCAAAAAGCCCCGGTTCGCGAGAACCGGGGCTTTTTTAT‑‑TGCTGTCGGGGCCTCTTCGCAGCACAAGGCTGCGCTTGAGCCTTGTGCTGTGAAACGAGTGGATCAACCGACCAGCGGATCACCGAc                                          >  2:12431/1‑149 (MQ=255)
                                                                                                             ccTTACCCCCCTGCACCATCAAAAAGCCCCGGTTCGCGAGAACCGGGGCTTTTTTAT‑‑TGCTGTCGGGGCCTCTTCGCAGCACAAGGCTGCGCTTGAGCCTTGTGCTGTGAAACGAGt                                                                     >  5:196360/1‑117 (MQ=255)
                                                                                                             ccTTACCCCCCTGCACCATCAAAAAGCCCCGGTTCGCGAGAACCGGGGCTTTTTTAT‑‑TGCTGTCGGGGCCTCTTCGCAGCACAAGGCTGCGCTTGAGCCTTGTGCTGTGAAACGAGt                                                                     <  6:196360/117‑1 (MQ=255)
                                                                                                                     cccTGCACCATCAAAAAGCCCCGGTTCGCGAGAACCGGGGCTTTTTTAT‑‑TGCTGTCGCGGCCTCTTCGCAGCACAAGGCTGCGCTTGAGCCTTGTGCTGTGAAACTAGTGGATCAAACGACCATCGGATCACCGACGTGCAGGATAtt                              >  7:351654/1‑148 (MQ=255)
                                                                                                                       cTGCACCATCAAAAAGCCCCGGTTCGCGAGAACCGGGGCTTTTTTAT‑‑TGCTGTCGGGGCCTCTTCGCAGCACAAGGCTGCGCTTGAGCCTTGTGCTGTGAAACGAGTGGATCAACCGACCAGCGGATCACCGACGTGCAGGATCTTCAt                           <  6:299218/149‑1 (MQ=255)
                                                                                                                               tCAAAAAGCCCCGGTTCGCGAGAACCGGGGCTTTTTTAT‑‑TGCTGTCGGGGCCTCTTCGCAGCACAAGGCTGCGCTTGAGCCTTGTGCTGTGAAACGAGTGGATCAACCGACCAGCGGATCACCGACGTGCAGGATCTTCATGCCATTgg                   >  3:94723/1‑149 (MQ=255)
                                                                                                                                  aaaaGCCCCGGTTCGCGAGAACCGGGGCTTTTTTAT‑‑TGCTGTCGGGGCCTCTTCGCAGCACAAGGCTGCGCTTGAGCCTTGTGCTGTGAAACGAGTGGATCAACCGACCAGCGGATCACCGACGTGCAGGATCTTCATGCCATTGGTGc                <  5:112946/149‑1 (MQ=255)
                                                                                                                                  aaaaGCCCCGGTTCGCGAGAACCGGGGCTTTTTTAT‑‑TGCTGTCGGGGCCTCTTCGCAGCACAAGGCTGCGCTTGAGCCTTGTGCTGTGAAACGAGTGGATCAACCGACCAGCGGATCACCGACGTGCAGGATCTTCATGCCATTGGTGc                >  2:239735/1‑149 (MQ=255)
                                                                                                                                   aaaGCCCCGGTTCGCGAGAACCGGGGCTTTTTTAT‑‑TGCTGTCGGGGCCTCTTCGCAGCACAAGGCTGCGCTTGAGCCTTGTGCTGTGAAACGAGTGGATCAACCGACCAGCGGATCACCGACGTGCAGGATCTTCATGCCATTGGTGcc               <  5:207155/149‑1 (MQ=255)
                                                                                                                                     aGCCCCGGTTCGCGAGAACCGGGGCTTTTTTAT‑‑TGCTGTCGGGGCCTCTTCGCAGCACAAGGCTGCGCTTGAGCCTTGTGCTGTGAAACGAGTGGATCAACCGACCAGCGGATCACCGACGTGCAGGATCTTCATGCCATTGGTGCCAc             >  6:245542/1‑149 (MQ=255)
                                                                                                                                      gCCCCGGTTCGCGAGAACCGGGGCTTTTTTAT‑‑TGCTGTCGGGGCCTCTTCGCAGCACAAGGCTGCGCTTGAGCCTTGTGCTGTGAAACGAGTGGATCAACCGACCAGCGGATCACCGACGTGCAGGATCTTCATGCCATTGGTGCCAcc            <  4:171035/149‑1 (MQ=255)
                                                                                                                                       ccccGGTTCTCGCGAACCGGGGCTTTTTGATGGTGCAGGGGGGCCTCTTCGCAGCACAAGGCTGCGCTTGAGCCTTGTGCTGTGAAACGAGTGGATCAACCGACCAGCGGATCACCGACGTGCAGGATCTTCATGCCATTGGt                  >  6:257365/1‑143 (MQ=37)
                                                                                                                                       ccccGGTTCTCGCGAACCGGGGCTTTTTGATGGTGCAGGGGGGCCTCTTCGCAGCACAAGGCTGCGCTTGAGCCTTGTGCTGTGAAACGAGTGGATCAACCGACCAGCGGATCACCGACGTGCAGGATCTTCATGCCATTGGt                  <  5:257365/143‑1 (MQ=37)
                                                                                                                                       ccccGGTTCGCGAGAACCGGGGCTTTTTTAT‑‑TGCTGTCGGGGCCTCTTCGCAGCACAAGGCTGCGCTTGAGCCTTGTGCTGTGAAACGAGTGGATCAACCGACCAGCGGATCACCGACGTGCAGGATCTTCAt                           <  5:283201/133‑1 (MQ=255)
                                                                                                                                       ccccGGTTCGCGAGAACCGGGGCTTTTTTAT‑‑TGCTGTCGGGGCCTCTTCGCAGCACAAGGCTGCGCTTGAGCCTTGTGCTGTGAAACGAGTGGATCAACCGACCAGCGGATCACCGACGTGCAGGATCTTCAt                           >  6:283201/1‑133 (MQ=255)
                                                                                                                                       ccccGGTTCGCGAGAACCGGGGCTTTTTTAT‑‑TGCTGTCGGGGCCTCTTCGCAGCACAAGGCTGCGCTTGAGCCTTGTGCTGTGAAACGAGTGGATCAACCGACCAGCGGATCACCGACGTGCAGGATCTTCATGCCATTGGTGCCACCg           >  2:245172/1‑149 (MQ=255)
                                                                                                                                        cccGTTTCTCGCGAACCGGGGCTTTTTGA‑‑TGGTGCAGAGGCCTCTTCGCAGCACAAGGCTGCGCTTGAGCCTTGTGCTGTGAAACGAGTGGATCAACCGACCAGCGGATCACCGACGTGCAGGATCTTCATGCCATTGGTGCCACCGAt         >  1:67278/1‑149 (MQ=255)
                                                                                                                                        cccGGTTCTCGCGAACCGGGGCTTTTTTAT‑‑TGCTGTCGGGGCCTCTTCGCAGCACAAGGCTGCGCTTGAGCCTTGTGCTGTGAAACGAGTGGATCAACCGACCAGCGGATCACCGACGTGCAGGATCTTCATGCCATTGGTGCCACCGa          >  5:301503/1‑149 (MQ=255)
                                                                                                                                        cccGGTTCGCGAGAACCGGGGCTTTTTTAT‑‑TGCTGTCGGGGCCTCTTCGCAGCACAAGGCTGCGCTTGAGCCTTGTGCTGTGAAACGAGTGGATCAACCGACCAGCGGATCACCGACGTGCAGGATCTTCATGCCATTGGTGCCACCGa          >  2:315684/1‑149 (MQ=255)
                                                                                                                                           cGGTTCGCGAGAACCGGGGCTTTTTTAT‑‑TGCTGTCGGGGCCTCTTCGCAGCACAAGGCTGCGCTTGAGCCTTGTGCTGTGAAACGAGTGGATCAACCGACCAGCGGATCACCGACGTGCAGGATCTTCATGCCATTGGTGCCACCGATg        >  4:277185/1‑149 (MQ=255)
                                                                                                                                             gTTCGCGAGAACCGGGGCTTTTTTAT‑‑TGCTGTCGGGGCCTCTTCGCAGCACAAGGCTGCGCTTGAGCCTTGTGCTGTGAAACGAGTGGATCAACCGACCAGCGGATCACCGACGTGCAGGATCTTCATGCCATTGGTGCCACCGATGgt      <  5:272145/149‑1 (MQ=255)
                                                                                                                                              ttCGCGTTATCCGGGGCTTTTTTAT‑‑TGCTGTCGGGGCCTCTTCGCAGCACAAGGCTGCGCTTGAGCCTTGTGCTGTGAAACGAGTGGATCAACCGACCAGCGGATCACCGACGTGCAGGATCTTCATGCCATCGGTGGCACCAATggca     <  1:128692/149‑3 (MQ=255)
                                                                                                                                                cgcgAGAACCGGGGCTTTTTTAT‑‑TGCTGTCGGGGCCTCTTCGCAGCACAAGGCTGCGCTTGAGCCTTGTGCTGTGAAACGAGTGGATCAACCGACCAGCGGATCACCGACGTGCAGGATCTTCATGCCATCGGTGGCACCAATggcat    <  6:317613/148‑4 (MQ=255)
                                                                                                                                                 gcgAGAACCGGGGCTTTTTTAT‑‑TGCTGTCGGGGCCTCTTCGCAGCACAAGGCTGCGCTTGAGCCTTGTGCTGTGAAACGAGTGGATCAACCGACCAGCGGATCACCGACGTGCAGGATCTTCATGCCATTGGTGCCACCGATGGTGTgg  >  2:143074/1‑149 (MQ=255)
                                                                                                                                                 |                                                                                                                                                      
TGATGATCACCACCGAAGCCATGGTTGCCGACCTGCCGGAAGACAAGCCAGCTGCTGGCATGCCTGATATGGGTGGCATGGGTGGCATGGGCGGCATGATGTAAGCCAGCCTTACCCCCCTGCACCATCAAAAAGCCCCGGTTCGCGAGAACCGGGGCTTTTTTAT‑‑TGCTGTCGGGGCCTCTTCGCAGCACAAGGCTGCGCTTGAGCCTTGTGCTGTGAAACGAGTGGATCAACCGACCAGCGGATCACCGACGTGCAGGATCTTCATGCCATTGGTGCCACCGATGGTGTGG  >  NC_002947/1551138‑1551430

Alignment Legend
Aligned base mismatch/match (shaded by quality score): ATCG/ATCG < 3 ≤ ATCG/ATCG < 14 ≤ ATCG/ATCG < 21 ≤ ATCG/ATCG < 32 ≤ ATCG/ATCG < 36 ≤ ATCG/ATCG
Unaligned base: atcg    Masked matching base: atcg    Alignment gap:     Deleted base: