Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 1,700,011 C→G 6.1% A329G (GCC→GGC)  wspC → biofilm formation methyltransferase WspC

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0029471,700,0110CG 6.1% 1606.6 / 65.6 475A329G (GCC→GGC) wspCbiofilm formation methyltransferase WspC
Reads supporting (aligned to +/- strand):  ref base C (193/253);  new base G (15/14);  total (208/267)
Fisher's exact test for biased strand distribution p-value = 4.41e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.91e-01

Only 100 of 484 total aligned reads displayed.

GCGTGCGTTCACGACGTTGCCCGCACCGGTTTATCCGCAGCCCAGTGTGCCGTTGCCTCGTTCGCGGCGTGTGCTGCCGGTGGCGGCCAGGCCAGCGCGAGCCCGCGAGCACAGCCACGAGGGGGCCAGCGAGTTGCTGG‑‑CCGGCATTGCCCGGCTGGCCAATGCCGGTGCCAGCGAGCAGGCCCGTAGCGAATGCCAACGCTACCTGAGCCAGTACCCGCCGTCGGCACAGGTGTACTACTGGCTGGGTTTGCTCAGTGATACCGAAGGCGATGCGCAGCAGGCGCTGAGCCACTATC  >  NC_002947/1699862‑1700160
                                                                                                                                                       |                                                                                                                                                     
  gtgcgtTCACGACGTTGCCCGCACCGGTTGATCCGCAGCCCAGTGTGCCGTTGCCTCGTTCGCGACGTGTGCTGCCGGTGGCGGCCAGGCCAGCGCGAGCCCGCGAGCACAGCCACGAGGGGGCCAGCGAGTTGCTGG‑‑CCGGCATTGccc                                                                                                                                                     <  2:6049844/150‑1 (MQ=255)
      gtTCACGACGTTGCCCGCACCGGTTTATCCGCAGCCCAGTGTGCCGTTGCCTCGTTCGCGGCGTGTGCTGCCGGTGGCGGCCAGGCCAGCGCGAGCCCGCGAGCACAGCCACGAGGGGGCCAGCGAGTTGCTGG‑‑CCGGCATTGCCCGGCt                                                                                                                                                 <  1:7325967/150‑1 (MQ=255)
      gtTCACGACGTTGCCCGCACCGGTTTATCCGCAGCCCAGTGTGCCGTTGCCTCGTTCGCGGCGTGTGCTGCCGGTGGCGGCCAGGCCAGCGCGAGCCCGCGAGCACAGCCACGAGGGGGCCAGCGAGTTGCTGG‑‑CCGGCATTGCCCGGCt                                                                                                                                                 >  1:7916905/1‑150 (MQ=255)
      gtTCACGACGTTGCCCGCACCGGTTTATCCGCAGCCCAGTGTGCCGTTGCCTCGTTCGCGGCGTGTGCTGCCGGTGGCGGCCAGGCCAGCGCGAGCCCGCGAGCACAGCCACGAGGGGGCCAGCGAGTTGCTGG‑‑CCGGCATTGCCCGGCt                                                                                                                                                 >  2:1599560/1‑150 (MQ=255)
       ttCACGACGTTGCCCGCACCGGTTTATCCGCAGCCCAGTGTGCCGTTGCCTCGTTCGCGGCGTGTGCTGCCGGTGGCGGCCAGGCCAGCGCGAGCCCGCGAGCACAGCCACGAGGGGGCCAGCGAGTTGCTGG‑‑CCGGCATTGCCCGGCTg                                                                                                                                                <  2:5140947/150‑1 (MQ=255)
       ttCACGACGTTGCCCGCACCGGTTTATCCGCAGCCCAGTGTGCCGTTGCCTCGTTCGCGGCGTGTGCTGCCGGTGGCGGCCAGGCCAGCGCGAGCCCGCGAGCACAGCCACGAGGGGGCCAGCGAGTTGCTGG‑‑CCGGCATTGCCCGGCTg                                                                                                                                                <  2:6106475/150‑1 (MQ=255)
          acgacgTTGCCCGCACCGGTTTATCCGCAGCCCAGTGTGCCGTTGCCTCGTTCGCGGCGTGTGCTGCCGGTGGCGGCCAGGCCAGCGCGAGCCCGCGAGCACAGCCACGAGGGGGCCAGCGAGTTGCTGG‑‑CCGGCATTGCCCGGCTGGcc                                                                                                                                             >  1:6004827/1‑150 (MQ=255)
           cgacgTTGCCCGCACCGGTTTATCCGCAGCCCAGTGTGCCGTTGCCTCGTTCGCGGCGTGTGCTGCCGGTGGCGGCCAGGCCAGCGCGAGCCCGCGAGCACAGCCACGAGGGGGCCAGCGAGTTGCTGG‑‑CCGGCATTGCCCGGCTGGCCa                                                                                                                                            <  1:5186911/150‑1 (MQ=255)
             acgTTGCCCGCACCGGTTTATCCGCAGCCCAGTGTGCCGTTGCCTCGTTCGCGGCGTGTGCTGCCGGTGGCGGCCAGGCCAGCGCGAGCCCGCGAGCACAGCCACGAGGGGGCCAGCGAGTTGCTGG‑‑CCGGCATTGCCCGGCTGGCCAAt                                                                                                                                          <  1:532073/150‑1 (MQ=255)
               gTTGCCCGCACCGGTTTATCCGCAGCCCAGTGTGCCGTTGCCTCGTTCGCGGCGTGTGCTGCCGGTGGCGGCCAGGCCAGCGCGAGCCCGCGAGCACAGCCACGAGGGGGCCAGCGAGTTGCTGG‑‑CCGGCATTGCCCGGCTGGCCAAAGc                                                                                                                                        >  1:1744355/1‑150 (MQ=255)
                 tGCCCGCACCGGTTTATCCGCCGCCCCGTGTGCCGTTGCCTCGTTCGCGGCGTGTGCTGCCGGTGGCGGCCAGGCCAGCGCGAGCCCGCGAGCACAGCCACGAGGGGGCCAGCGAGTTGCTGG‑‑CCGGCATTGCCCGGCTGGCCAATGCCg                                                                                                                                      <  2:5435604/150‑1 (MQ=255)
                 tGCCCGCACCGGTTTATCCGCAGCCCAGTGTGCCGTTGCCTCGTTCGCGGCGTGTGCTGCCGGTGGCGGCCAGGCCAGCGCGAGCCCGCGAGCACAGCCACGAGGGGGCCAGCGAGTTGCTGG‑‑CCGGCATTGCCCGGCTGGCCAATGCCg                                                                                                                                      <  1:5863423/150‑1 (MQ=255)
                 tGCCCGCACCGGTTTATCCGCAGCCCAGTGTGCCGTTGCCTCGTTCGCGGCGTGTGCTGCCGGTGGCGGCCAGGCCAGCGCGAGCCCGCGAGCACAGCCACGAGGGGGCCAGCGAGTTGCTGG‑‑CCGGCATTGCCCGGCTGGCCAATGCCg                                                                                                                                      <  2:5607514/150‑1 (MQ=255)
                    ccGCACCGGTTTATCCGCAGCCCAGTGTGCCGTTGCCTCGTTCGCGGCGTGTGCTGCCGGTGGCGGCCAGGCCAGCGCGAGCCCGCGAGCACAGCCACGAGGGGGCCAGCGAGTTGCTGG‑‑CCGGCATTGCCCGGCTGGCCAATGCCGGTg                                                                                                                                   >  1:1247061/1‑150 (MQ=255)
                     cGCACCGGTTTATCCGCAGCCCAGTGTGCCGTTGCCTCGTTCGCGGCGTGTGCTGCCGGTGGCGGCCAGGCCAGCGCGAGCCCGCGAGCACAGCCACGAGGGGGCCAGCGAGTTGCTGG‑‑CCGGCATTGCCCGGCTGGCCAATGCCGGtca                                                                                                                                  <  2:5499630/150‑3 (MQ=255)
                       cACCGGTTTATCCGCAGCCCAGTGTGCCGTTGCCTCGTTCGCGGCGTGTGCTGCCGGTGGCGGCCAGGCCAGCGCGAGCCCGCGAGCACAGCCACGAGGGGGCCAGCGAGTTGCTGG‑‑CCGGCATTGCCCGGCTGGCCAATGCCGGTGCCa                                                                                                                                <  1:7874692/150‑1 (MQ=255)
                       cACCGGTTTATCCGCAGCCCAGTGTGCCGTTGCCTCGTTCGCGGCGTGTGCTGCCGGTGGCGGCCAGGCCAGCGCGAGCCCGCGAGCACAGCCACGAGGGGGCCAGCGAGTTGCTGG‑‑CCGGCATTGCCCGGCTGGCCAATGCCGGTGCCa                                                                                                                                <  1:805913/150‑1 (MQ=255)
                           ggTTTATCCGCAGCCCAGTGTGCCGTTGCCTCGGTCGCGGCGTGTGCTGCCGGTGGCGGCCAGGCCAGCGCGAGCCCGCGAGCACAGCCACGAGGGGGCCAGCGAGTTGCTGG‑‑CCGGCATTGCCCGGCTGGCCAATGCCGGTGCCAGCGa                                                                                                                            <  2:8424593/150‑1 (MQ=255)
                            gTTTATCCGCAGCCCAGTGTGCCGTTGCCTCGTTCGCGGCGTGTGCTGCCGGTGGCGGCCAGGCCAGCGCGAGCCCGCGAGCACAGCCACGAGGGGGCCAGCGAGTCGCTGGCACCGGCATTGGCCAGCCGGGCAATGCCggccagcaa                                                                                                                              >  1:2274849/1‑142 (MQ=37)
                            gTTTATCCGCAGCCCAGTGTGCCGTTGCCTCGTTCGCGGCGTGTGCTGCCGGTGGCGGCCAGGCCAGCGCGAGCCCGCGAGCACAGCCACGAGGGGGCCAGCGAGTCGCTGGCACCGGCATTGGCCAGCCGGGCAATGCCggccagcaa                                                                                                                              <  2:2274849/149‑8 (MQ=37)
                             tttATCCGCAGCCCAGTGTGCCGTTGCCTCGTTCGCGGCGTGTGCTGCCGGTGGCGGCCAGGCCAGCGCGAGCCCGCGAGCACAGCCACGAGGGGGCCAGCGAGTTGCTGG‑‑CCGGCATTGCCCGGCTGGCCAATGCCGGTGCCAGCGAGc                                                                                                                          <  1:8441581/150‑1 (MQ=255)
                              ttATCCGCAGCCCAGTGTGCCGTTGCCTCGTTCGCGGCGTGTGCTGCCGGTGGCGGCCAGGCCAGCGCGAGCCCGCGAGCACAGCCACGAGGGGGCCAGCGAGTTGCTGG‑‑CCGGCATTGCCCGGCTGGCCAATGCCGGTGCCAGCGAGCa                                                                                                                         <  2:5823403/150‑1 (MQ=255)
                              ttATCCGCAGCCCAGTGTGCCGTTGCCTCGTTCGCGGCGTGTGCTGCCGGTGGCGGCCAGGCCAGCGCGAGCCCGCGAGCACAGCCACGAGGGGGCCAGCGAGTTGCTGG‑‑CCGGCATTGCCCGGCTGGCCAATGCCGGTGCCAGCGAGCa                                                                                                                         <  1:3353619/150‑1 (MQ=255)
                                aTCCGCAGCCCAGTGTGCCGTTGCCTCGTTCGCGGCGTGTGCTGCCGGTGGCGGCCAGGCCAGCGCGAGCCCGCGAGCACAGCCACGAGGGGGCCAGCGAGTTGCTGG‑‑CCGGCATTGCCCGGCTGGCCAATGCCGGTGCCAGCGAGCAgg                                                                                                                       <  1:1918929/150‑1 (MQ=255)
                                   cGCAGCCCAGTGTGCCGTTGCCTCGTTCGCGGCGTGTGCTGCCGGTGGCGGCCAGGCCAGCGCGAGCCCGCGAGCACAGCCACGAGGGGGCCAGCGAGTTGCTGG‑‑CCGGCATTGCCCGGCTGGCCAATGCCGGTGCCAGCGAGCAGGccc                                                                                                                    >  1:3689928/1‑150 (MQ=255)
                                   cGCAGCCCAGTGTGCCGTTGCCTCGTTCGCGGCGTGTGCTGCCGGTGGCGGCCAGGCCAGCGCGAGCCCGCGAGCACAGCCACGAGGGGGCCAGCGAGTTGCTGG‑‑CCGGCATTGCCCGGCTGGCCAATGCCGGTGCCAGCGAGCAGGccc                                                                                                                    <  1:1374235/150‑1 (MQ=255)
                                   cGCAGCCCAGTGTGCCGTTGCCTCGTTCGCGGCGTGTGCTGCCGGTGGCGGCCAGGCCAGCGCGAGCCCGCGAGCACAGCCACGAGGGGGCCAGCGAGTTGCTGG‑‑CCGGCATTGCCCGGCTGGCCAATGCCGGTGCCAGCGAGCAGGccc                                                                                                                    >  2:2079172/1‑150 (MQ=255)
                                    gCAGCCCAGTGTGCCGTTGCCTCGTTCGCGGCGTGTGCTGCCGGTGGCGGCCAGGCCAGCGCGAGCCCGCGAGCACAGCCACGAGGGGGCCAGCGAGTTGCTGG‑‑CCGGCATTGCCCGGCTGGCCAATGCCGGTGCCAGCGAGCAGGCCCg                                                                                                                   >  2:3154522/1‑150 (MQ=255)
                                        cccAGTGTGCCGTTGCCTCGTTCGCGGCGTGTGCTGCCGGTGGCGGCCAGGCCAGCGCGAGCCCGCGAGCACAGCCACGAGGGGGCCAGCGAGTTGCTGG‑‑CCGGCATTGCCCGGCTGGCCAATGCCGGTGCCAGCGAGCAGGCCCGTAGc                                                                                                               >  1:1656724/1‑150 (MQ=255)
                                        cccAGTGTGCCGTTGCCTCGTTCGCGGCGTGTGCTGCCGGTGGCGGCCAGGCCAGCGCGAGCCCGCGAGCACAGCCACGAGGGGGCCAGCGAGTTGCTGG‑‑CCGGCATTGCCCGGCTGGCCAATGCCGGTGCCAGCGAGCAGGCCCGTAGc                                                                                                               >  1:8794496/1‑150 (MQ=255)
                                         ccAGTGTGCCGTTGCCTCGTTCGCGGCGTGTGCTGCCGGTGGCGGCCAGGCCAGCGCGAGCCCGCGAGCACAGCCACGAGGGGGCCAGCGAGTTGCTGG‑‑CCGGCATTGCCCGGCTGGCCAATGCCGGTGCCAGCGAGCAGGCCCGTAGCg                                                                                                              >  2:4792351/1‑150 (MQ=255)
                                         ccAGAGTGCCGTTGCCTCGTTCGCGGCGGGTGCTGCCGGTGGCGGCCAGGCCAGCGCGAGCCCGCGAGCACAGCCACGAGGGGGCCAGCGAGTTGCTGG‑‑CCGGCATTGCCCGGCTGGCCAATGCCGGTGCCAGCGAGCAGGCCCGTAGCg                                                                                                              <  2:256203/150‑1 (MQ=255)
                                          cAGTGTGCCGTTGCCTCGTTCGCGGCGTGTGCTGCCGGTGGCGGCCAGGCCAGCGCGAGCCCGCGAGCACAGCCACGAGGGGGCCAGCGAGTTGCTGG‑‑CCGGCATTGCCCGGCTGGCCAATGCCGGTGCCAGCGAGCAGGCCCGTAGCGa                                                                                                             >  2:2000649/1‑150 (MQ=255)
                                             tgtgCCGTTGCCTCGTTCGCGGCGTGTGCTGCCGGTGGCGGCCAGGCCAGCGCGAGCCCGCGAGCACAGCCACGAGGGGGCCAGCGAGTTGCTGG‑‑CCGGCATTGCCCGGCTGGCCAATGCCGGTGCCAGCGAGCAGGCCCGTAGCGAATg                                                                                                          >  2:3704693/1‑150 (MQ=255)
                                             tgtgCCGTTGCCTCGTTCGCGGCGTGTGCTGCCGGTGGCGGCCAGGCCAGCGCGAGCCCGCGAGCACAGCCACGAGGGGGCCAGCGAGTTGCTGG‑‑CCGGCATTGCCCGGCTGGCCAATGCCGGTGCCAGCGAGCAGGCCCGTAGCGAATg                                                                                                          >  2:419439/1‑150 (MQ=255)
                                              gtgCCGTTGCCTCGTTCGCGGCGTGTGCTGCCGGTGGCGGCCAGGCCAGCGCGAGCCCGCGAGCACAGCCACGAGGGGGCCAGCGAGTTGCTGG‑‑CCGGCATTGCCCGGCTGGCCAATGCCGGTGCCAGCGAGCAGGCCCGTAGCGAATGc                                                                                                         <  2:7459018/150‑1 (MQ=255)
                                                gCCGTTGCCTCGTTCGCGGCGTGTGCTGCCGGTGGCGGCCAGGCCAGCGCGAGCCCGCGAGCACAGCCACGAGGGGGCCAGCGAGTTGCTGG‑‑CCGGCATTGCCCGGCTGGCCAATGCCGGTGCCAGCGAGCAGGCCCGTAGCGAATGCCa                                                                                                       <  2:6444260/150‑1 (MQ=255)
                                                 ccGTTGCCTCGTTCGCGGCGTGTGCTGCCGGTGGCGGCCAGGCCAGCGCGAGCCCGCGAGCACAGCCACGAGGGGGCCAGCGAGTTGCTGG‑‑CCGGCATTGCCCGGCTGGCCAATGCCGGTGCCAGCGAGCAGGCCCGTAGCGAATGCCaa                                                                                                      <  2:6590416/150‑1 (MQ=255)
                                                     tGCCTCGTTCGCGGCGTGTGCTGCCGGTGGCGGCCCGGCCCGCGCGAGCCCGCGAGCACAGCCACGAGGGGGCCAGCGAGTTGCTGG‑‑CCGGCATTGCCCGGCTGGCCAATGCCGGTGCCAGCGAGCAGGCCCGTAGCGAATGCCAACGCt                                                                                                  >  1:3021441/1‑150 (MQ=255)
                                                       ccTCGTTCGCGGCGTGTGCTGCCGGTGGCGGCCAGGCCAGCGCGAGCCCGCGAGCACAGCCACGAGGGGGCCAGCGAGTTGCTGG‑‑CCGGCATTGCCCGGCTGGCCAATGCCGGTGCCAGCGAGCAGGCCCGTAGCGAATGCCAACGCTAc                                                                                                <  2:883675/150‑1 (MQ=255)
                                                        cTCGTTCGCGGCGTGTGCTGCCGGTGGCGGCCAGGCCAGCGCGAGCCCGCGAGCACAGCCACGAGGGGGCCAGCGAGTTGCTGG‑‑CCGGCATTGCCCGGCTGGCCAATGCCGGTGCCAGCGAGCAGGCCCGTAGCGAATGCCAACGCTAcc                                                                                               >  2:1414962/1‑150 (MQ=255)
                                                        cTCGTTCGCGGCGTGTGCTGCCGGTGGCGGCCAGGCCAGCGCGAGCCCGCGAGCACAGCCACGAGGGGGCCAGCGAGTTGCTGG‑‑CCGGCATTGCCCGGCTGGCCAATGCCGGTGCCAGCGAGCAGGCCCGTAGCGAATGCCAACGCTAcc                                                                                               >  1:2002212/1‑150 (MQ=255)
                                                        cTCGTTCGCGGCGTGTGCTGCCGGTGGCGGCCAGGCCAGCGCGAGCCCGCGAGCACAGCCACGAGGGGGCCAGCGAGTTGCTGG‑‑CCGGCATTGCCCGGCTGGCCAATGCCGGTGCCAGCGAGCAGGCCCGTAGCGAATGCCAACGCTAcc                                                                                               >  2:8026308/1‑150 (MQ=255)
                                                         tCGTTCGCGGCGTGTGCTGCCGGTGGCGGCCAGGCCAGCGCGAGCCCGCGAGCACAGCCACGAGGGGGCCAGCGAGTTGCTGG‑‑CCGGCATTGCCCGGCTGGCCAATGCCGGTGCCAGCGAGCAGGCCCGTAGCGAATGCCAACGCTACCt                                                                                              >  1:3482090/1‑150 (MQ=255)
                                                         tCGTTCGCGGCGTGTGCTGCCGGTGGCGGCCAGGCCAGCGCGAGCCCGCGAGCACAGCCACGAGGGGGCCAGCGAGTTGCTGG‑‑CCGGCATTGCCCGGCTGGCCAATGCCGGTGCCAGCGAGCAGGCCCGTAGCGAATGCCAACGCTACCt                                                                                              >  1:2861783/1‑150 (MQ=255)
                                                            ttCGCGGCGTGTGCTGCCGGTGGCGGCCAGGCCAGCGCGAGCCCGCGAGCACAGCCACGAGGGGGCCAGCGAGTTGCTGG‑‑CCGGCATTGCCCGGCTGGCCAATGCCGGTGCCAGCGAGCAGGCCCGTAGCGATTGCCAACGCTACCTGAg                                                                                           <  1:717268/150‑1 (MQ=255)
                                                             tCGCGGCGTGTGCTGCCGGTGGCGGCCAGGCCAGCGCGAGCCCGCGAGCACAGCCACGAGGGGGCCAGCGAGTTGCTGG‑‑CCGGCATTGCCCGGCTGGCCAATGCCGGTGCCAGCGAGCAGGCCCGTAGCGAATGCCAACGCTACCTGAGc                                                                                          <  2:4497314/150‑1 (MQ=255)
                                                             tCGCGGCGTGTGCTGCCGGTGGCGGCCAGGCCAGCGCGAGCCCGCGAGCACAGCCACGAGGGGGCCAGCGAGTTGCTGG‑‑CCGGCATTGCCCGGCTGGCCAATGCCGGTGCCAGCGAGCAGGCCCGTAGCGAATGCCAACGCTACCTGAGc                                                                                          <  2:8679589/150‑1 (MQ=255)
                                                               gcggcgTGTGCTGCCGGTGGCGGCCAGGCCAGCGCGAGCCCGCGAGCACAGCCACGAGGGGGCCAGCGAGTTGCTGG‑‑CCGGCATTGCCCGGCTGGCCAATGCCGGTGCCAGCGAGCAGGCCCGTAGCGAATGCCAACGCTACCTGAGCCa                                                                                        <  1:8363069/150‑1 (MQ=255)
                                                                  gcgTGTGCTGCCGGTGGCGGCCAGGCCAGCGCGAGCCCGCGAGCACAGCCACGAGGGGGCCAGCGAGTTGCTGG‑‑CCGGCATTGCCCGGCTGGCCAATGCCGGTGCCAGCGAGCAGGCCCGTAGCGAATGCCAACGCTACCTGAGCCAGTa                                                                                     <  2:7060240/150‑1 (MQ=255)
                                                                   cgTGTGCTGCCGGTGGCGGCCAGGCCAGCGCGAGCCCGCGAGCACAGCCACGAGGGGGCCAGCGAGTTGCTGG‑‑CCGGCATTGCCCGGCTGGCCAATGCCGGTGCCAGCGAGCAGGCCCGTAGCGAATGCCAACGCTACCTGAGCCAGTAc                                                                                    <  1:504661/150‑1 (MQ=255)
                                                                     tgtgCTGCCGGTGGCGGCCAGGCCAGCGCGAGCCCGCGAGCACAGCCACGAGGGGGCCAGCGAGTTGCTGG‑‑CCGGCATTGCCCGGCTGGCCAATGCCGGTGCCAGCGAGCAGGCCCGTAGCGAATGCCAACGCTACCTGAGCCAGTAccc                                                                                  <  1:4384023/150‑1 (MQ=255)
                                                                     tgtgCTGCCGGTGGCGGCCAGGCCAGCGCGAGCCCGCGAGCACAGCCACGAGGGGGCCAGCGAGTTGCTGG‑‑CCGGCATTGCCCGGCTGGCCAATGCCGGTGCCAGCGAGCAGGCCCGTAGCGAATGCCAACGCTACCTGAGCCAGTAccc                                                                                  >  2:2914814/1‑150 (MQ=255)
                                                                     tgtgCTGCCGGTGGCGGCCAGGCCAGCGCGAGCCCGCGAGCACAGCCACGAGGGGGCCAGCGAGTTGCTGG‑‑CCGGCATTGCCCGGCTGGCCAATGCCGGTGCCAGCGAGCAGGCCCGTAGCGAATGCCAACGCTACCTGAGCCAGTAccc                                                                                  <  2:2451182/150‑1 (MQ=255)
                                                                            ccGGTGGCGGCCAGGACAGCGCGAGCCCGCGAGCACAGCCACGAGGGGGCCAGCGAGTTGCTGG‑‑CCGGCATTGCCCGGCTGGCCTATGCCGGTGCCAGCGAGCAGGCCCGTATAGAATGCCAACGCTACCTGAGCCATTTCCCGCCGTCg                                                                           <  2:3936901/150‑1 (MQ=255)
                                                                             cGGTGGCGGCCAGGCCAGCGCGAGCCCGCGAGCACAGCCACGAGGGGGCCAGCGAGTTGCTGG‑‑CCGGCATTGCCCGGCTGGCCAATGCCGGTGCCAGCGAGCAGGCCCGTAGCGAATGCCAACGCTACCTGAGCCAGTACCCGCCGTCgg                                                                          <  1:8599054/150‑1 (MQ=255)
                                                                               gTGGCGGCCAGGCCAGCGCGAGCCCGCGAGCACAGCCACGAGGGGGCCAGCGAGTTGCTGG‑‑CCGGCATTGCCCGGCTGGCCAATGCCGGTGCCAGCGAGCAGGCCCGTAGCGAATGCCAACGCTACCTGAGCCAGTACCCGCCGTCGGca                                                                        <  1:1575102/150‑1 (MQ=255)
                                                                                  gcggcCAGGCCAGCGCGAGCCCGCGAGCACAGCCACGAGGGGGCCAGCGAGTTGCTGG‑‑CCGGCATTGCCCGGCTGGCCAATGCCGGTGCCAGCGAGCAGGCCCGTAGCGAATGCCAACGCTCCCTGAGCCAGTACCCGCCGTCGGCACAg                                                                     >  1:2504580/1‑150 (MQ=255)
                                                                                   tggctgtgccagCGCGAGCCCGCGAGCACAGCCACGAGGGGGCCAGCGAGTTGCTGG‑‑CCGGCATTGCCCGGCTGGCCAATGCCGGTGCCAGCGAGCAGGCCCGTAGCGAATGCCAACGCTACCTGAGCCAGTACCCGCCGTCGGCACAgg                                                                    >  2:1015771/8‑150 (MQ=255)
                                                                                     gccaggccagCGCGAGCCCGCGAGCACAGCCACGAGGGGGCCAGCGAGTTGCTGG‑‑CCGGCATTGCCCGGCTGGCCAATGCCGGTGCCAGCGAGCAGGCCCGTAGCGAATGCCAACGCTACCTGAGCCAGTACCCGCCGTCGGCACAGgtg                                                                  >  2:8221028/1‑150 (MQ=255)
                                                                                     gccaggccagCGCGAGCCCGCGAGCACAGCCACGAGGGGGCCAGCGAGTTGCTGG‑‑CCGGCATTGCCCGGCTGGCCAATGCCGGTGCCAGCGAGCAGGCCCGTAGCGAATGCCAACGCTACCTGAGCCAGTACCCGCCGTCGGCACAGgtg                                                                  >  1:4090285/1‑150 (MQ=255)
                                                                                      ccaggccagCGCGAGCCCGCGAGCACAGCCACGAGGGGGCCAGCGAGTTGCTGG‑‑CCGGCATTGCCCGGCTGGCCAATGCCGGTGCCAGCGAGCAGGCCCGTAGCGAATGCCAACGCTACCTGAGCCAGTACCCGCCGTCGGCACAGgtgt                                                                 <  2:1802643/150‑1 (MQ=255)
                                                                                       caggccagCGCGAGCCCGCGAGCACAGCCACGAGGGGGCCAGCGAGTTGCTGG‑‑CCGGCATTGCCCGGCTGGCCAATGCCGGTGCCAGCGAGCAGGCCCGTAGCGAATGCCAACGCTACCTGAGCCAGTACCCGCCGTCGGCACAGGTGta                                                                <  2:6789920/150‑1 (MQ=255)
                                                                                          gccagCGCGAGCCCGCGAGCACAGCCACGAGGGGGCCAGCGAGTTGCTGG‑‑CCGGCATTGCCCGGCTGGCCAATGCCGGTGCCAGCGAGCAGGCCCGTAGCGAATGCCAACGCTACCTGAGCCAGTACCCGCCGTCGGCACAGGTGtacta                                                             <  1:4775646/150‑1 (MQ=255)
                                                                                          gccagCGCGAGCCCGCGAGCACAGCCACGAGGGGGCCAGCGAGTTGCTGG‑‑CCGGCATTGCCCGGCTGGCCAATGCCGGTGCCAGCGAGCAGGCCCGTAGCGAATGCCAACGCTACCTGAGCCAGTACCCGCCGTCGGCACAGGTGtacta                                                             >  1:752089/1‑150 (MQ=255)
                                                                                           ccagCGCGAGCCCGCGAGCACAGCCACGAGGGGGCCAGCGAGTTGCTGG‑‑CCGGCATTGCCCGGCTGGCCAATGCCGGTGCCAGCGAGCAGGCCCGTAGCGAATGCCAACGCTACCTGAGCCAGTACCCGCCGTCGGCACAGGTGtactac                                                            <  1:2190226/150‑1 (MQ=255)
                                                                                           ccagCGCGAGCCCGCGAGCACAGCCACGAGGGGGCCAGCGAGTTGCTGG‑‑CCGGCATTGCCCGGCTGGCCAATGCCGGTGCCAGCGAGCAGGCCCGTAGCGAATGCCAACGCTACCTGAGCCAGTACCCGCCGTCGGCACAGGTGtactac                                                            >  2:5720347/1‑150 (MQ=255)
                                                                                                    gCCCGCGAGCACAGCCACGAGGGGGCCAGCGAGTTGCTGG‑‑CCGGCATTGCCCGGCTGGCCAATGCCGGTGCCAGCGAGCAGGCCCGTAGCGAATGCCAACGCTACCTGAGCCAGTACCCGCCGTCGGCACAGGTGTACTACTGGCTGGGt                                                   >  2:8558683/1‑150 (MQ=255)
                                                                                                      ccGCGAGCACAGCCACGAGGGGGCCAGCGAGTTGCTGG‑‑CCGGCATTGCCCGGCTGGCCAATGCCGGTGCCAGCGAGCAGGCCCGTAGCGAATGCCAACGCTACCTGAGCCAGTACCCGCCGTCGGCACAGGTGTACTACTGGCTGGGttt                                                 >  1:4913006/1‑150 (MQ=255)
                                                                                                      ccGCGAGCACAGCCACGAGGGGGCCAGCGAGTTGCTGG‑‑CCGGCATTGCCCGGCTGGCCAATGCCGGTGCCAGCGAGCAGGCCCGTAGCGAATGCCAACGCTACCTGAGCCAGTACCCGCCGTCGGCACAGGTGTACTACTGGCTGGGttt                                                 >  2:6281991/1‑150 (MQ=255)
                                                                                                       cgcgAGCACAGCCACGAGGGGGCCAGCGAGTTGCTGG‑‑CCGGCATTGCCCGGCTGGCCAATGCCGGTGCCAGCGAGCAGGCCCGTAGCGAATGCCAACGCTACCTGAGCCAGTACCCGCCGTCGGCACAGGTGTACTACTGGCTGGGTTTg                                                <  2:7360317/150‑1 (MQ=255)
                                                                                                           aGCACAGCCACGAGGGGGCCAGCGAGTTGCTGG‑‑CCGGCATTGCCCGGCTGGCCAATGCCGGTGCCAGCGAGCAGGCCCGTAGCGAATGCCAACGCTACCTGAGCCAGTACCCGCCGTCGGCACAGGTGTACTACTGGCTGGGTTTGCTCa                                            <  1:2740176/150‑1 (MQ=255)
                                                                                                            gCACAGCCACGAGGGGGCCAGCGAGTTGCTGG‑‑CCGGCATTGCCCGGCTGGCCAATGCCGGTGCCAGCGAGCAGGCCCGTAGCGAATGCCAACGCTACCTGAGCCAGTACCCGCCGTCGGCACAGGTGTACTACTGGCTGGGTTTGCTCAg                                           <  2:773752/150‑1 (MQ=255)
                                                                                                              acaGCCACGAGGGGGCCAGCGAGTTGCTGG‑‑CCGGCATTGCCCGGCTGGCCAATGCCGGTGCCAGCGAGCAGGCCCGTAGCGAATGCCAACGCTACCTGAGCCAGTACCCGCCGTCGGCACAGGTGTACTACTGGCTGGGTTTGCTCAGTg                                         >  1:1079100/1‑150 (MQ=255)
                                                                                                               caGCCACGAGGGGGCCAGCGAGTTGCTGG‑‑CCGGCATTGCCCGGCTGGCCAATGCCGGTGCCAGCGAGCAGGCCCGTAGCGAATGCCAACGCTACCTGAGCCAGTACCCGCCGTCGGCACAGGTGTACTACTGGCTGGGTTTGCTCAGTGa                                        <  1:3173623/150‑1 (MQ=255)
                                                                                                                 gCCACGAGGGGGCCAGCGAGTTGCTGG‑‑CCGGCATTGCCCGGCTGGCCAATGCCGGTGCCAGCGAGCAGGCCCGTAGCGAATGCCAACGCTACCTGAGCCAGTACCCGCCGTCGGCACAGGTGTACTACTGGCTGGGTTTGCTCAGTGATa                                      <  2:4971722/150‑1 (MQ=255)
                                                                                                                  ccACGAGGGGGCCAGCGAGTTGCTGG‑‑CCGGCATTGCCCGGCTGGCCAATGCCGGTGCCAGCGAGCAGGCCCGTAGCGAATGCCAACGCTACCTGAGCCAGTACCCGCCGTCGGCACAGGTGTACTACTGGCTGGGTTTGCTCAGTGATAc                                     <  2:1247061/150‑1 (MQ=255)
                                                                                                                      gAGGGGGCCAGCGAGTTGCTGG‑‑CCGGCATTGCCCGGCTGGCCAATGCCGGTGCCAGCGAGCAGGCCCGTAGCGAATGCCAACGCTACCTGAGCCAGTACCCGCCGTCGGCACAGGTGTACTACTGGCTGGGTTTGCTCAGTGATACCGaa                                 <  1:4008030/150‑1 (MQ=255)
                                                                                                                      gAGGGGGCCAGCGAGTTGCTGG‑‑CCGGCATTGCCCGGCTGGCCAATGCCGGTGCCAGCGAGCAGGCCCGTAGCGAATGCCAACGCTACCTGAGCCAGTACCCGCCGTCGGCACAGGTGTACTACTGGCTGGGTTTGCTCAGTGATACCGaa                                 <  2:4112961/150‑1 (MQ=255)
                                                                                                                        gggggCCAGCGAGTTGCTGG‑‑CCGGCATTGCCCGGCTGGCCAATGCCGGTGCCAGCGAGCAGGCCCGTAGCGAATGCCAACGCTACCCGAGCCAGTACCCGCCGTCGGCACAGGTGTACTACTGGCTGGGTTTGCTCAGTGATACCGAAgg                               <  1:6209585/150‑1 (MQ=255)
                                                                                                                           ggCCAGCGAGTTGCTGG‑‑CCGGCATTGCCCGGCTGGCCAATGCCGGTGCCAGCGAGCAGGCCCGTAGCGAATGCCAACGCTACCTGAGCCAGTACCCGCCGTCGGCACAGGTGTACTACTGGCTGGGTTTGCTCAGTGATACCGAAGGCGa                            <  2:5232404/150‑1 (MQ=255)
                                                                                                                            gCCAGCGAGTTGCTGG‑‑CCGGCATTGCCCGGCTGGCCAATGCCGGTGCCAGCGAGCAGGCCCGTAGCGAATGCCAACGCTACCTGAGCCAGTACCCGCCGTCGGCACAGGTGTACTACTGGCTGGGTTTGCTCAGTGATACCGAAGGCGAt                           <  2:2843612/150‑1 (MQ=255)
                                                                                                                              cAGCGAGTTGCTGG‑‑CCGGCATTGCCCGGCTGGCCAATGCCGGTGCCAGCGAGCAGGCCCGTAGCGAATGCCAACGCTACCTGAGCCAGTACCCGCCGTCGGCACAGGTGTACTACTGGCTGGGTTTGCTCAGTGATACCGAAGGCGATgc                         >  2:4372762/1‑150 (MQ=255)
                                                                                                                                      tGCTGG‑‑CCGGCATTGCCCGGCTGGCCAATGCCGGTGCCAGCGAGCAGGCCCGTAGCGAATGCCAACGCTACCTGAGCCAGTACCCGCCGTCGGCACAGGTGTACTACTGGCTGGGTTTGCTCAGTGATACCGAAGGCGATGCGCAGCAgg                 <  1:5464812/150‑1 (MQ=255)
                                                                                                                                       gCTGG‑‑CCGGCATTGCCCGGCTGGCCAATGCCGGTGCCAGCGAGCAGGCCCGTAGCGAATGCCAACGCTACCTGAGCCAGTACCCGCCGTCGGCACAGGTGTACTACTGGCTGGGTTTGCTCAGTGATACCGAAGGCGATGCGCAGCAGgc                >  2:6214122/1‑150 (MQ=255)
                                                                                                                                        cTGG‑‑CCGGCATTGCCCGGCTGGCCAATGCCGGTGCCAGCGAGCAGGCCCGTAGCGAATGCCAACGCTACCTGAGCCAGTACCCGCCGTCGGCACAGGTGTACTACTGGCTGGGTTTGCTCAGTGATACCGAAGGCGATGCGCAGCAGgcg               <  1:373733/150‑1 (MQ=255)
                                                                                                                                         tggcaccGGCATTGGCCAGCCGGGCAATGCCGGCGCCAGCGAGCAGGCCCGTAGCGAATGCCAACGCTACCTGAGCCAGTACCCGCCGTCGGCACAGGTGTACTACTGGCTGGGTTTGCTCAGTGATACCGAAGGCGATGCGCAGCAGgc                >  2:303562/6‑150 (MQ=255)
                                                                                                                                         tggcaccGGCATTGGCCAGCCGGGCAATGCCG‑‑GCCAGCGAGCAGGCCCGTAGCGAATGCCAACGCTACCTGAGCCAGTACCCGCCGTCGGCACAGGTGTACTACTGGCTGGGTTTGCTCAGTGATACCGAAGGCGATGTGCAGCAGgcgc              >  1:5543849/6‑150 (MQ=38)
                                                                                                                                         tGG‑‑CCGGCATTGCCCGGCTGGCCAATGCCGGTGCCAGCGAGCAGGCCCGTAGCGAATGCCAACGCTACCTGAGCCAGTACCCGCCGTCGGCACAGGTGTACTACTGGCTGGGTTTGCTCAGTGATACCGAAGGCGATGCGCAGCAGgcgc              >  1:4181397/1‑150 (MQ=255)
                                                                                                                                           gcaccGGCATTGGCCAGCCGGGCAATGCCGGTGCCAGCGAGCAGGCCCGTAGCGAATGCCAACGCTACCTGAGCCAGTACCCGCCGTCGGCACAGGTGTACTACTGGCTGGGTTTGCTCAGTGATACCGAAGGCGATGCGCAGCAGgcgc              >  2:5395823/4‑150 (MQ=255)
                                                                                                                                           g‑‑CCGGCATTGCCCGGCTGGCCAATGCCGGTGCCAGCGAGCAGGCCCGTAGCGAATGGCAACGCTACCTGAGCCAGTACCCGCCGTCGGCACAGGTGTACTACTGGCTGGGTTTGCTCAGTGATACCGAAGGCGATGCGCAGCAGGCGCTg            >  1:2624416/1‑150 (MQ=255)
                                                                                                                                           g‑‑CCGGCATTGCCCGGCTGGCCAATGCCGGTGCCAGCGAGCAGGCCCGTAGCGAATGCCAACGCTACCTGAGCCAGTACCCGCCGTCGGCACAGGTGTACTACTGGCTGGGTTTGCTCAGTGATACCGAAGGCGATGCGCAGCAGGCGCTg            >  2:7136569/1‑150 (MQ=255)
                                                                                                                                              ccGGCATTGCCCGGCTGGCCAATGCCGGTGCCAGCGAGCAGGCCCGTAGCGAATGCCAACGCTACCTGAGCCAGTACCCGCCGTCGGCACAGGTGTACTACTGGCTGGGTTTGCTCAGTGATACCGAAGGCGATGCGCAGCAGGCGCTGa           <  1:3869813/150‑1 (MQ=255)
                                                                                                                                               cGGCATTGCCCGGCTGGCCAATGCCGGTGCCAGCGAGCAGGCCCGTAGCGAATGCCAACGCTACCTGAGCCAGTACCCGCCGTCGGCACAGGTGTACTACTGGCTGGGTTTGCTCAGTGATACCGAAGGCGATGCGCAGCAGGCGCTGAg          >  1:4590978/1‑150 (MQ=255)
                                                                                                                                               cGGCATTGCCCGGCTGGCCAATGCCGGTGCCAGCGAGCAGGCCCGTAGCGAATGCCAACGCTACCTGAGCCAGTACCCGCCGTCGGCACAGGTGTACTACTGGCTGGGTTTGCTCAGTGATACCGAAGGCGATGCGCAGCAGGCGCTGAg          >  1:6885307/1‑150 (MQ=255)
                                                                                                                                               cGGCATTGCCCGGCTGGCCAATGCCGGTGCCAGCGAGCAGGCCCGTAGCGAATGCCAACGCTACCTGAGCAAGTACCCGCCGTCGGCACAGGTGTACTACTGGCTGGGTTTGCTCAGTGATACCGAAGGCGATGCGCAGCAGGCGCTGAg          >  1:6471379/1‑150 (MQ=255)
                                                                                                                                                ggCATTGGCCAGCCGGGCAATGCCGGCGCCAGCGAGCAGGCCCGTAGCGAATGCCAACGCTACCTGAGCCAGTACCCGCCGTCGGCACAGGTGTACTACTGGCTGGGTTTGCTCAGTGATACCGAAGGCGATGCGCAGCAGGCGCTGAGc         >  2:788957/1‑150 (MQ=255)
                                                                                                                                                ggCATTGCCCGGCTGGCCAATGCCGGTGCCAGCGAGCAGGCCCGTAGCGAATGCCAACGCTACCTGAGCCAGTACCCGCCGTCGGCACAGGTGTACTACTGGCTGGGTTTGCTCAGTGATACCGAAGGCGATGCGCAGCAGGCGCTGAGc         >  1:6609919/1‑150 (MQ=255)
                                                                                                                                                 gCATTGGCCAGCCGGGCAATGCCGGTGCCAGCGAGCAGGCCCGTAGCGAATGCCAACGCTACCTGAGCCAGTACCCGCCGTCGGCACAGGTGTACTATTGGCTGGGTTTGCTCAGTGATACCGAAGGCGATGCGCAGCAGGCGCTGAGcc        >  1:6086519/1‑150 (MQ=255)
                                                                                                                                                       cccGGCTGGCCAATGCCGGTGCCAGCGAGCAGGCCCGTAGCGAATGCCAACGCTACCTGAGCCAGTACCCGCCGTCGGCACAGGTGTACTACTGGCTGGGTTTGCTCAGTGATACCGAAGGCGATGCGCAGCAGGCGCTGAGCCACTATc  >  2:1134348/1‑150 (MQ=255)
                                                                                                                                                       |                                                                                                                                                     
GCGTGCGTTCACGACGTTGCCCGCACCGGTTTATCCGCAGCCCAGTGTGCCGTTGCCTCGTTCGCGGCGTGTGCTGCCGGTGGCGGCCAGGCCAGCGCGAGCCCGCGAGCACAGCCACGAGGGGGCCAGCGAGTTGCTGG‑‑CCGGCATTGCCCGGCTGGCCAATGCCGGTGCCAGCGAGCAGGCCCGTAGCGAATGCCAACGCTACCTGAGCCAGTACCCGCCGTCGGCACAGGTGTACTACTGGCTGGGTTTGCTCAGTGATACCGAAGGCGATGCGCAGCAGGCGCTGAGCCACTATC  >  NC_002947/1699862‑1700160

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: