Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 1,551,285 A→C 8.7% 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.7% 199.2 / 9.3 69intergenic (+44/+82)groL/pykAchaperonin GroEL/pyruvate kinase II
Reads supporting (aligned to +/- strand):  ref base A (31/32);  new base C (2/4);  total (33/36)
Fisher's exact test for biased strand distribution p-value = 6.75e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.33e-01

ATGATCACCACCGAAGCCATGGTTGCCGACCTGCCGGAAGACAAGCCAGCTGCTGGCATGCCTGATATGGGTGGCATGGGTGGCATGGGCGGCATGATGTAAGCCAGCCTTACCCCCCTGCACCATCAAAAAGCCCCGGTTCGCGAGAACCGGGGCTTTTTTATTGCTGTCGGGGCCT‑‑CTTCGCAGCACAAGGCTGCGCTTGAGCCTTGTGCTGTGAAACGAGTGGATCAACCGACCAGCGGATCACCGACGTGCAGGATCTTCATGCCATTGGTGCCACCGATGGTGTGG  >  NC_002947/1551140‑1551430
                                                                                                                                                 |                                                                                                                                                    
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)
                                gCCGGAAGACAAGCCAGCTGCTGGCATGCCTGATATGGGTGGCATGGGTGGCATGGGCGGCATGATGTAAGCCAGCCTTACCCCCCTGCCCCATCAAAAAGCCCCGGTTCGCGAGAACCGGGGCTTTTTTATTGCTGTCGGGGCCT‑‑Ctt                                                                                                                >  2:194508/1‑149 (MQ=255)
                                 ccGGAAGACAAGCCAGCTGCTGGCATGCCTGATATGGGTGGCATGGGTGGCATGGGCGGCATGATGTAAGCCAGCCTTACCCCCCTGCACCATCAAAAAGCCCCGGTTCGCGAGAACCGGGGCTTTTTTATTGCTGTCGGGGCCT‑‑CTTc                                                                                                               >  6:298353/1‑149 (MQ=255)
                                   ggAAGACAAGCCAGCTGCTGGCATGCCTGATATGGGTGGCATGGGTGGCATGGGCGGCATGATGTAAGCCAGCCTTACCCCCCTGCACCATCAAAAAGCCCCGGTTCGCGAGAACCGGGGCTTTTTTATTGCTGTCGGGGCCT‑‑CTTCg                                                                                                              >  8:8540/1‑148 (MQ=255)
                                    gAAGACAAGCCAGCTGCTGGCATGCCTGATATGGGTGGCATGGGTGGCATGGGCGGCATGATGTAAGCCAGCCTTACCCCCCTGCACCATCAAAAAGCCCCGGTTCGCGAGAACCGGGGCTTTTTTATTGCTGTCGGGGCCT‑‑CTTCgca                                                                                                            <  8:138317/149‑1 (MQ=255)
                                      aGACAAGCCAGCTGCTGGCATGCCTGATATGGGTGGCATGGGTGGCATGGGCGGCATGATGTAAGCCAGCCTTACCCCCCTGCACCATCAAAAAGCCCCGGTTCGCGAGAACCGGGGCTTTTTTATTg                                                                                                                                  >  1:270232/1‑128 (MQ=255)
                                      aGACAAGCCAGCTGCTGGCATGCCTGATATGGGTGGCATGGGTGGCATGGGCGGCATGATGTAAGCCAGCCTTACCCCCCTGCACCATCAAAAAGCCCCGGTTCGCGAGAACCGGGGCTTTTTTATTg                                                                                                                                  <  2:270232/128‑1 (MQ=255)
                                          aaGCCAGCTGCTGGCATGCCTGATATGGGTGGCATGGGTGGCATGGGCGGCATGATGTAAGCCAGCCTTACCCCCCTGCACCATCAAAAAGCCCCGGTTCGCGAGAACCGGGGCTTTTTTATTGCTGTCGGGGCCT‑‑CTTCGCAGCACaa                                                                                                      <  3:194410/149‑1 (MQ=255)
                                                   gctgGCATGCCTGATATGGGTGGCATGGGTGGCATGGGCGGCATGATGTAAGCCAGCCTTACCCCCCTGCACCATCAAAAAGCCCCGGTTCGCGAGAACCGGGGCTTTTTTATTGCTGTCGGGGCCT‑‑CTTCGCAGCACAAGGCTGCGCt                                                                                             >  6:76876/1‑149 (MQ=255)
                                                     tgGCATGCCTGATATGGGTGGCATGGGTGGCATGGGCGGCATGATGTAAGCCAGCCTTACCCCCCTGCACCATCAAAAAGCCCCGGTTCGCGAGAACCGGGGCTTTTTTATTGCTGTCGGGGCCT‑‑CTTCGCAGCACAAGGCTGCGCtt                                                                                            >  3:277825/1‑148 (MQ=255)
                                                       gCATGCCTGATATGGGTGGCATGGGTGGCATGGGCGGCATGATGTAAGCCAGCCTTACCCCCCTGCACCATCAAAAAGCCCCGGTTCGCGAGAACCGGGGCTTTTTTATTGCTGTCGGGGCCT‑‑CTTCGCAGCACAAGGCTGCGCTTGa                                                                                          <  4:277825/148‑1 (MQ=255)
                                                           gCCTGATATGGGTGGCATGGGTGGCATGGGCGGCATGATGTAAGCCAGCCTTACCCCCCTGCACCATCAAAAAGCCCCGGTTCGCGAGAACCGGGGCTTTTTTATTGCTGTCGGGGCCT‑‑CTTCGCAGCACAAGGCTGCGCTTGAGCCtt                                                                                     >  7:227208/1‑149 (MQ=255)
                                                            ccTGATATGGGTGGCATGGGTGGCATGGGCGGCATGATGTAAGCCAGCCTTACCCCCCTGCACCATCAAAAAGCCCCGGTTCGCGAGAACCGGGGCTTTTTTATTGCTGTCGGGGCCT‑‑CTTCGCAGCACAAGGCTGCGCTTGAGCCTtg                                                                                    <  8:13548/149‑1 (MQ=255)
                                                                  ctgggtggcatgggtggcaTGGGCGGCATGATGTAAGCCAGCCTTACCCCCCTGCACCATCAAAAAGCCCCGGTTCGCGAGAACCGGGGCTTTTTTATTGCTGTCGGGGCCT‑‑CTTCGCAGCACAAGGCTGCGCTTGAGCCTTGTGCtg                                                                               <  4:66766/147‑1 (MQ=255)
                                                                    gggtggcatgggtggcatGGGCGGCATGATGTAAGCCAGCCTTACCCCCCTGCACCATCAAAAAGCCCCGGTTCGCGAGAACCGGGGCTTTTTTATTGCTGTCGGGGCCT‑‑CTTCGCAGCACAAgg                                                                                                    >  7:270093/1‑125 (MQ=255)
                                                                    gggtggcatgggtggcatGGGCGGCATGATGTAAGCCAGCCTTACCCCCCTGCACCATCAAAAAGCCCCGGTTCGCGAGAACCGGGGCTTTTTTATTGCTGTCGGGGCCT‑‑CTTCGCAGCACAAgg                                                                                                    <  8:270093/125‑1 (MQ=255)
                                                                          catgggtggcatgggCGGCATGATGTAAGCCAGCCTTACCCCCCTGCACCATCAAAAAGCCCCGGTTCGCGAGAACCGGGGCTTTTTTATTGCTGTCGGGGCCT‑‑CTTCGCAGcac                                                                                                        >  6:154210/1‑115 (MQ=255)
                                                                          catgggtggcatgggCGGCATGATGTAAGCCAGCCTTACCCCCCTGCACCATCAAAAAGCCCCGGTTCGCGAGAACCGGGGCTTTTTTATTGCTGTCGGGGCCT‑‑CTTCGCAGCACAAGGCTGCGCTTGAGCCTTGTGCTGTGAAACGAg                                                                      <  8:227208/149‑1 (MQ=255)
                                                                           atgggtggcatgggCGGCATGATGTAAGCCAGCCTTACCCCCCTGCACCATCAAAAAGCCCCGGTTCGCGAGAACCGGGGCTTTTTTATTGCTGTCGGGGCCT‑‑CTTCGCAGCACAAGGCTGCGCTTGAGCCTTGTGCTGTGAAACGAGt                                                                     <  7:146649/149‑1 (MQ=255)
                                                                            tgggtggcatgggCGGCATGATGTAAGCCAGCCTTACCCCCCTGCACCATCAAAAAGCCCCGGTTCGCGAGAACCGGGGCTTTTTTATTGCTGTCGGGGCCT‑‑CTTCGCAGCACAAGGCTGCGCTTGAGCCTTGTGCTGTGAAACGAGTg                                                                    <  3:75902/149‑1 (MQ=255)
                                                                                   catgggCGGCATTATGTAAGCCAGCCTTACCCCCCTGCACCATCAAAAAGCCCCGGTTCGCGAGAACCGGGGCTTTTTTATTGCTGTCGGGGCCT‑‑CTTCGCAGCACAAGGCTGCGCTTGAGCCTTGTGCTGTGAAACGAGTGGATCAAc                                                             <  5:146790/149‑1 (MQ=255)
                                                                                   catgggCGGCATGATGTAAGCCAGCCTTACCCCCCTGCACCATCAAAAAGCCCCGGTTCGCGAGAACCGGGGCTTTTTTATTGCTGTCGGGGCCT‑‑CTTCGCAGCACAAGGCTGCGCTTGAGCCTTGTGCTGTGAAACGAGTGGATCAAc                                                             <  3:272464/149‑1 (MQ=255)
                                                                                    atgggCGGCATGATGTAAGCCAGCCTTACCCCCCTGCACCATCAAAAAGCCCCGGTTCGCGAGAACCGGGGCTTTTTTATTGCTGTAGGGGCCATGGGGGCGCAGCACAAGCCTGCGCTTGAGAGTTGTGATGTGAAACTAGTGTATc                                                                >  1:102018/1‑148 (MQ=255)
                                                                                            cATGATGTAAGCCAGCCTTACCCCCCTGCACCATCAAAAAGCCCCGGTTCGCGAGAACCGGGGCTTTTTTATTGCTGTCGGGGCCT‑‑CTTCGCAGCACAAGGCTGCGCTTGAGCCTTGTGCTGTGAAGCGAGTGGATCAACCGACCAGCg                                                    >  4:199023/1‑149 (MQ=255)
                                                                                             atgatgTAAGCCAGCCTTACCCCCCTGCACCATCAAAAAGCCCCGGTTCGCGAGAACCGGGGCTTTTTTATTGCTGTCGGGGCCT‑‑CTTCGCAGCACAAGGCTGCGCTTGAGCCTTGTGCTGTGAAACGAGTGGATCAACCGACCAGCg                                                    >  8:12264/1‑148 (MQ=255)
                                                                                                 tgTAAGCCAGCCTTACCCCCCTGCACCATCAAAAAGCCCCGGTTCGCGAGAACCGGGGCTTTTTTATTGCTGTCGGGGCCT‑‑CTTCGCAGCACAAGGCTGCGCTTGAGCCTTGTGCTGTGAAACGAGTGGATCAACCGACCAGCGGATc                                                <  7:226809/148‑1 (MQ=255)
                                                                                                  gTAAGCCAGCCTTACCCCCCTGCACCATCAAAAAGCCCCGGTTCGCGAGAACCGGGGCTTTTTTATTGCTGTCGGGGCCT‑‑CTTCGCAGCACAAGGCTGCGCTTGAGCCTTGTGCTGTGAAACGAGTGGATCAACCGACCAGCGGATCAc                                              >  8:107560/1‑149 (MQ=255)
                                                                                                     agccagccTTACCCCCCTGCACCATCAAAACTCCCCGGTGCGCGAGAACCGGGGCTTTTTTATTGCTGTCGGGGCCT‑‑CTTCGCAGCACAAGGCTGCGCTTGAGCCTTGTGCTGTGAAACGAGTgg                                                                   <  2:160545/125‑1 (MQ=255)
                                                                                                     agccagccTTACCCCCCTGCACCATCAAAAAGCCCCGGTTCGCGAGAACCGGGGCTTTTTTATTGCTGTCGGGGCCT‑‑CTTCGCAGCACAAGGCTGCGCTTGAGCCTTGTGCTGTGAAACGAGTgg                                                                   >  1:160545/1‑125 (MQ=255)
                                                                                                         agccTTACCCCCCTGCACCATCAAAAAGCCCCGGTTCGCGAGAACCGGGGCTTTTTTATTGCTGTCGGGGCCT‑‑CTTCGCAGCACAAGGCTGCGCTTGAGCCTTGTGCTGTGAAACGAGTGGATCAACCGACCAGCgg                                                   >  6:237103/1‑137 (MQ=255)
                                                                                                         agccTTACCCCCCTGCACCATCAAAAAGCCCCGGTTCGCGAGAACCGGGGCTTTTTTATTGCTGTCGGGGCCT‑‑CTTCGCAGCACAAGGCTGCGCTTGAGCCTTGTGCTGTGAAACGAGTGGATCAACCGACCAGCgg                                                   <  5:237103/137‑1 (MQ=255)
                                                                                                           ccTTACCCCCCTGCACCATCAAAAAGCCCCGGTTCGCGAGAACCGGGGCTTTTTTATTGCTGTCGGGGCCT‑‑CTTCGCAGCACAAGGCTGCGCTTGAGCCTTGTGCTGTGAAACGAGTGGATCAACCGACCAGCGGATCACCGACGTg                                       >  5:297315/1‑147 (MQ=255)
                                                                                                                 cccccTGCACCATCAAAAAGCCCCGGTTCGCGAGAACCGGGGCTTTTTTATTGCTGTCGGGGCCT‑‑CTTCGCAGCACAAGGCTGCGCTTGAGCCTTGTGCTGTGAAACGAGTGGATCAACCGACCAGCGGATCACCGACGTGCAGGATCt                               >  8:2319/1‑149 (MQ=255)
                                                                                                                 cccccTGCACCATCAAAAAGCCCCGGTTCGCGAGAACCGGGGCTTTTTTATTGCTGTCGGGGCCT‑‑CTTCGCAGCACAAGGCTGCGCTTGAGCCTTGTGCTGTGAAACGAGTGGATCAACCGACCAGCGGATCACCGACGTGCAGGATCt                               >  6:303658/1‑149 (MQ=255)
                                                                                                                   cccTGCACCATCAAAAAGCCCCGGTTCGCGAGAACCGGGGCTTTTTTATTGCTGTCGGGGCCT‑‑CTTCGCAGCACAAGGCTGCGCTTGAGCCTTGTGCTGTGAAACGAGTGGATCAACCGCCCAGCGGATCACCGACGTGCAGGATCtt                              >  8:272482/1‑148 (MQ=255)
                                                                                                                   cccTGCACCATCAAAAAGCCCCGGTTCGCGAGAACCGGGGCTTTTTTATTGCTGTCGGGGCCT‑‑CTTCGCAGCACAAGGCTGCGCTTGAGCCTTGTGCTGTGAAACGAGTGGATCAACCGACCAGCGGATCACCGACGTGCAGGATCTTc                             <  6:67938/149‑1 (MQ=255)
                                                                                                                         gcaataaaaaaGCCCCGGTTCTCGCGAACCGGGGCTTTTTGATGGTGCAGGGGCCT‑‑CTTCGCAGCACAAGGCTGCGCTTGAGCCTTGTGCTGTGAAACGAGTGGATCAACCGACCAGCGGATCACCGACGTGCAGGATCTTCATGCCAt                      <  3:80699/143‑1 (MQ=25)
                                                                                                                           cATCAAATAGCCCAGGTTCGCGAGAAGCGGTGCTTTTTTATTGCTTTCGTGGCCT‑‑CTTCGCAGCACAAGGCTGCGCTTGAGCCTTGTGCTGTGAAACGAGTTGATCAACCGACCAGCGGATCACCGACGt                                        <  2:254514/130‑1 (MQ=255)
                                                                                                                           cATCAAAAAGCCCCTGTTCGCGAGAACCGGGGCTTTTTTATTGATGTCGGGGCCT‑‑CTTCGCAGCACAACGATGCGCTTGAGCCTTGTGCAGTGAAACGAGTGGATCAACAGACCAGCCGATCACCGACGt                                        >  1:254514/1‑130 (MQ=255)
                                                                                                                           cATCAAAAAGCCCCGGTTCGCGAGAACCGGGGCTTTTTTATTGCTGTCGGGGCCT‑‑CTTCGCAGCACAAGGCTGCGCTTGAGCCTTGTGCTGTGAAACGAGTGGATCAACCGACCAGCGGATCACCGACGTGCAGGATCTTCATGCCAtt                     <  6:297315/149‑1 (MQ=255)
                                                                                                                            aTCAAAAAGCCCCGGTTCGCGAGAACCGGGGCTTTTTTATTGCTGTCGGGGCCT‑‑CTTCGCAGCACAAGGCTGCGCTTGAGCCTTGTGCTGTGAAACGAGTGGATCAACCGACCAGCGGATCACCGACGTGCAg                                    >  7:11462/1‑133 (MQ=255)
                                                                                                                            aTCAAAAAGCCCCGGTTCGCGAGAACCGGGGCTTTTTTATTGCTGTCGGGGCCT‑‑CTTCGCAGCACAAGGCTGCGCTTGAGCCTTGTGCTGTGAAACGAGTGGATCAACCGACCAGCGGATCACCGACGTGCAg                                    <  8:11462/133‑1 (MQ=255)
                                                                                                                                aaaaGCCCCGGTTCTCGCGAACCGGGGCTTTTTTATTGCTGTCGGGGCCT‑‑CTTCGCAGCACAAGGCTGCGCTTGAGCCTTGTGCTGTGAAACGAGTGGATCAACCGACCAGCGGATCACCGACGTGCAGGATCTTCATGCCATTGGTGc                >  1:8173/1‑149 (MQ=255)
                                                                                                                                    gCCCCGGTTCTCGCGAACCGGGGCTTTTTTATTGCTGTCGGGGCCT‑‑CTTCGCAGCACAAGGCTGCGCTTGAGCCTTGTGCTGTGAAACGAGTGGATCAACCGACCAGCGGATCACCGACGTGCAGGATCTTc                             <  5:151312/132‑1 (MQ=255)
                                                                                                                                    gCCCCGGTTCTCGCGAACCGGGGCTTTTTTATTGCTGTCGGGGCCT‑‑CTTCGCAGCACAAGGCTGCGCTTGAGCCTTGTGCTGTGAAACGAGTGGATCAACCGACCAGCGGATCACCGACGTGCAGGATCTTc                             >  6:151312/1‑132 (MQ=255)
                                                                                                                                     ccccGGTTCGCGAGAACCGGGGCTTTTTTATTGCTGTCGGGGCCT‑‑CTTCGCAGCACAAGGCTGCGCTTGAGCCTTGTGCTGTGAAACGAGTGGATCAACCGACCAGCGGATCACCGACGTGCAGGATCTTCATGCCATTGGTGCCAcc            >  2:20915/1‑148 (MQ=255)
                                                                                                                                      cccGGTTCGCGAGAACCGGGGCTTTTTTATTGCTGTCGGGGCCT‑‑CTTCGCAGCACAAGGCTGCGCTTGAGCCTTGTGCTGTGAAACGAGTGGATCAACCGACCAGCGGATCACCGACGTGCAGGATCTTCATGCCATTGGTGCCAc             >  4:34929/1‑146 (MQ=255)
                                                                                                                                      cccGGTTCGCGAGAACCGGGGCTTTTTTATTGCTGTCGGGGCCT‑‑CTTCGCAGCACAAGGCTGCGCTTGAGCCTTGTGCTGTGAAACGAGTGGATCAACCGACCAGCGGATCACCGACGTGCAGGATCTTCATGCCATTGGTGCCAc             <  3:34929/146‑1 (MQ=255)
                                                                                                                                      cccGGTTCGCGAGAACCGGGGCTTTTTTATTGCTGTCGGGGCCT‑‑CTTCGCAGCACAAGGCTGCGCTTGAGCCTTGTGCTGTGAAACGAGTGGATCAACCGACCAGCGGATCACCGACGTGCAGGATCTTCATGCCATTGGTGCCACCGa          <  5:76876/149‑1 (MQ=255)
                                                                                                                                        cGGTTCGCGAGAACCGGGGCTTTTTTATTGCTGTCGGGGCCT‑‑CTTCGCAGCACAAGGCTGCGCTTGAGCCTTGTGCTGTGAAACGAGTGGATCAACCGACCAGCGGATCACCGACGTGCAGGATCTTCATGCCATTGGTGCCACCGATg        <  5:303658/149‑1 (MQ=255)
                                                                                                                                         ggTTCGCGAGAACCGGGGCTTTTTTATTGCTGTCGGGGCCT‑‑CTTCGCAGCACAAGGCTGCGCTTGAGCCTTGTGCTGTGAAACGAGTGGATCAACCGa                                                          <  2:273979/98‑1 (MQ=255)
                                                                                                                                         ggTTCGCGAGAACCGGGGCTTTTTTATTGCTGTCGGGGCCT‑‑CTTCGCAGCACAAGGCTGCGCTTGAGCCTTGTGCTGTGAAACGAGTGGATCAACCGa                                                          >  1:273979/1‑98 (MQ=255)
                                                                                                                                          gTTCGCGAGAACCGGGGCTTTTTTATTGCTGTCGGGGCCT‑‑CTTCGCAGCACAAGGCTGCGCTTGAGCCTTGTGCTGTGAAACGAGTGGATCAACCGACCAGCGGATCACCGACGTGCAGGATCTTCATGCCATTGGTGCCACCGATGgt      <  1:194508/149‑1 (MQ=255)
                                                                                                                                           ttCTCGCGAACCGGGGCTTTTTTATTGCTGTCGGGGCCT‑‑CTTCGCAGCACAAGGCTGCGCTTGAGCCTTGTGCTGTGAAACGAGTGGATCAACCGACCAGCGGATCACCGACGTGCAGGATCTTCATGCCATTGGTGCCACCGATGgtg     <  2:8173/149‑1 (MQ=255)
                                                                                                                                              gcgAGAACCGGGGCTTTTTTATTGCTGTCGGGGCCT‑‑CTTCGCAGCACAAGGCTGCGCTTGAGCCTTGTGCTGTGAAACGAGTGGATCAACCGACCAGCGGATCACCGACGTGCAGGATCTTCATGCCATTGGTGCCACCGATGGTGTgg  <  7:8540/149‑1 (MQ=255)
                                                                                                                                                 |                                                                                                                                                    
ATGATCACCACCGAAGCCATGGTTGCCGACCTGCCGGAAGACAAGCCAGCTGCTGGCATGCCTGATATGGGTGGCATGGGTGGCATGGGCGGCATGATGTAAGCCAGCCTTACCCCCCTGCACCATCAAAAAGCCCCGGTTCGCGAGAACCGGGGCTTTTTTATTGCTGTCGGGGCCT‑‑CTTCGCAGCACAAGGCTGCGCTTGAGCCTTGTGCTGTGAAACGAGTGGATCAACCGACCAGCGGATCACCGACGTGCAGGATCTTCATGCCATTGGTGCCACCGATGGTGTGG  >  NC_002947/1551140‑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: