Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 5,618,162 G→C 5.0% R181G (CGT→GGT)  wzy ← O‑antigen polymerase

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0029475,618,1620GC 5.0% 383.5 / 8.7 119R181G (CGT→GGT) wzyO‑antigen polymerase
Reads supporting (aligned to +/- strand):  ref base G (50/63);  new base C (4/2);  total (54/65)
Fisher's exact test for biased strand distribution p-value = 4.09e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 7.73e-01

Only 100 of 121 total aligned reads displayed.

AGAGAACACCCGGCTGTGGCGGTCACGGAACCACAACGGCGCCATCACCACGGTAATCACCAGTGCCAGAACCGCACCACGGCTCTGGCTGAGCATGGCGAAGGCCCCCAGGCAAGCCAGTGCCGCCAGCCACAACAACTGCAAGCCACGTCGCCGAGGTGGCTCGTAGAGCATCAAGAGCAAGGCCGAACCAATGACATAAGCCCCCAGGATGGGGTGGGATATCTCGCCGATACCCGCCAGCCTGAACAGCAAAGGCGCACGCTGCACCCCATAGAACTTGATGATGGAG  >  NC_002947/5618013‑5618304
                                                                                                                                                     |                                                                                                                                               
agagAACACCCGGCTGTGGCGGTCACGGAACCACAACGGCGCCATCACCACGGTAATCACCAGTGCCAGAACCGCACCACGGCTCTGGCTGAGCATGGCGAAGGCCCCCAGGCAAGCCAGTGCCGCCAGCCACAACAACTGCAAGCCACg                                                                                                                                                 >  2:566435/1‑150 (MQ=255)
agagAACACCCGGCTGTGGCGGTCACGGAACCACAACGGCGCCATCACCACGGTAATCACCAGTGCCAGAACCGCACCACGGCTCTGGCTGAGCATGGCGAAGGCCCCCAGGCAAGCCAGTGCCGCCAGCCACAACAACTGCAAGCCACg                                                                                                                                                 >  1:468911/1‑150 (MQ=255)
  agaACACCCGGCTGTGGCGGTCACGGAACCACAACGGCGCCATCACCACGGTAATCACCAGTGCCAGAACCGCACCACGGCTCTGGCTGAGCATGGCGAAGGCCCCCAGGCAAGCCAGTGCCGCCAGCCACAACAACTGCAAGCCACGTc                                                                                                                                               >  1:1072713/1‑150 (MQ=255)
       acCCGGCTGTGGCGGTCACGGAACCACAACGGCGCCATCACCACGGTAATCACCAGTGCCAGAACCGCACCACGGCTCTGGCTGAGCATGGCGAAGGCCCCCAGGCAAGCCAGTGCCGCCAGCCACAACAACTGCAAGCCACGTCGCCGa                                                                                                                                          >  1:788658/1‑150 (MQ=255)
       acCCGGCTGTGGCGGTCACGGAACCACAACGGCGCCATCACCACGGTAATCACCAGTGCCAGAACCGCACCACGGCTCTGGCTGAGCATGGCGAAGGCCCCCAGGCAAGCCAGTGCCGCCAGCCACAACAACTGCAAGCCACGTCGCCGa                                                                                                                                          <  1:229902/150‑1 (MQ=255)
          cGGCTGTGGCGGTCACGGAACCACAACGGCGCCATCACCACGGTAATCACCAGTGCCAGAACCGCACCACGGCTCTGGCTGAGCATGGCGAAGGCCCCCAGGCAAGCCAGTGCCGCCAGCCACAACAACTGCAAGCCACGTCGCCGAGGt                                                                                                                                       <  2:2212808/150‑1 (MQ=255)
          cGGCTGTGGCGGTCAAGGAACCGCTACAGCGCAATCACCACGGCACTCACAAGTGCCAGCACCTCATCACGGCTCTGGGTGAGCATGGCGAACGCCCCCAGGCCAGCCAGTGCCGCCAGCCACAACAACTGCAAGCCACGTCGCCGAGGt                                                                                                                                       <  2:389971/150‑1 (MQ=255)
            gCTGTGGCGGTCACGGAACCACAACGGCGCCATCACCACGGTAATCACCAGTGCCAGAACCGCACCACGGCTCTGGCTGAGCATGGCGAAGGCCCCCAGGCAAGCCAGTGCCGCCAGCCACAACAACTGCAAGCCACGTCGCCGAGGTgg                                                                                                                                     <  2:1197553/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)
                           gAACCACAACGGCCCCATCACCACGGTAATCACCAGTGCCAGAACCGCACCACGGCTCTGGCTGAGCATGGCGAAGGCCCCCAGGCAAGCCAGTGCCGCCAGCCACAACAACTGCAAGCCACGTCGCCGAGGTGGCTCGTAGAGCATCaa                                                                                                                     <  1:1432042/150‑1 (MQ=255)
                             aCCACAACGGCGCCATCACCACGGTAATCACCAGTGCCAGAACCGCACCACGGCTCTGGCTGAGCATGGCGAAGGCCCCCAGGCAAGCCAGTGCCGCCAGCCACAACAACTGCAAGCCACGTCGCCGAGGTGGCTCGTAGAGCATCAaga                                                                                                                   <  2:2251977/150‑1 (MQ=255)
                                acaACGGCGCCATCACCACGGTAATCACCAGTGCCAGAACCGCACCACGGCTCTGGCTGAGCATGGCGAAGGCCCCCAGGCAAGCCAGTGCCGCCAGCCACAACAACTGCAAGCCACGTCGCCGAGGTGGCTCGTAGAGCCTCAAGAGCa                                                                                                                >  1:609122/1‑150 (MQ=255)
                                acaACGGCGCCATCACCACGGTAATCACCAGTGCCAGAACCGCACCACGGCTCTGGCTGAGCATGGCGAAGGCCCCCAGGCAAGCCAGTGCCGCCAGCCACAACAACTGCAAGCCACGTCGCCGAGGTGGCTCGTAGAGCATCAAGAGCa                                                                                                                >  1:2085526/1‑150 (MQ=255)
                                acaACGGCGCCATCACCACGGTAATCACCAGTGCCAGAACCGCACCACGGCTCTGGCTGAGCATGGCGAAGGCCCCCAGGCAAGCCAGTGCCGCCAGCCACAACAACTGCAAGCCACGTCGCCGAGGTGGCTCGTAGAGCATCAAGAGCa                                                                                                                >  1:2322287/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                                                                                                        >  1:690097/1‑150 (MQ=255)
                                        gcCATCACCACGGTAATCACCAGTGCCAGAACCGCACCACGGCTCTGGCTGAGCATGGCGAAGGCCCCCAGGCAAGCCAGTGCCGCCAGCCACAACAACTGCAAGCCACGTCGCCGAGGTGGCTCGTAGAGCATCAAGAGCAAGGCCGaa                                                                                                        <  2:304300/150‑1 (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:2115640/1‑150 (MQ=255)
                                                    gTAATCACCAGTGCCAGAACCGCACCACGGCTCTGGCTGAGCATGGCGAAGGCCCCCAGGCAAGCCAGTGCCGCCAGCCACAACAACTGCAAGCCACGTCGCCGAGGTGGCTCGTAGAGCATCAAGAGCAAGGCCGAACCAATGACATaa                                                                                            <  1:1189110/150‑1 (MQ=255)
                                                    gTAATCACCAGTGCCAGAACCGCACCACGGCTCTGGCTGAGCATGGCGAAGGCCCCCAGGCAAGCCAGTGCCGCCAGCCACAACAACTGCAAGCCACGTCGCCGAGGTGGCTCGTAGAGCATCAAGAGCAAGGCCGAACCAATGACATaa                                                                                            <  2:1849840/150‑1 (MQ=255)
                                                    gTAATCACCAGTGCCAGAACCGCACCACGGCTCTGGCTGAGCATGGCGAAGGCCCCCAGGCAAGCCAGTGCCGCCAGCCACAACAACTGCAAGCCACGTCGCCGAGGTGGCTCGTAGAGCATCAAGAGCAAGGCCGAACCAATGACATaa                                                                                            <  2:690097/150‑1 (MQ=255)
                                                    ccaaTCATCGGTGCCAGAACCGCACTACGGCTCTGGCTGAGCAAGGCGAAGGCCCCCGGGCAAGCCAGTGCCCCCAGCCACAACAACTGCAATCCACGTCGCCGAGGTGGCTCGTAGAGCATCAAGAGCAAGGCCGAACCAATGACATaa                                                                                            <  2:468911/148‑1 (MQ=255)
                                                       aTCACCAGTGCCAGAACCGCACCACGGCTCTGGCTGAGCATGGCGAAGGCCCCCAGGCAAGCCAGTGCCGCCAGCCACAACAACTGCAAGCCACGTCGCCGAGGTGGCTCGTa                                                                                                                               >  1:640420/1‑113 (MQ=255)
                                                       aTCACCAGTGCCAGAACCGCACCACGGCTCTGGCTGAGCATGGCGAAGGCCCCCAGGCAAGCCAGTGCCGCCAGCCACAACAACTGCAAGCCACGTCGCCGAGGTGGCTCGTa                                                                                                                               <  2:640420/113‑1 (MQ=255)
                                                        tCACCAGTGCCAGAACCGCACCACGGCTCTGGCTGAGCATGGCGAAGGCCCCCAGGCAAGCCAGTGCCGCCAGCCACAACAACTGCAAGCCACGTCGCCGAGGTGGCTCGTCGAGCATCAAGAGCAAGGCCGAACCAATGACATAAGccc                                                                                        <  2:1391063/150‑1 (MQ=255)
                                                          aCCAGTGCCAGAACCGCACCACGGCTCTGGCTGAGCATGGCGAAGGCCCCCAGGCAAGCCAGTGCCGCCAGCCACAACAACTGCAAGCCACCTCGGCGACGTGGCttg                                                                                                                                 <  1:2202570/108‑3 (MQ=255)
                                                          aCCAGTGCCAGAACCGCACCACGGCTCTGGCTGAGCATGGCGAAGGCCCCCAGGCAAGCCAGTGCCGCCAGCCACAACAACTGCAAGCCACCTCGGCGACGTGGCttg                                                                                                                                 >  2:2202570/1‑106 (MQ=255)
                                                           ccAGTGCCAGAACCGCACCACGGCTCTGGCTGAGCATGGCGAAGGCCCCCAGGCAAGCCAGTGCCGCCAGCCACAACAACTGCAAGCCACGTCGCCGAGGTGGCTCg                                                                                                                                 <  1:1824057/107‑1 (MQ=255)
                                                             aGTGCCAGTACCGACCCACGGCGCTTGCTTAGCATGGCGAAGGCCCCCAGGCAAGCCAGTGCCGCCAGCCACAACAACTGCAAGCCACGTCGCCGAGGTGGCTCGTAGAGCATCAAGAGCAAGGCCGAACCAATGACATAAGCCCCCAgg                                                                                   <  1:1698283/150‑1 (MQ=255)
                                                             aGTGCCAGAACCGCACCACGGCTCTGGCTGAGCATGGCGAAGGCCCCCAGGCAAGCCAGTGCCGCCAGCCACAACAACTGCAAGCCACGTCGCCGAGGTGGCTCGTAGAGCATCAAGAGCAAGGCCGAACCAATGACATAAGCCCCCAgg                                                                                   >  2:42092/1‑150 (MQ=255)
                                                                gCCAGAACCGCACCACGGCTCTGGCTGAGCATGGCGAAGGCCCCCAGGCAAGCCAGTGCCGCCAGCCACAACAACTGCAAGCCACGTCGCCGAGGTGGCTCGTAGAGCATCAAGAGCAAGGCCGAACCAATGACATAAGCCCCCAGGATg                                                                                <  1:1674978/150‑1 (MQ=255)
                                                                  cAGAACCGCACCACGGCTCTGGCTGAGCATGGCGAAGGCCCCCAGGCAAGCCAGTGCCGCCAGCCACAACAACTGCAAGCCACGTCGCCGAGGTGGCTCGTAGAGc                                                                                                                          <  1:984538/106‑1 (MQ=255)
                                                                     aaCCGCACCACGGCTCTGGCTGAGCATGGCGAAGGCCCCCAGGCAAGCCAGTGCCGCCAGCCACAACAACTGCAAGCCACGTCGCCGAGGTGGCTCGTAGAGCATCAAGAGCAAGGCCGAACCAATGACATAAGCCCCCAGGATGGGGTg                                                                           <  2:1547599/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:1046541/1‑150 (MQ=255)
                                                                         gCACCACGGCTCTGGCTGAGCATGGCGAAGGCCCCCAGGCAAGCCAGTGCCGCCAGCCACAACAACTGCAAGCCACGTCGCCGAGGTGGCTCGTAGAGCATCAAGAGCAAGGCCGAACCAATGACATAAGCCCCCAGGATGGGGTGGGat                                                                       <  1:2044672/150‑1 (MQ=255)
                                                                                ggCTCTGGCTGAGCATGGCGAAGGCCCCCAGGCAAGCCAGTGCCGCCAGCCACAACAACTGCAAGCCACGTCGCCGAGGTGGCTCGTAGAGCATCAAGAGCAAGGCCGAACCAATGACATAAGCCCCCAGGATGGGGTGGGATATCTCGc                                                                <  1:1387186/150‑1 (MQ=255)
                                                                                ggCTCTGGCTGAGCATGGCGAAGGCCCCCAGGCAAGCCAGTGCCGCCAGCCACAACAACTGCAAGCCACGTCGCCGAGGTGGCTCGTAGAGCATCAAGAGCAAGGCCGAACCAATGACATAAGCCCCCAGGATGGGGTGGGATATCTCGc                                                                <  2:2404757/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)
                                                                                          gAGCATGGCGAAGGCCCCCAGGCAAGCCAGTGCCGCCAGCCACAACAACTGCAAGCCACGTCGCCGAGGTGGCTCGTAGAGCATCAAGAGCAAGGCCGAACCAATGACATAAGCCCCCAGGATGGGGTGGGATATCTCGCCGATACCCGc                                                      >  2:404542/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)
                                                                                                     aGGCCCCCAGGCAAGCCAGTGCCGCCAGCCACAACAACTGCAAGCCACGTCGCCGAGGTGGCTCGTAGAGCATCAAGAGCAAGGCCGAACCAATGACATAAGCCCCCAGGATGGGGTGGGATAtc                                                                    >  1:2346772/1‑125 (MQ=255)
                                                                                                      ggCCCCCAGGCAAGCCAGTGCCGCCAGCCACAACAACTGCAAGCCACGTCGCCGAGGTGGCTCGTAGAGCATCAAGAGCAAGGCCGAACCAATGACATAAGCCCCCAGGATGGGGTGGGATATCTCGCCGATACCCGCCAGCCTGAACAg                                          >  1:486768/1‑150 (MQ=255)
                                                                                                        cccccAGGCAAGCCAGTGCCGCCAGCCACAACAACTGCAAGCCACGTCGCCGAGGTGGCTCGTAGAGCATCAAGAGCAAgg                                                                                                             >  1:2391846/1‑81 (MQ=255)
                                                                                                        cccccAGGCAAGCCAGTGCCGCCAGCCACAACAACTGCAAGCCACGTCGCCGAGGTGGCTCGTAGAGCATCAAGAGCAAgg                                                                                                             <  2:2391846/81‑1 (MQ=255)
                                                                                                         ccccAGGCAAGCCAGTGCCGCCAGCCACAACAACTGCAAGCCACGTCGCCGAGGTGGCTCGTAGAGCATCAAGAGCAAGGCCGAACCAATGACATAAGCCCCCAGGATGGGGTGGGATATCTCGCCGATACCCGCCAGc                                                  <  1:290830/139‑1 (MQ=255)
                                                                                                         ccccAGGCAAGCCAGTGCCGCCAGCCACAACAACTGCAAGCCACGTCGCCGAGGTGGCTCGTAGAGCATCAAGAGCAAGGCCGAACCAATGACATAAGCCCCCAGGATGGGGTGGGATATCTCGCCGATACCCGCCAGc                                                  >  2:290830/1‑139 (MQ=255)
                                                                                                         ccccAGGCAAGCCAGTGCCGCCAGCCACAACAACTGCAAGCCACGTCGCCGAGGTGGCTCGTAGAGCATCAAGAGCAAGGCCGAACCAATGACATAAGCCCCCAGGATGGGGTGGGATATCTCGCCGATACCCGCCAGCCTGAACAGCaa                                       <  2:486768/150‑1 (MQ=255)
                                                                                                         ccccAGGCAAGCCAGTGCCGCCAGCCACAACAACTGCAAGCCACGTCGCCGAGGTGGCTCGTAGAGCATCAAGAGCAAGGCCGAACCAATGACATAAGCCCCCAGGATGGGGTGGGATATCTCGCCGATACCCGCCAGCCTGAACAGCaa                                       >  2:1866340/1‑150 (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                                                                                        <  1:360744/100‑1 (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)
                                                                                                                  aGCCAGTGCCGCCAGCCACAACAACTGCAAGCCACGTCGCCGAGGTGGCTCGTAGAGCATCAAGAGCAAGGCCGAACCAATGACATAAGCCCCCAGGATGGGGTGGGATATCTCGCCGATACCCGCCAGCCTGAACAGCAAAGGCGCACg                              >  1:465045/1‑150 (MQ=255)
                                                                                                                         gccgccAGCCACAACAACTGCAAGCCACGTCGCCGAGGTGGCTCGTAGAGCATCAAGAGCAAGGCCGAACCAATGACATAAGCCCCCAGGATGGGGTGGGATATCTCGCCGATACCCGCCAGCCTGAACAGCAAAGGCGCACGCTGCAcc                       <  2:2115640/150‑1 (MQ=255)
                                                                                                                            gccagccaCAACAACTGCAAGCCACGTCGCCGAGGTGGCTCGTAGAGCATCAAGAGCAAGGCCGAACCAATGACATAAGCCCCCAGGATGGGGTGGGATATCTCGCCGATACCCGCCAGCCTGAACAGCAAAGGCGCACGCTGCACCCCa                    >  1:451011/1‑150 (MQ=255)
                                                                                                                             ccagccaCAACAACTGCAAGCCACGTCGCCGAGGTGGCTCGTAGAGCATCAAGAGCAAGGCCGAACCAATGACATAAGCCCCCAGGATGGGGTGGGATATCTCGCCGATACCCGCCAGCCTGAACAGCAAAGGCGCACGCTGCACCCCAt                   <  2:1263276/150‑1 (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)
                                                                                                                                      acaacTGCAAACCACGTCGCCGAGGTGGCTCGTAGAGCATCAAGAGCAAGGCCGAACCAATGACATAAGCCCCCAGGATGGGGTGGGATATCTCGCCGATACCCGCCAGCCTGAACAGCAAAGGCGCACGCTGCACCCCATAGAACTtga          <  2:1622294/150‑1 (MQ=255)
                                                                                                                                          cTGCAAGCCACGTCGCCGAGGTGGCTCGTAGAGCATCAAGAGCAAGGCCGAACCAATGACATAAGCCCCCAGGATGGGGTGGGATATCTCGCCGATACCCGCCAGCCTGAACAGCAAAGGCGCACGCTGCACCCCATAGAACTTgatgat      <  1:1866340/150‑1 (MQ=255)
                                                                                                                                           tacgaGCCACCTCGGCGACGTGGCTCGTAGAGCATCAAGAGCAAGGCCGAACCAATGACATAAGCCCCCAGGATGGGGTGGGATATCTCGCCGATACCCGCCAGCCTGAACAGCAAAGGCGCACGCTGCACCCCATAGAACTTGatgatg     >  2:1906281/5‑150 (MQ=255)
                                                                                                                                           tGCAAGCCACGTCGCCGAGGTGGCTCGTAGAGCATCAAGAGCAAGGCCGAACCAATGACATAAGCCCCCAGGATGGGGTgg                                                                          >  1:1435802/1‑81 (MQ=255)
                                                                                                                                           tGCAAGCCACGTCGCCGAGGTGGCTCGTAGAGCATCAAGAGCAAGGCCGAACCAATGACATAAGCCCCCAGGATGGGGTgg                                                                          <  2:1435802/81‑1 (MQ=255)
                                                                                                                                           tGCAAGCCACGTCGCCGAGGTGGCTCGTAGAGCATCAAGAGCAAGGCCGAACCAATGACATAAGCCCCCAGGATGGGGTGGGATATCTCGCCGATACCCGCCAGCCTGAACAGc                                         >  1:1017197/1‑114 (MQ=255)
                                                                                                                                           tGCAAGCCACGTCGCCGAGGTGGCTCGTAGAGCATCAAGAGCAAGGCCGAACCAATGACATAAGCCCCCAGGATGGGGTGGGATATCTCGCCGATACCCGCCAGCCTGAACAGc                                         <  2:1017197/114‑1 (MQ=255)
                                                                                                                                           tGCAAGCCACGTCGCCGAGGTGGCTCGTAGAGCATCAAGAGCAAGGCCGAACCAATGACATAAGCCCCCAGGATGGGGTGGGATATCTCGCCGATACCCGCCAGCCTGAACAGCAAAGGCGCACGCTGCACCCCATAGAACTTgatgat      <  1:790753/149‑1 (MQ=255)
                                                                                                                                           tGCAAGCCACGTCGCCGAGGTGGCTCGTAGAGCATCAAGAGCAAGGCCGAACCAATGACATAAGCCCCCAGGATGGGGTGGGATATCTCGCCGATACCCGCCAGCCTGAACAGCAAAGGCGCACGCTGCACCCCATAGAACTTgatgat      >  2:790753/1‑149 (MQ=255)
                                                                                                                                             cAAGCCACGTCGCCGAGGTGGCTCGTAGAGCATCAAGAGCAAGGCCGAACCAATGACATAAGCCCCCAGGATGGGGTGGGATATCTCGCCGATACCCGCCAGCCTGAACAGCAAAGGCGCACGCTGCACCCCATAGAACTTGATGATGga   >  2:1418265/1‑150 (MQ=255)
                                                                                                                                              gaGCCACCTCGGCGACGTGGCTCGTAGAGCATCAAGAGCAAGGCCGAACCAATGACATAAGCCCCCAGGATGGGGTGGGATATCTCGCCGATACCCGCCAGCCTGAAACGGGCGCTGGCACGCTGCACCCCATAGAACTTGATGATGgag  >  1:126203/2‑150 (MQ=255)
                                                                                                                                              aaGCCACGTCGCCGAGGTGGCTCGTAGAGCATCAAGAGCAAGGCCGAACCAATGACATAAGCCCCCAGGATGGGGTGGGATATCTCGCCGATACCCGCCAGCCTGAACAGCAAAgg                                    >  1:912520/1‑116 (MQ=255)
                                                                                                                                              aaGCCACGTCGCCGAGGTGGCTCGTAGAGCATCAAGAGCAAGGCCGAACCAATGACATAAGCCCCCAGGATGGGGTGGGATATCTCGCCGATACCCGCCAGCCTGAACAGCAAAgg                                    <  2:912520/116‑1 (MQ=255)
                                                                                                                                                     |                                                                                                                                               
AGAGAACACCCGGCTGTGGCGGTCACGGAACCACAACGGCGCCATCACCACGGTAATCACCAGTGCCAGAACCGCACCACGGCTCTGGCTGAGCATGGCGAAGGCCCCCAGGCAAGCCAGTGCCGCCAGCCACAACAACTGCAAGCCACGTCGCCGAGGTGGCTCGTAGAGCATCAAGAGCAAGGCCGAACCAATGACATAAGCCCCCAGGATGGGGTGGGATATCTCGCCGATACCCGCCAGCCTGAACAGCAAAGGCGCACGCTGCACCCCATAGAACTTGATGATGGAG  >  NC_002947/5618013‑5618304

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: