Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 4,131,155 T→C 7.2% intergenic (+23/+89) PP_3634 → / ← PP_3635 GNAT family acetyltransferase/sulfonate ABC transporter permease

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0029474,131,1550TC 7.2% 423.4 / 11.3 149intergenic (+23/+89)PP_3634/PP_3635GNAT family acetyltransferase/sulfonate ABC transporter permease
Reads supporting (aligned to +/- strand):  ref base T (62/70);  new base C (3/7);  total (68/81)
Fisher's exact test for biased strand distribution p-value = 3.44e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.88e-01

Only 100 of 150 total aligned reads displayed.

CACCTATCTGCTTAACCAGCAGGGTGAGGATGCAATGATCGAAGCGGAGCAGATGCGCGATGAGCGGCTGGATGAGCTGACCGCTCAAGCCCAGGAATGTGTGCGCCAGCGCAAGAACGTGCATTGATCCGATGGCCTGTTCGCGGGCTTGCCCGCGCCCACAGGAATACCGTAAGGCCCAAGGCCGGTGATGTACCTGTGGGAGCGGGCGAGCCCGCGAAGAAGCCCCATTGGCCTTTCAGCCATGCACCCAGCGCCGATGCAGCCCACGCGCCAGAGCATCCAGCACAAACCCCA  >  NC_002947/4131006‑4131302
                                                                                                                                                     |                                                                                                                                                   
 acCTATCTGCTTAACCAGCAGGGTGAGGATGCAATGATCGAAGCGGAGCAGATGCGCGATGAGCGGCTGGATGAGCTGACCGCTCAAGCCCAGGAATGTGTGCGCCAGCGCAAGAACGTGCATTGATCCGATGGCCTGTTCGCGGGCTTg                                                                                                                                                    >  2:370054/1‑150 (MQ=255)
  ccTATCTGCTTAACCAGCAGGGTGAGGATGCAATGATCGAAGCGGAGCAGATGCGCGATGAGCGGCTGGATGAGCTGACCGCTCAAGCCCAGGAATGTGTGCGCCAGCGCAAGAACGTGCATTGATCCGATGGCCTGTTCGCGGGCTTGc                                                                                                                                                   <  1:2296903/150‑1 (MQ=255)
   cTATCTGCTTAACCAGCAGGGTGAGGATGCAATGATCGAAGCGGAGCAGATGCGCGATGAGCGGCTGGATGAGCTGACCGCTCAAGCCCAGGAATGTGTGCGCCAGCGCAAGAACGTGCATTGATCCGATGGCTTCTTCGCGGGCTCGcc                                                                                                                                                  >  2:760055/1‑150 (MQ=255)
   cTATCTGCTTAACCAGCAGGGTGAGGATGCAATGATCGAAGCGGAGCAGATGCGCGATGAGCGGCTGGATGAGCTGACCGCTCAAGCCCAGGAATGTGTGCGCCAGCGCAAGAACGTGCATTGATCCGATGGCCTGTTCGCGGGCTTGcc                                                                                                                                                  >  1:370413/1‑150 (MQ=255)
    tATCTGATTAACCAGCAGGGTGAGGATGCAATGATCGAAGCGGAGCAGATGCGCGATGAGCGGCTGGATGAGCTGACCGCTCAAGCCCAGGAATGAGTGCGCCAGCGCAAGAACGTGCATTGATCCGATGGCCTGTTCGCGGGCTTGccc                                                                                                                                                 >  2:1643885/1‑150 (MQ=255)
     aTCTGCTTAACCAGCAGGGTGAGGATGCAATGATCGAAGCGGAGCAGATGCGCGATGAGCGGCTGGATGAGCTGACCGCTCAAGCCCAGGAATGTGTGCGCCAGCGCAAGAACGTGCATTGATCCGATGGCCTGTTCGCGGGCTTGCCcg                                                                                                                                                >  1:607717/1‑150 (MQ=255)
      tCTGCTTAACCAGCAGGGTGAGGATGCAATGATCGAAGCGGAGCAGATGCGCGATGAGCGGCTGGATGAGCTGACCGCTCAAGCCCAGGAATGTGTGCGCCAGCGCAAGAACGTGCATTGATCCGATGGCCTGTTCGCGGGCTTGCCcgc                                                                                                                                               <  1:1051503/150‑1 (MQ=255)
        tGCTTAACCAGCAGGGTGAGGATGCAATGATCGAAGCGGAGCAGATGCGCGATGAGCGGCTGGATGAGCTGACCGCTCAAGCCCAGGAATGTGTGCGCCAGCGCAAGAACGTGCATTGATCCGATGGCCTGTTCGCGGGCTTGCCCGCGC                                                                                                                                             >  2:203579‑M1/1‑150 (MQ=255)
         gcataacCAGCAGGGCGAGGATGCAATGATCGAAGCGGAGCAGATGCGCGATGAGCGGCTGGATGAGCTGACGGCTCAAGCCCAGGAATGTGTGCGCCAACCCAAGAACGTGCATTGATCCGAAGGCCTGTTCTCGAGCTTGCCCGCGCC                                                                                                                                            >  2:1739081‑M1/4‑150 (MQ=255)
          cTTAACCAGCAGGGTGAGGATGCAATGATCGAAGCGGAGCAGATGCGCGATGAGCGGCTGGATGAGCTGACCGCTCAAGCCCAGGAATGTGTGCGCCAGCGCAAGAACGTGCATTGATCCGATGGCCTGTTCGCGGGCTTGCCCGCGCCC                                                                                                                                           >  2:1792635‑M1/1‑150 (MQ=255)
             aaCCAGCAGGGTGAGGATGCAATGATCGAAGCGGAGCAGATGCGCGATGAGCGGCTGGATGAGCTGACCGCTCAAGCCCAGGAATGTGTGCGCCAGCGCAAGAACGTGCATTGATCCGATGGCCTGTTCGCGGGCTTGCCCGCGCCCACA                                                                                                                                        <  1:2346876‑M1/150‑1 (MQ=255)
              aCCAGCAGGGTGAGGATGCAATGATCGAAGCGGAGCAGATGCGCGATGAGCGGCTGGATGAGCTGACCGCTCAAGCCCAGGAATGTGTGCGCCAGCGCAAGAACGTGCATTGATCCGATGGCCTGTTCGCGGGCTTGCCCGCGCCCACAG                                                                                                                                       <  2:243927‑M1/150‑1 (MQ=255)
               ccAGCAGGGTGAGGATGCAATGATCGAAGCGGAGCAGATGCGCGATGAGCGGCTGGATGAGCTGACCGCTCAAGCCCAGGAATGTGTGCGCCAGCGCAAGAACGTGCATTGATCCGATGGCCTGTTCGCGGGCTTGCCCGCGCCCACAGG                                                                                                                                      >  2:1882425‑M1/1‑150 (MQ=255)
               ccAGCAGGGTGAGGATGCAATGATCGAAGCGGAGCAGATGCGCGATGAGCGGCTGGATGAGCTGACCGCTCAAGCCCAGGAATGTGTGCGCCAGCGCAAGAACGTGCATTGATCCGATGGCCTGTTCGCGGGCTTGCCCGCGCCCACAGG                                                                                                                                      >  2:1415451‑M1/1‑150 (MQ=255)
                cagcagGGTGAGGATGCAATGATCGAAGCGGAGCAGATGCGCGATGAGCGGCTGGATGAGCTGACCGCTCAAGCCCAGGAATGTGTGCGCCAGCGCAAGAACGTGCATTGATCCGATGGCCTGTTCGCGGGCTTGCCCGCGCCCACAGGt                                                                                                                                     >  1:1647288‑M1/1‑149 (MQ=255)
                cagcagGGTGAGGATGCAATGATCGAAGCGGAGCAGATGCGCGATGAGCGGCTGGATGAGCTGACCGCTCAAGCCCAGGAATGTGTGCGCCAGCGCAAGAACGTGCATTGATCCGATGGCCTGTTCGCGGGCTTGCCCGCGCCCACAGGa                                                                                                                                     >  1:987844/1‑150 (MQ=255)
                 aatcgcgTGAGGATGCAATGATCGAAGCGGAGCAGATGCGCGATGAGCGGCTGGATGAGCTGACCGCTCAAGCCCAGGAATGTGTGCGCCAGCGCAAGAACGTGCATTGATCCGATGGCCTGTTCGCGGGCTTGCCCGCCCCCACAGGaa                                                                                                                                    >  2:1201591/7‑150 (MQ=255)
                        tGAGGATGCAATGATCGAAGCGGAGCAGATGCGCGATGAGCGGCTGGATGAGCTGACCGCTCAAGCCCAGGAATGTGTGCGCCAGCGCAAGAACGTGCATTGATCCGATGGCCTGTTCGCGGGCTCGCCCGCTCCCACAGGtacatcacc                                                                                                                             >  1:2391201‑M1/1‑141 (MQ=255)
                             aTGCAATGATCGAAGCGGAGCAGATGCGCGATGAGCGGCTGGATGAGCTGACCGCTCAAGCCCAGGAATGTGTGCGCCAGCGCAAGAACGTGCATTGATCCGATGGCCTGTTCGCGGGCAAGCCcgcgaa                                                                                                                                            >  1:953284/1‑128 (MQ=255)
                             aTGCAATGATCGAAGCGGAGCAGATGCGCGATGAGCGGCTGGATGAGCTGACCGCTCAAGCCCAGGAATGTGTGCGCCAGCGCAAGAACGTGCATTGATCCGATGGCCTGTTCGCGGGCAAGCCcgcgaa                                                                                                                                            <  2:953284/130‑3 (MQ=255)
                               gCAATGATCGAAGCGGAGCAGATGCGCGATGAGCGGCTGGATGAGCTGACCGCTCAAGCCCAGGAATGTGTGCGCCAGCGCAAGAACGTGCATTGATCCGATGGCCTGTTCGCGGGCTCGCCCGCTCCCACAGGtacatcaccggccttg                                                                                                                      <  2:2391201‑M1/150‑17 (MQ=255)
                                 aaTGATCGAAGCGGAGCAGATGCGCGATGAGCGGCTGGATGAGCTGACCGCTCAAGCCCAGGAATGTGTGCGCCAGCGCAAGAACGTGCATTGATCCGATGGCCTGTTCGCGGGCTTGCCCGCGCCCACAGGAATACCGTAAGGCCCAAg                                                                                                                    >  2:820810/1‑150 (MQ=255)
                                    gATCGAAGCGGAGCAGATGCGCGATGAGCGGCTGGATGAGCTGACCGCTCAAGCCCAGGAATGTGTGCGCCAGCGCAAGAACGTGCATTGATCCGATGGCCTGTTCGCGGGCTTGCCCGCGCCCACAGGAATACCGTAAGGCCCAAGGcc                                                                                                                 >  2:1567534/1‑150 (MQ=255)
                                      tCGAAGCGGAGCAGATGCGCGATGAGCGGCTGGATGAGCTGACCGCTCAAGCCCAGGAATGTGTGCGCCAGCGCAAGAACGTGCATTGATCCGATGGCCTGTTCGCGGGCTTGCCCGCGCCCACAGGAATACCGTAAGGCCCAGGCCGGt                                                                                                              <  1:245895/150‑1 (MQ=255)
                                       cGAAGCGGAGCAGATGCGCGATGAGCGGCTGGCTGAGCCCTGTTCGTAATGCCAGGAATGTGTGCGCCAGAGTTCAGACGTGCATTGATCCGATGGCCTGTTCGCGGGCTTGCCCGCGCCCACAGGAATACCGTAAGGCCCAAGGCCGGt                                                                                                              <  1:1739081/150‑1 (MQ=255)
                                          aGCGGAGCAGATGCGCGATGAGCGGCTGGATGAGCTGACCGCTCAAGCCCAGGAATGTGTGCTCCAGCGCAAGAACGTGCATTGATCCGATGGCCTGTTCGCGGGCTTGCCCGCGCCCACAGGAATACCGTAAGGCCCAAGGCCGGTGAt                                                                                                           <  1:1789604/150‑1 (MQ=255)
                                              gAGCAGATGCGCGATGAGCGGCTGGATGAGCTGACCGCTCAAGCCCAGGAATGTGTGCGCCAGCGCAAGAACGTGCATTGATCCGATGGCCTGTTCGCGGGCTTGCCCGCGCCCACAGGAATACCGTAAGGCCCAAGGCCGGTGATGTAc                                                                                                       >  1:894148‑M2/1‑149 (MQ=255)
                                              gAGCAGATGCGCGATGAGCGGCTGGATGAGCTGACCGCTCAAGCCCAGGAATGTGTGCGCCAGCGCAAGAACGTGCATTGATCCGATGGCCTGTTCGCCGGCTTGCCCGCGCCCACAGGAATACCGTAAGGCCCAAGGCCGGTGATGTAc                                                                                                       >  1:788701‑M2/1‑149 (MQ=255)
                                                  aGATGCGCGATGAGCGGCTGGATGAGCTGACCGCTCAAGCCCAGGAATGTGTGCGCCAGCGCAAGAACGTGCATTGATCCGATGGCCTGTTCGCGGGCTTGCCCGCGCCCACAGGtacatcaccggccttgggccttacggtattcctgt                                                                                                   <  2:1970067‑M1/150‑36 (MQ=255)
                                                   gATGCGCGATGAGCGGCTGGATGAGCTGACCGCTCAAGCCCAGGAATGTGTGCGCCAGCGCAAGAACGTGCATTGATCCGATGGCCTGTTCGCGGGCTTGCCCGCGCCCACAGGAATACCGTAAGGCCCAAGGCCGGTGATGTAcctgtg                                                                                                  >  2:1371187‑M2/1‑144 (MQ=255)
                                                   gATGCGCGATGAGCGGCTGGATGAGCTGACCGCTCAAGCCCAGGAATGTGTGCGCCAGCGCAAGAACGTGCATTGATCCGATGGCCTGTTCGCGGGCTCGCCCGCTCCCCCAGGtacatcaccggccttgggccttacggtattcctgtg                                                                                                  <  2:2271148‑M1/150‑37 (MQ=255)
                                                      gcgcGATGAGCGGCTGGATGAGCTGACCGCTCAAGCCCAGGAATGTGTGCGCCAGCGCAAGAACGTGCATTGATCCGATGGCCTGTTCGCGGGCTTGCCCGCGCCCACAGGAATACCGTAAGGCCCAAGGCCGGTGATGTACCTGTGGGa                                                                                               <  1:1882425/150‑1 (MQ=255)
                                                          gATGAGCGGCTGGATGAGCTGACCGCTCAAGCCCAGGAATGTGTGCGCCAGCGCAAGAACGTGCATTGATCCGATGGCCTGTTCGCGGGCTTGCCCGCGCCCACAGGtacatcaccggccttgggccttacggtattcctgtgggcgcgg                                                                                           <  1:2167945‑M1/150‑44 (MQ=255)
                                                            tGAGCGGCTGGATGAGCTGACCGCTCAAGCCCAGGAATGTGTGCGCCAGCGCAAGAACGTGCATTGATCCGATGGCCTGTTCGCGGGCTTGCCCGCGCCCACAGGAATACCGTAAGGCCCAAGGCCGGTGATGTACCTGTGGGAGCGGGc                                                                                         <  2:2263962/150‑1 (MQ=255)
                                                              aGCGGCTGGATGAGCTGACCGCTCAAGCCCAGGAATGTGTGCGCCAGCGCAAGAACGTGCATTGATCCGATGGCCTGTTCGCGGGCTTGCCCGCGCCCAGGGGAATACCGTAAGGCCCAAGGCCGGTGATGTACCTGTGGGAGCGGGCGa                                                                                       <  2:944295/150‑1 (MQ=255)
                                                               gCGGCTGGATGAGCTGACCGCTCAAGCCCAGGAATGTGTGCGCCAGCGCAAGAACGTGCATTGATCCGATGGCCTGTTCGCGGGCTTGCCCGCGCCCACAGGAATACCGTAAGGCCCAAGGCCGGTGATGTACCTGTGGGAGCGGGCGAg                                                                                      >  2:1266145/1‑150 (MQ=255)
                                                               gCGGCTGGATGAGCTGACCGCTCAAGCCCAGGAATGTGTGCGCCAGCGCAAGAACGTGCATTGATCCGATGGCCTGTTCGCGGGCTCGCCCGCTCCCACAGGtacatcaccggccttgggccttacggtattcctgtgggcgcgggcaag                                                                                      <  1:704180‑M1/150‑49 (MQ=255)
                                                                      gATGAGCTGACCGCTCAAGCCCAGGAATGTGTGCGCCAGCGCAAGAACGTGCATTGATCCGATGGCCTGTTCGCGGGCTTGCCCGCGCCCACAGGAATACCGTAAGGCCCAAGGCCGGTGAt                                                                                                           <  2:1118432/122‑1 (MQ=255)
                                                                      gATGAGCTGACCGCTCAAGCCCAGGAATGTGTGCGCCAGCGCAAGAACGTGCATTGATCCGATGGCCTGTTCGCGGGCTTGCCCGCGCCCACAGGAATACCGTAAGGCCCAAGGCCGGTGAt                                                                                                           >  1:1118432/1‑122 (MQ=255)
                                                                       aTGAGCTGACCGCTCAAGCCCAGGAATGTGTGCGCCAGCGCAAGAACGTGCATTGATCCGATGGCCTGTTCGCGGGCTTGCCCGCGCCCACAGGAATACCGTa                                                                                                                             <  2:2237481/103‑1 (MQ=255)
                                                                       aTGAGCTGACCGCTCAAGCCCAGGAATGTGTGCGCCAGCGCAAGAACGTGCATTGATCCGATGGCCTGTTCGCGGGCTTGCCCGCGCCCACAGGAATACCGTa                                                                                                                             >  1:2237481/1‑103 (MQ=255)
                                                                          aGCTGACCGCTCAAGCCCAGGAATGTGTGCGCCAGCGCAAGAACGTGCATTGATCCGATGGCCTGTTCGCGGGCTTGCCCGCGCCCACAGGAATACCGTAAGGCCCAAGGCCGGTGATGTACCTGTGGGAGCGGGCGAGCCCGCgaagaa                                                                           <  2:1556927/150‑1 (MQ=255)
                                                                          aGCTGACCGCTCAAGCCCAGGAATGTGTGCGCCAGCGCAAGAACGTGCATTGATCCGATGGCCTGTTCGCGGGCTTGCCCGCGCCCACAGGAATACCGTAAGGCCCAAGGCCGGTGATGTACCTGTGGGAGCGGGCGAGCCCGCgaagaa                                                                           <  2:776613/150‑1 (MQ=255)
                                                                          aGCTGACCGCTCAAGCCCAGGAATGTGTGCGCCAGCGCAAGAACGTGCATTGATCCGATGGCCTGTTCGCGGGCTTGCCCGCGCCCACAGGAATACCGTAAGGCCCAAGGCCGGTGATGTACCTGTGGGAGCGGGCGAGCCCGCgaagaa                                                                           <  2:776497/150‑1 (MQ=255)
                                                                           gCTGACCGCTCAAGCCCAGGAATGTGTGCGCCAGCGCAAGAACGTGCATTGATCCGATGGCCTGTTCGCGGGCTTGCCCGCGCCCACAGGAATACCGTAAGGCCCAAGGCCGGTGATGTACCTGTGGGAGCggg                                                                                          >  2:2088784/1‑134 (MQ=255)
                                                                           gCTGACCGCTCAAGCCCAGGAATGTGTGCGCCAGCGCAAGAACGTGCATTGATCCGATGGCCTGTTCGCGGGCTTGCCCGCGCCCACAGGAATACCGTAAGGCCCAAGGCCGGTGATGTACCTGTGGGAGCGGGCGAGCCCGCGaagaag                                                                          >  2:1302302/1‑150 (MQ=255)
                                                                            cTGACCGCTCAAGCCCAGGAATGTGTGCGCCAGCGCAAGAACGTGCATTGATCCGATGGCCTGTTCGCGGGCTTGCCCGCGCCCACAGGAATACCGTAAGGCCCAAGGCCGGTGATGTACCTGTGGGAGCGGGCGAGCCCGCGAAGAAGc                                                                         <  1:1792635/150‑1 (MQ=255)
                                                                              gACCGCTCAAGCCCAGGAATGTGTGCGCCAGCGCAAGAACGTGCATTGATCCGATGGCCTGTTCGCGGGCTTGCCCGCGCCCACAGGAATACCGTAAGGCCCAAGGCCGGTGATGTACCTGTGGGAGCGGGCGAGCCCGCGAAGAAGccc                                                                       >  1:281129/1‑150 (MQ=255)
                                                                               aCCGCTCAAGCCCAGGAATGTGTGCGCCAGCGCAAGAACGTGCATTGATCCGATGGCCTGTTCGCGGGCTTGCCCGCGCCCACAGGAATACCGTAAGGCCCAAGGCCGGTGATGTACCTGTGGGAGCGGGCGAGCCCGCGAAGAAGcccc                                                                      >  2:1256555/1‑150 (MQ=255)
                                                                               aCCGCTCAAGCCCAGGAATGTGTGCGCCAGCGCAAGAACGTGCATTGATCCGATGGCCTGTTCGCGGGCTTGCCCGCGCCCACAGGAATACCGTAAGGCCCAAGGCCGGTGATGTACCTGTGGGAGCGGGCGAGCCCGCGAAGAAGcccc                                                                      >  1:972490/1‑150 (MQ=255)
                                                                               aCCGCTCAAGCCCAGGAATGTGTGCGCCAGCGCAAGAACGTGCATTGATCCGATGGCCTGTTCGCGGGCTTGCCCGCGCCCACAGGAATACCGTAAGGCCCAAGGCCGGTGATGTACCTGTGGGAGCGGGCGAGCCCGCGAAGAAGcccc                                                                      >  1:1894390/1‑150 (MQ=255)
                                                                                  gCTCAAGCCCATGAATGTGTGCGCCAGCGCAAGAACGTGCATTGATCCGATGGCCTGTTCGCNGGCTTGCCCGCGCCCACAGGAATACCGTAAGGCCCAAGGCCGGTGAt                                                                                                           <  1:724102/110‑1 (MQ=255)
                                                                                  gCTCAAGCCCATGAATGTGTGCGCCAGCGCAAGAACGTGCATTGATCCGATGGCCTGTTCGCGGGCTTGCCCGCGCCCACAGGAATACCGTAAGGCCCAAGGCCGGTGAt                                                                                                           >  2:724102/1‑110 (MQ=255)
                                                                                  gCTCAAGCCCAGGAATGTGTGCGCCAGCGCAAGAACGTGCATTGATCCGATGGCCTGTTCGCGGGCTTGCCCGCGCCCACAGGAATACCGTAAGGCCCAAGGCCGGTGATGTACCTGTGGGAGCGGGCGAGCCCGCGAAGAAGCCCCAtt                                                                   <  2:403059/150‑1 (MQ=255)
                                                                                    tCAAGCCCAGGAATGTGTGCGCCAGCGCAAGAACGTGCATTGATCCGATGGCCTGTTCGCGGGCTTGCCCGCGCCCACAGGAATACCGTAAGGCCCAAGGCCGGTGATGTACCTGTGGGAGCGGGCGAGCCCGCGAAGAAGcc                                                                        <  2:854090/143‑1 (MQ=255)
                                                                                    tCAAGCCCAGGAATGTGTGCGCCAGCGCAAGAACGTGCATTGATCCGATGGCCTGTTCGCGGGCTTGCCCGCGCCCACAGGAATACCGTAAGGCCCAAGGCCGGTGATGTACCTGTGGGAGCGGGCGAGCCCGCGAAGAAGCCCCATTgg                                                                 >  1:2409849/1‑150 (MQ=255)
                                                                                         cccAGGAATGTGTGCGCCAGCGCAAGAACGTGCATTGATCCGATGGCTTCTTCGCGGGCTCGCCCGCTCCCACAGGtacatcaccggccttgggccttacggtattcctgtgggcgc                                                                                             <  1:661482‑M1/117‑42 (MQ=255)
                                                                                         cccAGGAATGTGTGCGCCAGCGCAAGAACGTGCATTGATCCGATGGCTTCTTCGCGGGCTCGCCCGCTCCCACAGGtacatcaccggccttgggccttacggtattcctgtgggcgc                                                                                             >  2:661482‑M1/1‑76 (MQ=255)
                                                                                           cAGGAATGTGTGCGCCAGCGCAAGAACGTGCATTGATCCGATGGCTTCTTCGCGGGCTCGCCCGCTCCCACAGGtacatcaccggcctggggccttacggtattcctgtgggcgcgggcaagcccgcgaacaggccatcggatcaatgca                                                          <  1:1519054‑M1/150‑77 (MQ=255)
                                                                                           cAGGAATGTGTGCGCCAGCGCAAGAACGTGCATTGATCCGATGGCCTGTTCGCGGGCTTGCCCGCGCCCACAGGAATACCGTAAGGCCCAAGGCCGGTGATGTACCTGTGGGAGCGGGCGAGCCCGCGAAGAAGCCCCATTGGCCTTTCa                                                          <  2:505482/150‑1 (MQ=255)
                                                                                            aGGAATGTGTGCGCCAGCGCAAGAACGTGCATTGATCCGATGGCTTCTCCGGGTGCTCGCCCGCTCCCACAGGtacatcaccggccttgggccttacggtattcctgtgggcgcgggcaagcccgcgaacaggccatcggatcaatgcac                                                         <  1:760055‑M1/150‑78 (MQ=255)
                                                                                            aGGAATGTGTGCGCCAGCGCAAGAACGTGCATTGATCCGATGGCCTGTTCGCGGGCTTGCCCGCGCCCACAGGAATACCGTAAGGCCCAAGGCCGGTGATGTACCTGTGGGAg                                                                                              >  1:2391281/1‑113 (MQ=255)
                                                                                            aGGAATGTGTGCGCCAGCGCAAGAACGTGCATTGATCCGATGGCCTGTTCGCGGGCTTGCCCGCGCCCACAGGAATACCGTAAGGCCCAAGGCCGGTGATGTACCTGTGGGAg                                                                                              <  2:2391281/113‑1 (MQ=255)
                                                                                            aGGAATGTGTGCGCCAGCGCAAGAACGTGCATTGATCCGATGGCCTGTTCGCGGGCTTGCCCGCGCCCACAGGAATACCGTAAGGCCCAAGGCCGGTGATGTACCTGTGGGAGCGGGCGAGCCCGCGAAGAAGCCCCATTGGCCTTTCAg                                                         <  1:1266145/150‑1 (MQ=255)
                                                                                                aTGTGTGCGCCAGCGCAAGAACGTGCATTGATCCGATGGCCTGTTCGCGGGCTTGCCCGCGCCCACAGGAATACCGTAAGGCCCAAGGCCGGTGATGTACCTGTGGGAGCGGGCGAGCCCGCGAAGAAGCCCCATTGGCCTTTCAGCCAt                                                     >  2:393132/1‑150 (MQ=255)
                                                                                                 tgtgtgCGCCAGCGCAAGAACGTGCATTGATCCGATGGCCTGTTCGCGGGCTTGCCCGCGCCCACAGGAATACCGTAAGGCCCAAGGCCGGTGATGTACCTGTGGGAGCGGGCGAGCCCGCGAAGAAGCCCCATTGGCCTTTCAGCCATg                                                    <  1:1256555/150‑1 (MQ=255)
                                                                                                         ccAGCGCAAGAACGTGCATTGATCCGATGGCCTGTTCGCGGGCTTGCCCGCGCCCACAGGtacatcaccggccttgggccttacggtat                                                                                                         >  1:705238‑M1/1‑60 (MQ=255)
                                                                                                              gcAAGAACTTGCATTGATCCGAGGCCCTGTTCGCAGCCTTGCCCGCGCCCACAGGAATACCGTAAGGCCCAAGGCCGGTGATGTACCTGTGGGAGCGGGCGAGCCCGCGAAGAAGCCCCATTGGCCTTTCAGCCATGCACCCAGCGCCGa                                       <  2:1083728/150‑1 (MQ=255)
                                                                                                                aaGAACGTGCATTGATCCGATGGCCTGTTCGCGGGCTTGCCCGCGCCCACAGGAATACCGTAAGGCCCAAGGCCGGt                                                                                                              <  1:270844/77‑1 (MQ=255)
                                                                                                                aaGAACGTGCATTGATCCGATGGCCTGTTCGCGGGCTTGCCCGCGCCCACAGGAATACCGTAAGGCCCAAGGCCGGTGATGTACCTGTGGGAGCGGGCGAGCCCGCGAAGAAGCCCCATTGGCCTTTCAGCCATGCACCCAGCGCCGATg                                     >  2:1578360/1‑150 (MQ=255)
                                                                                                                  gAACGTGCATTGATCCGATGGCCTGTTCGCGGGCTTGCCCGCGCCCACAGGAATACCGTAAGGCCCAAGGCCGGTGATGTACCTGTGGGAGCGGGCGAGCCCGCGAAGAAGCCCCATTGGCCTTTCa                                                          <  2:2028790/127‑1 (MQ=255)
                                                                                                                  gAACGTGCATTGATCCGATGGCCTGTTCGCGGGCTTGCCCGCGCCCACAGGAATACCGTAAGGCCCAAGGCCGGTGATGTACCTGTGGGAGCGGGCGAGCCCGCGAAGAAGCCCCATTGGCCTTTCa                                                          >  1:2028790/1‑127 (MQ=255)
                                                                                                                  gAACGTGCATTGATCCGATGGCCTGTTCGCGGGCTTGCCCGCGCCCACAGGAATACCGTAAGGCCCAAGGCCGGTGATGTACCTGTGGGAGCGGGCGAGCCCGCGAAGAAGCCCCATTGGCCTTTCAGCCATGCACCCAGCGCCGATGCa                                   <  1:370054/150‑1 (MQ=255)
                                                                                                                      gTGCATTGATCCGATGGCCTGTTCGCGGGCTTGCCCGCGCCCACAGGAATACCGTAAGGCCCAAGGCCGGTGATGTAcctgtgggagcgggcgagcacccggagaagccatcggatcaatgcacgttcttgcgctggcgcacacattcct                               <  1:760055‑M2/150‑74 (MQ=255)
                                                                                                                      gTGCATTGATCCGATGGCCTGTTCGCGGGCTTGCCCGCGCCCACAGGAATACCGTAAGGCCCAAGGCCGGTGATGTACCTGTGGGAGCGGGCGAGCCCGCGAAGAAGCCCCAtt                                                                   <  2:1225530/114‑1 (MQ=255)
                                                                                                                       tGCATTGATCCGATGGCCTGTTCGCGGGCTTGCCCGCGCCCACAGGAATACCGTAAGGCCCCAGGCCGGTGATGTAcctgtgggagcgggcgagcccgcgaagaagccatcggatcaatgcacgttcttgcgctggcgcacacattcctg                              <  1:1519054‑M2/150‑75 (MQ=255)
                                                                                                                             gATCCGATGGCCTGTTCGCGGGCTTGCCCGCGCCCACAGGAATACCGTAAGGCCCAAGGCCGGTGATGTACCTGTGGGAGCgg                                                                                           >  1:560803/1‑83 (MQ=255)
                                                                                                                                ccGATGGCCTGTTCGCGGGCTTGCCCGCGCCCACAGGAATACCGTAAGGCCCAAGGCCGGTGATGTACCTGTGGGAGCGGGCGAGCCCGCGAAGAAGCCCCATTGGCCTTTCAGCCATGCACCCAGCGCCGATGCAGCCCACGCGCCaga                     <  2:370413/150‑1 (MQ=255)
                                                                                                                                ccGATGGCCTGTTCGCGGGCTTGCCCGCGCCCACAGGAATACCGTAAGGCCCAAGGCCGGTGATGTACCTGTGGGAGCGGGCGAGCCCGCGAAGAAGCCCCATTGGCCTTTCAGCCATGCACCCAGCGCCGATGCAGCCCACGCGCCaga                     <  1:203579/150‑1 (MQ=255)
                                                                                                                                ccGATGGCCTGTTCGCGGGCTTGCCCGCGCCCACAGGAATACCGTAAGGCCCAAGGCCGGTGATGTACCTGTGGGAGCGGGCGAGCCCGCGAAGAAGCCCCATTGGCCTTTCAGCCATGCACCCAGCGCCGATGCAGCCCACGCGCCaga                     <  1:1415451/150‑1 (MQ=255)
                                                                                                                                ccGATGGCCTGTTCGCGGGCTTGCCCGCGCCCACAGGAATACCGTAAGGCCCAAGGCCGGTGATGTACCTGTGGGAGCGGGCGAGCCCGCGAAGAAGCCCCATTGGCCTTTCAGCCATGCACCCAGCGCCGATGCAGCCCACGCGCCaga                     <  1:1567534/150‑1 (MQ=255)
                                                                                                                                 cGATGGCCTGTTCGCGGGCTTGCCCGCGCCCACAGGAATACCGTAAGGCCCAAGGCCGGTGATGTACCTGTGGGAGCGGGCGAGCCCGCGAAGAAGCCCCATTGGCCTTTCAGCCATGCACCCAGCGCCGATGCAGCCCACGCGCCagag                    <  1:1302302/150‑1 (MQ=255)
                                                                                                                                 cGATGGCCTGTTCGCGGGCTTGCCCGCGCCCACAGGAATACCGTAAGGCCCAAGGCCGGTGATGTACCTGTGGGAGCGGGCGAGCCCGCGAAGAAGCCCCATTGGCCTTTCAGCCATGCACCCAGCGCCGATGCAGCCCACGCGCCagag                    >  1:2116946/1‑150 (MQ=255)
                                                                                                                                   aTGGCCTGTTCGCGGGCTTGCCCGCGCCCACAGGAATACCGTAAGGCCCAAGGCCGGTGATGTACCTGTGGGAGCGGGCGAGCCCGCGAAGAAGCCCCATTGGCCTTTCAGCCATGCTCCCAGCGCCGATGCAGCCCACgcgc                         >  1:89906/1‑143 (MQ=255)
                                                                                                                                   aTGGCCTGTTCGCGGGCTTGCCCGCGCCCACAGGAATACCGTAAGGCCCAAGGCCGGTGATGTACCTGTGGGAGCGGGCGAGCCCGCGAAGAAGCCCCATTGGCCTTTCAGCCATGCACCCAGCGCCGATGCAGTCCACGCGCCAGAGCa                  <  2:1415035/150‑1 (MQ=255)
                                                                                                                                    tGGCCTGTTCGCGGGCTTGCCCGCGCCCACAGGAATACCGTAAGGCCCAAGGCCGGt                                                                                                              <  1:1421920/57‑1 (MQ=255)
                                                                                                                                     ggCCTGTTCGCGGGCTTGCCCGCGCCCACAGGAATACCGTAAGGCCCAAGGCCGGTGATGTACCTGTGGGAGCGGGCGAGCCCGCGAAGAAGCCCCATTGGCCTTTCAGCCATGCACCCAGCGCCGATGCAGCCCACGCGCCAGAGCATc                <  1:1578360/150‑1 (MQ=255)
                                                                                                                                      gCCTGTTCGCGGGCTTGCCCGCGCCCACAGGAATACCGTAAGGCCCAAGGCCGGTGATGTACCTGTGGGAGCGGGCGAGCCCGCGAAGAAGCCCCATTGGCCTTTCAGCCATGCACCCAGCGCCGATGCAGCCCACGCGCCAGAGCATcc               <  1:820810/150‑1 (MQ=255)
                                                                                                                                      gCCTGTTCGCGGGCTTGCCCGCGCCCACAGGAATACCGTAAGGCCCAAGGCCGGTGATGTACCTGTGGGAGCGGGCGAGCCCGCGAAGAAGCCCCATTGGCCTTTCAGCCATGCACCCAGCGCCGATGCAGCCCACGCGCCAGAGCATcc               <  1:2220811/150‑1 (MQ=255)
                                                                                                                                       ccTGTTCGCGGGCTTGCCCGCGCCCACAGGAATACCGTAAGGCCCAAGGCCGGTGATGTACCTGTGGGAGCGGGCGAGCCCGCGAAGAAGCCCCATTGGCCTTTCAGCCATGCACCCAGCGCCGATGCAGCCCACGCGCCAGAGCATCCa              >  1:1343479/1‑150 (MQ=255)
                                                                                                                                        cTGTTCGCGGGCTTGCCCGCGCCCACAGGAATACCGTAAGGCCCAAGGCCGGTGATGTAcctgt                                                                                                   >  2:1839347‑M2/1‑59 (MQ=255)
                                                                                                                                        cTGTTCGCGGGCTTGCCCGCGCCCACAGGAATACCGTAAGGCCCAAGGCCGGTGATGTACCTGTGGGAGCGGGCGAGCCCGCGAAGAAGCCCCATTGGCCTTTCAGCCATGCACCCAGCGCCGATGCAGCCCACGCGCCAGAGCATCCAg             <  2:44490/150‑1 (MQ=255)
                                                                                                                                        cTGTTCGCGGGCTTGCCCGCGCCCACAGGAATACCGTAAGGCCCAAGGCCGGTGATGTACCTGTGGGAGCGGGCGAGCCCGCGAAGAAGCCCCATTGGCCTTTCAGCCATGCACCCAGCGCCGATGCAGCCCACGCGCCAGAGCATCCAg             <  2:281129/150‑1 (MQ=255)
                                                                                                                                         tGTTCGCGGGCTTGCCCGCGCCCACAGGAATACCGTAAGGCCCAAGGCCGGTGATGTAcctgtgggagcgggcgagcccgcgaagaggccatcggatcaatgcacgttcttgcgctggcgcacacattcctgggcttgagcggtcagctc            <  2:1100783‑M2/150‑93 (MQ=255)
                                                                                                                                         tGTTCGCGGGCTTGCCCGCGCCCACAGGAATACCGTAAGGCCCAAGGCCGGTGATGTACCTGTGGGAGCGGGCGAGCCCGCGAAGAAGCCCCATTGGCCTTTCAGCCATGCACCCAGCGCCGATGCAGCCCACGCGCCAGAGCATCCAGc            >  2:191475/1‑150 (MQ=255)
                                                                                                                                            gcgcgGGCAAGCCCGCGCCCACAGGAATACCGTAAGGCCCAAGGCCGGTGATGTACCTGTGGGAGCGGGCGAGCCCGCGAAGAAGCCCCATTGGCCTTTCAGCCATGCACCCAGCGCCGATGCAGcc                                <  2:1767060/126‑1 (MQ=255)
                                                                                                                                            gcgcgGGCAAGCCCGCGCCCACAGGAATACCGTAAGGCCCAAGGCCGGTGATGTACCTGTGGGAGCGGGCGAGCCCGCGAAGAAGCCCCATTGGCCTTTCAGCCATGCACCCAGCGCCGATGCAGcc                                >  1:1767060/2‑127 (MQ=255)
                                                                                                                                            gcgcgGGCAAGCCCGCGCCCACAGGAATACCGTAAGGCCCAAGGCCGGTGATGTACCTGTGGGAGCGGGCGAGCCCGCGAAGAAGCCCCATTGGCCTTTCAGCCATGCACCCAGCGCCGATGCAGCCCACGCGCCAGAGCATCCAGca           <  2:2014707/147‑1 (MQ=255)
                                                                                                                                                gggCTTGCCCGCGCCCACAGGAATACCGTAAGGCCCAAGGCCGGTGATGTACCTGTGGGAGCGGGCGAGCCCGCGAAGAAGCCCCATTGGCCttt                                                            <  2:653044/95‑1 (MQ=255)
                                                                                                                                                  gCTTGCCCGCGCCCACAGGAATAACGTAGGGCCCAAGGCCAGTGATGTACCTGTGGGTGCGGGCGAGCCCGCGAAGAAGCCCCATTGGCCTTTCAGCCATGCACCCAGCGCCGATGCAGCCCACGCGCCAGAGCATCCAGCACAAAcccc   >  1:1529936/1‑150 (MQ=255)
                                                                                                                                                     |                                                                                                                                                   
CACCTATCTGCTTAACCAGCAGGGTGAGGATGCAATGATCGAAGCGGAGCAGATGCGCGATGAGCGGCTGGATGAGCTGACCGCTCAAGCCCAGGAATGTGTGCGCCAGCGCAAGAACGTGCATTGATCCGATGGCCTGTTCGCGGGCTTGCCCGCGCCCACAGGAATACCGTAAGGCCCAAGGCCGGTGATGTACCTGTGGGAGCGGGCGAGCCCGCGAAGAAGCCCCATTGGCCTTTCAGCCATGCACCCAGCGCCGATGCAGCCCACGCGCCAGAGCATCCAGCACAAACCCCA  >  NC_002947/4131006‑4131302

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