Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 1,551,282 G→T 8.8% 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.8% 205.1 / 9.4 68intergenic (+41/+85)groL/pykAchaperonin GroEL/pyruvate kinase II
Reads supporting (aligned to +/- strand):  ref base G (31/31);  new base T (2/4);  total (33/35)
Fisher's exact test for biased strand distribution p-value = 6.74e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 2.22e-01

TGATGATCACCACCGAAGCCATGGTTGCCGACCTGCCGGAAGACAAGCCAGCTGCTGGCATGCCTGATATGGGTGGCATGGGTGGCATGGGCGGCATGATGTAAGCCAGCCTTACCCCCCTGCACCATCAAAAAGCCCCGGTTCGCGAGAACCGGGGCTTTTTTATTGCTGTCGGGGCC‑T‑‑CTTCGCAGCACAAGGCTGCGCTTGAGCCTTGTGCTGTGAAACGAGTGGATCAACCGACCAGCGGATCACCGACGTGCAGGATCTTCATGCCATTGGTGCCACCGATGGTGTGG  >  NC_002947/1551138‑1551430
                                                                                                                                                |                                                                                                                                                       
tgatgaTCACCACCGAAGCCATGGTTGCCGACCTGCCGGAAGACAAGCCAGCTGCTGGCATGCCTGATATGGGTGGCATGGGTGGCATGGGCGGCATGATGTAAGCCAGCCTTACCCCCCTGCACCATCAAAAAGCCCCGGTTcgc                                                                                                                                                        >  5:283092/1‑146 (MQ=255)
  atgatCACCACCGAAGCCATGGTTGCCGACCTGCCGGAAGACAAGCCAGCTGCTGGCATGCCTGATATGGGTGGCATGGGTGGCATGGGCGGCATGATGTAAGCCAGCCTTACCCCCCTGCACCATCAAAAAGCCCCGGTTCGCGAGaa                                                                                                                                                   >  6:205505/1‑149 (MQ=255)
    gatCACCACCGAAGCCATGGTTGCCGACCTGCCGGAAGACAAGCCAGCTGCTGGCATGCCTGATATGGGTGGCATGGGTGGCATGGGCGGCATGATGTAAGCCAGCCTTACCCCCCTGCACCATCAAAAAGCCCCGGTTCGCGAGAAc                                                                                                                                                  <  3:19702/148‑1 (MQ=255)
      tCACCACCGAAGCCATGGTTGCCGACCTGCCGGAAGACAAGCCAGCTGCTGGCATGCCTGATATGGGTGGCATGGGTGGCATGGGCGGCATGATGTAAGCCAGCCTTACCCCCCTGCACCATCAAAAAGCCCCGGTTCGCGAGAACCgg                                                                                                                                               <  1:123247/149‑1 (MQ=255)
        accaccGAAGCCATGGTTGCCGACCTGCCGGAAGACAAGCCAGCTGCTGGCATGCCTGATATGGGTGGCATGGGTGGCATGGGCGGCATGATGTAAGCCAGCCTTACCCCCCTGCACCATCAAAAAGCCCCGGTTCGCGAGAACCgggg                                                                                                                                             >  7:138317/1‑149 (MQ=255)
        accaccGAAGCCATGGTTGCCGACCTGCCGGAAGACAAGCCAGCTGCTGGCATGCCTGATATGGGTGGCATGGGGGGCATGGGCGGCATGATGTAAGCCAGCCTTACCCCCCTGCACCATCAAAAAGCCCCGGTTCTCGCGAACCgggg                                                                                                                                             <  5:250862/149‑1 (MQ=255)
         ccacaGAAGCCATGTTTGCCGACCTTCCGGAAGACAAGCCAGCTGCTGGCATGCCTGATATGGGTGGCATGGGTGGCATGGGCGTCATTATGTAAGCCAGCCTTACCCCCCTGCACCATCAAAAAGCCCCGGTTCGCGAGAACCGGGGc                                                                                                                                            <  3:174787/149‑1 (MQ=255)
          caccGAAGCCATGGTTGCCGACCTGCCGGAAGACAAGCCAGCTGCTGGCATGCCTGATATGGGTGGCATGGGTGGCATGGGCGGCATGATGTAAGCCAGCCTTACCCCCCTGCACCATCAAAAAGCCCCGGTTCGCGAGAACCGGGGCt                                                                                                                                           >  4:194410/1‑149 (MQ=255)
             cGAAGCCATGGTTGCCGACCTGCCGGAAGACAAGCCAGCTGCTGGCATGCCTGATATGGGTGGCATGGGTGGCATGGGCGGCATGATGTAAGCCAGCCTTACCCCCCTGCACCATCAAAAAGCCCCGGTTCGCGAGAACCGGGGCtttt                                                                                                                                        <  5:205505/149‑1 (MQ=255)
              gAAGCCATGGTTGCCGACCTGCCGGAAGACAAGCCAGCTGCTGGCATGCCTGATATGGGTGGCATGGGTGGCATGGGCGGCATGATGTAAGCCAGCCTTACCCCCCTGCACCATCAAAAAGCCCCGGTTCGCGAGAACCGGGGCttttt                                                                                                                                       >  4:75902/1‑149 (MQ=255)
                 gCCAGGGTTGCCGACCTGCCGGAAGACCAGCCAGCTGCTGGCATGCCTGATATGGGTGGCATGGGTGGCATGGGCGGCATGATGTAAGCCAGCCTTACCCCCCTGCACCATCAAAAAGCCCCGGGTCGCGAGAACCGGGGCTTTTTTAt                                                                                                                                    <  4:259584/149‑1 (MQ=255)
                      ggTTGCCGACCTGCCGGAAGACAAGCCAGCTGCTGGCATGCCTGATATGGGTGGCATGGGTGGCATGGGCGGCATGATGTAAGCCAGCCTTACCCCCCTGCACCATCAAAAAGCCCCGGTTCGCGAGAACCGGGGCTTTTTTATTGCTg                                                                                                                               >  5:67938/1‑149 (MQ=255)
                       gTTGCCGACCTGCCGGAAGACAAGCCAGCTGCTGGCATGCCTGATATGGGTGGCATGGGTGGCATGGGCGGCATGATGTAAGCCAGCCTTACCCCCCTGCACCATCAAAAAGCCCCGGTTCGCGAGAACCGGGGCTTTTTTATTGCTGt                                                                                                                              >  6:303137/1‑149 (MQ=255)
                          gCCGACCTGCCGGAAGACAAGCCAGCTGCTGGCATGCCTGATATGGGTGGCATGGGTGGCATGGGCGGCATGATGTAAGCCAGCCTTACCCCCCTGCACCATCAAAAAGCCCCGGTTCGCGAGAACCGGGGCTTTTTTATTGCTGTCgg                                                                                                                           >  8:146649/1‑149 (MQ=255)
                           ccGACCTGCCGGAAGACAAGCCAGCTGCTGGCATGCCTGATATGGGTGGCATGGGTGGCATGGGCGGCATGATGTAAGCCAGCCTTACCCCCCTGCACCATCAAAAAGCCCCGGTTCGCGAGAACCGGGGCTTTTTTATTGCTGTCggg                                                                                                                          <  6:283092/149‑1 (MQ=255)
                           ccGACCTGCCGGAAGACAAGCCAGCTGCTGGCATGCCTGATATGGGTGGCATGGGTGGCATGGGCGGCATGATGTAAGCCAGCCTTACCCCCCTGCACCATCAAAAAGCCCCGGTTCGCGAGAACCGGGGCTTTTTTATTGCTGTCggg                                                                                                                          <  5:303137/149‑1 (MQ=255)
                                  gCCGGAAGACAAGCCAGCTGCTGGCATGCCTGATATGGGTGGCATGGGTGGCATGGGCGGCATGATGTAAGCCAGCCTTACCCCCCTGCCCCATCAAAAAGCCCCGGTTCGCGAGAACCGGGGCTTTTTTATTGCTGTCGGGGCC‑T‑‑Ctt                                                                                                                >  2:194508/1‑149 (MQ=255)
                                   ccGGAAGACAAGCCAGCTGCTGGCATGCCTGATATGGGTGGCATGGGTGGCATGGGCGGCATGATGTAAGCCAGCCTTACCCCCCTGCACCATCAAAAAGCCCCGGTTCGCGAGAACCGGGGCTTTTTTATTGCTGTCGGGGCC‑T‑‑CTTc                                                                                                               >  6:298353/1‑149 (MQ=255)
                                     ggAAGACAAGCCAGCTGCTGGCATGCCTGATATGGGTGGCATGGGTGGCATGGGCGGCATGATGTAAGCCAGCCTTACCCCCCTGCACCATCAAAAAGCCCCGGTTCGCGAGAACCGGGGCTTTTTTATTGCTGTCGGGGCC‑T‑‑CTTCg                                                                                                              >  8:8540/1‑148 (MQ=255)
                                      gAAGACAAGCCAGCTGCTGGCATGCCTGATATGGGTGGCATGGGTGGCATGGGCGGCATGATGTAAGCCAGCCTTACCCCCCTGCACCATCAAAAAGCCCCGGTTCGCGAGAACCGGGGCTTTTTTATTGCTGTCGGGGCC‑T‑‑CTTCgca                                                                                                            <  8:138317/149‑1 (MQ=255)
                                        aGACAAGCCAGCTGCTGGCATGCCTGATATGGGTGGCATGGGTGGCATGGGCGGCATGATGTAAGCCAGCCTTACCCCCCTGCACCATCAAAAAGCCCCGGTTCGCGAGAACCGGGGCTTTTTTATTg                                                                                                                                  <  2:270232/128‑1 (MQ=255)
                                        aGACAAGCCAGCTGCTGGCATGCCTGATATGGGTGGCATGGGTGGCATGGGCGGCATGATGTAAGCCAGCCTTACCCCCCTGCACCATCAAAAAGCCCCGGTTCGCGAGAACCGGGGCTTTTTTATTg                                                                                                                                  >  1:270232/1‑128 (MQ=255)
                                            aaGCCAGCTGCTGGCATGCCTGATATGGGTGGCATGGGTGGCATGGGCGGCATGATGTAAGCCAGCCTTACCCCCCTGCACCATCAAAAAGCCCCGGTTCGCGAGAACCGGGGCTTTTTTATTGCTGTCGGGGCC‑T‑‑CTTCGCAGCACaa                                                                                                      <  3:194410/149‑1 (MQ=255)
                                                     gctgGCATGCCTGATATGGGTGGCATGGGTGGCATGGGCGGCATGATGTAAGCCAGCCTTACCCCCCTGCACCATCAAAAAGCCCCGGTTCGCGAGAACCGGGGCTTTTTTATTGCTGTCGGGGCC‑T‑‑CTTCGCAGCACAAGGCTGCGCt                                                                                             >  6:76876/1‑149 (MQ=255)
                                                       tgGCATGCCTGATATGGGTGGCATGGGTGGCATGGGCGGCATGATGTAAGCCAGCCTTACCCCCCTGCACCATCAAAAAGCCCCGGTTCGCGAGAACCGGGGCTTTTTTATTGCTGTCGGGGCC‑T‑‑CTTCGCAGCACAAGGCTGCGCtt                                                                                            >  3:277825/1‑148 (MQ=255)
                                                         gCATGCCTGATATGGGTGGCATGGGTGGCATGGGCGGCATGATGTAAGCCAGCCTTACCCCCCTGCACCATCAAAAAGCCCCGGTTCGCGAGAACCGGGGCTTTTTTATTGCTGTCGGGGCC‑T‑‑CTTCGCAGCACAAGGCTGCGCTTGa                                                                                          <  4:277825/148‑1 (MQ=255)
                                                             gCCTGATATGGGTGGCATGGGTGGCATGGGCGGCATGATGTAAGCCAGCCTTACCCCCCTGCACCATCAAAAAGCCCCGGTTCGCGAGAACCGGGGCTTTTTTATTGCTGTCGGGGCC‑T‑‑CTTCGCAGCACAAGGCTGCGCTTGAGCCtt                                                                                     >  7:227208/1‑149 (MQ=255)
                                                              ccTGATATGGGTGGCATGGGTGGCATGGGCGGCATGATGTAAGCCAGCCTTACCCCCCTGCACCATCAAAAAGCCCCGGTTCGCGAGAACCGGGGCTTTTTTATTGCTGTCGGGGCC‑T‑‑CTTCGCAGCACAAGGCTGCGCTTGAGCCTtg                                                                                    <  8:13548/149‑1 (MQ=255)
                                                                    ctgggtggcatgggtggcaTGGGCGGCATGATGTAAGCCAGCCTTACCCCCCTGCACCATCAAAAAGCCCCGGTTCGCGAGAACCGGGGCTTTTTTATTGCTGTCGGGGCC‑T‑‑CTTCGCAGCACAAGGCTGCGCTTGAGCCTTGTGCtg                                                                               <  4:66766/147‑1 (MQ=255)
                                                                      gggtggcatgggtggcatGGGCGGCATGATGTAAGCCAGCCTTACCCCCCTGCACCATCAAAAAGCCCCGGTTCGCGAGAACCGGGGCTTTTTTATTGCTGTCGGGGCC‑T‑‑CTTCGCAGCACAAgg                                                                                                    >  7:270093/1‑125 (MQ=255)
                                                                      gggtggcatgggtggcatGGGCGGCATGATGTAAGCCAGCCTTACCCCCCTGCACCATCAAAAAGCCCCGGTTCGCGAGAACCGGGGCTTTTTTATTGCTGTCGGGGCC‑T‑‑CTTCGCAGCACAAgg                                                                                                    <  8:270093/125‑1 (MQ=255)
                                                                            catgggtggcatgggCGGCATGATGTAAGCCAGCCTTACCCCCCTGCACCATCAAAAAGCCCCGGTTCGCGAGAACCGGGGCTTTTTTATTGCTGTCGGGGCC‑T‑‑CTTCGCAGcac                                                                                                        >  6:154210/1‑115 (MQ=255)
                                                                            catgggtggcatgggCGGCATGATGTAAGCCAGCCTTACCCCCCTGCACCATCAAAAAGCCCCGGTTCGCGAGAACCGGGGCTTTTTTATTGCTGTCGGGGCC‑T‑‑CTTCGCAGCACAAGGCTGCGCTTGAGCCTTGTGCTGTGAAACGAg                                                                      <  8:227208/149‑1 (MQ=255)
                                                                             atgggtggcatgggCGGCATGATGTAAGCCAGCCTTACCCCCCTGCACCATCAAAAAGCCCCGGTTCGCGAGAACCGGGGCTTTTTTATTGCTGTCGGGGCC‑T‑‑CTTCGCAGCACAAGGCTGCGCTTGAGCCTTGTGCTGTGAAACGAGt                                                                     <  7:146649/149‑1 (MQ=255)
                                                                              tgggtggcatgggCGGCATGATGTAAGCCAGCCTTACCCCCCTGCACCATCAAAAAGCCCCGGTTCGCGAGAACCGGGGCTTTTTTATTGCTGTCGGGGCC‑T‑‑CTTCGCAGCACAAGGCTGCGCTTGAGCCTTGTGCTGTGAAACGAGTg                                                                    <  3:75902/149‑1 (MQ=255)
                                                                                     catgggCGGCATTATGTAAGCCAGCCTTACCCCCCTGCACCATCAAAAAGCCCCGGTTCGCGAGAACCGGGGCTTTTTTATTGCTGTCGGGGCC‑T‑‑CTTCGCAGCACAAGGCTGCGCTTGAGCCTTGTGCTGTGAAACGAGTGGATCAAc                                                             <  5:146790/149‑1 (MQ=255)
                                                                                     catgggCGGCATGATGTAAGCCAGCCTTACCCCCCTGCACCATCAAAAAGCCCCGGTTCGCGAGAACCGGGGCTTTTTTATTGCTGTCGGGGCC‑T‑‑CTTCGCAGCACAAGGCTGCGCTTGAGCCTTGTGCTGTGAAACGAGTGGATCAAc                                                             <  3:272464/149‑1 (MQ=255)
                                                                                      atgggCGGCATGATGTAAGCCAGCCTTACCCCCCTGCACCATCAAAAAGCCCCGGTTCGCGAGAACCGGGGCTTTTTTATTGCTGTAGGGGCCATGGGGGCGCAGCACAAGCCTGCGCTTGAGAGTTGTGATGTGAAACTAGTGTATc                                                                >  1:102018/1‑148 (MQ=255)
                                                                                              cATGATGTAAGCCAGCCTTACCCCCCTGCACCATCAAAAAGCCCCGGTTCGCGAGAACCGGGGCTTTTTTATTGCTGTCGGGGCC‑T‑‑CTTCGCAGCACAAGGCTGCGCTTGAGCCTTGTGCTGTGAAGCGAGTGGATCAACCGACCAGCg                                                    >  4:199023/1‑149 (MQ=255)
                                                                                               atgatgTAAGCCAGCCTTACCCCCCTGCACCATCAAAAAGCCCCGGTTCGCGAGAACCGGGGCTTTTTTATTGCTGTCGGGGCC‑T‑‑CTTCGCAGCACAAGGCTGCGCTTGAGCCTTGTGCTGTGAAACGAGTGGATCAACCGACCAGCg                                                    >  8:12264/1‑148 (MQ=255)
                                                                                                   tgTAAGCCAGCCTTACCCCCCTGCACCATCAAAAAGCCCCGGTTCGCGAGAACCGGGGCTTTTTTATTGCTGTCGGGGCC‑T‑‑CTTCGCAGCACAAGGCTGCGCTTGAGCCTTGTGCTGTGAAACGAGTGGATCAACCGACCAGCGGATc                                                <  7:226809/148‑1 (MQ=255)
                                                                                                    gTAAGCCAGCCTTACCCCCCTGCACCATCAAAAAGCCCCGGTTCGCGAGAACCGGGGCTTTTTTATTGCTGTCGGGGCC‑T‑‑CTTCGCAGCACAAGGCTGCGCTTGAGCCTTGTGCTGTGAAACGAGTGGATCAACCGACCAGCGGATCAc                                              >  8:107560/1‑149 (MQ=255)
                                                                                                       agccagccTTACCCCCCTGCACCATCAAAACTCCCCGGTGCGCGAGAACCGGGGCTTTTTTATTGCTGTCGGGGCC‑T‑‑CTTCGCAGCACAAGGCTGCGCTTGAGCCTTGTGCTGTGAAACGAGTgg                                                                   <  2:160545/125‑1 (MQ=255)
                                                                                                       agccagccTTACCCCCCTGCACCATCAAAAAGCCCCGGTTCGCGAGAACCGGGGCTTTTTTATTGCTGTCGGGGCC‑T‑‑CTTCGCAGCACAAGGCTGCGCTTGAGCCTTGTGCTGTGAAACGAGTgg                                                                   >  1:160545/1‑125 (MQ=255)
                                                                                                           agccTTACCCCCCTGCACCATCAAAAAGCCCCGGTTCGCGAGAACCGGGGCTTTTTTATTGCTGTCGGGGCC‑T‑‑CTTCGCAGCACAAGGCTGCGCTTGAGCCTTGTGCTGTGAAACGAGTGGATCAACCGACCAGCgg                                                   >  6:237103/1‑137 (MQ=255)
                                                                                                           agccTTACCCCCCTGCACCATCAAAAAGCCCCGGTTCGCGAGAACCGGGGCTTTTTTATTGCTGTCGGGGCC‑T‑‑CTTCGCAGCACAAGGCTGCGCTTGAGCCTTGTGCTGTGAAACGAGTGGATCAACCGACCAGCgg                                                   <  5:237103/137‑1 (MQ=255)
                                                                                                             ccTTACCCCCCTGCACCATCAAAAAGCCCCGGTTCGCGAGAACCGGGGCTTTTTTATTGCTGTCGGGGCC‑T‑‑CTTCGCAGCACAAGGCTGCGCTTGAGCCTTGTGCTGTGAAACGAGTGGATCAACCGACCAGCGGATCACCGACGTg                                       >  5:297315/1‑147 (MQ=255)
                                                                                                                   cccccTGCACCATCAAAAAGCCCCGGTTCGCGAGAACCGGGGCTTTTTTATTGCTGTCGGGGCC‑T‑‑CTTCGCAGCACAAGGCTGCGCTTGAGCCTTGTGCTGTGAAACGAGTGGATCAACCGACCAGCGGATCACCGACGTGCAGGATCt                               >  6:303658/1‑149 (MQ=255)
                                                                                                                   cccccTGCACCATCAAAAAGCCCCGGTTCGCGAGAACCGGGGCTTTTTTATTGCTGTCGGGGCC‑T‑‑CTTCGCAGCACAAGGCTGCGCTTGAGCCTTGTGCTGTGAAACGAGTGGATCAACCGACCAGCGGATCACCGACGTGCAGGATCt                               >  8:2319/1‑149 (MQ=255)
                                                                                                                     cccTGCACCATCAAAAAGCCCCGGTTCGCGAGAACCGGGGCTTTTTTATTGCTGTCGGGGCC‑T‑‑CTTCGCAGCACAAGGCTGCGCTTGAGCCTTGTGCTGTGAAACGAGTGGATCAACCGCCCAGCGGATCACCGACGTGCAGGATCtt                              >  8:272482/1‑148 (MQ=255)
                                                                                                                     cccTGCACCATCAAAAAGCCCCGGTTCGCGAGAACCGGGGCTTTTTTATTGCTGTCGGGGCC‑T‑‑CTTCGCAGCACAAGGCTGCGCTTGAGCCTTGTGCTGTGAAACGAGTGGATCAACCGACCAGCGGATCACCGACGTGCAGGATCTTc                             <  6:67938/149‑1 (MQ=255)
                                                                                                                           gcaataaaaaaGCCCCGGTTCTCGCGAACCGGGGCTTTTTGA‑TGGTGCAGGGGCC‑T‑‑CTTCGCAGCACAAGGCTGCGCTTGAGCCTTGTGCTGTGAAACGAGTGGATCAACCGACCAGCGGATCACCGACGTGCAGGATCTTCATGCCAt                      <  3:80699/143‑1 (MQ=25)
                                                                                                                             cATCAAATAGCCCAGGTTCGCGAGAAGCGGTGCTTTTTTATTGCTTTCGTGGCC‑T‑‑CTTCGCAGCACAAGGCTGCGCTTGAGCCTTGTGCTGTGAAACGAGTTGATCAACCGACCAGCGGATCACCGACGt                                        <  2:254514/130‑1 (MQ=255)
                                                                                                                             cATCAAAAAGCCCCTGTTCGCGAGAACCGGGGCTTTTTTATTGATGTCGGGGCC‑T‑‑CTTCGCAGCACAACGATGCGCTTGAGCCTTGTGCAGTGAAACGAGTGGATCAACAGACCAGCCGATCACCGACGt                                        >  1:254514/1‑130 (MQ=255)
                                                                                                                             cATCAAAAAGCCCCGGTTCGCGAGAACCGGGGCTTTTTTATTGCTGTCGGGGCC‑T‑‑CTTCGCAGCACAAGGCTGCGCTTGAGCCTTGTGCTGTGAAACGAGTGGATCAACCGACCAGCGGATCACCGACGTGCAGGATCTTCATGCCAtt                     <  6:297315/149‑1 (MQ=255)
                                                                                                                              aTCAAAAAGCCCCGGTTCGCGAGAACCGGGGCTTTTTTATTGCTGTCGGGGCC‑T‑‑CTTCGCAGCACAAGGCTGCGCTTGAGCCTTGTGCTGTGAAACGAGTGGATCAACCGACCAGCGGATCACCGACGTGCAg                                    >  7:11462/1‑133 (MQ=255)
                                                                                                                              aTCAAAAAGCCCCGGTTCGCGAGAACCGGGGCTTTTTTATTGCTGTCGGGGCC‑T‑‑CTTCGCAGCACAAGGCTGCGCTTGAGCCTTGTGCTGTGAAACGAGTGGATCAACCGACCAGCGGATCACCGACGTGCAg                                    <  8:11462/133‑1 (MQ=255)
                                                                                                                                  aaaaGCCCCGGTTCTCGCGAACCGGGGCTTTTTTATTGCTGTCGGGGCC‑T‑‑CTTCGCAGCACAAGGCTGCGCTTGAGCCTTGTGCTGTGAAACGAGTGGATCAACCGACCAGCGGATCACCGACGTGCAGGATCTTCATGCCATTGGTGc                >  1:8173/1‑149 (MQ=255)
                                                                                                                                      gCCCCGGTTCTCGCGAACCGGGGCTTTTTTATTGCTGTCGGGGCC‑T‑‑CTTCGCAGCACAAGGCTGCGCTTGAGCCTTGTGCTGTGAAACGAGTGGATCAACCGACCAGCGGATCACCGACGTGCAGGATCTTc                             >  6:151312/1‑132 (MQ=255)
                                                                                                                                      gCCCCGGTTCTCGCGAACCGGGGCTTTTTTATTGCTGTCGGGGCC‑T‑‑CTTCGCAGCACAAGGCTGCGCTTGAGCCTTGTGCTGTGAAACGAGTGGATCAACCGACCAGCGGATCACCGACGTGCAGGATCTTc                             <  5:151312/132‑1 (MQ=255)
                                                                                                                                       ccccGGTTCGCGAGAACCGGGGCTTTTTTATTGCTGTCGGGGCC‑T‑‑CTTCGCAGCACAAGGCTGCGCTTGAGCCTTGTGCTGTGAAACGAGTGGATCAACCGACCAGCGGATCACCGACGTGCAGGATCTTCATGCCATTGGTGCCAcc            >  2:20915/1‑148 (MQ=255)
                                                                                                                                        cccGGTTCGCGAGAACCGGGGCTTTTTTATTGCTGTCGGGGCC‑T‑‑CTTCGCAGCACAAGGCTGCGCTTGAGCCTTGTGCTGTGAAACGAGTGGATCAACCGACCAGCGGATCACCGACGTGCAGGATCTTCATGCCATTGGTGCCAc             >  4:34929/1‑146 (MQ=255)
                                                                                                                                        cccGGTTCGCGAGAACCGGGGCTTTTTTATTGCTGTCGGGGCC‑T‑‑CTTCGCAGCACAAGGCTGCGCTTGAGCCTTGTGCTGTGAAACGAGTGGATCAACCGACCAGCGGATCACCGACGTGCAGGATCTTCATGCCATTGGTGCCAc             <  3:34929/146‑1 (MQ=255)
                                                                                                                                        cccGGTTCGCGAGAACCGGGGCTTTTTTATTGCTGTCGGGGCC‑T‑‑CTTCGCAGCACAAGGCTGCGCTTGAGCCTTGTGCTGTGAAACGAGTGGATCAACCGACCAGCGGATCACCGACGTGCAGGATCTTCATGCCATTGGTGCCACCGa          <  5:76876/149‑1 (MQ=255)
                                                                                                                                          cGGTTCGCGAGAACCGGGGCTTTTTTATTGCTGTCGGGGCC‑T‑‑CTTCGCAGCACAAGGCTGCGCTTGAGCCTTGTGCTGTGAAACGAGTGGATCAACCGACCAGCGGATCACCGACGTGCAGGATCTTCATGCCATTGGTGCCACCGATg        <  5:303658/149‑1 (MQ=255)
                                                                                                                                           ggTTCGCGAGAACCGGGGCTTTTTTATTGCTGTCGGGGCC‑T‑‑CTTCGCAGCACAAGGCTGCGCTTGAGCCTTGTGCTGTGAAACGAGTGGATCAACCGa                                                          >  1:273979/1‑98 (MQ=255)
                                                                                                                                           ggTTCGCGAGAACCGGGGCTTTTTTATTGCTGTCGGGGCC‑T‑‑CTTCGCAGCACAAGGCTGCGCTTGAGCCTTGTGCTGTGAAACGAGTGGATCAACCGa                                                          <  2:273979/98‑1 (MQ=255)
                                                                                                                                            gTTCGCGAGAACCGGGGCTTTTTTATTGCTGTCGGGGCC‑T‑‑CTTCGCAGCACAAGGCTGCGCTTGAGCCTTGTGCTGTGAAACGAGTGGATCAACCGACCAGCGGATCACCGACGTGCAGGATCTTCATGCCATTGGTGCCACCGATGgt      <  1:194508/149‑1 (MQ=255)
                                                                                                                                             ttCTCGCGAACCGGGGCTTTTTTATTGCTGTCGGGGCC‑T‑‑CTTCGCAGCACAAGGCTGCGCTTGAGCCTTGTGCTGTGAAACGAGTGGATCAACCGACCAGCGGATCACCGACGTGCAGGATCTTCATGCCATTGGTGCCACCGATGgtg     <  2:8173/149‑1 (MQ=255)
                                                                                                                                                gcgAGAACCGGGGCTTTTTTATTGCTGTCGGGGCC‑T‑‑CTTCGCAGCACAAGGCTGCGCTTGAGCCTTGTGCTGTGAAACGAGTGGATCAACCGACCAGCGGATCACCGACGTGCAGGATCTTCATGCCATTGGTGCCACCGATGGTGTgg  <  7:8540/149‑1 (MQ=255)
                                                                                                                                                |                                                                                                                                                       
TGATGATCACCACCGAAGCCATGGTTGCCGACCTGCCGGAAGACAAGCCAGCTGCTGGCATGCCTGATATGGGTGGCATGGGTGGCATGGGCGGCATGATGTAAGCCAGCCTTACCCCCCTGCACCATCAAAAAGCCCCGGTTCGCGAGAACCGGGGCTTTTTTATTGCTGTCGGGGCC‑T‑‑CTTCGCAGCACAAGGCTGCGCTTGAGCCTTGTGCTGTGAAACGAGTGGATCAACCGACCAGCGGATCACCGACGTGCAGGATCTTCATGCCATTGGTGCCACCGATGGTGTGG  >  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: ‑