Predicted mutation
evidence seq id position mutation annotation gene description
RA NC_000913 4,148,681 C→T G529S (GGC→AGC)  eptC ← LPS heptose I phosphoethanolamine transferase

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009134,148,6810CT100.0% 273.2 / NA 85G529S (GGC→AGC) eptCLPS heptose I phosphoethanolamine transferase
Reads supporting (aligned to +/- strand):  ref base C (0/0);  new base T (44/41);  total (44/41)

CCATAGGGCAGTGTGTCGTAATCGATCAGTGCGTTTTTCTTATACGGGTTACCAATCCAGCGGGTAGTTTCTTTGAACTGCGGATTCACCACTGAACGGGTTGGATCGTAACCGTCGTAAGATAAGCCCGCCAAATCTGACCAGGTGTGGATCAGTTCCGCCAGGCTGTATTTACGATCAACATCCTGCGAGAAATCACGGGGATGAGTCGCTTGCCATTTTTCTGACGTCCACAGCAGGAACGGAATGGTG  >  NC_000913/4148554‑4148805
                                                                                                                               |                                                                                                                            
ccATAGGGCAGTGTGTCGTAATCGATCAGTGCGTTTTTCTTATACGGGTTACCAATCCAGCGGGTAGTTTCTTTGAACTGCGGATTCACCACTGAACGGGTTGGATCGTAACCGTCGTAAGATAAGCTCGCCAAATCTg                                                                                                                   >  1:1506403/1‑139 (MQ=255)
  aTAGGGCAGTGTGTCGTAATCGATCAGTGCGTTTTTCTTATACGGGTTACCAATCCAGCGGGTAGTTTCTTTGAACTGCGGATTCACCACTGAACGGGTTGGATCGTAACCGTCGTAAGATAAGCTCGCCAAATCTGAc                                                                                                                 <  2:838265/139‑1 (MQ=255)
   tAGGGCAGTGTGTCGTAATCGATCAGTGCGTTTTTCTTATACGGGTTACCAATCCAGCGGGTAGTTTCTTTGAACTGCGGATTCACCACTGAACGGGTTGGATCGTAACCGTCGTAAGATAAGCTCGCCAAATCTGAcc                                                                                                                <  1:1341348/139‑1 (MQ=255)
   tAGGGCAGTGTGTCGTAATCGATCAGTGCGTTTTTCTTATACGGGTTACCAATCCAGCGGGTAGTTTCTTTGAACTGCGGATTCACCACTGAACGGGTTGGATCGTAACCGTCGTAAGATAAGCTCGCCAAATCTGAcc                                                                                                                >  1:1490107/1‑139 (MQ=255)
        cAGTGTGTCGTAATCGATCAGTGCGTTTTTCTTATACGGGTTACCAATCCAGCGGGTAGTTTCTTTGAACTGCGGATTCACCACTGAACGGGTTGGATCGTAACCGTCGTAAGATAAGCTCGCCAAATCTGAcc                                                                                                                >  2:587065/1‑134 (MQ=255)
        cAGTGTGTCGTAATCGATCAGTGCGTTTTTCTTATACGGGTTACCAATCCAGCGGGTAGTTTCTTTGAACTGCGGATTCACCACTGAACGGGTTGGATCGTAACCGTCGTAAGATAAGCTCGCCAAATCTGAcc                                                                                                                <  1:587065/134‑1 (MQ=255)
           tgtgtCGTAATCGATCAGTGCGTTTTTCTTATACGGGTTACCAATCCAGCGGGTAGTTTCTTTGAACTGCGGATTCACCACTGAACGGGTTGGATCGTAACCGTCGTAAGATAAGCTCGCCAAATCTGACCAGGTGTgg                                                                                                        >  2:816855/1‑139 (MQ=255)
           tgtgtCGTAATCGATCAGTGCGTTTTTCTTATACGGGTTACCAATCCAGCGGGTAGTTTCTTTGAACTGCGGATTCACCACTGAACGGGTTGGATCGTAACCGTCGTAAGATAAGCTCGCCAAATCTGACCAGGTGTgg                                                                                                        >  2:777857/1‑139 (MQ=255)
           tgtgtCGTAATCGATCAGTGCGTTTTTCTTATACGGGTTACCAATCCAGCGGGTAGTTTCTTTGAACTGCGGATTCACCACTGAACGGGTTGGATCGTAACCGTCGTAAGATAAGCTCGCCAAATCTGACCAGGTGTgg                                                                                                        >  1:1237848/1‑139 (MQ=255)
           tgtgtCGTAATCGATCAGTGCGTTTTTCTTATACGGGTTACCAATCCAGCGGGTAGTTTCTTTGAACTGCGGATTCACCACTGAACGGGTTGGATCGTAACCGTCGTAAGATAAGCTCGCCAAATCTGACCAGGTGTgg                                                                                                        <  2:1490107/139‑1 (MQ=255)
             tgtCGTAATCGATCAGTGCGTTTTTCTTATACGGGTTACCAATCCAGCGGGTAGTTTCTTTGAACTGCGGATTCACCACTGAACGGGTTGGATCGTAACCGTCGTAAGATAAGCTCGCCAAATCTGACCAGGTGTGGAt                                                                                                      <  1:816855/139‑1 (MQ=255)
               tCGTAATCGATCAGTGCGTTTTTCTTATACGGGTTACCAATCCAGCGGGTAGTTTCTTTGAACTGCGGATTCACCACTGAACGGGTTGGATCGTAACCGTCGTAAGATAAGCTCGCCAAATCTGACCAGGTGTGGATCa                                                                                                    >  2:1474008/1‑139 (MQ=255)
                        aTCAGTGCGTTTTTCTTATACGGGTTACCAATCCAGCGGGTAGTTTCTTTGAACTGCGGATTCACCACTGAACGGGTTGGATCGTAACCGTCGTAAGATAAGCTCGCCAAATCTGACCAGGTGTGGATCAGTTCCGCCa                                                                                           >  2:1561433/1‑139 (MQ=255)
                        aTCAGTGCGTTTTTCTTATACGGGTTACCAATCCAGCGGGTAGTTTCTTTGAACTGCGGATTCACCACTGAACGGGTTGGATCGTAACCGTCGTAAGATAAGCTCGCCAAATCTGACCAGGTGTGGATCAGTTCCGCCa                                                                                           >  2:1494599/1‑139 (MQ=255)
                          cAGTGCGTTTTTCTTATACGGGTTACCAATCCAGCGGGTAGTTTCTTTGAACTGCGGATTCACCACTGAACGGGTTGGATCGTAACCGTCGTAAGATAAGCTCGCCAAATCTGACCAGGTGTGGATCAGTTCCGCCAgg                                                                                         <  2:1636613/139‑1 (MQ=255)
                           aGTGCGTTTTTCTTATACGGGTTACCAATCCAGCGGGTAGTTTCTTTGAACTGCGGATTCACCACTGAACGGGTTGGATCGTAACCGTCGTAAGATAAGCTCGCCAAATCTGACCAGGTGTGGATCAGTTCCGCCAGGc                                                                                        <  1:1561433/139‑1 (MQ=255)
                             tGCGTTTTTCTTATACGGGTTACCAATCCAGCGGGTAGTTTCTTTGAACTGCGGATTCACCACTGAACGGGTTGGATCGTAACCGTCGTAAGATAAGCTCGCCAAATCTGACCAGGTGTGGATCAGTTCCGCCAGGCTg                                                                                      >  2:990669/1‑139 (MQ=255)
                             tGCGTTTTTCTTATACGGGTTACCAATCCAGCGGGTAGTTTCTTTGAACTGCGGATTCACCACTGAACGGGTTGGATCGTAACCGTCGTAAGATAAGCTCGCCAAATCTGACCAGGTGTGGATCAGTTCCGCCAGGCTg                                                                                      <  1:990669/139‑1 (MQ=255)
                                   tttCTTATACGGGTTACCAATCCAGCGGGTAGTTTCTTTGAACTGCGGATTCACCACTGAACGGGTTGGATCGTAACCGTCGTAAGATAAGCTCGCCAAATCTGACCAGGTGTGGATCAGTTCCGCCAGGCTGTATTTa                                                                                <  1:1494599/139‑1 (MQ=255)
                                    ttCTTATACGGGTTACCAATCCAGCGGGTAGTTTCTTTGAACTGCGGATTCACCACTGAACGGGTTGGATCGTAACCGTCGTAAGATAAGCTCGCCAAATCTGACCAGGTGTGGATCAGTTCCGCCAGGCTGTATTTAc                                                                               <  2:1237848/139‑1 (MQ=255)
                                       ttATACGGGTTACCAATCCAGCGGGTGGTTTCTTTGAACTGCGGATTCACCACTGAACGGGTTGGATCGTAACCGTCGTAAGATAAGCTCGCCAAATCTGACCAGGTGTGGATCAGTTCCGCCAGGCTGTAtt                                                                                  >  1:648443/1‑133 (MQ=255)
                                       ttATACGGGTTACCAATCCAGCGGGTAGTTTCTTTGAACTGCGGATTCACCACTGAACGGGTTGGATCGTAACCGTCGTAAGATAAGCTCGCCAAATCTGACCAGGTGTGGATCAGTTCCGCCAGGCTGTAtt                                                                                  <  2:648443/133‑1 (MQ=255)
                                       ttATACGGGTTACCAATCCAGCGGGTAGTTTCTTTGAACTGCGGATTCACCACTGAACGGGTTGGATCGTAACCGTCGTAAGATAAGCTCGCCAAATCTGACCAGGTGTGGATCAGTTCCGCCAGGCTGTATTTACGAt                                                                            <  2:1506403/139‑1 (MQ=255)
                                       ttATACGGGTTACCAATCCAGCGGGTAGTTTCTTTGAACTGCGGATTCACCACTGAACGGGTTGGATCGTAACCGTCGTAAGATAAGCTCGCCAAATCTGACCAGGTGTGGATCAGTTCCGCCAGGCTGTATTTACGAt                                                                            >  2:891690/1‑139 (MQ=255)
                                         ataCGGGTTACCAATCCAGCGGGTAGTTTCTTTGAACTGCGGATTCACCACTGAACGGGTTGGATCGTAACCGTCGTAAGATAAGCTCGCCAAATCTGAc                                                                                                                 >  2:498161/1‑100 (MQ=255)
                                         ataCGGGTTACCAATCCAGCGGGTAGTTTCTTTGAACTGCGGATTCACCACTGAACGGGTTGGATCGTAACCGTCGTAAGATAAGCTCGCCAAATCTGAc                                                                                                                 <  1:498161/100‑1 (MQ=255)
                                               gTTACCAATCCAGCGGGTAGTTTCTTTGAACTGCGGATTCACCACTGAACGGGTTGGATCGTAACCGTCGTAAGATAAGCTCGCCAAATCTGACCAGGTGTGGATCAg                                                                                                   <  1:1627837/108‑1 (MQ=255)
                                               gTTACCAATCCAGCGGGTAGTTTCTTTGAACTGCGGATTCACCACTGAACGGGTTGGATCGTAACCGTCGTAAGATAAGCTCGCCAAATCTGACCAGGTGTGGATCAg                                                                                                   >  2:1627837/1‑108 (MQ=255)
                                               gTTACCAATCCAGCGGGTAGTTTCTTTGAACTGCGGATTCACCACTGAACGGGTTGGATCGTAACCGTCGTAAGATAAGCTCGCCAAATCTGACCAGGTGTGGATCAGTTCCGCCAGGCTGTAttt                                                                                 <  1:1264539/126‑1 (MQ=255)
                                               gTTACCAATCCAGCGGGTAGTTTCTTTGAACTGCGGATTCACCACTGAACGGGTTGGATCGTAACCGTCGTAAGATAAGCTCGCCAAATCTGACCAGGTGTGGATCAGTTCCGCCAGGCTGTAttt                                                                                 >  2:1264539/1‑126 (MQ=255)
                                                     aaTCCAGCGGGTAGTTTCTTTGAACTGCGGATTCTCCACTGAACGGGTTGGATCGTAACCGTCGTAAGATAAGCTCGCCAAATCTGACCAGGTGTGGATCAGTTCCGCCAGGCTGTATTTACGATCAACATCCTGCgag                                                              >  1:1370836/1‑139 (MQ=255)
                                                       tCCAGCGGGTAGTTTCTTTGAACTGCGGATTCACCACTGAACGGGTTGGATCGTAACCGTCGTAAGATAAGCTCGCCAAATCTGAcc                                                                                                                <  2:150978/87‑1 (MQ=255)
                                                       tCCAGCGGGTAGTTTCTTTGAACTGCGGATTCACCACTGAACGGGTTGGATCGTAACCGTCGTAAGATAAGCTCGCCAAATCTGAcc                                                                                                                >  1:150978/1‑87 (MQ=255)
                                                         cAGCGGGTAGTTTCTTTGAACTGCGGATTCACCACTGAACGGGTTGGATCGTAACCGTCGTAAGATAAGCTCGCCAAATCTGACCAGGTGTGGATCAGTTCCGCCAGGCTGTATTTACGATCAACATCCTGCGAGAAAt                                                          >  2:1155888/1‑139 (MQ=255)
                                                         cAGCGGGTAGTTTCTTTGAACTGCGGATTCACCACTGAACGGGTTGGATCGTAACCGTCGTAAGATAAGCTCGCCAAATCTGACCAGGTGTGGATCAGTTCCGCCAGGCTGTATTTACGATCAACATCCTGCGAGAAAt                                                          >  1:782191/1‑139 (MQ=255)
                                                           gCGGGTAGTTTCTTTGAACTGCGGATTCACCACTGAACGGGTTGGATCGTAACCGTCGTAAGATAAGCTCGCCAAATCTGAcc                                                                                                                >  1:848437/1‑83 (MQ=255)
                                                           gCGGGTAGTTTCTTTGAACTGCGGATTCACCACTGAACGGGTTGGATCGTAACCGTCGTAAGATAAGCTCGCCAAATCTGAcc                                                                                                                <  2:848437/83‑1 (MQ=255)
                                                            cGGGTAGTTTCTTTGAACTGCGGATTCACCACTGAACGGGTTGGATCGTAACCGTCGTAAGATAAGCTCGCCAAATCTGACCAGGTGTGGATCAGTTCCGCCAGGCTGTATTTAc                                                                               >  1:1227101/1‑115 (MQ=255)
                                                            cGGGTAGTTTCTTTGAACTGCGGATTCACCACTGAACGGGTTGGATCGTAACCGTCGTAAGATAAGCTCGCCAAATCTGACCAGGTGTGGATCAGTTCCGCCAGGCTGTATTTAc                                                                               <  2:1227101/115‑1 (MQ=255)
                                                             gggTAGTTTCTTTGAACTGCGGATTCACCACTGAACGGGTTGGATCGTAACCGTCGTAAGATAAGCTCGcc                                                                                                                          >  2:190908/1‑71 (MQ=255)
                                                             gggTAGTTTCTTTGAACTGCGGATTCACCACTGAACGGGTTGGATCGTAACCGTCGTAAGATAAGCTCGcc                                                                                                                          <  1:190908/71‑1 (MQ=255)
                                                              ggTAGTTTCTTTGAACTGCGGATTCACCACTGAACGGGTTGGATCGTAACCGTCGTAAGATAAGCTCGCCAAATCTGACCAGGTGTGGATCAGTTCCGCCAGGc                                                                                        >  2:349736/1‑104 (MQ=255)
                                                              ggTAGTTTCTTTGAACTGCGGATTCACCACTGAACGGGTTGGATCGTAACCGTCGTAAGATAAGCTCGCCAAATCTGACCAGGTGTGGATCAGTTCCGCCAGGc                                                                                        <  1:349736/104‑1 (MQ=255)
                                                              ggTAGTTTCTTTGAACTGCGGATTCACCACTGAACGGGTTGGATCGTAACCGTCGTAAGATAAGCTCGCCAAATCTGACCAGGTGTGGATCAGTTCCGCCAGGCTGTATTTACGATCAACATCCTGCGAGAAATCACgg                                                     <  1:777857/139‑1 (MQ=255)
                                                                  gTTTCTTTGAACTGCGGATTCACCACTGAACGGGTTGGATCGTAACCGTCGTAAGATAAGCTCGCCAAATCTGACCAGGTGTGGATCAGTTCCGCCAGGCTGTATTTACGATCAACATCCTGCGAGaaa                                                           >  1:1212948/1‑129 (MQ=255)
                                                                  gTTTCTTTGAACTGCGGATTCACCACTGAACGGGTTGGATCGTAACCGTCGTAAGATAAGCTCGCCAAATCTGACCAGGTGTGGATCAGTTCCGCCAGGCTGTATTTACGATCAACATCCTGCGAGaaa                                                           <  2:1212948/129‑1 (MQ=255)
                                                                        ttGAACTGCGGATTCACCACTGAACGGGTTGGATCGTAACCGTCGTAAGATAAGCTCGCCAAATCTGACCAGGTGTGGATCAGTTCCGCCAGGCTGTATTTAc                                                                               <  1:1131757/103‑1 (MQ=255)
                                                                        ttGAACTGCGGATTCACCACTGAACGGGTTGGATCGTAACCGTCGTAAGATAAGCTCGCCAAATCTGACCAGGTGTGGATCAGTTCCGCCAGGCTGTATTTAc                                                                               >  2:1131757/1‑103 (MQ=255)
                                                                             cTGCGGATTCACCACTGAGCGGGTTGGATCGTAACCGTCGTAAGATAAGTTCGCCAAATCTGACCAGGTGTGGATCAGTTCCGCCAGGCTGTAttt                                                                                 <  2:1029034/96‑1 (MQ=255)
                                                                             cTGCGGATTCACCACTGAACGGGTTGGATCGTAACCGTCGTAAGATAAGTTAGCCAAATCTGACCAGGTGTGGATCAGTTCCGCCAGGCTGTAttt                                                                                 >  1:1029034/1‑96 (MQ=255)
                                                                                cGGATTCACCACTGAACGGGTTGGATCGTAACCGTCGTAAGATAAGCTCGCCAAATCTGACCAGgt                                                                                                            >  1:914649/1‑66 (MQ=255)
                                                                                cGGATTCACCACTGAACGGGTTGGATCGTAACCGTCGTAAGATAAGCTCGCCAAATCTGACCAGgt                                                                                                            <  2:914649/66‑1 (MQ=255)
                                                                                cGGATTCACCACTGAACGGGTTGGATCGTAACCGTCGTAAGATAAGCTCGCCAAATCTGACCAGGTGTGGATCAGtt                                                                                                 >  2:561593/1‑77 (MQ=255)
                                                                                cGGATTCACCACTGAACGGGTTGGATCGTAACCGTCGTAAGATAAGCTCGCCAAATCTGACCAGGTGTGGATCAGtt                                                                                                 <  1:561593/77‑1 (MQ=255)
                                                                                cGGATTCACCACTGAACGGGTTGGATCGTAACCGTCGTAAGATAAGCTCGCCAAATCTGACCAGGTGTGGATCAGTTCCGcc                                                                                            >  2:677647/1‑82 (MQ=255)
                                                                                cGGATTCACCACTGAACGGGTTGGATCGTAACCGTCGTAAGATAAGCTCGCCAAATCTGACCAGGTGTGGATCAGTTCCGcc                                                                                            <  1:677647/82‑1 (MQ=255)
                                                                                cGGATTCACCACTGAACGGGTTGGATCGTAACCGTCGTAAGATAAGCTCGCCAAATCTGACCAGGTGTGGATCAGTTCCGCCAGGCTGTATTTACGATCAACATCCTGCGAGAAATCACGGGGATGAGTCGCTTGCCAt                                   >  1:1559097/1‑139 (MQ=255)
                                                                                cGGATTCACCACTGAACGGGTTGGATCGTAACCGTCGTAAGATAAGCTCGCCAAATCTGACCAGGTGTGGATCAGTTCCGCCAGGCTGTATTTACGATCAACATCCTGCGAGAAATCACGGGGATGAGTCGCTTGCCAt                                   >  2:1663160/1‑139 (MQ=255)
                                                                                cGGATTCACCACTGAACGGGTTGGATCGTAACCGTCGTAAGATAAGCTCGCCAAATCTGACCAGGTGTGGATCAGTTCCGCCAGGCTGTATTTACGATCAACATCCTGCGAGAAATCACGGGGATGAGTCGCTTGCCAt                                   >  2:200323/1‑139 (MQ=255)
                                                                                cGGATTCACCACTGAACGGGTTGGATCGTAACCGTCGTAAGATAAGCTCGCCAAATCTGACCAGGTGTGGATCAGTTCCGCCAGGCTGTATTTACGATCAACATCCTGCGAGAAATCACGGGGATGAGTCGCTAGCCAt                                   >  1:1542257/1‑139 (MQ=255)
                                                                                     tCACCACTGAACGGGTTGGATCGTAACCGTCGTAAGATAAGCTCGCCAAATCTGAcc                                                                                                                <  2:799963/57‑1 (MQ=255)
                                                                                     tCACCACTGAACGGGTTGGATCGTAACCGTCGTAAGATAAGCTCGCCAAATCTGAcc                                                                                                                >  1:799963/1‑57 (MQ=255)
                                                                                     tCACCACTGAACGGGTTGGATCGTAACCGTCGTAAGATAAGCTCGCCAAATCTGACCAGGTGTGGATCAGTTCCGCCAGGCTGTATTTACGATCAACATCCTGCGAGAAATCACGGGGATGAGTCGCTTGCCATTTTTc                              >  2:1198011/1‑139 (MQ=255)
                                                                                            tGAACGGGTTGGATCGTAACCGTCGTAAGATAAGCTCGCCAAATCTGACCAGGTGTGGATCAGTTCCGCCAGGCTGTATTTACGATCAACATCCTGCGAGAAATCACGGGGATGAGTCGCTTGCCATTTTTCTGACGTc                       <  1:1474008/139‑1 (MQ=255)
                                                                                              aaCGGGTTGGATCGTAACCGTCGTAAGATAAGCTCGCCAAATCTGACCAGGTGTGGATCAGTTCCGCCAGGCTGTATTTACGATCAACATCCTGCGAGAAATCACGGGGATGAGTCGCTTGCCATTTTTCTGACGTCca                     >  1:179467/1‑139 (MQ=255)
                                                                                                cGGGTTGGATCGTAACCGTCGTAAGATAAGCTCGCCAAATCTGACCAGGTGTGGATCAGTTCCGCCAGGCTGTATTTACGATCAACATCCTGCGAGAAATCACGGGGATGAGTCGCTTGCCATTTTTCTGACGTCcaca                   >  2:396674/1‑139 (MQ=255)
                                                                                                  ggTTGGATCGTAACCGTCGTAAGATAAGCTCGCCAAATCTGACCAGGTGTGGATCAGTTCCGCCAGGCTGTATTTACGATCAACATCCTGCGAGAAATCACGGGGATGAGTCGCTTGCCATTTTTCTGACGTCCAcagc                 <  2:782191/139‑1 (MQ=255)
                                                                                                   gTTGGATCGTAACCGTCGTAAGATAAGCTCGCCAAATCTGACCAGGTGTGGATCAGTTCCGCCAGGCTGTATTTACGATCAACATCCTGCGAGAAATCACGGGGATGAGTCGCTTGCCATTTTTCTGACGTCCAcagca                >  2:46124/1‑139 (MQ=255)
                                                                                                            tAACCGTCGTAAGATAAGCTCGCCAAATCTGACCAGGTGTGGATCAGTTCCGCCAGGCTGTATTTACGATCAACATCCTGCGAGAAATCACGGGGATGAGTCGCTTGCCATTTTTCTGACGTCCACAGCAGGAACGGaa       <  2:1370836/139‑1 (MQ=255)
                                                                                                             aaCCGTCGTAAGATAAGCTCGCCAAATCTGACCAGGTGTGGATCAGTTCCGCCAGGCTGTATTTACGATCAACATCCTGCGAGAAATCACGGGGATGAGTCGCTTGCCATTTTTCTGACGTCCACAGCAGGAACGGAAt      >  2:1060924/1‑139 (MQ=255)
                                                                                                             aaCCGTCGTAAGATAAGCTCGCCAAATCTGACCAGGTGTGGATCAGTTCCGCCAGGCTGTATTTACGATCAACATCCTGCGAGAAATCACGGGGATGAGTCGCTTGCCATTTTTCTGACGTCCACAGCAGGAACGGAAt      <  2:1559097/139‑1 (MQ=255)
                                                                                                             aaCCGTCGTAAGATAAGCTCGCCAAATCTGACCAGGTGTGGATCAGTTCCGCCAGGCTGTATTTACGATCAACATCCTGCGAGAAATCACGGGGATGAGTCGCTTGCCATTTTTCTGACGTCCACAGCAGGAACGGAAt      <  2:1542257/139‑1 (MQ=255)
                                                                                                                 gtcgtAAGATAAGCTCGCCAAATCTGACCAGGTGTGGATCAGTTCCGCCAGGCTGTATTTACGATCAACATCCTGCGAGAAATCACGGGGATGAGTCGCTTGCCATTTTTCTGACGTCCACAGCAGGAACGGAATGgtg  <  1:1060924/139‑1 (MQ=255)
                                                                                                                 gtcgtAAGATAAGCTCGCCAAATCTGACCAGGTGTGGATCAGTTCCGCCAGGCTGTATTTACGATCAACATCCTGCGAGAAATCACGGGGATGAGTCGCTTGCCATTTTTCTGACGTCCACAGCAGGAACGGAATGgtg  <  1:396674/139‑1 (MQ=255)
                                                                                                                  tcgtAAGATAAGCTCGCCAAATCTGACCAGGTGTGGATCAGTTCCGCCAGGc                                                                                        >  2:1280859/1‑52 (MQ=255)
                                                                                                                  tcgtAAGATAAGCTCGCCAAATCTGACCAGGTGTGGATCAGTTCCGCCAGGc                                                                                        <  1:1280859/52‑1 (MQ=255)
                                                                                                                  tcgtAAGATAAGCTCGCCAAATCTGACCAGGTGTGGATCAGTTCCGCCAGGCTGTAttt                                                                                 >  2:1075032/1‑59 (MQ=255)
                                                                                                                  tcgtAAGATAAGCTCGCCAAATCTGACCAGGTGTGGATCAGTTCCGCCAGGCTGTAttt                                                                                 <  1:1075032/59‑1 (MQ=255)
                                                                                                                  tcgtAAGATAAGCTCGCCAAATCTGACCAGGTGTGGATCAGTTCCGCCAGGCTGTATTTACGATCAACATCCTGCGAGAAATCACGGGGATGAGTCGCTTGCCAttttt                               >  1:94401/1‑109 (MQ=255)
                                                                                                                  tcgtAAGATAAGCTCGCCAAATCTGACCAGGTGTGGATCAGTTCCGCCAGGCTGTATTTACGATCAACATCCTGCGAGAAATCACGGGGATGAGTCGCTTGCCAttttt                               <  2:94401/109‑1 (MQ=255)
                                                                                                                   cgtAAGATAAGCTCGCCAAATCTGACCAGGTGTGGATCAGTTCCGCCAGGCt                                                                                       >  2:1557086/1‑52 (MQ=255)
                                                                                                                   cgtAAGATAAGCTCGCCAAATCTGACCAGGTGTGGATCAGTTCCGCCAGGCt                                                                                       <  1:1557086/52‑1 (MQ=255)
                                                                                                                            aGCTCGCCAAATCTGACCAGGTGTGGATCAGTTCCGCCAg                                                                                          <  1:1509616/40‑1 (MQ=255)
                                                                                                                            aGCTCGCCAAATCTGACCAGGTGTGGATCAGTTCCGCCAGGCTGTATTTACGATCAACATc                                                                     <  2:414100/61‑1 (MQ=255)
                                                                                                                            aGCTCGCCAAATCTGACCAGGTGTGGATCAGTTCCGCCAGGCTGTATTTACGATCAACATc                                                                     >  1:414100/1‑61 (MQ=255)
                                                                                                                               |                                                                                                                            
CCATAGGGCAGTGTGTCGTAATCGATCAGTGCGTTTTTCTTATACGGGTTACCAATCCAGCGGGTAGTTTCTTTGAACTGCGGATTCACCACTGAACGGGTTGGATCGTAACCGTCGTAAGATAAGCCCGCCAAATCTGACCAGGTGTGGATCAGTTCCGCCAGGCTGTATTTACGATCAACATCCTGCGAGAAATCACGGGGATGAGTCGCTTGCCATTTTTCTGACGTCCACAGCAGGAACGGAATGGTG  >  NC_000913/4148554‑4148805

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