Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 1,700,014 G→A 5.9% R330Q (CGG→CAG)  wspC → biofilm formation methyltransferase WspC

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0029471,700,0140GA 5.9% 1544.5 / 51.6 486R330Q (CGG→CAG) wspCbiofilm formation methyltransferase WspC
Reads supporting (aligned to +/- strand):  ref base G (202/255);  new base A (15/14);  total (217/269)
Fisher's exact test for biased strand distribution p-value = 4.47e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

Only 100 of 489 total aligned reads displayed.

TGCGTTCACGACGTTGCCCGCACCGGTTTATCCGCAGCCCAGTGTGCCGTTGCCTCGTTCGCGGCGTGTGCTGCCGGTGGCGGCCAGGCCAGCGCGAGCCCGCGAGCACAGCCACGAGGGGGCCAGCGAGTTGCTGGCCGGCATTGCCCGGCTGGCCAATGCCGGTGCCAGCGAGCAGGCCCGTAGCGAATGCCAACGCTACCTGAGCCAGTACCCGCCGTCGGCACAGGTGTACTACTGGCTGGGTTTGCTCAGTGATACCGAAGGCGATGCGCAGCAGGCGCTGAGCCACTATCGC  >  NC_002947/1699865‑1700162
                                                                                                                                                     |                                                                                                                                                    
tgcgtTCACGACGTTGCCCGCACCGGTTTATCCGCAGCCCAGTGTGCCGTTGCCTCGTTCGCGGCGTGTGCTGCCGGTGGCGGCCAGGCCAGCGCGAGCCCGCGAGCACAGCCACGAGGGGGCCAGCGAGTTGCTGGCCGGCATTGCCCg                                                                                                                                                      <  1:2070393/150‑1 (MQ=255)
  cgtTCACGACGTTGCCCGCACCGGTTTATCCGCAGCCCAGTGTGCCGTTGCCTCGTTCGCGGCGTGTGCTGCCGGTGGCGTCCAGGCCAGCGCGAGCCCGCGAGCACAGCCACGAGGGGGCCAGCGAGTTGCTGGCCGGCATTGCCCGGc                                                                                                                                                    <  2:3827647/150‑1 (MQ=255)
   gtTCACGACGTTGCCCGCACCGGTTTATCCGCAGCCCAGTGTGCCGTTGCCTCGTTCGCGGCGTGTGCTGCCGGTGGCGGCCAGGCCAGCGCGAGCCCGCGAGCACAGCCACGAGGGGGCCAGCGAGTTGCTGGCCGGCATTGCCCGGCt                                                                                                                                                   <  1:7325967/150‑1 (MQ=255)
   gtTCACGACGTTGCCCGCACCGGTTTATCCGCAGCCCAGTGTGCCGTTGCCTCGTTCGCGGCGTGTGCTGCCGGTGGCGGCCAGGCCAGCGCGAGCCCGCGAGCACAGCCACGAGGGGGCCAGCGAGTTGCTGGCCGGCATTGCCCGGCt                                                                                                                                                   >  1:5213470/1‑150 (MQ=255)
   gtTCACGACGTTGCCCGCACCGGTTTATCCGCAGCCCAGTGTGCCGTTGCCTCGTTCGCGGCGTGTGCTGCCGGTGGCGGCCAGGCCAGCGCGAGCCCGCGAGCACAGCCACGAGGGGGCCAGCGAGTTGCTGGCCGGCATTGCCCGGCt                                                                                                                                                   >  2:5446327/1‑150 (MQ=255)
   gtTCACGACGTTGCCCGCACCGGTTTATCCGCAGCCCAGTGTGCCGTTGCCTCGTTCGCGGCGTGTGCTGCCGGTGGCGGCCAGGCCAGCGCGAGCCCGCGAGCACAGCCACGAGGGGGCCAGCGAGTTGCTGGCCGGCATTGCCCGGCt                                                                                                                                                   <  2:1852916/150‑1 (MQ=255)
   gtTCACGACGTTGCCCGCACCGGTTTATCCGCAGCCCAGTGTGCCATTGCCTCGTTCGCGGCGTGTGCTGCCGGTGGCGGCCAGGCCAGCGCGAGCCCGCGAGCACAGCCACGAGGGGGCCAGCGAGTTGCTGGCCGGCATTGCCCGGCt                                                                                                                                                   <  2:5741331/150‑1 (MQ=255)
    gtCACGACGGTGCCCGCACCGGTTTATCCGCAGCCCAGTGTGCCGTTGCCTCGTTCGCGGCGTGTGCTGCCGGTGGCGGCCAGGCCAGCGCGAGCCCGCGAGCACAGCCACGAGGGGGCCAGCGAGTTGCTGGCCGGCATTGCCCGGCTg                                                                                                                                                  <  2:6358246/149‑1 (MQ=255)
      cACGACGTTGCCCGCACCGGTTTATCCGCAGCCCAGTGTGCCGTTGCCTCGTTCGCGGCGTGTGCTGCCGGTGGCGGCCAGGCCAGCGCGAGCCCGCGAGCACAGCCACGAGGGGGCCAGCGAGTTGCTGGCCGGCATTGCCCGGCTGGc                                                                                                                                                <  2:4062851/150‑1 (MQ=255)
          acgTTGCCCGCACCGGTTTATCCGCAGCCCAGTGTGCCGTTGCCTCGTTCGCGGCGTGTGCTGCCGGTGGCGGCCAGGCCAGCGCGAGCCCGCGAGCACAGCCACGAGGGGGCCAGCGAGTTGCTGGCCGGCATTGCCCGGCTGGCCAAt                                                                                                                                            <  2:3146793/150‑1 (MQ=255)
            gTTGCCCGCACCGGTTTATCCGCAGCCCAGTGTGCCGTTGCCTCGTTCGCGGCGTGTGCTGCCGGTGGCGGCCAGGCCAGCGCGAGCCCGCGAGCACAGCCACGAGGGGGCCAGCGAGTTGCTGGCCGGCATTGCCCGGCTGGCCAATGc                                                                                                                                          <  2:7681680/150‑1 (MQ=255)
            gTTGCCCGCACCGGTTTATCCGCAGCCCAGTGTGCCGTTGCCTCGTTCGCGGCGTGTGCTGCCGGTGGCGGCCAGGCCAGCGCGAGCCCGCGAGCACAGCCACGAGGGGGCCAGCGAGTTGCTGGCCGGCATTGCCCGGCTGGCCAAAGc                                                                                                                                          >  1:1744355/1‑150 (MQ=255)
            gTTGCCCGCACCGGTTTATCCGCAGCCCAGTGTGCCGTTGCCTCGTTCGCGGCGTGTGCTGCCGGTGGCGGCCAGGCCAGCGCGAGCCCGCGAGCACAGCCACGAGGGGGCCAGCGAGTTGCTGACCGGCATTGGCCAGCCGGGCAATGc                                                                                                                                          <  2:1583838/150‑1 (MQ=255)
             ttGCCCGCACCGGTTTATCCGCAGCCCAGTGTGCCGTTGCCTCGTTCGCGGCGTGTGCTGCCGGTGGCGGCCAGGCCAGCGCGAGCCCGCGAGCACAGCCACGAGGGGGCCAGCGAGTTGCTGGCCGGCATTGCCCGGCTGGCCAATGcc                                                                                                                                         <  1:5972669/150‑1 (MQ=255)
              tGCCCGCACCGGTTTATCCGCAGCCCAGTGTGCCGTTGCCTCGTTCGCGGCGTGTGCTGCCGGTGGCGGCCAGGCCAGCGCGAGCCCGCGAGCACAGCCACGAGGGGGCCAGCGAGTTGCTGGCCGGCATTGGCCAGCCGGGCAATGCCg                                                                                                                                        <  1:7967525/150‑1 (MQ=255)
              tGCCCGCACCGGTTTATCCGCAGCCCAGTGTGCCGTTGCCTCGTTCGCGGCGTGTGCTGCCGGTGGCGGCCAGGCCAGCGCGAGCCCGCGAGCACAGCCACGAGGGGGCCAGCGAGTTGCTGGCCGGCATTGCCCGGCTGGCCAATGCCg                                                                                                                                        <  2:418389/150‑1 (MQ=255)
              tGCCCGCACCGGTTTATCCACAGCCCAGTGTGCCGTTGCCTCGTTCGCGGCGTGTGCTGCCGGTGGCGGCCAGGCCAGCGCGAGCCCGCGAGCACAGCCACGAGGGGGCCAGCGAGTTGCTGGCCGGCATTGCCCGGCTGGCCAATGCCg                                                                                                                                        <  2:6183825/150‑1 (MQ=255)
              tGCCCGCACCGGTTTATCCACAGCCCAGTGTGCCGTTGCCTCGTTCGCGGCGTGTGCTGCCGGTGGCGGCCAGGCCAGCGCGAGCCCGCGAGCACAGCCACGAGGGGGCCAGCGAGTTGCTGGCCGGCATTGCCCGGCTGGCCAATGCCg                                                                                                                                        <  2:5965549/150‑1 (MQ=255)
                cccGCACCGGTTTATCCGCAGCCCAGTGTGCCGTTGCCTCGTTCGCGGCGTGTGCTGCCGGTGGCGGCCAGGCCAGCGCGAGCCCGCGAGCACAGCCACGAGGGGGCCAGCGAGTTGCTGACCGGCATTGGCCAGCCGGGCAATGCCggc                                                                                                                                      <  1:8168136/150‑2 (MQ=255)
                cccGCACCGGTTTATCCGCAGCCCAGTGTGCCGTTGCCTCGTTCGCGGCGTGTGCTGCCGGTGGCGGCCAGGCCAGCGCGAGCCCGCGAGCACAGCCACGAGGGGGCCAGCGAGTTGCTGACCGGCATTGGCCAGCCGGGCAATGCCggc                                                                                                                                      <  2:220615/150‑2 (MQ=255)
                 ccGCACCGGTTTATCCGCAGCCCAGTGTGCCGTTGCCTCGTTCGCGGCGTGTGCTGCCGGTGGCGGCCAGGCCAGCGCGAGCCCGCGAGCACAGCCACGAGGGGGCCAGCGAGTTGCTGGCCGGCATTGCCCGGCTGGCCAATGCCGGTg                                                                                                                                     >  1:3743695/1‑150 (MQ=255)
                 ccGCACCGGTTTATCCGCAGCCCAGTGTGCCGTTGCCTCGTTCGCGGCGTGTGCTGCCGGTGGCGGCCAGGCCAGCGCGAGCCCGCGAGCACAGCCACGAGGGGGCCAGCGAGTTGCTGGCCGGCATTGCCCGGCTGGCCAATGCCGGTg                                                                                                                                     >  1:1247061/1‑150 (MQ=255)
                    cACCGGTTTATCCGCAGCCCAGTGTGCCGTTGCCTCGTTCGCGGCGTGTGCTGCCGGTGGCGGCCAGGCCAGCGCGAGCCCGCGAGCACAGCCACGAGGGGGCCAGCGAGTTGCTGGCCGGCATTGCCCGGCTGGCCAATGCCGGTGCCa                                                                                                                                  <  1:7874692/150‑1 (MQ=255)
                        ggTTTATCCGCAGCCCAGTGTGCCGTTGCCTCGGTCGCGGCGTGTGCTGCCGGTGGCGGCCAGGCCAGCGCGAGCCCGCGAGCACAGCCACGAGGGGGCCAGCGAGTTGCTGGCCGGCATTGCCCGGCTGGCCAATGCCGGTGCCAGCGa                                                                                                                              <  2:8424593/150‑1 (MQ=255)
                          tttATCCGCAGCCCAGTGTGCCGTTGCCTCGTTCGCGGCGTGTGCTGCCGGTGGCGGCCAGGCCAGCGCGAGCCCGCGAGCACAGCCACGAGGGGGCCAGCGAGTTGCTGGCCGGCATTGCCCGGCTGGCCAATGCCGGTGCCAGCGAGc                                                                                                                            <  1:8441581/150‑1 (MQ=255)
                           ttATCCGCAGCCCAGTGTGCCGTTGCCTCGTTCGCGGCGTGTGCTGCCGGTGGCGGCCAGGCCAGCGCGAGCCCGCGAGCACAGCCACGAGGGGGCCAGCGAGTTGCTGGCCGGCATTGCCCGGCTGGCCAATGCCGGTGCCGGCGAgac                                                                                                                           <  2:2065482/150‑3 (MQ=255)
                           ttATCCGCAGCCCAGTGCGCCGTTGCCTCGTTCGCGGCGTGTGCTGCCGGTGGCGGCCAGGCCAGCGCGAGCCCGCGAGCACAGCCACGAGGGGGCCAGCGAGTTGCTGGCCGGCATTGCCCGGCTGGCCAATGCCGGTGCCAGCGCGCa                                                                                                                           <  2:4935317/150‑1 (MQ=255)
                            tATCCGCAGCCCAGTGTGCCGTTGCCTCGTTCGCGGCGTGTGCTGCCGGTGGCGGCCAGGCCAGCGCGAGCCCGCGAGCACAGCCACGAGGGGGCCAGCGAGTTGCTGGCCGGCATTGCCCGGCTGGCCAATGCCGGTGCCAGCGAGCAg                                                                                                                          <  1:5499905/150‑1 (MQ=255)
                             aTCCGCAGCCCAGTGTGCCGTTGCCTCGTTCGCGGCGTGTGCTGCCGGTGGCGGCCAGGCCAGCGCGAGCCCGCGAGCACAGCCACGAGGGGGCCAGCGAGTTGCTGGCCGGCATTGCCCGGCTGGCCAATGCCGGTGCCAGCGAGCAgg                                                                                                                         <  1:1918929/150‑1 (MQ=255)
                                cGCAGCCCAGTGTGCCGTTGCCTCGTTCGCGGCGTGTGCTGCCGGTGGCGGCCAGGCCAGCGCGAGCCCGCGAGCACAGCCACGAGGGGGCCAGCGAGTTGCTGGCCGGCATTGCCCGGCTGGCCAATGCCGGTGCCAGCGAGCAGGccc                                                                                                                      <  1:1374235/150‑1 (MQ=255)
                                cGCAGCCCAGTGTGCCGTTGCCTCGTTCGCGGCGTGTGCTGCCGGTGGCGGCCAGGCCAGCGCGAGCCCGCGAGCACAGCCACGAGGGGGCCAGCGAGTTGCTGGCCGGCATTGCCCGGCTGGCCAATGCCGGTGCCAGCGAGCAGGccc                                                                                                                      >  1:516493/1‑150 (MQ=255)
                                 gCAGCCCAGTGTGCCGTTGCCTCGTTCGCGGCGTGTGCTGCCGGTGGCGGCCAGGCCAGCGCGAGCCCGCGAGCACAGCCACGAGGGGGCCAGCGAGTTGCTGGCCGGCATTGCCCGGCTGGCCAATGCCGGTGCCAGCGAGCAGGCCCg                                                                                                                     >  2:28276/1‑150 (MQ=255)
                                 gCAGCCCAGTGTGCCGTTGCCTCGTTCGCGGCGTGTGCTGCCGGTGGCGGCCAGGCCAGCGCGAGCCCGCGAGCACAGCCACGAGGGGGCCAGCGAGTTGCTGGCCGGCATTGCCCGGCTGGCCAATGCCGGTGCCAGCGAGCAGGCCCg                                                                                                                     >  2:7429774/1‑150 (MQ=255)
                                   aGCCCAGTGTGCCGTTGCCTCGTTCGCGGCGTGTGCTGCCGGTGGCGGCCAGGCCAGCGCGAGCCCGCGAGCACAGCCACGAGGGGGCCAGCGAGTTGCTGGCCGGCATTGCCCGGCTGGCCAATGCCGGTGCCAGCGAGCAGGCCCGTa                                                                                                                   <  1:570940/150‑1 (MQ=255)
                                    gCCCAGTGTGCCGTTGCCTCGTTCGCGGCGTGTGCTGCCGGTGGCGGCCAGGCCAGCGCGAGCCCGCGAGCACAGCCACGAGGGGGCCAGCGAGTTGCTGGCCGGCATTGCCCGGCTGGCCAATGCCGGTGCCAGCGAGCAGGCCCGTAg                                                                                                                  >  2:2549788/1‑150 (MQ=255)
                                     cccAGTGTGCCGTTGCCTCGTTCGCGGCGTGTGCTGCCGGTGGCGGCCAGGCCAGCGCGAGCCCGCGAGCACAGCCACGAGGGGGCCAGCGAGTTGCTGGCCGGCATTGCCCGGCTGGCCAATGCCGGTGCCAGCGAGCAGGCCCGTAGc                                                                                                                 >  1:8794496/1‑150 (MQ=255)
                                     cccAGTGTGCCGTTGCCTCGTTCGCGGCGTGTGCTGCCGGTGGCGGCCAGGCCAGCGCGAGCCCGCGAGCACAGCCACGAGGGGGCCAGCGAGTTGCTGGCCGGCATTGCCCGGCTGGCCAATGCCGGTGCCAGCGAGCAGGCCCGTAGc                                                                                                                 >  1:1656724/1‑150 (MQ=255)
                                      ccAGTGTGCCGTTGCCTCGTTCGCGGCGTGTGCTGCCGGTGGCGGCCAGGCCAGCGCGAGCCCGCGAGCACAGCCACGAGGGGGCCAGCGAGTTGCTGGCCGGCATTGCCCGGCTGGCCAATGCCGGTGCCAGCGAGCAGGCCCGTAGCg                                                                                                                >  2:4792351/1‑150 (MQ=255)
                                       cAGTGTGCCGTTGCCTCGTTCGCGGCGTGTGCTGCCGGTGGCGGCCAGGCCAGCGCGAGCCCGCGAGCACAGCCACGAGGGGGCCAGCGAGTTGCTGGCCGGCATTGCCCGGCTGGCCAATGCCGGTGCCAGCGAGCAGGCCCGTAGCGa                                                                                                               <  2:7833918/150‑1 (MQ=255)
                                        aGTGTGCCGTTGCCTCGTTCGCGGCGTGTGCTGCCGGTGGCGGCCAGGCCAGCGCGAGCCCGCGAGCACAGCCACGAGGGGGCCAGCGAGTTGCTGGCCGGCATTGCCCGGCTGGCCAATGCCGGTGCCAGCGAGCAGGCCCGTAGCGaa                                                                                                              <  1:4011399/150‑1 (MQ=255)
                                         gtgtGCCGTTGCCTCGTTCGCGGCGTGTGCTGCCGGTGGCGGCCAGGCCAGCGCGAGCCCGCGAGCACAGCCACGAGGGGGCCAGCGAGTTGCTGGCCGGCATTGCCCGGCTGGCCAATGCCGGTGCCAGCGAGCAGGCCCGTAGCGAAt                                                                                                             <  2:5206221/150‑1 (MQ=255)
                                         gtgtGCCGTTGCCTCGTTCGCGGCGTGTGCTGCCGGTGGCGGCCAGGCCAGCGCGAGCCCGCGAGCACAGCCACGAGGGGGCCAGCGAGTTGCTGGCCGGCATTGCCCGGCTGGCCAATGCCGGTGCCAGCGAGCAGGCCCGTAGCGAAt                                                                                                             >  1:4191964/1‑150 (MQ=255)
                                                    ccTCGTTCGCGGCGTGTGCTGCCGGTGGCGGCCAGGCCAGCGCGAGCCCGCGAGCACAGCCACGAGGGGGCCAGCGAGTTGCTGGCCGGCATTGCCCGGCTGGCCAATGCCGGTGCCAGCGAGCAGGCCCGTAGCGAATGCCAACGCTAc                                                                                                  <  2:883675/150‑1 (MQ=255)
                                                     cTCGTTCGCGGCGTGTGCTGCCGGTGGCGGCCAGGCCAGCGCGAGCCCGCGAGCACAGCCACGAGGGGGCCAGCGAGTTGCTGGCCGGCATTGCCCGGCTGGCCAATGCCGGTGCCAGCGAGCAGGCCCGTAGCGAATGCCAACGCTAcc                                                                                                 >  2:8026308/1‑150 (MQ=255)
                                                      tCGTTCGCGGCGTGTGCTGCCGGTGGCGGCCAGGCCAGCGCGAGCCCGCGAGCACAGCCACGAGGGGGCCAGCGAGTTGCTGGCCGGCATTGCCCGGCTGGCCAATGCCGGTGCCAGCGAGCAGGCCCGTAGCGAATGCCAACGCTACCt                                                                                                >  1:3365611/1‑150 (MQ=255)
                                                         ttCGCGGCGTGTGCTGCCGGTGGCGGCCAGGCCAGCGCGAGCCCGCGAGCACAGCCACGAGGGGGCCAGCGAGTTGCTGGCCGGCATTGCCCGGCTGGCCAATGCCGGTGCCAGCGAGCAGGCCCGTAGCGATTGCCAACGCTACCTGAg                                                                                             <  1:717268/150‑1 (MQ=255)
                                                          tCGCGGCGTGTGCTGCCGGTGGCGGCCAGGCCAGCGCGAGCCCGCGAGCACAGCCACGAGGGGGCCAGCGAGTTGCTGGCCGGCATTGCCCGGCTGGCCAATGCCGGTGCCAGCGAGCAGGCCCGTAGCGAATGCCAACGCTACCTGAGc                                                                                            <  2:8679589/150‑1 (MQ=255)
                                                          tCGCGGCGTGTGCTGCCGGTGGCGGCCAGGCCAGCGCGAGCCCGCGAGCACAGCCACGAGGGGGCCAGCGAGTTGCTGGCCGGCATTGCCCGGCTGGCCAATGCCGGTGCCAGCGAGCAGGCCCGTAGCGAATGCCAACGCTACCTGAGc                                                                                            <  1:3515899/150‑1 (MQ=255)
                                                          tCGCGGCGTGTGCTGCCGGTGGCGGCCAGGCCAGCGCGAGCCCGCGAGCACAGCCACGAGGGGGCCAGCGAGTTGCTGGCCGGCATTGCCCGGCTGGCCAATGCCGGTGCCAGCGAGCAGGCCCGTAGCGAATGCCAACGCTACCTGAGc                                                                                            <  2:4497314/150‑1 (MQ=255)
                                                            gcggcgTGTGCTGCCGGTGGCGGCCAGGCCAGCGCGAGCCCGCGAGCACAGCCACGAGGGGGCCAGCGAGTTGCTGGCCGGCATTGCCCGGCTGGCCAATGCCGGTGCCAGCGAGCAGGCCCGTAGCGAATGCCAACGCTACCTGAGCCa                                                                                          <  2:6565498/150‑1 (MQ=255)
                                                            gcggcgTGTGCTGCCGGTGGCGGCCAGGCCAGCGCGAGCCCGCGAGCACAGCCACGAGGGGGCCAGCGAGTTGCTGGCCGGCATTGCCCGGCTGGCCAATGCCGGTGCCAGCGAGCAGGCCCGTAGCGAATGCCAACGCTACCTGAGCCa                                                                                          <  1:8363069/150‑1 (MQ=255)
                                                               gcgTGTGCTGCCGGTGGCGGCCAGGCCAGCGCGAGCCCGCGAGCACAGCCACGAGGGGGCCAGCGAGTTGCTGGCCGGCATTGCCCGGCTGGCCAATGCCGGTGCCAGCGAGCAGGCCCGTAGCGAATGCCAACGCTACCTGAGCCAGTa                                                                                       <  2:7060240/150‑1 (MQ=255)
                                                                  tgtgCTGCCGGTGGCGGCCAGGCCAGCGCGAGCCCGCGAGCACAGCCACGAGGGGGCCAGCGAGTTGCTGGCCGGCATTGCCCGGCTGGCCAATGCCGGTGCCAGCGAGCAGGCCCGTAGCGAATGCCAACGCTACCTGAGCCAGTAccc                                                                                    >  2:2914541/1‑150 (MQ=255)
                                                                       tgcCGGTGGCGGCCAGGCCAGCGCGAGCCCGCGAGCACAGCCACGAGGGGGCCAGCGAGTTGCTGGCCGGCATTGCCCGGCTGGCCAATGCCGGTGCCAGCGAGCAGGCCCGTAGCGAATGCCAACGCTACCTGAGCCAGTACCCGCCGt                                                                               <  2:1685850/150‑1 (MQ=255)
                                                                         ccGGTGGCGGCCAGGCCAGCGCGAGCCCGCGAGCACAGCCACGAGGGGGCCAGCGAGTTGCTGGCCGGCATTGCCCGGCTGGCCAATGCCGGTGCCAGCGAGCAGGCCCGTAGCGAATGCCAACGCTACCTGAGCCAGTACCCGCCGTCg                                                                             <  1:2865835/150‑1 (MQ=255)
                                                                          cGGTGGCGGCCAGGCCAGCGCGAGCCCGCGAGCACAGCCACGAGGGGGCCAGCGAGTTGCTGGCCGGCATTGCCCGGCTGGCCAATGCCGGTGCCAGCGAGCAGGCCCGTAGCGAATGCCAACGCTACCTGAGCCAGTACCCGCCGTCgg                                                                            <  1:8599054/150‑1 (MQ=255)
                                                                            gTGGCGGCCAGGCCAGCGCGAGCCCGCGAGCACAGCCACGAGGGGGCCAGCGAGTTGCTGGCCGGCATTGCCCGGCTGGCCAATGCCGGTGCCAGCGAGCAGGCCCGTAGCGAATGCCAACGCTACCTGAGCCAGTACCCGCCGTCGGca                                                                          <  1:1575102/150‑1 (MQ=255)
                                                                             tGGCGGCCAGGCCAGCGCGAGCCCGCGAGCACAGCCACGAGGGGGCCAGCGAGTTGCTGGCCGGCATTGCCCGGCTGGCCAATGCCGGTGCCAGCGAGCAGGCCCGTAGCGAATGCCAACGCTACCTGAGCCAGTACCCGCCGTCGGcac                                                                         >  1:2211042/1‑150 (MQ=255)
                                                                               gcggcCAGGCCAGCGCGAGCCCGCGAGCACAGCCACGAGGGGGCCAGCGAGTTGCTGGCCGGCATTGCCCGGCTGGCCAATGCCGGTGCCAGCGAGCAGGCCCGTAGCGAATGCCAACGCTCCCTGAGCCAGTACCCGCCGTCGGCACAg                                                                       >  1:2504580/1‑150 (MQ=255)
                                                                                 ggccaggccaGCGCGAGCCCGCGAGCACAGCCACGAGGGGGCCAGCGAGTTGCTGGCCGGCATTGCCCGGCTGGCCAATGCCGGTGCCAGCGAGCAGGCCCGTAGCGAATGCCAACGCTACCTGAGCCAGTACCCGCCGTCGGCACAGgt                                                                     <  1:4114603/150‑1 (MQ=255)
                                                                                  gccaggccagCGCGAGCCCGCGAGCACAGCCACGAGGGGGCCAGCGAGTTGCTGGCCGGCATTGCCCGGCTGGCCAATGCCGGTGCCAGCGAGCAGGCCCGTAGCGAATGCCAACGCTACCTGAGCCAGTACCCGCCGTCGGCACAGgtg                                                                    >  2:8221028/1‑150 (MQ=255)
                                                                                     aggccagCGCGAGCCCGCGAGCACAGCCACGAGGGGGCCAGCGAGTTGCTGGCCGGCATTGCCCGGCTGGCCAATGCCGGTGCCAGCGAGCAGGCCCGTAGCGAATGCCAACGCTACCTGAGCCAGTACCCGCCGTCGGCACAGGTGtac                                                                 >  2:6105209/1‑150 (MQ=255)
                                                                                       gccagCGCGAGCCCGCGAGCACAGCCACGAGGGGGCCAGCGAGTTGCTGGCCGGCATTGCCCGGCTGGCCAATGCCGGTGCCAGCGAGCAGGCCCGTAGCGAATGCCAACGCTACCTGAGCCAGTACCCGCCGTCGGCACAGGTGtacta                                                               <  1:4775646/150‑1 (MQ=255)
                                                                                       gccagCGCGAGCCCGCGAGCACAGCCACGAGGGGGCCAGCGAGTTGCTGGCCGGCATTGCCCGGCTGGCCAATGCCGGTGCCAGCGAGCAGGCCCGTAGCGAATGCCAACGCTACCTGAGCCAGTACCCGCCGTCGGCACAGGTGtacta                                                               >  1:752089/1‑150 (MQ=255)
                                                                                           gcgcGAGCCCGCGAGCACAGCCACGAGGGGGCCAGCGAGTTGCTGGCCGGCATTGCCCGGCTGGCCAATGCCGGTGCCAGCGAGCAGGCCCGTAGCGAATGCCAACGCTACCTGAGCCAGTACCCGCCGTCGGCACAGGTGTACTActgg                                                           <  1:2279281/150‑1 (MQ=255)
                                                                                           gcgcGAGCCCGCGAGCACAGCCACGAGGGGGCCAGCGAGTTGCTGGCCGGCATTGCCCGGCTGGCCAATGCCGGTGCCAGCGAGCAGGCCCGTAGCGAATGCCAACGCTACCTGAGCCAGTACCCGCCGTCGGCACAGGTGTACTActgg                                                           >  2:3584002/1‑150 (MQ=255)
                                                                                               gAGCCCGCGAGCACAGCCACGAGGGGGCCAGCGAGTTGCTGGCCGGCATTGCCCGGCTGGCCAATGCCGGTGCCAGCGAGCAGGCCCGTAGCGAATGCCAACGCTACCTGAGCCAGTACCCGCCGTCGGCACAGGTGTACTActggctgg                                                       <  2:5397941/150‑1 (MQ=255)
                                                                                                 gCCCGCGAGCACAGCCACGAGGGGGCCAGCGAGTTGCTGGCCGGCATTGCCCGGCTGGCCAATGCCGGTGCCAGCGAGCAGGCCCGTAGCGAATGCCAACGCTACCTGAGCCAGTACCCGCCGTCGGCACAGGTGTACTACTGGCTGGGt                                                     >  2:8558683/1‑150 (MQ=255)
                                                                                                 gCCCGCGAGCACAGCCACGAGGGGGCCAGCGAGTTGCTGGCCGGCATTGCCCGGCTGGCCAATGCCGGTGCCAGCGAGCAGGCCCGTAGCGAATGCCAACGCTACCTGAGCCAGTACCCGCCGTCGGCACAGGTGTACTACTGGCTGGGt                                                     <  2:5589009/150‑1 (MQ=255)
                                                                                                   ccGCGAGCACAGCCACGAGGGGGCCAGCGAGTTGCTGGCCGGCATTGCCCGGCTGGCCAATGCCGGTGCCAGCGAGCAGGCCCGTAGCGAATGCCAACGCTACCTGAGCCAGTACCCGCCGTCGGCACAGGTGTACTACTGGCTGGGttt                                                   >  1:4913006/1‑150 (MQ=255)
                                                                                                    cgcgAGCACAGCCACGAGGGGGCCAGCGAGTTGCTGGCCGGCATTGCCCGGCTGGCCAATGCCGGTGCCAGCGAGCAGGCCCGTAGCGAATGCCAACGCTACCTGAGCCAGTACCCGCCGTCGGCACAGGTGTACTACTGGCTGGGTTTg                                                  <  2:7360317/150‑1 (MQ=255)
                                                                                                        aGCACAGCCACGAGGGGGCCAGCGAGTTGCTGGCCGGCATTGCCCGGCTGGCCAATGCCGGTGCCAGCGAGCAGGCCCGTAGCGAATGCCAACGCTACCTGAGCCAGTACCCGCCGTCGGCACAGGTGTACTACTGGCTGGGTTTGCTCa                                              <  1:2740176/150‑1 (MQ=255)
                                                                                                          cacaGCCACGAGGGGGCCAGCGAGTTGCTGGCCGGCATTGCCCGGCTGGCCAATGCCGGTGCCAGCGAGCAGGCCCGTAGCGAATGCCAACGCTACCTGAGCCAGTACCCGCCGTCGGCACAGGTGTACTACTGGCTGGGTTTGCTCAGt                                            <  1:97299/150‑1 (MQ=255)
                                                                                                           acaGCCACGAGGGGGCCAGCGAGTTGCTGGCCGGCATTGCCCGGCTGGCCAATGCCGGTGCCAGCGAGCAGGCCCGTAGCGAATGCCAACGCTACCTGAGCCAGTACCCGCCGTCGGCACAGGTGTACTACTGGCTGGGTTTGCTCAGTg                                           >  1:1079100/1‑150 (MQ=255)
                                                                                                              gCCACGAGGGGGCCAGCGAGTTGCTGGCCGGCATTGCCCGGCTGGCCAATGCCGGTGCCAGCGAGCAGGCCCGTAGCGAATGCCAACGCTACCTGAGCCAGTACCCGCCGTCGGCACAGGTGTACTACTGGCTGGGTTTGCTCAGTGATa                                        <  2:6281760/150‑1 (MQ=255)
                                                                                                               ccACGAGGGGGCCAGCGAGTTGCTGGCCGGCATTGCCCGGCTGGCCAATGCCGGTGCCAGCGAGCAGGCCCGTAGCGAATGCCAACGCTACCTGAGCCAGTACCCGCCGTCGGCACAGGTGTACTACTGGCTGGGTTTGCTCAGTGATAc                                       <  2:1247061/150‑1 (MQ=255)
                                                                                                                   gAGGGGGCCAGCGAGTTGCTGGCCGGCATTGCCCGGCTGGCCAATGCCGGTGCCAGCGAGCAGGCCCGTAGCGAATGCCAACGCTACCTGAGCCAGTACCCGCCGTCGGCACAGGTGTACTACTGGCTGGGTTTGCTCAGTGATACCGaa                                   <  1:3193951/150‑1 (MQ=255)
                                                                                                                     gggggCCAGCGAGTTGCTGGCCGGCATTGCCCGGCTGGCCAATGCCGGTGCCAGCGAGCAGGCCCGTAGCGAATGCCAACGCTACCCGAGCCAGTACCCGCCGTCGGCACAGGTGTACTACTGGCTGGGTTTGCTCAGTGATACCGAAgg                                 <  1:6209585/150‑1 (MQ=255)
                                                                                                                           cAGCGAGTTGCTGGCCGGCATTGCCCGGCTGGCCAATGCCGGTGCCAGCGAGCAGGCCCGTAGCGAATGCCAACGCTACCTGAGCCAGTACCCGCCGTCGGCACAGGTGTACTACTGGCTGGGTTTGCTCAGTGATACCGAAGGCGATgc                           >  2:4372762/1‑150 (MQ=255)
                                                                                                                           cAGCGAGTTGCTGGCCGGCATTGCCCGGCTGGCCAATGCCGGTGCCAGCGAGCAGGCCCGTAGCGAATGCCAACGCTACCTGAGCCAGTACCCGCCGTCGGCACAGGTGTACTACTGGCTGGGTTTGCTCAGTGATACCGAAGGCGATgc                           >  1:6078767/1‑150 (MQ=255)
                                                                                                                             gCGAGTTGCTGGCCGGCATTGCCCGGCTGGCCAATGCCGGTGCCAGCGAGCAGGCCCGTAGCGAATGCCAACGCTACCTGAGCCAGTACCCGCCGTCGGCACAGGTGTACTACTGGCTGGGTTTGCTCAGTGATACCGAAGGCGATgcgc                         >  2:5681832/1‑150 (MQ=255)
                                                                                                                                  ttGCTGGCCGGCATTGCCCGGCTGGCCAATGCCGGTGCCAGCGAGCAGGCCCGTAGCGAATGCCAACGCTACCTGAGCCAGTACCCGCCGTCGGCACAGGTGTACTACTGGCTGGGTTTGCTCAGTGATACCGAAGGCGATGCGcagcag                    >  2:2430622/1‑150 (MQ=255)
                                                                                                                                    tggcaccGGCATTGGCCAGCCGGGCAATGCCGGCGCCAGCGAGCAGGCCCGTAGCGAATGCCAACGCTACCTGAGCCAGTACCCGCCGTCGGCACAGGTGTACTACTGGCTGGGTTTGCTCAGTGATACCGAAGGCGATGCGCAGCAGgc                  >  2:303562/6‑150 (MQ=255)
                                                                                                                                    gCTGGCCGGCATTGCCCGGCTGGCCAATGCCGGTGCCAGCGAGCAGGCCCGTAGCGAATGCCAACGCTACCTGAGCCAGTACCCGCCGTCGGCACAGGTGTACTACTGGCTGGGTTTGCTCAGTGATACCGAAGGCGATGCGCAGCAGgc                  >  2:5098761/1‑150 (MQ=255)
                                                                                                                                    gCTGGCCGGCATTGCCCGGCTGGCCAATGCCGGTGCCAGCGAGCAGGCCCGTAGCGAATGCCAACGCTACCTGAGCCAGTACCCGCCGTCGGCACAGGTGTACTACTGGCTGGGTTTGCTCAGTGATACCGAAGGCGATGCGCAGCAGgc                  >  1:3909542/1‑150 (MQ=255)
                                                                                                                                        gCCGGCATTGCCCGGCTGGCCAATGCCGGTGCCAGCGAGCAGGCCCGTAGCGAATGGCAACGCTACCTGAGCCAGTACCCGCCGTCGGCACAGGTGTACTACTGGCTGGGTTTGCTCAGTGATACCGAAGGCGATGCGCAGCAGGCGCTg              >  1:2624416/1‑150 (MQ=255)
                                                                                                                                        gCCGGCATTGCCCGGCTGGCCAATGCCGGTGCCAGCGAGCAGGCCCGTAGCGAATGCCAACGCTACCTGAGCCAGTACCCGCCGTCGGCACAGGTGTACTACTGGCTGGGTTTGCTCAGTGATACCGAAGGCGATGCGCAGCAGGCGCTg              >  2:7136569/1‑150 (MQ=255)
                                                                                                                                         ccGGCATTGGCCAGCCGGGCAATGCCGGCGCCAGCGAGCAGGCCCGTAGCGAATGCCAACGCTACCTGAGCCAGTACCCGCCGTCGGCACAGGTGTACTACTGGCTGGGTTTGCTCAGTGATACCGAAGGCGATGCGCAGCAGGCGCTGa             <  1:303562/150‑1 (MQ=255)
                                                                                                                                         ccGGCATTGCCCGGCTGGCCAATGCCGGTGCCAGCGAGCAGGCCCGTAGCGAATGCCAACGCTACCTGAGCCAGTACCCGCCGTCGGCACAGGTGTACTACTGGCTGGGTTTGCTCAGTGATACCGAAGGCGATGCGCAGCAGGCGCTGa             <  2:6774360/150‑1 (MQ=255)
                                                                                                                                         ccGGCATTGCCCGGCTGGCCAATGCCGGTGCCAGCGAGCAGGCCCGTAGCGAATGCCAACGCTACCTGAGCCAGTACCCGCCGTCGGCACAGGTGTACTACTGGCTGGGTTTGCTCAGTGATACCGAAGGCGATGCGCAGCAGGCGCTGa             <  1:5445915/150‑1 (MQ=255)
                                                                                                                                          cGGCATTGCCCGGCTGGCCAATGCCGGTGCCAGCGAGCAGGCCCGTAGCGAATGCCAACGCTACCTGAGCCAGTACCCGCCGTCGGCACAGGTGTACTACTGGCTGGGTTTGCTCAGTGATACCGAAGGCGATGCGCAGCAGGCGCTGAg            >  1:6885307/1‑150 (MQ=255)
                                                                                                                                          cGGCATTGCCCGGCTGGCCAATGCCGGTGCCAGCGAGCAGGCCCGTAGCGAATGCCAACGCTACCTGAGCCAGTACCCGCCGTCGGCACAGGTGTACTACTGGCTGGGTTTGCTCAGTGATACCGAAGGCGATGCGCAGCAGGCGCTGAg            >  1:4590957/1‑150 (MQ=255)
                                                                                                                                          cGGCATTGCCCGGCTGGCCAATGCCGGTGCCAGCGAGCAGGCCCGTAGCGAATGCCAACGCTACCTGAGCAAGTACCCGCCGTCGGCACAGGTGTACTACTGGCTGGGTTTGCTCAGTGATACCGAAGGCGATGCGCAGCAGGCGCTGAg            >  1:6471379/1‑150 (MQ=255)
                                                                                                                                           ggCATTGCCCGGCTGGCCAATGCCGGTGCCAGCGAGCAGGCCCGTAGCGAATGCCAACGCTACCTGAGCCAGTACCCGCCGTCGGCACAGGTGTACTACTGGCTGGGTTTGCTCAGTGATACCGAAGGCGATGCGCAGCAGGCGCTGAGc           >  1:6609919/1‑150 (MQ=255)
                                                                                                                                                 gCCCGGCTGGCCAATGCCGGTGCCAGCGAGCAGGCCCGTAGCGAATGCCAACGCTACCTGAGCCAGTACCCGCCGTCGGCACAGGTGTACTACTGGCTGGGTTTGCTCAGTGATACCGAAGGCGATGCGCAGCAGGCGCTGAGCCACTAt     >  2:5276957/1‑150 (MQ=255)
                                                                                                                                                  cccGGCTGGCCAATGCCGGTGCCAGCGAGCAGGCCCGTAGCGAATGCCAACGCTACCTGAGCCAGTACCCGCCGTCGGCACAGGTGTACTACTGGCTGGGTTTGCTCAGTGATACCGAAGGCGATGCGCAGCAGGCGCTGAGCCACTATc    >  2:1134348/1‑150 (MQ=255)
                                                                                                                                                  cccGGCTGGCCAAGGCCGGTGCCAGCGAGCAGGCCCGTAGCGAATGCCAACGCTACCTGAGCCAGTACCCGCCGTCGGCACAGGTGTACTACTGGCTGGGTTTGCTCAGTGATACCGAAGGCGATGCGCAGCAGGCGCTGAGCCACTATc    <  1:3697747/150‑1 (MQ=255)
                                                                                                                                                   ccGGCTGGCCAATGCCGGTGCCAGCGAGCAGGCCCGTAGCGAATGCCAACGCTACCTGAGCCAGTACCCGCCGTCGGCACAGGTGTACTACTGGCTGGGTTTGCTCAGTGATACCGAAGGCGATGCGCAGCAGGCGCTGAGCCACTATCg   <  1:4974153/150‑1 (MQ=255)
                                                                                                                                                   ccGGCTGGCCAATGCCGGTGCCAGCGAGCAGGCCCGTAGCGAATGCCAACGCTACCTGAGCCAGTACCCGCCGTCGGCACAGGTGTACTACTGGCTGGGTTTGCTCAGTGATACCGAAGGCGATGCGCAGCAGGCGCTGAGCCACTATCg   >  1:4379590/1‑150 (MQ=255)
                                                                                                                                                   ccGGCTGGCCAATGCCGGTGCCAGCGAGCAGGCCCGTAGCGAATGCCAACGCTACCTGAGCCAGTACCCGCCGTCGGCACAGGTGTACTACTGGCTGGGTTTGCTCAGTGATACCGAAGGCGATGCGCAGCAGGCGCTGAGCCACTATCg   >  2:1391237/1‑150 (MQ=255)
                                                                                                                                                     |                                                                                                                                                    
TGCGTTCACGACGTTGCCCGCACCGGTTTATCCGCAGCCCAGTGTGCCGTTGCCTCGTTCGCGGCGTGTGCTGCCGGTGGCGGCCAGGCCAGCGCGAGCCCGCGAGCACAGCCACGAGGGGGCCAGCGAGTTGCTGGCCGGCATTGCCCGGCTGGCCAATGCCGGTGCCAGCGAGCAGGCCCGTAGCGAATGCCAACGCTACCTGAGCCAGTACCCGCCGTCGGCACAGGTGTACTACTGGCTGGGTTTGCTCAGTGATACCGAAGGCGATGCGCAGCAGGCGCTGAGCCACTATCGC  >  NC_002947/1699865‑1700162

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: