Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 5,367,686 G→C 9.0% L280L (CTC→CTG ftsH ← integral membrane ATP‑dependent zinc metallopeptidase

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0029475,367,6860GC 9.0% 576.2 / 36.5 178L280L (CTC→CTGftsHintegral membrane ATP‑dependent zinc metallopeptidase
Reads supporting (aligned to +/- strand):  ref base G (86/76);  new base C (11/5);  total (97/81)
Fisher's exact test for biased strand distribution p-value = 2.96e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.47e-01

Only 100 of 178 total aligned reads displayed.

GGATGTCCGGCAGGCCAACCACCACCTGGCGGTCGAAGCGGCCGGGGCGCAGCAGCGCAGGGTCGAGCACATCCGGACGGTTGGTGGCGGCGATGACGATGATGCCATCATTCATCTCGAAGCCGTCCATCTCCACCAGCAACTGGTTGAGGGTTTGCTCACGCTCGTCGTGACCGCCGCCCATGCCGGCGCCACGATGGCGACCAACGGCGTCGATTTCGTCGATGAAGATGATGCACGGGGCGTGCTTCTTGGCCTGCTCGAACATGTCGCGCACACGGCTGGCACCGACACCCA  >  NC_002947/5367538‑5367834
                                                                                                                                                    |                                                                                                                                                    
ggATGTCCGGCAGGCCAACCACCACCTGGCGGTCGAAGCGGCCGGGGCGCAGCAGCGCAGGGTCGAGCACATCCGGACGGTTGGTGGCGGCGATGACGATGATGCCATCATTCATCTCGAAGCCGTCCATCTCCACCAGCAACTGGTTGa                                                                                                                                                     >  2:917174/1‑150 (MQ=255)
ggATGTCCGGCAGGCCAACCACCACCTGGCGGTCGAAGCGGCCGGGGCGCAGCAGCGCAGGGTCGAGCACATCCGGACGGTTGGTGGCGGCGATGACGATGATGCCATCATTCATCTCGAAGCCGTCCATCTCCACCAGCAACTGGTTGa                                                                                                                                                     >  1:1612980/1‑150 (MQ=255)
 gATGTCCGGCAGGCCAACCACCACCTGGCGGTCGAAGCGGCCGGGGCGCAGCAGCGCAGGGTCGAGCACATCCGGACGGTTGGTGGCGGCGATGACGATGATGCCATCATTCATCTCGAAGCCGTCCATCTCCACCAGCAACTGGTTGAg                                                                                                                                                    >  2:1454551/1‑150 (MQ=255)
   tGTCCGGCAGGCCAACCACCACCTGGCGGTCGAAGCGGCCGGGGCGCAGCAGCGCAGGGTCGAGCACATCCGGACGGTTGGTGGCGGCGATGACGATGATGCCATCATTCATCTCGAAGCCGTCCATCTCCACCAGCAACTGGTTGAggg                                                                                                                                                  >  2:1054362/1‑150 (MQ=255)
     tCCGGCAGGCCAACCACCACCTGGCGGTCGAAGCGGCCGGGGCGCAGCAGCGCAGGGTCGAGCACATCCGGACGGTTGGTGGCGGCGATGACGATGATGCCATCATTCATCTCGAAGCCGTCCATCTCCACCAGCAACTGGTTGAGGGtt                                                                                                                                                >  1:3252053/1‑150 (MQ=255)
     tCCGGCAGGCCAACCACCACCTGGCGGTCGAAGCGGCCGGGGCGCAGCAGCGCAGGGTCGAGCACATCCGGACGGTTGGTGGCGGCGATGACGATGATGCCATCATTCATCTCGAAGCCGTCCATCTCCACCAGCAACTGGTTGAGGGtt                                                                                                                                                <  2:3525885/150‑1 (MQ=255)
      ccGGCAGGCCAACCACCACCTGGCGGTCGAAGCGGCCGGGGCGCAGCAGCGCAGGGTCGAGCACATCCGGACGGTTGGTGGCGGCGATGACGATGATGCCATCATTCATCTCGAAGCCGTCCATCTCCACCAGCAACTGGTTGAGGGttt                                                                                                                                               <  2:3942505/150‑1 (MQ=255)
       cGGCAGGCCAACCACCACCTGGCGGTCGAAGCGGCCGGGGCGCAGCAGCGCAGGGTCGAGCACATCCGGACGGTTGGTGGCGGCGATGACGATGATGCCATCATTCATCTCGAAGCCGTCCATCTCCACCAGCAACTGGTTGAGGGTTTg                                                                                                                                              <  2:1900579/150‑1 (MQ=255)
        ggCAGGCCAACCACCACCTGGCGGTCGAAGCGGCCGGGGCGCAGCAGCGCAGGGTCGAGCACATCCGGACGGTTGGTGGCGGCGATGACGATGATGCCATCATTCATCTCGAAGCCGTCCATCTCCACCAGCAACTGGTTGAGGGTTTGc                                                                                                                                             >  2:3945296/1‑150 (MQ=255)
             gCCAACCACCACCTGGCGGTCGAAGCGGCCGGGGCGCAGCTGCGCAGGGTCGAGCACATCCGGACGGTTGGTGGCGGCGATGACGATGATGCCATCATTCATCTCGAAGCCGTCCATCTCCACCAGCAACTGGTTGAGGGTTTGCTCACg                                                                                                                                        >  2:3802804/1‑150 (MQ=255)
              ccAACCACCACCTGGCGGTCGAAGCGGCCGGGGCGCAGCAGCGCAGGGTCGAGCACATCCGGACGGTTGGTGGCGGCGATGACGATGATGCCATCATTCATCTCGAAGCCGTCCATCTCCACCAGCAACTGGTTGAGGGTTTGCTCACGc                                                                                                                                       >  1:2114336/1‑150 (MQ=255)
                aaCCACCACCTGGCGGTCGAAGCGGCCGGGGCGCAGCAGCGCAGGGTCGAGCACATCCGGACGGTTGGTGGCGGCGATGACGATGATGCCATCATTCATCTCGAAGCCGTCCATCTCCACCAGCAACTGGGTGAGGGTTTGCTCACGCtc                                                                                                                                     >  2:3652638/1‑150 (MQ=255)
                  ccaccaccTGGCGGTCGAAGCGGCCGGGGCGCAGCAGCGCAGGGTCGAGCACATCCGGACGGTTGGTGGCGGCGATGACGATGATGCCATCATTCATCTCGAAGCCGTCCATCTCCACCAGCAACTGGTTGAGGGTTTGCTCACGCtcgt                                                                                                                                   <  2:4216326/150‑1 (MQ=255)
                   caccaccTGGCGGTCGAAGCGGCCGGGGCGCAGCAGCGCAGGGTCGAGCACATCCGGACGGTTGGTGGCGGCGATGACGATGATGCCATCATTCATCTCGAAGCCGTCCATCTCCACCAGCAACTGGTTGAGGGTTTGCTCACGCtcgtc                                                                                                                                  <  1:1876802/150‑1 (MQ=255)
                     ccaccTGGCGGTCGAAGCGGCCGGGGCGCAGCAGCGCAGGGTCGAGCACATCCGGACGGTTGGTGGCGGCGATGACGATGATGCCATCATTCATCTCGAAGCCGTCCATCTCCACCAGCAACTGGTTGAGGGt                                                                                                                                                 >  1:3862370/1‑133 (MQ=255)
                     ccaccTGGCGGTCGAAGCGGCCGGGGCGCAGCAGCGCAGGGTCGAGCACATCCGGACGGTTGGTGGCGGCGATGACGATGATGCCATCATTCATCTCGAAGCCGTCCATCTCCACCAGCAACTGGTTGAGGGt                                                                                                                                                 <  2:3862370/133‑1 (MQ=255)
                     ccaccTGGCGGTCGAAGCGGCCGGGGCGCAGCAGCGCAGGGTCGAGCACATCCGGACGGTTGGTGGCGGCGATGACGATGATGCCATCATTCATCTCGAAGCCGTCCATCTCCACCAGCAACTGGTTGAGGGTTTGCTCACGCTcgtcgt                                                                                                                                <  1:3209983/150‑1 (MQ=255)
                      caccTGGCGGTCGAAGCGGCCGGGGCGCAGCAGCGCAGGGTCGAGCACATCCGGACGGTTGGTGGCGGCGATGACGATGATGCCATCATTCATCTCGAAGCCGTCCATCTCCACCAGCAACTGGTTGAGGGTTTGCTCACGCTCGTCGTg                                                                                                                               >  1:2448675/1‑150 (MQ=255)
                      caccTGGCGGTCGAAGCGGCCGGGGCGCAGCAGCGCAGGGTCGAGCACATCCGGACGGTTGGTGGCGGCGATGACGATGATGCCATCATTCATCTCGAAGCCGTCCATCTCCACCAGCAACTGGTTGAGGGTTTGCTCACGCTCGTCGTg                                                                                                                               >  2:280955/1‑150 (MQ=255)
                      caccTGGCGGTCGAAGCGGCCGGGGCGCAGCAGCGCAGGGTCGAGCACATCCGGACGGTTGGTGGCGGCGATGACGATGATGCCATCATTCATCTCGAAGCCGTCCATCTCCACCAGCAACTGGTTGAGGGTTTGCTCACGCTCGTCGTg                                                                                                                               >  1:570294/1‑150 (MQ=255)
                      cacaTGGCCGTCCAAGCGGCCGGCGCGCAGCAGCGCACGGTCGAGCACATCCGGACGGTTGGTGGCGGCGATGACGATGATGCCATCATTCATATCGAAGCCGACCAAATCCACCAGCAACAGGTTGAGGGATTGCTCACGCTCGTCGTg                                                                                                                               >  2:281692/1‑150 (MQ=255)
                              ggTCGAAGCGGCCGGGGCGCAGCAGCGCAGGGTCGAGCACATCCGGACGGTTGGTGGCGGCGATGACGATGATGCCATCATTCATCTCGAAGCCGTCCATCTCCACCAGCAACTGGTTGAGGGTTTGCTCACGCTCGTCGTGACcgccgc                                                                                                                       >  2:2155473/1‑150 (MQ=255)
                                tCGAAGCGGCCGGGGCGCAGCAGCGCAGGGTCGAGCACATCCGGACGGTTGGTGGCGGCGATGACGATGATGCCATCATTCATCTCGAAGCCGTCCATCTCCACCAGCAACTGGTTGAGGGTTTGCTCACGCTCGTCGTGACCGCCGccc                                                                                                                     >  2:2144033/1‑150 (MQ=255)
                                tCGAAGCGGCCGGGGCGCAGCAGCGCAGGGTCGAGCACATCCGGACGGTTGGTGGCGGCGATGACGATGATGCCATCATTCATCTCGAAGCCGTCCATCTCCACCAGCAACTGGTTGAGGGTTTGCTCACGCTCGTCGTGACCGCCGccc                                                                                                                     >  2:4724920/1‑150 (MQ=255)
                                  gAAGCGGCCGGGGCGCAGCAGCGCAGGGTCGAGCACATCCGGACGGTTGGTGGCGGCGATGACGATGATGCCATCATTCATCTCGAAGCCGTCCATCTCCACCAGCAACTGGTTGAGGGTTTGCTCACGCTCGTCGTGACCGCCGCCCAt                                                                                                                   >  1:2549723/1‑150 (MQ=255)
                                   aaGCGGCCGGGGCGCAGCAGCGCAGGGTCGAGCACATCCGGACGGTTGGTGGCGGCGATGACGATGATGCCATCATTCATCTCGAAGCCGTCCATCTCCACCAGCAACTGGTTGAGGGTTTGCTCACGCTCGTCGTGACCGCCGCCCATg                                                                                                                  <  2:2241554/150‑1 (MQ=255)
                                    aGCGGCCGGGGCGCAGCAGCGCAGGGTCGAGCACATCCGGACGGTTGGTGGCGGCGATGACGATGATGCCATCATTCATCTCGAAGCCGTCCATCTCCACCAGCAACTGGTTGAGGGTTTGCTCACGCTCGTCGTGACCGCCGCCCATGc                                                                                                                 <  2:4322964/150‑1 (MQ=255)
                                        gaccctgcgcTGCTGCGCAGGGTCGAGCACATCCGGACGGTTGGTGGCGGCGATGACGATGATGCCATCATTCATCTCGAAGCCGTCCATCTCCACCAGCAACTGGTTGAGGGTTTGCTCACGCTCGTCGTGACCGCCGCCCATGCCGgc                                                                                                             >  2:2531581/7‑150 (MQ=255)
                                               cgcAGCAGCGCAGGGTCGAGCACATCCGGACGGTTGGTGGCGGCGATGACGATGATGCCATCATTCATCTCGAAGCCGTCCATCTCCACCAGCAACTGGTTGAGGGTTTGCTCACGCTCGTCGTGACCGCCGCCCATGCCGGCGCCACGa                                                                                                      <  2:1169112/150‑1 (MQ=255)
                                               cgcAGCAGCGCAGGGTCGAGCACATCCGGACGGTTGGTGGCGGCGATGACGATGATGCCATCATTCATCTCGAAGCCGTCCATCTCCACCAGCAACTGGTTGAGGGTTTGCTCACGCTCGTCGTGACCGCCGCCCATGCCGGCGCCACGa                                                                                                      <  2:4111980/150‑1 (MQ=255)
                                                gcagcaGCGCAGGGTCGAGCACATCCGGACGGTTGGTGGCGGCGATGACGATGATGCCATCATTCATCTCGAAGCCGTCCATCTCCACCAGCAACTGGTTGAGGGTTTGCTCACGCTCGTCGTGACCGCCGCCCATGCCGGCGCCACGAt                                                                                                     <  1:726114/150‑1 (MQ=255)
                                                  agcagcGCAGGGTCGAGCACATCCGGACGGTTGGTGGCGGCGATGACGATGATGCCATCATTCATCTCGAAGCCGTCCATCTCCACCAGCAACTGGTTGAGGGTTTGCTCACGCTCGTCGTGACCGCCGCCCATGCCGGCGCCACGATgg                                                                                                   >  2:1475714/1‑150 (MQ=255)
                                                  agcagcGCAGGGTCGAGCACATCCGGACGGTTGGTGGCGGCGATGACGATGATGCCATCATTCATCTCGAAGCCGTCCATCTCCACCAGCAACTGGTTGAGGGTTTGCTCACGCTCGTCGTGACCGCCGCCCATGCCGGCGCCACGATgg                                                                                                   >  2:360793/1‑150 (MQ=255)
                                                   gcagcGCAGGGTCGAGCACATCCGGACGGTTGGTGGCGGCGATGACGATGATGCCATCATTCATCTCGAAGCCGTCCATCTCCACCAGCAACTGGTTGAGGGTTTGCTCACGCTCGTCGTGACCGCCGCCCATGCCGGCGCCACGATGGc                                                                                                  <  1:1258730/150‑1 (MQ=255)
                                                      gcgcAGGGTCGAGCACATCCGGACGGTTGGTGGCGGCGATGACGATGATGCCATCATTCATCTCGAAGCCGTCCATCTCCACCAGCAACTGGTTGAGGGTTTGCTCACGCTCGTCGTGACCGCCGCCCATGCCGGCGCCACGATGGCGAc                                                                                               <  1:4622710/150‑1 (MQ=255)
                                                           gggTCGAGCACATCCGGACGGTTGGTGGCGGCGATGACGATGATGCCATCATTCATCTCGAAGCCGTCCATCTCCACCAGCAACTGGTTGAGGGTTTGCTCACGCTCGTCGTGACCGCCGCCCATGCCGGCGCCACGATGGCGACCAACg                                                                                          <  1:1974723/150‑1 (MQ=255)
                                                                 aGCACATCCGGACGGTTGGTGGCGGCGATGACGATGATGCCATCATTCATCTCGAAGCCGTCCATCTCCACCAGCAACTGGTTGAGGGTTTGCTCACGCTCGTCGTGACCGCCGCCCATGCCGGCGCCACGATGGCGACCAACGGCGTCg                                                                                    <  1:3054059/150‑1 (MQ=255)
                                                                   cacaTCCGGACGGTTGGTGGCGGCGATGACGATGATGCCATCATTCATCTCGAAGCCGTCCATCTCCACCAGCAACTGGTTGAGGGTTTGCTCACGCTCGTCGTGACCGCCGCCCATGCCGGCGCCACGATGGCGACCAACGGCGTCGAt                                                                                  <  1:1583569/150‑1 (MQ=255)
                                                                      aTCCGGACGGTTGGTGGCGGCGATGACGATGATGCCATCATTCATCTCGAAGCCGTCCATCTCCACCAGCAACTGGTTGAGGGTTTGCTCACGCTCGTCGTGACCGCCGCCCATGCCGGCGCCACGATGGCGACCAACGGCGTCGATTtc                                                                               <  2:3739730/150‑1 (MQ=255)
                                                                      aTCCGGACGGTTGGTGGCGGCGATGACGATGATGCCATCATTCATCTCGAAGCCGTCCATCTCCACCAGCAACTGGTTGAGGGTTTGCTCACGCTCGTCGTGACCGCCGCCCATGCCGGCGCCACGATGGCGACCAACGGCGTCGATTtc                                                                               <  1:997296/150‑1 (MQ=255)
                                                                                ttGGTGGCGGCGATGACGATGATGCCATCATTCATCTCGAAGCCGTCCATCTCCACCAGCAACTGGTTGAGGGTTTGCTCACGCTCGTCGTGACCGCCGCCCATGCCGGCGCCACGATGGCGACCAACGGCGTCGATTTCGTCGATGAAg                                                                     <  1:2144033/150‑1 (MQ=255)
                                                                                   gtggCGGCGATGACGATGATGCCATCATTCATCTCGAAGCCGTCCATCTCCACCAGCAACTGGTTGAGGGTTTGCTCACGCTCGTCGTGACCGCCGCCCATGCCGGCGCCACGATGGCGACCAACGGCGTCGATTTCGTCGATGAAgatg                                                                  <  2:1918926/150‑1 (MQ=255)
                                                                                    tggCGGCGATGACGATGATGCCATCATTCATCTCGAAGCCGTCCATCTCCACCAGCAACTGGTTGAGGGTTTGCTCACGCTCGTCGTGACCGCCGCCCATGCCGGCGCCACGATGGCGACCAACGGCGTCGATTTCGTCGATGAAgatga                                                                 >  1:3006194/1‑150 (MQ=255)
                                                                                    tggCGGCGATGACGATGATGCCATCATTCATCTCGAAGCCGTCCATCTCCACCAGCAACTGGTTGAGGGTTTGCTCACGCTCGTCGTGACCGCCGCCCATGCCGGCGCCACGATGGCGACCAACGGCGTCGATTTCGTCGATGAAgatga                                                                 <  1:2838218/150‑1 (MQ=255)
                                                                                    tggCGGCGATGACGATGATGCCATCATTCATCTCGAAGCCGTCCATCTCCACCAGCAACTGGTTGAGGGTTTGCTCACGCTCGTCGTGACCGCCGCCCATGCCGGCGCCACGATGGCGACCAACGGCGTCGATTTCGTCGATGAAgatga                                                                 >  2:1428559/1‑150 (MQ=255)
                                                                                            atgacgatgaTGCCATCATTCATCTCGAAGCCGTCCATCTCCACCAGCAACTGGTTGAGGGTTTGCTCACGCTCGTCGTGACCGCCGCCCATGCCGGCGCCACGATGGCGACCAACGGCGTCGATTTCGTCGATGAAGATGATGCACggg                                                         <  1:4771162/150‑1 (MQ=255)
                                                                                            atgacgatgaTGCCATCATTCATCTCGAAGCCGTCCATCTCCACCAGCAACTGGTTGAGGGTTTGCTCACGCTCGTCGTGACCGCCGCCCATGCCGGCGCCACGATGGCGACCAACGGCGTCGATTTCGTCGATGAAGATGATGCACggg                                                         <  2:1212993/150‑1 (MQ=255)
                                                                                                 gatgatGCCATCATTCATCTCGAAGCCGTCCATCTCCACCAGCAACTGGTTGAGGGTTTGCTCACGCTCGTCGTGACCGCCGCCCATGCCGGCGCCACGATGGCGACCAACGGCGTCGATTTCGTCGATGAAGATGATGCACGGGGCGTg                                                    >  1:4982645/1‑150 (MQ=255)
                                                                                                 gatgatGCCATCATTCATCTCGAAGCCGTCCATCTCCACCAGCAACTGGTTGAGGGTTTGCTCACGCTCGTCGTGACCGCCGCCCATGCCGGCGCCACGATGGCGACCAACGGCGTCGATTTCGTCGATGAAGATGATGCACGGGGCGTg                                                    >  2:1580652/1‑150 (MQ=255)
                                                                                                    gatgCCATCATTCATCTCGAAGCCGTCCATCTCCACCAGCAACTGGTTGAGGGTTTGCTCACGCTCGTCGTGACCGCCGCCCATGCCGGCGCCACGATGGCGACCAACGGCGTCGATTTCGTCGATGAAGATGATGCACGGGGCGTGctt                                                 >  1:2742235/1‑150 (MQ=255)
                                                                                                      tgCCATCATTCATCTCGAAGCCGTCCATCTCCACCAGCAACTGGTTGAGGGTTTGCTCACGCTCGTCGTGACCGCCGCCCATGCCGGCGCCACGATGGCGACCAACGGCGTCGATTTCGTCGATGAAGATGATGCACGGGGCGTGcttct                                               >  1:5036089/1‑150 (MQ=255)
                                                                                                      tgCCATCATTCATCTCGAAGCCGTCCATCTCCACCAGCAACTGGTTGAGGGTTTGCTCACGCTCGTCGTGACCGCCGCCCATGCCGGCGCCACGATGGCGACCAACGGCGTCGATTTCGTCGATGAAGATGATGCACGGGGCGTGcttct                                               <  2:3994282/150‑1 (MQ=255)
                                                                                                       gCCATCATTCATCTCGAAGCCGTCCATCTCCACCAGCAACTGGTTGAGGGTTTGCTCACGCTCGTCGTGACCGCCGCCCATGCCGGCGCCACGATGGCGACCAACGGCGTCGATTTCGTCGATGAAGATGATGCACGGGGCGTGcttctt                                              <  2:2434689/150‑1 (MQ=255)
                                                                                                       gCCATCATTCATCTCGAAGCCGTCCATCTCCACCAGCAACTGGTTGAGGGTTTGCTCACGCTCGTCGTGACCGCCGCCCATGCCGGCGCCACGATGGCGACCAACGGCGTCGATTTCGTCGATGAAGATGATGCACGGGGCGTGcttctt                                              <  2:3730254/150‑1 (MQ=255)
                                                                                                        ccATCATTCATCTCGAAGCCGTCCATCTCCACCAGCAACTGGTTGAGGGTTTGCTCACGCTCGTCGTGACCGCCGCCCATGCCGGCGCCACGATGGCGACCAACGGCGTCGATTTCGTCGATGAAGATGATGCACGGGGCGTGCTTCTTg                                             <  1:1063050/150‑1 (MQ=255)
                                                                                                        ccATCATTCATCTCGAAGCCGTCCATCTCCACCAGCAACTGGTTGAGGGTTTGCTCACGCTCGTCGTGACCGCCGCCCATGCCGGCGCCACGATGGCGACCAACGGCGTCGATTTCGTCGATGAAGATGATGCACGGGGCGTGCTTCTTg                                             <  1:917174/150‑1 (MQ=255)
                                                                                                        ccATCATTCATCTCGAAGCCGTCCATCTCCACCAGCAACTGGTTGAGGGTTTGCTCACGCTCGTCGTGACCGCCGCCCATGCCGGCGCCACGATGGCGACCAACGGCGTCGATTTCGTCGATGAAGATGATGCACGGGGCGTGCTTCTTg                                             >  2:2569921/1‑150 (MQ=255)
                                                                                                        ccATCATTCATCTCGAAGCCGTCCATCTCCACCAGCAACTGGTTGAGGGTTTGCTCACGCTCGTCGTGACCGCCGCCCATGCCGGCGCCACGATGGCGACCAACGGCGTCGATTTCGTCGATGAAGATGATGCACGGGGCGTGCTTCTTg                                             <  1:5084572/150‑1 (MQ=255)
                                                                                                         catcatTCATCTCGAAGCAGTCCATCTCCACCAGCAACTGGTTGAGGGTTTGCTCACTCTCGTCGTGACCTCCGCCCATGCCTGCGCCACTATGGCTACCAACGGCGTCGATTTCGTCGATGAAGATGATGCACGGTGCGTGCTTCTTgg                                            <  2:1584039/150‑1 (MQ=255)
                                                                                                          atcatTCATCTCGAAGCCGTCCATCTCCACCAGCAACTGGTTGAGGGTTTGCTCACGCTCGTCGTGACCGCCGCCCATGCCGGCGCCACGATGGCGACCAACGGCGTCGATTTCGTCGATGAAGATGATGCACGGGGCGTGCTTCTTGGc                                           >  2:4392713/1‑150 (MQ=255)
                                                                                                          atcatTCATCTCGAAGCCGTCCATCTCCACCAGCAACTGGTTGAGGGTTTGCTCACGCTCGTCGTGACCGCCGCCCATGCCGGCGCCACGATGGCGACCAACGGCGTCGATTTCGTCGATGAAGATGATGCACGGGGCGTGCTGCTTGGc                                           >  2:4946581/1‑150 (MQ=255)
                                                                                                              ttcatCTCGAAGCCGTCCATCTCCACCAGCAACTGGTTGAGGGTTTGCTCACGCTCGTCGTGACCGCCGCCCATGCCGGCGCCACGATGGCGACCAACGGCGTCGATTTCGTCGATGAAGATGATGCACGGGGCGTGcttcttcttcatc                                       <  1:4487547/150‑8 (MQ=255)
                                                                                                              ttcatCTCGAAGCCGTCCATCTCCACCAGCAACTGGTTGAGGGTTTGCTCACGCTCGTCGTGACCGCCGCCCATGCCGGCGCCACGATGGCGACCAACGGCGTCGATTTCGTCGATGAAGATGATGCACGGGGCGTGCTTCTTGGCCTGc                                       <  2:4792486/150‑1 (MQ=255)
                                                                                                               tcatCTCGAAGCCGTCCATCTCCACCAGCAACTGGTTGAGGGTTTGCTCACGCTCGTCGTGACCGCCGCCCATGCCGGCGCCACGATGGCGACCAACGGCGTCGATTTCGTCGATGAAGATGATGCACGGGGCGTGCTTCTTGGCCTGCt                                      >  1:3718315/1‑150 (MQ=255)
                                                                                                                 atCTCGAAGCCGTCCATCTCCACCAGCAACTGGTTGAGGGTTTGCTCACGCTCGTCGTGACCGCCGCCCATGCCGGCGCCACGATGGCGACCAACGGCGTCGATTTCGTCGATGAAGATGATGCACGGGGCGTGCTTCTTGGCCTGCTCg                                    >  1:421296/1‑150 (MQ=255)
                                                                                                                  tctcGAAGCCGTCCATCTCCACCAGCAACTGGTTGAGGGTTTGCTCACGCTCGTCGTGACCGCCGCCCATGCCGGCGCCACGATGGCGACCAACGGCGTCGATTTCGTCGATGAAGATGATGCACGGGGCGTGCTTCTTGGCCTGCTCGa                                   >  2:4522051/1‑150 (MQ=255)
                                                                                                                   ctcGAAGCCGTCCATCTCCACCAGCAACTGGTTGAGGGTTTGCTCACGCTCGTCGTGACCGCCGCCCATGCCGGCGCCACGATGGCGACCAACGGCGTCGATTTCGTCGATGAAGATGATGCACGGGGCGTGCTTCTTGGCCTGCTCGaa                                  <  1:3662294/150‑1 (MQ=255)
                                                                                                                    tcGAAGCCGTCCATCTCCACCAGCAACTGGTTGAGGGTTTGCTCACGCTCGTCGTGACCGCCGCCCATGCCGGCGCCACGATGGCGACCAACGGCGTCGATTTCGTCGATGAAGATGATGCACGGGGCGTGCTTCTTGGCCTGCTCGAAc                                 >  1:3915792/1‑150 (MQ=255)
                                                                                                                    tcGAAGCCGTCCATCTCCACCAGCAACTGGTTGAGGGTTTGCTCACGCTCGTCGTGACCGCCGCCCATGCCGGCGCCACGATGGCGACCAACGGCGTCGATTTCGTCGATGAAGATGATGCACGGGGCGTGCTTCTTGGCCTGCTCGAAc                                 >  2:3389958/1‑150 (MQ=255)
                                                                                                                    tcGAAGCCGTCCATCTCCACCAGCAACTGGTTGAGGGTTTGCTCACGCTCGTCGTGACCGCCGCCCATGCCGGCGCCACGATGGCGACCAACGGCGTCGATTTCGTCGATGAAGATGATGCACGGGGCGTGCTTCTTGGCCTGCTCGAAc                                 >  2:3391952/1‑150 (MQ=255)
                                                                                                                           catCCATTTCGACCAGAAGCTGGTTCAGCGTCTGTTCGCGTTCGTCGTTGCCGCCGCCAAGACCGGCACCACGGTGGCGACCGACGGCGTCGATTTCGTCGATGAAGATGATACAAGGCGCATTCTTCTTGGCCTGTTCGAACATGTcac                          <  1:1454048/148‑3 (MQ=255)
                                                                                                                           cGTCCATCTCCACCAGCAACTGGTTGAGGGTTTGCTCACGCTCGTCGTGACCGCCGCCCATGCCGGCGCCACGATGGCGACCAACGGCGTCGATTTCGTCGATGAAGATGATGCACGGGGCGTGCTTCTTGGCCTGCTCGAACATGTcgc                          <  1:1329732/150‑1 (MQ=255)
                                                                                                                           cGTCCATCTCCACCAGCAACTGGTTGAGGGTTTGCTCACGCTCGTCGTGACCGCCGCCCATGCCGGCGCCACGATGGCGACCAACGGCGTCGATTTCGTCGATGAAGATGATGCACGGGGCGTGCTTCTTGGCCTGCTCGAACATGTcgc                          <  1:1884641/150‑1 (MQ=255)
                                                                                                                            gTCCATCTCCACCAGCAACTGGTTGAGGGTTTGCTCACGCTCGTCGTGACCGCCGCCCATGCCGGCGCCACGATGGCGACCAACGGCGTCGATTTCGTCGATGAAGATGATGCACGGGGCGTGCTTCTTGGCCTGCTCGAACATGTcgcg                         >  1:1838527/1‑150 (MQ=255)
                                                                                                                            gTCCATCTCCACCAGCAACTGGTTGAGGGTTTGCTCACGCTCGTCGTGACCGCCGCCCATGCCGGCGCCACGATGGCGACCAACGGCGTCGATTTCGTCGATGAAGATGATGCACGGGGCGTGCTTCTTGGCCTGCTCGAACATGTcgcg                         >  2:4949851/1‑150 (MQ=255)
                                                                                                                            gTCCATCTCCACCAGCAACTGGTTGAGGGTTTGCTCACGCTCGTCGTGACCGCCGCCCATGCCGGCGCCACGATGGCGACCAACGGCGTCGATTTCGTCGATGAAGATGATGCACGGGGCGTGCTTCTTGGCCTGCTCGAACATGTcgcg                         >  2:3452989/1‑150 (MQ=255)
                                                                                                                             tCCATTTCGACCAGAAGCTGGTTCAGCGTCTGTTCGCGTTCGTCGTTGCCGCCGCCAAGACCGGCACCACGGTGGCGACCGACGGCGTCGATTTCGTCGATGAAGATGATGCAAGGCGCATTCTTCTTGGCCTGTTCGAACATGTcacgg                        >  2:1757082/1‑149 (MQ=255)
                                                                                                                             tCCATCTCCACCAGCAACTGGTTGAGGGTTTGCTCACGCTCGTCGTGACCGCCGCCCATGCCGGCGCCACGATGGCGACCAACGGCGTCGATTTCGTCGATGAAGATGATGCACGGGGCGTGCTTCTTGGCCTGCTCGAACATGTCgcgc                        >  1:3177378/1‑150 (MQ=255)
                                                                                                                               cATTTCGACCAGAAGCTGGTTCAGCGTCTGTTCGCGTTCGTCGTTGCCGCCGCCAAGACCGGCACCACGGTGGCGACCGACGGCGTCGATTTCGTCGATGAAGATGATGCAAGGCGCATTCTTCTTGGCCTGTTCGAACATGTCACGgac                      >  1:4182869/1‑150 (MQ=255)
                                                                                                                               cATTTCGACCAGAAGCTGGTTCAGCGTCTGTTCGCGTTCGTCGTTGCCGCCGCCAAGACCGGCACCACGGTGGCGACCGACGGCGTCGATTTCGTCGATGAAGATGATGCAAGGCGCATTCTTCTTGGCCTGTTCGAACATGTCACGgac                      >  1:4785506/1‑150 (MQ=255)
                                                                                                                                 tctcCACCAGCAACTGGTTGAGGGTTTGCTCACGCTCGTCGTGACCGCCGCCCATGCCGGCGCCTCGATGGCGACCAACGGCGTCGTTTTCGTCGATGAAGATGATGCACGGGGCGTGCTTCTTGGCCTGCTCGAACATGTCGCGCacac                    <  1:3802804/150‑1 (MQ=255)
                                                                                                                                  ctcCACCAGCAACTGGTTGAGGGTTTGCTCACGCTCGTCGTGACCGCCGCCCATGCCGGCGCCACGATGGCGACCAACGGCGTCGATTTCGTCGATGAAGATGATGCACGGGGCGTGCTTCTTGGCCTGCTCGAACATGTCGCGCACACg                   >  1:2623926/1‑150 (MQ=255)
                                                                                                                                   tcCACCAGCAACTGGTTGAGGGTTTGCTCACGCTCGTCGTGACCGCCGCCCATGCCGGCGCCACGATGGCGACCAACGGCGTCGATTTCGTCGATGAAGATGATGCACGGGGCGTGCTTCTTGGCCTGCTCGAACATGTCGCGCACACgg                  <  1:3605761/150‑1 (MQ=255)
                                                                                                                                    cgaccaGAAGCTGGTTCAGCGTCTGTTCGCGTTCGTCGTTGCCGCCGCCAAGACCGGCACCACGGTGGCGACCGACGGCGTCGATTTCGTCGATGAAGATGATGCAAGGCGCATTCTTCTTGGCCTGTTCGAACATGTCACGGACGCGGc                 <  2:531864/148‑1 (MQ=255)
                                                                                                                                     caccaGCAACTGGTTGAGGGTTTGCTCACGCTCGTCGTGACCGCCGCCCATGCCGGCGCCACGATGGCGACCAACGGCGTCGATTTCGTCGATGAAGATGATGCACggg                                                         <  2:1925412/109‑1 (MQ=255)
                                                                                                                                     caccaGCAACTGGTTGAGGGTTTGCTCACGCTCGTCGTGACCGCCGCCCATGCCGGCGCCACGATGGCGACCAACGGCGTCGATTTCGTCGATGAAGATGATGCACGGGGCGTGCTTCTTGGCCTGCTCGAACATGTCGCGCACACGGCt                >  2:3091246/1‑150 (MQ=255)
                                                                                                                                     caccaGCAACTGGTTGAGGGTTTGCTCACGCTCGTCGTGACCGCCGCCCATGCCGGCGCCACGATGGCGACCAACGGCGTCGATTTCGTCGATGAAGATGATGCACGGGGCGTGCTTCTTGGCCTGCTCGAACATGTCGCGCACACGGCt                >  1:4082998/1‑150 (MQ=255)
                                                                                                                                     caccaGCAACTGGTTGAGGGTTTGCTCACGCTCGTCGTGACCGCCGCCCATGCCGGCGCCACGATGGCGACCAACGGCGTCGATTTCGTCGATGAAGATGATGCACGGGGCGTGCTTCTTGGCCTGCTCGAACATGTCGCGCACACGGCt                <  1:3945296/150‑1 (MQ=255)
                                                                                                                                      accaGCAACTGGTTGAGGGTTTGCTCACGCTCGTCGTGACCGCCGCCCATGCCGGCGCCACGATGGCGACCAACGGCGTCGATTTCGTCGATGAAGATGATGCACGGGGCGTGCTTCTTGGCCTGCTCGAACATGTCGCGCACACGGCTg               <  1:1454551/150‑1 (MQ=255)
                                                                                                                                      accaGCAACTGGTTGAGGGTTTGCTCACGCTCGTCGTGACCGCCGCCCATGCCGGCGCCACGATGGCGACCAACGGCGTCGATTTCGTCGATGAAGATGATGCACGGGGCGTGCTTCTTGGCCTGCTCGAACATGTCGCGCACACGGCTg               <  2:1505752/150‑1 (MQ=255)
                                                                                                                                      accaGAAGCTGGTTCAGCGTCTGTTCGCGTTCGTCGTTGCCGCCGCCAAGACCGGCACCACGGTGGCGACCGACGGCGTCGATTTCGTCGATGAAGATGATGCAAGGCGCATTCTTCTTGGCCTGTTCGAACATGTCACGGACGCGGCtt               <  1:1107363/150‑2 (MQ=255)
                                                                                                                                       ccaGAAGCTGGTTCAGCGTCTGTTCGCGTTCGTCGTTGCCGCCGCCAAGACCGGCACCACGGTGGCGACCGACGGCGTCGATTTCGTCGATGAAGATGATGCAAGGCGCATTctt                                                 >  2:2353317/1‑115 (MQ=255)
                                                                                                                                       ccaGAAGCTGGTTCAGCGTCTGTTCGCGTTCGTCGTTGCCGCCGCCAAGACCGGCACCACGGTGGCGACCGACGGCGTCGATTTCGTCGATGAAGATGATGCAAGGCGCATTctt                                                 <  1:2353317/115‑1 (MQ=255)
                                                                                                                                        cagaagcTGGTTCAGCGTCTGTTCGCGTTCGTCGTTGCCGCCGCCAAGACCGGCACCACGGTGGCGACCGACGGCGTCGATTTCGTCGATGAAGATGATGCAAGGCGCATTCTTCTTGGCCTGTTCGAACATGTCACGGACGCGGCTTGc             >  2:5022872/7‑150 (MQ=255)
                                                                                                                                        cagaagcTGGTTCAGCGTCTGTTCGCGTTCGTCGTTGCCGCCGCCAAGACCGGCACCACGGTGGCGACCGACGGCGTCGATTTCGTCGATGAAGATGATGCAAGGCGCATTCTTCTTGGCCTGTTCGAACATGTCACGGACGCGGCTTGc             <  1:202594/144‑1 (MQ=255)
                                                                                                                                        caGCAACTGGTTGAGGGTTTGCTCACGCTCGTCGTGACCGCCGCCCATGCCGGCGCCACGATGGCGACCAACGGCGTCGATTTCGTCGAAGAAGATGATGCACGGGGCGTGCTTCTTGGCCTGCTCGAACATGTCGCGCACACGGCTGGc             >  1:4561798/1‑150 (MQ=255)
                                                                                                                                         agaagcTGGTTCAGCGTCTGTTCGCGTTCGTCGTTGCCGCCGCCAAGACCGGCACCACGGTGGCGACCGACGGCGTCGATTTCGTCGATGAAGATGATGCAAGGCGCATTCTTCTTGGCCTGTTCGAACATGTCACGGACGCGGCTTGCa            >  2:1063628/6‑150 (MQ=255)
                                                                                                                                           cAACTGGTTGAGGGTTTGCTCACGCTCGTCGTGACCGCCGCCCATGCCGGCGCCACGATGGCGACCAACGGCGTCGATTTCGTCGATGAAGATGATGCACGGGGCGTGCTTCTTGGCCTGCTCGAACATGTCGCGCACACGGCTGGCAcc          >  2:1447651/1‑150 (MQ=255)
                                                                                                                                             aCTGGTTGAGGGTTTGCTCACGCTCGTCGTGACCGCCGCCCATGCCGGCGCCACGATGGCGACCAACGGCGTCGATTTCGTCGATGAAGATGATGCACGGGGCGTGCTTCTTGGCCTGCTCGAACATGTCGCGCACACGGCTGGCACCGa        <  1:616508/150‑1 (MQ=255)
                                                                                                                                                   tGAGGGTTTGCTCACGCTCGTCGTGACCGCCGCCCATGCCGGCGCCACGATGGCGACCAACGGCGTCGATTTCGTCGATGAAGATGATGCACGGGGCGTGCTTCTTGGCCTGCTCGAACATGTCGCGCACACGGCTGGCACCGACACCCa  >  1:364114/1‑150 (MQ=255)
                                                                                                                                                    |                                                                                                                                                    
GGATGTCCGGCAGGCCAACCACCACCTGGCGGTCGAAGCGGCCGGGGCGCAGCAGCGCAGGGTCGAGCACATCCGGACGGTTGGTGGCGGCGATGACGATGATGCCATCATTCATCTCGAAGCCGTCCATCTCCACCAGCAACTGGTTGAGGGTTTGCTCACGCTCGTCGTGACCGCCGCCCATGCCGGCGCCACGATGGCGACCAACGGCGTCGATTTCGTCGATGAAGATGATGCACGGGGCGTGCTTCTTGGCCTGCTCGAACATGTCGCGCACACGGCTGGCACCGACACCCA  >  NC_002947/5367538‑5367834

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: