Predicted mutation
evidence seq id position mutation annotation gene description
RA NC_000913 3,789,047 C→T A316V (GCT→GTT) 
*345* (TAG→TAA
yibQ →
waaH ←
putative polysaccharide deacetylase
LPS(HepIII)‑glucuronic acid glycosyltransferase

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009133,789,0470CT100.0% 267.8 / NA 82A316V (GCT→GTT) 
*345* (TAG→TAA
yibQ
waaH
putative polysaccharide deacetylase
LPS(HepIII)‑glucuronic acid glycosyltransferase
Reads supporting (aligned to +/- strand):  ref base C (0/0);  new base T (41/41);  total (41/41)

ACCTGTTTATGCTAATCGCTTCTTTGAAGTATTAAGCGAAAGCATCAGCCAGAGCACGCTGATCGTTTACTTCCAGCATCAGTGGCAAGGCTGGGGCAAACAGCCCGAAGCGGCGAAGTTTAACGCTAGCGCAAATTAAAGGCGCTACGCGTAATGCGATGGTGCGACCGCGTTTTGTCGCGCCATTGCCATAAGCGGAAAGACCACAGTAACGCCTGGTAACCGACCTTCACGCTGCGCACATTGGTAATCA  >  NC_000913/3788922‑3789174
                                                                                                                             |                                                                                                                               
aCCTGTTTATGCTAATCGCTTCTTTGAAGTATTAAGCGAAAGCATCAGCCAGAGCACGCTGATCGTTTACTTCCAGCATCAGTGGCAAGGCTGGGGCAAACAGCCCGAAGCGGCGAAGTTTAACGTTAGCGCAAATTaa                                                                                                                    <  2:1032219/139‑1 (MQ=255)
aCCTGTTTATGCTAATCGCTTCTTTGAAGTATTAAGCGAAAGCATCAGCCAGAGCACGCTGATCGTTTACTTCCAGCATCAGTGGCAAGGCTGGGGCAAACAGCCCGAAGCGGCGAAGTTTAACGTTAGCGCAAATTaa                                                                                                                    >  1:586745/1‑139 (MQ=255)
      ttATGCTAATCGCTTCTTTGAAGTATTAAGCGAAAGCATCAGCCAGAGCACGCTGATCGTTTACTTCCAGCATCAGTGGCAAGGCTGGGGCAAACAGCCCGAAGCGGCGAAGTTTAACGTTAGCGCAAATTAAAGgcgc                                                                                                              >  2:1505371/1‑139 (MQ=255)
       tATGCTAATCACTTCTTTGAAGTATTAAGCGAAAGCATCAGCCAGAGCACGCTGATCGTTTACTTCCAGCATCAGTGGCAAGGCTGGGGCAAACAGCCCGAAGCGGCGAAGTGTAACGTTAGCGCAAATTAAAGGCGCt                                                                                                             <  1:1125649/139‑1 (MQ=255)
             aaTCGCTTCTTTGAAGTATTAAGCGAAAGCATCAGCCAGAGCACGCTGATCGTTTACTTCCAGCATCAGTGGCAAGGCTGGGGCAAACAGCCCGAAGCGGCGAAGTTTAACGTTAGCGCAAATTAAAGGCGCTACGCGt                                                                                                       >  2:849313/1‑139 (MQ=255)
               tCGCTTCTTTGAAGTATTAAGCGAAAGCATCAGCCAGAGCACGCTGATCGTTTACTTCCAGCATCAGTGGCAAGGCTGGGGCAAACAGCCCGAAGCGGCGAAGTTTAACGTTAGCGCAAATTAAAGGCGCTACGCGTaa                                                                                                     >  2:353836/1‑139 (MQ=255)
               tCGCTTCTTTGAAGTATTAAACGAAAGCATCAGCCAGAGCACGCTGATCGTTTACTTCCAGCATCAGTGGCAAGGCTGGGGCAAACAGCCCGAAGCGGCGAAGTTTAACGTTAGCGCAAATTAAAGGCGCTACGCGTaa                                                                                                     >  2:124060/1‑139 (MQ=255)
                            gTATTAAGCGAAAGCATCAGCCAGAGCACGCTGATCGTTTACTTCCAGCATCAGTGGCAAGGCTGGGGCAAACAGCCCGAAGCGGCGAAGTTTAACGTTAGCGCAAATTAAAGGCGCTACGCGTAATGCGATGGTGCGa                                                                                        <  1:1505371/139‑1 (MQ=255)
                             tATTAAGCGAAAGCATCAGCCAGAGCACGCTGATCGTTTACTTCCAGCATCAGTGGCAAGGCTGGGGCAAACAGCCCGAAGCGGCGAAGTTTAACGTTAGCGCAAATTAAAGGCGCTACGCGTAATGCGATGGTGCGAc                                                                                       <  1:438748/139‑1 (MQ=255)
                               ttAAGCGAAAGCATCAGCCCGAGCACGCTGATCGTTTACTTCCAGCATCAGTGGCAAGGCTGGGGCAAACAGCCCGAAGCGGCGAAGTTTAACGTTAgcg                                                                                                                            >  1:1004544/1‑100 (MQ=255)
                               ttAAGCGAAAGCATCAGCCAGAGCACGCTGATCGTTTACTTCCAGCATCAGTGGCAAGGCTGGGGCAAACAGCCCGAAGCGGCGAAGTTTAACGTTAgcg                                                                                                                            <  2:1004544/100‑1 (MQ=255)
                               ttAAGCGAAAGCATCAGCCAGAGCACGCTGATCGTTTACTTCCAGCATCAGTGGCAAGGCTGGGGCAAACAGCCCGAAGCGGCGAAGTTTAACGTTAGCGCAAATTAAAg                                                                                                                  <  1:808693/110‑1 (MQ=255)
                               ttAAGCGAAAGCATCAGCCAGAGCACGCTGATCGTTTACTTCCAGCATCAGTGGCAAGGCTGGGGCAAACAGCCCGAAGCGGCGAAGTTTAACGTTAGCGCAAATTAAAg                                                                                                                  >  2:808693/1‑110 (MQ=255)
                               ttAAGCGAAAGCATCAGCCAGAGCACGCTGATCGTTTACTTCCAGCATCAGTGGCAAGGCTGGGGCAAACAGCCCGAAGCGGCGAAGTTTAACGTTAGCGCAAATTAAAGGCGCTACGCGTAATGCGATGGTGCGAc                                                                                       <  2:1417030/137‑1 (MQ=255)
                               ttAAGCGAAAGCATCAGCCAGAGCACGCTGATCGTTTACTTCCAGCATCAGTGGCAAGGCTGGGGCAAACAGCCCGAAGCGGCGAAGTTTAACGTTAGCGCAAATTAAAGGCGCTACGCGTAATGCGATGGTGCGAc                                                                                       >  1:1417030/1‑137 (MQ=255)
                               ttAAGCGAAAGCATCAGCCAGAGCACGCTGATCGTTTACTTCCAGCATCAGTGGCAAGGCTGGGGCAAACAGCCCGAAGCGGCGAAGTTTAACGTTAGCGCAAATTAAAGGCGCTACGCGTAATGCGATGGTGCGACcg                                                                                     >  2:1073386/1‑139 (MQ=255)
                               ttAAGCGAAAGCATCAGCCAGAGCACGCTGATCGTTTACTTCCAGCATCAGTGGCAAGGCTGGGGCAAACAGCCCGAAGCGGCGAAGTTTAACGTTAGCGCAAATTAAAGGCGCTACGCGTAATGCGATGGTGCGACcg                                                                                     >  1:641181/1‑139 (MQ=255)
                               ttAAGCGAAAGCATCAGCCAGAGCACGCTGATCGTTTACTTCCAGCATCAGTGGCAAGGCTGGGGCAAACAGCCCGAAGCGGCGAAGTTTAACGTTAGCGCAAATTAAAGGCGCTACGCGTAATGCGATGGTGCGACcg                                                                                     >  1:1502180/1‑139 (MQ=255)
                                tAAGCGAAAGCATCAGCCAGAGCACGCTGATCGTTTACTTCCAGCATCAGTGGCAAGGCTGGGGCAAACAGCCCGAAGCGGCGAAGTTTAACGTTAGCGCAAATTAAAGGCGCTACGCGTAATGCGATGGTGCGACcgc                                                                                    >  1:994976/1‑139 (MQ=255)
                                  aGCGAAAGCATCAGCCAGAGCACGCTGATCGTTTACTTCCAGCATCAGTGGCAAGGCTGGGGCAAACAGCCCGAAGCGGCGAAGTTTAACGTTAGCGCAAATTAAAGGCGCTACGCGTAATGCGATGGTGCGACCGCGt                                                                                  <  1:1073386/139‑1 (MQ=255)
                                       aaGCATCAGCCAGAGCACGCTGATCGTTTACTTCCAGCATCAGTGGCAAGGCTGGGGCAAACAGCCCGAAGCGGCGAAGTTTAACGTTAGCGCAAATTAAAGGCGCt                                                                                                             >  1:1067849/1‑107 (MQ=255)
                                       aaGCATCAGCCAGAGCACGCTGATCGTTTACTTCCAGCATCAGTGGCAAGGCTGGGGCAAACAGCCCGAAGCGGCGAAGTTTAACGTTAGCGCAAATTAAAGGCGCt                                                                                                             <  2:1067849/107‑1 (MQ=255)
                                       aaGCATCAGCCAGAGCACGCTGATCGTTTACTTCCAGCATCAGTGGCAAGGCTGGGGCAAACAGCCCGAAGCGGCGAAGTTTAACGTTAGCGCAAATTAAAGGCGCTACGCGTAATGCGATGGTGCGACCGCGt                                                                                  <  2:1180571/134‑1 (MQ=255)
                                       aaGCATCAGCCAGAGCACGCTGATCGTTTACTTCCAGCATCAGTGGCAAGGCTGGGGCAAACAGCCCGAAGCGGCGAAGTTTAACGTTAGCGCAAATTAAAGGCGCTACGCGTAATGCGATGGTGCGACCGCGt                                                                                  >  1:1180571/1‑134 (MQ=255)
                                       aaGCATCAGCCAGAGCACGCTGATCGTTTACTTCCAGCATCAGTGGCAAGGCTGGGGCAAACAGCCCGAAGCGGCGAAGTTTAACGTTAGCGCAAATTAAAGGCGCTACGCGTAATGCGATGGTGCGACCGCGTTTTGt                                                                             >  1:753363/1‑139 (MQ=255)
                                           aTCAGCCAGAGCACGCTGATCGTTTACTTCCAGCATCAGTGGCAAGGCTGGGGCAAACAGCCCGAAGCGGCGAAGTTTAACGTTAGCGCAAATTAAAGGCGCTACGCGTAATGCGATGGTGCGACCGCGTTTTGTcgcg                                                                         <  1:353836/139‑1 (MQ=255)
                                           aTCAGCCAGAGCACGCTGATCGTTTACTTCCAGCATCAGTGGCAAGGCTGGGGCAAACAGCCCGAAGCGGCGAAGTTTAACGTTAGCGCAAATTAAAGGCGCTACGCGTAATGCGATGGTGCGACCGCGTTTTGTcgcg                                                                         <  2:580559/139‑1 (MQ=255)
                                            tCAGCCAGAGCACGCTGATCGTTTACTTCCAGCATCAGTGGCAAGGCTGGGGCAAACAGCCCGAAGCGGCGAAGTTTAACGTTAGCGCAAATTAAAGGCGCTACGCGTAATGCGATGGTGCGACCGCGTTTTGTCgcgc                                                                        >  1:897587/1‑139 (MQ=255)
                                                    agCACGCTGATCGTTTACTTCCAGCATCAGTGGCAAGGCTGGGGCAAACAGCCCGAAGTGGCGAAGTTTAACGTTAGCGCAAATTAAAGGCGCTACGCGTAATGCGATGGTGCGACCGCGTTTTGTCGCGCCATTGCCa                                                                >  1:1461162/1‑139 (MQ=255)
                                                    agCACGCTGATCGTTTACTTCCAGCATCAGTGGCAAGGCTGGGGCAAACAGCCCGAAGCGGCGAAGTTTAACGTTAg                                                                                                                              <  2:392184/77‑1 (MQ=255)
                                                    agCACGCTGATCGTTTACTTCCAGCATCAGTGGCAAGGCTGGGGCAAACAGCCCGAAGCGGCGAAGTTTAACGTTAg                                                                                                                              >  1:392184/1‑77 (MQ=255)
                                                    agCACGCTGATCGTTTACTTCCAGCATCAGTGGCAAGGCTGGGGCAAACAGCCCGAAGCGGCGAAGTTTAACGTTAGCGCAAATTAAAGGCGCTACGCGTAATGCGATGGTGCGACCGCGTTTTGTcgcg                                                                         <  1:501048/130‑1 (MQ=255)
                                                    agCACGCTGATCGTTTACTTCCAGCATCAGTGGCAAGGCTGGGGCAAACAGCCCGAAGCGGCGAAGTTTAACGTTAGCGCAAATTAAAGGCGCTACGCGTAATGCGATGGTGCGACCGCGTTTTGTcgcg                                                                         >  2:501048/1‑130 (MQ=255)
                                                    agCACGCTGATCGTTTACTTCCAGCATCAGTGGCAAGGCTGGGGCAAACAGCCCGAAGCGGCGAAGTTTAACGTTAGCGCAAATTAAAGGCGCTACGCGTAATGCGATGGTGCGACCGCGTTTTGTcg                                                                           <  1:187978/128‑1 (MQ=255)
                                                    agCACGCTGATCGTTTACTTCCAGCATCAGTGGCAAGGCTGGGGCAAACAGCCCGAAGCGGCGAAGTTTAACGTTAGCGCAAATTAAAGGCGCTACGCGTAATGCGATGGTGCGACCGCGTTTTGTcg                                                                           >  2:187978/1‑128 (MQ=255)
                                                     gCACGCTGATCGTTTACTTCCAGCATCAGTGGCAAGGCTGGGGCAAACAGCCCGAAGCGGCGAAGTTTAACGTTAGCGCAAATTAAAGGCGCTACGCGTAATGCGATGGTGCGACCGCGTTTTGTCGCGCCATTGCCAt                                                               <  2:994976/139‑1 (MQ=255)
                                                     gCACGCTGATCGTTTACTTCCAGCATCAGTGGCAAGGCTGGGGCAAACAGCCCGAAGCGGCGAAGTTTAACGTTAGCGCAAATTAAAGGCGCTACGCGTAATGCGATGGTGCGACCGCGTTTTGTCGCGCCATTGCCAt                                                               >  1:998132/1‑139 (MQ=255)
                                                         gCTGATCGTTTACTTCCAGCATCAGTGGCAAGGCTGGGGCAAACAGCCCGAAGCGGCGAAGTTTAACGTTAGCGCAAATTAAAGGCGCTACGCGTAATGCGATGGTGCGACCGCGTTTTGTCGCGCCATTGCCATAAGc                                                           >  2:461837/1‑139 (MQ=255)
                                                          cTGATCGTTTACTTCCAGCATCAGTGGCAAGGCTGGGGCAAACAGCCCGAAGCGGCGAAGTTTAACGTTAGCGCAAATTAAAGGCGCTACGCGTAATGCGATGGTGCGACCGCGTTTTGTCGCGCCATTGCCATAAGCg                                                          >  2:126055/1‑139 (MQ=255)
                                                              tCGTTTACTTCCAGCATCAGTGGCAAGGCTGGGGCAAACAGCCCGAAGCGGCGAAGTTTAACGTTAgcg                                                                                                                            >  2:1452017/1‑69 (MQ=255)
                                                              tCGTTTACTTCCAGCATCAGTGGCAAGGCTGGGGCAAACAGCCCGAAGCGGCGAAGTTTAACGTTAgcg                                                                                                                            <  1:1452017/69‑1 (MQ=255)
                                                                 tttACTTCCAGCATCAGTGGCAAGGCTGGGGCAAACAGCCCGAAGCGGCGAAGTTTAACGTTAGCGCAAATTAAAGGCGCTACGCGTAATGCGATGGTGCGACCGCGTTTTGTCGCGCCATTGCCATAAGCGGAAAGAc                                                   <  2:586745/139‑1 (MQ=255)
                                                                   tACTTCCAGCATCAGTGGCAAGGCTGGGGCAAACAGCCCGAAGCGGCGAAGTTTAACGTTAGCGCAAATTAAAGGCGCTACGCGTAATGCGATGGTGCGACCGCGTTTTGTcg                                                                           >  2:1150466/1‑113 (MQ=255)
                                                                   tACTTCCAGCATCAGTGGCAAGGCTGGGGCAAACAGCCCGAAGCGGCGAAGTTTAACGTTAGCGCAAATTAAAGGCGCTACGCGTAATGCGATGGTGCGACCGCGTTTTGTcg                                                                           <  1:1150466/113‑1 (MQ=255)
                                                                       tCCAGCATCAGTGGCAAGGCTGGGGCAAACAGCCCGAAGCGGCGAAGTTTAACGTTAGCGCAAATTAAAGGCGCTACGCGTAATGCGATGGTGCGACCGCGTTTTGTCGCGcc                                                                       >  1:37594/1‑113 (MQ=255)
                                                                       tCCAGCATCAGTGGCAAGGCTGGGGCAAACAGCCCGAAGCGGCGAAGTTTAACGTTAGCGCAAATTAAAGGCGCTACGCGTAATGCGATGGTGCGACCGCGTTTTGTCGCGcc                                                                       <  2:37594/113‑1 (MQ=255)
                                                                       tCCAGCATCAGTGGCAAGGCTGGGGCAAACAGCCCGAAGCGGCGAAGTTTAACGTTAGCGCAAATTAAAGGCGCTACGCGTAATGCGATGCTGCGACCGCGTTTTGTCGCGCCATTGCCATAAGCGGAAAGACCACAGt                                             >  2:469179/1‑139 (MQ=255)
                                                                         cAGCATCAGTGGCAAGGCTGGGGCAAACAGCCCGAAGCGGCGAAGTTTAACGTTAGCGCAAATTAAAGGCGCTACGCGTAATGCGATGGTGCGACCGCGTTTTGTCGCGCCATTGCCATAAGCGGAAAGACCACAGTaa                                           <  2:641181/139‑1 (MQ=255)
                                                                          aGCATCAGTGGCAAGGCTGGGGCAAACAGCCCGAAGCGGCGAAGTTTAACGTTAGCGCAAATTAAAGGCGCTACGCGTAATGCGATGGTGCGACCGCGTTTTGTCGCGCCATTGCCATAAGCGGAAAGACCACAGTAAc                                          <  1:126055/139‑1 (MQ=255)
                                                                           gCATCAGTGGCAAGGCTGGGGCAAACAGCCCGAAGCGGCGAAGTTTAACGTTAGCGCaa                                                                                                                         >  2:400493/1‑59 (MQ=255)
                                                                           gCATCAGTGGCAAGGCTGGGGCAAACAGCCCGAAGCGGCGAAGTTTAACGTTAGCGCaa                                                                                                                         <  1:400493/59‑1 (MQ=255)
                                                                            cATCAGTGGCAAGGCTGGGGCAAACAGCCCGAAGCGGCGAAGTTTAACGTTAGCGCAAATTAAAGGCGCTACGCGTAATGCGATGGTGCGACCGCGTTTTGTCGCGCCATTGCCATAAGCGGAAAGACCACAGTAACGc                                        >  2:737478/1‑139 (MQ=255)
                                                                                aGTGGCAAGGCTGGGGCAAACAGCCCGAAGCGGCGAAGTTTAACGTTAGCGCAAATTAAAGGCGCTACGCGTAATGCGATGGTGCGACCGCGTTTTGTCGCGCCATTGCCATAAGCGGAAAGACCACAGTAACGCCTgg                                    >  1:7064/1‑139 (MQ=255)
                                                                                  tGGCAAGGCTGGGGCAAACAGCCCGAAGCGGCGAAGTTTAACGTTAGCGCAAATTAAAGGCGCTACGCGTAATGCGATGGTGCGAc                                                                                       >  2:243271/1‑86 (MQ=255)
                                                                                  tGGCAAGGCTGGGGCAAACAGCCCGAAGCGGCGAAGTTTAACGTTAGCGCAAATTAAAGGCGCTACGCGTAATGCGATGGTGCGAc                                                                                       <  1:243271/86‑1 (MQ=255)
                                                                                   ggCAAGGCTGGGGCAAACAGCCCGAAGCGGCGAAGTTTAACGTTAGCGCAAATTAAAGGCGCTACGCGTAATGCGATGGTGCGACCGCGTTTTGTCGCGCCATTGTCATAAGCGGAAAGACCACAGTAACGCCTGGTaa                                 >  2:1474727/1‑139 (MQ=255)
                                                                                   ggCAAGGCTGGGGCAAACAGCCCGAAGCGGCGAAGTTTAACGTTAGCGCAAATTAAAGGCGCTACGCGTAATGCGATGGTGCGACCGCGTTTTGTCGCGCCATTGCCATAAGCGGAAAGACCACAGTAACGCCTGGTaa                                 <  2:897587/139‑1 (MQ=255)
                                                                                         gCTGGGGCAAACAGCCCGAAGCGGCGAAGTTTAACGTTAGCGCAAATTAAAGGCGCTAcg                                                                                                          >  2:1139489/1‑60 (MQ=255)
                                                                                         gCTGGGGCAAACAGCCCGAAGCGGCGAAGTTTAACGTTAGCGCAAATTAAAGGCGCTAcg                                                                                                          <  1:1139489/60‑1 (MQ=255)
                                                                                               gCAAACAGCCCGAAGTGGCGAAGTTTAACGTTAGCGCAAATTAAAGGCGCTACGCGTAATGCGATGGTGCGACCGCGTTTTGTCGCGCCATTGCCATAAGCGGAAAGACCACAGTAACGCCTGGTAACCGACCTTCACg                     <  2:1461162/139‑1 (MQ=255)
                                                                                                 aaaCAGCCCGAAGCGGCGAAGTTTAACGTTAGCGCAAATTAAAGGCGCTACGCGTAATGCGATGGTGCGACCGCGTTTTGTCGCGCCATTGCCATAAGCGGAAAGACCACAGTAACGCCTGGTAACCGACCTTCACGCt                   <  1:737478/139‑1 (MQ=255)
                                                                                                  aaCAGCCCGAAGCGGCGAAGTTTAACGTTAGCGCAAATTAAAGGCGCTACGCGTAATGCGATGGTGCGACCGCGTTTTGTCGCGCCATTGTCATAAGCGGAAAGACCACAGTAACGCCTGGTAACCGACCTTCACGCTg                  <  1:1474727/139‑1 (MQ=255)
                                                                                                   aCAGCCCGAAGCGGCGAAGTTTAACGTTAGCGCAAATTAAAGGCGCTACGCGTAATGCGATGGTGCGACCGCGTTTTGTCGCGCCATTGCCATAAGCGGAAAGACCACAGTAACGCCTGGTAACCGACCTTCACGCTgc                 >  1:825847/1‑139 (MQ=255)
                                                                                                    cAGCCCGAAGCGGCGAAGTTTAACGTTAGCGCAAATTAAAGGCGCTACGCGTAATGCGATGGTGCGACCGCGTTTTGTCGCGCCATTGCCATAAGCGGAAAGACCACAGTAACGCCTGGTAACCGACCTTCACGCTgcg                <  2:753363/139‑1 (MQ=255)
                                                                                                      gCCCGAAGCGGCGAAGTTTAACGTTAGCGCAAATTAAAGGCGCTACGCGTAATGCGATGGTGCGACCGCGTTTTGTCGCGCCATTGCCATAAGCGGAAAGACCACAGTAACGCCTGGTAACCGACCTTCACGCTGCGca              <  2:998132/139‑1 (MQ=255)
                                                                                                        ccGAAGCGGCGAAGTTTAACGTTAGCGCAAATTAAAGGCGCTa                                                                                                            >  2:594455/1‑43 (MQ=255)
                                                                                                        ccGAAGCGGCGAAGTTTAACGTTAGCGCAAATTAAAGGCGCTa                                                                                                            <  1:594455/43‑1 (MQ=255)
                                                                                                         cGAAGCGGCGAAGTTTAACGTTAGCGCAAATTAAAGGCGCTACGCGTAATGCGATGGTGCGACCGCGTTTTGTCGCGCCATTGCCATAAGCGGAAAGACCACAGTAACGCCTGGTAAcc                               <  2:1003597/119‑1 (MQ=255)
                                                                                                         cGAAGCGGCGAAGTTTAACGTTAGCGCAAATTAAAGGCGCTACGCGTAATGCGATGGTGCGACCGCGTTTTGTCGCGCCATTGCCATAAGCGGAAAGACCACAGTAACGCCTGGTAAcc                               >  1:1003597/1‑119 (MQ=255)
                                                                                                            aGCGGCGAAGTTTAACGTTAGCGCAAATTAAAGGCGCTACGCGTAATGCGATGGTGCGACCGCGTTTTGTCGCGCCATTGCCATAAGCGGAAAGACCACAGTAACGCCTGGTAACCGACCTTCACGCTGCGCACATTgg        >  2:773671/1‑139 (MQ=255)
                                                                                                                 cgAAGTTTAACGTTAGCGCAAATTAAAGGCGCTACGCGTAATGCGATGGTGCGACCGCGTTTTGTCGCGCCATTGCCATAAGc                                                           >  2:590565/1‑83 (MQ=255)
                                                                                                                 cgAAGTTTAACGTTAGCGCAAATTAAAGGCGCTACGCGTAATGCGATGGTGCGACCGCGTTTTGTCGCGCCATTGCCATAAGc                                                           <  1:590565/83‑1 (MQ=255)
                                                                                                                 cgAAGTTTAACGTTAGCGCAAATTAAAGGCGCTACGCGTAATGCGATGGTGCGACCGCGTTTTGTCGCGCCATTGCCATAAGCGGaaa                                                      >  2:1443440/1‑88 (MQ=255)
                                                                                                                 cgAAGTTTAACGTTAGCGCAAATTAAAGGCGCTACGCGTAATGCGATGGTGCGACCGCGTTTTGTCGCGCCATTGCCATAAGCGGaaa                                                      <  1:1443440/88‑1 (MQ=255)
                                                                                                                 cgAAGTTTAACGTTAGCGCAAATTAAAGGCGCTACGCGTAATGCGATGGTGCGACCGCGTTTTGTCGCGCCATTGCCATAAGCGGAAAGACCACAGTAACGCCTGGTAACCGACCTTCACGCTGCGCACATTg         <  1:356871/133‑1 (MQ=255)
                                                                                                                 cgAAGTTTAACGTTAGCGCAAATTAAAGGCGCTACGCGTAATGCGATGGTCCGACCGCGTTTTGTCGCGCCATTGCCATAAGCGGAAAGACCACAGTAACGCCTGGTAACCGACCTTCACGCTGCGCACATTg         >  2:356871/1‑133 (MQ=255)
                                                                                                                  gAAGTTTAACGTTAGCGCAAATTAAAGGCGCTACGCGTAATGCGATGGTGCGACCGCGTTTTGTCGCGCCATTGCCATAAGCGGAAAGACCACAGTAACGCCTGGTAACCGACCTTCACGCTGCGCACATTGGTAATCa  <  2:1502180/139‑1 (MQ=255)
                                                                                                                  gAAGTTTAACGTTAGCGCAAATTAAAGGCGCTACGCGTAATGCGATGGTGCGACCGCGTTTTGTCGCGCCATTGCCATAAGCGGAAAGACCACAGTAACGCCTGGTAACCGACCTTCACGCTGCGCACATTGGTAATCa  <  1:469179/139‑1 (MQ=255)
                                                                                                                  gAAGTTTAACGTTAGCGCAAATTAAAGGCGCTACGCGTAATGCGATGGTGCGACCGCGTTTTGTCGCGCCATTGCCATAAGCGGAAAGACCACAGTAACGCCTGGTAACCGACCTTCACGCTGCGCACATTGGTAATCa  <  1:773671/139‑1 (MQ=255)
                                                                                                                  gAAGTTTAACGTTAGCGCAAATTAAAGGCGCTACGCGTAATGCGATGGTGCGACCGCGTTTTGTCGCGCCATTGCCATAAGCGGAAAGACCACAGTAACGCCTGGTAACCGACCTTCACGCTGCGCACATTGGTAATCa  <  1:849313/139‑1 (MQ=255)
                                                                                                                      tttAACGTTAGCGCAAATTAAAGGCGCTACGCGTAATGCGATGGTGCGACCGCGTTTTGTCGCGCCATTGCCATAAGCgg                                                         >  1:70483/1‑80 (MQ=255)
                                                                                                                      tttAACGTTAGCGCAAATTAAAGGCGCTACGCGTAATGCGATGGTGCGACCGCGTTTTGTCGCGCCATTGCCATAAGCgg                                                         <  2:70483/80‑1 (MQ=255)
                                                                                                                             |                                                                                                                               
ACCTGTTTATGCTAATCGCTTCTTTGAAGTATTAAGCGAAAGCATCAGCCAGAGCACGCTGATCGTTTACTTCCAGCATCAGTGGCAAGGCTGGGGCAAACAGCCCGAAGCGGCGAAGTTTAACGCTAGCGCAAATTAAAGGCGCTACGCGTAATGCGATGGTGCGACCGCGTTTTGTCGCGCCATTGCCATAAGCGGAAAGACCACAGTAACGCCTGGTAACCGACCTTCACGCTGCGCACATTGGTAATCA  >  NC_000913/3788922‑3789174

Alignment Legend
Aligned base mismatch/match (shaded by quality score): ATCG/ATCG < 3 ≤ ATCG/ATCG < 16 ≤ ATCG/ATCG < 34 ≤ ATCG/ATCG < 38 ≤ ATCG/ATCG < 40 ≤ ATCG/ATCG
Unaligned base: atcg    Masked matching base: atcg    Alignment gap:     Deleted base: