Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 1,700,011 C→G 6.0% 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.0% 1010.7 / 39.1 298A329G (GCC→GGC) wspCbiofilm formation methyltransferase WspC
Reads supporting (aligned to +/- strand):  ref base C (129/151);  new base G (9/9);  total (138/160)
Fisher's exact test for biased strand distribution p-value = 8.10e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.98e-01

Only 100 of 307 total aligned reads displayed.

CGTGCGTTCACGACGTTGCCCGCACCGGTTTATCCGCAGCCCAGTGTGCCGTTGCCTCGTTCGCGGCGTGTGCTGCCGGTGGCGGCCAGGCCAGCGCGAGCCCGCGAGCACAGCCACGAGGGGGCCAGCGAGTTGCTGGCCGGCATTGCCCGGCTGGCCAATGCCGGTGCCAGCGAGCAGGCCCGTAGCGAATGCCAACGCTACCTGAGCCAGTACCCGCCGTCGGCACAGGTGTACTACTGGCTGGGTTTGCTCAGTGATACCGAAGGCGATGCGCAGCAGGCGCTGAGCCACTATC  >  NC_002947/1699863‑1700160
                                                                                                                                                    |                                                                                                                                                     
 gtgtgtTCACGACGTTGCCCGCACCGGTTTATCCGCAGCCCAGTGTGCCGTTGCCTCGTTCGCGGCGTGTGCTGCCGGTGGCGGCCAGGCCAGCGCGAGCCCGCGAGCACAGCCACGAGGGGGCCAGCGAGTTGCTGGCCGGCATTGccc                                                                                                                                                     <  1:1004249/150‑1 (MQ=255)
   gcgtTCACGACGTTGCCCGCACCGGTTTATCCGCAGCCCAGTGTGCCGTTGCCTCGTTCGCGGCGTGTGCTGCCGGTGGCGGCCAGGCCAGCGCGAGCCCGCGAGCACAGCCACGAGGGGGCCAGCGAGTTGCTGGCCGGCATTGCCCgg                                                                                                                                                   <  1:1391807/150‑1 (MQ=255)
     gtTCACGACGTTGCCCGCACCGGTTTATCCGCAGCCCAGTGTGCCGTTGCCTCGTTCGCGGCGTGTGCTGCCGGTGGCGGCCAGGCCAGCGCGAGCCCGCGAGCACAGCCACGAGGGGGCCAGCAAGTTGCTGGCCGGCATTGCCCGGCt                                                                                                                                                 <  2:1509000/150‑1 (MQ=255)
       tCACGACGTTGCCCGCACCGGTTTATCCGCAGCCCAGTGTGCCGTTGCCTCGTTCGCGGCGTGTGCTGCCGGTGGCGGCCAGGCCAGCGCGAGCCCGCGAGCACAGCCACGAGGGGGCCAGCGAGTTGCTGGCCGGCATTGCCCGGCTgg                                                                                                                                               <  2:2248292/150‑1 (MQ=255)
        cACGACGTTGCCCGCACCGGTTTATCCGCAGCCCAGTGTGCCGTTGCCTCGTTCGCGGCGTGTGCTGCCGGTGGCGGCCAGGCCAGCGCGAGCCCGCGAGCACAGCCACGAGGGGGCCAGCGAGTTGCTGGCCGGCATTGCCCGGCTGGc                                                                                                                                              <  1:4634122/150‑1 (MQ=255)
         acgacgTTGCCCGCACCGGTTTATCCGCAGCCCAGTGTGCCGTTGCCTCGTTCGCGGCGTGTGCTGCCGGTGGCGGCCAGGCCAGCGCGAGCCCGCGAGCACAGCCACGAGGGGGCCAGCGAGTTGCTGGCCGGCATTGCCCGGCTGGcc                                                                                                                                             >  2:1310740/1‑150 (MQ=255)
            acgTTGCCCGCACCGGTTTATCCGCAGCCCAGTGTGCCGTTGCCTCGTTCGCGGCGTGTGCTGCCGGTGGCGGCCAGGCCAGCGCGAGCCCGCGAGCACAGCCACGAGGGGGCCAGCGAGTTGCTGGCCGGCATTGCCCGGCTGGCCAAt                                                                                                                                          <  1:1068095/150‑1 (MQ=255)
             cgTTGCCCGCACCGGTTTATCCGCAGCCCAGTGTGCCGTTGCCTCGTTCGCGGCGTGTGCTGCCGGTGGCGGCCAGGCCAGCGCGAGCCCGCGAGCACAGCCACGAGGGGGCCAGCGAGTTGCTGACCGGCATTGGCCAGCCGGGCAATg                                                                                                                                         <  2:2352670/150‑1 (MQ=255)
               ttGCCCGCACCGGTTTATCCGCAGCCCAGTGTGCCGTTGCCTCGTTCGCGGCGTGTGCTGCCGGTGGCGGCCAGGCCAGCGCGAGCCCGCGAGCACAGCCACGAGGGGGCCAGCGAGTTGCTGGCCGGCATTGCCCGGCTGGCCAATGcc                                                                                                                                       <  1:1800306/150‑1 (MQ=255)
               ttGCCCGCACCGGTTTATCCGCAGCCCAGTGTGCCGTTGCCTCGTTCGCGGCGTGTGCTGCCGGTGGCGGCCAGGCCAGCGCGAGCCCGCGAGCACAGCCACGAGGGGGCCAGCGAGTTGCTGGCCGGCATTGCCCGGCTGGCCAATGcc                                                                                                                                       >  1:417229/1‑150 (MQ=255)
               ttGCCCGCACCGGTTTATCCGCAGCCCAGTGTGCCGTTGCCTCGTTCGCGGCGTGTGCTGCCGGTGGCGGCCAGGCCAGCGCGAGCCCGCGAGCACAGCCACGAGGGGGCCAGCGAGTTGCTGGCCGGCATTGCCCGGCTGGCCAATGcc                                                                                                                                       >  1:2170473/1‑150 (MQ=255)
                tGCCCGCACCGGTTTATCCGCAGCCCAGTGTGCCGTTGCCTCGTTCGCGGCGTGTGCTGCCGGTGGCGGCCAGGCCAGCGCGAGCCCGCGAGCACAGCCACGAGGGGGCCAGCGAGTTGCTGGCCGGCATTGGCCAGCCGGGCAATgact                                                                                                                                      <  2:2438866/150‑4 (MQ=255)
                   ccGCACCGGTTTATCCGCAGCCCAGTGTGCCGTTGCCTCGTTCGCGGCGTGTGCTGCCGGTGGCGGCCAGGCCAGCGCGAGCCCGCGAGCACAGCCACGAGGGGGCCAGCGAGTTGCTGGCCGGCATTGCCCGGCTGGCCAATGCCGGTg                                                                                                                                   >  1:1744414/1‑150 (MQ=255)
                   ccGCACCGGTTTATCCGCAGCCCAGTGTGCCGTTGCCTCGTTCGCGGCGTGTGCTGCCGGTGGCGGCCAGGCCAGCGCGAGCCCGCGAGCACAGCCACGAGGGGGCCAGCGAGTTGCTGGCCGGCATTGCCCGGCTGGCCAATGCCGGTg                                                                                                                                   >  1:4868519/1‑150 (MQ=255)
                    cGCACCGGTTTATCCGCAGCCCAGTGTGCCGTTGCCTCGTTCGCGGCGTGTGCTGCCGGTGGCGGCCAGGCCAGCGCGAGCCCGCGAGCACAGCCACGAGGGGGCCAGCGAGTTGCTGGCCGGCATTGCCCGGCTGGCCAATGCCGGTGc                                                                                                                                  <  1:1544669/150‑1 (MQ=255)
                       aCCGGTTTATCCGCAGCCCAGTGTGCCGTTGCCTCGTTCGCGGCGTGTGCTGCCGGTGGCGGCCAGGCCAGCGCGAGCCCGCGAGCACAGCCACGAGGGGGCCAGCGAGTTGCTGGCCGGCATTGCCCGGCTGGCCAATGCCGGTGCCAg                                                                                                                               >  1:1685481/1‑150 (MQ=255)
                        ccGGTTTATCCGCAGCCCAGTGTGCCGTTGCCTCGTTCGCGGCGTGTGCTGCCGGTGGCGGCCAGGCCAGCGCGAGCCCGCGAGCACAGCCACGAGGGGGCCAGCGAGTTGCTGGCCGGCATTGCCCGGCTGGCCAATGCCGGTGCCAGc                                                                                                                              <  1:2986515/150‑1 (MQ=255)
                        ccGGTTTATCCGCAGCCCAGTGTGCCGTTGCCTCGTTCGCGGCGTGTACTGCCGGTGGCGGCCAGGCCAGCGCGAGCCCGCGAGCACAGCCACGAGGGGGCCAGCGAGTTGCTGGCCGGCATTGCCCGGCTGGCCAATGCCGGTGCCAGc                                                                                                                              <  1:5460233/150‑1 (MQ=255)
                        ccGGTTTATCCGCAGCCCAGTGTGCCGGTGCCTCGTTCGCGGCGTGTGCTGCCGGTGGCGGCCAGGCCAGCGCGAGCCCGCGAGCACAGCCACGAGGGGGCCAGCGAGGTGCTGGCCGGCATTGCCCGGCTGGCCAATGCCGGTGCCAGc                                                                                                                              <  2:3091505/150‑1 (MQ=255)
                          ggTTTATCCGCAGCCCAGTGTGCCGTTGCCTCGTTCGCGGCGTGTGCTGCCGGTGGCGGCCAGGCCAGCGCGAGCCCGCGAGCACAGCCACGAGGGGGCCAGCGAGTTGCTGGCCGGCATTGCCCGGCTGGCCAATGCCGGTGCCAGCGa                                                                                                                            <  1:2245386/150‑1 (MQ=255)
                           gTTTATCCGCAGCCCAGTGTGCCGTTGCCTCGTTCGCGGCGTGTGCTGCCGGTGGCGGCCAGGCCAGCGCGAGCCCGCGAGCACAGCCACGAGGGGGCCAGCGAGTTGCTGGCCGGCATTGCCCGGCTGGCCAATGCCGGTGCCAGccgg                                                                                                                           <  1:982150/150‑4 (MQ=255)
                              tATCCGCAGCCCAGTGTGCCGTTGCCTCGTTCGCGGCGTGTGCTGCCGGTGGCGGCCAGGCCAGCGCGAGCCCGCGAGCACAGCCACGAGGGGGCCAGCGAGTTGCTGGCCGGCATTGGCCAGCCGGGCAATGCCg                                                                                                                                      <  2:2434155/136‑1 (MQ=255)
                              tATCCGCAGCCCAGTGTGCCGTTGCCTCGTTCGCGGCGTGTGCTGCCGGTGGCGGCCAGGCCAGCGCGAGCCCGCGAGCACAGCCACGAGGGGGCCAGCGAGTTGCTGGCCGGCATTGGCCAGCCGGGCAATGCCg                                                                                                                                      >  1:2434172/1‑136 (MQ=255)
                                tCCGCAGCCCAGTGTGCCGTTGCCTCGTTCGCGGCGTGTGCTGCCGGTGGCGGCCAGGCCAGCGCGAGCCCGCGAGCACAGCCACGAGGGGGCCAGCGAGTTGCTGGCCGGCATTGCCCGGCTGGCCAATGCCGGTGCCAGCGAGCAGGc                                                                                                                      <  2:4868523/150‑1 (MQ=255)
                                  cGCAGCCCAGTGTGCCGTTGCCTCGTTCGCGGCGTGTGCTGCCGGTGGCGGCCAGGCCAGCGCGAGCCCGCGAGCACAGCCACGAGGGGGCCAGCGAGTTGCTGGCCGGCATTGCCCGGCTGGCCAATGCCGGTGCCAGCGAGCAGGccc                                                                                                                    <  1:277702/150‑1 (MQ=255)
                                   gCAGCCCAGTGTGCCGTTGCCTCGTTCGCGGCGTGTGCTGCCGGTGGCGGCCAGGCCAGCGCGAGCCCGCGAGCACAGCCACGAGGGGGCCAGCGAGTTGCTGGCCGGCATTGCCCGGCTGGCCAATGCCGGTGCCAGCGAGCAGGCCCg                                                                                                                   >  2:103561/1‑150 (MQ=255)
                                     aGCCCAGTGTGCCGTTGCCTCGTTCGCGGCGTGTGCTGCCGGTGGCGGCCAGGCCAGCGCGAGCCCGCGAGCACAGCCACGAGGGGGCCAGCGAGTTGCTGGCCGGCATTGCCCGGCTGGCCAATGCCGGTGCCAGCGAGCAGGCCCGTa                                                                                                                 <  2:5244905/150‑1 (MQ=255)
                                     aGCCCAGTGTGCCGTTGCCTCGTTCGCGGCGTGTGCTGCCGGTGGCGGCCAGGCCAGCGCGAGCCCGCGAGCACAGCCACGAGGGGGCCAGCGAGTTGCTGGCCGGCATTGCCCGGCTGGCCAATGCCGGTGCCAGCGAGCAGGCCCGTa                                                                                                                 <  1:4921275/150‑1 (MQ=255)
                                      gCCCAGTGTGCCGTTGCCTCGTTCGCGGCGTGTGCTGCCGGTGGCGGCCAGGCCAGCGCGAGCCCGCGAGCACAGCCACGAGGGGGCCAGCGAGTTGCTGGCCGGCATTGCCCGGCTGGCCAATGCCGGTGCCAGCGAGCAGGCCCGTAg                                                                                                                <  1:1510851/150‑1 (MQ=255)
                                       cccAGTGTGCCGTTGCCTCGTTCGCGGCGTGTGCTGCCGGTGGCGGCCAGGCCAGCGCGAGCCCGCGAGCACAGCCACGAGGGGGCCAGCGAGTTGCTGGCCGGCATTGCCCGGCTGGCCAATGCCGGTGCCAGCGAGCAGGCCCGTAGc                                                                                                               <  2:4547567/150‑1 (MQ=255)
                                       cccAGTGTGCCGTTGCCTCGTTCGCGGCGTGTGCTGCCGGTGGCGGCCAGGCCAGCGCGAGCCCGCGAGCACAGCCACGAGGGGGCCAGCGAGTTGCTGGCCGGCATTGCCCGGCTGGCCAATGCCGGTGCCAGCGAGCAGGCCCGTAGc                                                                                                               <  1:3727432/150‑1 (MQ=255)
                                            tgtgCCGTTGCCTCGTTCGCGGCGTGTGCTGCCGGTGGCGGCCAGGCCAGCGCGAGCCCGCGAGCACAGCCACGAGGGGGCCAGCGAGTTGCTGGCCGGCATTGCCCGGCTGGCCAATGCCGGTGCCAGCGAGCAGGCCCGTAGCGAATg                                                                                                          >  2:2876731/1‑150 (MQ=255)
                                                ccGGTGCCTCGTTCGCGGCGTGTGCTGCCGGTGGCGGCCAGGCCAGCGCGAGCCCGCGAGCACAGCCACGAGGGGGCCAGCGAGTTGCTGGCCGGCAGTGCCCGGCTGGCCAATGCCGGTGCCAGCGAGCAGGCCCGTAGCGAATGCCaa                                                                                                      >  1:1146638/1‑150 (MQ=255)
                                                 cGTTGCCTCGTTCGCGGCGTGTGCTGCCGGTGGCGGCCAGGCCAGCGCGAGCCCGCGAGCACAGCCACGAGGGGGCCAGCGAGTTGCTGGCCGGCATTGCCCGGCTGGCCAATGCCGGTGCCAGCGAGCAGGCCCGTAGCGAATGCCAAc                                                                                                     <  1:4762569/150‑1 (MQ=255)
                                                 cGTTGCCTCGTTCGCGGCGTGTGCTGCCGGTGGCGGCCAGGCCAGCGCGAGCCCGCGAGCACAGCCACGAGGGGGCCAGCGAGTTGCTGGCCGGCATTGCCCGGCTGGCCAATGCCGGTGCCAGCGAGCAGGCCCGTAGCGAATGCCAAc                                                                                                     >  1:2565919/1‑150 (MQ=255)
                                                   ttGCCTCGTTCGCGGCGTGTGCTGCCGGTGGCGGCCAGGCCAGCGCGAGCCCGCGAGCACAGCCACGAGGGGGCCAGCGAGTTGCTGGCCGGCATTGCCCGGCTGGCCAATGCCGGTGCCAGCGAGCAGGCCCGTAGCGAATGCCAACGc                                                                                                   >  1:1933577/1‑150 (MQ=255)
                                                       cTCGTTCGCGGCGTGTGCTGCCGGTGGCGGCCAGGCCAGCGCGAGCCCGCGAGCACAGCCACGAGGGGGCCAGCGAGTTGCTGGCCGGCATTGCCCGGCTGGCCAATGCCGGTGCCAGCGAGCAGGCCCGTAGCGAATGCCAACGCTAcc                                                                                               >  1:2368750/1‑150 (MQ=255)
                                                        tCGTTCGCGGCGTGTGCTGCCGGTGGCGGCCAGGCCAGCGCGAGCCCGCGAGCACAGCCACGAGGGGGCCAGCGAGTTGCTGGCCGGCATTGCCCGGCTGGCCAATGCCGGTGCCAGCGAGCAGGCCCGTAGCGAATGCCAACGCTACCt                                                                                              >  2:3741039/1‑150 (MQ=255)
                                                          gTTCGCGGCGTGTGCTGCCGGTGGCGGCCAGGCCAGCGCGAGCCCGCGAGCACAGCCACGAGGGGGCCAGCGAGTTGCTGGCCGGCATTGCCCGGCTGGCCAATGCCGGTGCCAGCGAGCAGGCCCGTAGCGAATGCCAACGCTACCTGa                                                                                            <  2:5054705/150‑1 (MQ=255)
                                                          gTTCGCGGCGTGTGCTGCCGGTGGCGGCCAGGCCAGCGCGAGCCCGCGAGCACAGCCACGAGGGGGCCAGCGAGTTGCTGGCCGGCATTGCCCGGCTGGCCAATGCCGGTGCCAGCGAGCAGGCCCGTAGCGAATGCCAACGCTACCTGa                                                                                            >  2:1832476/1‑150 (MQ=255)
                                                          gTTCGCGGCGTGTGCTGCCGGTGGCGGCCAGGCCAGCGCGAGCCCGCGAGCACAGCCACGAGGGGGCCAGCGAGTTGCTGGCCGGCATTGCCCGGCTGGCCAATGCCGGTGCCAGCGAGCAGGCCCGTAGCGAATGCCAACGCTACCTGa                                                                                            <  2:196139/150‑1 (MQ=255)
                                                           ttCGCGGCGTGTGCTGCCGGTGGCGGCCAGGCCAGCGCGAGCCCGCGAGCACAGCCACGAGGGGGCCAGCGAGTTGCTGGCCGGCATTGCCCGGCTGGCCAATGCCGGTGCCAGCGAGCAGGCCCGTAGCGAATGCCAACGCTACCTGAg                                                                                           >  1:437519/1‑150 (MQ=255)
                                                            tCGCGGCGTGTGCTGCCGGTGGCGGCCAGGCCAGCGCGAGCCCGCGAGCACAGCCACGAGGGGGCCAGCGAGTTGCTGGCCGGCATTGCCCGGCTGGCCAATGCCGGTGCCAGCGAGCAGGCCCGTAGCGAATGCCAACGCTACCTGAGc                                                                                          <  2:1950094/150‑1 (MQ=255)
                                                            tCGCGGCGTGTGCTGCCGGTGGCGGCCAGGCCAGCGCGAGCCCGCGAGCACAGCCACGAGGGGGCCAGCGAGTTGCTGGCCGGCATTGCCCGGCTGGCCAATGCCGGTGCCAGCGAGCAGGCCCGTAGCGAATGCCAACGCTACCTGAGc                                                                                          <  2:3640957/150‑1 (MQ=255)
                                                            tCGCGGCGTGTGCTGCCGGTGGCGGCCAGGCCAGCGCGAGCCCGCGAGCACAGCCACGAGGGGGCCAGCGAGTTGCTGGCCGGCACTGCCCGGCTGGCCAATGCCGGTGCCAGCGAGCAGGCCCGTAGCGAATGCCAACGCTACCTGAGc                                                                                          <  1:3200607/150‑1 (MQ=255)
                                                              gcggcgTGTGCTGCCGGTGGCGGCCAGGCCAGCGCGGGCCCGCGAGCACAGCCGCGAGGGGGCCCGCGAGTTGCTGGCCGGCCTGGCCCGGCTGGCCCATGCCGGTGCCCGCGAGCCGGCCCGTCGCGCATGCCCACGCTACCTGAGccc                                                                                        >  1:3570257/1‑149 (MQ=255)
                                                              gcggcgTGTGCTGCCGGTGGCGGCCAGGCCAGCGCGAGCCCGCGAGCACAGCCACGAGGGGGCCAGCGAGTTGCTGGCCGGCATTGCCCGGCTGGCCAATGCCGGTGCCAGCGAGCAGGCCCGTAGCGAATGCCAACGCTACCTGAGCCa                                                                                        >  2:3398549/1‑150 (MQ=255)
                                                                  cgTGTGCTGCCGGTGGCGGCCAGGCCAGCGCGAGCCCGCGAGCACAGCCACGAGGGGGCCAGCGAGTTGCTGGCCGGCATTGCCCGGCTGGCCAATGCCGGTGCCAGCGAGCAGGCCCGTAGCGAATGCCAACGCTACCTGAGCCAGTAc                                                                                    <  2:3868199/150‑1 (MQ=255)
                                                                  cgTGTGCTGCCGGTGGCGGCCAGGCCAGCGCGAGCCCGCGAGCACAGCCACGAGGGGGCCAGCGAGTTGCTGGCCGGCATTGCCCGGCTGGCCAATGCCGGTGCCAGCGAGCAGGCCCGTAGCGAATGCCAACGCTACCTGAGCCAGTAc                                                                                    >  2:2108271/1‑150 (MQ=255)
                                                                     gtgCTGCCGGTGGCGGCCAGGCCAGCGCGAGCCCGCGAGCACAGCCACGAGGGGGCCAGCGAGTTGCTGGCCGGCATTGCCCGGCTGGCCAATGCCGGTGCCAGCGAGCAGGCCCGTAGCGAATGCCAACGCTACCTGAGCCAGTACccg                                                                                 >  2:1166833/1‑150 (MQ=255)
                                                                      tgctgcCGGTGGCGGCCAGGCCAGCGCGAGCCCGCGAGCACAGCCACGAGGGGGCCAGCGAGTTGCTGGCCGGCATTGCCCGGCTGGCCAATGCCGGTGCCAGCGAGCAGGCCCGTAGCGAATGCCAACGCTACCTGAGCCAGTACccgc                                                                                <  1:385641/150‑1 (MQ=255)
                                                                      tgctgcCGGTGGCGGCCAGGCCAGCGCGAGCCCGCGAGCACAGCCACGAGGGGGCCAGCGAGTTGCTGGCCGGCATTGCCCGGCTGGCCAATGCCGGTGCCAGCGAGCAGGCCCGTAGCGAATGCCAACGCTACCTGAGCCAGTACccgc                                                                                <  2:4746607/150‑1 (MQ=255)
                                                                      tgctgcCGGTGGCGGCCAGGCCAGCGCGAGCCCGCGAGCACAGCCACGAGGGGGCCAGCGAGTTGCTGGCCGGCATTGCCCGGCTGGCCAATGCCGGTGCCAGCGAGCAGGCCCGTAGCGAATGCCAACGCTACCTGAGCCAGTACccgc                                                                                <  1:3406106/150‑1 (MQ=255)
                                                                        ctgcCGGTGGCGGCCAGGCCAGCGCGAGCCCGCGAGCACAGCCACGAGGGGGCCAGCGAGTTGCTGGCCGGCATTGCCCGGCTGGCCAATGCCGGTGCCAGCGAGCAGGCCCGTAGCGAATGCCAACGCTACCTGAGCCAGTACccgccg                                                                              >  1:3473272/1‑150 (MQ=255)
                                                                            cGGTGGCGGCCAGGCCAGCGCGAGCCCGCGAGCACAGCCACGAGGGGGCCAGCGAGTTGCTGGCCGGCATTGCCCGGCTGGCCAATGCCGGTGCCAGCGAGCAGGCCCGTAGCGAATGCCAACACTACCTGAGCCAGTACCCGCCGTCgg                                                                          <  2:4728374/150‑1 (MQ=255)
                                                                                ggcggcCAGGCCAGCGCGAGCCCGCGAGCACAGCCACGAGGGGGCCAGCGAGTTGCTGGCCGGCATTGCCCGGCTGGCCAATGCCGGTGCCAGCGAGCAGGCCCGTAGCGAATGCCAACGCTACCTGAGCCAGTACCCGCCGTCGGcaca                                                                      <  1:568801/150‑1 (MQ=255)
                                                                                    gccaggccagCGCGAGCCCGCGAGCACAGCCACGAGGGGGCCAGCGAGTTGCTGGCCGGCATTGCCCGGCTGGCCAATGCCGGTGCCAGCGAGCAGGCCCGTAGCGAATGCCAACGCTACCTGAGCCAGTACCCGCCGTCGGCACAGgtg                                                                  >  2:4359949/1‑150 (MQ=255)
                                                                                      caggccagCGCGAGCCCGCGAGCACAGCCACGAGGGGGCCAGCGAGTTGCTGGCCGGCATTGCCCGGCTGGCCAATGCCGGTGCCAGCGAGCAGGCCCGTAGCGAATGCCAACGCTACCTGAGCCAGTACCCGCCGTCGGCACAGGTGta                                                                >  1:4959759/1‑150 (MQ=255)
                                                                                         gccagCGCGAGCCCGCGAGCACAGCCACGAGGGGGCCAGCGAGTTGCTGGCCGGCATTGCCCGGCTGGCCAATGCCGGTGCCAGCGAGCAGGCCCGTAGCGAATGCCAACGCTACCTGAGCCAGTACCCGCCGTCGGCACAGGTGtacta                                                             >  1:3885477/1‑150 (MQ=255)
                                                                                           cagCGCGAGCCCGCGAGCACAGCCACGAGGGGGCCAGCGAGTTGCTGGCCGGCATTGCCCGGCTGGCCAATGCCGGTGCCAGCGAGCAGGCCCGTAGCGAATGCCAACGCTACCTGAGCCAGTACCCGCCGTCGGCACAGGTGTactact                                                           <  1:5328727/150‑1 (MQ=255)
                                                                                            agCGCGAGCCCGCGAGCACAGCCACGAGGGGGCCAGCGAGTTGCTGGCCGGCATTGCCCGGCTGGCCAATGCCGGTGCCAGCGAGCAGGCCCGTAGCGAATGCCAACGCTACCTGAGCCAGTACCCGACGTCGGCACAGGTGTACTActg                                                          >  2:905070/1‑150 (MQ=255)
                                                                                             gcgcGAGCCCGCGAGCACAGCCACGAGGGGGCCAGCGAGTTGCTGGCCGGCATTGCCCGGCTGGCCAATGCCGGTGCCAGCGAGCAGGCCCGTAGCGAATGCCAACGCTACCTGAGCCAGTACCCGCCGTCGGCACAGGTGTACTActgg                                                         >  2:1095916/1‑150 (MQ=255)
                                                                                              ctggcccccGCGAGCACAGCCACGAGGGGGCCAGCGAGTTGCTGGCCGGCATTGCCCGGCTGGCCAATGCCGGTGCCAGCGAGCAGGCCCGTAGCGAATGCCAACGCTACCTGAGCCAGTACCCGCCGTCGGCACAGGTGTACTActggc                                                        >  1:521800/7‑150 (MQ=255)
                                                                                               gcgAGCCCGCGAGCACAGCCACGAGGGGGCCAGCGAGTTGCTGGCCGGCATTGCCCGGCTGGCCAATGCCGGTGCCAGCGAGCAGGCCCGTAGCGAATGCCAACGCTACCTGAGCCAGTCCCCGCCGTCGGCACCGGTGTACTActggct                                                       >  2:4182308/1‑150 (MQ=255)
                                                                                                 gAGCCCGCGAGCACAGCCACGAGGGGGCCAGCGAGTTGCTGGCCGGCATTGCCCGGCTGGCCAATGCCGGTGCCAGCGAGCAGGCCCGTAGCGAATGCCAACGCTACCTGAGCCAGTACCCGCCGTCGGCACAGGTGTACTActggctgg                                                     <  1:50990/150‑1 (MQ=255)
                                                                                                     ccGCGAGCACAGCCACGAGGGGGCCAGCGAGTTGCTGGCCGGCATTGCCCGGCTGGCCAATGCCGGTGCCAGCGAGCAGGCCCGTAGCGAATGCCAACGCTACCTGAGCCAGTACCCGCCGTCGGCACAGGTGTACTACTGGCTGGGttt                                                 >  2:4999612/1‑150 (MQ=255)
                                                                                                     ccGCGAGCACAGCCACGAGGGGGCCAGCGAGTTGCTGGCCGGCATTGCCCGGCTGGCCAATGCCGGTGCCAGCGAGCAGGCCCGTAGCGAATGCCAACGCTACCTGAGCCAGTACCCGCCGTCGGCACAGGTGTACTACTGGCTGGGttt                                                 <  2:2773666/150‑1 (MQ=255)
                                                                                                      cgcgAGCACAGCCACGAGGGGGCCAGCGAGTTGCTGGCCGGCATTGCCCGGCTGGCCAATGCCGGTGCCAGCGAGCAGGCCCGTAGCGAATGCCAACGCTACCTGAGCCAGTACCCGCCGTCGGCACAGGTGTACTACTGGCTGGGTTTg                                                <  2:2331120/150‑1 (MQ=255)
                                                                                                      cgcgAGCACAGCCACGAGGGGGCCAGCGAGTTGCTGGCCGGCATTGCCCGGCTGGCCAATGCCGGTGCCAGCGAGCAGGCCCGTAGCGAATGCCAACGCTACCTGAGCCAGTACCCGCCGTCGGCACAGGTGTACTACTGGCTGGGTTTg                                                >  1:402998/1‑150 (MQ=255)
                                                                                                       gcgAGCACAGCCACGAGGGGGCCAGCGAGTTGCTGGCCGGCATTGCCCGGCTGGCCAATGCCGGTGCCAGCGAGCAGGCCCGTAGCGAATGCCAACGCTACCTGAGCCAGTACCCGCCGTCGGCACAGGTGTACTACTGGCTGGGTTTGc                                               >  2:2135754/1‑150 (MQ=255)
                                                                                                        cgAGCACAGCCACGAGGGGGCCAGCGAGTTGCTGGCCGGCATTGCCCGGCTGGCCAATGCCGGTGCCAGCGAGCAGGCCCGTAGCGAATGCCAACGCTACCTGAGCCAGTACCCGCCGTCGGCACAGGTGTACTACTGGCTGGGTTTGCt                                              <  2:1666686/150‑1 (MQ=255)
                                                                                                         gAGCACAGCCACGAGGGGGCCAGCGAGTTGCTGGCCGGCATTGCCCGGCTGGCCAATGCCGGTGCCAGCGAGCAGGCCCGTAGCGAATGCCAACGCTACCTGAGCCAGTACCCGCCGTCGGCACAGGTGTACTACTGGCTGGGTTTGCTc                                             >  2:5394121/1‑150 (MQ=255)
                                                                                                          aGCACAGCCACGAGGGGGCCAGCGAGTTGCTGGCCGGCATTGCCCGGCTGGCCAATGCCGGTGCCAGCGAGCAGGCCCGTAGCGAATGCCAACGCTACCTGAGCCAGTACCCGCCGTCGGCACAGGTGTACTACTGGCTGGGTTTGCTCa                                            <  1:2907731/150‑1 (MQ=255)
                                                                                                           gCACAGCCACGAGGGGGCCAGCGAGTTGCTGGCCGGCATTGCCCGGCTGGCCAATGCCGGTGCCAGCGAGCAGGCCCGTAGCGAATGCCAACGCTACCTGAGCCAGTACCCGCCGTCGGCACAGGTGTACTACTGGCTGGGTTTGCTcgg                                           >  1:1185142/1‑148 (MQ=255)
                                                                                                            cacaGCCACGAGGGGGCCAGCGAGTTGCTGGCCGGCATTGCCCGGCTGGCCAATGCCGGTGCCAGCGAGCAGGCCCGTAGCGAATGCCAACGCTACCTGAGCCAGTACCCGCCGTCGGCACAGGTGTACTACTGGCTGGGTTTGCTCAGt                                          <  2:4959629/150‑1 (MQ=255)
                                                                                                                 ccACGAGGGGGCCAGCGAGTTGCTGGCCGGCATTGCCCGGCTGGCCAATGCCGGTGCCAGCGAGCAGGCCCGTAGCGAATGCCAACGCTACCTGAGCCAGTACCCGCCGTCGGCACAGGTGTACTACTGGCTGGGTTTGCTCAGTGATAc                                     >  1:1279982/1‑150 (MQ=255)
                                                                                                                 ccACGAGGGGGCCAGCGAGTTGCTGGCCGGCATTGCCCGGCTGGCCAATGCCGGTGCCAGCGAGCAGGCCCGTAGCGAATGCCAACGCTACCTGAGCCAGTACCCGCCGTCGGCACAGGTGTACTACTGGCTGGGTTTGCTCAGTGATAc                                     >  2:3970837/1‑150 (MQ=255)
                                                                                                                 ccACGAGGGGGCCAGCGAGTTGCTGGCCGGCATTGCCCGGCTGGCCAATGCCGGTGCCAGCGAGCAGGCCCGTAGCGAATGCCAACGCTACCTGAGCCAGTACCCGCCGTCGGCACAGGTGTACTACTGGCTGGGTTTGCTCAGTGATAc                                     >  1:4481100/1‑150 (MQ=255)
                                                                                                                          ggCCAGCGAGTTGCTGGCCGGCATTGCCCGGCTGGCCAATGCCGGTGCCAGCGAGCAGGCCCGTAGCGAATGCCAACGCTACCTGAGCCAGTACCCGCCGTCGGCACAGGTGTACTACTGGCTGGGTTTGCTCAGTGATACCGAAGGCGa                            >  2:4015678/1‑150 (MQ=255)
                                                                                                                           gCCAGCGAGTTGCTGGCCGGCATTGCCCGGCTGGCCAATGCCGGTGCCAGCGAGCAGGCCCGTAGCGAATGCCAACGCTACCTGAGCCAGTACCCGCCGTCGGCACAGGTGTACTACTGGCTGGGTTTGCTCAGTGATACCGAAGGCGAt                           >  1:784294/1‑150 (MQ=255)
                                                                                                                           gCCAGCGAGGTGCTGGCCGGCATTGCCCGGCTGGCCAATGCCGGTGCCAGCGAGCAGGCCCGTAGCGAATGCCAACGCTACCTGAGCCAGTACCCGCCGTCGGCACAGGTGTACTACTGGCTGGGTTTGCTCAGTGATACCGAAGGCGAt                           <  2:5617492/150‑1 (MQ=255)
                                                                                                                             cAGCGAGTTGCTGGCCGGCATTGCCCGGCTGGCCAATGCCGGTGCCAGCGAGCAGGCCCGTAGCGAATGCCAACGCTACCTGAGCCAGTACCCGCCGTCGGCACAGGTGTACTACTGGCTGGGTTTGCTCAGTGATACCGAAGGCGATgc                         >  2:1398789/1‑150 (MQ=255)
                                                                                                                             cAGCGAGTTGCTGGCCGGCATTGCCCGGCTGGCCAATGCCGGTGCCAGCGAGCAGGCCCGTAGCGAATGCCAACGCTACCTGAGCCAGTACCCGCCGTCGGCACAGGTGTACTACTGGCTGGGTTTGCTCAGTGATACCGAAGGCGATgc                         >  1:4580164/1‑150 (MQ=255)
                                                                                                                              aGCGAGTTGCTGGCCGGCATTGCCCGGCTGGCCAATGCCGGTGCCAGCGAGCAGGCCCGTAGCGAATGCCAACGCTACCTGAGCCAGTACCCGCCGTCGGAACAGGTGTACTACTGGCTGGGTTTGCTCAGTGATACCGAAGGCGATgcg                        >  1:5583976/1‑150 (MQ=255)
                                                                                                                               gCGAGTTGCTGGCCGGCATTGCCCGGCTGGCCAATGCCGGTGCCAGCGAGCAGGCCCGTAGCGAATGCCAACGCTACCTGAGCCAGTACCCGCCGTCGGCACAGGTGTACTACTGGCTGGGTTTGCTCAGTGATACCGAAGGCGATgcgc                       <  2:3195569/150‑1 (MQ=255)
                                                                                                                                 gAGTTGCTGGCCGGCATTGCCCGGCTGGCCAATGCCGGTGCCAGCGAGCAGGCCCGTAGCGAATGCCAACGCTACCTGAGCCAGTACCCGCCGTCGGCACAGGTGTACTACTGGCTGGGTTTGCTCAGTGATACCGAAGGCGATGCgcag                     <  1:4412145/150‑1 (MQ=255)
                                                                                                                                   gTTGCTGGCCGGCATTGCCCGGCTGGCCAATGCCGGTGCCAGCGAGCAGGCCCGTAGCGAATGCCAACGCTACCTGAGCCAGTACCCGCCGTCGGCACAGGTGTACTACTGGCTGGGTTTGCTCAGTGATACCGAAGGCGATGCgcagca                   <  1:4478387/150‑1 (MQ=255)
                                                                                                                                    ttGCTGGCCGGCATTGCCCGGCTGGCCAATGCCGGTGCCAGCGAGCAGGCCCGTAGCGAATGCCAACGCTACCTGAGCCAGTACCCGCCGTCGGCACAGGTGTACTACTGGCTGGGTTTGCTCAGTGATACCGAAGGCGATGCGcagcag                  >  2:3248549/1‑150 (MQ=255)
                                                                                                                                     tGCTGGCCGGCATTGCCCGGCTGGCCAATGCCGGTGCCAGCGAGCAGGCCCGTAGCGAATGCCAACGCTACCTGAGCCAGTACCCGCCGTCGGCACAGGTGTACTACTGGCTGGGTTTGCTCAGTGATACCGAAGGCGATGCGCAGCAgg                 <  2:4315627/150‑1 (MQ=255)
                                                                                                                                      tggcaccGGCATTGGCCAGCCGGGCAATGCCG‑‑GCCAGCAAGCAGGCCCGTAGCGAATGCCAACGCTACCTGAGCCAGTACCCGCCGTCGGCACAGGTGTACTACTGGCTGGGTTTGCTCAGTGATACCGAAGGCGATGCGCAGCAGgcgc              >  2:686827/6‑150 (MQ=38)
                                                                                                                                      gCTGGCCGGCATTGCCCGGCTGGCCAATGCCGGTGCCAGCGAGCAGGCCCGTAGCGAATGCCAACGCTACCTGAGCCAGTACCCGCCGTCGGCACAGGTGTACTACTGGCTGGGTTTGCTCAGTGATACCGAAGGCGATGCGCAGCAGgc                <  2:2516543/150‑1 (MQ=255)
                                                                                                                                        tGGCCGGCATTGCCCGGCTGGCCAATGCCGGTGCCAGCGAGCAGGCCCGTAGCGAATGCCAACGCTACCTGAGCCAGTACCCGCCGTCGGCACAGGTGTACTACTGGCTGGGTTTGCTCAGTGATACCGAAGGCGATGCGCAGCAGgcgc              >  2:4800685/1‑150 (MQ=255)
                                                                                                                                        tGGCCGGCATTGCCCGGCTGGCCAATGCCGGTGCCAGCGAGCAGGCCCGTAGCGAATGCCAACGCTACCTGAGCCAGTACCCGCCGTCGGCACAGGTGTACTACTGGCTGGGTTTGCTCAGTGATACCGAAGGCGATGCGCAGCAGgcgc              >  1:2124792/1‑150 (MQ=255)
                                                                                                                                        tGGCCGGCATTGCCCGGCTGGCCAATGCCGGTGCCAGCGAGCAGGCCCGTAGCGAATGCCAACGCTACCTGAGCCAGTACCCGCCGTCGGCACAGGTGTACTACTGGCTGGGTTTGCTCAGTGATACCGAAGGCGATGCGCAGCAGgcgc              >  1:684850/1‑150 (MQ=255)
                                                                                                                                            cGGCATTGGCCAGCCGGGCAATGCCGGTGCCAGCGAGCAGGCCCGTAGCGAATGCCAACGCTACCTGAGCCAGTACCCGCCGTCGGCACAGGTGTACTACTGGCTGGGTTTGCTCAGTGATACCGAAGGCGATGCGCAGCAGGCGCTGAg          >  1:272239/1‑150 (MQ=255)
                                                                                                                                             ggCATTGCCCGGCTGGCCAATGCCGGTGCCAGCGAGCAGGCCCGTAGCGAATGCCAACGCTACCTGAGCCAGTACCCGCCGTCGGCACAGGTGTACTACTGGCTGGGTTTGCTCAGTGATACCGAAGGCGATGCGCAGCAGGCGCTGAGc         >  2:4246280/1‑150 (MQ=255)
                                                                                                                                              gCATTGGCCAGCCGGGCAATGCCGGTGCCAGCGAGCAGGCCCGTAGCGAATGCCAAc                                                                                                     >  2:5467441/1‑57 (MQ=25)
                                                                                                                                               cATTGCCCGGCTGGCCAATGCCGGTGCCAGCGAGCAGGCCCGTAGCGAATGCCAACGCTACCTGAGCCAGTACCCGCCGTCGGCACAGGTGTACTACTGGCTGGGTTTGCTCAGTGATACCGAAGGCGATGCGCAGCAGGCGCTGAGCCa       <  1:5311326/150‑1 (MQ=255)
                                                                                                                                                    cccGGCTGGCCAATGCCGGTGCCAGCGAGCAGGCCCGTAGCGAATGCCAACGCTACCTGAGCCAGTACCCGCCGTCGGCACAGGTGTACTACTGGCTGGGTTTGCTCAGTGATACCGAAGGCGATGCGCAGCAGGCGCTGAGCCACTATc  >  2:5513024/1‑150 (MQ=255)
                                                                                                                                                    cccGGCTGGCCAATGCCGGTGCCAGCGAGCAGGCCCGTAGCGAATGCCAACGCTACCTGAGCCAGTACCCGCCATCGGCACAGGTGTACTACTGGCTGGGTTTGCTCAGTGATACCGAAGGCGATGCGCAGCAGGCGCTGAGCCACTATc  <  2:662155/150‑1 (MQ=255)
                                                                                                                                                    |                                                                                                                                                     
CGTGCGTTCACGACGTTGCCCGCACCGGTTTATCCGCAGCCCAGTGTGCCGTTGCCTCGTTCGCGGCGTGTGCTGCCGGTGGCGGCCAGGCCAGCGCGAGCCCGCGAGCACAGCCACGAGGGGGCCAGCGAGTTGCTGGCCGGCATTGCCCGGCTGGCCAATGCCGGTGCCAGCGAGCAGGCCCGTAGCGAATGCCAACGCTACCTGAGCCAGTACCCGCCGTCGGCACAGGTGTACTACTGGCTGGGTTTGCTCAGTGATACCGAAGGCGATGCGCAGCAGGCGCTGAGCCACTATC  >  NC_002947/1699863‑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: