Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 5,618,170 G→C 6.2% P178R (CCT→CGT)  wzy ← O‑antigen polymerase

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0029475,618,1700GC 6.2% 371.3 / 9.6 117P178R (CCT→CGT) wzyO‑antigen polymerase
Reads supporting (aligned to +/- strand):  ref base G (47/62);  new base C (5/3);  total (52/65)
Fisher's exact test for biased strand distribution p-value = 4.64e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

Only 100 of 119 total aligned reads displayed.

CCGGCTGTGGCGGTCACGGAACCACAACGGCGCCATCACCACGGTAATCACCAGTGCCAGAACCGCACCACGGCTCTGGCTGAGCATGGCGAAGGCCCCCAGGCAAGCCAGTGCCGCCAGCCACAACAACTGCAAGCCACGTCGCCGAGGTGGCTCGTAGAGCATCAAGAGCAAGGCCGAACCAATGACATAAGCCCCCAGGATGGGGTGGGATATCTCGCCGATACCCGCCAGCCTGAACAGCAAAGGCGCACGCTGCACCCCATAGAACTTGATGATGGAGACCAACGCAGCGAT  >  NC_002947/5618022‑5618318
                                                                                                                                                    |                                                                                                                                                      
ccGGCTGTGGCGGTCACGGAACCACAACGGCGCCATCACCACGGTAATCACCAGTGCCAGAACCGCACCACGGCTCTGGCTGAGCATGGCGAAGGCCCCCAGGCAAGCCAGTGCCGCCAGCCACAACAACTGCAAGCCACGTCGCCGAgg                                                                                                                                                       <  1:2014685/150‑1 (MQ=255)
ccGGCTGTGGCGGTCACGGAACCACAACGGCGCCATCACCACGGTAATCACCAGTGCCAGAACCGCACCACGGCTCTGGCTGAGCATGGCGAAGGCCCCCAGGCAAGCCAGTGCCGCCAGCCACAACAACTGCAAGCCACGTCGCCGAgg                                                                                                                                                       >  2:435027/1‑150 (MQ=255)
 cGGCTGTGGCGGTCACGGAACCACAACGGCGCCATCACCACGGTAATCACCAGTGCCAGAACCGCACCACGGCTCTGGCTGAGCATGGCGAAGGCCCCCAGGCAAGCCAGTGCCGCCAGCCACAACAACTGCAAGCCACGTCGCCGAGGt                                                                                                                                                      <  2:2212808/150‑1 (MQ=255)
 cGGCTGTGGCGGTCAAGGAACCGCTACAGCGCAATCACCACGGCACTCACAAGTGCCAGCACCTCATCACGGCTCTGGGTGAGCATGGCGAACGCCCCCAGGCCAGCCAGTGCCGCCAGCCACAACAACTGCAAGCCACGTCGCCGAGGt                                                                                                                                                      <  2:389971/150‑1 (MQ=255)
   gCTGTGGCGGTCACGGAACCACAACGGCGCCATCACCACGGTAATCACCAGTGCCAGAACCGCACCACGGCTCTGGCTGAGCATGGCGAAGGCCCCCAGGCAAGCCAGTGCCGCCAGCCACAACAACTGCAAGCCACGTCGCCGAGGTgg                                                                                                                                                    <  2:1197553/150‑1 (MQ=255)
    cTGTGGCGGTCACGGAACCACAACGGCGCCATCACCACGGTAATCACCAGTGCCAGAACCGCACCACGGCTCTGGCTGAGCATGGCGAAGGCCCCCAGGCAAGCCAGTGCCGCCAGCCACAACAACTGCAAGCCACGTCGCCGAGGTGGc                                                                                                                                                   <  1:220594/150‑1 (MQ=255)
        ggCGGTCACGGAACCACAACGGCGCCATCACCACGGTAATCACCAGTGCCAGAACCGCACCACGGCTCTGGCTGAGCATGGCGAAGGCCCCCAGGCAAGCCAGTGCCGCCAGCCACAACAACTGCAAGCCACGTCGCCGAGGTGGCTCGt                                                                                                                                               <  1:1997347/150‑1 (MQ=255)
             tCACGGAACCACAACGGCGCCATCACCACGGTAATCACCAGTGCCAGAACCGCACCACGGCTCTGGCTGAGCATGGCGAAGGCCCCCAGGCAAGCCAGTGCCGCCAGCCACAACAACTGCAAGCCACGTCGCCGAGGTGGCTCGTAGAGc                                                                                                                                         <  1:566435/150‑1 (MQ=255)
               aCGGAACCACAACGGCGCCATCACCACGGTAATCACCAGTGCCAGAACCGCACCACGGCTCTGGCTGAGCATGGCGAAGGCCCCCAGGCAAGCCAGTGCTCTCAGCCACAACAACTGCAAGCCACGTCGCCGAGGTGGCTCGTAGAGCAt                                                                                                                                       <  1:1042523/150‑1 (MQ=255)
               aCGGAACCACAACGGCGCCATCACCACGGTAATCACCAGTGCCAGAACCGCACCACGGCTCTGGCTGAGCATGGCGAAGGCCCCCAGGCAAGCCAGTGCCGCCAGCCACAACAACTGCAAGCCACGTCGCCGAGGTGGCTCGTAGAGCAt                                                                                                                                       <  1:258121/150‑1 (MQ=255)
               aCGGAACCACAACGGCGCCATCACCACGGTAATCACCAGTGCCAGAACCGCACCACGGCTCTGACTGAGCATGGCGAAGGCCCCCAGGCAAGCCAGTGCCGCCAGCCACAACAACTGCAAGCCACGTCGCCGGGGTGGATCGTagagaac                                                                                                                                       >  2:1097674/1‑147 (MQ=255)
                 ggAACCACAACGGCGCCATCACCACGGTAATCACCAGTGCCAGAACCGCACCACGGCTCTTGCTGAGCATGGCGAAGGCCCCCAGGCAAGCCAGTGCCGCCAGCCACAACAACTGCAAGCCACGTCGCCGAGGTGGCTCGTAGAGCATc                                                                                                                                      <  1:81740/149‑1 (MQ=255)
                 ggAACCACAACGGCGCCATCACCACGGTAATCACCAGTGCCAGAACCGCACCACGGCTCTTGCTGAGCATGGCGAAGGCCCCCAGGCAAGCCAGTGCCGCCAGCCACAACAACTGCAAGCCACGTCGCCGAGGTGGCTCGTAGAGCATc                                                                                                                                      >  2:81740/1‑149 (MQ=255)
                  gAACCACAACGGCCCCATCACCACGGTAATCACCAGTGCCAGAACCGCACCACGGCTCTGGCTGAGCATGGCGAAGGCCCCCAGGCAAGCCAGTGCCGCCAGCCACAACAACTGCAAGCCACGTCGCCGAGGTGGCTCGTAGAGCATCaa                                                                                                                                   <  1:1432042/150‑1 (MQ=255)
                    aCCACAACGGCGCCATCACCACGGTAATCACCAGTGCCAGAACCGCACCACGGCTCTGGCTGAGCATGGCGAAGGCCCCCAGGCAAGCCAGTGCCGCCAGCCACAACAACTGCAAGCCACGTCGCCGAGGTGGCTCGTAGAGCATCAaga                                                                                                                                 <  2:2251977/150‑1 (MQ=255)
                       acaACGGCGCCATCACCACGGTAATCACCAGTGCCAGAACCGCACCACGGCTCTGGCTGAGCATGGCGAAGGCCCCCAGGCAAGCCAGTGCCGCCAGCCACAACAACTGCAAGCCACGTCGCCGAGGTGGCTCGTAGAGCCTCAAGAGCa                                                                                                                              >  1:609122/1‑150 (MQ=255)
                           cGGCGCCATCACCACGGTAATCACCAGTGCCAGAACCGCACCACGGCTCTGGCTGAGCATGGCGAAGGCCCCCAGGCAAGCCAGTGCCGCCAGCCACAACAACTGCAAGCCACGTCGCCGAGGTGGCTCGTAGAGCATCAAGAGCAAGGc                                                                                                                          <  1:435027/150‑1 (MQ=255)
                            ggCGCCATCACCACGGTAATCACCAGTGCCAGAACCGCACCACGGCTCTGGCTGAGCATGGCGAAGGCCCCCAGGCAAGCCAGTGCCGCCAGCCACAACAACTGCAAGCCACGTCGCCGAGGTGGCTCGTAGAGCATCAAGAGCAAGGcc                                                                                                                         <  2:479781/150‑1 (MQ=255)
                             tctcgaccaccacTGTAATCACCAGTGCCAGAACCGCACCACGGCTCTGGCTGAGCATGGCGAAGGCCTCAAGGAAAGCCAGTGCCGCCAGCCACAACAACTGCAAGCCACGTCGCCGAGGTGGCTCGTAGAGCATCAAGAGCAAGGCCg                                                                                                                        <  1:158775/143‑1 (MQ=255)
                               gcCATCACCACGGTAATCACCAGTGCCAGAACCGCACCACGGCTCTGGCTGAGCATGGCGAAGGCCCCCAGGCAAGCCAGTGCCGCCAGCCACAACAACTGCAAGCCACGTCGCCGAGGTGGCTCGTAGAGCATCAAGAGCAAGGCCGaa                                                                                                                      <  2:304300/150‑1 (MQ=255)
                               gcCATCACCACGGTAATCACCAGTGCCAGAACCGCACCACGGCTCTGGCTGAGCATGGCGAAGGCCCCCAGGCAAGCCAGTGCCGCCAGCCACAACAACTGCAAGCCACGTCGCCGAGGTGGCTCGTAGAGCATCAAGAGCAAGGCCGaa                                                                                                                      >  1:690097/1‑150 (MQ=255)
                                 cATCACCACGGTAATCACCAGTGCCAGAACCGCACCACGGCTCTGGCTGAGCATGGCGAAGGCCCCCAGGCAAGCCAGTGCCGCCAGCCACAACAACTGCAAGCCACGTCGCCGAGGTGGCTCGTAGAGCATCAAGAGCAAGGCCGAAcc                                                                                                                    >  2:1207695/1‑150 (MQ=255)
                                   tCACCACGGTAATCACCAGTGCCAGAACCGCACCACGGCTCTGGCTGAGCATGGCGAAGGCCCCCAGGCAAGCCAGTGCCGCCAGCCACAACAACTGCAAGCCACGTCGCCAAGGTGGCTCGTAGAGCATCAAGAGCAAGGCCGAACCaa                                                                                                                  >  1:701334/1‑150 (MQ=255)
                                        acGGTAATCACCAGTGCCAGAACCGCACCACGGCTCTGGCTGAGCATGGCGAAGGCCCCCAGGCAAGCCAGTGCCGCCAGCCACAACAACTGCAAGCCACGTCGCCGAGGTGGCTCGTag                                                                                                                                            <  2:2211842/120‑1 (MQ=255)
                                        acGGTAATCACCAGTGCCAGAACCGCACCACGGCTCTGGCTGAGCATGGCGAAGGCCCCCAGGCAAGCCAGTGCCGCCAGCCACAACAACTGCAAGCCACGTCGCCGAGGTGGCTCGTag                                                                                                                                            >  1:2211842/1‑120 (MQ=255)
                                        acGGTAATCACCAGTGCCAGAACCGCACCACGGCTCTGGCTGAGCATGGCGAAGGCCCCCAGGCAAGCCAGTGCCGCCAGCCACAACAACTGCAAGCCACGTCGCCGAGGTGGCTCGTAGAGCATCAAGAGCAAGGCCGAACCAATGACa                                                                                                             >  1:1849840/1‑150 (MQ=255)
                                        acGGTAATCACCAGTGCCAGAACCGCACCACGGCTCTGGCTGAGCATGGCGAAGGCCCCCAGGCAAGCCAGTGCCGCCAGCCACAACAACTGCAAGCCACGTCGCCGAGGTGGCTCGTAGAGCATCAAGAGCAAGGCCGAACCAATGACa                                                                                                             >  1:2115640/1‑150 (MQ=255)
                                           gTAATCACCAGTGCCAGAACCGCACCACGGCTCTGGCTGAGCATGGCGAAGGCCCCCAGGCAAGCCAGTGCCGCCAGCCACAACAACTGCAAGCCACGTCGCCGAGGTGGCTCGTAGAGCATCAAGAGCAAGGCCGAACCAATGACATaa                                                                                                          <  2:1849840/150‑1 (MQ=255)
                                           gTAATCACCAGTGCCAGAACCGCACCACGGCTCTGGCTGAGCATGGCGAAGGCCCCCAGGCAAGCCAGTGCCGCCAGCCACAACAACTGCAAGCCACGTCGCCGAGGTGGCTCGTAGAGCATCAAGAGCAAGGCCGAACCAATGACATaa                                                                                                          <  2:690097/150‑1 (MQ=255)
                                           ccaaTCATCGGTGCCAGAACCGCACTACGGCTCTGGCTGAGCAAGGCGAAGGCCCCCGGGCAAGCCAGTGCCCCCAGCCACAACAACTGCAATCCACGTCGCCGAGGTGGCTCGTAGAGCATCAAGAGCAAGGCCGAACCAATGACATaa                                                                                                          <  2:468911/148‑1 (MQ=255)
                                              aTCACCAGTGCCAGAACCGCACCACGGCTCTGGCTGAGCATGGCGAAGGCCCCCAGGCAAGCCAGTGCCGCCAGCCACAACAACTGCAAGCCACGTCGCCGAGGTGGCTCGTa                                                                                                                                              <  2:640420/113‑1 (MQ=255)
                                              aTCACCAGTGCCAGAACCGCACCACGGCTCTGGCTGAGCATGGCGAAGGCCCCCAGGCAAGCCAGTGCCGCCAGCCACAACAACTGCAAGCCACGTCGCCGAGGTGGCTCGTa                                                                                                                                              >  1:640420/1‑113 (MQ=255)
                                               tCACCAGTGCCAGAACCGCACCACGGCTCTGGCTGAGCATGGCGAAGGCCCCCAGGCAAGCCAGTGCCGCCAGCCACAACAACTGCAAGCCACGTCGCCGAGGTGGCTCGTCGAGCATCAAGAGCAAGGCCGAACCAATGACATAAGccc                                                                                                      <  2:1391063/150‑1 (MQ=255)
                                                 aCCAGTGCCAGAACCGCACCACGGCTCTGGCTGAGCATGGCGAAGGCCCCCAGGCAAGCCAGTGCCGCCAGCCACAACAACTGCAAGCCACCTCGGCGACGTGGCttg                                                                                                                                                <  1:2202570/108‑3 (MQ=255)
                                                 aCCAGTGCCAGAACCGCACCACGGCTCTGGCTGAGCATGGCGAAGGCCCCCAGGCAAGCCAGTGCCGCCAGCCACAACAACTGCAAGCCACCTCGGCGACGTGGCttg                                                                                                                                                >  2:2202570/1‑106 (MQ=255)
                                                  ccAGTGCCAGAACCGCACCACGGCTCTGGCTGAGCATGGCGAAGGCCCCCAGGCAAGCCAGTGCCGCCAGCCACAACAACTGCAAGCCACGTCGCCGAGGTGGCTCg                                                                                                                                                >  2:1824057/1‑107 (MQ=255)
                                                  ccAGTGCCAGAACCGCACCACGGCTCTGGCTGAGCATGGCGAAGGCCCCCAGGCAAGCCAGTGCCGCCAGCCACAACAACTGCAAGCCACGTCGCCGAGGTGGCTCg                                                                                                                                                <  1:1824057/107‑1 (MQ=255)
                                                  ccAGTGCCAGAACCGCACCACGGCTCTGGCTGAGCATGGCGAAGGCCCCCAGGCAAGCCAGTGCCGCCAGCCACAACAACTGCAAGCCACGTCGCCGAGGTGGCTCGTAGAGCATCAAGAGCAAGGCCGAACCAATGACATAAGCCCCCa                                                                                                   >  2:1127143/1‑150 (MQ=255)
                                                    aGTGCCAGTACCGACCCACGGCGCTTGCTTAGCATGGCGAAGGCCCCCAGGCAAGCCAGTGCCGCCAGCCACAACAACTGCAAGCCACGTCGCCGAGGTGGCTCGTAGAGCATCAAGAGCAAGGCCGAACCAATGACATAAGCCCCCAgg                                                                                                 <  1:1698283/150‑1 (MQ=255)
                                                    aGTGCCAGAACCGCACCACGGCTCTGGCTGAGCATGGCGAAGGCCCCCAGGCAAGCCAGTGCCGCCAGCCACAACAACTGCAAGCCACGTCGCCGAGGTGGCTCGTAGAGCATCAAGAGCAAGGCCGAACCAATGACATAAGCCCCCAgg                                                                                                 >  2:42092/1‑150 (MQ=255)
                                                     gTGCCAGAACCGCACCACGGCTCTGGCTGAGCATGGCGAAGGCCCCCAGGCAAGCCAGTGCCGCCAGCCACAACAACTGCAAGCCACGTCGCCGAGGTGGCTCGTAGAGCATCAAGAGCAAGGCCGAACCAATGACATAAGCCCCCAGGa                                                                                                >  1:2404757/1‑150 (MQ=255)
                                                        ccAGAACCGCACCACGGCTCTGGCTGAGCATGGCGAAGGCCCCCAGGCAAGCCAGTGCCGCCAGCCACAACAACTGCAAGCCACGTCGCCGAGGTGGCTCGTAGAGCATCAAGAGCAAGGCCGAACCAATGACATAAGCCCCCAGGATgg                                                                                             <  2:2085526/150‑1 (MQ=255)
                                                            aaCCGCACCACGGCTCTGGCTGAGCATGGCGAAGGCCCCCAGGCAAGCCAGTGCCGCCAGCCACAACAACTGCAAGCCACGTCGCCGAGGTGGCTCGTAGAGCATCAAGAGCAAGGCCGAACCAATGACATAAGCCCCCAGGATGGGGTg                                                                                         <  2:1547599/150‑1 (MQ=255)
                                                             aCCGCACCACGGCTCTGGCTGAGCATGGCGAAGGCCCCCAGGCAAGCCAGTGCCGCCAGCCACAACAACTGCAAGCCACGTCGCCGAGGTGGCTCGTAGAGCATCAAGAGCAAGGCCGAACCAATGACATAAGCCCCCAGGATGGGGTgg                                                                                        <  1:1127143/150‑1 (MQ=255)
                                                             aCCGCACCACGGCTCTGGCTGAGCATGGCGAAGGCCCCCAGGCAAGCCAGTGCCGCCAGCCACAACAACTGCAAGCCACGTCGCCGAGGTGGCTCGTAGAGCATCAAGAGCAAGGCCGAACCAATGACATAAGCCCCCAGGATGGGGTgg                                                                                        <  2:1264255/150‑1 (MQ=255)
                                                             aCCGCACCACGGCTCTGGCTGAGCATGGCGAAGGCCCCCAGGCAAGCCAGTGCCGCCAGCCACAACAACTGCAAGCCACGTCGCCGACGTGGCTCGTAGAGCACCAAGAGCAAGGCCGAACCAATGCCCTAAGCCCCCAGCATGGGGTgg                                                                                        <  2:788658/150‑1 (MQ=255)
                                                             aCCGCACCACGGCTCTGGCTGAGCATGGCGAAGGCCCCCAGCCAAGCCAGTGCCGCCAGCCACAACAACTGCAAGCCACGTCGCCGAGGTGGCTCGCAGAGCATCAAGAGCAAGGCCGAACCAATGACATAAGCCCCCAGGATGGGGTgg                                                                                        <  1:1774815/150‑1 (MQ=255)
                                                                gCACCACGGCTCTGGCTGAGCATGGCGAAGGCCCCCAGGCAAGCCAGTGCCGCCAGCCACAACAACTGCAAGCCACGTCGCCGAGGTGGCTCGTAGAGCATCAAGAGCAAGGCCGAACCAATGACATAAGCCCCCAGGATGGGGTGGGat                                                                                     <  1:2044672/150‑1 (MQ=255)
                                                                gCACCACGGCTCTGGCTGAGCATGGCGAAGGCCCCCAGGCAAGCCAGTGCCGCCAGCCACAACAACTGCAAGCCACGTCGCCGAGGTGGCTCGTAGAGCATCAAGAGCAAGGCCGAACCAATGACATAAGCCCCCAGGATGGGGTGGGat                                                                                     >  1:1046541/1‑150 (MQ=255)
                                                                       ggCTCTGGCTGAGCATGGCGAAGGCCCCCAGGCAAGCCAGTGCCGCCAGCCACAACAACTGCAAGCCACGTCGCCGAGGTGGCTCGTAGAGCATCAAGAGCAAGGCCGAACCAATGACATAAGCCCCCAGGATGGGGTGGGATATCTCGc                                                                              <  2:2404757/150‑1 (MQ=255)
                                                                       ggCTCTGGCTGAGCATGGCGAAGGCCCCCAGGCAAGCCAGTGCCGCCAGCCACAACAACTGCAAGCCACGTCGCCGAGGTGGCTCGTAGAGCATCAAGAGCAAGGCCGAACCAATGACATAAGCCCCCAGGATGGGGTGGGATATCTCGc                                                                              <  1:1387186/150‑1 (MQ=255)
                                                                          tctGGGCGAGCACGGCGAAGGCCACCAGGGAAGCCAGTGCCGCCCGCCACAACACCTGCAAGCCACGTCGCCGAGGTGGCTCGTAGAGCATCAAGAGCAAGGCCGAACCAATGACATAGGACCCCAGGATGGGGTGGGATATCTCGCCGa                                                                           <  2:867844/150‑1 (MQ=255)
                                                                          tctGGCTGAGCATGGCGAAGGCCCCCAGGCAAGCCAGTGCCGCCAGCCACAACAACCGCAAGCCACGTCGCCGAGGTGGCTCGTAGAGCATCAAGAGCAAGGCCGAACCAATGACATAAGCCCCCAGGATGGGGTGGGATATCTCGCCGa                                                                           <  1:1948428/150‑1 (MQ=255)
                                                                              gCTGAGCATGGCGAAGGCCCCCAGGCAAGCCAGTGCCGCCAGCCACAACAACTGCAAGCCACGTCGCCGAGGTGGCTCGTAGAGCATCAAGAGCAAGGCCGAACCAATGACATAAGCCCCCAGGATGGGGTGGGATATCTCGCCGATAcc                                                                       >  2:954071/1‑150 (MQ=255)
                                                                                   gCATGGCGAAGGCCCCCAGGCAAGCCAGTGCCGCCAGCCACAACAACTGCAAGCCACGTCGCCGAGGTGGCTCGTAGGTGCAACAGGGCAAGGCCGAACCAATGACATAAGCCCCCAGGATGGGGAGGGATATCTTGCCGATACCCGCCa                                                                  >  1:150195/1‑150 (MQ=255)
                                                                                   gCATGGCGAAGGCCCCCAGGCAAGCCAGTGCCGCCAGCCACAACAACTGCAAGCCACGTCGCCGAGGTGGCTCGTAGAGCATCAAGAGCAAGGCCGAACCAATGACATAAGCCCCCAGGATGGGGTGGGATATCTCGCCGATACCCGCCa                                                                  >  2:841343/1‑150 (MQ=255)
                                                                                       ggCGAAGGCCCCCAGGCAAGCCAGTGCCGCCAGCCACAACAACTGCAAGCCACGTCGCCGAGGTGGCTCGTAGAGCATCAAGAGCAAGGCCGAAc                                                                                                                     >  2:455398/1‑95 (MQ=255)
                                                                                       ggCGAAGGCCCCCAGGCAAGCCAGTGCCGCCAGCCACAACAACTGCAAGCCACGTCGCCGAGGTGGCTCGTAGAGCATCAAGAGCAAGGCCGAAc                                                                                                                     <  1:455398/95‑1 (MQ=255)
                                                                                            aGGCCCCCAGGCAAGCCAGTGCCGCCAGCCACAACAACTGCAAGCCACGTCGCCGAGGTGGCTCGTAGAGCATCAAGAGCAAGGCCGAACCAATGACATAAGCCCCCAGGATGGGGTGGGATAtc                                                                                  >  1:2346772/1‑125 (MQ=255)
                                                                                            aGGCCCCCAGGCAAGCCAGTGCCGCCAGCCACAACAACTGCAAGCCACGTCGCCGAGGTGGCTCGTAGAGCATCAAGAGCAAGGCCGAACCAATGACATAAGCCCCCAGGATGGGGTGGGATAtc                                                                                  <  2:2346772/125‑1 (MQ=255)
                                                                                             ggCCCCCAGGCAAGCCAGTGCCGCCAGCCACAACAACTGCAAGCCACGTCGCCGAGGTGGCTCGTAGAGCATCAAGAGCAAGGCCGAACCAATGACATAAGCCCCCAGGATGGGGTGGGATATCTCGCCGATACCCGCCAGCCTGAACAg                                                        >  1:486768/1‑150 (MQ=255)
                                                                                               cccccAGGCAAGCCAGTGCCGCCAGCCACAACAACTGCAAGCCACGTCGCCGAGGTGGCTCGTAGAGCATCAAGAGCAAgg                                                                                                                           >  1:2391846/1‑81 (MQ=255)
                                                                                                ccccAGGCAAGCCAGTGCCGCCAGCCACAACAACTGCAAGCCACGTCGCCGAGGTGGCTCGTAGAGCATCAAGAGCAAGGCCGAACCAATGACATAAGCCCCCAGGATGGGGTGGGATATCTCGCCGATACCCGCCAGc                                                                <  1:290830/139‑1 (MQ=255)
                                                                                                ccccAGGCAAGCCAGTGCCGCCAGCCACAACAACTGCAAGCCACGTCGCCGAGGTGGCTCGTAGAGCATCAAGAGCAAGGCCGAACCAATGACATAAGCCCCCAGGATGGGGTGGGATATCTCGCCGATACCCGCCAGc                                                                >  2:290830/1‑139 (MQ=255)
                                                                                                ccccAGGCAAGCCAGTGCCGCCAGCCACAACAACTGCAAGCCACGTCGCCGAGGTGGCTCGTAGAGCATCAAGAGCAAGGCCGAACCAATGACATAAGCCCCCAGGATAGGTTGGGATATCTCGCCGATACCCGCCAGCCTGAACAGCaa                                                     >  1:1834220/1‑150 (MQ=255)
                                                                                                ccccAGGCAAGCCAGTGCCGCCAGCCACAACAACTGCAAGCCACCTCGGCGACGTGGCt                                                                                                                                                  <  1:663944/59‑1 (MQ=37)
                                                                                                ccccAGGCAAGCCAGTGCCGCCAGCCACAACAACTGCAAGCCACCTCGGCGACGTGGCt                                                                                                                                                  >  2:663944/1‑59 (MQ=37)
                                                                                                 cccAGGCAAGCCAGTGCCGCCAGCCACAACAACTGCAAGCCACGTCGCCGAGGTGGGTc                                                                                                                                                 <  1:2010974/59‑1 (MQ=255)
                                                                                                 cccAGGCAAGCCAGTGCCGCCAGCCACAACAACTGCAAGCCACGTCGCCGAGGTGGCTc                                                                                                                                                 >  2:2010974/1‑59 (MQ=255)
                                                                                                 cccAGGCAAGCCAGTGCCGCCAGCCACAACAACTGCAAGCCACGTCGCCGAGGTGGCTCGTAGAGCATCAAGCGCTGGACCGAACCAATGACATAAGCCCCCAGGATGGGGTGGGATATCTCGCCGATACCCGCCAGCCTGAACAGCaaa                                                    >  1:2099070/1‑150 (MQ=255)
                                                                                                 cccAGGCAAGCCAGTGCCGCCAGCCACAACAACTGCAAGCCACGTCGCCGAGGTGGCTCGTAGAGCATCAAGAGCAAGGCCGAACCAATGACATAAGccc                                                                                                      >  2:360744/1‑100 (MQ=255)
                                                                                                  ccAGGCAAGCCAGTGCCGCCAGCCACAACAACTGCAAGCCACGTCGCCGAGGTGGCTCGTAGAGCATCAAGAGCAAGGCCGAACCAATGACATAAGCCCCCAGGATGGGGTGGGATATCTCGCCGATACCCGCCAGCCTGAACAGCAAAg                                                   >  1:1263276/1‑150 (MQ=255)
                                                                                                    aGGCAAGCCAGTGGAGCCAGCCACAACAACTGCAAGGCACGTCCCCGAGGGGGCTCTTAGAGCATCAAGAGCAAGGCCGAACCAACGATATAAGCTCCCAGGATGGGGTGGGATATCTCGCCTATACCCGCCAGCCTGAACAGCAAAGgc                                                 <  1:1097674/150‑1 (MQ=255)
                                                                                                    aGGCAAGCCAGTGCCGCCAGCCACAACAACTGCAAGCCACGTCGCCGAGGTGGCTCGTAGAGCATCAAGAGCAAGGCCGAACCAATGACATAAGCCCCCAGGATGGGGTGGGATATCTCGCCGATACCCGCCAGCCTGAACAGCAAAGgc                                                 <  2:1072713/150‑1 (MQ=255)
                                                                                                    aGGCAAGCCAGTGCCGCCAGCCACAACAACTGCAAGCCACGTCGCCGAGGTGGCTCGTAGAGCATCAAGAGCAAGGCCGAACCAATGACATAAGCCCCCAGGATGGGGTGGGATATCTCGCCGATACCCGCCAGCCTGAACAGCAAAGgc                                                 <  1:1207695/150‑1 (MQ=255)
                                                                                                      gCAAGCCAGTGCCGCCAGCCACAACAACTGCAAGCCACGTCGCCGAGGTGGCTCGTAGAGCATCAAGAGCAAGGCCGAACCAATGACATAAGCCCCCAGGATGGGGTGGGATATCTCGCCGATACCCGCCAGCCTGAACAGCAAAGgcgc                                               >  1:1159007/1‑150 (MQ=255)
                                                                                                        aaGCCAGTGCCGCCAGCCACAACAACTGCAAGCCACGTCGCCGAGGTGGCTCGTAGAGCATCAAGAGCAAGGCCGAACCAATGACATAAGCCCCCAGGATGGGGTGGGATATCTCGCCGATACCCGCCAGCCTGAACAGCAAAGGCGCAc                                             >  2:1338177/1‑150 (MQ=255)
                                                                                                         aGCCAGTGCCGCCAGCCACAACAACTGCAAGCCACGTCGCCGAGGTGGCTCGTAGAGCATCAAGAGCAAGGCCGAACCAATGACATAAGCCCCCAGGATGGGGTGGGATATCTCGCCGATACCCGCCAGCCTGAACAGCAAAGGCGCACg                                            >  1:465045/1‑150 (MQ=255)
                                                                                                                   gccagccaCAACAACTGCAAGCCACGTCGCCGAGGTGGCTCGTAGAGCATCAAGAGCAAGGCCGAACCAATGACATAAGCCCCCAGGATGGGGTGGGATATCTCGCCGATACCCGCCAGCCTGAACAGCAAAGGCGCACGCTGCACCCCa                                  >  1:451011/1‑150 (MQ=255)
                                                                                                                    ccagccaCAACAACTGCAAGCCACGTCGCCGAGGTGGCTCGTAGAGCATCAAGAGCAAGGCCGAACCAATGACATAAGCCCCCAGGATGGGGTGGGATATCTCGCCGATACCCGCCAGCCTGAACAGCAAAGGCGCACGCTGCACCCCAt                                 <  2:1263276/150‑1 (MQ=255)
                                                                                                                    ccagccaCAACAACTGCAAGCCACGTAGCCGACATGGCTCGAACAAGTCACCTAGGAAGGCCGAACCAAGTGCATAAGCCCCCAGGATGGGGTGGGATATCTCGCCGATACCCGCCAGCCTGAACAGCAAAGGCGCACGCTGCACCCCAt                                 >  1:1622294/1‑150 (MQ=255)
                                                                                                                     cagccaCAACAACTGCAAGCCACGTCGCCGAGGTGGCTCGTAGAGCATCAAGAGCAAGGCCGAACCAATGACATAAGCCCCCAGGATGGGGTGGGATATCTCGCCGATACCCGCCAGCCTGAACAGCAAAGGCGCACGCTGCACCCCATa                                <  1:841343/150‑1 (MQ=255)
                                                                                                                        ccaCAACAACTGCAAGCCACGTCGCCGAGGTGGCTCGGAGAGCATCAAGAGCAAGGCCGAACCAATGACATAAGCCCCCAGGATGGGGTGGGATATCTCGCCGATACCCGCCAGCCTGAACAGCAAAGGCGCACGCTGCACCCCATAGaa                             <  2:1046541/150‑1 (MQ=255)
                                                                                                                         cacaACAACTGCAAGCCACGTCGCCGAGGTGGCTCGTAGAGCATCAAGAGCAAGGCCGAACCAATGACATAAGCCCCCAGGATGGGGTGGGATATCTCGCCGATACCCGCCAGCCTGAACAGCAAAGGCGCACGCTGCACCCCATAGAAc                            >  2:360121/1‑150 (MQ=255)
                                                                                                                            aacaacTGCAAGCCACGTCGCCGAGGTGGCTCGTAGAGCATCAAGAGCAAGGCCGAACCAATGACATAAGCCCCCAGGATGGGGTGGGATATCTCGCCGATACCCGCCAGCCTGAACAGCAAAGGCGCACGCTGCACCCCATAGAACTtg                         >  2:2006301/1‑150 (MQ=255)
                                                                                                                            aacaacTGCAAGCCACGTCGCCGAGGTGGCTCGTAGAGCATCAAGAGCAAGGCCGAACCAATGACATAAGCCCCCAGGATGGGGTGGGATATCTCGCCGATACCCGCCAGCCTGAACAGCAAAGGCGCACGCTGCACCCCATAGAACTtg                         >  1:1379483/1‑150 (MQ=255)
                                                                                                                             acaacTGCAAGCCACGTCGCCGAGGTGGCTCGTAGAGCATCAAGAGCAAGGCCGAACCAATGACATAAGCCCCCAGGATGGGGTGGGATATCTCGCCGATACCCGCCAGCCTGAACAGCAAAGGCGCACGCTGCACCCCATAGAACTtga                        <  1:954071/150‑1 (MQ=255)
                                                                                                                             acaacTGCAAACCACGTCGCCGAGGTGGCTCGTAGAGCATCAAGAGCAAGGCCGAACCAATGACATAAGCCCCCAGGATGGGGTGGGATATCTCGCCGATACCCGCCAGCCTGAACAGCAAAGGCGCACGCTGCACCCCATAGAACTtga                        <  2:1622294/150‑1 (MQ=255)
                                                                                                                              caacTGCAAGCCACGTCGCCGAGGTGGCTCGTAGAGCATCAAGAGCAAGGCCGAACCAATGACATAAGCCCCCAGGATGGGGTgg                                                                                        <  2:596635/85‑1 (MQ=255)
                                                                                                                                  tacgaGCCACCTCGGCGACGTGGCTCGTAGAGCATCAAGAGCAAGGCCGAACCAATGACATAAGCCCCCAGGATGGGGTGGGATATCTCGCCGATACCCGCCAGCCTGAACAGCAAAGGCGCACGCTGCACCCCATAGAACTTGatgatg                   >  2:1906281/5‑150 (MQ=255)
                                                                                                                                  tGCAAGCCACGTCGCCGAGGTGGCTCGTAGAGCATCAAGAGCAAGGCCGAACCAATGACATAAGCCCCCAGGATGGGGTgg                                                                                        <  2:1435802/81‑1 (MQ=255)
                                                                                                                                  tGCAAGCCACGTCGCCGAGGTGGCTCGTAGAGCATCAAGAGCAAGGCCGAACCAATGACATAAGCCCCCAGGATGGGGTgg                                                                                        >  1:1435802/1‑81 (MQ=255)
                                                                                                                                  tGCAAGCCACGTCGCCGAGGTGGCTCGTAGAGCATCAAGAGCAAGGCCGAACCAATGACATAAGCCCCCAGGATGGGGTGGGATATCTCGCCGATACCCGCCAGCCTGAACAGc                                                       <  2:1017197/114‑1 (MQ=255)
                                                                                                                                  tGCAAGCCACGTCGCCGAGGTGGCTCGTAGAGCATCAAGAGCAAGGCCGAACCAATGACATAAGCCCCCAGGATGGGGTGGGATATCTCGCCGATACCCGCCAGCCTGAACAGc                                                       >  1:1017197/1‑114 (MQ=255)
                                                                                                                                  tGCAAGCCACGTCGCCGAGGTGGCTCGTAGAGCATCAAGAGCAAGGCCGAACCAATGACATAAGCCCCCAGGATGGGGTGGGATATCTCGCCGATACCCGCCAGCCTGAACAGCAAAGGCGCACGCTGCACCCCATAGAACTTgatgat                    <  1:790753/149‑1 (MQ=255)
                                                                                                                                     gaGCCACCTCGGCGACGTGGCTCGTAGAGCATCAAGAGCAAGGCCGAACCAATGACATAAGCCCCCAGGATGGGGTGGGATATCTCGCCGATACCCGCCAGCCTGAAACGGGCGCTGGCACGCTGCACCCCATAGAACTTGATGATGgag                >  1:126203/2‑150 (MQ=255)
                                                                                                                                     aaGCCACGTCGCCGAGGTGGCTCGTAGAGCATCAAGAGCAAGGCCGAACCAATGACATAAGCCCCCAGGATGGGGTGGGATATCTCGCCGATACCCGCCAGCCTGAACAGCAAAgg                                                  >  1:912520/1‑116 (MQ=255)
                                                                                                                                     aaGCCACGTCGCCGAGGTGGCTCGTAGAGCATCAAGAGCAAGGCCGAACCAATGACATAAGCCCCCAGGATGGGGTGGGATATCTCGCCGATACCCGCCAGCCTGAACAGCAAAgg                                                  <  2:912520/116‑1 (MQ=255)
                                                                                                                                                 cGAGGTGGCTCGTAGAGCATCAAGAGCAAGGCCGAACCAATGACATAAGCCCCCAGGATGGGGTGGGATATCTCGCCGATACCCGCCAGCCTGAACAGCAAAAGCGCACGCTGCACCCCATAGAACTTGATGATGGAGACCAACGCAGCg    >  1:1774953/1‑150 (MQ=255)
                                                                                                                                                   aGGTGGCTCGTAGAGCATCAAGAGCAAGGCAGAACCAATGACATAAGCCCCCAGGATGGGGTGGGATATCTCGCCGATACCCGCCAGCCTGAACAGCAAAGGCGCACGCTGCACCACATAGAACTTGATGATGGAGACCAACGCAGCGAt  <  2:2099070/150‑1 (MQ=255)
                                                                                                                                                    |                                                                                                                                                      
CCGGCTGTGGCGGTCACGGAACCACAACGGCGCCATCACCACGGTAATCACCAGTGCCAGAACCGCACCACGGCTCTGGCTGAGCATGGCGAAGGCCCCCAGGCAAGCCAGTGCCGCCAGCCACAACAACTGCAAGCCACGTCGCCGAGGTGGCTCGTAGAGCATCAAGAGCAAGGCCGAACCAATGACATAAGCCCCCAGGATGGGGTGGGATATCTCGCCGATACCCGCCAGCCTGAACAGCAAAGGCGCACGCTGCACCCCATAGAACTTGATGATGGAGACCAACGCAGCGAT  >  NC_002947/5618022‑5618318

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