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% 236.2 / NA 72A316V (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 (38/34);  total (38/34)

TGTTTATGCTAATCGCTTCTTTGAAGTATTAAGCGAAAGCATCAGCCAGAGCACGCTGATCGTTTACTTCCAGCATCAGTGGCAAGGCTGGGGCAAACAGCCCGAAGCGGCGAAGTTTAACGCTAGCGCAAATTAAAGGCGCTACGCGTAATGCGATGGTGCGACCGCGTTTTGTCGCGCCATTGCCATAAGCGGAAAGACCACAGTAACGCCTGGTAACCGACCTTCACGCTGCGCACATTGGTAATCAGGCG  >  NC_000913/3788925‑3789178
                                                                                                                          |                                                                                                                                   
tGTTTATCCTAATCGCTTCTTTGAAGTATTAAGCGAAAGCATCAGCCAGAGCACGCTGATCGTTTACTTCCAGCATCAGTGGCAAGGCTGGGGCAAACAGCCCGAAGCGGCGAAGTTTAACGTTAGCGCAAATTAAAgg                                                                                                                     >  1:1036294/1‑139 (MQ=255)
          aaTCGCTTCTTTGAAGTATTAAGCGAAAGCATCAGCCAGAGCACGCTGATCGTTTACTTCCAGCATCAGTGGCAAGGCTGGGGCAAACAGCCCGAAGCGGCGAAGTTTAACGTTAGCGCAAATTAAAGGCGCTACGCGt                                                                                                           >  1:1119878/1‑139 (MQ=255)
            tCGCTTCTTTGAAGTATTAAGCGAAAGCATCAGCCAGAGCACGCTGATCGTTTACTTCCAGCATCAGTGGCAAGGCTGGGGCAAACAGCCCGAAGCGGCGAAGTTTAACGTTAGCGCAAATTAAAGGCGCTACGCGTaa                                                                                                         >  1:63877/1‑139 (MQ=255)
            tCGCTTCTTTGAAGTATTAAGCGAAAGCATCAGCCAGAGCACGCTGATCGTTTACTTCCAGCATCAGTGGCAAGGCTGGGGCAAACAGCCCGAAGCGGCGAAGTTTAACGTTAGCGCAAATTAAAGGCGCTACGCGTaa                                                                                                         >  2:787653/1‑139 (MQ=255)
                 tcttTGAAGTATTAAGCGAAAGCATCAGCCAGAGCACGCTGATCGTTTACTTCCAGCATCAGTGGCAAGGCTGGGGCAAACAGCCCGAAGCGGCGAAGTTTAACGTTAGCGCAAATTAAAGGCGCTACGCGTAATGCGa                                                                                                    >  2:675959/1‑139 (MQ=255)
                            ttAAGCGAAAGCATCAGCCAGAGCACGCTGATCGTTTACTTCCAGCATCAGTGGCAAGGCTGGGGCAAACAGCCCGAAGCGGCGAAGTTTAACGTTAgcg                                                                                                                                >  2:1556085/1‑100 (MQ=255)
                            ttAAGCGAAAGCATCAGCCAGAGCACGCTGATCGTTTACTTCCAGCATCAGTGGCAAGGCTGGGGCAAACAGCCCGAAGCGGCGAAGTTTAACGTTAgcg                                                                                                                                <  1:1556085/100‑1 (MQ=255)
                              aaGCGAAAGCATCAGCCAGAGCACGCTGATCGTTTACTTCCAGCATCAGTGGCAAGGCTGGGCCAAACAGCCCGAAGCGGCGAAGTTTAACGTTAGCGCAAATTAAAGGCGCTACGCGTAATGCGATGGTACGACcgcg                                                                                       >  2:959462/1‑139 (MQ=255)
                                    aaGCATCAGCCAGAGCACGCTGATCGTTTACTTCCAGCATCAGTGGCAAGGCTGGGGCAAACAGCCCGAAGCGGCGAAGTTTAACGTTAGCGCAAATTAAAGGCGCTACGCGTAATGCGATGGTGCGACCGCGTTTTGt                                                                                 <  1:807470/139‑1 (MQ=255)
                                          cAGCCAGAGCACGCTGATCGTTTACTTCCAGCATCAGTGGCAAGGCTGGGGCAAACAGCCCGAAGCGGCGAAGTTTAACGTTAGCGCAAATTAAAGGCGCTACGCGTAATGCGATGGTGCGACCGCGTTTTGTCGCGcc                                                                           <  1:1545495/139‑1 (MQ=255)
                                           aGCCAGAGCACGCTGATCGTTTACTTCCAGCATCAGTGGCAAGGCTGGGGCAAACAGCCCGAAGCGGCGAAGTTTAACGTTAGCGCAAATTAAAg                                                                                                                      <  2:916284/95‑1 (MQ=255)
                                           aGCCAGAGCACGCTGATCGTTTACTTCCAGCATCAGTGGCAAGGCTGGGGCAAACAGCCCGAAGCGGCGAAGTTTAACGTTAGCGCAAATTAAAg                                                                                                                      >  1:916284/1‑95 (MQ=255)
                                                 agCACGCTGATCGTTTACTTCCAGCATCAGTGGCAAGGCTGGGGCAAACAGCCCGAAGCGGCGAAGTTTAACGTTAGCGCAAATTAAAGGCGCTACGCGTAATGCGATGGTGCGACCGCGTTTTGTCGCGCCATTGCCa                                                                    <  2:1036294/139‑1 (MQ=255)
                                                 agCACGCTGATCGTTTACTTCCAGCATCAGTGGCAAGGCTGGGGCAAACAGCCCGAAGCGGCGAAGTTTAACGTTAGCGCAAATTAAAGGCGCTACGCGTAATGCGATGGTGCGACCGCGTTTTGTCGCGCCATTGCCa                                                                    >  2:1225308/1‑139 (MQ=255)
                                                   cACGCTGATCGTTTACTTCCAGCATCAGTGGCAAGGCTGGGGCAAACAGCCCGAAGCGGCGAAGTTTAACGTTAGCGCAAATTAAAGGCGCTACGCGTAATGCGATGGTGCGACCGCGTTTTGTCGCGCCATTGCCATa                                                                  >  2:1441328/1‑139 (MQ=255)
                                                      gCTGATCGTTTACTTCCAGCATCAGTGGCAAGGCTGGGGCAAACAGCCCGAAGCGGCGAAGTTTAACGTTAGCGCAAATTAAAGGCGCTACGCGTAATGCGATGGTGCGACCGCGt                                                                                      <  1:984005/116‑1 (MQ=255)
                                                      gCTGATCGTTTACTTCCAGCATCAGTGGCAAGGCTGGGGCAAACAGCCCGAAGCGGCGAAGTTTAACGTTAGCGCAAATTAAAGGCGCTACGCGTAATGCGATGGTGCGACCGCGt                                                                                      >  2:984005/1‑116 (MQ=255)
                                                      gCTGATCGTTTACTTCCAGCATCAGTGGCAAGGCTGGGGCAAACAGCCCGAAGCGGCGAAGTTTAACGTTAGCGCAAATTAAAGGCGCTACGCGTAATGCGATGGTGCGACCGCGTTTTGTCGCGCCATTGCCATaa                                                                 <  2:772872/137‑1 (MQ=255)
                                                      gCTGATCGTTTACTTCCAGCATCAGTGGCAAGGCTGGGGCAAACAGCCCGAAGCGGCGAAGTTTAACGTTAGCGCAAATTAAAGGCGCTACGCGTAATGCGATGGTGCGACCGCGTTTTGTCGCGCCATTGCCATaa                                                                 >  1:772872/1‑137 (MQ=255)
                                                      gCTGATCGTTTACTTCCAGCATCAGTGGCAAGGCTGGGGCAAACAGCCCGAAGCGGCGAAGTTTAACGTTAGCGCAAATTAAAGGCGCTACGCGTAATGCGATGGTGCGACCGCGTTTTGTCGCGCCATTGCCATAAGc                                                               >  2:1541109/1‑139 (MQ=255)
                                                      gCTGATCGTTTACTTCCAGCATCAGTGGCAAGGCTGGGGCAAACAGCCCGAAGCGGCGAAGTTTAACGTTAGCGCAAATTAAAGGCGCTACGCGTAATGCGATGGTGCGACCGCGTTTTGTCGCGCCATTGCCATAAGc                                                               >  1:891267/1‑139 (MQ=255)
                                                      gCTGATCGTTTACTTCCAGCATCAGTGGCAAGGCTGGGGCAAACAGCCCGAAGCGGCGAAGTTTAACGTTAGCGCAAATTAAAGGCGCTACGCGTAATGCGATGGTGCGACCGCGTTTTGTCGCGCCATTGCCATAAGc                                                               >  1:1046198/1‑139 (MQ=255)
                                                       cTGATCGTTTACTTCCAGCATCAGTGGCAAGGCTGGGGCAAACAGCCCGAAGCGGCGAAGTTTAACGTTAGCGCAAATTaaa                                                                                                                       <  2:1408951/82‑1 (MQ=255)
                                                       cTGATCGTTTACTTCCAGCATCAGTGGCAAGGCTGGGGCAAACAGCCCGAAGCGGCGAAGTTTAACGTTAGCGCAAATTaaa                                                                                                                       >  1:1408951/1‑82 (MQ=255)
                                                       cTGATCGTTTACTTCCAGCATCAGTGGCAAGGCTGGGGCAAACAGCCCGAAGCGGCGAAGTTTAACGTTAGCGCAAATTAAAGGCGCTACGCGTAATGCGATGGTGCGa                                                                                            <  2:773676/109‑1 (MQ=255)
                                                       cTGATCGTTTACTTCCAGCATCAGTGGCAAGGCTGGGGCAAACAGCCCGAAGCGGCGAAGTTTAACGTTAGCGCAAATTAAAGGCGCTACGCGTAATGCGATGGTGCGa                                                                                            <  1:231276/109‑1 (MQ=255)
                                                       cTGATCGTTTACTTCCAGCATCAGTGGCAAGGCTGGGGCAAACAGCCCGAAGCGGCGAAGTTTAACGTTAGCGCAAATTAAAGGCGCTACGCGTAATGCGATGGTGCGa                                                                                            >  2:231276/1‑109 (MQ=255)
                                                       cTGATCGTTTACTTCCAGCATCAGTGGCAAGGCTGGGGCAAACAGCCCGAAGCGGCGAAGTTTAACGTTAGCGCAAATTAAAGGCGCTACGCGTAATGCGATGGTGCGa                                                                                            >  1:773676/1‑109 (MQ=255)
                                                       cTGATCGTTTACTTCCAGCATCAGTGGCAAGGCTGGGGCAAACAGCCCGAAGCGGCGAAGTTTAACGTTAGCGCAAATTAAAGGCGCTACGCGTAATGCGATGGTGCGA‑CGCGTTTTGTCGCGCCATTGCCATAAGCgg                                                             >  2:944737/1‑139 (MQ=255)
                                                         gATCGTTTACTTCCAGCATCAGTGGCAAGGCTGGGGCAAACAGCCCGAAGCGGCGAAGTTTAACGTTAGCGCAAATTAAAGGCGCTACGCGTAATGCGATGGTGCGACCGCGTTTTGTCGCGCCATTGCCATAAGCGGa                                                            <  1:258361/139‑1 (MQ=255)
                                                         gATCGTTTACTTCCAGCATCAGTGGCAAGGCTGGGGCAAACAGCCCGAAGCGGCGAAGTTTAACGTTAGCGCAAATTAAAGGCGCTACGCGTAATGCGATGGTGCGACCGCGTTTTGTCGCGCCATTGCCATAAGCGGa                                                            <  2:1119878/139‑1 (MQ=255)
                                                         gATCGTTTACTTCCAGCATCAGTGGCAAGGCTGGGGCAAACAGCCCGAAGCGGCGAAGTTTAACGTTAGCGCAAATTAAAGGCGCTACGCGTAATGCGATGGTGCGACCGCGTTTTGTCGCGCCATTGCCATAAGCGGa                                                            <  1:1441328/139‑1 (MQ=255)
                                                          aTCGTTTACTTCCAGCATCAGTGGCAAGGCTGGGGCAAACAGCCCGAAGCGGCGAAGTTTAACGTTAGCGCAAATTa                                                                                                                         >  1:1391259/1‑77 (MQ=255)
                                                          aTCGTTTACTTCCAGCATCAGTGGCAAGGCTGGGGCAAACAGCCCGAAGCGGCGAAGTTTAACGTTAGCGCAAATTa                                                                                                                         <  2:1391259/77‑1 (MQ=255)
                                                              tttACTTCCAGCATCAGTGGCAAGGCTGGGGCAAACAGCCCGAAGCGGCGAAGTTTAACGTTAGCGCAAATTAAAGGCGCTACGCGTAATGCGATGGTGCGACCGCGTTTTGTCGCGCCATTGCCATAAGCGGAAAGAc                                                       <  2:1389641/139‑1 (MQ=255)
                                                               ttACTTCCAGCATCAGTGGCAAGGCTGGGGCAAACAGCCCGAAGCGGCGAAGTTTAACGTTAGCGCAAATTAAAGGCGCTACGCGTAATGCGATGGTGCGACCGCGTTTTGTCGCGCCATTGCCATAAGCGGAAAGAcc                                                      <  1:649828/139‑1 (MQ=255)
                                                                tACTTCCAGCATCAGTGGCAAGGCTGGGGCAAACAGCCCGAAGCGGCGAAGTTTAACGTTAGCGCAAATTAAAGGCGCTACGCGTAATGCGATGGTGCGACCGCGTTTTGTCGCGCCATTGCCATAAGCGGAAAGACca                                                     >  2:355383/1‑139 (MQ=255)
                                                                    tCCAGCATCAGTGGCAAGGCTGGGGCAAACAGCCCGAAGCGGCGAAGTTTAACGTTAGCGCAAATTAAAGgcgc                                                                                                                  >  1:1306737/1‑74 (MQ=255)
                                                                    tCCAGCATCAGTGGCAAGGCTGGGGCAAACAGCCCGAAGCGGCGAAGTTTAACGTTAGCGCAAATTAAAGgcgc                                                                                                                  <  2:1306737/74‑1 (MQ=255)
                                                                    tCCAGCATCAGTGGCAAGGCTGGGGCAAACAGCCCGAAGCGGCGAAGTTTAACGTTAGCGCAAATTAAAGGCGCTACGCGTAATGCGATGGTGCGACCGCGTTTTGTCGCGCCATTGCCATAAGCGGAAAGACCACAGt                                                 <  2:891267/139‑1 (MQ=255)
                                                                    tCCAGCATCAGTGGCAAGGCTGGGGCAAACAGCCCGAAGCGGCGAAGTTTAACGTTAGCGCAAATTAAAGGCGCTACGCGTAATGCGATGGTGCGACCGCGTTTTGTCGCGCCATTGCCATAAGCGGAAAGACCACAGt                                                 >  2:1198092/1‑139 (MQ=255)
                                                                     ccAGCATCAGTGGCAAGGCTGGGGCAAACAGCCCGAAGCGGCGAAGTTTAACGTTAGCGCAAATTAAAGGCGCTACGCGTAATGCGATGGTGCGACCGCGTTTTGTCGCGCCATTGCCATAAGCGGAAAGACCACAGTa                                                >  2:1474515/1‑139 (MQ=255)
                                                                      cAGCATCAGTGGCAAGGCTGGGGCAAACAGCCCGAAGCGGCGAAGTTTAACGTTAGCGCAAATTAAAGGCGCTACGCGTAATGCGATGGTGCGACCGCGTTTTGTCGCGCCATTGCCATAAGCGGAAAGACCACAGTaa                                               >  1:1610564/1‑139 (MQ=255)
                                                                               tGGCAAGGCTGGGGCAAACAGCCCGAAGCGGCGAAGTTTAACGTTAGCGCAAATTAAAGGCGCTACGCGTAATGCGATGGTGCGACCGCGTTTTGTcgcg                                                                             >  2:304814/1‑100 (MQ=255)
                                                                               tGGCAAGGCTGGGGCAAACAGCCCGAAGCGGCGAAGTTTAACGTTAGCGCAAATTAAAGGCGCTACGCGTAATGCGATGGTGCGACCGCGTTTTGTcgcg                                                                             <  1:304814/100‑1 (MQ=255)
                                                                                ggCAAGGCTGGGGCAAACAGCCCGAAGCGGCGAAGTTTAACGTTAGCGCAAATTAAAGGCGCTACGCGTAATGCGATGGTGCGACCGCGTTTTGTCGCGCCATTGCCATAAGCGGAAAGACCACAGTAACGCCTGGTaa                                     >  1:1192263/1‑139 (MQ=255)
                                                                                  cAAGGCTGGGGCAACCAGCTCGATTCCGCTCAGTTTCACGTTAGCGCAAATTAACGGCGCTACGCGTAATGCGCTGGTGCGACCGCGTTTTGTCGCGCCATTGCCACAAGCGGAAAGACCACAGTAACGCCTGGTAAcc                                   <  1:208097/139‑1 (MQ=255)
                                                                                  cAAGGCTGGGGCAAACAGCCCGAAGCGGCGAAGTTTAACGTTAGCGCAAATTAAAGGCGCTACGCGTAATGCGATGGTGCGACCGCGTTTTGTCGCGCCATTGCCATAAGCGGAAAGACCACAGTAACGCCTGGTAAcc                                   >  2:1057674/1‑139 (MQ=255)
                                                                                  cAAGGCTGGGGCAAACAGCCCGAAGCGGCGAAGTTTAACGTTAGCGCAAATTAAAGGCGCTACGCGTAATGCGATGGTGCGACCGCGTTTTGTCGCGCCATTGCCATAAGCGGAAAGACCACAGTAACGCCTGGTAAcc                                   <  1:675959/139‑1 (MQ=255)
                                                                                    aGGCTGGGGCAAACAGCCCGAAGCGGCGAAGTTTAACGTTAGCGCaa                                                                                                                             <  1:971373/47‑1 (MQ=255)
                                                                                    aGGCTGGGGCAAACAGCCCGAAGCGGCGAAGTTTAACGTTAGCGCaa                                                                                                                             >  2:971373/1‑47 (MQ=255)
                                                                                         ggggCAAACAGCCCGAAGCGGCGAAGTTTAACGTTAGCGCAAATTAAAGGCGCTACGCGTAATGCGATGGTGCGACCGCGTTTTGTCGCGCCATTGCCATAAGCGGAAAGACCACAGTAACGcc                                           >  1:357082/1‑124 (MQ=255)
                                                                                         ggggCAAACAGCCCGAAGCGGCGAAGTTTAACGTTAGCGCAAATTAAAGGCGCTACGCGTAATGCGATGGTGCGACCGCGTTTTGTCGCGCCATTGCCATAAGCGGAAAGACCACAGTAACGcc                                           <  2:357082/124‑1 (MQ=255)
                                                                                         ggggCAAACAGCCCGAAGCGGCGAAGTTTAACGTTAGCGCAAATTAAAGGCGCTACGCGTAATGCGATGGTGCGACCGCGTTTTGTCGCGCCATTGCCATAAGCGGAAAGACCACAGTAACGCCTGGTAAcc                                   <  1:1571479/132‑1 (MQ=255)
                                                                                         ggggCAAACAGCCCGAAGCGGCGAAGTTTAACGTTAGCGCAAATTAAAGGCGCTACGCGTAATGCGATGGTGCGACCGCGTTTTGTCGCGCCATTGCCATAAGCGGAAAGACCACAGTAACGCCTGGTAAcc                                   >  2:1571479/1‑132 (MQ=255)
                                                                                          gggCAAACAGCCCGAAGCGGCGAAGTTTAACGTTAGCGCAAATTAAAGGCGCTACGCGTAATGCGATGGTGCGACCGCGTTTTGTCGCGCCATTGCCATAAGCGGAAAGACCACAGTaa                                               >  2:1156122/1‑119 (MQ=255)
                                                                                          gggCAAACAGCCCGAAGCGGCGAAGTTTAACGTTAGCGCAAACTAAAGGCGCTACGCGTAATGCGATGGTGCGACCGCGTTTTGTCGCGCCATTGCCATAAGCGGAAAGACCACAGTaa                                               <  1:1156122/119‑1 (MQ=255)
                                                                                           ggCAAACAGCCCGAAGCGGCGAAGTTTAACGTTAGCGCAAATTAAAGGCGCTAcgcg                                                                                                            <  2:1394842/57‑1 (MQ=255)
                                                                                           ggCAAACAGCCCGAAGCGGCGAAGTTTAACGTTAGCGCAAATTAAAGGCGCTAcgcg                                                                                                            >  1:1394842/1‑57 (MQ=255)
                                                                                           ggCAAACAGCCCGAAGCGGCGAAGTTTAACGTTAGCGCAAATTAAAGGCGCTACGCGTAATGCGa                                                                                                    >  2:711855/1‑65 (MQ=255)
                                                                                           ggCAAACAGCCCGAAGCGGCGAAGTTTAACGTTAGCGCAAATTAAAGGCGCTACGCGTAATGCGa                                                                                                    <  1:711855/65‑1 (MQ=255)
                                                                                           ggCAAACAGCCCGAAGCGGCGAAGTTTAACGTTAGCGCAAATTAAAGGCGCTACGCGTAATGCGATGGTGCGACCGCGTTTTGTCGCGCCATTGCCATAAGCGGAAAGACCACAGTAACGCCTGGTAACCGACCTTCAc                          >  1:1552950/1‑139 (MQ=255)
                                                                                                aCAGCCCGAAGCGGCGAAGTTTAACGTTAGCGCAAATTAAAGGCGCTACGCGTAATGCGATGGTGCGACCGCGTTTTGTCGCGCCATTGCCATAAGCGGAAAGACCACAGTAACGCCTGGTAACCGACCTTCACGCTgc                     >  1:342698/1‑139 (MQ=255)
                                                                                                   gCCCGAAGCGGCGAAGTTTAACGTTAGCGCAAATTAAAGGCGCTACGCGTAATGCGATGGTGCGACCGCGTTTTGTCGCGCCATTGCCATAAGCGGAAAGACCACAGTAACGCCTGGTAACCGACCTTCACGCTGCGca                  <  1:1225308/139‑1 (MQ=255)
                                                                                                    cccGAAGCGGCGAAGTTTAACGTTAGCGCAAATTAAAGGCGCTACGCGTAATGCGATGGTGCGACCGCGTTTTGTcg                                                                               <  1:1464857/77‑1 (MQ=255)
                                                                                                    cccGAAGCGGCGAAGTTTAACGTTAGCGCAAATTAAAGGCGCTACGCGTAATGCGATGGTGCGACCGCGTTTTGTcg                                                                               >  2:1464857/1‑77 (MQ=255)
                                                                                                     ccGAAGCGGCGAAGTTTAACGTTAGCGCAAATTAAAGGCGCTACGCGTAATGCGATGGTGCGACCGCGTTTTGTCGCGCCATTGCCATAAGCGGAAAGACCACAGTAACGCCTGGTAACCGACCTTCACGCTGCGcaca                <  1:1198092/139‑1 (MQ=255)
                                                                                                      cGAAGCGGCGACGTTTAACGTTAGCGCAAATTAAAGGCGCTACGCGTAATGCGATGGTGCGACCGCGTTTTGTCGCGCCATTGCCATAAGCGGAAAGACCACAGTAACGCCTGGTAACCGACCTTCACGCTGCGCACAt               <  2:1046198/139‑1 (MQ=255)
                                                                                                      cGAAGCGGCGAACTTTAACGTTAGCGCAAATTAAAGGCGCTACGCGTAATGCGATGGTGCGACCGCGTTTTGTCGCGCCATTGCCATAAGCGGAAAGACCACAGTAGCGCCTGGTAACCGACCTTCACGCTGCGCACAt               <  1:959462/139‑1 (MQ=255)
                                                                                                         aGCGGCGAAGTTTAACGTTAGCGCAAATTAAAGGCGCTACGCGTAATGCGATGGTGCGACCGCGTTTTGTCGCGCCATTGCCATAAGCGGAAAGACCACAGTAACGCCTGGTAACCGACCTTCACg                         >  1:167532/1‑126 (MQ=255)
                                                                                                         aGCGGCGAAGTTTAACGTTAGCGCAAATTAAAGGCGCTACGCGTAATGCGATGGTGCGACCGCGTTTTGTCGCGCCATTGCCATAAGCGGAAAGACCACAGTAACGCCTGGTAACCGACCTTCACg                         <  2:167532/126‑1 (MQ=255)
                                                                                                             gcgAAGTTTAACGTTAGCGCAAATTAAAGGCGCTACGCGTAATGCGATGGTGCGACCGCGTTTTGTCGCGCCATTGCCATAAGCGGAAAGACCACAGTAACGCCTGGTAACCGACCTTCACGCTGCGCACATTGGTAAt        <  1:1057674/139‑1 (MQ=255)
                                                                                                                 aGTTTAACGTTAGCGCAAATTAAAGGCGCTACGCGTAATGCGATGGTGCGACCGCGTTTTGTCGCGCCATTGCCATAAGCGGAAAGACCACAGTAACGCCTGGTAACCGACCTTCACGCTGCGCACATTGGTAATCAgg    <  2:1552950/139‑1 (MQ=255)
                                                                                                                 aGTTTAACGGTAGCGGAAAATAAAGGCGCTACGCGTAATGGGATTGTGCCACCGGGTTTTGTCGCCCCATTTGCATAAGCGGAAAGACCACAGTAACGCCTGGTAACCGACCTTCACGCTGCGCACATTGGTAATCAgg    <  1:355383/139‑1 (MQ=255)
                                                                                                                  gTTTAACGTTAGCGCAAATTAAAGGCGCTACGCGTAATGCGATGGGGCGA‑CGCGTTTTGTCGCGCCATTGCCATAAGCGGAAAGACCACAGTAACGCCTGGTAACCGACCTTCACGCTGCGCACATTGGTAATCAGgcg  <  1:944737/139‑1 (MQ=255)
                                                                                                                          |                                                                                                                                   
TGTTTATGCTAATCGCTTCTTTGAAGTATTAAGCGAAAGCATCAGCCAGAGCACGCTGATCGTTTACTTCCAGCATCAGTGGCAAGGCTGGGGCAAACAGCCCGAAGCGGCGAAGTTTAACGCTAGCGCAAATTAAAGGCGCTACGCGTAATGCGATGGTGCGACCGCGTTTTGTCGCGCCATTGCCATAAGCGGAAAGACCACAGTAACGCCTGGTAACCGACCTTCACGCTGCGCACATTGGTAATCAGGCG  >  NC_000913/3788925‑3789178

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