Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 1,695,068 T→C 5.2% L331P (CTG→CCG)  wspC → biofilm formation methyltransferase WspC

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0029471,695,0680TC 5.2% 460.1 / 8.8 155L331P (CTG→CCG) wspCbiofilm formation methyltransferase WspC
Reads supporting (aligned to +/- strand):  ref base T (76/70);  new base C (3/5);  total (79/76)
Fisher's exact test for biased strand distribution p-value = 4.87e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 8.33e-01

Only 100 of 159 total aligned reads displayed.

GTTCACGACGTTGCCCGCACCGGTTTATCCGCAGCCCAGTGTGCCGTTGCCTCGTTCGCGGCGTGTGCTGCCGGTGGCGGCCAGGCCAGCGCGAGCCCGCGAGCACAGCCACGAGGGGGCCAGCGAGTTGCTGGCCGGCATTGCCCGGCTGGCCAATGCCGGTGCCAGCGAGCAGGCCCGTAGCGAATGCCAACGCTACCTGAGCCAGTACCCGCCGTCGGCACAGGTGTACTACTGGCTGGGTTTGCTCAGTGATACCGAAGGCGATGCGCAGCAGGCGCTGAGCCACTATCGCAAGG  >  NC_002947/1694919‑1695217
                                                                                                                                                     |                                                                                                                                                     
gtTCACGACGTTGCCCGCACCGGTTTATCCGCAGCCCAGTGTGCCGTTGCCTCGTTCGCGGCGTGTGCTGCCGGTGGCGGCCAGGCCAGCGCGAGCCCGCGAGCACAGCCACGAGGGGGCCAGCGAGTTGCTGGCCGGCATTGCCCGGCt                                                                                                                                                       <  1:1306469/150‑1 (MQ=255)
    acgacgTTGCCCGCACCGGTTTATCCGCAGCCCAGTGTGCCGTTGCCTCGTTCGCGGCGTGTGCTGCCGGTGGCGGCCAGGCCAGCGCGAGCCCGCGAGCACAGCCACGAGGGGGCCAGCGAGTTGCTGGCCGGCATTGCCCGGCTGGcc                                                                                                                                                   >  2:1249077/1‑150 (MQ=255)
         gTTGCCCGCACCGGTTTATCCGCAGCCCAGTGTGCCGTTGCCTCGTTCGCGGCGTGTGCTGCCGGTGGCGGCCAGGCCAGCGCGAGCCCGCGAGCACAGCCACGAGGGGGCCAGCGAGTTGCTGGCCGGCCTTGCCCGGCTGGCCAATGc                                                                                                                                              <  2:821279/150‑1 (MQ=255)
         gTTGCCCGCACCGGTTTATCCGCAGCCCAGTGTGCCGTTGCCTCGTTCGCGGCGTGTGCTGCCGGTGGCGGCCAGGCCAGCGCGAGCCCGCGAGCACAGCCACGAGGGGGCCAGCGAGTTGCTGGCCGGCATTGCCCGGCTGGCCAATGc                                                                                                                                              >  2:685289/1‑150 (MQ=255)
          ttGCCCGCACCGGTTTATCCGCAGCCCAGTGTGCCGTTGCCTCGTTCGCGGCGTGTGCTGCCGGTGGCGGCCAGGCCAGCGCGAGCCCGCGAGCACAGCCACGAGGGGGCCAGCGAGTTGCTGGCCGGCATTGCCCGGCTGGCCAATGcc                                                                                                                                             <  1:2903256/150‑1 (MQ=255)
          ttGCCCGCACCGGTTTATCCGCAGCCCAGTGTGCCGTTGCCTCGTTCGCGGCGTGTGCTGCCGGTGGCGGCCAGGCCAGCGCGAGCCCGCGAGCACAGCCACGAGGGGGCCAGCGAGTTGCTGGCCGGCATTGCCCGGCTGGCCAATGcc                                                                                                                                             >  2:2168087/1‑150 (MQ=255)
           tGCCCGCACCGGTTTATCCGCAGCCCAGTGTGCCGTTGCCTCGTTCGCGGCGTGTGCTGCCGGTGGCGGCCAGGCCAGCCCGAGCCCGCGAGCACAGCCACGAGGGGGCCAGCGAGTTGCTGGCCGGCATTGCCCGGCTGGCCAATGCCg                                                                                                                                            <  2:2539028/150‑1 (MQ=255)
             cccGCACCGGTTTATCCGCAGCCCAGTGTGCCGTTGCCTCGTTCGCGGCGTGTGCTGCCGGTGGCGGCCAGGCCAGCGCGAGCCCGCGAGCACAGCCACGAGGGGGCCAGCGAGTTGCTGACCGGCATTGGCCAGCCGGGCAATGCCggc                                                                                                                                          <  1:2991172/150‑2 (MQ=255)
              ccGCACCGGTTTATCCGCATCCCAGTGTGCCGTTGCCTCGTTCGCGGCGTGTGCTGCCGGTGGCGGCCAGGCCAGCGCGAGCCCGCGAGCACAGCCACGAGGGGGCCAGCGAGTTGCTGGCCGGCATTGCCCGGCTGGCCAATGCCGGTg                                                                                                                                         >  1:2257493/1‑150 (MQ=255)
               cGCACCGGTTTATCCGCAGCCCAGCGTGCCGTTGCCTCGTTCGCGGCGTGTGCTGCCGGTGGCGGCCAGGCCAGCGCGAGCCCGCGCGCACAGCCACGAGGGGGCCAGCGAGTTGCTGGCCGGCATTGCCCGGCTGGCCAATGCCGGTGc                                                                                                                                        >  2:1681369/1‑150 (MQ=255)
                gCACCGGTTTATCCGCAGCCCAGTGTGCCGTTGCCTCGTTCGCGGCGTGTGCTGCCGGTGGCGGCCAGGCCAGCGCGAGCCCGCGAGCACAGCCACGAGGGGGCCAGCGAGTTGCTGGCCGGCATTGCCCGGCTGGCCAATGCCGGTGcc                                                                                                                                       <  2:2862012/150‑1 (MQ=255)
                           tCCGCAGCCCAGTGTGCCGTTGCCTCGTTCGCGGCGTGTGCTGCCGGTGGCGGCCAGGCCAGCGCGAGCCCGCGAGCACAGCCACGAGGGGGCCAGCGAGTTGCTGGCCGGCATTGCCCGGCTGGCCAATGCCGGTGCCAGCGAGCAGGc                                                                                                                            >  1:260452/1‑150 (MQ=255)
                             cGCAGCCCAGTGTGCCGTTGCCTCGTTCGCGGCGTGTGCTGCCGGTGGCGGCCAGGCCAGCGCGAGCCCGCGAGCACAGCCACGAGGGGGCCAGCGAGTTGGTGGCCGGCATTGCCCGGCTGGCCAATGCCGGTGCCAGCGAGCAGGccc                                                                                                                          >  1:2968349/1‑150 (MQ=255)
                                aGCCCAGTGTGCCGTTGCCTCGTTCGCGGCGTGTGCTGCCGGTGGCGGCCAGGCCAGCGCGAGCCCGCGAGCACAGCCACGAGGGGGCCAGCGAGTTGCTGGCCGGCATTGCCCGGCTGGCCAATGCCGGTGCCAGCGAGCAGGCCCGTa                                                                                                                       <  1:1444932/150‑1 (MQ=255)
                                  cccAGTGTGCCGTTGCCTCGTTCGCGGCGTGTGCTGCCGGTGGCGGCCAGGCCAGCGCGAGCCCGCGAGCACAGCCCCGAGGGGGCCAGCGAGTTGCTGGCCGGCATTGCCCGGCTGGCCAATGCCGGTGCCAGCGAGCAGGGCCGTAgg                                                                                                                     >  1:1413415/1‑149 (MQ=255)
                                  cccAGTGTGCCGTTGCCTCGTTCGCGGCGTGTGCTGCCGGTGGCGGCCAGGCCAGCGCGAGCCCGCGAGCACAGCCACGAGGGGGCCAGCGAGTTGCTGGCCGGCATTGCCCGGCTGGCCAATGCCGGTGCCAGCGAGCAGGCCCGTAGc                                                                                                                     >  2:1285710/1‑150 (MQ=255)
                                  cccAGTGTGCCGTTGCCTCGTTCGCGGCGTGTGCTGCCGGTGGCGGCCAGGCCAGCGCGAGCCCGCGAGCACAGCCACGAGGGGGCCAGCGAGTTGCTGGCCGGCATTGCCCGGCTGGCCAATGCCGGTGCCAGCGAGCAGGCCCGTAGc                                                                                                                     >  2:2407529/1‑150 (MQ=255)
                                  cccAGTGTGCCGTTGCCTCGTTCGCGGCGTGTGCTGCCGGTGGCGGCCAGGCCAGCGCGAGCCCGCGAGCACAGCCACGAGGGGGCCAGCGAGTTGCTGGCCGGCATTGCCCGGCTGGCCAATGCCGGTGCCAGCGAGCAGGCCCGTAGc                                                                                                                     >  1:322960/1‑150 (MQ=255)
                                    cAGTGTGTCGTTGCCTCGTTCGCGGCGTGTGCTGCCGGTGGCGGCCAGGCCAGCGCGAGCCCGCGAGCACAGCCACGAGGGGGCCAGCGAGTTGCTGGCCGGCATTGCCCGGCTGGCCAATGCCGGTGCCAGCGAGCAGGCCCGTAGCGa                                                                                                                   >  2:1219301/1‑150 (MQ=255)
                                     aGTGTGCCGTTGCCTCGTTCGCGGCGTGTGCTGCCGGTGGCGGCCAGGCCAGCGCGAGCCCGCGAGCACAGCCACGAGGGGGCCAGCGAGTTGCTGGCCGGCATTGCCCGGCTGGCCAATGCCGGTGCCAGCGAGCAGGCCCGTAGCGaa                                                                                                                  <  2:1484670/150‑1 (MQ=255)
                                      tatcactgagcaaacccagccagtagtacacctgtgccgacggcgggtactggctcaggtagcgttggcattcgctacGGGCCAGCGAGTTGCTGGCCGGCATTGCCCGGCTGGCCAATGCCGGTGCCAGCGAGCAGGCCCGTAGCGAAt                                                                                                                 >  1:1503501‑M1/79‑150 (MQ=255)
                                         tgCCGTTGCCTCGTTCGCGGCGTGTGCTGCCGGTGGCGGCCAGGCCAGCGCGAGCCCGCGAGCACAGCCACGAGGGGGCCAGCGAGTTGCTGGCCGGCATTGCCCGGCTGGCCAATGCCGGTGTCAGCGAGCAGGCCCGTAGCGAATGcc                                                                                                              >  1:478387/1‑150 (MQ=255)
                                         tgCCGTTGCCTCGTTCGCGGCGTGTGCTGCCGGTGGCGGCCAGGCCAGCGCGAGCCCGCGAGCACAGCCACGAGGGGGCCAGCGAGTTGCTGGCCGGCATTGCCCGGCTGGCCAATGCCGGTGCCAGCGAGCAGGCCCGTAGCGAATGcc                                                                                                              >  2:285654/1‑150 (MQ=255)
                                         tgCCGTTGCCTCGTTCGCGGCGTGTGCTGCCGGTGGCGGCCAGGCCAGCGCGAGCCCGCGAGCACAGCCACGAGGGGGCCAGCGAGTTGCTGGCCGGCATTGCCCGGCTGGCCAATGCCGGTGCCAGCGAGCAGGCCCGTAGCGAATGcc                                                                                                              >  2:283074/1‑150 (MQ=255)
                                          actgagcaaacccagccagtagtacacctgtgccgacggcgggtactggctcaggtagcgttggcattcgctacGGGCCAGCGAGTTGCTGGCCGGCATTGCCCGGCTGGCCAATGCCGGTGCCAGCGAGCAGGCCCg                                                                                                                         >  2:1396534‑M1/75‑138 (MQ=255)
                                          actgagcaaacccagccagtagtacacctgtgccgacggcgggtactggctcaggtagcgttggcattcgctacGGGCCAGCGAGTTGCTGGCCGGCATTGCCCGGCTGGCCAATGCCGGTGCCAGCGAGCAGGCCCg                                                                                                                         <  1:1396534‑M1/64‑1 (MQ=255)
                                              ttGCCTCGTTCGCGGCGTGTGCTGCCGGTGGCGGCCAGGCCAGCGCGAGCCCGCGAGCACAGCCACGAGGGGGCCAGCGAGTTGCTGGCCGGCATTGCCCGGCTGGcc                                                                                                                                                   >  2:1861626/1‑108 (MQ=255)
                                               tGCCTCGTTCGCGGCGTGTGCTGCCGGTGGCGGCCAGGCCAGCGCGAGCCCGCGAGCACAGCCACGAGGGGGCCAGCGAGTTGCTGGCCGGCATTGCCCGGCTGGCCAATGCCGGTGCCAGCGAGCAGGCCCGTAGCGAATGCCAACGCt                                                                                                        <  2:3060095/150‑1 (MQ=255)
                                               tGCCGCGGTCGGGGCGTGGGCTGCCGGGGGCGGCCAGGCCAGCGCGAGCCCGCGAGCACAGCCACGAGGGGGCCAGCGAGTTGCTGGCCGGCATGGCCCGGCTGGCCAATGCCGGTGCCAGCGAGCAGGCCCGTAGCGAATGCCAActct                                                                                                        <  2:121074/150‑1 (MQ=255)
                                                 ccTCGTTCGCGGCGTGTGCTGCCGGTGGCGGCCAGGCCAGCGCGAGCCCGCGAGCACAGCCACGAGGGGGCCAGCGAGTTGCTGGCCGGCATTGCCCGGCTGGCCAATGCCGGTGCCAGc                                                                                                                                    <  2:329006/120‑1 (MQ=255)
                                                 ccTCGTTCGCGGCGTGTGCTGCCGGTGGCGGCCAGGCCAGCGCGAGCCCGCGAGCACAGCCACGAGGGGGCCAGCGAGTTGCTGGCCGGCATTGCCCGGCTGGCCAATGCCGGTGCCAGCGAGCAGGCCCGTAGCGAATGCCAACGCTAc                                                                                                      >  1:2713726/1‑150 (MQ=255)
                                                  cTCGTTCGCGGCGTGTGCTGCCGGTGGGGTCCAGGCCAGCGCGAGCCCGCGAGCACAGCCACGAGGGGGCCAGCGAGTTGCTGGCCGGCATTGCCCGGCTGGCCAATGCCGGTGCCAGCGAGCAGGCCCGTAGCGAATGCCAACGCTAcc                                                                                                     >  1:1330264/1‑150 (MQ=255)
                                                  cTCGTTCGCGGCGTGTGCTGCCGGTGGCGGCCAGGCCAGCGCGAGCCCGCGAGCACAGCCACGAGGGGGCCAGCGAGTTGCTGGCCGGCATTGCCCGGCTGGCCAATGCCGGTGCCAGCGAGCAGGCCCGTAGCGAATGCCAACGCTAcc                                                                                                     >  1:668095/1‑150 (MQ=255)
                                                  cTCGTTCGCGGCGTGTGCTGCCGGTGGCGGCCAGGCCAGCGCGAGCCCGCGAGCACAGCCACGAGGGGGCCAGCGAGTTGCTGGCCGGCATTGCCCGGCTGGCCAATGCCGGTGCCAGCGAGCAGGCCCGTAGCGAATGCCAACGCTAcc                                                                                                     >  1:501870/1‑150 (MQ=255)
                                                   tCGTTCGCGGCGTGTGCTGCCGGTGGCGGCCAGGCCAGCGCGAGCCCGCGAGCACAGCCACGAGGGGGCCAGCGAGTTGCTGGCCGGCATTGCCCGGCTGGCCAATGCCGGTGCCAGCGa                                                                                                                                  <  1:2359681/120‑1 (MQ=255)
                                                   tCGTTCGCGGCGTGTGCTGCCGGTGGCGGCCAGGCCAGCGCGAGCCCGCGAGCACAGCCACGAGGGGGCCAGCGAGTTGCTGGCCGGCATTGCCCGGCTGGCCAATGCCGGTGCCAGCGa                                                                                                                                  >  2:2359681/1‑120 (MQ=255)
                                                   tCGTTCGCGGCGTGTGCTGCCGGTGGCGGCCAGGCCAGCGCGAGCCCGCGAGCACAGCCACGAGGGGGCCAGCGAGTTGCTGGCCGGCATTGCCCGGCTGGCCAATGCCGGTGCCAGCGAGCAGGCCCGTAGCGAATGCCAACGCTACCt                                                                                                    >  2:2770547/1‑150 (MQ=255)
                                                      ttCGCGGCGTGTGCTGCCGGTGGCGGCCAGGCCAGCGCGAGCCCGCGAGCACAGCCACGAGGGGGCCAGCGAGTTGCTGGCCGGCATTGCCCGGCTGGCCAATGCCGGTGCCAGCGAGCAGGCCCGTAGCGaa                                                                                                                  >  2:1372008/1‑133 (MQ=255)
                                                      ttCGCGGCGTGTGCTGCCGGTGGCGGCCAGGCCAGCGCGAGCCCGCGAGCACAGCCACGAGGGGGCCAGCGAGTTGCTGGCCGGCATTGCCCGGCTGGCCAATGCCGGTGCCAGCGAGCAGGCCCGTAGCGaa                                                                                                                  <  1:1372008/133‑1 (MQ=255)
                                                         gcggcgTGTGCTGCCGGTGGCGGCCAGGCCAGCGCGAGCCCGCGAGCACAGCCACGAGGGGGCCAGCGAGTTGCTGGCCGGCATTGCCCGGCTGGCCAATGCCGGTGCCAGCGAGCAGGCCCGTAGCGAATGCCAACGCTACCTGAGCCa                                                                                              <  2:501870/150‑1 (MQ=255)
                                                            gcgTGTGCTGCCGGTGGCGGCCAGGCCAGCGCGAGCCCGCGAGCACAGCCACGAGGGGGCCAGCGAGTTGCTGGCCGGCATTGCCCGGCTGGCCAATGCCGGTGCCAGCGAGCAGGCCCGTAGCGAATGCCAACGCTACCTGAGCCAGTa                                                                                           >  2:52871/1‑150 (MQ=255)
                                                              gtgtGCTGCCGGTGGCGGCCAGGCCAGCGCGAGCCCGCGAGCACAGCCACGAGGGGGCCAGCGAGTTGCTGGCCGGCATTGCCCGGCTGGCCAATGCCGGTGCCAGCGAGCCGGCCCGTAGCGAATGCCAACGCTACCTGAGCCAGTAcc                                                                                         <  1:242233/150‑1 (MQ=255)
                                                                 tgctgcCGGTGGCGGCCAGGCCAGCGCGAGCCCGCGAGCACAGCCACGAGGGGGCCAGCGAGTTGCTGGCCGGCATTGCCCGGCTGGCCAATGCCGGTGCCAGCGAGCAGGCCCGTAGCGAATGCCAACGCTACCTGAGCCAGTACccgc                                                                                      >  1:954401/1‑150 (MQ=255)
                                                                    tgcCGGTGGCGGCCAGGCCAGCGCGAGCCCGCGAGCACAGCCACGAGGGGGCCAGCGAGTTGCTGGCCGGCATTGCCCGGCTGGCCAATGCCGGTGCCAGCGAGCAGGCCCGTAGCGAATGCCAACGCTACCTGAGCCAGTACCCGCCGt                                                                                   >  2:2628858/1‑150 (MQ=255)
                                                                      ccGGTGGCGGCCAGGCCAGCGCGAGCCCGCGAGCACAGCCACGAGGGGGCCAGCGAGTTGCTGGCCGGCATTGCCCGGCTGGCCAATGCCGGTGCCAGCGAGCAGGCCCGTAGCGAATg                                                                                                                <  1:1697888/119‑1 (MQ=255)
                                                                         gTGGCGGCCAGGCCAGCGCGAGCCCGCGAGCACAGCCACGAGGGGGCCAGCGAGTTGCTGGCCGGCATTGCCCGGCTGGCCAATGCCGGTGCCAGCGAGCAGGCCCGTAGCGAATGCCAACGCTACCTGAGCCAGTACCCGCCGTCGGca                                                                              >  2:2583543/1‑150 (MQ=255)
                                                                               gccaggccagCGCGAGCCCGCGAGCACAGCCACGAGGGGGCCAGCGAGTTGCTGGCCGGCATTGCCCGGCTGGCCAATGCCGGTGCCAGCGAGCa                                                                                                                               <  1:2860942/95‑1 (MQ=255)
                                                                               gccaggccagCGCGAGCCCGCGAGCACAGCCACGAGGGGGCCAGCGAGTTGCTGGCCGGCATTGCCCGGCTGGCCAATGCCGGTGCCAGCGAGCAGGCCCGTAGCGAATGCCAACGCTACCTGAGCCAGTACCCGCCGTCGGCACAGgtg                                                                        <  2:554946/150‑1 (MQ=255)
                                                                                  aggccagCGCGAGCCCGCGAGCACAGCCACGAGGGGGCCAGCGAGTTACTGGCCGGCATTGCCCGGCTGGCCAATGCCGGTGCCAGCGAGCAGGCCCGTAGCGAATg                                                                                                                <  1:821018/107‑1 (MQ=255)
                                                                                  aggccagCGCGAGCCCGCGAGCACAGCCACGAGGGGGCCAGCGAGTTACTGGCCGGCATTGCCCGGCTGGCCAATGCCGGTGCCAGCGAGCAGGCCCGTAGCGAATg                                                                                                                >  2:821018/1‑107 (MQ=255)
                                                                                     ccagCGCGAGCCCGCGAGCACAGCCACGAGGGGGCCAGCGAGTTGCTGGCCGGCATTGCCCGGCTGGCCAATGCCGGTGCCAGCGAGCAGGCCCGTAGCGAATGCCAACGCTACCTGAGCCAGTACCCGCCGTCGGCACAGGTGtactgc                                                                  >  1:1651097/1‑148 (MQ=255)
                                                                                      cagCGCGAGCCCGCGAGCACAGCCACGCGGGGGCCGGCGCGTCGCTGGCAGGCATTGCCCGGCGGGCCAATGCCGGTGCCAGCGAGCAGGCCCGTAGCGCATGCCCACCCCACCTGCGCCAGTCACCGCCGTCGGCACAGGTGTactact                                                                 <  2:2968349/150‑1 (MQ=255)
                                                                                      cagCGCGAGCCCGCGAGCACAGCCACGAGGGGGCCAGCGAGTTGCTGGCCGGCATTGCCCGGCTGGCCAATGCCGGTGCCAGCGAGCAGGCCCGTAGCGaa                                                                                                                  <  2:2653763/101‑1 (MQ=255)
                                                                                      cagCGCGAGCCCGCGAGCACAGCCACGAGGGGGCCAGCGAGTTGCTGGCCGGCATTGCCCGGCTGGCCAATGCCGGTGCCAGCGAGCAGGCCCGTAGCGaa                                                                                                                  >  1:2653763/1‑101 (MQ=255)
                                                                                      cagCGCGAGCCCGCGAGCACAGCCACGAGGGGGCCAGCGAGTTGCTGGCCGGCATTGCCCGGCTGGCCAATGCCGGTGCCAGCGAGCAGGCCCGTAGCGAATGCCAACGCTACCTGAGCCAGTACCCGCCGTCGGCACAGGTGTactact                                                                 <  2:322960/150‑1 (MQ=255)
                                                                                        gcgcGAGCCCGCGAGCACAGCCACGAGGGGGCCAGCGAGTTGCTGGCCGGCATTGCCCGGCTGGCCAATGCCGGTGCCAGCGAGCAGGCCCGTAGCGAATGCCAACGCTACCTGAGCCAGTACCCGCCGTCGGCACAGGTGTACTActgg                                                               >  2:2251807/1‑150 (MQ=255)
                                                                                          gcgAGCCCGCGAGCACAGCCACGAGGGGGCCAGCGAGTTGCTGGCCGGCATTGCCCGGCTGGCCAATGCCGGTGCCAGCGAGCAGGCCCGTAGCGAATGCCAACGCTACCTGATCCAGTACCCGCCGTCGGCACAGGTGTACTActggct                                                             <  1:283074/150‑1 (MQ=255)
                                                                                          gcgAGCCCGCGAGCACAGCCACGAGGGGGCCAGCGAGTTGCTGGCCGGCATTGCCCGGCTGGCCAATGCCGGTGCCAGCGAGCAGGCCCGTAGCGAATGCCAACGCTACCTGAGCCAGTACCCGCCGTCGGCACAg                                                                           >  1:1222343/1‑136 (MQ=255)
                                                                                            gAGCCCGCGAGCACAGCCACGAGGGGGCCAGCGAGTTGCTGGCCGGCATTGCCCGGCTGGCCAATGCCGGTGCCAGCGAGCAGGCCCGTAGCGAATGCCAACGCTACCTGAGCCAGTACCCGCCGTCGGCACAGGTGTACTActggctgg                                                           <  2:260452/150‑1 (MQ=255)
                                                                                            gAGCCCGCGAGCACAGCCACGAGGGGGCCAGCGAGTTGCTGGCCGGCATTGCCCGGCTGGCCAATGCCGGTGCCAGCGAGCAGGCCCGTAGCGAATGCCAACGCTACCTGAGCCAGTACCCGCCGTCGGCACAGGTGTACTActggctgg                                                           <  1:2583543/150‑1 (MQ=255)
                                                                                             aGCCCGCGAGCACAGCCACGAGGGGGCCAGCGAGTTGCTGGCCGGCATTGCCCGGCTGGCCAATGCCGGTGCCAGCGAGCAGGCCCGTAGCGAATGc                                                                                                               >  1:708719/1‑97 (MQ=255)
                                                                                             aGCCCGCGAGCACAGCCACGAGGGGGCCAGCGAGTTGCTGGCCGGCATTGCCCGGCTGGCCAATGCCGGTGCCAGCGAGCAGGCCCGTAGCGAATGc                                                                                                               <  2:708719/97‑1 (MQ=255)
                                                                                               cccGCGAGCACAGCCACGAGGGGGCCAGCGAGTTGCTGGCCGGCATTGCCCGGCTGGCCAATGCCGGTGCCAGCGAGCAGGCTCGTAGCGAATGCCAACGCTACCTGAGCCAGTACCCGCCGTCGGCACAGGTGTACTACTGGCTGGGtt                                                        <  1:685289/150‑1 (MQ=255)
                                                                                               cccGCGAGCACAGCCACGAGGGGGCCAGCGAGTTGCTGGCCGGCATTGCCCGGCTGGCCAATGCCGGTGCCAGCGAGCAGGCCCGTAGCGAATGCCAACGCTACCTGAGCCAGTACCCGCCGTCGGCACAGGTGTACTACTGGCTGGGtt                                                        >  1:1789730/1‑150 (MQ=255)
                                                                                                ccGCGAGAACAGCCACGAGGGGGCCAGCGAGTTGCTGGCCGGCATTGCCCGGCTGGCCAATGCCGGTGCCAGCGAGCAGGCCCGTAGCGAATGCCAACGCTACCTGAGCCAGTACCCGCCGTCGGCACAGGTGTACTACTGGCTGGGttt                                                       >  1:1952105/1‑150 (MQ=255)
                                                                                                      gCACAGCCACGAGGGGGCCAGCGAGTTGCTGGCCGGCATTGCCCGGCTGGCCAATGCCGGTGCCAGCGAGCAGGCCCGTAGCGAATGCCAACGCTACCTGAGCCAGTACCCGCCGTCGGCACAGGTGTACTACTGGCTGGGTTTGCTCAg                                                 >  1:1254155/1‑150 (MQ=255)
                                                                                                       cacaGCCACGAGGGGGCCAGCGAGTTGCTGGCCGGCATTGCCCGGCTGGCCAATGCCGGTGCCAGCGAGCAGGCCCGTAGCGAATGCCAACGCTACCTGAGCCAGTACCCGCCGTCGGCACAGGTGTACTActggctgg                                                           <  1:2063052/139‑1 (MQ=255)
                                                                                                       cacaGCCACGAGGGGGCCAGCGAGTTGCTGGCCGGCATTGCCCGGCTGGCCAATGCCGGTGCCAGCGAGCAGGCCCGTAGCGAATGCCAACGCTACCTGAGCCAGTACCCGCCGTCGGCACAGGTGTACTActggctgg                                                           >  2:2063052/1‑139 (MQ=255)
                                                                                                       cacaGCCACGAGGGGGCCAGCGAGTTGCTGGCCGGCATTGCCCGGCTGGCCAATGCCGGTGCCAGCGAGCAGGCCCGTAGCGAATGCCAACGCTACCTGAGCCAGTACCCGCCGTCGGCACAGGTGTACTACTGGCTGGGTTTGCTCAGt                                                <  2:954401/150‑1 (MQ=255)
                                                                                                       cacaGCCACGAGGGGGCCAGCGAGTTGCTGGCCGGCATTGCCCGGCTGGCCAATGCCGGTGCCAGCGAGCAGGCCCGTAGCGAATGCCAACGCTACCTGAGCCAGTACCCGCCGTCGGCACAGGTGTACTACTGGCTGGGTTTGCTCAGt                                                <  1:2843235/150‑1 (MQ=255)
                                                                                                            ccACGAGGGGGCCAGCGAGTTGCTGGCCGGCATTGCCCGGCTGGCCAATGCCGGTGTCAGCGAGCAGGCCCGTAGCGAATGCCAACGCTACCTGAGCCAGTACCCGCCGTCGGCACAGGTGTACTACTGGCTGGGTTTGCTCAGTGATAc                                           <  2:478387/150‑1 (MQ=255)
                                                                                                            ccACGAGGGGGCCAGCGAGTTGCTGGCCGGCATTGCCCGGCTGGCCAATGCCGGTGCCAGCGAGCAGGCCCGTa                                                                                                                       <  1:2171195/74‑1 (MQ=255)
                                                                                                            ccACGAGGGGGCCAGCGAGTTGCTGGCCGGCATTGCCCGGCTGGCCAATGCCGGTGCCAGCGAGCAGGCCCGTAGCGAATGCCAACGCTACCTGAGCCAGTACCCGCCGTCGGCACAGGTGTACTACTGGCTGGGTTTGCTCAg                                                 >  2:1610443/1‑144 (MQ=255)
                                                                                                            ccACGAGGGGGCCAGCGAGTTGCTGGCCGGCATTGCCCGGCTGGCCAATGCCGGTGCCAGCGAGCAGGCCCGTAGCGAATGCCAACGCTACCTGAGCCAGTACCCGCCGTCGGCACAGGTGTACTACTGGCTGGGTTTGCTCAg                                                 <  1:1610443/144‑1 (MQ=255)
                                                                                                                gAGGGGGCCAGCGAGTTGCTGGCCGGCATTGCCCGGCTGGCCAATGCCGGTGCCAGCGAGCAGGCCCGTAGCGAATGCCAACGCTACCTGAGCCAGTACCCGCCGTCGGCACAGGTGTACTACTGGCTGGGTTTGCTCAGTGATACCGaa                                       <  2:2713726/150‑1 (MQ=255)
                                                                                                                  gggggCCAGCGAGTTGCTGGCCGGCATTGCCCGGCTGGCCAATGCCGGTGCCAGCGAGCAGGCCCGTAGCGAATGCCAACGCTACCTGAGCCAGTACCCGCCGTCGGCACAGGTGTACTACTGGCTGGGTTTGCTCAGTGATACCGAAgg                                     <  2:2016183/150‑1 (MQ=255)
                                                                                                                    cGGCCAGCAACTCGCTGGCCGGCATTGCCCGGCTGGCCAATGCCGGTGCCAGCGAGCAGGCCCGTAGCGAATGCCAACGCTACCTGAGCCAGTACCCGCCGTCGGCACAGGTGTACTACTGGCTGGGTTTGCTCAGTGATACCGaa                                       >  1:2343997‑M1/2‑146 (MQ=255)
                                                                                                                    cGGCCAGCAACTCGCTGGCCGGCATTGCCCGGCTGGCCAATGCCGGTGCCAGCGAGCAGGCCCGTAGCGAATGCCAACGCTACCTGAGCCAGTACCCGCCGTCGGCACAGGTGTACTACTGGCTGGGTTTGCTCAGTGATACCGaa                                       <  2:2343997‑M1/145‑1 (MQ=255)
                                                                                                                     GGCCAGCGAGTTGCTGGCCGGCATTGCCCGGCTGGCCAATGCCGGTGCCAGCGAGCAGGCCCGTAGCGAATGCCAACGCTACCTGAGCCAGTACCCGCCGTCGGCACAGGTGTACTACTGGCTGGGTTTGCTCAGTGATACCGAAGGCGa                                  <  2:1666409‑M1/150‑1 (MQ=255)
                                                                                                                      cCCAGCGAGGTGCTGGACGGCATTGCCCGGCTGGCCAATGCCGGTGCCAGCGAGCAGGACCGTAGCGAATGCAAACGCTACATGAGCCAGTACCCGCCGTCGGCACAGGTGTACTACTGGCTGGGTTTGCTCAGTGATACCGAAGGCGAt                                 >  2:1413415‑M1/2‑150 (MQ=255)
                                                                                                                      GCCAGCGAGTTGCTGGCCGGCATTGCCCGGCTGGCCAATGCCGGTGCCAGCGAGCAGGCCCGTAGCGAATGCCAACGCTACCTGAGCCAGTACCCGCCGTCGGCACAGGTGTACTACTGGCTGGGTTTGCTCAGTGATACCGAAGGCGAt                                 >  1:1219301‑M1/1‑150 (MQ=255)
                                                                                                                       CCAGCGAGTTGCTGGCCGGCATTGCCCGGCTGGCCAATGCCGGTGCCAGCGAGCAGGCCCGTAGCGAATGCCAACGCTACCTGAGCCAGTACCCGCCGTCGGCACAGGTGtactac                                                                  >  1:1892731‑M1/1‑116 (MQ=255)
                                                                                                                       CCAGCGAGTTGCTGGCCGGCATTGCCCGGCTGGCCAATGCCGGTGCCAGCGAGCAGGCCCGTAGCGAATGCCAACGCTACCTGAGCCAGTACCCGCCGTCGGCACAGGTGtactac                                                                  <  2:1892731‑M1/116‑1 (MQ=255)
                                                                                                                        CAGCGAGTTGCTGGCCGGCATTGCCCGGCTGGCCAATGCCGGTGCCAGCGAGCAGGCCCGTAGCGAATGCCAACGCTACCTGAGCCAGTACCCGCCGTCGGCACAGGTGTACTActggc                                                              >  1:2266719‑M1/1‑119 (MQ=255)
                                                                                                                          gCGAGTTGCTGGCCGGCATTGCCCGGCTGGCCAATGCCGGTGCCAGCGAGCATGACCGTAGCGAATGCCAACGCTACCTGAGCCAGTACCCGCCGTCGGCACAGGTGTAGTACTGGCTGGGTTTGCTCAGTGATACCGAAGGCGATgcgc                             >  1:363466/1‑150 (MQ=255)
                                                                                                                            gAGTTGCTGGCCGGCATTGCCCGGCTGGCCAATGCCGGTGCCAGCGAGCAGGCCCGTAGCGAATGCCAAAGCTACCTGAGCCAGTACCCGCCGTCGGCACAGGTGTACTACTGGCTGGGTTTGCTCAGTGATACCGAAGGCGATGCgcag                           <  2:2257493/150‑1 (MQ=255)
                                                                                                                              gTTGCTGGCCGGCATTGCCCGGCTGGCCAATGCCGGTGCCAGCGAGCAGGCCCGTAGCGAATGCCAACGCTACCTGAGCCAGTACCCGCCGTCGGCACAGGTGTACTACTGGCTGGGTTTGCTCAGTGATACCGAAGGCGATGCgcagca                         >  1:2785336/1‑150 (MQ=255)
                                                                                                                                   tGGCCGGCATTGCCCGGCTGGCCAATGCCGGTGCCAGCGAGCAGGCCCGTAGCGAATGCCAACGCTACCTGAGCCAGTACCCGCCGTCGGCACAGGTGTACTACTGGCTGGGTTTGCTCAGTGATACCGAAGGCGATGCGCAGCAGgcgc                    <  1:2381226/150‑1 (MQ=255)
                                                                                                                                   tGGCCGGCATTGCCCGGCTGGCCAATGCCGGTGCCAGCGAGCAGGCCCGTAGCGAATGCCAACGCTACCTGAGCCAGTACCCGCCGTCGGCACAGGTGTACTACTGGCTGGGTTTGCTCAGTGATACCGAAGGCGATGCGCAGCAGgcgc                    <  1:1167424/150‑1 (MQ=255)
                                                                                                                                     tttggCATTGCCCGGCTGGCCAATGCCGGTGCCAGCGAGCAGGCCCGTAGCGAATGCCAACGCTACCTGAGCCAGTACCCGCCGTCGGCACAGGTGTACTACTGGCTGGGTTTGCTCAGTGATACCGAAGGCGATGCGCAGCAGGCGCTg                  >  2:443654/4‑150 (MQ=255)
                                                                                                                                     accGGCATTGGCCAGCCGGGCAATGCCGGTGCCAGCGAGCAGGCCCGTAGCGAATGCCAACGCTACCTGAGCCAGTACCCGCCGTCGGca                                                                              <  1:2029110/89‑1 (MQ=38)
                                                                                                                                     accGGCATTGGCCAGCCGGGCAATGCCGGTGCCAGCGAGCAGGCCCGTAGCGAATGCCAACGCTACCTGAGCCAGTACCCGCCGTCGGCACAGGTGta                                                                      >  2:330992/2‑98 (MQ=38)
                                                                                                                                     accGGCATTGGCCAGCCGGGCAATGCCGGTGCCAGCGAGCAGGCCCGTAGCGAATGCCAACGCTACCTGAGCCAGTACCCGCCGTCGGCACAGGTGta                                                                      <  1:330992/97‑1 (MQ=38)
                                                                                                                                      ccGGCATTGGCCAGCCGGGCAATGCCGGTGCCAGCGAGCAGGCCCGTAGCGAATGCCAACGCTACCTGAGCCAGTACCCGCCGTCGGCACAGGTGTACTActgg                                                               >  1:330480/1‑104 (MQ=38)
                                                                                                                                           tttGCGCGGCTGGCCAATGCCGGTGCGCGGGACCGGGCCCGTAGCGAATGGCAACGCTACCTGAGCCAGGACCCGCCGTCGGCACAGGTGTACTACTGGCTGGGTTCGGTCAGTGATACCGAAGGCGATGCCCAGCGGCCGGTGAGTCAc            <  2:125926/149‑1 (MQ=255)
                                                                                                                                               cccGGCTGGCCAATGCCGGTGCCAGCGAGCAGGCCCGTAGCGAATGCCAACGCTACCTGAGCCAGTACCCGCCGTCGGCACAGGTGTACTACTGGCTGGGTTTGCTCACTGATACCGAAGGCGATGCGCAGCAGGCGCTGAGcagatcgg        >  1:2187876/1‑143 (MQ=255)
                                                                                                                                                ccGGCTGGCCAATGCCGGTGCCAGCGAGCAGGCCCGTAGCGAATGCCAACGCTACCTGAGCCAGTACCCGCCGTCGGCACAGGTGTACTACTGGCTGGGTTTGCTCAg                                                 >  2:1059601/1‑108 (MQ=255)
                                                                                                                                                ccGGCTGGCCAATGCCGGTGCCAGCGAGCAGGCCCGTAGCGAATGCCAACGCTACCTGAGCCAGTACCCGCCGTCGGCACAGGTGTACTACTGGCTGGGTTTGCTCAg                                                 <  1:1059601/108‑1 (MQ=255)
                                                                                                                                                  ggCTGGCCAATGCCGGTGCCAGCGAGCAGGCCCGTAGCGAATGCCAACGCTACCTGAGCCAGTACCCGCCGTCGGCACAGGTGTACTACTGGCTGGGTTTGCTCAGTGATACCGAAGGCGATGCGCAGCAGGCGCTGAGCCACTATCGCa     >  1:1979023/1‑150 (MQ=255)
                                                                                                                                                  ggCTGGCCAATGCCGGTGCCAGCGAGCAGGCCCGTAGCGAATGCCAACGCTACCTGAGCCAGTACCCGCCGTCGGCACAGGTGTACTACTGGCTGGGTTTGCTCAGTGATACCGAAGGCGATGCGCAGCAGGCGCTGAGCCACTATCGCa     <  1:732033/150‑1 (MQ=255)
                                                                                                                                                     |                                                                                                                                                     
GTTCACGACGTTGCCCGCACCGGTTTATCCGCAGCCCAGTGTGCCGTTGCCTCGTTCGCGGCGTGTGCTGCCGGTGGCGGCCAGGCCAGCGCGAGCCCGCGAGCACAGCCACGAGGGGGCCAGCGAGTTGCTGGCCGGCATTGCCCGGCTGGCCAATGCCGGTGCCAGCGAGCAGGCCCGTAGCGAATGCCAACGCTACCTGAGCCAGTACCCGCCGTCGGCACAGGTGTACTACTGGCTGGGTTTGCTCAGTGATACCGAAGGCGATGCGCAGCAGGCGCTGAGCCACTATCGCAAGG  >  NC_002947/1694919‑1695217

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