Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 5,558,386 C→A 5.8% R406L (CGC→CTC)  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,3860CA 5.8% 658.9 / 19.0 197R406L (CGC→CTC) PP_4888methyl‑accepting chemotaxis transducer
Reads supporting (aligned to +/- strand):  ref base C (98/87);  new base A (7/5);  total (105/92)
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 213 total aligned reads displayed.

GCCCGGCAAGGTTCTCAATCAGGCTCGCCAGTTGCGTCGAGGAACGGCTGTTCTCATCCGTCGCCTGGCGAATGTCCACCACCTGCTGCAACGCCGCCTGCATGGCATGGCTTTCAGCCTGAGCCTCGTCAGCCATGCTCGACAGGTCGCGCAGGCTGGCGGCCACTTCGTCACGCTGCAACGCCGCCGCGGCATCGGCACCGGCATTGCGCTGGGCCATGGCGCCGATCTCGACGCCGGTACGCTGGGCCACCTCGCCGGCCTCGCGGACGATGGGCTGCAGCTTGTCGACAAAGCGA  >  NC_002947/5558237‑5558535
                                                                                                                                                     |                                                                                                                                                      
gCCCGGCAAGGTTCTCAATCAGGCTCGCCAGTTGCGTCGAGGAACGGCTGTTCTCATCCGTCGCCTGGCGAATGTCCACCACCTGCTGCAACGCCGCCTGCATGGCATGGCTTTCAGCCTGAGCCTCGTCAGCCATGCTCGACAGGTcgc                                                                                                                                                        <  1:2062947/150‑1 (MQ=255)
gCCCGGCAAGGTTCTCAATCAGGCTCGCCAGTTGCGTCGAGGAACGGCTGTTCTCATCCGTCGCCTGGCGAATGTCCACCACCTGCTGCAACGCCGCCTGCATGGCATGGCTTTCAGCCTGAGCCTCGTCAGCCATGCTCGACAGGTcgc                                                                                                                                                        >  2:415967/1‑150 (MQ=255)
    ggCAAGGTTCTCAATCAGGCTCGCCAGTTGCGTCGAGGAACGGCTGTTCTCATCCGTCGCCTGGCGAATGTCCACCACCTGCTGCAACGCCGCCTGCATGGCATGGCTTTCAGCCTGAGCCTCGTCAGCCATGCTCGACAGGTCGCGCAg                                                                                                                                                   <  1:3716362/150‑1 (MQ=255)
      cAAGGTTCTCCATCAGGCTCGCCAGTTGCGTCGAGGAACGGCTGTTCTCATCCGTCGCCTGGCGAATGTCCACCACCTGCTGCAACGCCGCCTGCATGGCACGGCTTTCAGCCTGAGCCTCGTCAGCCATGCGCGACCTGTCGAGCatgg                                                                                                                                                  <  2:1724631/150‑4 (MQ=255)
        aGGTTCTCAATCAGGCTCGCCAGTTGCGTCGAGGAACGGCTGTTCTCATCCGTCGCCTGGCGAATGTCCACCACCTGCTGCAACGCCGCCTGCATGGCATGGCTTTCAGCCTGAGCCTCGTCAGCCATGCTCGACAGGTCGCGCAGGCTg                                                                                                                                               <  2:885157/150‑1 (MQ=255)
          gTTCTCAATCAGGCTCGCCAGTTGCGTCGAGGAACGGCTGTTCTCATCCGTCGCCTGGCGAATGTCCACCACCTGCTGCAACGCCGCCTGCATGGCATGGCTTTCAGCCTGAGCCTCGTCAGCCATGCTCGACAGGTCGCGCAGGCTggc                                                                                                                                             >  2:779130/1‑150 (MQ=255)
          gTTCTCAATCAGGCTCGCCAGTTGCGTCGAGGAACGGCTGTTCTCATCCGTCGCCTGGCGAATGTCCACCACCTGCTGCAACGCCGCCTGCATGGCATGGCTTTCAGCCTGAGCCTCGTCAGCCATGCTCGACAGGTCGCGCAGGCTggc                                                                                                                                             >  2:2126124/1‑150 (MQ=255)
            tctcAATCAGGCTCGCCAGTTGCGTCGAGGAACGGCTGTTCTCATCCGTCGCCTGGCGAATGTCCACCACCTGCTGCAACGCCGCCTGCATGGCATGGCTTTCAGCCTGAGCCTCGCCAGCCTGCGCGACCTGTCGAGCATGGCTgacga                                                                                                                                           <  1:1592420/150‑2 (MQ=38)
             ctcAATCAGGCTCGCCAGTTGCGTCGAGGAACGGCTGTTCTCATCCGTCGCCTGGCGAATGTCCACCACCTGCTGCAACGCCGCCTGCATGGCATGGCTTTCAGCCTGAGCCTCGTCAGCCATGCTCGACAGGTCGCGCAGGCTggcggc                                                                                                                                          <  2:1259967/150‑1 (MQ=255)
              tcAATCAGGCTCGCCAGTTGCGTCGAGGAACGGCTGTTCTCATCCGTCGCCTGGCGAATGTCCACCACCTGCTGCAACGCCGCCTGCATGGCATGGCTTTCAGCCTGAGCCTCGTCAGCCATGCTCGACAGGTCGCGCAGGCTGGCGGcc                                                                                                                                         >  1:3140168/1‑150 (MQ=255)
              tcAATCAGGCTCGCCAGTTGCGTCGAGGAACGGCTGTTCTCATCCGTCGCCTGGCGAATGTCCACCACCTGCTGCAACGCCGCCTGCATGGCATGGCTTTCAGCCTGAGCCTCGTCAGCCATGCGCGACCTGTCGAGCatg                                                                                                                                                   >  2:404339/1‑139 (MQ=255)
                aaTCAGGCTCGCCAGTTGCGTCGAGGAACGGCTGTTCTCATCCGTCGCCTGGCGAATGTCCACCACCTGCTGCAACGCCGCCTGCATGGCATGGCTTTCAGCCTGAGCCTCGTCAGCCATGCTCGACAGGTCGCGCAGGCTGGCGGCCAc                                                                                                                                       >  1:4255402/1‑150 (MQ=255)
                aaTCAGGCTCGCCAGTTGCGTCGAGGAACGGCTGTTCTCATCCGTCGCCTGGCGAATGTCCACCACCTGCTGCAACGCCGCCTGCATGGCATGGCTTTCAGCCTGAGCCTCGTCAGCCATGCTCGACAGGTCGCGCAGGCTGGCGGCCAc                                                                                                                                       >  2:2309394/1‑150 (MQ=255)
                      gCTCGCCAGTTGCGTCGAGGAACGGCTGTTCTCATCCGTCGCCTGGCGAATGTCCACCACCTGCTGCAACGCCGCCTGCATGGCATGGCTTTCAGCCTGAGCCTCGTCAGCCATGCTCGACAGGTCGCGCAGGCTGGCGGCCACTTCGTc                                                                                                                                 >  1:539700/1‑150 (MQ=255)
                      gCTCGCCAGTTGCGTCGAGGAACGGCTGTTCTCATCCGTCGCCTGGCGAATGTCCACCACCTGCTGCAACGCCGCCTGCATGGCATGGCTTTCAGCCTGAGCCTCGTCAGCCATGCTCGACAGGTCGCGCAGGCTGGCGGCCACTTCGTc                                                                                                                                 >  1:1493236/1‑150 (MQ=255)
                        tCGCCAGTTGCGTCGAGGAACGGCTGTTCTCATCCGTCGCCTGGCGAATGTCCACCACCTGCTGCAACGCCGCCTGCATGGCATGGCTTTCAGCCTGAGCCTCGTCAGCCATGCTCGACAGGTCGCGCAGGCt                                                                                                                                                >  2:4671292/1‑133 (MQ=255)
                        tCGCCAGTTGCGTCGAGGAACGGCTGTTCTCATCCGTCGCCTGGCGAATGTCCACCACCTGCTGCAACGCCGCCTGCATGGCATGGCTTTCAGCCTGAGCCTCGTCAGCCATGCTCGACAGGTCGCGCAGGCt                                                                                                                                                <  1:4671292/133‑1 (MQ=255)
                            cAGTTGCGTCGAGGAACGGCTGTTCTCATCCGTCGCCTGGCGAATGTCCACCACCTGCTGCAACGCCGCCTGCATGGCATGGCTTTCAGCCTGAGCCTCGTCAGCCATGCTCGACAGGTCGCGCAGGCTGGCGGCCACTTCGTCACGCTg                                                                                                                           <  1:1499426/150‑1 (MQ=255)
                               ttGCGTCGAGGAACGGCTGTTCTCATCCGTCGCCTGGCGAATGTCCACCACCTGCTGCAACGCCGCCTGCATGGCATGGCTTTCAGCCTGAGCCTCGTCAGCCATGCTCGACAGGTCGCGCAGGCTGGCGGCCACTTCGTCACGCTGCaa                                                                                                                        >  1:748812/1‑150 (MQ=255)
                               ttGCGTCGAGGAACGGCTGTTCTCATCCGTCGCCTGGCGAATGTCCACCACCTGCTGCAACGCCGCCTGCATGGCATGGCTTTCAGCCTGAGCCTCGTCAGCCATGCTCGACAGGTCGCGCAGGCTGGCGGCCACTTCGTCACGCTGCaa                                                                                                                        >  1:4035041/1‑150 (MQ=255)
                                  cGTCGAGGAACGGCTGTTCTCATCCGTCGCCTGGCGAATGTCCACCACCTGCTGCAACGCCGCCTGCATGGCATGGCTTTCAGCCTGAGCCTCGTCAGCCATGCTCGACAGGTCGCGCAGGCTGGCGGCCACTTCGTCACGCTGCAAcgc                                                                                                                     <  1:2320971/150‑1 (MQ=255)
                                        ggAACGGCTGTTCTCATCCGTCGCCTGGCGAATGTCCACCACCTGCTGCAACGCCGCCTGCATGGCATGGCTTTCAGCCTGAGCCTCGTCAGCCATGCTCGACAGGTCGCGCAGGCTGGCGGCCACTTCGTCACGCTGCAAcgccgccgc                                                                                                               >  1:910352/1‑150 (MQ=255)
                                             ggCTGTTCTCATCCGTCGCCTGGCGAATGTCCACCACCTGCTGCAACGCCGCCTGCATGGCATGGCTTTCAGCCTGAGCCTCGTCAGCCATGCTCGACAGGTCGCGCAGGCTGGCGGCCACTTCGTCACGCTGCAACGCCGCCGCGGCAt                                                                                                          <  2:5078931/150‑1 (MQ=255)
                                             ggCTGTTCTCATCCGTCGCCTGGCGAATGTCCACCACCTGCTGCAACGCCGCCTGCATGGCATGGCTTTCAGCCTGAGCCTCGTCAGCCATGCTCGACAGGTCGCGCAGGCTGGCGGCCACTTCGTCACGCTGCAACGCCGCCGCGGCAt                                                                                                          <  2:529483/150‑1 (MQ=255)
                                             ggCTGTTCTCATCCGTCGCCTGGCGAATGTCCACCACCTGCTGCAACGCCGCCTGCATGGCATGGCTTTCAGCCTGAGCCTCGTCAGCCATGCTCGACAGGTCGCGCAGGCTGGCGGCCACTTCGTCACGCTGCAACGCCGCCGCGGCAt                                                                                                          >  2:5260008/1‑150 (MQ=255)
                                               cTGTTCTCATCCGTCGCCTGGCGAATGTCCACCACCTGCTGCAACGCCGCCTGCATGGCATGGCTTTCAGCCTGAGCCTCGTCAGCCATGCTCGACAGGTCGCGCAGGCTGGCGGCCACTTCGTCACGCTGCAACGCCGCCGCGGCATCg                                                                                                        <  2:3049308/150‑1 (MQ=255)
                                                tGTTCTCATCCGTCGCCTGGCGAATGTCCACCACCTGCTGCAACGCCGCCTGCATGGCATGGCTTTCAGCCTGAGCCTCGTCAGCCATGCTCGACAGGTCGCGCAGGCTGGCGGCCACTTCGTCACGCTGCAACGCCGCCGCGGCATCgg                                                                                                       <  1:396226/150‑1 (MQ=255)
                                                 gTTCTCATCCGTCGCCTGGCGAATGTCCACCACCTGCTGCAACGCCGCCTGCATGGCATGGCTTTCAGCCTGAGCCTCGTCAGCCATGCTCGACAGGTCGCGCAGGCTGGCGGCCACTTCGTCACGCTGCAACGCCGCCGCGGCATCGGc                                                                                                      >  1:4048907/1‑150 (MQ=255)
                                                  ttCTCATCCGTCGCCTGGCGAATGTCCACCACCTGCTGCAACGCCGCCTGCATGGCATGGCTTTCAGCCTGAGCCTCGTCAGCCATGCTCGACAGGTCGCGCAGGCTGGCGGCCACTTCGTCACGCTGCAACGCCGCCGCGGCATCGGCa                                                                                                     >  1:4383614/1‑150 (MQ=255)
                                                   tctcATCCGTCGCCTGGCGAATGTCCACCACCTGCTGCAACGCCGCCTGCATGGCATGACTTTCAGCCTGAGCCTCGTCAGCCATGCTCGACAGGTCGCGCAGGCTGGCGGCCACTTCGTCACGCTGCAACGCCGCCGCGGCATCGGCAc                                                                                                    <  1:4550101/150‑1 (MQ=255)
                                                     tcATCCGTCGCCTGGCGAATGTCCACCACCTGCTGCAACGCCGCCTGCATGGCATGGCTTTCAGCCTGAGCCTCGTCAGCCATGCTCGACAGGTCGCGCAGGCTGGCGGCCACTTCGTCACGCTGCAACGCCGCCGCGGCATCGGCACCg                                                                                                  <  2:2246165/150‑1 (MQ=255)
                                                     tcATCCGTCGCCTGGCGAATGTCCACCACCTGCTGCAACGCCGCCTGCATGGCATGGCTTTCAGCCTGAGCCTCGTCAGCCATGCTCGACAGGTCGCGCAGGCTGGCGGCCACTTCGTCACGCTGCAACGCCGCCGCGGCATCGGCACCg                                                                                                  >  1:598180/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)
                                                                 tGGCGAATGTCCACCACCTGCTGCAACGCCGCCTGCATGGCATGGCTTTCAGCCTGAGCCTCGTCAGCCAGGCTCGACAGGTCGCGCAGGCTGGCGGCCACTTCGTCACGCTGCAACGCCGCCGCGGCATCGGCACCGGCATTGCGCTgg                                                                                      <  2:1834149/150‑1 (MQ=255)
                                                                  ggCGAATGTCCACCACCTGCTGCAACGCCGCCTGCATGGCATGGCTTTCAGCCTGAGCCTCGTCAGCCATGCTCGACAGGTCGCGCAGGCTGGCGGCCACTTCGTCACGCTGCAACGCCGCCGCGGCATCGGCACCGGCATTGCGCTggg                                                                                     <  1:967415/150‑1 (MQ=255)
                                                                     gAATGTCCACCACCTGCTGCAACGCCGCCTGCATGGCATGGCTTTCAGCCTGAGCCTCGTCAGCCATGCTCGACAGGTCGCGCAGGCTGGCGGCCACTTCGTCACGCTGCAACGCCGCCGCGGCATCGGCACCGGCATTGCGCTGGGCCa                                                                                  >  2:1997588/1‑150 (MQ=255)
                                                                      aaTGTCCACCACCTGCTGCAACGCCGCCTGCATGGCATGGCTTTCAGCCTGAGCCTCGTCAGCCATGCTCGACAGGTCGCGCAGGCTGGCGGCCACTTCGTCACGCTGCAACGCCGCCGCGGCATCGGCACCGGCATTGCGCTGGGCCAt                                                                                 <  2:3742527/150‑1 (MQ=255)
                                                                            caccacCTGCTGCAACGCCGCCTGCATGGCATGGCTTTCAGCCTGAGCCTCGTCAGCCATGCTCGACAGGTCGCGCAGGCTGGCGGCCACTTCGTCACGCTGCAAGGCCGCCGCGGCATCGGCACCGGCATTGCGCTGGGCCATGGCGcc                                                                           <  1:2025803/150‑1 (MQ=255)
                                                                            caccacCTGCTGCAACGCCGCCTGCATGGCATGGCTTTCAGCCTGAGCCTCGTCAGCCATGCTCGACAGGTCGCGCAGGCTGGCGGCCACTTCGTCACGCTGCAACGCCGCCGCGGCATCGGCACCGGCATTGCGCTGGGCCATGGCGcc                                                                           >  1:3827662/1‑150 (MQ=255)
                                                                                  ctgctgCAACGCCGCCTGCATGGCATGGCTTTCAGCCTGAGCCTCGTCAGCCATGCTCGACAGGTCGCGCAGGCTGGCGGCCACTTCGTCACGCTGCAACGCCGCCGCGGCATCGGCACCGGCATTGCGCTGGGCCATGGCGCCGAtctc                                                                     <  2:3071351/150‑1 (MQ=255)
                                                                                   tgctgcAACGCCGCCTGCATGGCATGGCTTTCAGCCTGAGCCTCGTCAGCCATGCTCGACAGGTCGCGCAGGCTGGCGGCCACTTCGTCACGCTGCAACGCCGCCGCGGCATCGGCACCGGCATTGCGCTGGGCCATGGCGCCGATCTCg                                                                    <  2:2723781/150‑1 (MQ=255)
                                                                                    gctgcAACGCCGCCTGCATGGCATGGCTTTCAGCCTGAGCCTCGTCAGCCATGCTCGACAGGTCGCGCAGGCTGGCGGCCACTTCGTCACGCTGCAACGCCGCCGCGGCATCGGCACCGGCATTGCGCTGGGCCATGGCGCCGATCTCGa                                                                   >  1:148729/1‑150 (MQ=255)
                                                                                     ctgcAACGCCGCCTGCATGGCATGGCTTTCAGCCTGAGCCTCGTCAGCCATGCTCGACAGGTCGCGCAGGCTGGCGGCCACTTCGTCACGCTGCAACGCCGCCGCGGCATCGGCACCGGCATTGCGCTGGGCCATGGCGCCGATCTCGAc                                                                  >  2:4121855/1‑150 (MQ=255)
                                                                                           cgccgcCTGCATGGCATGGCTTTCAGCCTGAGCCTCGTCAGCCATGCTCGACAGGTCGCGCAGGCTGGCGGCCACTTCGTCACGCTGCAACGCCGCCGCGGCATCGGCACCGGCATTGCGCTGGGCCATGGCGCCGATCTCGACGCCGGt                                                            >  1:511495/1‑150 (MQ=255)
                                                                                             ccgccTGCATGGCATGGCTTTCAGCCTGAGCCTCGTCAGCCATGCTCGACAGGTCGCGCAGGCTGGCGGCCACTTCGTCACGCTGCAACGCCGCCGCGGCATCGGCACCGGCATTGCGCTGGGCCATGGCGCCGATCTCGACGCCGGTAc                                                          >  2:757220/1‑150 (MQ=255)
                                                                                                 cTGCATGGCATGGCTTTCAGCCTGAGCCTCGTCAGCCATGCTCGACAGGTCGCGCAGGCTGGCGGCCACTTCGTCACGCTGCAACGCCGCCGCGGCATCGGCACCGGCATTGCGCTGGGCCATGGCGCCGATCTCGACGCCGGTACGCTg                                                      >  1:2972151/1‑150 (MQ=255)
                                                                                                  tGCATGGCATGGCTTTCAGCCTGAGCCTCGTCAGCCATGCTCGACAGGTCGCGCAGGCTGGCGGCCACTTCGTc                                                                                                                                 <  1:2235445/74‑1 (MQ=255)
                                                                                                  tGCATGGCATGGCTTTCAGCCTGAGCCTCGTCAGCCATGCTCGACAGGTCGCGCAGGCTGGCGGCCACTTCGTCACGCTGCAACGCCGCCGCGGCATCGGCACCGGCATTGCGCTGGGCCATGGCGCCGATCTCGACGCCGGTACGCTgg                                                     <  1:5260008/150‑1 (MQ=255)
                                                                                                    catggcatggCTTTCAGCCTGAGCCTCGTCAGCCATGCTCGACAGGTCGCGCAGGCTGGCGGCCACTTCGTCACGCTGCAACGCCGCCGCGGCATCGGCACCGGCATTGCGCTGGGCCATGGCGCCGATCTCGACGCCGATACGCTGGGc                                                   >  2:4937269/1‑150 (MQ=255)
                                                                                                     atggcatggcTTTCAGCCTGAGCCTCGTCAGCCATGCTCGACAGGTCGCGCAGGCTGGCGGCCACTTCGTCACGCTGCAACGCCGCCGCGGCATCGGCACCGGCATTGCGCTGGGCCATGGCGCCGATCTCGACGCCGGTACGCTgg                                                     <  2:1348001/147‑1 (MQ=255)
                                                                                                         catggcTTTCAGCCTGAGCCTCGTCAGCCATGCTCGACAGGTCGCGCAGGCTGGCGGCCACTTCGTCACGCTGCAACGCCGCCGCGGCATCGGCACCGGCATTGCGCTGGGCCATGGCGCCGATCTCGACGCCGGTACGCTGGGCCACCt                                              >  2:100281/1‑150 (MQ=255)
                                                                                                          atggcTTTCAGCCTGAGCCTCGTCAGCCATGCTCGACAGGTCGCGCAGGCTGGCGGCCACTTCGTCACGCTGCAACGCCGCCGCGGCATCGGCACCGGCATTGCGCTGGGCCATGGCGCCGATCTCGACGCCGGTACGCTGGGCCACCTc                                             >  2:3201710/1‑150 (MQ=255)
                                                                                                             gcTTTCAGCCTGAGCCTCGTCAGCCATGCTCGACAGGTCGCGCAGGCTGGCGGCCACTTCGTCACGCTGCAACGCCGCCGCGGCATCGGCACCGGCATTGCGCTGGGCCATGGCGCCGATCTCGACGCCGGTACGCTGGGCCACCTCGcc                                          <  2:2886718/150‑1 (MQ=255)
                                                                                                              cTTTCAGCCTGAGCCTCGTCAGCCATGCTCGACAGGTCGCGCAGGCTGGCGGCCACTTCGTCACGCTGCAACGCCGCCGCGGCATCGGCACCGGCATTGCGCTGGGCCATGGCGCCGATCTCGACGCCGGTACGCTGGGCCACCTCGCCg                                         <  2:3488642/150‑1 (MQ=255)
                                                                                                              cTTTCAGCCTGAGCCTCGTCAGCCATGCTCGACAGGTCGCGCAGGCTGGCGGCCACTTCGTCACGCTGCAACGCCGCCGCGGCATCGGCACCGGCATTGCGCTGGGCCATGGCGCCGATCTCGACGCCGGTACGCTGGGCCACCTCGCCg                                         >  1:2106774/1‑150 (MQ=255)
                                                                                                              cTTTCAGCCTGAGCCTCGTCAGCCATGCTCGACAGGTCGCGCAGGCTGGCGGCCACTTCGTCACGCTGCAACGCCGCCGCGGCATCGGCACCGGCATTGCGCTGGGCCATGGCGCCGATCTCGACGCCGGTACGCTGGGCCACCTCGCCg                                         >  2:45915/1‑150 (MQ=255)
                                                                                                               tttCAGCCTGAGCCTCGTCAGCCATGCTCGACAGGTCGCGCAGGCTGGCGGCCACTTCGTCACGCTGCAACGCCGCCGCGGCATCGGCACCGGCAtt                                                                                             >  2:1519960/1‑97 (MQ=255)
                                                                                                               tttCAGCCTGAGCCTCGTCAGCCATGCTCGACAGGTCGCGCAGGCTGGCGGCCACTTCGTCACGCTGCAACGCCGCCGCGGCATCGGCACCGGCAtt                                                                                             <  1:1519960/97‑1 (MQ=255)
                                                                                                               tttCAGCCTGAGCCTCGTCAGCCATGCTCGACAGGTCGCGCAGGCTGGCGGCCACTTCGTCACGCTGCAACGCCGCCGCGGCATCGGCACCGGCATTGCGCTGGGCCATGGCGCCGATCTCGACGCCGGTACGCTGGGCCACCTCGCCgg                                        >  1:1423512/1‑150 (MQ=255)
                                                                                                                 tCAGCCTGAGCCTCGTCAGCCATGCTCGACAGGTCGCGCAGGCTGGCGGCCACTTCGTCACGCTGCAACGCCGCCGCGGCATCGGCACCGGCATTGCGCTGGGCCATGGCGCCGATCTCGACGCCGGTACGCTGGGCCACCTCGCCGGcc                                      >  2:1432591/1‑150 (MQ=255)
                                                                                                                   aGCCTGAGCCTCGTCAGCCATGCTCGACAGGTCGCGCAGGCTGGCGGCCACTTCGTCACGCTGCAACGCCGCCGCGGCATCGGCACCGGCATTGCGCTGGGCCATGGCGCCGATCTCGACGCCGGTACGCTGGGCCACCTCGCCGGCCTc                                    >  1:1673678/1‑150 (MQ=255)
                                                                                                                             tCGTCAGCCATGCTCGACAGGTCGCGCAGGCTGGCGGCCACTTCGTCACGCTGCAACGCCGCCGCGGCATCGGCACCGGCATTGCGCTGGGCCATGGCGCCGATCTCGACGCCGGTACGCTGGGCCACCTCGCCGGCCTCGCGGACGATg                          <  2:1868703/150‑1 (MQ=255)
                                                                                                                             tCGTCAGCCATGCTCGACAGGTCGCGCAGGCTGGCGGCCACGTCGTCACGCTGCAACGCCGCCGCGGCATCGGCACCGGCATTGCGCTGGGCCATGGCGCCGATCTCGACGCCGGTACGCTGGGCCACCTCGCCGGCCTCGCGGACGATg                          <  2:3523595/150‑1 (MQ=255)
                                                                                                                              cGTCAGCCATGCTCGACAGGTCGCGCAGGCTGGCGGCCACTTCGTCACGCTGCAACGCCGCCGCGGCATCGGCACCGGCATTGCGCTGGGCCATGGCGCCGATCTCGACGCCGGTACGCTGGGCCACCTCGCCGGCCTCGCGGACGATgg                         <  2:4261715/150‑1 (MQ=255)
                                                                                                                               gTCAGCCATGCTCGACAGGTCGCGCAGGCTGGCGGCCACTTCGTCACGCTGCAACGCCGCCGCGGCATCGGCACCGGCATTGCGCTGGGCCATGGCGCCGATCTCGACGCCGGTACGCTGGGCCACCTCGCCGGCCTCGCGGACGATggg                        <  2:26514/150‑1 (MQ=255)
                                                                                                                               gTCAGCCATGCTCGACAGGTCGCGCAGGCTGGCGGCCACTTCGTCACGCTGCAACGCCGCCGCGGCATCGGCACCGGCATTGCGCTGGGCCATGGCGCCGATCTCGACGCCGGTACGCTGGGCCACCTCGCCGGCCTCGCGGACGATggg                        >  2:590490/1‑150 (MQ=255)
                                                                                                                                  aGCCATGCTCGACAGGTCGCGCAGGCTGGCGGCCACTTCGTCACGCTGCAACGCCGCCGCGGCATCGGCACCGGCATTGCGCTGGGCCATGGCGCCGAtctc                                                                     >  2:5312634/1‑102 (MQ=255)
                                                                                                                                  aGCCATGCTCGACAGGTCGCGCAGGCTGGCGGCCACTTCGTCACGCTGCAACGCCGCCGCGGCATCGGCACCGGCATTGCGCTGGGCCATGGCGCCGATCTCGACGCCGGTACGCTGGGCCACCTCTCCGGCCTCGCGGACGATGGGCTg                     <  2:1190964/150‑1 (MQ=255)
                                                                                                                                  aGCCATGCTCGACAGGTCGCGCAGGCTGGCGGCCACTTCGTCACGCTGCAACGCCGCCGCGGCATCGGCACCGGCATTGCGCTGGGCCATGGCGCCGATCTCGACGCCGGTACGCTGGGCCACCTCGCCGGCCTCGCGGACGATGGGCTg                     >  2:574360/1‑150 (MQ=255)
                                                                                                                                   gCCATGCTCGACAGGTCGCGCAGGCTGGCGGCCACTTCGTCACGCTGCAACGCCGCCGCGGCATCGGCACCGGCATTGCGCTGGGCCATGGCGCCGATCTCGACGCCGGTACGCTGGGCCACCTCGCCGGCCTCGCGGACGATGGGCTGc                    >  1:3157853/1‑150 (MQ=255)
                                                                                                                                    gcctGCGCGACCTGTCGAGCATGGCTGACGGCCACTTCGTCACGCTGCAACGCCGCCGCGGCATCGGCACCGGCATTGCGCTGGGCCATGGCGCCGATCTCGACGCCGGTACGCTGGGCCACCTCGCCGGCCTCGCGGACGATGGGCTGc                    >  2:4373309/4‑150 (MQ=255)
                                                                                                                                    ccATGCTCGACAGGTCGCGCAGGCTGGCGGCCACTTCGTCACGCTGCAACGCCGCCGCGGCATCGGCACCGGCATTGCGCTGGGCCATGGCGCCGATCTCGACGCCGGTACGCTGGGCCACCTCGCCGGCCTCGCGGACGATGGGCTGCa                   <  2:2679922/150‑1 (MQ=255)
                                                                                                                                     cctGCGCGACCTGTCGAGCAGGCTGGCGGCCACTTCGTCACGCTGCAACGCCGCCGCGGCATCGGCACCGGCATTGCGCTGGGCCATGGCGCCGATCTCGACGCCGGTACGCTGGGCCACCTCGCCGGCCTCGCGGACGATGGGCTGCAg                  >  1:4851020/3‑150 (MQ=255)
                                                                                                                                     cATGCTCGACAGGTCGCGCAGGCTGGCGGCCACTTCGTCACGCTGCAACGCCGCCGCGGCATCGGCACCGGCATTGCGCTGGGCCATgg                                                                               <  1:1881759/89‑1 (MQ=255)
                                                                                                                                      ctGCGCGACCTGTCGAGCATGGCTGGCGGCCACTTCGTCACGCTGCAACGCCGCCGCGGCATCGGCACCGGCATTGCGCTGGGCCATGGCGCCGATCTCGACGCCGGTACGCTGGGCCACCTCGCCGGCCTCGCGGACGATGGGCTGCAg                  >  2:3885519/2‑150 (MQ=255)
                                                                                                                                      aTGCTCGACAGGTCGCGCAGGCTGGCGGCCACTTCGTCACGCTGCAACGCCGCCGCGGCATCGGCACCGGCATTGCGCTGGGCCATGGCGCCGATCTCGACGCCGGTACGCTGGGCCACCTCGCCGGCCTCGCGGACGATGGGCTGCAGc                 >  2:808329/1‑150 (MQ=255)
                                                                                                                                      aTGCTCGACAGGTCGCGCAGGCTGGCGGCCACTTCGTCACGCTGCAACGCCGCCGCGGCATCGGCACCGGCATTGCGCTGGGCCATGGCGCCGATCTCGACGCCGGTACGCTGGGCCACCTCGCCGGCCTCGCGGACGATGGGCTGCAGc                 <  2:2861730/150‑1 (MQ=255)
                                                                                                                                      aTGCTCGACAGGTCGCGCAGGCTGGCGGCCACTTCGTCACGCTGCAACGCCGCCGCGGCATCGGCACCGGCATTGCGCTGGGCCATGGCGCCGATCTCGACGCCGGTACGCTGGGACACCTCGCCGGCCTCGCGGACGATGGGCTGCAGc                 <  1:5020749/150‑1 (MQ=255)
                                                                                                                                       tGCTCGACAGGTCGCGCAGGCTGGCGGCCACTTCGTCACGCTGCAACGCCGCCGCGGCATCGGCACCGGCATTGCGCTGGGCCATGGCGCCGATCTCGACGCCGGTACGCTGGGCCACCTCGCCGGCCTCGCGGACGATggg                        <  2:4786778/142‑1 (MQ=255)
                                                                                                                                       tGCTCGACAGGTCGCGCAGGCTGGCGGCCACTTCGTCACGCTGCAACGCCGCCGCGGCATCGGCACCGGCATTGCGCTGGGCCATGGCGCCGATCTCGACGCCGGTACGCTGGGCCACCTCGCCGGCCTCGCGGACGATggg                        >  1:4786778/1‑142 (MQ=255)
                                                                                                                                         cTCGACAGGTCGCGCAGGCTGGCGGCCACTTCGTCACGCTGCAACGCCGCCGCGGCATCGGCACCGGCATTGCGCTGGGCCATGGCGCCGATCTCGACGCCGGTACGCTGGGCCACCTCGCCGGCCTCGCGGACGATGGGCTGCAGCTTg              <  2:4051957/150‑1 (MQ=255)
                                                                                                                                          tCGACAGGTCGCGCAGGCTGGCGGCCACTTCGTCACGCTGCAACGCCGCCGCGGCATCGGCACCGGCATTGCGCTGGGCCATGGCGCCGATCTCGACGCCGGTACGCTGGGCCACCTCGCCGGCCTCGCGGACGATGGGCTGCAGCTTGt             >  2:3595401/1‑150 (MQ=255)
                                                                                                                                           cGACAGGTCGCGCAGGCTGGCGGCCACTTCGTCACGCTGCAACGCCGCCGCGGCATCGGCACCGGCATTGCGCTGGGCCATGGCGCCGATCTCGACGCCGGTACGCTGGGCCACCTCGCCGGCCTCGCGGACGATGGGCTGCAGCTTGTc            <  1:937702/150‑1 (MQ=255)
                                                                                                                                           cGACAGGTCGCGCAGGCTGGCGGCCACTTCGTCACGCTGCAACGCCGCCGCGGCATCGGCACCGGCATTGCGCTGGGCCATGGCGCCGATCTCGACGCCGGTACGCTGGGCCACCTCGCCGGCCTCGCGGACGATGGGCTGCAGCTTGTc            <  2:4240068/150‑1 (MQ=255)
                                                                                                                                           cGACAGGTCGCGCAGGCTGGCGGCCACTTCGTCACGCTGCAACGCCGCCGCGGCATCGGCACCGGCATTGCGCTGGGCCATGGCGCCGATCTCGACGCCGGTACGCTGGGCCACCTCGCCGGCCTCGCGGACGATGGGCTGCAGCTTGTc            <  1:747065/150‑1 (MQ=255)
                                                                                                                                           cGACAGGTCGCGCAGGCTGGCGGCCACTTCGTCACGCTGCAACGCCGCCGCGGCATCGGCACCGGCATTGCGCTGGGCCATGGCGCCGATCTCGACGCCGGTACGCTGGGCCACCTCGCCGGCCTCGCGGACGATGGGCTGCAGCTTGTc            <  2:1493236/150‑1 (MQ=255)
                                                                                                                                           cGACAGGTCGCGCAGGCTGGCGGCCACTTCGTCACGCTGCAACGCCGCCGCGGCATCGGCACCGGCATTGCGCTGGGCCATGGCGCCGATCTCGACGCCGGTACGCTGGGCCACCTCGCCGGCCTCGCGGACGATGGGCTGCAGCTTGTc            <  1:5054146/150‑1 (MQ=255)
                                                                                                                                           cGACAGGTCGCGCAGGCTGGCGGCCACTTCGTCACGCTGCAACGCCGCCGCGGCATCGGCACCGGCATTGCGCTGGGCCATGGCGCCGATCTCGACGCCGGTACGCTGGGCCACCTCGCCGGCCTCGCGGACGATGGGCTGCAGCTTGTc            <  2:539700/150‑1 (MQ=255)
                                                                                                                                            gACAGGTCGCGCAGGCTGGCGGCCACTTCGTCACGCTGCAACGCCGCCGCGGCATCGGCACCGGCATTGCGCTGGGCCATGGCGCCGATCTCGACGCCGGTACGCTGGGCCACCTCGCCGGCCTCGCGGACGATGGGCTGCAGCTTGTCg           >  2:1039437/1‑150 (MQ=255)
                                                                                                                                            gACAGGTCGCGCAGGCTGGCGGCCACTTCGTCACGCTGCAACGCCGCCGCGGCATCGGCACCGGCATTGCGCTGGGCCATGGCGCCGATCTCGACGCCGGTACGCTGGGCCACCTCGCCGGCCTCGCGGACGATGGGCTGCAGCTTGTCg           >  1:3416403/1‑150 (MQ=255)
                                                                                                                                              cAGGTCGCGCAGGCTGGCGGCCACTTCGTCACGCTGCAACGCCGCCGCGGCATCGGCACCGGCATTGCGCTGGGCCATGGCGCCGAACTCGACGCCGGTACGCTGGGCCACCTCGCCGGCCTCGCGGACGATGGGCTGCAGCTTGTCCAc         <  2:1674959/150‑1 (MQ=255)
                                                                                                                                               aGGTCGCGCAGGCTGGCGGCCACTTCGTCACGCTGCAACGCCGCCGCGGCATCGGCACCGGCATTGCGCTGGGCCATGGCGCCGATCTCGACGCCGGTACGCTGGGCCACCTCGCCGGCCTCGCGGACGATGGGCTGCAGCTTGTCGACa        >  2:297546/1‑150 (MQ=255)
                                                                                                                                                 gTCGCGCAGGCTGGCGGCCACTTCGTCACGCTGCAACGCCGCCGCGGCATCGGCACCGGCATTGCGCTGGGCCATGGCGCCGATCTCGACGCCGGTACGCTGGGCCACCTCGCCGGCCTCGCGGACGATGGGCTGCAGCTTGTCGACaaa      <  1:1039437/150‑1 (MQ=255)
                                                                                                                                                 gTCGCGCAGGCTGGCGGCCACTTCGTCACGCTGCAACGCCGCCGCGGCATCGGCACCGGCATTGCGCTGGGCCATGGCGCCGATCTCGACGCCGGTACGCTGGGCCACCTCACCGGCCTCGCGGACGATGGGCTGCAGCTTGTCGACaaa      >  1:1300413/1‑150 (MQ=255)
                                                                                                                                                    gcgcAGGCTGGCGGCCACTTCGTCACGCTGCAACGCCGCCGCGGCATCGGCACCGGCATTGCGCTGGGCCATGGCGCCGATCTCGACGCCGGTACGCTGGGCCACCTCGCCGGCCTCGCGGACGATGGGCTGCAGCTTGTCGACAAAGCg   >  1:2992737/1‑150 (MQ=255)
                                                                                                                                                     cgcAGGCTGGCGGCCACTTCGTCACGCTGCAACGCCGCCGCGGCATCGGCACCGGCATTGCGCTGGGCCATGGCGCCGATCTCGACGCCGGTACGCTGGGCCACCTCGCCGGCCTCGCGGACGATGGGCTGCAGCTTGTCGACAAAGCGa  <  1:1343862/150‑1 (MQ=255)
                                                                                                                                                     cgcAGGCTGGCGGCCACTTCGTCACGCTGCAACGCCGCCGCGGCATCGGCACCGGCATTGCGCTGGGCCATGGCGCCGATCTCGACGCCGGTACGCTGGGCCACCTCGCCGGCCTCGCGGACGATGGGCTGCAGCTTGTCGACAAAGCGa  <  1:5146308/150‑1 (MQ=255)
                                                                                                                                                     |                                                                                                                                                      
GCCCGGCAAGGTTCTCAATCAGGCTCGCCAGTTGCGTCGAGGAACGGCTGTTCTCATCCGTCGCCTGGCGAATGTCCACCACCTGCTGCAACGCCGCCTGCATGGCATGGCTTTCAGCCTGAGCCTCGTCAGCCATGCTCGACAGGTCGCGCAGGCTGGCGGCCACTTCGTCACGCTGCAACGCCGCCGCGGCATCGGCACCGGCATTGCGCTGGGCCATGGCGCCGATCTCGACGCCGGTACGCTGGGCCACCTCGCCGGCCTCGCGGACGATGGGCTGCAGCTTGTCGACAAAGCGA  >  NC_002947/5558237‑5558535

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: