Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 5,558,380 A→C 5.8% L408R (CTG→CGG) ‡ 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,3800AC 5.8% 706.4 / 23.0 207L408R (CTG→CGG) ‡PP_4888methyl‑accepting chemotaxis transducer
Reads supporting (aligned to +/- strand):  ref base A (103/92);  new base C (7/5);  total (110/97)
Fisher's exact test for biased strand distribution p-value = 7.73e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.42e-01

Only 100 of 213 total aligned reads displayed.

CAACCTGCCCGGCAAGGTTCTCAATCAGGCTCGCCAGTTGCGTCGAGGAACGGCTGTTCTCATCCGTCGCCTGGCGAATGTCCACCACCTGCTGCAACGCCGCCTGCATGGCATGGCTTTCAGCCTGAGCCTCGTCAGCCATGCTCGACAGGTCGCGCAGGCTGGCGGCCACTTCGTCACGCTGCAACGCCGCCGCGGCATCGGCACCGGCATTGCGCTGGGCCATGGCGCCGATCTCGACGCCGGTACGCTGGGCCACCTCGCCGGCCTCGCGGACGATGGGCTGCAGCTTGTCGACA  >  NC_002947/5558231‑5558529
                                                                                                                                                     |                                                                                                                                                      
cAACCTGCCCGGCAAGGTTCTCAATCAGGCTCGCCAGTTGCGTCGAGGAACGGCTGTTCTCATCCGTCGCCTGGCGAATGTCCACCACCTGCTGCAACGCCGCCTGCATGGCATGGCTTTCAGCCTGAGCCTCGTCAGCCATGCTCGACa                                                                                                                                                        <  1:4576006/150‑1 (MQ=255)
 aaCCTGCCCGGCAAGGTTCTCAATCAGGCTCGCCAGTTGCGTCGAGGAACGGCTGTTCTCATCCGTCGCCTGGCGAATGTCCACCACCTGCTGCAACGCCGCCTGCATGGCATGGCTTTCAGCCTGAGCCTCGTCAGCCATGCTCGACAg                                                                                                                                                       >  1:20525/1‑150 (MQ=255)
  aCCTGCCCGGCAAGGTTCTCAATCAGGCTCGCCAGTTGCGTCGAGGAACGGCTGTTCTCATCCGTCGCCTGGCGAATGTCCACCACCTGCTGCAACGCCGCCTGCATGGCATGGCTTTCAGCCTGAGCCTCGTCAGCCATGCTCGACAgg                                                                                                                                                      >  2:3825521/1‑150 (MQ=255)
    cTGCCCGGCAAGGTTCTCAATCAGGCTCGCCAGTTGCGTCGAGGAACGGCTGTTCTCATCCGTCGCCTGGCGAATGTCCACCACCTGCTGCAACGCCGCCTGCATGGCATGGCTTTCAGCCTGAGCCTCGTCAGCCATGCTCGACAGGTc                                                                                                                                                    >  2:4490809/1‑150 (MQ=255)
      gCCCGGCAAGGTTCTCAATCAGGCTCGCCAGTTGCGTCGAGGAACGGCTGTTCTCATCCGTCGCCTGGCGAATGTCCACCACCTGCTGCAACGCCGCCTGCATGGCATGGCTTTCAGCCTGAGCCTCGTCAGCCATGCTCGACAGGTcgc                                                                                                                                                  >  1:435189/1‑150 (MQ=255)
      gCCCGGCAAGGTTCTCAATCAGGCTCGCCAGTTGCGTCGAGGAACGGCTGTTCTCATCCGTCGCCTGGCGAATGTCCACCACCTGCTGCAACGCCGCCTGCATGGCATGGCTTTCAGCCTGAGCCTCGTCAGCCATGCTCGACAGGTcgc                                                                                                                                                  <  1:2062947/150‑1 (MQ=255)
        ccGGCAAGGTTCTCAATCAGGCTCGCCAGTTGCGTCGAGGAACGGCTGTTCTCATCCGTCGCCTGGCGAATGTCCACCACCTGCTGCAACGCCGCCTGCATGGCATGGCTTTCAGCCTGAGCCTCGTCAGCCATGCTCGACAGGTCgcgc                                                                                                                                                <  1:3148683/150‑1 (MQ=255)
        ccGGCAAGGTTCTCAATCAGGCTCGCCAGTTGCGTCGAGGAACGGCTGTTCTCATCCGTCGCCTGGCGAATGTCCACCACCTGCTGCAACGCCGCCTGCATGGCATGGCTTTCAGCCTGAGCCTCGTCAGCCATGCTCGACAGGTCgcgc                                                                                                                                                <  1:94175/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)
                   ctcAATCAGGCTCGCCAGTTGCGTCGAGGAACGGCTGTTCTCATCCGTCGCCTGGCGAATGTCCACCACCTGCTGCAACGCCGCCTGCATGGCATGGCTTTCAGCCTGAGCCTCGTCAGCCATGCTCGACAGGTCGCGCAGGCTggcggc                                                                                                                                    <  2:1259967/150‑1 (MQ=255)
                    tcAATCAGGCTCGCCAGTTGCGTCGAGGAACGGCTGTTCTCATCCGTCGCCTGGCGAATGTCCACCACCTGCTGCAACGCCGCCTGCATGGCATGGCTTTCAGCCTGAGCCTCGTCAGCCATGCGCGACCTGTCGAGCatg                                                                                                                                             <  1:404339/141‑3 (MQ=255)
                      aaTCAGGCTCGCCAGTTGCGTCGAGGAACGGCTGTTCTCATCCGTCGCCTGGCGAATGTCCACCACCTGCTGCAACGCCGCCTGCATGGCATGGCTTTCAGCCTGAGCCTCGTCAGCCATGCTCGACAGGTCGCGCAGGCTGGCGGCCAc                                                                                                                                 >  1:2680411/1‑150 (MQ=255)
                      aaTCAGGCTCGCCAGTTGCGTCGAGGAACGGCTGTTCTCATCCGTCGCCTGGCGAATGTCCACCACCTGCTGCAACGCCGCCTGCATGGCATGGCTTTCAGCCTGAGCCTCGTCAGCCATGCTCGACAGGTCGCGCAGGCTGGCGGCCAc                                                                                                                                 >  2:2309394/1‑150 (MQ=255)
                       aTCAGGCTCGCCAGTTGCGTCGAGGAACGGCTGTTCTCATCCGTCGCCTGGCGAATGTCCACCACCTGCTGCAACGCCGCCTGCATGGCATGGCTTTCAGCCTGAGCCTCGTCAGCCATGCTCGACAGGTCGCGCAGGCTGGCGGCCACt                                                                                                                                >  1:5083854/1‑150 (MQ=255)
                        tCAGGCTCGCCAGTTGCGTCGAGGAACGGCTGTTCTCATCCGTCGCCAGGCGAATGTCCACCACCTGCTGCAACGCCGCCTGCATGGCATGGCTTTCAGCCTGAGCCTCGTCAGCCATGCTCGACAGGTCGCGCAGGCTGGCGGCCACtt                                                                                                                               >  1:1913998/1‑150 (MQ=255)
                          aGGCTCGCCAGTTGCGTCGAGGAACGGCTGTTCTCATCCGTCGCCTGGCGAATGTCCACCACCTGCTGCAACGCCGCCTGCATGGCATGGCTTTCAGCCTGAGCCTCGTCAGCCATGCTCGACAGGTCGCGCAGGCTGGCGGCCACTTCg                                                                                                                             <  2:4221477/150‑1 (MQ=255)
                            gCTCGCCAGTTGCGTCGAGGAACGGCTGTTCTCATCCGTCGCCTGGCGAATGTCCACCACCTGCTGCAACGCCGCCTGCATGGCATGGCTTTCAGCCTGAGCCTCGTCAGCCATGCTCGACAGGTCGCGCAGGCTGGCGGCCACTTCGTc                                                                                                                           >  1:3574912/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)
                                      tGCGTCGAGGAACGGCTGTTCTCATCCGTCGCCTGGCGAATGTCCACCACCTGCTGCAACGCCGCCTGCATGGCATGGCTTTCAGCCTGAGCCTCGTCAGCCATGCTCGACAGGTCGCGCAGGCTGGCGGCCACTTCGTCACGCTGCAAc                                                                                                                 <  2:4048953/150‑1 (MQ=255)
                                        cGTCGAGGAACGGCTGTTCTCATCCGTCGCCTGGCGAATGTCCACCACCTGCTGCAACGCCGCCTGCATGGCATGGCTTTCAGCCTGAGCCTCGTCAGCCATGCTCGACAGGTCGCGCAGGCTGGCGGCCACTTCGTCACGCTGCAAcgc                                                                                                               <  1:2320971/150‑1 (MQ=255)
                                              ggAACGGCTGTTCTCATCCGTCGCCTGGCGAATGTCCACCACCTGCTGCAACGCCGCCTGCATGGCATGGCTTTCAGCCTGAGCCTCGTCAGCCATGCTCGACAGGTCGCGCAGGCTGGCGGCCACTTCGTCACGCTGCAAcgccgccgc                                                                                                         >  1:910352/1‑150 (MQ=255)
                                              ggAACGGCTGTTCTCATCCGTCGCCTGGCGAATGTCCACCACCTGCTGCAACGCCGCCTGCATGGCATGGCTTTCAGCCTGAGCCTCGTCAGCCATGCTCGACAGGTCGCGCAGGCTGGCGGCCACTTCGTCACGCTGCAAcgccgccgc                                                                                                         <  2:3022688/150‑1 (MQ=255)
                                               gAACGGCTGTTCTCATCCGTCGCCTGGCGAATGTCCACCACCTGCTGCAACGCCGCCTGCATGGCATGGCTTTCAGCCTGAGCCTCGTCAGCCATGCTCGACAGGTCGCGCAGGCTGGCGGCCACGTCGTCACGCTGCAACGCCGCcgcg                                                                                                        >  1:1148624/1‑150 (MQ=255)
                                                aaCGGCTGTTCTCATCCGTCGCCTGGCGAATGTCCACCACCTGCTGCAACGCCGCCTGCATGGCATGGCTTTCAGCCTGAGCCTCGTCAGCCATGCTCGACAGGTCGCGCAGGCTGGCGGCCACTTCGTCACGCTGCAACGCCGCCGCgg                                                                                                       <  2:2680411/150‑1 (MQ=255)
                                                   ggCTGTTCTCATCCGTCGCCTGGCGAATGTCCACCACCTGCTGCAACGCCGCCTGCATGGCATGGCTTTCAGCCTGAGCCTCGTCAGCCATGCTCGACAGGTCGCGCAGGCTGGCGGCCACTTCGTCACGCTGCAACGCCGCCGCGGCAt                                                                                                    >  2:3068468/1‑150 (MQ=255)
                                                   ggCTGTTCTCATCCGTCGCCTGGCGAATGTCCACCACCTGCTGCAACGCCGCCTGCATGGCATGGCTTTCAGCCTGAGCCTCGTCAGCCATGCTCGACAGGTCGCGCAGGCTGGCGGCCACTTCGTCACGCTGCAACGCCGCCGCGGCAt                                                                                                    <  2:529483/150‑1 (MQ=255)
                                                   ggCTGTTCTCATCCGTCGCCTGGCGAATGTCCACCACCTGCTGCAACGCCGCCTGCATGGCATGGCTTTCAGCCTGAGCCTCGTCAGCCATGCTCGACAGGTCGCGCAGGCTGGCGGCCACTTCGTCACGCTGCAACGCCGCCGCGGCAt                                                                                                    >  2:5260008/1‑150 (MQ=255)
                                                   ggCTGTTCTCATCCGTCGCCTGGCGAATGTCCACCACCTGCTGCAACGCCGCCTGCATGGCATGGCTTTCAGCCTGAGCCTCGTCAGCCATGCTCGACAGGTCGCGCAGGCTGGCGGCCACTTCGTCACGCTGCAACGCCGCCGCGGCAt                                                                                                    <  2:5078931/150‑1 (MQ=255)
                                                      tGTTCTCATCCGTCGCCTGGCGAATGTCCACCACCTGCTGCAACGCCGCCTGCATGGCATGGCTTTCAGCCTGAGCCTCGTCAGCCATGCTCGACAGGTCGCGCAGGCTGGCGGCCACTTCGTCACGCTGCAACGCCGCCGCGGCATCgg                                                                                                 <  1:396226/150‑1 (MQ=255)
                                                          ctcATCCGTCGCCTGGCGAATGTCCACCACCTGCTGCAACGCCGCCTGCATGGCATGGCTTTCAGCCTGAGCCTCGTCAGCCATGCTCGACAGGTCGCGCAGGCTGGCGGCCACTTCGTCACGCTGCAACGCCGCCGCGGCATCGGCAcc                                                                                             <  2:935859/150‑1 (MQ=255)
                                                           tcATCCGTCGCCTGGCGAATGTCCACCACCTGCTGCAACGCCGCCTGCATGGCATGGCTTTCAGCCTGAGCCTCGTCAGCCATGCTCGACAGGTCGCGCAGGCTGGCGGCCACTTCGTCACGCTGCAACGCCGCCGCGGCATCGGCACCg                                                                                            >  1:1496077/1‑150 (MQ=255)
                                                           tcATCCGTCGCCTGGCGAATGTCCACCACCTGCTGCAACGCCGCCTGCATGGCATGGCTTTCAGCCTGAGCCTCGTCAGCCATGCTCGACAGGTCGCGCAGGCTGGCGGCCACTTCGTCACGCTGCAACGCCGCCGCGGCATCGGCACCg                                                                                            <  2:2246165/150‑1 (MQ=255)
                                                           tcATCCGTCGCCTGGCGAATGTCCACCACCTGCTGCAACGCCGCCTGCATGGCATGGCTTTCAGCCTGAGCCTCGTCAGCCATGCTCGACAGGTCGCGCAGGCTGGCGGCCACTTCGTCACGCTGCAACGCCGCCGCGGCATCGGCACCg                                                                                            >  1:598180/1‑150 (MQ=255)
                                                                      cTGGCGAATGTCCACCACCTGCTGCAACGCCGCCTGCATGGCATGGCTTTCAGCCTGAGCCTCGTCAGCCATGCTCGACAGGTCGCGCAGGCTGGCGGCCACTTCGTCACGCTGCAACGCCGCCGCGGCATCGGCACCGGCATTGCGCTg                                                                                 >  1:3730691/1‑150 (MQ=255)
                                                                       tGGCGAATGTCCACCACCTGCTGCAACGCCGCCTGCATGGCATGGCTTTCAGCCTGAGCCTCGTCAGCCATGCTCGACAGGTCGCGCAGGCTGGCGGCCACTTCGTCACGCTGCAACGCCGCCGCGGCATCGGCACCGGCATTGCGCTgg                                                                                <  1:2588449/150‑1 (MQ=255)
                                                                       tGGCGAATGTCCACCACCTGCTGCAACGCCGCCTGCATGGCATGGCTTTCAGCCTGAGCCTCGTCAGCCAGGCTCGACAGGTCGCGCAGGCTGGCGGCCACTTCGTCACGCTGCAACGCCGCCGCGGCATCGGCACCGGCATTGCGCTgg                                                                                <  2:1834149/150‑1 (MQ=255)
                                                                           gAATGTCCACCACCTGCTGCAACGCCGCCTGCATGGCATGGCTTTCAGCCTGAGCCTCGTCAGCCATGCTCGACAGGTCGCGCAGGCTGGCGGCCACTTCGTCACGCTGCAACGCCGCCGCGGCATCTGCACCGGCATTGCGCTGGGCCa                                                                            >  2:787333/1‑150 (MQ=255)
                                                                           gAATGTCCACCACCTGCTGCAACGCCGCCTGCATGGCATGGCTTTCAGCCTGAGCCTCGTCAGCCATGCTCGACAGGTCGCGCAGGCTGGCGGCCACTTCGTCACGCTGCAACGCCGCCGCGGCATCGGCACCGGCATTGCGCTGGGCCa                                                                            >  2:1997588/1‑150 (MQ=255)
                                                                           gAATGTCCACCACCTGCTGCAACGCCGCCTGCATGGCATGGCTTTCAGCCTGAGCCTCGTCAGCCATGCTCGACAGGTCGCGCAGGCTGGCGGCCACTTCGTCACGCTGCAACGCCGCCGCGGCATCGGCACCGGCATTGCGCTGGGCCa                                                                            >  2:1018740/1‑150 (MQ=255)
                                                                            aaTGTCCACCACCTGCTGCAACGCCGCCTGCATGGCATGGCTTTCAGCCTGAGCCTCGTCAGCCATGCTCGACAGGTCGCGCAGGCTGGCGGCCACTTCGTCACGCTGCAACGCCGCCGCGGCATCGGCACCGGCATTGCGCTGGGCCAt                                                                           <  2:3742527/150‑1 (MQ=255)
                                                                              tGTCCACCACCTGCTGCAACGCCGCCTGCATGGCATGGCTTTCAGCCTGAGCCTCGTCAGCCATGCTCGACAGGTCGCGCAGGCTGGCGGCCACTTCGTCACGCTGCAACGCCGCCGCGGCATCGGCACCGGCATTGCGCTGGGCCATgg                                                                         <  1:1475827/150‑1 (MQ=255)
                                                                              tGTCCACCACCTGCTGCAACGCCGCCTGCATGGCATGGCTTTCAGCCTGAGCCTCGTCAGCCATGCTCGACAGGTCGCGCAGGCTGGCGGCCACTTCGTCACGCTGCAACGCCGCCGCGGCATCGGCACCGGCATTGCGCTGGGCCATgg                                                                         >  2:669522/1‑150 (MQ=255)
                                                                                 ccaccaCCTGCTGCAACGCCGCCTGCATGGCATGGCTTTCAGCCTGAGCCTCGTCAGCCATGCTCGACAGGTCGCGCAGGCTGGCGGCCACTTCGTCACGCTGCAACGCCGCCGCGGCATCGGCACCGGCATTGCGCTGGGCCATGgcgc                                                                      >  1:1868703/1‑150 (MQ=255)
                                                                                  caccacCTGCTGCAACGCCGCCTGCATGGCATGGCTTTCAGCCTGAGCCTCGTCAGCCATGCTCGACAGGTCGCGCAGGCTGGCGGCCACTTCGTCACGCTGCAACGCCGCCGCGGCATCGGCACCGGCATTGCGCTGGGCCATGGCGcc                                                                     >  1:3827662/1‑150 (MQ=255)
                                                                                  caccacCTGCTGCAACGCCGCCTGCATGGCATGGCTTTCAGCCTGAGCCTCGTCAGCCATGCTCGACAGGTCGCGCAGGCTGGCGGCCACTTCGTCACGCTGCAACGCCGCCGCGGCATCGGCACCGGCATTGCGCTGGGCCATGGCGcc                                                                     >  1:5125054/1‑150 (MQ=255)
                                                                                  caccacCTGCTGCAACGCCGCCTGCATGGCATGGCTTTCAGCCTGAGCCTCGTCAGCCATGCTCGACAGGTCGCGCAGGCTGGCGGCCACTTCGTCACGCTGCAACGCCGCCGCGGCATCGGCACCGGCATTGCGCTGGGCCATGGCGcc                                                                     >  2:5382150/1‑150 (MQ=255)
                                                                                  caccacCTGCTGCAACGCCGCCTGCATGGCATGGCTTTCAGCCTGAGCCTCGTCAGCCATGCTCGACAGGTCGCGCAGGCTGGCGGCCACTTCGTCACGCTGCAACGCCGCCGCGGCATCGGCACCGGCATTGCGCTGGGCCATGGCGcc                                                                     >  1:4895326/1‑150 (MQ=255)
                                                                                    ttttcagctgcAACGCCGCCTGCATGGCATGGCTTTCAGCCTGAGCCTCGTCAGCCATGCTCGACAGGTCGCGCAGGCTGGCGGCCACTTCGTCACGCTGCAACGCCGCCGCGGCATCGGCACCGGCATTGCGCTGGGCCATGGCGCCGa                                                                   >  2:584466/7‑150 (MQ=255)
                                                                                         tgctgcAACGCCGCCTGCATGGCATGGCTTTCAGCCTGAGCCTCGTCAGCCATGCTCGACAGGTCGCGCAGGCTGGCGGCCACTTCGTCACGCTGCAACGCCGCCGCGGCATCGGCACCGGCATTGCGCTGGGCCATGGCGCCGATCTCg                                                              >  1:3074925/1‑150 (MQ=255)
                                                                                         tgctgcAACGCCGCCTGCATGGCATGGCTTTCAGCCTGAGCCTCGTCAGCCATGCTCGACAGGTCGCGCAGGCTGGCGGCCACCTCGTCACGCTGCAACGCCGCCGCGGCATCGGCACCGGCATTGCGCTGGGCCATGGCGCCGATCTCg                                                              <  2:4255402/150‑1 (MQ=255)
                                                                                          gctgcAACGCCGCCTGCATGGCATGGCTTTCAGCCTGAGCCTCGTCAGCCATGCTCGACAGGTCGCGCAGGCTGGCGGCCACTTCGTCACGCTGCAACGCCGCCGCGGCATCGGCACCGGCATTGCGCTGGGCCATGGCGCCGATCTCGa                                                             <  2:1496077/150‑1 (MQ=255)
                                                                                           ctgcAACGCCGCCTGCATGGCATGGCTTTCAGCCTGAGCCTCGTCAGCCATGCTCGACa                                                                                                                                                        <  1:2876776/59‑1 (MQ=255)
                                                                                           ctgcAACGCCGCCTGCATGGCATGGCTTTCAGCCTGAGCCTCGTCAGCCATGCTCGACa                                                                                                                                                        >  2:2876776/1‑59 (MQ=255)
                                                                                           ctgcAACGCCGCCTGCATGGCATGGCTTTCAGCCTGAGCCTCGTCAGCCATGCTCGACAGGTCGCGCAGGCTGGCGGCCACTTCGTCACGCTGCAACGCCGCCGCGGCATCGGCACCGGCATTGCGCTGGGCCATGGCGCCGATCTCGAc                                                            >  2:4121891/1‑150 (MQ=255)
                                                                                                aCGCCGCCTGCATGGCATGGCTTTCAGCCTGAGCCTCGTCAGCCATGCTCGACAGGTCGCGCAGGCTGGCGGCCACTTCGTCACGCTGCAACGCCGCCGCGGCATCGGCACCGGCATTGCGCTGGGCCATGGCGCCGATCTCGACGCCgg                                                       >  2:157869/1‑150 (MQ=255)
                                                                                                       cTGCATGGCATGGCTTTCAGCCTGAGCCTCGTCAGCCATGCTCGACAGGTCGCGCAGGCTGGCGGCCACTTCGTCACGCTGCAACGCCGCCGCGGCATCGGCACCGGCATTGCGCTGGGCCATGGCGCCGATCTCGACGCCGGTACGCTg                                                >  2:2663743/1‑150 (MQ=255)
                                                                                                        tGCATGGCATGGCTTTCAGCCTGAGCCTCGTCAGCCATGCTCGACAGGTCGCGCAGGCTGGCGGCCACTTCGTCACGCTGCAACGCCGCCGCGGCATCGGCACCGGCATTGCGCTGGGCCATGGCGCCGATCTCGACGCCGGTACGCTgg                                               <  2:144729/150‑1 (MQ=255)
                                                                                                        tGCATGGCATGGCTTTCAGCCTGAGCCTCGTCAGCCATGCTCGACAGGTCGCGCAGGCTGGCGGCCACTTCGTCACGCTGCAACGCCGCCGCGGCATCGGCACCGGCATTGCGCTGGGCCATGGCGCCGATCTCGACGCCGGTACGCTgg                                               >  2:4578700/1‑150 (MQ=255)
                                                                                                        tGCATGGCATGGCTTTCAGCCTGAGCCTCGTCAGCCATGCTCGACAGGTCGCGCAGGCTGGCGGCCACTTCGTCACGCTGCAACGCCGCCGCGGCATCGGCACCGGCATTGCGCTGGGCCATGGCGCCGATCTCGACGCCGGTACGCTgg                                               <  1:5260008/150‑1 (MQ=255)
                                                                                                          catggcatggCTTTCAGCCTGAGCCTCGTCAGCCATGCTCGACAGGTCGCGCAGGCTGGCGGCCACTTCGTCACGCTGCAACGCCGCCGCGGCATCGGCACCGGCATTGCGCTGGGCCATGGCGCCGATCTCGACGCCGATACGCTGGGc                                             >  2:4937269/1‑150 (MQ=255)
                                                                                                           atggcatggcTTTCAGCCTGAGCCTCGTCAGCCATGCTCGACAGGTCGCGCAGGCTGGCGGCCACTTCGTCACGCTGCAACGCCGCCGCGGCATCGGCACCGGCATTGCGCTGGGCCATGGCGCCGATCTCGACGCCGGTACGCTgg                                               >  1:1348001/1‑147 (MQ=255)
                                                                                                           atggcatggcTTTCAGCCTGAGCCTCGTCAGCCATGCTCGACAGGTCGCGCAGGCTGGCGGCCACTTCGTCACGCTGCAACGCCGCCGCGGCATCGGCACCGGCATTGCGCTGGGCCATGGCGCCGATCTCGACGCCGGTACGCTgg                                               <  2:1348001/147‑1 (MQ=255)
                                                                                                                atggcTTTCAGCCTGAGCCTCGTCAGCCATGCTCGACAGGTCGCGCAGGCTGGCGGCCACTTCGTCACGCTGCAACGCCGCCGCGGCATCGGCACCGGCATTGCGCTGGGCCATGGCGCCGATCTCGACGCCGGTACGCTGGGCCACCTc                                       <  1:5382150/150‑1 (MQ=255)
                                                                                                                atggcTTTCAGCCTGAGCCTCGTCAGCCATGCTCGACAGGTCGCGCAGGCTGGCGGCCACTTCGTCACGCTGCAACGCCGCCGCGGCATCGGCACCGGCATTGCGCTGGGCCATGGCGCCGATCTCGACGCCGGTACGCTGGGCCACCTc                                       >  2:3201710/1‑150 (MQ=255)
                                                                                                                   gcTTTCAGCCTGAGCCTCGTCAGCCATGCTCGACAGGTCGCGCAGGCTGGCGGCCACTTCGTCACGCTGCAACGCCGCCGCGGCATCGGCACCGGCATTGCGCTGGGCCATGGCGCCGATCTCGACGCCGGTACGCTGGGCCACCTCGcc                                    <  2:2886718/150‑1 (MQ=255)
                                                                                                                    cTTTCAGCCTGAGCCTCGTCAGCCATGCTCGACAGGTCGCGCAGGCTGGCGGCCACTTCGTCACGCTGCAACGCCGCCGCGGCATCGGCACCGGCATTGCGCTGGGCCATGGCGCCGATCTCGACGCCGGTACGCTGGGCCACCTCGCCg                                   >  2:2773298/1‑150 (MQ=255)
                                                                                                                    cTTTCAGCCTGAGCCTCGTCAGCCATGCTCGACAGGTCGCGCAGGCTGGCGGCCACTTCGTCACGCTGCAACGCCGCCGCGGCATCGGCACCGGCATTGCGCTGGGCCATGGCGCCGATCTCGACGCCGGTACGCTGGGCCACCTCGCCg                                   >  1:2106774/1‑150 (MQ=255)
                                                                                                                    cTTTCAGCCTGAGCCTCGTCAGCCATGCTCGACAGGTCGCGCAGGCTGGCGGCCACTTCGTCACGCTGCAACGCCGCCGCGGCATCGGCACCGGCATTGCGCTGGGCCATGGCGCCGATCTCGACGCCGGTACGCTGGGCCACCTCGCCg                                   >  2:4107571/1‑150 (MQ=255)
                                                                                                                    cTTTCAGCCTGAGCCTCGTCAGCCATGCTCGACAGGTCGCGCAGGCTGGCGGCCACTTCGTCACGCTGCAACGCCGCCGCGGCATCGGCACCGGCATTGCGCTGGGCCATGGCGCCGATCTCGACGCCGGTACGCTGGGCCACCTCGCCg                                   <  2:3488642/150‑1 (MQ=255)
                                                                                                                     tttCAGCCTGAGCCTCGTCAGCCATGCTCGACAGGTCGCGCAGGCTGGCGGCCACTTCGTCACGCTGCAACGCCGCCGCGGCATCGGCACCGGCATTGCGCTGGGCCATGGCGCCGATCTCGACGCCGGTACGCTGGGCCACCTCGCCgg                                  <  1:4387326/150‑1 (MQ=255)
                                                                                                                       tCAGCCTGAGCCTCGTCAGCCATGCTCGACAGGTCGCGCAGGCTGGCGGCCACTTCGTCACGCTGCAACGCCGCCGCGGCATCGGCACCGGCATTGCGCTGGGCCATGGCGCCGATCTCGACGCCGGTACGCTGGGCCACCTCGCCGGcc                                >  2:1432591/1‑150 (MQ=255)
                                                                                                                               aGCCTCGTCAGCCATGCTCGACAGGTCGCGCAGGCTGGCGGCCACTTCGTCACGCTGCAACGCCGCCGCGGCATCGGCACCGGCATTGCGCTGGGCCATGGCGCCGATCTCGACGCCGGTACGCTGGGCCACCTCGCCGGCCTCGCGGAc                        <  2:2096203/150‑1 (MQ=255)
                                                                                                                                 ccTCGTCAGCCATGCTCGACAGGTCGCGCAGGCTGGCGGCCACTTCGTCACGCTGCAACGCCGCCGCGGCATCGGCACCGGCATTGCGCTGGGCCATGGCGCCGATCTCGACGCCGGTACGCTGGGCCACCTCGCCGGCCTCGCGGACGa                      >  1:5351371/1‑150 (MQ=255)
                                                                                                                                   tCGTCAGCCATGCTCGACAGGTCGCGCAGGCTGGCGGCCACTTCGTCACGCTGCAACGCCGCCGCGGCATCGGCACCGGCATTGCGCTGGGCCATGGCGCCGATCTCGACGCCGGTACGCTGGGCCACCTCGCCGGCCTCGCGGACGATg                    <  1:1576315/150‑1 (MQ=255)
                                                                                                                                   tCGTCAGCCATGCTCGACAGGTCGCGCAGGCTGGCGGCCACTTCGTCACGCTGCAACGCCGCCGCGGCATCGGCACCGGCATTGCGCTGGGCCATGGCGCCGATCTCGACGCCGGTACGCTGGGCCACCTCGCCGGCCTCGCGGACGATg                    <  2:1868703/150‑1 (MQ=255)
                                                                                                                                   tCGTCAGCCATGCTCGACAGGTCGCGCAGGCTGGCGGCCACGTCGTCACGCTGCAACGCCGCCGCGGCATCGGCACCGGCATTGCGCTGGGCCATGGCGCCGATCTCGACGCCGGTACGCTGGGCCACCTCGCCGGCCTCGCGGACGATg                    <  2:3523595/150‑1 (MQ=255)
                                                                                                                                    cGTCAGCCATGCTCGACAGGTCGCGCAGGCTGGCGGCCACTTCGTCACGCTGCAACGCCGCCGCGGCATCGGCACCGGCATTGCGCTGGGCCATGGCGCCGATCTCGACGCCGGTACGCTGGGCCACCTCGCCGGCCTCGCGGACGATgg                   <  1:5030107/150‑1 (MQ=255)
                                                                                                                                      tCAGCCATGCTCGACAGGTCGCGCAGGCTGGCGGCCACTTCGTCACGCTGCAACGCCGCCGCGGCATCGGCACTGGCATTGCGCTGGGCCATGGCGCCGATCTCGACGCCGGTACGCTGGGCCACCTCGCCGGCCTCGCGGACGATGGGc                 >  1:1491297/1‑150 (MQ=255)
                                                                                                                                      tCAGCCATGCTCGACAGGTCGCGCAGGCTGGCGGCCACTTCGTCACGCTGCAACGCCGCCGCGGCATCGGCACCGGCATTGCGCTGGGCCATGGCGCCGATCTCGACGCCGGTACGCTGGGCCACCTCGCCGGCCTCGCGGACGATGGGc                 <  2:435189/150‑1 (MQ=255)
                                                                                                                                       cAGCCATGCTCGACAGGTCGCGCAGGCTGGCGGCCACTTCGTCACGCTGCAACGCCGCCGCGGCATCGGCACCGGCATTGCGCTGGGCCATGGCGCCGATCTCGACGCCGGTACGCTGGGCCACCTCGCCGGCCTCGCGGACGATGGGCt                <  1:971195/150‑1 (MQ=255)
                                                                                                                                        aGCCATGCTCGACAGGTCGCGCAGGCTGGCGGCCACTTCGTCACGCTGCAACGCCGCCGCGGCATCGGCACCGGCATTGCGCTGGGCCATGGCGCCGAtctc                                                               >  2:5312634/1‑102 (MQ=255)
                                                                                                                                        aGCCATGCTCGACAGGTCGCGCAGGCTGGCGGCCACTTCGTCACGCTGCAACGCCGCCGCGGCATCGGCACCGGCATTGCGCTGGGCCATGGCGCCGATCTCGACGCCGGTACGCTGGGCCACCTCTCCGGCCTCGCGGACGATGGGCTg               <  2:1190964/150‑1 (MQ=255)
                                                                                                                                        aGCCATGCTCGACAGGTCGCGCAGGCTGGCGGCCACTTCGTCACGCTGCAACGCCGCCGCGGCATCGGCACCGGCATTGCGCTGGGCCATGGCGCCGATCTCGACGCCGGTACGCTGGGCCACCTCGCCGGCCTCGCGGACGATGGGCTg               >  2:574360/1‑150 (MQ=255)
                                                                                                                                           cctGCGCGACCTGTCGAGCATGCTGGCGGCCACTTCGTCACGCTGCAACGCCGCCGCGGCATCGGCACCGGCATTGCGCTGGGCCATGGCGCCGATCTCGACGCCGGTACGCTGGGCCACCTCGCCGGCCTCGCGGACGATGGGCTGCAg            >  2:2136992/3‑150 (MQ=255)
                                                                                                                                            ctGCGCGACCTGTCGAGCATGGCTGGCGGCCACTTCGTCACGCTGCAACGCCGCCGCGGCATCGGCACCGGCATTGCGCTGGGCCATGGCGCCGATCTCGACGCCGGTACGCTGGGCCACCTCGCCGGCCTCGCGGACGATGGGCTGCAg            >  2:3885519/2‑150 (MQ=255)
                                                                                                                                             tGCTCGACAGGTCGCGCAGGCTGGCGGCCACTTCGTCACGCTGCAACGCCGCCGCGGCATCGGCACCGGCATTGCGCTGGGCCATGGCGCCGATCTCGACGCCGGTACGCTGGGCCACCTCGCCGGCCTCGCGGACGATggg                  <  2:4786778/142‑1 (MQ=255)
                                                                                                                                             tGCTCGACAGGTCGCGCAGGCTGGCGGCCACTTCGTCACGCTGCAACGCCGCCGCGGCATCGGCACCGGCATTGCGCTGGGCCATGGCGCCGATCTCGACGCCGGTACGCTGGGCCACCTCGCCGGCCTCGCGGACGATggg                  >  1:4786778/1‑142 (MQ=255)
                                                                                                                                             tGCGCGACCTGTCGAGCAGGCTGGCGGCCACTTCGTCACGCTGCAACGCCGCCGCGGCATCGGCACCGGCAtt                                                                                       <  1:5425273/73‑1 (MQ=255)
                                                                                                                                                tCGACAGGTCGCGCAGGCTGGCGGCCACTTCGTCACGCTGCAACGCCGCCGCGGCATCGGCACCGGCATTGCGCTGGGCCATGGCGCCGATCTCGACGCCGGTACGCTGGGCCACCTCGCCGGCCTCGCGGACGATGGGCTGCAGCTTGt       >  1:3364041/1‑150 (MQ=255)
                                                                                                                                                tCGACAGGTCGCGCAGGCTGGCGGCCACTTCGTCACGCTGCAACGCCGCCGCGGCATCGGCACCGGCATTGCGCTGGGCCATGGCGCCGATCTCGACGCCGGTACGCTGGGCCACCTCGCCGGCCTCGCGGACGATGGGCTGCAGCTTGt       >  2:3595401/1‑150 (MQ=255)
                                                                                                                                                 cGACAGGTCGCGCAGGCTGGCGGCCACTTCGTCACGCTGCAACGCCGCCGCGGCATCGGCACCGGCATTGCGCTGGGCCATGGCGCCGATCTCGACGCCGGTACGCTGGGCCACCTCGCCGGCCTCGCGGACGATGGGCTGCAGCTTGTc      <  1:937702/150‑1 (MQ=255)
                                                                                                                                                 cGACAGGTCGCGCAGGCTGGCGGCCACTTCGTCACGCTGCAACGCCGCCGCGGCATCGGCACCGGCATTGCGCTGGGCCATGGCGCCGATCTCGACGCCGGTACGCTGGGCCACCTCGCCGGCCTCGCGGACGATGGGCTGCAGCTTGTc      <  1:747065/150‑1 (MQ=255)
                                                                                                                                                     aGGTCGCGCAGGCTGGCGGCCACTTCGTCACGCTGCAACGCCGCCGCGGCATCGGCACCGGCATTGCGCTGGGCCATGGCGCCGATCTCGACGCCGGTACGCTGGGCCACCTCGCCGGCCTCGCGGACGATGGGCTGCAGCTTGTCGACa  >  2:297546/1‑150 (MQ=255)
                                                                                                                                                     |                                                                                                                                                      
CAACCTGCCCGGCAAGGTTCTCAATCAGGCTCGCCAGTTGCGTCGAGGAACGGCTGTTCTCATCCGTCGCCTGGCGAATGTCCACCACCTGCTGCAACGCCGCCTGCATGGCATGGCTTTCAGCCTGAGCCTCGTCAGCCATGCTCGACAGGTCGCGCAGGCTGGCGGCCACTTCGTCACGCTGCAACGCCGCCGCGGCATCGGCACCGGCATTGCGCTGGGCCATGGCGCCGATCTCGACGCCGGTACGCTGGGCCACCTCGCCGGCCTCGCGGACGATGGGCTGCAGCTTGTCGACA  >  NC_002947/5558231‑5558529

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: