Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 4,888,157 A→T 6.0% intergenic (+69/+54) PP_4295 → / ← PP_4296 TetR family transcriptional regulator/hypothetical protein

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0029474,888,1570AT 6.0% 337.0 / 4.1 109intergenic (+69/+54)PP_4295/PP_4296TetR family transcriptional regulator/hypothetical protein
Reads supporting (aligned to +/- strand):  ref base A (50/50);  new base T (2/4);  total (55/54)
Fisher's exact test for biased strand distribution p-value = 6.79e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

Only 100 of 114 total aligned reads displayed.

CTGGCCGACAGTGATTATGACGCGGCGGCCGAGACCATCATTCGTATGGTCTTGAAAGGTTGTGAGCCTGAGGCAGCCTGAACCATCTCGAAGGGCCATTCGCGGGCAAGCCCGCTCCCACAGGATCACCACATGCCTGAAGGTTGTGGAGTCCCTGTGGGAGCGGGCTTGCCCGCGAATGGGCCAGAACAGGCAATAGAGATTTAGGCCGCTACCCCGGCATCTGCCCGCAGCCCCACCTCCTCGATCGCACTGATCGCACACTGCTCGTCGATATCCGACGTATCGCCGCTGATCCC  >  NC_002947/4888008‑4888306
                                                                                                                                                     |                                                                                                                                                     
cTGGCCGACAGTGATTATGACGCGGCGGCCGAGACCATCATTCGTATGGTCTTGAAAGGTTGTGAGCCTGAGGCAGCCTGAACCATCTCGAAGGGCCATTCGCGGGCAAGCCCGCTCCCACAGGATCACCACATGCCTGAAGGTTGTGGa                                                                                                                                                       <  1:85320/150‑1 (MQ=255)
cTGGCCGACAGTGATTATGACGCGGCGGCCGAGACCATCATTCGTATGGTCTTGAAAGGTTGTGAGCCTGAGGCAGCCTGAACCATCTCGAAGGGCCATTCGCGGGCAAGCCCGCTCCCACAGGATCACCACATGCCTGAAGGTTGTGGa                                                                                                                                                       >  2:3691528/1‑150 (MQ=255)
cTGGCCGACAGTGATTATGACGCGGCGGCCGAGACCATCATTCGTATGGTCTTGAAAGGTTGTGAGCCTGAGGCAGCCTGAACCATCTCGAAGGGCCATTCGCGGGCAAGCCCGCTCCCACAGGATCACCACATGCCTGAAGGTTGTGGa                                                                                                                                                       <  2:1578181/150‑1 (MQ=255)
  ggCCGACAGTGATTATGCCGCGGCGGCCGAGACCATCATTCGTATGGTCTTGAAAGGTTGTGAGCCTGAGGCAGCCTGAACCATCTCGAAGGGCCATTCGCGGGCAAGCCCGCTCCCACAGGATCACCACATGCCTGAAGGTTGTGGAGt                                                                                                                                                     <  1:3936716/150‑1 (MQ=255)
  ggCCGACAGTGATTATGACGCGGCGGCCGAGACCATCATTCGTATGGTCTTGAAAGGTTGTGAGCCTGAGGCAGCCTGAACCATCTCGAAGGGCCATTCGCGGGCAAGCCCGCTCCCACAGGATCACCACATGCCTGAAGGTTGTGGAGt                                                                                                                                                     <  1:3936703/150‑1 (MQ=255)
    cgacgagcagtgtgcgatcagtgcgatcgaggaggtggggctgcgggcagatgccggggtagcggcctaaatctctattgcctgttctggcCCATTCGCGGGCAAGCCCGCTCCCACAGGATCACCACATGCCTGAAGGTTGTGGAGTcc                                                                                                                                                   >  2:3930884‑M1/92‑150 (MQ=255)
    ccGACAGTGATTATGACGCGGCGGCCGAGACCATCATTCGTATGGTCTTGAAAGGTTGTGAGCCTGAGGCAGCCTGAACCATCTCGAAGGGCCATTCGCGGGCAAGCCCGCTCCCCCCCGCTCACCCCAGGCCGGAAGGTTGTGGGGTcc                                                                                                                                                   >  1:137901/1‑150 (MQ=255)
      gACAGTGATTATGACGCGGCGGCCGAGACCATCATTCGTATGGTCTTGAAAGGTTGTGAGCCTGAGGCAGCCTGAACCATCTCGAAGGGCCATTCGCGGGCAATCCCGCTCCCACAGGATCACCACATGCCTGAAGGTTGTGGAGTCCCt                                                                                                                                                 <  2:4181848/150‑1 (MQ=255)
      gACAGTGATTATGACGCGGCGGCCGAGACCATCATTCGTATGGTCTTGAAAGGTTGTGAGCCTGAGGCAGCCTGAACCATCTCGAAGGGCCATTCGCGGGCAAGCCCGCTCCCACAGGATCACCACATGCCTGAAGGTTGTGGAGTCCCt                                                                                                                                                 <  1:3961470/150‑1 (MQ=255)
        gagcagtgtgcgatcagtgcgatcgaggaggtggggctgcgggcagatgccggggtagcggcctaaatctctattgcctgttctggcCCATTCGCGGGCAAGCCCGCTCCCACAGGATCACCACATGCCTGAAGGTTGTGGAGTCCCtgt                                                                                                                                               >  2:2456735‑M1/88‑150 (MQ=255)
        gagcagtgtgcgatcagtgcgatcgaggaggtggggctgcgggcagatgccggggtagcggcctaaatctctattgcctgttctggcCATTCGCGGGCAAGCCCGCTCCCACAGGATCACCACATGCCTGAAGGTTGTGGAGTCCCtgtg                                                                                                                                              >  1:4422771‑M1/88‑150 (MQ=255)
        cAGTGATTATGACGCGGCGGCCGAGACCATCATTCGTATGGTCTTGAAAGGTTGTGAGCCTGAGGCAGCCTGAACCATCTCGAAGGGCCATTCGCGGGCAAGCCCGCTCCCACAGGATCACCACATGCCTGAAGGTTGTGGAGTCCCtgt                                                                                                                                               <  1:1393991/150‑1 (MQ=255)
        cAGTGATTATGACGCGGCGGCCGAGACCATCATTCGTATGGTCTTGAAAGGTTGTGAGCCTGAGGCAGCCTGAACCATCTCGAAGGGCCATTCGCGGGCAAGCCCGCTCCCACAGGATCACCACATGCCTGAAGGTTGTGGAGTCCCtgt                                                                                                                                               <  1:1751736/150‑1 (MQ=255)
        cAGTGATTATGACGCGGCGGCCGAGACCATCATTCGTATGGTCTTGAAAGGTTGTGAGCCTGAGGCAGCCTGAACCATCTCGAAGGGCCATTCGCGGGCAAGCCCGCTCCCACAGGATCACCACATGCCTGAAGGTTGTGGAGTCCCtgt                                                                                                                                               <  2:1967818/150‑1 (MQ=255)
         aGTGATTATGACGCGGCGGCCGAGACCATCATTCGTATGGTCTTGAAAGGTTGTGAGCCTGAGGCAGCCTGAACCATCTCGAAGGGCCATTCGCGGGCAAGCCCGCTCCCACAGGATCACCACATGCCTGAAGGTTGTGGAGTCCCtgtg                                                                                                                                              >  2:3952095/1‑150 (MQ=255)
          gTGATTATGACGCGGCGGCCGAGACCATCATTCGTATGGTCTTGAAAGGTTGTGAGCCTGAGGCAGCCTGAACCATCTCGAAGGGCCATTCGCGGGCAAGCCCGCTCCCACAGGATCACCACATGCCTGAAGGTTGTGGAGTCCCTGTgg                                                                                                                                             <  1:2439413/150‑1 (MQ=255)
          gTGATTATGACGCGGCGGCCGAGACCATCATTCGTATGGTCTTGAAAGGTTGTGAGCCTGAGGCAGCCTGAACCATCTCGAAGGGCCATTCGCGGGAAAGCCCGCTCCCACAGGATCACCACATGCCTGAAGGTTGTGGAGTCCCTGTgg                                                                                                                                             <  2:3016049/150‑1 (MQ=255)
            agtgtgcgatcagtgcgatcgaggaggtggggctgcgggcagatgccggggtagcggcctaaatctctattgcctgttctggcCCATTCGCGGGCAAGCCCGCTCCCACAGGATCACCACATGCCTGAAGGTTGTGGAGTCCCTGTGGGa                                                                                                                                           >  1:2196364‑M1/84‑150 (MQ=255)
              ttATGACGCGGCGGCCGAGACCATCATTCGTATGGTCTTGAAAGGTTGTGAGCCTGAGGCAGCCTGAACCATCTCGAAGGGCCATTCGCGGGCAAGCCCGCTCCCACAGGATCACCACATGCCTGAAGGTTGTGGAGTCCCTGTGGGAGc                                                                                                                                         <  2:3823260/150‑1 (MQ=255)
              ttATGACGCGGCGGCCGAGACCATCATTCGTATGGTCTTGAAAGGTTGTGAGCCTGAGGCAGCCTGAACCATCTCGAAGGGCCATTCGCGGGCAAGCCCGCTCCCACAGGATCACCACATGCCTGAAGGTTGTGGAGTCCCTGTGGGAGc                                                                                                                                         <  1:2395328/150‑1 (MQ=255)
              ttATGACGCGGCGGCCGAGACCATCATTCGTATGGTCTTGAAAGGTTGTGAGCCTGAGGCAGCCTGAACCATCTCGAAGGGCCATTCGCGGGCAAGCCCGCTCCCACAGGATCACCACATGCCTGAAGGTTGTGGAGTCCCTGTGGGAGc                                                                                                                                         <  2:3052631/150‑1 (MQ=255)
              ttATGACGCGGCGGCCGAGACCATCATTCGTATGGTCTTGAAAGGTTGTGAGCCTGAGGCAGCCTGAACCATCTCGAAGGCCCATTCGCGGGCAAGCCCGCTCCCACAGGGACTCCACAACCTTCAGGCATGTGGTGATCCTGTGGGAGc                                                                                                                                         <  2:4205495/150‑1 (MQ=18)
                aTGACGCGGCGGCCGAGACCATCATTCGTATGGTCTTGAAAGGTTGTGAGCCTGAGGCAGCCTGAACCATCTCGAAGGGCCATTCGCGGGCAAGCCCGCTCCCACAGGATCACCACATGCCTGAAGGTTGTGGAGTCCCTGTGGGAGCgg                                                                                                                                       <  1:2657700/150‑1 (MQ=255)
                aTGACGCGGCGGCCGAGACCATCATTCGTATGGTCTTGAAAGGTTGTGAGCCTGAGGCAGCCTGAACCATCTCGAAGGGCCATTCGCGGGCAAGCCCGCTCCCACAGGATCACCACATGCCTGAAGGTTGTGGAGTCCCTGTGGGAGCgg                                                                                                                                       <  1:2657760/150‑1 (MQ=255)
                 tGACGCGGCGGCCGAGACCAGCATTCGTATGGTCTTGAAAGGTTGTGAGCCTGAGGCAGCCTGAACCATCTCGAAGGGCCATTCGCGGGCAAGCCCGCTCCCACAGGATCACCACATGCCTGAAGGTTGTGGAGTCCCTGTGGGAGCggg                                                                                                                                      <  1:1788204/150‑1 (MQ=255)
                  gACGCGGCGGCCGAGACCATCATTCGTATGGTCTTGAAAGGTTGTGAGCCTGAGGCAGCCTGAACCATCTCGAAGGGCCATTCGCGGGCAAGCCCGCTCCCACAGGATCACCACATGCCTGAAGGTTGTGGAGTCCCTGTGGGAGCGGGc                                                                                                                                     >  2:3486941/1‑150 (MQ=255)
                  gACGCGGCGGCCGAGACCATCATTCGTATGGTCTTGAAAGGTTGTGAGCCTGAGGCAGCCTGAACCATCTCGAAGGGCCATTCGCGGGCAAGCCCGCTCCCACAGGATCACCACATGCCTGAAGGTTGTGGAGTCCCTGTGGGAGCGGGc                                                                                                                                     <  1:1041874/150‑1 (MQ=255)
                    cgcgGCGGCCGAGACCATCATTCGTATGGTCTTGAAAGGTTGTGAGCCTGAGGCAGCCTGAACCATCTCGAAGGGCCATTCGCGGGCAAGCCCGCTCCCACAGGGACTCCACAACCTTCAGGCATGTGGTGATCCTGTGGGAGCGGGCtt                                                                                                                                   <  1:3095520/150‑1 (MQ=18)
                    cgcgGCGGCCGAGACCATCATTCGTATGGTCTTGAAAGGTTGTGAGCCTGAGGCAGCCTGAACCATCTCGAAGGGCCATTCGCGGGCAAGCCCGCTCCCACAGGATCACCACATGCCTGAAGGTTGTGGAGTCCCTGTGGGAGCGGGCtt                                                                                                                                   <  1:4446205/150‑1 (MQ=255)
                       ggcggcCGAGACCATCATTCGTATGGTCTTGAAAGGTTGTGAGCCTGAGGCAGCCTGAACCATCTCGAAGGGCCATTCGCGGGCAAGCCCGCTCCCACAGGATCACCACATGCCTGAAGGTTGTGGAGTCCCTGTGGGAGCGGGCTGGcc                                                                                                                                >  2:668238/1‑150 (MQ=255)
                        gcggcCGAGACCATCATTCGTATGGTCTTGAAAGGTTGTGAGCCTGAGGCAGCCTGAACCATCTCGAAGGGCCATTCGCGGGCAAGCCCGCGCCCACAGGATCACCACATGCCTGAAGGTTGTGGAGTCCCTGTGGGAGCGGGCTTGccc                                                                                                                               >  1:4283188/1‑150 (MQ=255)
                          ggcCGAGACCATCATTCGTATGGTCTTGAAAGGTTGTGAGCCTGAGGCAGCCTGAACCATCTCGAAGGGCCATTCGCGGGCAAGCCCGCTCCCACAGGATCACCACATGCCTGAAGGTTGTGGAGTCCCTGTGGGAGCGGGCTTGCCcgc                                                                                                                             <  2:4057087/150‑1 (MQ=255)
                                   catcatTCGTATGGTCTTGAAAGGTTGTGAGCCTGAGGCAGCCTGAACCATCTCGAAGGGCCATTCGCGGGCAAGCCCGCTCCCACAGGATCACCACATGCCAGAAGGTTGTGGAGTCCCTGTGGGAGCCGGCATGACCGCGAAAGGGcc                                                                                                                    >  2:563576/1‑150 (MQ=39)
                                    atcatTCGTATGGTCTTGAAAGGTTGTGAGCCTGTGGCAGCCTGAACCATCTCGAAGGGCCATTCGCGGGCAAGCCCGCTCCCACAGGATCACCACATGCCTGATGGTTGTGGAGTCCCTGTGGGAGGGGGCTTGCCCGCGAATGGGCCa                                                                                                                   <  2:533688/150‑1 (MQ=255)
                                                   ttGAAAGGTTGTGAGCCTGAGGCAGCCTGAACCATCTCGAAGGGCCATTCGCGGGCAAGCCCGCTCCCACACGCTCACCACAAGCCTGAAGGTTGTGGTGTCCCTGTGGGAGCGGGCTTGCCCGCGAATGGGCCAGAACACGCACTagag                                                                                                    >  1:3115711/1‑150 (MQ=38)
                                                   ttGAAAGGTTGTGAGCCTGAGGCAGCCTGAACCATCTCGAAGGGCCATTCGCGGGCAAGCCCGCGCCCACACGCTAACCACATGCCTGAAGGTTGTGGAGTCCCAGTGGGAGCGGGCTTGCCCGCGAAAGGGCCAGAACAGGCCAAagag                                                                                                    >  1:950007/1‑150 (MQ=255)
                                                     gAAAGGTTGTGAGCCTGAGGCAGCCTGAACCATCTCGAAGGGCCATTCGCGGGCAAGCCCGCTCCCACAGGATCACCACATGCCTGAAGGTTGTGGAGTCCCtgtg                                                                                                                                              >  2:518357/1‑106 (MQ=255)
                                                     gAAAGGTTGTGAGCCTGAGGCAGCCTGAACCATCTCGAAGGGCCATTCGCGGGCAAGCCCGCTCCCACAGGATCACCACATGCCTGAAGGTTGTGGAGTCCCtgtg                                                                                                                                              <  1:518357/106‑1 (MQ=255)
                                                     gAAAGGTTGTGAGCCTGAGGCAGCCTGAACCATCTCGAAGGGCCATTCGCGGGCAAGCCCGCTCCCACAGGATCACCACAGGCCTGAAGGGTGTGGAGTCCCTGTGGGAGCGGGCTTGCCCGCGAATGGGCCAGAACAGGCAAGAGAGAt                                                                                                  >  2:1076718/1‑150 (MQ=255)
                                                      aaaGGTTGTGAGCCTGAGGCAGCCTGAACCATCTCGAAGGGCCATTCGCGGGCAAGCCCGCTCCCACAGGATCACCACATGCCTGAAGGTTGTGGAGTCCCTGTGGGAGCGGGCTTGCCCGCGAATGGGCCAGAACAGGCAATAGAGAtt                                                                                                 >  2:1147613/1‑150 (MQ=255)
                                                        aGGTTGTGAGCCTGAGGCAGCCTGAACCATCTCGACGGGCCATTCGCGGGCAAGCCCGCTCCCACAGGATCACCACTTGCCTTAAGTTTGTGGAGTCCCTGTGGGAGCGGGCTTGCCCGCGAATGGGCCAGAACAGGCAATAGAGATTTa                                                                                               <  2:1633814/150‑1 (MQ=255)
                                                         ggTTGTGAGCCTGAGGCAGCCTGAACCATCTCGAAGGGCCATTCGCGGGCAAGCCCGCTCCCACAGGATCACCACATGCCTGAAGGTTGTGGAGTCCCTGTGGGAGCGGGCTTGCCCCCGAATGGGCCAGAACACGCAAAAGAGATTTAg                                                                                              >  1:3124323/1‑150 (MQ=255)
                                                         ggTTGTGAGCCTGAGGCAGCCTGAACCATCTCGAAGGGCCATTCGCGGGCAAGCCCGCGCCCACAGGATCACCACCTGCCGGAAGGGTGTGGTGTGCCAGTGGGAGCGGGCTTGCCCGCGAATGGGCCAGAACAGGCAATAGAGATTTat                                                                                              >  1:3124342/1‑149 (MQ=38)
                                                          gTTGTGAGCCTGAGGCAGCCTGAACCACCACGAAGGGCCATCCGCGGGCAAGCCCGCTCCCACAGGATCACCACATCCCTCAAGTTTGTGGAGACCCTGTGGGAGCGGGCTTGCCCGCGAATGGGCCAGAACAGGCAATAGAGATTTAgg                                                                                             <  2:907122/150‑1 (MQ=37)
                                                             gtgAGCCTGAGGCAGCCTGAACCATCTCGAAGGGCCATTCGCGGGCAAGCCCGCTCCCACAGGATCACCACATGCCTGAAGGTTGTGGAGTCCCTGTGGGAGCGGGCTTGCCCGCGAATGGGCCAGAACAGGCAATAGAGATTTAGGCCg                                                                                          <  2:355711/150‑1 (MQ=255)
                                                                   ctgtgccAGCCTGAACCATCTCGACGGGCCATTCGCGGGCACGCCCGCTCCCACAGTAACACAACTTGCCTGATGGTTTGGGAGTCCCTGTGGGCGGGGGCTTGCCCGCGAATGGGCCAGAACAGGCAATAGAGATTTAGGCCGCTAccc                                                                                    <  1:3920253/144‑1 (MQ=255)
                                                                    tGAGGCAGCCTGAACCATCTCGAAGGGCCATTCGCGGGCAAGCCCGCTCCCACAGGATCACCACATGCCTGAAGGTTGTGGAGTCCCTGTGGGAGCGGGCTTGCCCGCGAATGGGCCAGAACAGGCAATAGAGCTTTAGGCCGCTAcccc                                                                                   >  2:3880151/1‑150 (MQ=255)
                                                                            ccTGAACCATCTCGAAGGGCCATTCGCGGGCAAGCCCGCTCCCACAGGATCACCACATGCCTGAAGGTTGTGGAGTCCCTGTGGGAGCGGGCTTGCCCGCGAATGGGCCAGAACAGGCAATAGAGATTTAGGCCGCTACCCCGGCATCTg                                                                           >  2:1376974/1‑150 (MQ=255)
                                                                               gAACCATCTCGAAGGGCCATTCGCGGGCAAGCCCGCTCCCACAGGATCACCACATGCCTGAAGGTTGTGGAGTCCCTGTGGGAGCGGGCTTGCCCGCGAATGGGCCAGAACAGGCAATAGAGATTTAGGCCGCTACCCCGGCATCTGccc                                                                        <  2:760255/150‑1 (MQ=255)
                                                                                     tctcGAAGGGCCATTCGCGGGCAAGCCCGCTCCCACAGGATCACCACATGCCTGAAGGTTGTGGAGTCCCTGTGGGAGCGGGCTTGCCCGCGAATGGGCCAGAACAGGCAATAGAGATTTAGGCCGCTACCCCGGCATCTGCCCGCAGcc                                                                  >  2:2709441/1‑150 (MQ=255)
                                                                                         gAAGGGCCATTCGCGGGCAAGCCCGCTCCCACAGGATCACCACATGCCAGAAGGTTGTGGAGTCCCGGTGGGAGCGGGCTTGCCCGCGAATGGGCCACAACAGGCAATAGAGATTTAGGCCGCTACCCCCGCATCTGCCCGCACCCCCAc                                                              >  1:2969922/1‑150 (MQ=39)
                                                                                               CCATTCGCGGGCACGCCCGCTCCCACAGCATCACCCACTGCCTGAAGGTGGGGGAGTCCCTGTGGGAGCGGGCGTGCCCGCGAATGGGCCAGAACAGGCAATAGAGATTTAGGCCGCTACCCCGGCATCTGCCCGCAGCCCCACctcctc                                                        >  1:3214610‑M1/1‑150 (MQ=255)
                                                                                               CCATTCGCGGGCAAGCCCGCTCCCACAGGATCACCACATGCCTGAAGGTTGTGGAGTCCCTGTGGGAGCGGGCTTGCCCGCGAATGGGCCAGAACAGGCAATAGAGATTTAGGCCGCTACCCCGGCATCTGCCCGCAGCCCCACctcctc                                                        <  2:3245180‑M1/150‑1 (MQ=255)
                                                                                               CCATTCGCGGGCAAGCCCGCTCCCACAGGATCACCACATGCCTGAAGGTTGTGGAGTCCCTGTGGGAGCGGGCTTGCCCGCGAATGGGCCAGAACAGGCAATAGAGATTTAGGCCGCTACCCCGGCATCTGCCCGCAGCCCCACctcctc                                                        <  2:2786599‑M1/150‑1 (MQ=255)
                                                                                               CCATTCGCGGGCAAGCCCGCTCCCACAGGATCACCACATGCCTGAAGGTTGTGGAGTCCCTGTGGGAGCGGGCTTGCCCGCGAATGGGCCAGAACAGGCAATAGAGATTTAGGCCGCTACCCCGGCATCTGCCCGCAGCCCCACctcctc                                                        >  2:2230453‑M1/1‑150 (MQ=255)
                                                                                                CATTCGCGGGCAAGCCCGCTCCCACAGGATCACCACATGCCTGAAGGTTGTGGAGTACCTGTGGGAGCGGGCTTGCCCGCGAATGGGCCACAACAGGCAATAGAGATTTAGGCCGCTAACCCCGCATCTGCCCGCAGCCCCACCTCCTCg                                                       <  1:718399‑M1/150‑1 (MQ=255)
                                                                                                  TTCGCGGGCAAGCCCGCTCCCACAGGATCACCACATGCCTGAAGGTTGTGGAGTCCCTGTGGGAGCGGGCTTGCCCGCGAATGGGCCAGAACAGGCAATAGAGATTTAGGCCGCTACCCCGGCATCTGCCCGCAGCCCCACCTCCTCGAt                                                     >  2:3507237‑M1/1‑150 (MQ=255)
                                                                                                      CGGGCAAGCCCGCTCCCACAGTATCACCACTTGCCTGATGGTTGTGGAGTCCCTGTGGGCGCGGGCTTGCCCGCGAATGGGCCAGAACAGGCAATAGAGATTTAGGCCGCTACCCCGGCATCTGCCCGCAGCCCCACCTCCTCGATCGCa                                                 >  2:137901‑M1/1‑150 (MQ=255)
                                                                                                      CGGGCAAGCCCGCTCCCACAGGATCACCACATGCCTGAAGGTTGTGGAGTCCCTGTGGGCGCGGGCTTGCCCGCGAATGGGCCAGAACTGGCAATAGAGATTTAGGCCGCTACCCCGGCATCTGCCCGCAGCCCCACCTCCTCGATCGCa                                                 >  1:3691528‑M1/1‑150 (MQ=255)
                                                                                                       GGGCAAGCCCGCTCCCACAGGATCACCACATGCCTGAAGGTTGTGGAGTCCCTGTGGGAGCGGGCTTGCCCGCGAATGGGCCAGAACAGGCAATAGAGATTTAGGCCGCTACCCCGGCATCTGCCCGCAGCCCCACCTCCTCGATCGCAc                                                <  2:188336‑M1/150‑1 (MQ=255)
                                                                                                         GCAAGCCCGCTCCCACAGGATCACCACATGCCTGAAGGTTGTGGAGTCCCTGTGGGAGCGGGCTTGCCCGCGAATGGGCCAGAACAGGCAATAGAGAGTTAGGCCGCTACCCCGGCATCTGCCCGCAGCCCCACCTCCTCGATCGCACTg                                              <  2:3502475‑M1/150‑1 (MQ=255)
                                                                                                           AAGCCCGCTCCCACAGGATCACCACATGCCTGAAGGTTGTGGAGTCCCTGTGGGAGCGGGCTTGCCCGCGAATGGGCCAGAACAGGCAATAGAGATTTAGGCCGCTACCCCGGCATCTGCCCGCAGCCCCACCTCCTCGATCGCACTGAt                                            <  2:2275023‑M1/150‑1 (MQ=255)
                                                                                                           AAGCCCGCTCCCACAGGATCACCACATGCCTGAAGGTTGTGGAGTCCCTGTGGGAGCGGGCTGGCCCGCGAATGGGCCAGAACAGGCAATAGAGATTTAGGCCGCTACCCCGGCATCTGCCCGCAGCCCCACCTCCTCGATCGCACTGAt                                            >  1:668238‑M1/1‑150 (MQ=255)
                                                                                                               CCGCTCCCACAGGGACTCCACAACCTTCAGGCATGTGGTGATCCTGTGGGAGCGGGCTTGCCCGCGAATGGGCCAGAACAGGCAATAGAGATTTAGGCCGCTACCCCGGCATCTGCCCGCAGCCCCACCTCCTCGATCGCACTGATCGca                                        <  2:885426‑M1/150‑1 (MQ=255)
                                                                                                               CCGCTCCCACAGGATCACCACATGCCTGAAGGTTGTGGAGTCCCTGTGGGAGCGGGCTTGCCCGCGAATGGGCCAGAACAGGCAATAGAGATTTAGGCCGCTACCCCGGCATCTGCCCGCAGCCCCACCTCCTCGATCGCACTGATCGca                                        <  1:914258‑M1/150‑1 (MQ=255)
                                                                                                                 GCTCCCACTGGATCACCACATGCCTGACGTTTGTGTAGGCCCGGGGGGCGCGGGCTTGCCCGCGAATGGGCCAGAACAGGCAATAGAGATTTAGGCCGCTACCCCGGCATCTGCCCGCAGCCCCACCTCCTCGATCGCACTGATCGcaca                                      >  1:2257390‑M1/1‑150 (MQ=255)
                                                                                                                      CACAGGATCACCACATGCCTGAAGGTTGTGGAGTCCCTGTGGGAGCGGGCTTGCCCGCGAATGGGCCAGAACAGGCAATAGAGATTTAGGCCGCTACCCCGGCATCTGCCCGCAGCCCCACCTCCTCGATCGCACTGATCGCACACTGCt                                 <  2:872960‑M1/150‑1 (MQ=255)
                                                                                                                      CACAGGATCACCACATGCCTGAAGGTTGTGGAGTCCCTGTGGGAGCGGGCTTGCCCGCGAATGGGCCAGAACAGGCAATAGAGATTTAGGCCGCTACCCCGGCATCTGCCCGCAGCCCCACCTCCTCGATCGCACTGATCGCACACTGCt                                 <  1:1008876‑M1/150‑1 (MQ=255)
                                                                                                                       ACAGGATCACCACATGCCTGAAGGTTGTGGAGTCCCTGTGGGAGCGGGCTTGCCCGCGAATGGGCCAGAACAGGCAATAGAGATTTAGGCCGCTACCCCGGCATCTGCCCGCAGCCCCACCTCCTCGATCGCACTGATCGCACACTGCtc                                <  2:3436981‑M1/150‑1 (MQ=255)
                                                                                                                         AGGATCACCACATGCCTGAAGGTTGTGGAGTCCCTGTGGGAGCGGGCTTGCCCGCGAATGGGCCAGAACAGGCAATAGAGATTTAGGCCGCTACCCCGGCATCTGCCCGCAGCCCCACCTCCTCGATCGCACTGATCGCACACTGCtcgt                              <  1:1838717‑M1/150‑1 (MQ=255)
                                                                                                                            aTCACCACATGCCTGAAGGTTGTGGAGTCCCTGTGGGAGCGGGCTTGCCCGCGAATGGGCCAGAACAGGCAATAGAGATTTAGGCCGCTACCCCGGCATCTGCCCGCAGCCCCACCTCCTCGATCGCACTGATCGCACACTGCTCGTCGa                           >  1:4326187/1‑150 (MQ=255)
                                                                                                                            aTCACCACATGCCTGAAGGTTGTGGAGTCCCTGTGGGAGCGGGCTTGCCCGCGAATGGGCCAGAACAGGCAATAGAGATTTAGGCCGCTACCCCGGCATCTGCCCGCAGCCCCACCTCCTCGATCGCACTGATCGCACACTGCTCGTCGa                           >  1:4326173/1‑150 (MQ=255)
                                                                                                                            aTCACCACATGCCTGAAGGTTGTGGAGTCCCTGTGGGAGCGGGCTTGCCCGCGAATGGGCCAGAACAGGCAATAGAGATTTAGGCCGCTACCCCGGCATCTGCCCGCAGCCCCACCTCCTCGATCGCACTGATCGCACACTGCTCGTCGa                           >  1:1525154/1‑150 (MQ=255)
                                                                                                                             tCACCACATGCCTGAAGGTTGTGGAGTCCCTGTGGGAGCGGGCTTGCCCGCGAATGGGCCAGAACAGGCAATAGAGATTTAGGCCGCTACCCCGGCATCTGCCCGCAGCCCCACCTCCTCGATCGCACTGATCGCACACTGCTCGTCGat                          >  1:1999981/1‑150 (MQ=255)
                                                                                                                                ccacATGCCTGAAGGTTGTGGAGTCCCTGTGGGAGCGGGCTTGCCCGCGAATGGGCCAGAACAGGCAATAGAGATTTAGGCCGCTACCCCGGCATCTGCCCGCAGCCCCACCTCCTCGATCGCACTGATCGCACACTGCTCGTCGATATc                       >  2:1400701/1‑150 (MQ=255)
                                                                                                                                ccacATGCCTGAAGGTTGTGGAGTCCCTGTGGGAGCGGGCTTGCCCGCGAATGGGCCAGAACAGGCAATAGAGATTTAGGCCGCTACCCCGGCATCTGCCCGCAGCCCCACCTCCTCGATCGCACTGATCGCACACTGCTCGTCGATATc                       >  2:3128609/1‑150 (MQ=255)
                                                                                                                                 cacaTGCCTGAAGGTTGTGGAGTCCCTGTGGGAGCGGGCTTGCCCGCGAATGGGCCAGAACAGGCAATAGAGATTTAGGCCGCTACCCCGGCATCTGCCCGCAGCCCCACCTCCTCGATCGCACTGATCGCACACTGCTCGTCGATATcc                      >  1:1254962/1‑150 (MQ=255)
                                                                                                                                        cTGAAGGTTGTGGAGTCCCTGTGGGAGCGGGCTTGCCCGCGAATGGGCCAGAACAGGCAATAGAGATTTAGGCCGCTACCCCGGCATCTGCCCGCAGCCCCACCTCCTCGATCGCACTGATCGCACACTGCTCGTCGATATCCGACGTAt               >  1:1119774/1‑150 (MQ=255)
                                                                                                                                         tGAAGGTTGTGGAGTCCCTGTGGGAGCGGGCTTGCCCGCGAATGGGCCAGAACAGGCAATAGAGATTTAGGCCGCTACCCCGGCATCTGCCCGCAGCCCCACCTCCTCGATCGCACTGATCGCACACTGCTCGTCGATATCCGACGTATc              >  2:2950420/1‑150 (MQ=255)
                                                                                                                                         tGAAGGTTGTGGAGTCCCTGTGGGAGCGGGCTTGCCCGCGAATGGGCCAGAACAGGCAATAGAGATTTAGGCCGCTACCCCGGCATCTGCCCGCAGCCCCACCTCCTCGATCGCACTGATCGCACACTGCTCGTCGATATCCGACGTATc              >  2:1502057/1‑150 (MQ=255)
                                                                                                                                         tGAAGGTTGTGGAGTCCCTGTGGGAGCGGGCTTGCCCGCGAATGGGCCAGAACAGGCAATAGAGATTTAGGCCGCTACCCCGGCATCTGCCCGCAGCCCCACCTCCTCGATCGCACTGATCGCACACTGCTCGTCGATATCCGACGTATc              >  1:3865652/1‑150 (MQ=255)
                                                                                                                                         tGAAGGTTGTGGAGTCCCTGTGGGAGCGGGCTTGCCCGCGAATGGGCCAGAACAGGCAATAGAGATTTAGGCCGCTACCCCGGCATCTGCCCGCAGCCCCACCTCCTCGATCGCACTGATCGCACACTGCTCGTCGATATCCGACGTATc              >  2:4411478/1‑150 (MQ=255)
                                                                                                                                          gAAGGTGGTGGAGTCCCTGTGGGAGCGGGCTTGCCCGCGAATGGGCCAGAACAGGCAATAGAGATTTAGGCCGCTACCCCGGCATCTGCCCGCAGCCCCACCTCCTCGATCGCACTGATCGCACACTGCTCGTCGATATCCGACGTATcg             <  1:2950179/150‑1 (MQ=255)
                                                                                                                                           aaGGTTGTGGAGTCCCTGTGGGAGCGGGCTTGCCCGCGAATGGGCCAGAACAGGCAATAGAGATTTAGGCCGCTACCCCGGCATCTGCCCGCAGCCCCACCTCCTCGATCGCACTGATCGCACACTGCTCGTCGATATCCGACGTATcgc            >  1:1305795/1‑150 (MQ=255)
                                                                                                                                           aaGGTTGTGGAGTCCCTGTGGGAGCGGGCTTGCCCGCGAATGGGCCAGAACAGGCAATAGAGATTTAGGCCGCTACCCCGGCATCTGCCCGCAGCCCCACCTCCTCGATCGCACTGATCGCACACTGCTCGTCGATATCCGACGTATcgc            <  2:4283196/150‑1 (MQ=255)
                                                                                                                                           aaGGTTGTGGAGTCCCTGTGGGAGCGGGCTTGCCCGCGAATGGGCCAGAACAGGCAATAGAGATTTAGGCCGCTACCCCGGCATCTGCCCGCAGCCCCACCTCCTCGATCGCACTGATCGCACACTGCTCGTCGATATCCGACGTATcgc            <  2:4283188/150‑1 (MQ=255)
                                                                                                                                              gTTGTGGAGTCCCTGTGGGAGCGGGCTTGCCCGCGAATGGGCCAGAACAGGCAATAGAGATTTAGGCCGCTACCCCGGCATCTGCCCGCAGCCCCACCTCCTCGATCGCACTGATCGCACACTGCTCGTCGATATCCGACGTATcgccgc         <  1:1754347/150‑1 (MQ=255)
                                                                                                                                              gTTGTGGAGTCCCTGTGGGAGCGGGCTTGCCCGCGAATGGGCCAGAACAGGCAATAGAGATTTAGGCCGCTACCCCGGCATCTGCCCGCAGCCCCACCTCCTCGATCGCACTGATCGCACACTGCTCGTCGATATCCGACGTATcgccgc         <  2:3849034/150‑1 (MQ=255)
                                                                                                                                              gTTGTGGAGTCCCTGTGGGAGCGGGCTTGCCCGCGAATGGGCCAGAACAGGCAATAGAGATTTAGGCCGCTACCCCGGCATCTGCCCGCAGCCCCACCTCCTCGATCGCACTGATCGCACACTGCTCGTCGATATCCGACGTATcgccgc         <  2:1999981/150‑1 (MQ=255)
                                                                                                                                               ttGTGGAGTCCCTGTGGGAGCGGGCTTGCCCGCGAATGGGCCAGAACAGGCAATAGAGATTTAGGCCGCTACCCCGGCATCTGCCCGCAGCCCCACCTCCTCGATCGCACTGATCGCACACTGCTCGTCGATATCCGAc                   <  2:2340052/139‑1 (MQ=255)
                                                                                                                                               ttGTGGAGTCCCTGTGGGAGCGGGCTTGCCCGCGAATGGGCCAGAACAGGCAATAGAGATTTAGGCCGCTACCCCGGCATCTGCCCGCAGCCCCACCTCCTCGATCGCACTGATCGCACACTGCTCGTCGATATCCGAc                   >  1:2340052/1‑139 (MQ=255)
                                                                                                                                               ttGTGGAGTCCCTGTGGGAGCGGGCTTGCCCGCGAATGGGCCAGAACAGGCAATAGAGATTTAGGCCGCTACCCCGGCATCTGCCCGCAGCCCCACCTCCTCGATCGCACTGATCGCACACTGCTCGTCGATATCCGACGTATCGCCGCt        >  2:4114062/1‑150 (MQ=255)
                                                                                                                                               ttGTGGAGTCCCTGTGGGAGCGGGCTTGCCCGCGAATGGGCCAGAACAGGCAATAGAGATTTAGGCCGCTACCCCGGCATCTGCCCGCAGCCCCACCTCCTCGATCGCACTGATCGCACACTGCTCGTCGATATCCGACGTATCGCCGCt        >  2:335129/1‑150 (MQ=255)
                                                                                                                                                tgtgGAGTCCCTGTGGGAGCGGGCTTGCCCGCGAATGGGCCAGAACAGGCAATAGAGATTTAGGCCGCTACCCCGGCATCTGCCCGCAGCCCCACCTCCTCGATCGCACTGATCGCACACTGCTCGTCGATATCCGACGTATCGCCGCTg       >  2:3569515/1‑150 (MQ=255)
                                                                                                                                                  tgGAGTCCCTGTGGGAGCGGGCTTGCCCGCGAATGGGCCAGAACAGGCAATAGAGATTTAGGCCGCTACCCCGGCATCTGCCCGCAGCCCCACCTCCTCGATCGCACTGATCGCACACTGCTCGTCGATATCCGACGTATCGCCGCTGAt     <  2:608737/150‑1 (MQ=255)
                                                                                                                                                  tgGAGTCCCTGTGGGAGCGGGCTTGCCCGCGAATGGGCCAGAACAGGCAATAGAGATTTAGGCCGCTACCCCGGCATCTGCCCGCAGCCCCACCTCCTCGATCGCACTGATCGCACACTGCTCGTCGATATCCGACGTATCGCCGCTGAt     <  2:2059186/150‑1 (MQ=255)
                                                                                                                                                   ggAGTCCCTGTGGGAGCGGGCTTGCCCGCGAATGGGCCAGAACAGGCAATAGAGATTTAGGCCGCTACCCCGGCATCTGCCCGCAGCCCCACCTCCTCGATCGCACTGATCGCACACTGCTCGTCGATATCCGACGTATCGCCGCTGATc    >  2:3366720/1‑150 (MQ=255)
                                                                                                                                                    gAGTCCCTGTGGGAGCGGGCTTGCCCGCGAATGGGCCAGAACAGGCAATAGAGATTTAGGCCGCTACCCCGGCATCTGCCCGCAGCCCCACCTCCTCGATCGCACTGATCGCACACTGCTCGTCGATATCCGACGTATCGCCGCTGATcc   >  2:1623475/1‑150 (MQ=255)
                                                                                                                                                     aGTCCCTGTGGGAGCGGGCTTGCCCGCGACTGGGCCAGAACAGGCAATAGAGATTTAGGCCGCTACCCCGGCATCTGCCCGCAGCCCCACCTCCTCGATCGCACTGATCGCACACTGCTCGTCGATATCCGACGTATCGCCGCTGATccc  >  1:2417874/1‑150 (MQ=255)
                                                                                                                                                     aGTCCCTGTGGGAGCGGGCTTGCCCGCGAATGGGCCAGAACAGGCAATAGAGATTTAGGCCGCTACCCCGGCATCTGCCCGCAGCCCCACCTCCTCGATCGCACTGATCGCACACTGCTCGTCGATATCCGACGTATCGCCGCTGATccc  >  1:2417791/1‑150 (MQ=255)
                                                                                                                                                     |                                                                                                                                                     
CTGGCCGACAGTGATTATGACGCGGCGGCCGAGACCATCATTCGTATGGTCTTGAAAGGTTGTGAGCCTGAGGCAGCCTGAACCATCTCGAAGGGCCATTCGCGGGCAAGCCCGCTCCCACAGGATCACCACATGCCTGAAGGTTGTGGAGTCCCTGTGGGAGCGGGCTTGCCCGCGAATGGGCCAGAACAGGCAATAGAGATTTAGGCCGCTACCCCGGCATCTGCCCGCAGCCCCACCTCCTCGATCGCACTGATCGCACACTGCTCGTCGATATCCGACGTATCGCCGCTGATCCC  >  NC_002947/4888008‑4888306

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