Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 1,551,285 A→C 8.6% intergenic (+44/+82) 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,2850AC 8.6% 230.4 / 11.6 82intergenic (+44/+82)groL/pykAchaperonin GroEL/pyruvate kinase II
Reads supporting (aligned to +/- strand):  ref base A (39/35);  new base C (4/3);  total (43/39)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 6.74e-01

TCACCACCGAAGCCATGGTTGCCGACCTGCCGGAAGACAAGCCAGCTGCTGGCATGCCTGATATGGGTGGCATGGGTGGCATGGGCGGCATGATGTAAGCCAGCCTTACCCCCCTGCACCATCAAAAAGCCCCGGTTCGCGAGAACCGGGGCTTTTTTAT‑‑TGCTGTCGGGGCCTCTTCGCAGCACAAGGCTGCGCTTGAGCCTTGTGCTGTGAAACGAGTGGATCAACCGACCAGCGGATCACCGACGTGCAGGATCTTCATGCCATTGGTGCCACCGATGGTGTGG  >  NC_002947/1551144‑1551430
                                                                                                                                              |                                                                                                                                                   
tCACCACCGAAGCCATGGTTGCCGACCTGCCGGAAGACAAGCCAGCTGCTGGCATGCCTGATATGGGTGGCATGGGTGGCATGGGCGGCATGATGTAAGCCAGCCTTACCCCCCTGCACCATCAAAAAGCCCCGGTTCGCGAGAACCgg                                                                                                                                              <  6:58812/149‑1 (MQ=255)
tCACCACCGAAGCCATGGTTGCCGACCTGCCGGAAGACAAGCCAGCTGCTGGCATGCCTGATATGGGTGGCATGGGTGGCATGGGCGGCATGATGTAAGCCAGCCTTACCCCCCTGCACCATCAAAAAGCCCCGGTTCGCGAGAACCgg                                                                                                                                              <  1:204417/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                                                                                                                                         >  4:253737/1‑122 (MQ=255)
                                gAAGACAAGCCAGCTGCTGGCATGCCTGATATGGGTGGCATGGGTGGCATGGGCGGCATGATGTAAGCCAGCCTTACCCCCCTGCACCATCAAAAAGCCCCGGTTCGCGAGAACCGGGGCtt                                                                                                                                         <  3:253737/122‑1 (MQ=255)
                                   gACAAGCCAGCTGCTGGCATGCCTGATATGGGTGGCATGGGTGGCAAGGGCGGCATGATGTAAGCCAGCCTTACCCCCCTGCACCATCAAAAAGCCCCGGTTCGCGAGAACCGGGGCTTTTTTAT‑‑TGCTGTCGGGGCCTCTTCgcagca                                                                                                         >  5:292673/1‑149 (MQ=255)
                                      aaGCCAGCTGCTGGCATGCCTGATATGGGTGGCATGGGTGGCATGGGCGGCATGATGTAAGCCAGCCTTACCCCCCTGCACCATCAAAAAGCCCCGGTTCGCGAGAACCGGGGCTTTTTTAT‑‑TGCTGTCGGGGCCTCTTCGCAGcac                                                                                                        >  5:188242/1‑147 (MQ=255)
                                      aaGCCAGCTGCTGGCATGCCTGATATGGGTGGCATGGGTGGCATGGGCGGCATGATGTAAGCCAGCCTTACCCCCCTGCACCATCAAAAAGCCCCGGTTCGCGAGAACCGGGGCTTTTTTAT‑‑TGCTGTCGGGGCCTCTTCGCAGcac                                                                                                        >  1:267873/1‑147 (MQ=255)
                                      aaGCCAGCTGCTGGCATGCCTGATATGGGTGGCATGGGTGGCATGGGCGGCATGATGTAAGCCAGCCTTACCCCCCTGCACCATCAAAAAGCCCCGGTTCGCGAGAACCGGGGCTTTTTTAT‑‑TGCTGTCGGGGCCTCTTCGCAGcac                                                                                                        >  3:216461/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                                                                         >  5:299218/1‑149 (MQ=255)
                                                                   tggcatgggtggcatgggCGGCATGATGTAAGCCAGCCTTACCCCCCTGCACCATCAAAAAGCCCCGGTTCGCGAGAACCGGGGCTTTTTTAT‑‑TGCTGTCGGGGCCTCTTCGCAGCACAAGGCTGCGCTTGAGCCTTGTGCTGTGAAAc                                                                         <  4:331965/149‑1 (MQ=255)
                                                                       atgggtggcatgggCGGCATGATGTAAGCCAGCCTTACCCCCCTGCACCATCAAAAAGCCCCGGTTCGCGAGAACCGGGGCTTTTTTAT‑‑TGCTGTCGGGGCCTCTTCGCAGCACAAGGCTGCGCTTGAGCCTTGTGCTGTGAAACGAGt                                                                     >  3:106706/1‑149 (MQ=255)
                                                                       atgggtggcatgggCGGCATGATGTAAGCCAGCCTTACCCCCCTGCACCATCAAAAAGCCCCGGTTCGCGAGAACCGGGGCTTTTTTAT‑‑TGCTGTCGGGGCCTCTTCGCAGCACAAGGCTGCGCTTGAGCCTTGTGCTGTGAAACGAGt                                                                     >  5:317613/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                                                                     <  6:196360/117‑1 (MQ=255)
                                                                                                       ccTTACCCCCCTGCACCATCAAAAAGCCCCGGTTCGCGAGAACCGGGGCTTTTTTAT‑‑TGCTGTCGGGGCCTCTTCGCAGCACAAGGCTGCGCTTGAGCCTTGTGCTGTGAAACGAGt                                                                     >  5:196360/1‑117 (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                  <  5:257365/143‑1 (MQ=37)
                                                                                                                                 ccccGGTTCTCGCGAACCGGGGCTTTTTGATGGTGCAGGGGGGCCTCTTCGCAGCACAAGGCTGCGCTTGAGCCTTGTGCTGTGAAACGAGTGGATCAACCGACCAGCGGATCACCGACGTGCAGGATCTTCATGCCATTGGt                  >  6:257365/1‑143 (MQ=37)
                                                                                                                                 ccccGGTTCGCGAGAACCGGGGCTTTTTTAT‑‑TGCTGTCGGGGCCTCTTCGCAGCACAAGGCTGCGCTTGAGCCTTGTGCTGTGAAACGAGTGGATCAACCGACCAGCGGATCACCGACGTGCAGGATCTTCAt                           >  6:283201/1‑133 (MQ=255)
                                                                                                                                 ccccGGTTCGCGAGAACCGGGGCTTTTTTAT‑‑TGCTGTCGGGGCCTCTTCGCAGCACAAGGCTGCGCTTGAGCCTTGTGCTGTGAAACGAGTGGATCAACCGACCAGCGGATCACCGACGTGCAGGATCTTCAt                           <  5:283201/133‑1 (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)
                                                                                                                                             gagaACCGGGGCTTTTTTAT‑‑TGCTGTCGGGGCCTCTTCGCAGCACAAGGCTGCGCTTGAGCCTTGTGCTGTGAAACGAGTGGATCAACCGACCAGCGGATCACCGACGTGCAGGATCTTCATGCCATTGGTGCCACCGATGGTGTgg  >  4:274315/1‑147 (MQ=255)
                                                                                                                                              |                                                                                                                                                   
TCACCACCGAAGCCATGGTTGCCGACCTGCCGGAAGACAAGCCAGCTGCTGGCATGCCTGATATGGGTGGCATGGGTGGCATGGGCGGCATGATGTAAGCCAGCCTTACCCCCCTGCACCATCAAAAAGCCCCGGTTCGCGAGAACCGGGGCTTTTTTAT‑‑TGCTGTCGGGGCCTCTTCGCAGCACAAGGCTGCGCTTGAGCCTTGTGCTGTGAAACGAGTGGATCAACCGACCAGCGGATCACCGACGTGCAGGATCTTCATGCCATTGGTGCCACCGATGGTGTGG  >  NC_002947/1551144‑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: