Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 5,558,375 T→G 5.8% S410R (AGC→CGC)  PP_4888 ← methyl‑accepting chemotaxis transducer

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0029475,558,3750TG 5.8% 682.1 / 20.8 200S410R (AGC→CGC) PP_4888methyl‑accepting chemotaxis transducer
Reads supporting (aligned to +/- strand):  ref base T (100/88);  new base G (7/5);  total (107/93)
Fisher's exact test for biased strand distribution p-value = 7.74e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

Only 100 of 202 total aligned reads displayed.

GTTTCAACCTGCCCGGCAAGGTTCTCAATCAGGCTCGCCAGTTGCGTCGAGGAACGGCTGTTCTCATCCGTCGCCTGGCGAATGTCCACCACCTGCTGCAACGCCGCCTGCATGGCATGGCTTTCAGCCTGAGCCTCGTCAGCCATGCTCGACAGGTCGCGCAGGCTGGCGGCCACTTCGTCACGCTGCAACGCCGCCGCGGCATCGGCACCGGCATTGCGCTGGGCCATGGCGCCGATCTCGACGCCGGTACGCTGGGCCACCTCGCCGGCCTCGCGGACGATGGGCTGCAGCTTGT  >  NC_002947/5558227‑5558524
                                                                                                                                                    |                                                                                                                                                      
gTTTCAACCTGCCCGGCAAGGTTCTCAATCAGGCTCGCCAGTTGCGTCGAGGAACGGCTGTTCTCATCCGTCGCCTGGCGAATGTCCACCACCTGCTGCAACGCCGCCTGCATGGCATGGCTTTCAGCCTGAGCCTCGTCAGCCATGCTc                                                                                                                                                       <  1:496710/150‑1 (MQ=255)
gTTTCAACCTGCCCGGCAAGGTTCTCAATCAGGCTCGCCAGTTGCGTCGAGGAACGGCTGTTCTCATCCGTCGCCTGGCGAATGTCCACCACCTGCTGCAACGCCGCCTGCATGGCATGGCTTTCAGCCTGAGCCTCGTCAGCCATGCTc                                                                                                                                                       <  1:3218983/150‑1 (MQ=255)
gTTTCAACCTGCCCGGCAAGGTTCTCAATCAGGCTCGCCAGTTGCGTCGAGGAACGGCTGTTCTCATCCGTCGCCTGGCGAATGTCCACCACCTGCTGCAACGCCGCCTGCATGGCATGGCTTTCAGCCTGAGCCTCGTCAGCCATGCTc                                                                                                                                                       <  2:4906731/150‑1 (MQ=255)
    cAACCTGCCCGGCAAGGTTCTCAATCAGGCTCGCCAGTTGCGTCGAGGAACGGCTGTTCTCATCCGTCGCCTGGCGAATGTCCACCACCTGCTGCAACGCCGCCTGCATGGCATGGCTTTCAGCCTGAGCCTCGTCAGCCATGCTCGACa                                                                                                                                                   <  1:4576006/150‑1 (MQ=255)
          gCCCGGCAAGGTTCTCAATCAGGCTCGCCAGTTGCGTCGAGGAACGGCTGTTCTCATCCGTCGCCTGGCGAATGTCCACCACCTGCTGCAACGCCGCCTGCATGGCATGGCTTTCAGCCTGAGCCTCGTCAGCCATGCTCGACAGGTcgc                                                                                                                                             >  1:435189/1‑150 (MQ=255)
          gCCCGGCAAGGTTCTCAATCAGGCTCGCCAGTTGCGTCGAGGAACGGCTGTTCTCATCCGTCGCCTGGCGAATGTCCACCACCTGCTGCAACGCCGCCTGCATGGCATGGCTTTCAGCCTGAGCCTCGTCAGCCATGCTCGACAGGTcgc                                                                                                                                             <  1:2062971/150‑1 (MQ=255)
          gCCCGGCAAGGTTCTCAATCAGGCTCGCCAGTTGCGTCGAGGAACGGCTGTTCTCATCCGTCGCCTGGCGAATGTCCACCACCTGCTGCAACGCCGCCTGCATGGCATGGCTTTCAGCCTGAGCCTCGTCAGCCATGCTCGACAGGTcgc                                                                                                                                             <  1:2062643/150‑1 (MQ=255)
            ccGGCAAGGTTCTCAATCAGGCTCGCCAGTTGCGTCGAGGAACGGCTGTTCTCATCCGTCGCCTGGCGAATGTCCACCACCTGCTGCAACGCCGCCTGCATGGCATGGCTTTCAGCCTGAGCCTCGTCAGCCATGCTCGACAGGTCgcgc                                                                                                                                           <  1:94175/150‑1 (MQ=255)
            ccGGCAAGGTTCTCAATCAGGCTCGCCAGTTGCGTCGAGGAACGGCTGTTCTCATCCGTCGCCTGGCGAATGTCCACCACCTGCTGCAACGCCGCCTGCATGGCATGGCTTTCAGCCTGAGCCTCGTCAGCCATGCTCGACAGGTCgcgc                                                                                                                                           <  1:3148683/150‑1 (MQ=255)
             cGGCAAGGTTCTCAATCAGGCTCGCCAGTTGCGTCGAGGAACGGCTGTTCTCATCCGTCGCCTGGCGAATGTCCACCACCTGCTGCAACGCCGCCTGCATGGCATGGCTTTCAGCCTGAGCCTCGTCAGCCATGCTCGACAGGTCGCGCa                                                                                                                                          <  2:4517157/150‑1 (MQ=255)
              ggCAAGGTTCTCAATCAGGCTCGCCAGTTGCGTCGAGGAACGGCTGTTCTCATCCGTCGCCTGGCGAATGTCCACCACCTGCTGCAACGCCGCCTGCATGGCATGGCTTTCAGCCTGAGCCTCGTCAGCCATGCTCGACAGGTCGCGCAg                                                                                                                                        <  1:4650134/150‑1 (MQ=255)
              ggCAAGGTTCTCAATCAGGCTCGCCAGTTGCGTCGAGGAACGGCTGTTCTCATCCGTCGCCTGGCGAATGTCCACCACCTGCTGCAACGCCGCCTGCATGGCATGGCTTTCAGCCTGAGCCTCGTCAGCCATGCTCGACAGGTCGCGCAg                                                                                                                                        <  2:777053/150‑1 (MQ=255)
                cAAGGTTCTCCATCAGGCTCGCCAGTTGCGTCGAGGAACGGCTGTTCTCATCCGTCGCCTGGCGAATGTCCACCACCTGCTGCAACGCCGCCTGCATGGCACGGCTTTCAGCCTGAGCCTCGTCAGCCATGCGCGACCTGTCGAGCatgg                                                                                                                                       <  2:1724631/150‑4 (MQ=255)
                  aGGTTCTCAATCAGGCTCGCCAGTTGCGTCGAGGAACGGCTGTTCTCATCCGTCGCCTGGCGAATGTCCACCACCTGCTGCAACGCCGCCTGCATGGCATGGCTTTCAGCCTGAGCCTCGTCAGCCATGCTCGACAGGTCGCGCAGGCTg                                                                                                                                    <  2:885157/150‑1 (MQ=255)
                      tctcAATCAGGCTCGCCAGTTGCGTCGAGGAACGGCTGTTCTCATCCGTCGCCTGGCGAATGTCCACCACCTGCTGCAACGCCGCCTGCATGGCATGGCTTTCAGCCTGAGCCTCGCCAGCCTGCGCGACCTGTCGAGCATGGCTgacga                                                                                                                                <  1:1592420/150‑2 (MQ=38)
                        tcAATCAGGCTCGCCAGTTGCGTCGAGGAACGGCTGTTCTCATCCGTCGCCTGGCGAATGTCCACCACCTGCTGCAACGCCGCCTGCATGGCATGGCTTTCAGCCTGAGCCTCGTCAGCCATGCTCGACAGGTCGCGCAGGCTGGCGGcc                                                                                                                              >  2:3839941/1‑150 (MQ=255)
                        tcAATCAGGCTCGCCAGTTGCGTCGAGGAACGGCTGTTCTCATCCGTCGCCTGGCGAATGTCCACCACCTGCTGCAACGCCGCCTGCATGGCATGGCTTTCAGCCTGAGCCTCGTCAGCCATGCGCGACCTGTCGAGCatg                                                                                                                                        <  1:404339/141‑3 (MQ=255)
                          aaTCAGGCTCGCCAGTTGCGTCGAGGAACGGCTGTTCTCATCCGTCGCCTGGCGAATGTCCACCACCTGCTGCAACGCCGCCTGCATGGCATGGCTTTCAGCCTGAGCCTCGTCAGCCATGCTCGACAGGTCGCGCAGGCTGGCGGCCAc                                                                                                                            >  2:2309394/1‑150 (MQ=255)
                          aaTCAGGCTCGCCAGTTGCGTCGAGGAACGGCTGTTCTCATCCGTCGCCTGGCGAATGTCCACCACCTGCTGCAACGCCGCCTGCATGGCATGGCTTTCAGCCTGAGCCTCGTCAGCCATGCTCGACAGGTCGCGCAGGCTGGCGGCCAc                                                                                                                            >  1:3966571/1‑150 (MQ=255)
                            tCAGGCTCGCCAGTTGCGTCGAGGAACGGCTGTTCTCATCCGTCGCCTGGCGAATGTCCACCACCTGCTGCAACGCCGCCTGCATGGCATGGCTTTCAGCCTGAGCCTCGTCAGCCATGCTCGACAGGTCGCGCAGGCTGGCGGCCACtt                                                                                                                          >  1:3860033/1‑150 (MQ=255)
                              aGGCTCGCCAGTTGCGTCGAGGAACGGCTGTTCTCATCCGTCGCCTGGCGAATGTCCACCACCTGCTGCAACGCCGCCTGCATGGCATGGCTTTCAGCCTGAGCCTCGTCAGCCATGCTCGACAGGTCGCGCAGGCTGGCGGCCACTTCg                                                                                                                        <  2:4221477/150‑1 (MQ=255)
                                gCTCGCCAGTTGCGTCGAGGAACGGCTGTTCTCATCCGTCGCCTGGCGAATGTCCACCACCTGCTGCAACGCCGCCTGCATGGCATGGCTTTCAGCCTGAGCCTCGTCAGCCATGCTCGACAGGTCGCGCAGGCTGGCGGCCACTTCGTc                                                                                                                      >  1:3574912/1‑150 (MQ=255)
                                gCTCGCCAGTTGCGTCGAGGAACGGCTGTTCTCATCCGTCGCCTGGCGAATGTCCACCACCTGCTGCAACGCCGCCTGCATGGCATGGCTTTCAGCCTGAGCCTCGTCAGCCATGCTCGACAGGTCGCGCAGGCTGGCGGCCACTTCGTc                                                                                                                      >  2:3796017/1‑150 (MQ=255)
                                gCTCGCCAGTTGCGTCGAGGAACGGCTGTTCTCATCCGTCGCCTGGCGAATGTCCACCACCTGCTGCAACGCCGCCTGCATGGCATGGCTTTCAGCCTGAGCCTCGTCAGCCATGCTCGACAGGTCGCGCAGGCTGGCGGCCACTTCGTc                                                                                                                      >  1:539700/1‑150 (MQ=255)
                                gCTCGCCAGTTGCGTCGAGGAACGGCTGTTCTCATCCGTCGCCTGGCGAATGTCCACCACCTGCTGCAACGCCGCCTGCATGGCATGGCTTTCAGCCTGAGCCTCGTCAGCCATGCTCGACAGGTCGCGCAGGCTGGCGGCCACTTCGTc                                                                                                                      >  1:4240068/1‑150 (MQ=255)
                                  tCGCCAGTTGCGTCGAGGAACGGCTGTTCTCATCCGTCGCCTGGCGAATGTCCACCACCTGCTGCAACGCCGCCTGCATGGCATGGCTTTCAGCCTGAGCCTCGTCAGCCATGCTCGACAGGTCGCGCAGGCt                                                                                                                                     >  2:4671292/1‑133 (MQ=255)
                                      cAGTTGCGTCGAGGAACGGCTGTTCTCATCCGTCGCCTGGCGAATGTCCACCACCTGCTGCAACGCCGCCTGCATGGCATGGCTTTCAGCCTGAGCCTCGTCAGCCATGCTCGACAGGTCGCGCAGGCTGGCGGCCACTTCGTCACGCTg                                                                                                                <  1:1499539/150‑1 (MQ=255)
                                         ttGCGTCGAGGAACGGCTGTTCTCATCCGTCGCCTGGCGAATGTCCACCACCTGCTGCAACGCCGCCTGCATGGCATGGCTTTCAGCCTGAGCCTCGTCAGCCATGCTCGACAGGTCGCGCAGGCTGGCGGCCACTTCGTCACGCTGCaa                                                                                                             >  1:748812/1‑150 (MQ=255)
                                         ttGCGTCGAGGAACGGCTGTTCTCATCCGTCGCCTGGCGAATGTCCACCACCTGCTGCAACGCCGCCTGCATGGCATGGCTTTCAGCCTGAGCCTCGTCAGCCATGCTCGACAGGTCGCGCAGGCTGGCGGCCACTTCGTCACGCTGCaa                                                                                                             <  1:5170109/150‑1 (MQ=255)
                                          tGCGTCGAGGAACGGCTGTTCTCATCCGTCGCCTGGCGAATGTCCACCACCTGCTGCAACGCCGCCTGCATGGCATGGCTTTCAGCCTGAGCCTCGTCAGCCATGCTCGACAGGTCGCGCAGGCTGGCGGCCACTTCGTCACGCTGCAAc                                                                                                            <  2:4048953/150‑1 (MQ=255)
                                                  ggAACGGCTGTTCTCATCCGTCGCCTGGCGAATGTCCACCACCTGCTGCAACGCCGCCTGCATGGCATGGCTTTCAGCCTGAGCCTCGTCAGCCATGCTCGACAGGTCGCGCAGGCTGGCGGCCACTTCGTCACGCTGCAAcgccgccgc                                                                                                    >  1:910352/1‑150 (MQ=255)
                                                   gAACGGCTGTTCTCATCCGTCGCCTGGCGAATGTCCACCACCTGCTGCAACGCCGCCTGCATGGCATGGCTTTCAGCCTGAGCCTCGTCAGCCATGCTCGACAGGTCGCGCAGGCTGGCGGCCACGTCGTCACGCTGCAACGCCGCcgcg                                                                                                   >  1:1148624/1‑150 (MQ=255)
                                                       ggCTGTTCTCATCCGTCGCCTGGCGAATGTCCACCACCTGCTGCAACGCCGCCTGCATGGCATGGCTTTCAGCCTGAGCCTCGTCAGCCATGCTCGACAGGTCGCGCAGGCTGGCGGCCACTTCGTCACGCTGCAACGCCGCCGCGGCAt                                                                                               <  2:529483/150‑1 (MQ=255)
                                                       ggCTGTTCTCATCCGTCGCCTGGCGAATGTCCACCACCTGCTGCAACGCCGCCTGCATGGCATGGCTTTCAGCCTGAGCCTCGTCAGCCATGCTCGACAGGTCGCGCAGGCTGGCGGCCACTTCGTCACGCTGCAACGCCGCCGCGGCAt                                                                                               >  2:5260008/1‑150 (MQ=255)
                                                         cTGTTCTCATCCGTCGCCTGGCGAATGTCCACCACCTGCTGCAACGCCGCCTGCATGGCATGGCTTTCAGCCTGAGCCTCGTCAGCCATGCTCGACAGGTCGCGCAGGCTGGCGGCCACTTCGTCACGCTGCAACGCCGCCGCGGCATCg                                                                                             <  2:3049308/150‑1 (MQ=255)
                                                         cTGTTCTCATCCGTCGCCTGGCGAATGTCCACCACCTGCTGCAACGCCGCCTGCATGGCATGGCTTTCAGCCTGAGCCTCGTCAGCCATGCTCGACAGGTCGCGCAGGCTGGCGGCCACTTCGTCACGCTGCAACGCCGCCGCGGCATCg                                                                                             <  2:3966571/150‑1 (MQ=255)
                                                          tGTTCTCATCCGTCGCCTGGCGAATGTCCACCACCTGCTGCAACGCCGCCTGCATGGCATGGCTTTCAGCCTGAGCCTCGTCAGCCATGCTCGACAGGTCGCGCAGGCTGGCGGCCACTTCGTCACGCTGCAACGCCGCCGCGGCATCgg                                                                                            >  2:5054146/1‑150 (MQ=255)
                                                              ctcATCCGTCGCCTGGCGAATGTCCACCACCTGCTGCAACGCCGCCTGCATGGCATGGCTTTCAGCCTGAGCCTCGTCAGCCATGCTCGACAGGTCGCGCAGGCTGGCGGCCACTTCGTCACGCTGCAACGCCGCCGCGGCATCGGCAcc                                                                                        <  2:935859/150‑1 (MQ=255)
                                                               tcATCCGTCGCCTGGCGAATGTCCACCACCTGCTGCAACGCCGCCTGCATGGCATGGCTTTCAGCCTGAGCCTCGTCAGCCATGCTCGACAGGTCGCGCAGGCTGGCGGCCACTTCGTCACGCTGCAACGCCGCCGCGGCATCGGCACCg                                                                                       <  2:2246165/150‑1 (MQ=255)
                                                               tcATCCGTCGCCTGGCGAATGTCCACCACCTGCTGCAACGCCGCCTGCATGGCATGGCTTTCAGCCTGAGCCTCGTCAGCCATGCTCGACAGGTCGCGCAGGCTGGCGGCCACTTCGTCACGCTGCAACGCCGCCGCGGCATCGGCACCg                                                                                       >  1:598180/1‑150 (MQ=255)
                                                               tcATCCGTCGCCTGGCGAATGTCCACCACCTGCTGCAACGCCGCCTGCATGGCATGGCTTTCAGCCTGAGCCTCGTCAGCCATGCTCGACAGGTCGCGCAGGCTGGCGGCCACTTCGTCACGCTGCAACGCCGCCGCGGCATCGGCACCg                                                                                       >  1:1496077/1‑150 (MQ=255)
                                                               tcATCCGTCGCCTGGCGAATGTCCACCACCTGCTGCAACGCCGCCTGCATGGCATGGCTTTCAGCCTGAGCCTCGTCAGCCATGCTCGACAGGTCGCGCAGGCTGGCGGCCACTTCGTCACGCTGCAACGCCGCCGCGGCATCGGCACCg                                                                                       <  1:5359876/150‑1 (MQ=255)
                                                                cATCCGTCGCCTGGCGAATGTCCACCACCTGCTGCAACGCCGCCTGCATGGCATGGCTTTCAGCCTGAGCCTCGTCAGCCATGCTCGACAGGTCGCGCAGGCTGGCGGCCACTTCGTCACGCTGCAACGCCGCCGCGGCATCGGCACCgg                                                                                      >  1:2611293/1‑150 (MQ=255)
                                                                 aTCCGTCGCCTGGCGAATGTCCACCACCTGCTGCAACGCCGCCTGCATGGCATGGCTTTCAGCCTGAGCCTCGTCAGCCATGCTCGACAGGTCGCGCAGGCTGGCGGCCACTTCGTCACGCTGCAACGCCGCCGCGGCATCGGCACCGGc                                                                                     >  1:5442014/1‑150 (MQ=255)
                                                                   ccGTCGCCTGGCGAATGTCCACCACCTGCTGCAACGCCGCCTGCATGGCATGGCTTTCAGCCTGAGCCTCGTCAGCCATGCTCGACAGGTCGCGCAGGCTGGCGGCCACTTCGTCACGCTGCAACGCCGCCGCGGCATCGGCACCGGCAt                                                                                   >  1:2752026/1‑150 (MQ=255)
                                                                    cGTCGCCTGGCGAATGTCCACCACCTGCTGCAACGCCGCCTGCATGGCATGGCTTTCAGCCTGAGCCTCGTCAGCCATGCTCGACAGGTCGCGCAGGCTGGCGGCCACTTCGTCACGCTGCAACGCCGCCGCGGCATCGGCACCGGCAtt                                                                                  <  2:3720062/150‑1 (MQ=255)
                                                                          cTGGCGAATGTCCACCACCTGCTGCAACGCCGCCTGCATGGCATGGCTTTCAGCCTGAGCCTCGTCAGCCATGCTCGACAGGTCGCGCAGGCTGGCGGCCACTTCGTCACGCTGCAACGCCGCCGCGGCATCGGCACCGGCATTGCGCTg                                                                            >  1:3730691/1‑150 (MQ=255)
                                                                           tGGCGAATGTCCACCACCTGCTGCAACGCCGCCTGCATGGCATGGCTTTCAGCCTGAGCCTCGTCAGCCAGGCTCGACAGGTCGCGCAGGCTGGCGGCCACTTCGTCACGCTGCAACGCCGCCGCGGCATCGGCACCGGCATTGCGCTgg                                                                           <  2:1834149/150‑1 (MQ=255)
                                                                               gAATGTCCACCACCTGCTGCAACGCCGCCTGCATGGCATGGCTTTCAGCCTGAGCCTCGTCAGCCATGCTCGACAGGTCGCGCAGGCTGGCGGCCACTTCGTCACGCTGCAACGCCGCCGCGGCATCTGCACCGGCATTGCGCTGGGCCa                                                                       >  2:787333/1‑150 (MQ=255)
                                                                               gAATGTCCACCACCTGCTGCAACGCCGCCTGCATGGCATGGCTTTCAGCCTGAGCCTCGTCAGCCATGCTCGACAGGTCGCGCAGGCTGGCGGCCACTTCGTCACGCTGCAACGCCGCCGCGGCATCGGCACCGGCATTGCGCTGGGCCa                                                                       >  2:1018740/1‑150 (MQ=255)
                                                                               gAATGTCCACCACCTGCTGCAACGCCGCCTGCATGGCATGGCTTTCAGCCTGAGCCTCGTCAGCCATGCTCGACAGGTCGCGCAGGCTGGCGGCCACTTCGTCACGCTGCAACGCCGCCGCGGCATCGGCACCGGCATTGCGCTGGGCCa                                                                       >  2:1997588/1‑150 (MQ=255)
                                                                                  tGTCCACCACCTGCTGCAACGCCGCCTGCATGGCATGGCTTTCAGCCTGAGCCTCGTCAGCCATGCTCGACAGGTCGCGCAGGCTGGCGGCCACTTCGTCACGCTGCAACGCCGCCGCGGCATCGGCACCGGCATTGCGCTGGGCCATgg                                                                    <  1:1475827/150‑1 (MQ=255)
                                                                                  tGTCCACCACCTGCTGCAACGCCGCCTGCATGGCATGGCTTTCAGCCTGAGCCTCGTCAGCCATGCTCGACAGGTCGCGCAGGCTGGCGGCCACTTCGTCACGCTGCAACGCCGCCGCGGCATCGGCACCGGCATTGCGCTGGGCCATgg                                                                    >  2:669522/1‑150 (MQ=255)
                                                                                      caccacCTGCTGCAACGCCGCCTGCATGGCATGGCTTTCAGCCTGAGCCTCGTCAGCCATGCTCGACAGGTCGCGCAGGCTGGCGGCCACTTCGTCACGCTGCAAGGCCGCCGCGGCATCGGCACCGGCATTGCGCTGGGCCATGGCGcc                                                                <  1:2025803/150‑1 (MQ=255)
                                                                                      caccacCTGCTGCAACGCCGCCTGCATGGCATGGCTTTCAGCCTGAGCCTCGTCAGCCATGCTCGACAGGTCGCGCAGGCTGGCGGCCACTTCGTCACGCTGCAACGCCGCCGCGGCATCGGCACCGGCATTGCGCTGGGCCATGGCGcc                                                                >  2:5382150/1‑150 (MQ=255)
                                                                                       accaccTGCTGCAACGCCGCCTGCATGGCATGGCTTTCAGCCTGAGCCTCGTCAGCCATGCTCGACAGGTCGCGCAGGCTGGCGGCCACTTCGTCACGCTGCAACGCCGCCGCGGCATCGGCACCGGCATTGCGCTGGGCCATGGCGCCg                                                               >  1:2377685/1‑150 (MQ=255)
                                                                                        ttttcagctgcAACGCCGCCTGCATGGCATGGCTTTCAGCCTGAGCCTCGTCAGCCATGCTCGACAGGTCGCGCAGGCTGGCGGCCACTTCGTCACGCTGCAACGCCGCCGCGGCATCGGCACCGGCATTGCGCTGGGCCATGGCGCCGa                                                              >  2:584466/7‑150 (MQ=255)
                                                                                            ctgctgCAACGCCGCCTGCATGGCATGGCTTTCAGCCTGAGCCTCGTCAGCCATGCTCGACAGGTCGCGCAGGCTGGCGGCCACTTCGTCACGCTGCAACGCCGCCGCGGCATCGGCACCGGCATTGCGCTGGGCCATGGCGCCGAtctc                                                          <  2:3071351/150‑1 (MQ=255)
                                                                                             tgctgcAACGCCGCCTGCATGGCATGGCTTTCAGCCTGAGCCTCGTCAGCCATGCTCGACAGGTCGCGCAGGCTGGCGGCCACTTCGTCACGCTGCAACGCCGCCGCGGCATCGGCACCGGCATTGCGCTGGGCCATGGCGCCGATCTCg                                                         >  1:3074925/1‑150 (MQ=255)
                                                                                             tgctgcAACGCCGCCTGCATGGCATGGCTTTCAGCCTGAGCCTCGTCAGCCATGCTCGACAGGTCGCGCAGGCTGGCGGCCACTTCGTCACGCTGCAACGCCGCCGCGGCATCGGCACCGGCATTGCGCTGGGCCATGGCGCCGATCTCg                                                         <  2:2723781/150‑1 (MQ=255)
                                                                                              gctgcAACGCCGCCTGCATGGCATGGCTTTCAGCCTGAGCCTCGTCAGCCATGCTCGACAGGTCGCGCAGGCTGGCGGCCACTTCGTCACGCTGCAACGCCGCCGCGGCATCGGCACCGGCATTGCGCTGGGCCATGGCGCCGATCTCGa                                                        <  2:1496077/150‑1 (MQ=255)
                                                                                               ctgcAACGCCGCCTGCATGGCATGGCTTTCAGCCTGAGCCTCGTCAGCCATGCTCGACa                                                                                                                                                   <  1:2876776/59‑1 (MQ=255)
                                                                                               ctgcAACGCCGCCTGCATGGCATGGCTTTCAGCCTGAGCCTCGTCAGCCATGCTCGACAGGTCGCGCAGGCTGGCGGCCACTTCGTCACGCTGCAACGCCGCCGCGGCATCGGCACCGGCATTGCGCTGGGCCATGGCGCCGATCTCGAc                                                       >  2:1218853/1‑150 (MQ=255)
                                                                                               ctgcAACGCCGCCTGCATGGCATGGCTTTCAGCCTGAGCCTCGTCAGCCATGCTCGACAGGTCGCGCAGGCTGGCGGCCACTTCGTCACGCTGCAACGCCGCCGCGGCATCGGCACCGGCATTGCGCTGGGCCATGGCGCCGATCTCGAc                                                       >  2:4121891/1‑150 (MQ=255)
                                                                                                    aCGCCGCCTGCATGGCATGGCTTTCAGCCTGAGCCTCGTCAGCCATGCTCGACAGGTCGCGCAGGCTGGCGGCCACTTCGTCACGCTGCAACGCCGCCGCGGCATCGGCACCGGCATTGCGCTGGGCCATGGCGCCGATCTCGACGCCgg                                                  >  2:157869/1‑150 (MQ=255)
                                                                                                     cgccgcCTGCATGGCATGGCTTTCAGCCTGAGCCTCGTCAGCCATGCTCGACAGGTCGCGCAGGCTGGCGGCCACTTCGTCACGCTGCAACGCCGCCGCGGCATCGGCACCGGCATTGCGCTGGGCCATGGCGCCGATCTCGACGCCGGt                                                 >  1:511495/1‑150 (MQ=255)
                                                                                                            tGCATGGCATGGCTTTCAGCCTGAGCCTCGTCAGCCATGCTCGACAGGTCGCGCAGGCTGGCGGCCACTTCGTc                                                                                                                      <  1:2235445/74‑1 (MQ=255)
                                                                                                            tGCATGGCATGGCTTTCAGCCTGAGCCTCGTCAGCCATGCTCGACAGGTCGCGCAGGCTGGCGGCCACTTCGTCACGCTGCAACGCCGCCGCGGCATCGGCACCGGCATTGCGCTGGGCCATGGCGCCGATCTCGACGCCGGTACGCTgg                                          <  2:144729/150‑1 (MQ=255)
                                                                                                               atggcatggcTTTCAGCCTGAGCCTCGTCAGCCATGCTCGACAGGTCGCGCAGGCTGGCGGCCACTTCGTCACGCTGCAACGCCGCCGCGGCATCGGCACCGGCATTGCGCTGGGCCATGGCGCCGATCTCGACGCCGGTACGCTgg                                          >  1:1348001/1‑147 (MQ=255)
                                                                                                                    atggcTTTCAGCCTGAGCCTCGTCAGCCATGCTCGACAGGTCGCGCAGGCTGGCGGCCACTTCGTCACGCTGCAACGCCGCCGCGGCATCGGCACCGGCATTGCGCTGGGCCATGGCGCCGATCTCGACGCCGGTACGCTGGGCCACCTc                                  >  2:3201710/1‑150 (MQ=255)
                                                                                                                       gcTTTCAGCCTGAGCCTCGTCAGCCATGCTCGACAGGTCGCGCAGGCTGGCGGCCACTTCGTCACGCTGCAACGCCGCCGCGGCATCGGCACCGGCATTGCGCTGGGCCATGGCGCCGATCTCGACGCCGGTACGCTGGGCCACCTCGcc                               <  2:2886650/150‑1 (MQ=255)
                                                                                                                        cTTTCAGCCTGAGCCTCGTCAGCCATGCTCGACAGGTCGCGCAGGCTGGCGGCCACTTCGTCACGCTGCAACGCCGCCGCGGCATCGGCACCGGCATTGCGCTGGGCCATGGCGCCGATCTCGACGCCGGTACGCTGGGCCACCTCGCCg                              >  2:45915/1‑150 (MQ=255)
                                                                                                                        cTTTCAGCCTGAGCCTCGTCAGCCATGCTCGACAGGTCGCGCAGGCTGGCGGCCACTTCGTCACGCTGCAACGCCGCCGCGGCATCGGCACCGGCATTGCGCTGGGCCATGGCGCCGATCTCGACGCCGGTACGCTGGGCCACCTCGCCg                              >  2:4107571/1‑150 (MQ=255)
                                                                                                                        cTTTCAGCCTGAGCCTCGTCAGCCATGCTCGACAGGTCGCGCAGGCTGGCGGCCACTTCGTCACGCTGCAACGCCGCCGCGGCATCGGCACCGGCATTGCGCTGGGCCATGGCGCCGATCTCGACGCCGGTACGCTGGGCCACCTCGCCg                              <  2:3488644/150‑1 (MQ=255)
                                                                                                                         tttCAGCCTGAGCCTCGTCAGCCATGCTCGACAGGTCGCGCAGGCTGGCGGCCACTTCGTCACGCTGCAACGCCGCCGCGGCATCGGCACCGGCATTGCGCTGGGCCATGGCGCCGATCTCGACGCCGGTACGCTGGGCCACCTCGCCgg                             <  1:4387326/150‑1 (MQ=255)
                                                                                                                           tCAGCCTGAGCCTCGTCAGCCATGCTCGACAGGTCGCGCAGGCTGGCGGCCACTTCGTCACGCTGCAACGCCGCCGCGGCATCGGCACCGGCATTGCGCTGGGCCATGGCGCCGATCTCGACGCCGGTACGCTGGGCCACCTCGCCGGcc                           >  2:1432591/1‑150 (MQ=255)
                                                                                                                             aGCCTGAGCCTCGTCAGCCATGCTCGACAGGTCGCGCAGGCTGGCGGCCACTTCGTCACGCTGCAACGCCGCCGCGGCATCGGCACCGGCATTGCGCTGGGCCATGGCGCCGATCTCGACGCCGGTACGCTGGGCCACCTCGCCGGCCTc                         >  1:1673678/1‑150 (MQ=255)
                                                                                                                              gCCTGAGCCTCGTCAGCCATGCTCGACAGGTCGCGCAGGCTGGCGGCCACGTCGTCACGCTGCAACGCCGCCGCGGCATCGGCACCGGCATTGCGCTGGGCCATGGCGCCGATCTCGACGCCGGTACGCTGGGCCACCTCGCCGGCCTcg                        >  2:2945349/1‑150 (MQ=255)
                                                                                                                                 tGAGCCTCGTCAGCCATGCTCGACAGGTCGCGCAGGCTGGCGGCCACTTCGTCACGCTGCAACGCCGCCGCGGCATCGGCACCGGCATTGCGCTGGGCCATGGCGCCGATCTCGACGCCGGTACGCTGGGCCACCTCGCCGGCCTCGCgg                     <  2:3574912/150‑1 (MQ=255)
                                                                                                                                   aGCCTCGTCAGCCATGCTCGACAGGTCGCGCAGGCTGGCGGCCACTTCGTCACGCTGCAACGCCGCCGCGGCATCGGCACCGGCATTGCGCTGGGCCATGGCGCCGATCTCGACGCCGGTACGCTGGGCCACCTCGCCGGCCTCGCGGAc                   <  2:2096203/150‑1 (MQ=255)
                                                                                                                                       tCGTCAGCCATGCTCGACAGGTCGCGCAGGCTGGCGGCCACTTCGTCACGCTGCAACGCCGCCGCGGCATCGGCACCGGCATTGCGCTGGGCCATGGCGCCGATCTCGACGCCGGTACGCTGGGCCACCTCGCCGGCCTCGCGGACGATg               <  2:1868703/150‑1 (MQ=255)
                                                                                                                                       tCGTCAGCCATGCTCGACAGGTCGCGCAGGCTGGCGGCCACTTCGTCACGCTGCAACGCCGCCGCGGCATCGGCACCGGCATTGCGCTGGGCCATGGCGCCGATCTCGACGCCGGTACGCTGGGCCACCTCGCCGGCCTCGCGGACGATg               <  1:1576315/150‑1 (MQ=255)
                                                                                                                                        cGTCAGCCATGCTCGACAGGTCGCGCAGGCTGGCGGCCACTTCGTCACGCTGCAACGCCGCCGCGGCATCGGCACCGGCATTGCGCTGGGCCATGGCGCCGATCTCGACGCCGGTACGCTGGGCCACCTCGCCGGCCTCGCGGACGATgg              <  2:4261715/150‑1 (MQ=255)
                                                                                                                                        cGTCAGCCATGCTCGACAGGTCGCGCAGGCTGGCGGCCACTTCGTCACGCTGCAACGCCGCCGCGGCATCGGCACCGGCATTGCGCTGGGCCATGGCGCCGATCTCGACGCCGGTACGCTGGGCCACCTCGCCGGCCTCGCGGACGATgg              <  1:5030107/150‑1 (MQ=255)
                                                                                                                                         gTCAGCCATGCTCGACAGGTCGCGCAGGCTGGCGGCCACTTCGTCACGCTGCAACGCCGCCGCGGCATCGGCACCGGCATTGCGCTGGGCCATGGCGCCGATCTCGACGCCGGTACGCTGGGCCACCTCGCCGGCCTCGCGGACGATggg             <  2:26514/150‑1 (MQ=255)
                                                                                                                                          tCAGCCATGCTCGACAGGTCGCGCAGGCTGGCGGCCACTTCGTCACGCTGCAACGCCGCCGCGGCATCGGCACTGGCATTGCGCTGGGCCATGGCGCCGATCTCGACGCCGGTACGCTGGGCCACCTCGCCGGCCTCGCGGACGATGGGc            >  1:1491297/1‑150 (MQ=255)
                                                                                                                                           cAGCCATGCTCGACAGGTCGCGCAGGCTGGCGGCCACTTCGTCACGCTGCAACGCCGCCGCGGCATCGGCACCGGCATTGCGCTGGGCCATGGCGCCGATCTCGACGCCGGTACGCTGGGCCACCTCGCCGGCCTCGCGGACGATGGGCt           <  1:971195/150‑1 (MQ=255)
                                                                                                                                            aGCCATGCTCGACAGGTCGCGCAGGCTGGCGGCCACTTCGTCACGCTGCAACGCCGCCGCGGCATCGGCACCGGCATTGCGCTGGGCCATGGCGCCGAtctc                                                          <  1:5312634/102‑1 (MQ=255)
                                                                                                                                            aGCCATGCTCGACAGGTCGCGCAGGCTGGCGGCCACTTCGTCACGCTGCAACGCCGCCGCGGCATCGGCACCGGCATTGCGCTGGGCCATGGCGCCGAtctc                                                          >  2:5312634/1‑102 (MQ=255)
                                                                                                                                            aGCCATGCTCGACAGGTCGCGCAGGCTGGCGGCCACTTCGTCACGCTGCAACGCCGCCGCGGCATCGGCACCGGCATTGCGCTGGGCCATGGCGCCGATCTCGACGCCGGTACGCTGGGCCACCTCTCCGGCCTCGCGGACGATGGGCTg          <  2:1190964/150‑1 (MQ=255)
                                                                                                                                            aGCCATGCTCGACAGGTCGCGCAGGCTGGCGGCCACTTCGTCACGCTGCAACGCCGCCGCGGCATCGGCACCGGCATTGCGCTGGGCCATGGCGCCGATCTCGACGCCGGTACGCTGGGCCACCTCGCCGGCCTCGCGGACGATGGGCTg          >  2:574360/1‑150 (MQ=255)
                                                                                                                                              gcctGCGCGACCTGTCGAGCATGGCTGACGGCCACTTCGTCACGCTGCAACGCCGCCGCGGCATCGGCACCGGCATTGCGCTGGGCCATGGCGCCGATCTCGACGCCGGTACGCTGGGCCACCTCGCCGGCCTCGCGGACGATGGGCTGc         >  2:4373309/4‑150 (MQ=255)
                                                                                                                                              ccATGCTCGACAGGTCGCGCAGGCTGGCGGCCACTTCGTCACGCTGCAACGCCGCCGCGGCATCGGCACCGGCATTGCGCTGGGCCATGGCGCCGATCTCGACGCCGGTACGCTGGGCCACCTCGCCGGCCTCGCGGACGATGGGCTGCa        <  2:2679922/150‑1 (MQ=255)
                                                                                                                                               cctGCGCGACCTGTCGAGCATGCTGGCGGCCACTTCGTCACGCTGCAACGCCGCCGCGGCATCGGCACCGGCATTGCGCTGGGCCATGGCGCCGATCTCGACGCCGGTACGCTGGGCCACCTCGCCGGCCTCGCGGACGATGGGCTGCAg       >  2:2136992/3‑150 (MQ=255)
                                                                                                                                               cctGCGCGACCTGTCGAGCAGGCTGGCGGCCACTTCGTCACGCTGCAACGCCGCCGCGGCATCGGCACCGGCATTGCGCTGGGCCATGGCGCCGATCTCGACGCCGGTACGCTGGGCCACCTCGCCGGCCTCGCGGACGATGGGCTGCAg       >  1:4851020/3‑150 (MQ=255)
                                                                                                                                               cATGCTCGACAGGTCGCGCAGGCTGGCGGCCACTTCGTCACGCTGCAACGCCGCCGCGGCATCGGCACCGGCATTGCGCTGGGCCATgg                                                                    <  1:1881759/89‑1 (MQ=255)
                                                                                                                                                aTGCTCGACAGGTCGCGCAGGCTGGCGGCCACTTCGTCACGCTGCAACGCCGCCGCGGCATCGGCACCGGCATTGCGCTGGGCCATGGCGCCGATCTCGACGCCGGTACGCTGGGCCACCTCGCCGGCCTCGCGGACGATGGGCTGCAGc      <  2:2861730/150‑1 (MQ=255)
                                                                                                                                                 tGCTCGACAGGTCGCGCAGGCTGGCGGCCACTTCGTCACGCTGCAACGCCGCCGCGGCATCGGCACCGGCATTGCGCTGGGCCATGGCGCCGATCTCGACGCCGGTACGCTGGGCCACCTCGCCGGCCTCGCGGACGATggg             >  1:4786778/1‑142 (MQ=255)
                                                                                                                                                 tGCTCGACAGGTCGCGCAGGCTGGCGGCCACTTCGTCACGCTGCAACGCCGCCGCGGCATCGGCACCGGCATTGCGCTGGGCCATGGCGCCGATCTCGACGCCGGTACGCTGGGCCACCTCGCCGGCCTCGCGGACGATggg             <  2:4786778/142‑1 (MQ=255)
                                                                                                                                                 tGCGCGACCTGTCGAGCAGGCTGGCGGCCACTTCGTCACGCTGCAACGCCGCCGCGGCATCGGCACCGGCATTGCGCTGGGCCATGGCGCCGATCTCGACGCCGGTACGCTGGGCCACCTCGCCGGCCTCGCGGACGATGGGCTGCAGCt     >  2:1275855/1‑150 (MQ=255)
                                                                                                                                                    |                                                                                                                                                      
GTTTCAACCTGCCCGGCAAGGTTCTCAATCAGGCTCGCCAGTTGCGTCGAGGAACGGCTGTTCTCATCCGTCGCCTGGCGAATGTCCACCACCTGCTGCAACGCCGCCTGCATGGCATGGCTTTCAGCCTGAGCCTCGTCAGCCATGCTCGACAGGTCGCGCAGGCTGGCGGCCACTTCGTCACGCTGCAACGCCGCCGCGGCATCGGCACCGGCATTGCGCTGGGCCATGGCGCCGATCTCGACGCCGGTACGCTGGGCCACCTCGCCGGCCTCGCGGACGATGGGCTGCAGCTTGT  >  NC_002947/5558227‑5558524

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