Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 1,700,020 C→G 6.2% A332G (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,0200CG 6.2% 1578.0 / 65.9 468A332G (GCC→GGC) wspCbiofilm formation methyltransferase WspC
Reads supporting (aligned to +/- strand):  ref base C (197/242);  new base G (15/14);  total (212/256)
Fisher's exact test for biased strand distribution p-value = 5.64e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.78e-01

Only 100 of 475 total aligned reads displayed.

CACGACGTTGCCCGCACCGGTTTATCCGCAGCCCAGTGTGCCGTTGCCTCGTTCGCGGCGTGTGCTGCCGGTGGCGGCCAGGCCAGCGCGAGCCCGCGAGCACAGCCACGAGGGGGCCAGCGAGTTGCTGG‑‑CCGGCATTGCCCGGCTGGCCAATGCCGGTGCCAGCGAGCAGGCCCGTAGCGAATGCCAACGCTACCTGAGCCAGTACCCGCCGTCGGCACAGGTGTACTACTGGCTGGGTTTGCTCAGTGATACCGAAGGCGATGCGCAGCAGGCGCTGAGCCACTATCGCAAGGCCC  >  NC_002947/1699871‑1700169
                                                                                                                                                       |                                                                                                                                                     
  cgacgTTGCCCGCACCGGTTTATCCGCAGCCCAGTGTGCCGTTGCCTCGTTCGCGGCGTGTGCTGCCGGTGGCGGCCAGGCCAGCGCGAGCCCGCGAGCACAGCCACGAGGGGGCCAGCGAGTTGCTGG‑‑CCGGCATTGCCCGGCTGGCCa                                                                                                                                                     <  1:5186911/150‑1 (MQ=255)
   gacgTTGCCCGCACCGGTTTATCCGCAGCCCAGTGTGCCGTTGCCTCGTTCGCGGCGTGTGCTGCCGGTGGCGGCCAGGCCAGCGCGAGCCCGCGAGCACAGCCACGAGGGGGCCAGCGAGTTGCTGG‑‑CCGGCATTGCCCGGCTGGCCaa                                                                                                                                                    >  1:6177909/1‑150 (MQ=255)
   gacgTTGCCCGCACCGGTTTATCCGCAGCCCAGTGTGCCGGTGCCTCGTTCGCGGCGTGTGCTGCCGGTGGCGGCCAGGCCAGCGCGAGCCCGCGAGCACAGCCACGAGGGGGCCAGCGAGTTGCTGG‑‑CCGGCATTGCCCGGCTGGCCaa                                                                                                                                                    <  2:7813177/150‑1 (MQ=255)
    acgTTGCCCGCACCGGTTTATCCGCAGCCCAGTGTGCCGTTGCCTCGTTCGCGGCGTGTGCTGCCGGTGGCGGCCAGGCCAGCGCGAGCCCGCGAGCACAGCCACGAGGGGGCCAGCGAGTTGCTGG‑‑CCGGCATTGCCCGGCTGGCCAAt                                                                                                                                                   <  2:3061325/150‑1 (MQ=255)
      gTTGCCCGCACCGGTTTATCCGCAGCCCAGTGTGCCGTTGCCTCGTTCGCGGCGTGTGCTGCCGGTGGCGGCCAGGCCAGCGCGAGCCCGCGAGCACAGCCACGAGGGGGCCAGCGAGTTGCTGG‑‑CCGGCATTGCCCGGCTGGCCAATGc                                                                                                                                                 <  2:7681680/150‑1 (MQ=255)
      gTTGCCCGCACCGGTTTATCCGCAGCCCAGTGTGCCGTTGCCTCGTTCGCGGCGTGTGCTGCCGGTGGCGGCCAGGCCAGCGCGAGCCCGCGAGCACAGCCACGAGGGGGCCAGCGAGTTGCTGA‑‑CCGGCATTGGCCAGCCGGGCAATGc                                                                                                                                                 <  2:1583838/150‑1 (MQ=255)
       ttGCCCGCACCGGTTTATCCGCAGCCCAGTGTGCCGTTGCCTCGTTCGCGGCGTGTGCTGCCGGTGGCGGCCAGGCCAGCGCGAGCCCGCGAGCACAGCCACGAGGGGGCCAGCGAGTTGCTGG‑‑CCGGCATTGCCCGGCTGGCCAATGcc                                                                                                                                                <  1:6325642/150‑1 (MQ=255)
        tGCCCGCACCGGTTTATCCGCAGCCCAGTGTGCCGTTGCCTCGTTCGCGGCGTGTGCTGCCGGTGGCGGCCAGGCCAGCGCGAGCCCGCGAGCACAGCCACGAGGGGGCCAGCGAGTTGCTGG‑‑CCGGCATTGCCCGGCTGGCCAATGCCg                                                                                                                                               <  1:8493241/150‑1 (MQ=255)
           ccGCACCGGTTTATCCGCAGCCCAGTGTGCCGTTGCCTCGTTCGCGGCGTGTGCTGCCGGTGGCGGCCAGGCCAGCGCGAGCCCGCGAGCACAGCCACGAGGGGGCCAGCGAGTTGCTGG‑‑CCGGCATTGCCCGGCTGGCCAATGCCGGTg                                                                                                                                            >  1:1247061/1‑150 (MQ=255)
           ccGCACCGGTTTATCCGCAGCCCAGTGTGCCGTTGCCTCGTTCGCGGCGTGTGCTGCCGGTGGCGGCCAGGCCAGCGCGAGCCCGCGAGCACAGCCACGAGGGGGCCAGCGAGTTGCTGG‑‑CCGGCATTGCCCGGCTGGCCAATGCCGGTg                                                                                                                                            >  2:3173623/1‑150 (MQ=255)
           ccGCACCGGTTTATCCGCAGCCCAGTGTGCCGTTGCCTCGTTCGCGGCGTGTGCTGCCGGTGGCGGCCAGGCCAGCGCGAGCCCGCGAGCACAGCCACGAGGGGGCCAGCGAGTTGCTGG‑‑CCGGCATTGCCCGGCTGGCCAATGCCGGTg                                                                                                                                            >  2:3961060/1‑150 (MQ=255)
                  ggTTTATCCGCAGCCCAGTGTGCCGTTGCCTCGTTCGCGGCGTGTGCTGCCGGTGGCGGCCAGGCCAGCGCGAGCCCGCGAGCACAGCCACGAGGGGGCCAGCGAGTTGCTGG‑‑CCGGCATTGCCCGGCTGGCCAATGCCGGTGCCAGCGa                                                                                                                                     <  2:5937797/150‑1 (MQ=255)
                  ggTTTATCCGCAGCCCAGTGTGCCGTTGCCTCGGTCGCGGCGTGTGCTGCCGGTGGCGGCCAGGCCAGCGCGAGCCCGCGAGCACAGCCACGAGGGGGCCAGCGAGTTGCTGG‑‑CCGGCATTGCCCGGCTGGCCAATGCCGGTGCCAGCGa                                                                                                                                     <  2:8424593/150‑1 (MQ=255)
                  ggTATATCCGCAGCCCAGTGTGCCGTTGCCTCGTTCGCGGCGTGTGCTGCCGGTGGCGGCCAGGCCAGCGCGAGCCCGCGAGCACAGCCACGAGGGGGCCAGCGAGTTGCTGG‑‑CCGGCATTGCCCGGCTGGCCAATGCCGGTGCCAGCGa                                                                                                                                     <  1:1573173/150‑1 (MQ=255)
                   gTTTATCCGCAGCCCAGTGTGCCGTTGCCTCGTTCGCGGCGTGTGCTGCCGGTGGCGGCCAGGCCAGCGCGAGCCCGCGAGCACAGCCACGAGGGGGCCAGCGAGTCGCTGGCACCGGCATTGGCCAGCCGGGCAATGCCggccagcaa                                                                                                                                       >  1:2274849/1‑142 (MQ=37)
                     ttATCCGCAGCCCAGTGCGCCGTTGCCTCGTTCGCGGCGTGTGCTGCCGGTGGCGGCCAGGCCAGCGCGAGCCCGCGAGCACAGCCACGAGGGGGCCAGCGAGTTGCTGG‑‑CCGGCATTGCCCGGCTGGCCAATGCCGGTGCCAGCGCGCa                                                                                                                                  <  2:4935317/150‑1 (MQ=255)
                          cGCAGCCCAGTGTGCCGTTGCCTCGTTCGCGGCGTGTGCTGCCGGTGGCGGCCAGGCCAGCGCGAGCCCGCGAGCACAGCCACGAGGGGGCCAGCGAGTTGCTGG‑‑CCGGCATTGCCCGGCTGGCCAATGCCGGTGCCAGCGAGCAGGccc                                                                                                                             >  2:2499984/1‑150 (MQ=255)
                           gCAGCCCAGTGTGCCGTTGCCTCGTTCGCGGCGTGTGCTGCCGGTGGCGGCCAGGCCAGCGCGAGCCCGCGAGCACAGCCACGAGGGGGCCAGCGAGTTGCTGG‑‑CCGGCATTGCCCGGCTGGCCAATGCCGGTGCCAGCGAGCAGGCCCg                                                                                                                            >  2:28699/1‑150 (MQ=255)
                           gCAGCCCAGTGTGCCGTTGCCTCGTTCGCGGCGTGTGCTGCCGGTGGCGGCCAGGCCAGCGCGAGCCCGCGAGCACAGCCACGAGGGGGCCAGCGAGTTGCTGG‑‑CCGGCATTGCCCGGCTGGCCAATGCCGGTGCCAGCGAGCAGGCCCg                                                                                                                            >  1:603107/1‑150 (MQ=255)
                           gCAGCCCAGTGTGCCGTTGCCTCGTTCGCGGCGTGTGCTGCCGGTGGCGGCCAGGCCAGCGCGAGCCCGCGAGCACAGCCACGAGGGGGCCAGCGAGTTGCTGG‑‑CCGGCATTGCCCGGCTGGCCAATGCCGGTGCCAGCGAGCAGGCCCg                                                                                                                            >  2:7429774/1‑150 (MQ=255)
                               cccAGTGTGCCGTTGCCTCGTTCGCGGCGTGTGCTGCCGGTGGCGGCCAGGCCAGCGCGAGCCCGCGAGCACAGCCACGAGGGGGCCAGCGAGTTGCTGG‑‑CCGGCATTGCCCGGCTGGCCAATGCCGGTGCCAGCGAGCAGGCCCGTAGc                                                                                                                        >  2:5728674/1‑150 (MQ=255)
                               cccAGTGTGCCGTTGCCTCGTTCGCGGCGTGTGCTGCCGGTGGCGGCCAGGCCAGCGCGAGCCCGCGAGCACAGCCACGAGGGGGCCAGCGAGTTGCTGG‑‑CCGGCATTGCCCGGCTGGCCAATGCCGGTGCCAGCGAGCAGGCCCGTAGc                                                                                                                        >  1:1656724/1‑150 (MQ=255)
                               ccaaGTGTGCCGTTGCCTCGTTCGCGGCGTGTGCTGCCGGTGGCGGCCAGGCCAGCGCGAGCCCGCGAGCACAGCCACGAGGGGGCCAGCGAGTTGCTGG‑‑CCGGCATTGCCCGGCTGGCCAATGCCGGTGCCAGCGAGCAGGCCCGTAGc                                                                                                                        >  1:3640972/4‑150 (MQ=255)
                                 cAGTGTGCCGTTGCCTCGTTCGCGGCGTGTGCTGCCGGTGGCGGCCAGGCCAGCGCGAGCCCGCGAGCACAGCCACGAGGGGGCCAGCGAGTTGCTGG‑‑CCGGCATTGCCCGGCTGGCCAATGCCGGTGCCAGCGAGCAGGCCCGTAGCGa                                                                                                                      >  2:2000649/1‑150 (MQ=255)
                                  aGTGTGCCGTTGCCTCGTTCGCGGCGTGTGCTGCCGGTGGCGGCCAGGCCAGCGCGAGCCCGCGAGCACAGCCACGAGGGGGCCAGCGAGTTGCTGG‑‑CCGGCATTGCCCGGCTGGCCAATGCCGGTGCCAGCGAGCAGGCCCGTAGCGaa                                                                                                                     <  1:1844270/150‑1 (MQ=255)
                                    tgtgCCGTTGCCTCGTTCGCGGCGTGTGCTGCCGGTGGCGGCCAGGCCAGCGCGAGCCCGCGAGCACAGCCACGAGGGGGCCAGCGAGTTGCTGG‑‑CCGGCATTGCCCGGCTGGCCAATGCCGGTGCCAGCGAGCAGGCCCGTAGCGAATg                                                                                                                   >  2:419439/1‑150 (MQ=255)
                                     gtgCCGTTGCCTCGTTCGCGGCGTGTGCTGCCGGTGGCGGCCAGGCCAGCGCGAGCCCGCGAGCACAGCCACGAGGGGGCCAGCGAGTTGCTGG‑‑CCGGCATTGCCCGGCTGGCCAATGCCGGTGCCAGCGAGCAGGCCCGTAGCGAATGc                                                                                                                  <  2:4810630/150‑1 (MQ=255)
                                       gCCGTTGCCTCGTTCGCGGCGTGTGCTGCCGGTGGCGGCCAGGCCAGCGCGAGCCCGCGAGCACAGCCACGAGGGGGCCAGCGAGTTGCTGG‑‑CCGGCATTGCCCGGCTGGCCAATGCCGGTGCCAGCGAGCAGGCCCGTAGCGAATGCCa                                                                                                                <  2:6444260/150‑1 (MQ=255)
                                        ccGTTGCCTCGTTCGCGGCGTGTGCTGCCGGTGGCGGCCAGGCCAGCGCGAGCCCGCGAGCACAGCCACGAGGGGGCCAGCGAGTTGCTGG‑‑CCGGCATTGCCCGGCTGGCCAATGCCGGTGCCAGCGAGCAGGCCCGTAGCGAATGCCaa                                                                                                               <  2:6590416/150‑1 (MQ=255)
                                        ccGTTGCCTCGTTCGCGGCGTGTGCTGCCGGTGGCGGCCAGGCCAGCGCGAGCCCGCGAGCACAGCCACGAGGGGGCCAGCGAGTTGCTGG‑‑CCGGCATTGCCCGGCTGGCCAATGCCGGTGCCAGCGAGCAGGCCCGTAGCGAATGCCaa                                                                                                               <  1:8304351/150‑1 (MQ=255)
                                          gTTGCCTCGTTCGCGGCGTGTGCTGCCGGTGGCGGCCAGGCCAGCGCGAGCCCGCGAGCACAGCCACGAGGGGGCCAGCGAGTTGCTGG‑‑CCGGCATTGCCCGGCTGGCCAATGCCGGTGCCAGCGAGCAGGCCCGTAGCGAATGCCAACg                                                                                                             >  1:4469148/1‑150 (MQ=255)
                                              ccTCGTTCGCGGCGTGTGCTGCCGGTGGCGGCCAGGCCAGCGCGAGCCCGCGAGCACAGCCACGAGGGGGCCAGCGAGTTGCTGG‑‑CCGGCATTGCCCGGCTGGCCAATGCCGGTGCCAGCGAGCAGGCCCGTAGCGAATGCCAACGCTAc                                                                                                         <  2:883675/150‑1 (MQ=255)
                                               cTCGTTCGCGGCGTGTGCTGCCGGTGGCGGCCAGGCCAGCGCGAGCCCGCGAGCACAGCCACGAGGGGGCCAGCGAGTTGCTGG‑‑CCGGCATTGCCCGGCTGGCCAATGCCGGTGCCAGCGAGCAGGCCCGTAGCGAATGCCAACGCTAcc                                                                                                        >  2:1414962/1‑150 (MQ=255)
                                               cTCGTTCGCGGCGTGTGCTGCCGGTGGCGGCCAGGCCAGCGCGAGCCCGCGAGCACAGCCACGAGGGGGCCAGCGAGTTGCTGG‑‑CCGGCATTGCCCGGCTGGCCAATGCCGGTGCCAGCGAGCAGGCCCGTAGCGAATGCCAACGCTAcc                                                                                                        >  2:8026308/1‑150 (MQ=255)
                                               cTCGTTCGCGGCGTGTGCTGCCGGTGGCGGCCAGGCCAGCGCGAGCCCGCGAGCACAGCCACGAGGGGGCCAGCGAGTTGCTGG‑‑CCGGCATTGCCCGGCTGGCCAATGCCGGTGCCAGCGAGCAGGCCCGTAGCGAATGCCAACGCTAcc                                                                                                        >  1:2002212/1‑150 (MQ=255)
                                                tCGTTCGCGGCGTGTGCTGCCGGTGGCGGCCAGGCCAGCGCGAGCCCGCGAGCACAGCCACGAGGGGGCCAGCGAGTTGCTGG‑‑CCGGCATTGCCCGGCTGGCCAATGCCGGTGCCAGCGAGCAGGCCCGTAGCGAATGCCAACGCTACCt                                                                                                       >  1:2861783/1‑150 (MQ=255)
                                                tCGTTCGCGGCGTGTGCTGCCGGTGGCGGCCAGGCCAGCGCGAGCCCGCGAGCACAGCCACGAGGGGGCCAGCGAGTTGCTGG‑‑CCGGCATTGCCCGGCTGGCCAATGCCGGTGCCAGCGAGCAGGCCCGTAGCGAATGCCAACGCTACCt                                                                                                       >  1:868339/1‑150 (MQ=255)
                                                tCGTTCGCGGCGTGTGCTGCCGGTGGCGGCCAGGCCAGCGCGAGCCCGCGAGCACAGCCACGAGGGGGCCAGCGAGTTGCTGG‑‑CCGGCATTGCCCGGCTGGCCAATGCCGGTGCCAGCGAGCAGGCCCGTAGCGAATGCCAACGCTACCt                                                                                                       >  2:5637315/1‑150 (MQ=255)
                                                tCGTTCGCGGCGTGTGCTGCCGGTGGCGGCCAGGCCAGCGCGAGCCCGCGAGCACAGCCACGAGGGGGCCAGCGAGTTGCTGG‑‑CCGGCATTGCCCGGCTGGCCAATGCCGGTGCCAGCGAGCAGGCCCGTAGCGAATGCCAACGCTACCt                                                                                                       >  2:412697/1‑150 (MQ=255)
                                                tCGTTCGCGGCGTGTGCTGCCGGTGGCGGCCAGGCCAGCGCGAGCCCGCGAGCACAGCCACGAGGGGGCCAGCGAGTTGCTGG‑‑CCGGCATTGCCCGGCTGGCCAATGCCGGTGCCAGCGAGCAGGCCCGTAGCGAATGCCAACGCTACCt                                                                                                       <  1:3147714/150‑1 (MQ=255)
                                                tCGTTCGCGGCGTGTGCTGCCGGTGGCGGCCAGGCCAGCGCGAGCCCGCGAGCACAGCCACGAGGGGGCCAGCGAGTTGCTGG‑‑CCGGCATTGCCCGGCTGGCCAATGCCGGTGCCAGCGAGCAGGCCCGTAGCGAATGCCAACGCTACCt                                                                                                       >  1:8383972/1‑150 (MQ=255)
                                                    tCGCGGCGTGTGCTGCCGGTGGCGGCCAGGCCAGCGCGAGCCCGCGAGCACAGCCACGAGGGGGCCAGCGAGTTGCTGG‑‑CCGGCATTGCCCGGCTGGCCAATGCCGGTGCCAGCGAGCAGGCCCGTAGCGAATGCCAACGCTACCTGAGc                                                                                                   <  2:8679589/150‑1 (MQ=255)
                                                    tCGCGGCGTGTGCTGCCGGTGGCGGCCAGGCCAGCGCGAGCCCGCGAGCACAGCCACGAGGGGGCCAGCGAGTTGCTGG‑‑CCGGCATTGCCCGGCTGGCCAATGCCGGTGCCAGCGAGCAGGCCCGTAGCGAATGCCAACGCTACCTGAGc                                                                                                   <  1:7876078/150‑1 (MQ=255)
                                                      gcggcgTGTGCTGCCGGTGGCGGCCAGGCCAGCGCGAGCCCGCGAGCACAGCCACGAGGGGGCCAGCGAGTTGCTGG‑‑CCGGCATTTCCCGGCTGGCCAATGCCGGTGCCAGCGAGCAGGCCCGTAGCGAATGCCAACGCTACCTGAGCCa                                                                                                 <  2:126628/150‑1 (MQ=255)
                                                      gcggcgTGTGCTGCCGGTGGCGGCCAGGCCAGCGCGAGCCCGCGAGCACAGCCACGAGGGGGCCAGCGAGTTGCTGG‑‑CCGGCATTGCCCGGCTGGCCAATGCCGGTGCCAGCGAGCAGGCCCGTAGCGAATGCCAACGCTACCTGAGCCa                                                                                                 <  1:5672383/150‑1 (MQ=255)
                                                       cggcgTGTGCTGCCGGTGGCGGCCAGGCCAGCGCGAGCCCGCGAGCACAGCCACGAGGGGGCCAGCGAGTTGCTGG‑‑CCGGCATTGCCCGGCTGGCCAATGCCGGTGCCAGCGAGCAGGCCCGTAGCGAATGCCAACGCTACCTGAGCCAg                                                                                                >  2:5515360/1‑150 (MQ=255)
                                                          cgTGTGCTGCCGGTGGCGGCCAGGCCAGCGCGAGCCCGCGAGCACAGCCACGAGGGGGCCAGCGAGTTGCTGG‑‑CCGGCATTGCCCGGCTGGCCAATGCCGGTGCCAGCGAGCAGGCCCGTAGCGAATGCCAACGCTACCTGAGCCAGTAc                                                                                             <  1:504661/150‑1 (MQ=255)
                                                          cgTGTGCTGCCGGTGGCGGCCAGGCCAGCGCGAGCCCGCGAGCACAGCCACGAGGGGGCCAGCGAGTTGCTGG‑‑CCGGCATTGCCCGGCTGGCCAATCCCGGTGCCAGCGAGCAGGCCCGTAGCGACAGCCAACGCTACCTGAGCCAGTAc                                                                                             >  1:4943237/1‑150 (MQ=255)
                                                              tgctgcCGGTGGCGGCCAGGCCAGCGCGAGCCCGCGAGCACAGCCACGAGGGGGCCAGCGAGTTGCTGG‑‑CCGGCATTGCCCGGCTGGCCAATGCCGGTGCCAGCGAGCAGGCCCGTAGCGAATGCCAACGCTACCTGAGCCAGTACccgc                                                                                         >  1:5340769/1‑150 (MQ=255)
                                                                     ggTGGCGGCCAGGCCAGCGCGAGCCCGCGAGCACAGCCCCGAGGGGGCCAGCGAGTTGCTGG‑‑CCGGCATTGCCCGGCTGGCCAATGCCGGTGCCAGCGAGCAGGCCCGTAGCGAATGCCAACGCTACCTGAGCCAGTACCCGCCGTCGGc                                                                                  <  2:7091244/150‑1 (MQ=255)
                                                                       tGGCGGCCAGGCCAGCGCGAGCCCGCGAGCACAGCCACGAGGGGGCCAGCGAGTTGCTGG‑‑CCGGCATTGCCCGGCTGGCCAATGCCGGTGCCAGCGAGCAGGCCCGTAGCGAATGCCAACGCTACCTGAGCCAGTACCCGCCGTCGGcac                                                                                >  1:2211042/1‑150 (MQ=255)
                                                                         gcggcCAGGCCAGCGCGAGCCCGCGAGCACAGCCACGAGGGGGCCAGCGAGTTGCTGG‑‑CCGGCATTGCCCGGCTGGCCAATGCCGGTGCCAGCGAGCAGGCCCGTAGCGAATGCCAACGCTCCCTGAGCCAGTACCCGCCGTCGGCACAg                                                                              >  1:2504580/1‑150 (MQ=255)
                                                                         gcggcCAGGCCAGCGCGAGCCCGCGAGCACAGCCACGAGGGGGCCAGCGAGTTGCTGG‑‑CCGGCATTGCCCGGCTGGCCAATGCCGGTGCCAGCGAGCAGGCCCGTAGCGAATGCCAACGCTACCTGAGCCAGTACCCGCCGTCGGCACAg                                                                              >  2:5126921/1‑150 (MQ=255)
                                                                           ggccaggccaGCGCGAGCCCGCGAGCACAGCCACGAGGGGGCCAGCGAGTTGCTGG‑‑CCGGCATTGCCCGGCTGGCCAATGCCGGTGCCAGCGAGCAGGCCCGTAGCGAATGCCAACGCTACCTGAGCCAGTACCCGCCGTCGGCACAGgt                                                                            >  1:5495806/1‑150 (MQ=255)
                                                                           ggccaggccaGCGCGAGCCCGCGAGCACAGCCACGAGGGGGCCAGCGAGTTGCTGG‑‑CCGGCATTGCCCGGCTGGCCAATGCCGGTGCCAGCGAGCAGGCCCGTAGCGAATGCCAACGCTACCTGAGCCAGTACCCGCCGTCGGCACAGgt                                                                            >  1:6779910/1‑150 (MQ=255)
                                                                            gccaggccagCGCGAGCCCGCGAGCACAGCCACGAGGGGGCCAGCGAGTTGCTGG‑‑CCGGCATTGCCCGGCTGGCCAATGCCGGTGCCAGCGAGCAGGCCCGTAGCGAATGCCAACGCTACCTGAGCCAGTACCCGCCGTCGGCACAGgtg                                                                           >  1:4594294/1‑150 (MQ=255)
                                                                            gccaggccagCGCGAGCCCGCGAGCACAGCCACGAGGGGGCCAGCGAGTTGCTGG‑‑CCGGCATTGCCCGGCTGGCCAATGCCGGTGCCAGCGAGCAGGCCCGTAGCGAATGCCAACGCTACCTGAGCCAGTACCCGCCGTCGGCACAGgtg                                                                           >  1:7299426/1‑150 (MQ=255)
                                                                            gccaggccagCGCGAGCCCGCGAGCACAGCCACGAGGGGGCCAGCGAGTTGCTGG‑‑CCGGCATTGCCCGGCTGGCCAATGCCGGTGCCAGCGAGCAGGCCCGTAGCGAATGCCAACGCTACCTGAGCCAGTACCCGCCGTCGGCACAGgtg                                                                           >  2:8221028/1‑150 (MQ=255)
                                                                             ccaggccagCGCGAGCCCGCGAGCACAGCCACGAGGGGGCCAGCGAGTTGCTGG‑‑CCGGCATTGCCCGGCTGGCCAATGCCGGTGCCAGCGAGCAGGCCCGTAGCGAATGCCAACGCTACCTGAGCCAGTACCCGCCGTCGGCACAGgtgt                                                                          <  2:1802643/150‑1 (MQ=255)
                                                                              caggccagCGCGAGCCCGCGAGCACAGCCACGAGGGGGCCAGCGAGTTGCTGG‑‑CCGGCATTGCCCGGCTGGCCAATGCCGGTGCCAGCGAGCAGGCCCGTAGCGAATGCCAACGCTACCTGAGCCAGTACCCGCCGTCGGCACAGGTGta                                                                         <  2:6789920/150‑1 (MQ=255)
                                                                              cagcccagCGCCAGCCCGCCAGCACAGCCCCGAGGGGGCCAGCGAGTTGCCGG‑‑CCGGCATTGCCCGGCTCGCCTATGCCGGTGCCAGCGAGCAGGCCCGTAGCGAATGCCAACGCAACCTGAGCCAGTACCCGCCGTCGGCACAGGTGta                                                                         <  2:8520284/150‑1 (MQ=255)
                                                                               aggccagCGCGAGCCCGCGAGCACAGCCACGAGGGGGCCAGCGAGTTGCTGG‑‑CCGGCATTGCCCGGCTGGCCAATGCCGGTGCCAGCGAGCAGGCCCGTAGCGAATGCCAACGCTACCTGAGCCAGTACCCGCCGTCGGCACAGGTGtac                                                                        <  2:1037450/150‑1 (MQ=255)
                                                                               aggccagCGCGAGCCCGCGAGCACAGCCACGAGGGGGCCAGCGAGTTGCTGG‑‑CCGGCATTGCCCGGCTGGCCAATGCCGGTGCCAGCGAGCAGGCCCGTAGCGAATGCCAACGCTACCTGAGCCAGTACCCGCCGTCGGCACAGGTGtac                                                                        >  2:6105209/1‑150 (MQ=255)
                                                                                     gcgcGAGCCCGCGAGCACAGCCACGAGGGGGCCAGCGAGTTGCTGG‑‑CCGGCATTGCCCGGCTGGCCAATGCCGGTGCCAGCGAGCAGGCCCGTAGCGAATGCCAACGCTACCTGAGCCAGTACCCGCCGTCGGCACAGGTGTACTActgg                                                                  >  2:3584002/1‑150 (MQ=255)
                                                                                       gcgAGCCCGCGAGCACAGCCACGAGGGGGCCAGCGAGTTGCTGG‑‑CCGGCATTGCCCGGCTGGCCAATGCCGGTGCCAGCGAGCAGGCCCGTAGCGAATGCCAACGCTACCTGAGCCAGTACCCGCCGTCGGCACAGGTGTACTActggct                                                                <  1:1719620/150‑1 (MQ=255)
                                                                                             ccGCGAGCACAGCCACGAGGGGGCCAGCGAGTTGCTGG‑‑CCGGCATTGCCCGGCTGGCCAATGCCGGTGCCAGCGAGCAGGCCCGTAGCGAATGCCAACGCTACCTGAGCCAGTACCCGCCGTCGGCACAGGTGTACTACTGGCTGGGttt                                                          <  1:7514435/150‑1 (MQ=255)
                                                                                             ccGCGAGCACAGCCACGAGGGGGCCAGCGAGTTGCTGG‑‑CCGGCATTGCCCGGCTGGCCAATGCCGGTGCCAGCGAGCAGGCCCGTAGCGAATGCCAACGCTACCTGAGCCAGTACCCGCCGTCGGCACAGGTGTACTACTGGCTGGGttt                                                          >  2:6281991/1‑150 (MQ=255)
                                                                                               gcgAGCACAGCCACGAGGGGGCCAGCGAGTTGCTGG‑‑CCGGCATTGCCCGGCTGGCCAATGCCGGTGCCAGCGAGCAGGCCCGTAGCGAATGCCAACGCTACCTGAGCCAGTACCCGCCGTCGGCACAGGTGTACTACTGGCTGGGTTTGc                                                        <  2:5426508/150‑1 (MQ=255)
                                                                                                  aGCACAGCCACGAGGGGGCCAGCGAGTTGCTGG‑‑CCGGCATTGCCCGGCTGGCCAATGCCGGTGCCAGCGAGCAGGCCCGTAGCGAATGCCAACGCTACCTGAGCCAGTACCCGCCGTCGGCACAGGTGTACTACTGGCTGGGTTTGCTCa                                                     <  1:2740176/150‑1 (MQ=255)
                                                                                                   gCACAGCCACGAGGGGGCCAGCGAGTTGCTGG‑‑CCGGCATTGCCCGGCTGGCCAATGCCGGTGCCAGCGAGCAGGCCCGTAGCGAATGCCAACGCTACCTGAGCCAGTACCCGCCGTCGGCACAGGTGTACTACTGGCTGGGTTTGCTCAg                                                    <  1:965294/150‑1 (MQ=255)
                                                                                                   gCACAGCCACGAGGGGGCCAGCGAGTTGCTGG‑‑CCGGCATTGCCCGGCTGGCCAATGCCGGTGCCAGCGAGCAGGCCCGTAGCGAATGCCAACGCTACCTGAGCCAGTACCCGCCGTCGGCACAGGTGTACTACTGGCTGGGTTTGCTCAg                                                    <  1:4369346/150‑1 (MQ=255)
                                                                                                    cacaGCCACGAGGGGGCCAGCGAGTTGCTGG‑‑CCGGCATTGCCCGGCTGGCCAATGCCGGTGCCAGCGAGCAGGCCCGTAGCGAATGCCAACGCTACCTGAGCCAGTACCCGCCGTCGGCACAGGTGTACTACTGGCTGGGTTTGCTCAGt                                                   <  1:4086575/150‑1 (MQ=255)
                                                                                                     acaGCCACGAGGGGGCCAGCGAGTTGCTGG‑‑CCGGCATTGCCCGGCTGGCCAATGCCGGTGCCAGCGAGCAGGCCCGTAGCGAATGCCAACGCTACCTGAGCCAGTACCCGCCGTCGGCACAGGTGTACTACTGGCTGGGTTTGCTCAGTg                                                  >  1:1079100/1‑150 (MQ=255)
                                                                                                       aGCCACGAGGGGGCCAGCGAGTTGCTGG‑‑CCGGCATTGCCCGGCTGGCCAATGCCGGTGCCAGCGAGCAGGCCCGTAGCGAATGCCAACGCTACCTGAGCCAGTACCCGCCGTCGGCACAGGTGTACTACTGGCTGGGTTTGCTCAGTGAt                                                >  2:4651524/1‑150 (MQ=255)
                                                                                                       aGCCACGAGGGGGCCAGCGAGTTGCTGG‑‑CCGGCATTGCCCGGCTGGCCAATGCCGGTGCCAGCGAGCAGGCCCGTAGCGAATGCCAACGCTACCTGAGCCAGTACCCGCCGTCGGCACAGGTGTACTACTGGCTGGGTTTGCTCAGTGAt                                                >  2:2137670/1‑150 (MQ=255)
                                                                                                        gCCACGAGGGGGCCAGCGAGTTGCTGG‑‑CCGGCATTGCCCGGCTGGCCAATGCCGGTGCCAGCGAGCAGGCCCGTAGCGAATGCCAACGCTACCTGAGCCAGTACCCGCCGTCGGCACAGGTGTACTACTGGCTGGGTTTGCTCAGTGATa                                               <  2:2393425/150‑1 (MQ=255)
                                                                                                           aCGAGGGGGCCAGCGAGTTGCTGG‑‑CCGGCATTGCCCGGCTGGCCAATGCCGGTGCCAGCGAGCAGGCCCGTAGCGAATGCCAACGCTACCTGAGCCAGTACCCGCCGTCGGCACAGGTGTACTACTGGCTGGGTTTGCTCAGTGATACCg                                            >  1:5891982/1‑150 (MQ=255)
                                                                                                             gAGGGGGCCAGCGAGTTGCTGG‑‑CCGGCATTGCCCGGCTGGCCAATGCCGGTGCCAGCGAGCAGGCCCGTAGCGAATGCCAACGCTACCTGAGCCAGTACCCGCCGTCGGCACAGGTGTACTACTGGCTGGGTTTGCTCAGTGATACCGaa                                          <  1:7686427/150‑1 (MQ=255)
                                                                                                             gAGGGGGCCAGCGAGTTGCTGG‑‑CCGGCATTGCCCGGCTGGCCAATGCCGGTGCCAGCGAGCAGGCCCGTAGCGAATGCCAACGCTACCTGAGCCAGTACCCGCCGTCGGCACAGGTGTACTACTGGCTGGGTTTGCTCAGTGATACCGaa                                          <  1:4804415/150‑1 (MQ=255)
                                                                                                             gAGGGGGCCAGCGAGTTGCTGG‑‑CCGGCATTGCCCGGCTGGCCAATGCCGGTGCCAGCGAGCAGGCCCGTAGCGAATGCCAACGCTACCTGAGCCAGTACCCGCCGTCGGCACAGGTGTACTACTGGCTGGGTTTGCTCAGTGATACCGaa                                          <  1:4008007/150‑1 (MQ=255)
                                                                                                             gAGGGGGCCAGCGAGTTGCTGG‑‑CCGGCATTGCCCGGCTGGCCAATGCCGGTGCCAGCGAGCAGGCCCGTAGCGAATGCCAACGCTACCTGAGCCAGTACCCGCCGTCGGCACAGGTGTACTACTGGCTGGGTTTGCTCAGTGATACCGaa                                          >  2:1197805/1‑150 (MQ=255)
                                                                                                               gggggCCAGCGAGTTGCTGG‑‑CCGGCATTGCCCGGCTGGCCAATGCCGGTGCCAGCGAGCAGGCCCGTAGCGAATGCCAACGCTACCTGAGCCAGTACCCGCCGTCGGCACAGGTGTACTACTGGCTGGGTTTGCTCAGTGATACCGAAgg                                        <  2:2976812/150‑1 (MQ=255)
                                                                                                               gggggCCAGCGAGTTGCTGG‑‑CCGGCATTGCCCGGCTGGCCAATGCCGGTGCCAGCGAGCAGGCCCGTAGCGAATGCCAACGCTACCCGAGCCAGTACCCGCCGTCGGCACAGGTGTACTACTGGCTGGGTTTGCTCAGTGATACCGAAgg                                        <  1:4164253/150‑1 (MQ=255)
                                                                                                                 gggCCAGCGAGTTGCTGG‑‑CCGGCATTGCCCGGCTGGCCAATGCCGGTGCCAGCGAGCAGGCCCGTAGCGAATGCCAACGCTACCTGAGCCAGTACCCGCCGTCGGCACAGGTGTACTACTGGCTGGGTTTGCTCAGTGATACCGAAGGCg                                      <  2:4373449/150‑1 (MQ=255)
                                                                                                                  ggCCAGCGAGTTGCTGG‑‑CCGGCATTGCCCGGCTGGCCAATGCCGGTGCCAGCGAGCAGGCCCGTAGCGAATGCCAACGCTACCTGAGCCAGTACCCGCCGTCGGCACAGGTGTACTACTGGCTGGGTTTGCTCAGTGATACCGAAGGCGa                                     <  2:5232404/150‑1 (MQ=255)
                                                                                                                     cAGCGAGTTGCTGG‑‑CCGGCATTGCCCGGCTGGCCAATGCCGGTGCCAGCGAGCAGGCCCGTAGCGAATGCCAACGCTACCTGAGCCAGTACCCGCCGTCGGCACAGGTGTACTACTGGCTGGGTTTGCTCAGTGATACCGAAGGCGATgc                                  >  1:3325002/1‑150 (MQ=255)
                                                                                                                     cAGCGAGTTGCTGG‑‑CCGGCATTGCCCGGCTGGCCAATGCCGGTGCCAGCGAGCAGGCCCGTAGCGAATGCCAACGCTACCTGAGCCAGTACCCGCCGTCGGCACAGGTGTACTACTGGCTGGGTTTGCTCAGTGATACCGAAGGCGATgc                                  >  1:8058440/1‑150 (MQ=255)
                                                                                                                     cAGCGAGTTGCTGG‑‑CCGGCATTGCCCGGCTGGCCAATGCCGGTGCCAGCGAGCAGGCCCGTAGCGAATGCCAACGCTACCTGAGCCAGTACCCGCCGTCGGCACAGGTGTACTACTGGCTGGGTTTGCTCAGTGATACCGAAGGCGATgc                                  >  1:3012305/1‑150 (MQ=255)
                                                                                                                          aGTTGCTGG‑‑CCGGCATTGCCCGGCTGGCCAATGCCGGTGCCAGCGAGCAGGCCCGTAGCGAATGCCAACGCTACCTGAGCCAGTACCCGCCGTCGGCACAGGTGTACTACTGGCTGGGTTTGCTCAGTGATACCGAAGGCGATGCgcagc                             <  1:3736652/150‑1 (MQ=255)
                                                                                                                              gCTGG‑‑CCGGCATTGCCCGGCTGGCCAATGCCGGTGCCAGCGAGCAGGCCCGTAGCGAATGCCAACGCTACCTGAGCCAGTACCCGCCGTCGGCACAGGTGTACTACTGGCTGGGTTTGCTCAGTGATACCGAAGGCGATGCGCAGCAGgc                         >  2:6214122/1‑150 (MQ=255)
                                                                                                                              gCTGG‑‑CCGGCATTGCCCGGCTGGCCAATGCCGGTGCCAGCGAGCAGGCCCGTAGCGAATGCCAACGCTACCTGAGCCAGTACCCGCCGTCGGCACAGGTGTACTACTGGCTGGGTGTGCTCAGTGATACCGAAGGCGATGCGCAGCAGgc                         >  1:3866998/1‑150 (MQ=255)
                                                                                                                                tggcaccGGCATTGGCCAGCCGGGCAATGCCG‑‑GCCAGCGAGCAGGCCCGTAGCGAATGCCAACGCTACCTGAGCCAGTACCCGCCGTCGGCACAGGTGTACTACTGGCTGGGTTTGCTCAGTGATACCGAAGGCGATGCGCAGCAGgcgc                       >  1:7089161/6‑150 (MQ=38)
                                                                                                                                tggcaccGGCATTGGCCAGCCGGGCAATGCCG‑‑GCCAGCAAGCAGGCCCGTAGCGAATGCCAACGCTACCTGAGCCAGTACCCGCCGTCGGCACAGGTGTACTACTGGCTGGGTTTGCTCAGTGATACCGAAGGCGATGCGCAGCAGgcgc                       >  1:1300548/6‑150 (MQ=38)
                                                                                                                                tGG‑‑CCGGCATTGCCCGGCTGGCCAATGCCGGTGCCAGCGAGCAGGCCCGTAGCGAATGCCAACGCTACCTGAGCCAGTACCCGCCGTCGGCACAGGTGTACTACTGGCTGGGTTTGCTCAGTGATACCGAAGGCGATGCGCAGCAGgcgc                       >  2:2655173/1‑150 (MQ=255)
                                                                                                                                  g‑‑CCGGCATTGCCCGGCTGGCCAATGCCGGTGCCAGCGAGCAGGCCCGTAGCGAATGGCAACGCTACCTGAGCCAGTACCCGCCGTCGGCACAGGTGTACTACTGGCTGGGTTTGCTCAGTGATACCGAAGGCGATGCGCAGCAGGCGCTg                     >  1:2624416/1‑150 (MQ=255)
                                                                                                                                  g‑‑CCGGCATTGCCCGGCTGGCCAATGCCGGTGCCAGCGAGCAGGCCCGTAGCGAATGCCAACGCTACCTGAGCCAGTACCCGCCGTCGGCACAGGTGTACTACTGGCTGGGTTTGCTCAGTGATACCGAAGGCGATGCGCAGCAGGCGCTg                     >  2:7136569/1‑150 (MQ=255)
                                                                                                                                     ccGGCATTGCCCGGCTGGCCAATGCCGGTGCCAGCGAGCAGGCCCGTAGCGAATGCCAACGCTACCTGAGCCAGTACCCGCCGTCGGCACAGGTGTACTACTGGCTGGGTTTGCTCAGTGATACCGAAGGCGATGCGCAGCAGGCGCTGa                    <  1:6603293/150‑1 (MQ=255)
                                                                                                                                              cccGGCTGGCCAATGCCGGTGCCAGCGAGCAGGCCCGTAGCGAATGCCAACGCTACCTGAGCCAGTACCCGCCGTCGGCACAGGTGTACTACTGGCTGGGTTTGCTCAGTGATACCGAAGGCGATGCGCAGCAGGCGCTGAGCCACTATc           >  2:7346562/1‑150 (MQ=255)
                                                                                                                                               ccGGCTGGCCAATGCCGGTGCCAGCGAGCAGGCCCGTAGCGAATGCCAACGCTACCTGAGCCAGTACCCGCCGTCGGCACAGGTGTACTACTGGCTGGGTTTGCTCAGTGATACCGAAGGCGATGCGCAGCAGGCGCTGAGCCACTATCg          >  2:5372816/1‑150 (MQ=255)
                                                                                                                                                      gCCAATGCCGGTGCCAGCGAGCAGGCCCGTAGCGAATGCCAACGCTACCTGAGCCAGTACCCGCCGTCGGCACAGGTGTACTACTGGCTGGGTTTGCTCAGTGATACCGAAGGCGATGCGCAGCAGGCGCTGAGCCACTATCGCAAGGcc   <  1:3411313/150‑1 (MQ=255)
                                                                                                                                                       |                                                                                                                                                     
CACGACGTTGCCCGCACCGGTTTATCCGCAGCCCAGTGTGCCGTTGCCTCGTTCGCGGCGTGTGCTGCCGGTGGCGGCCAGGCCAGCGCGAGCCCGCGAGCACAGCCACGAGGGGGCCAGCGAGTTGCTGG‑‑CCGGCATTGCCCGGCTGGCCAATGCCGGTGCCAGCGAGCAGGCCCGTAGCGAATGCCAACGCTACCTGAGCCAGTACCCGCCGTCGGCACAGGTGTACTACTGGCTGGGTTTGCTCAGTGATACCGAAGGCGATGCGCAGCAGGCGCTGAGCCACTATCGCAAGGCCC  >  NC_002947/1699871‑1700169

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: