Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 3,705,984 T→G 33.7% intergenic (‑194/+114) dppB ← / ← dppA dipeptide/heme ABC transporter permease/dipeptide/heme ABC transporter periplasmic binding protein; dipeptide chemotaxis receptor

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009133,705,9840TG33.7% 492.6 / inf 291intergenic (‑194/+114)dppB/dppAdipeptide/heme ABC transporter permease/dipeptide/heme ABC transporter periplasmic binding protein; dipeptide chemotaxis receptor
Reads supporting (aligned to +/- strand):  ref base T (82/111);  new base G (92/6);  total (174/117)
Fisher's exact test for biased strand distribution p-value = 2.63e-19
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 3.34e-02

Only 100 of 295 total aligned reads displayed.

TTGCTGCCGGATGGCGCTGCACGCAGAGGCACAGCGCCACCGGGTATACAGATTGTTGACAAAGTGCGCTTTGTTTATGCCGGATGCGGCGTTAACGCCTTATCCGGCCT‑‑ACAAAATCGTGCAAATTCAATACATTGCAGAATTTTC‑‑GTAGGCCTGATAAGCGGAGCGCATCAGGCATTTTTGCCTTTGCCATCAGTCT  >  NC_000913/3705885‑3706083
                                                                                                   |                                                                                                           
  gctgcCGGATGGCGCTGCACGCAGAGGCACAGCGCCACCGGGTATACAGATTGTTGACAAAGTGCGCTTTGTTTATGCCGGATGCGGCGTTAACGCCTTa                                                                                                           <  1:1774437/100‑1 (MQ=255)
   ctgcCGGATGGCGCTGCACGCAGAGGCACAGCGCCACCGGGTATACAGATTGTTGACAAAGTGCGCTTTGTTTATGCCGGATGCGGCGTTAACGCCTTAt                                                                                                          >  1:2797547/1‑100 (MQ=255)
    tgcCGGATGGCGCTGCACGCAGAGGCACAGCGCCACCGGGTATACAGATTGTTGACAAAGTGCGCTTTGTTTATGCCGGATGCGGCGTTAACGCCTTATc                                                                                                         >  2:1239567/1‑100 (MQ=255)
     gcCGGATGGCGCTGCACGCAGAGGCACAGCGCCACCGGGTATACAGATTGTTGACAAAGTGCGCTTTGTTTATGCCGGATGCGGCGTTAACGCCTTATcc                                                                                                        >  2:1896680/1‑100 (MQ=255)
      ccGGATGGCGCTGCACGCAGAGGCACAGCGCCACCGGGTATACAGATTGTTGACAAAGTGCGCTTTGTTTATGCCGGATGCGGCGTTAACGCCTTATCCg                                                                                                       >  2:634613/1‑100 (MQ=255)
      caggATGGCGCTGCACGCAGAGGCACAGCGCCACCGGGTATACAGATTGTTGACAAAGTGCGCTTTGTTTATGCCGGATGCGGCGTTAACGCCTTATCCg                                                                                                       <  1:2549909/98‑1 (MQ=255)
      caggATGGCGCTGCACGCAGAGGCACAGCGCCACCGGGTATACAGATTGTTGACAAAGTGCGCTTTGTTTATGCCGGATGCGGCGTTAACGCCTTATCCg                                                                                                       <  1:3440299/98‑1 (MQ=255)
        ggATGGCGCTGCACGCAGAGGCACAGCGCCACCGGGTATACAGATTGTTGACAAAGTGCGCTTTGTTTATGCCGGATAAGGCGTTAACGCCGCATCCGGc                                                                                                     <  1:1981352/100‑1 (MQ=37)
        ggATGGCGCTGCACGCAGAGGCACAGCGCCACCGGGTATACAGATTGTTGACAAAGTGCGCTTTGTTTATGCCGGATAAGGCGTTAACGCCGCATCCGGc                                                                                                     <  2:1005494/100‑1 (MQ=37)
        ggATGGCGCTGCACGCAGACGCACAGCGCCACCGGGTATACAGATTGTTGACAAAGTGCGCTTGGTTTATGCCGTTTGCGGCTTTAACGCCtttttcg                                                                                                       >  2:1585974/1‑93 (MQ=25)
         gATGGCGCTGCACGCAGAGGCACAGCGCCACCGGGTATACAGATTGTTGACAAAGTGCGCTTTGTTTATGCCGGATGCGGCGTTAACGCCTTATCCGGcc                                                                                                    <  1:552597/100‑1 (MQ=255)
         gATGGCGCTGCACGCAGAGGCACAGCGCCACCGGGTATACAGATTGTTGACAAAGTGCGCTTTGTGTATGCCGGATGCGGCGTTAACGCCTTATCCGGcc                                                                                                    <  2:295458/100‑1 (MQ=255)
          aTGGCGCTGCACGCAGAGGCACAGCGCCACCGGGTATACAGATTGTTGACAAAGTGCGCTTTGTTTATGCCGGATGCGGCGTTAACGCCTTATCCGGcc                                                                                                    <  2:3391859/99‑1 (MQ=255)
          aTGGCGCTGCACGCAGAGGCACAGCGCCACCGGGTATACAGATTGTTGACAAAGTGCGCTTTGTTTATGCCGGATGCGGCGTTAACGCCTTATCCGGCCt                                                                                                   <  2:1868115/100‑1 (MQ=255)
           tGGCGCTGCACGCAGAGGCACAGCGCCACCGGGTATACAGATTGTTGACAAAGTGCGCTTTGTTTATGCCGGATGCGGCGTTAACGCCTTATCCGGCCT‑‑a                                                                                                >  1:2932017/1‑100 (MQ=255)
              cgcTGCACGCAGAGGCACAGCGCCACCGGGTATACAGATTGTTGACAAAGTGCGCTTTGTTTATGCCGGATGCGGCGTTAACGCCTTATCCGGCCT‑‑ACaa                                                                                            >  1:3822560/1‑100 (MQ=255)
                cTGCACGCAGAGGCACAGCGCCACCGGGTATACAGATTGTTGACAAAGTGCGCTTTGTTTATACCGGATCCGGCGTTAACGCCTTATCCGCCCT‑‑ACaaa                                                                                           >  2:2724628/1‑99 (MQ=255)
                 tGCACGCAGAGGCACAGCGCCACCGGGTATACAGATTGTTGACAAAGTGCGCTTTGTTTATGCCGGATGCGGCGTTAACGCCTTATCCGGCCT‑‑ACAAAAt                                                                                         >  1:1325645/1‑100 (MQ=255)
                 tGCACGCAGAGGCACAGCGCCACCGGGTATACAGATTGTTGACAAAGTGCGCTTTGTTTATGCCGGATGCGGCGTTAACGCCTTATCCGGCCT‑‑ACAAAAt                                                                                         >  1:3256633/1‑100 (MQ=255)
                  gCACGCAGAGGCACAGCGCCACCGGGTATACAGATTGTTGACAAAGTGCGCTTTGTTTATGCCGGATGCGGCGTTAACGCCTTATCCGGCCT‑‑ACAAAATc                                                                                        >  1:1916829/1‑100 (MQ=255)
                  gCACGCAGAGGCACAGCGCCACCGGGTATACAGATTGTTGACAAAGTGCGCTTTGTTTATGCCGGATGCGGCGTTAACGCCTTATCCGGCCT‑‑ACAAAATc                                                                                        >  1:3011338/1‑100 (MQ=255)
                   gaCGCAGAGGCACAGCGCCACCGGGTATACAGATTGTTGACACAGTACGCTATGTTTATGCCGAATGCGTCTTTAACGCCTTATCCGTCCT‑‑Acgca                                                                                            >  2:1485222/2‑93 (MQ=25)
                   cACGCAGAGGCACAGCGCCACCGGGTATACAGATTGTTGACAAAGTGCGCTTTGTTTATGCCGGATGCGGCGTTAACGCCTTATCCGGCCT‑‑ACAAAATCg                                                                                       <  1:3510278/100‑1 (MQ=255)
                   cACGCAGAGGCACAGCGCCACCGGGTATACAGATTGTTGACAAAGTGCGCTTTGTTTATGCCGGATGCGGCGTTAACGCCTTATCCGGCCT‑‑ACAAAATCg                                                                                       <  2:3814777/100‑1 (MQ=255)
                    aCGCAGAGGCACAGCGCCACCGGGTATACAGATTGTTGACAAAGTGCGCTTTGTTTATGCCGGATGCGGCGTTAACGCCTTATCCGGCCT‑‑ACAAAATCGt                                                                                      <  2:834424/100‑1 (MQ=255)
                     cGCAGAGGCACAGCGCCACCGGGTATACAGATTGTTGACAAAGTGCGCTTTGTTTATGCCGGATGCGGCGTTAACGCCTTATCCGGCCT‑‑ACAAAATCGt                                                                                      >  2:1686935/1‑99 (MQ=255)
                     cGCAGAGGCACAGCGCCACCGGGTATACAGATTGTTGACAAAGTGCGCTTTGTTTATGCCGGATGCGGCGTTAACGCCTTATCCGGCCT‑‑ACAAAATCGTg                                                                                     <  2:814803/100‑1 (MQ=255)
                     cGCAGAGGCACAGCGCCACCGGGTATACAGATTGTTGACAAAGTGCGCTTTGTTTATGCCGGATGCGGCGTTAACGCCTTATCCGGCCT‑‑ACAAAATCGTg                                                                                     <  1:3787153/100‑1 (MQ=255)
                          agGCACAGCGCCACCGGGTATACAGATTGTTGACAAAGTGCGCTTTGTTTATGCCGGATGCGGCGTTAACGCCTTATCCGGCCT‑‑ACAAAATCGTGCAAAt                                                                                <  1:3620351/100‑1 (MQ=255)
                            gCACAGCGCTACCGGGTATTCAGATGGTTGACAAAGTGCGCTTTGTTTATGCCGGATGCGGCGTTAACGCCTTATCCGGCCT‑‑ACAAAATCGTGCAAATTc                                                                              <  2:1924443/100‑1 (MQ=255)
                            gCACAGCGCCACCGGGTATACAGATTGTTGACAAAGTGCGCTTTGTTTATGCCGGATGCGGCGTTAACGCCTTATCCGGCCT‑‑ACAAAATCGTGCAAATTc                                                                              <  1:1653585/100‑1 (MQ=255)
                             cacaGCGCCACCGGGTATACAGATTGTTGACAAAGTGCGCTTTGTTTATGCCGGATGCGGCGTTAACGCCTTATCCGGCCT‑‑ACAAAATCGTGCAAATTCa                                                                             >  1:1399735/1‑100 (MQ=255)
                              acaGCGCCACCGGGTATACAGATTGTTGACAAAGTGCGCTTTGTTTATGCCGGATGCGGCGTTAACGCCTTATCCGGCCT‑‑ACAAAATCGTGCAAATTCaa                                                                            <  1:1145316/100‑1 (MQ=255)
                               caGCGCCACCGGGTATACAGATTGTTGACTAAGTGCGCTTTGTTTATGCCGGATGCGGCGTTAACGTCTTATCCGGCCT‑‑ACAAAATCGTGCAAATTCAAt                                                                           <  2:2364817/100‑1 (MQ=255)
                                aGCGCCACCGGGTATACAGATTGTTGACAAAGTGCGCTTTGTTTATGCCGGATGCGGCGTTAACGCCTTATCCGGCCT‑‑ACAAAATCGTGCAAATTCAATa                                                                          <  2:128994/100‑1 (MQ=255)
                                     cACCGGGTATACAGATTGTTGACAAAGTGCGCTTTGTTTATGCCGGATGCGGCGTTAACGCCTTATCCGGCCT‑‑ACAAAATCGTGCAAATTCAATACATTg                                                                     <  2:1773465/100‑1 (MQ=255)
                                     cACCGGGTATACAAATTGTTGACAAAGTGCGATTTGTTTATGGCGGATGCGGCGTTAACGCCTTATCCGGCCT‑‑ACAAAATCGTGCAAATTCAATACATTg                                                                     >  1:2662340/1‑100 (MQ=255)
                                         gggTATACAGATTGTTGACAAAGTGCGCTTTGTTTATGCCGGATGCGGCGTTAACGCCTTATCCGGCCT‑‑ACAAAATCGTGCAAATTCAATACATTGCAGa                                                               >  1:819259/1‑100 (MQ=255)
                                             ataCAGATTGTTGACAAAGTGCGCTTTGTTTATGCCGGATGCGGCGTTAACGCCTTATCCGGCCT‑‑ACAAAATCGTGCAAATTCAATACATTGCAGAAttt                                                          <  2:433215/100‑1 (MQ=255)
                                              taCAGATTGTTGACAAAGTGCGCTTTGTTTATGCCGGATGCGGCGTTAACGCCTTATCCGGCCT‑‑ACAAAATCGTGCAAATTCAATACATTGCAGAAtttt                                                         >  1:1104764/1‑100 (MQ=255)
                                              taCAGATTGTTGACAAAGTGCGCTTTGTTTATGCCGGATGCGGCGTTAACGCCTTATCCGGCCT‑‑ACAAAATCGTGCAAATTCAATACATTGCAGAAtttt                                                         >  1:163557/1‑100 (MQ=255)
                                                cAGATTGTTGACAAAGTGCGCTTTGTTTATGCCGGATGCGGCGTTAACGCCTTATCCGGCCT‑‑ACAAAATCGTGCAAATTCAATACATTGCAGAATTTTC‑‑g                                                     <  2:2290613/100‑1 (MQ=255)
                                                cAGATTGTTGACAAAGTGCGCTTTGTTTATGCCGGATGCGGCGTTAACGCCTTATCCGGCCT‑‑ACAAAATCGTGCAAATTCAATACATTGCAGAATTTTC‑‑g                                                     <  2:3715525/100‑1 (MQ=255)
                                                    ttgttgACAAAGTGCGCTTTGTTTATGCCGGATGCGGCGTTAACGCCTTATCCGGCCT‑‑ACAAAATCGTGCAAATTCAATACATTGCAGAATTTTC‑‑GTAgg                                                 <  1:698438/100‑1 (MQ=255)
                                                     cataaaCAAAGCGCACTTTGTTTATGCCGGATGCGGCGTTAACGCCTTATCCGGCCT‑‑ACAAAATCGTGCAAATTCAATACATTGCAGAATTTTC‑‑GTAGGc                                                >  2:3227649/6‑100 (MQ=18)
                                                      gttgACAAAGTGCGCTTTGTTTATGCCGGATGCGGCGTTAACGCCTTATCCGGCCT‑‑ACAAAATCGTGCAAATTCAATACATTGCAGAATTTTC‑‑GTAGGcc                                               <  2:2211694/100‑1 (MQ=255)
                                                       taaaCAAAGCGCACTTTGTTTATGCCGGATGCGGCGTTAACGCCTTATCCGGCCT‑‑ACAAAATCGTGCAAATTCAATACATTGCAGAATTTTC‑‑GTAGGCCt                                              >  2:1983579/4‑100 (MQ=18)
                                                         gACAAAGTGCGCTTTGTTTATGCCGGATGCGGCGTTAACGCCTTATCCGGCCT‑‑ACAAAATCGTGCAAATTCAATACATTGCAGAATTTTC‑‑GTAGGCCTGa                                            >  2:2556248/1‑100 (MQ=255)
                                                             aaGTGCGCTTTGTTTATGCCGGATGCGGCGTTAACGCCTTATCCGGCCT‑‑ACAAAATCGTGCAAATTCAATACATTGCAGAATTTTC‑‑GTAGGCCTGATAAg                                        <  1:1216033/100‑1 (MQ=255)
                                                              aGTGCGCTTTGTTTATGCCGGATGCGGCGTTAACGCCTTATCCGGCCT‑‑ACAAAATCGTGCAAATTCAATACATTGCAGAATTTTC‑‑GTAGGCCTGATAAGc                                       <  1:723624/100‑1 (MQ=255)
                                                                 gcgcTTTGTTTATGCCGGATGCGGCGTTAACGCCTTATCCGGCCT‑‑ACAAAATCGTGCAAATTCAATACATTGCAGAATTTTC‑‑GTAGGCCTGATAAGCGGa                                    >  2:1394693/1‑100 (MQ=255)
                                                                      ttGTTTATGCCGGATGCGGCGTTAACGCCTTATCCGGCCT‑‑ACAAAATCGTGCAAATTCAATACATTGCAGAATTTTC‑‑GTAGGCCTGATAAGCGGAGCGCa                               <  2:1123426/100‑1 (MQ=255)
                                                                        gTTTATGCCGGATGCGGCGTTAACGCCTTATCCGGCCT‑‑ACAAAATCGTGCAAATTCAATACATTGCAGAATTTTC‑‑GTAGGCCTGATAAGCGGAGCGCATc                             <  2:3528445/100‑1 (MQ=255)
                                                                        gTTTATGCCGGATGCGGCGTTAACGCCTTATCCGGCCT‑‑ACAAAATCGTGCAAATTCAATACATTGCAGAATTTTC‑‑GTAGGCCTGATAAGCGGAGCGCATc                             >  2:2122604/1‑100 (MQ=255)
                                                                          gtaggCCGGATAAGGCGTTTACGCCGCATCCGGCCT‑‑ACAAAATCGTGCAAATTCAACATATTGCAATTCTCTT‑‑GTAGGCCTGATAAGCGTAGCGCATCAg                           >  1:2104335/5‑100 (MQ=18)
                                                                          gtaggCCGGATAAGGCGTTCACGCCGCATCCGGCCT‑‑ACAAAAATCTTGCCAATTCAATATATTGCAGGACCAT‑‑GTAGGCCTGATAAGCGTAGCGCATCAg                           <  1:2324286/96‑1 (MQ=18)
                                                                           taggCCGGATAAGGCGTTCACGCCGCATCCGGCCT‑‑ACAAAAATCTTGCCAATTCAATATATTGCAGGACCAT‑‑GTAGGCCTGATAAGCGTAGCGCATCAgg                          >  1:890024/4‑100 (MQ=18)
                                                                           taggCCGGATAAGGCGTTCACGCCGCATCCGGCCT‑‑ACAAAGTCGTATAAATTCAATGTATTGCAGTTTATTG‑‑GTAGGCCTGATAAGCGCAGCGCATCAgg                          >  1:450411/4‑100 (MQ=25)
                                                                           tATGCCGGATGCGGCGTTAACGCCTTATCCGGCCT‑‑ACAAAATCGTGCAAATTCAATACATTGCAGAATTTTC‑‑GTAGGCCTGATAAGCGGAGCGCATCAgg                          <  2:3018724/100‑1 (MQ=255)
                                                                            aggCCGGATAAGGAGATCACGCCGCATCCGGCCT‑‑ACAATACTGTGCAAATTCAATATATTGCAGGAAACAC‑‑GTAGGCCTGATAAGCCAAGCGCATCAGGc                         >  2:2434693/3‑100 (MQ=17)
                                                                             ggCCGGATAAGGGGCTTACGCCGCATCCGGCAT‑‑ACAAGAACGTGCAAATTCAACAAATTGTGATGCATCCT‑‑GTAGGCCTGATAAGCGTAGCGCATCAGGc                         >  2:746507/2‑100 (MQ=18)
                                                                             ggCCGGATAAGGCGTTTACGCCGCATCCGGCCT‑‑ACAAATTCTGCAAATTCAACAAATTGCGAATCCCTT‑‑GTAGGCCTGATAAGCGTAGCGCATCAGGCAt                       >  1:2442379/2‑100 (MQ=18)
                                                                             ggCCGGATAAGGCGTTTACGCCGCATCCGGCCT‑‑ACAAAAGCATGCAAATTCAATATATTGCAGAGTTCAC‑‑GTAGGCCTGATAAGCGTAGCGCATCAGGCa                        >  2:1631156/2‑100 (MQ=17)
                                                                             ggCCGGATAAGGCGTTTACGCCGCATCCGGCCT‑‑ACAAAAGCATGCAAATTCAATATATTGCAGAGTTCAC‑‑GTAGGCCTGATAAGCGTAGCGCATCAGGCa                        >  1:840017/2‑100 (MQ=17)
                                                                             ggCCGGATAAGGCGTTCACGCCGCATCCGGCCT‑‑ACAAAAATCTTGCCAATTCAATATATTGCAGGACCAT‑‑GTAGGCCTGATAAGCGTAGCGCATCAGGCa                        >  1:599739/2‑100 (MQ=18)
                                                                             ggCCGGATAAGGCGTTCACGCCGCATCCGGCCT‑‑ACAAGATCGTGCAAATTCAATATATTGCATGAGCCAC‑‑GTAGGCCTGATAAGCGTAGCGCATCAGGCa                        >  2:135189/2‑100 (MQ=17)
                                                                             ggCCGGATAAGGCGTTCACGCCGCATCCGGCCT‑‑ACAAAACTTTGCAAATTCAATATATTGCAATCTCCGTGTAGGCCTGATAAGCGTAGCGCATCAGGCa                        >  1:2194927/2‑100 (MQ=18)
                                                                             ggCCGGATAAGGCGTTCACGCCGCATCCGGCCT‑‑ACAAAACTTTGCAAATTCAATATATTGCAATCTCCGTGTAGGCCTGATAAGCGTAGCGCATCAGGCa                        >  1:1823785/2‑100 (MQ=18)
                                                                             ggCCGGATAAGGCGTTCACGCCGCATCCGGCCT‑‑ACAAAACCTTGCTAATTCAATATATTGCAGGGACTAT‑‑GTAGGCCTGATAAGCATAGCGCATCAGGCa                        >  1:2888249/2‑100 (MQ=18)
                                                                              gCCGGATAAGGCGTTTACGCCGCATCCGGCCT‑‑ACAAAATCGTGCTAATTCAATATATTGCAGAAACCTT‑‑GTAGGCCTGATAAGCGTAGCGCATCAGGCag                       >  1:2328186/1‑99 (MQ=18)
                                                                              gCCGGATAAGGCGTTCACGCCGCATCCGGCCT‑‑ACAAAAATCTTGCCAATTCAATATATTGCAGGACCAT‑‑GTAGGCCTGATAAGCGTAGCGCATCAGGCag                       >  2:3093156/1‑99 (MQ=18)
                                                                              gCCGGATAAGGCGTTCACGCCGCATCCGGCCT‑‑ACAAAAGCATGAAAATTCAATATATTGCAGGAGCTGC‑‑GTAGGCCTGATAAGCGTAGCGCATCAGGCag                       >  1:3892116/1‑99 (MQ=17)
                                                                              gCCGGATAAGGCGTTAACGCCTTATCCGGCCT‑‑ACAAAATCGTGCAAATTCAATACATTGCAGAATTTTC‑‑GTAGGCCTGATAAGCGGAGCGCATCAGGCAt                       >  2:202301/1‑100 (MQ=255)
                                                                               ctggataaggcgTTCACGCCGCATCCAGCCT‑‑ACAAAATTGTGCAAAATCAATGGATTGCACAGCTAGC‑‑GTAGGCCTGATAAGCGTAGCGCATCAGGCaga                      >  1:2728497/9‑98 (MQ=25)
                                                                               ccggataaggcgTTTACGCCGCATCCGGCCT‑‑ACAAAAGCATGCAAATTCAATATATTGCAGAGTTCAC‑‑GTAGGCCTGATAAGCGTAGCGCATCAGGCagt                      >  1:1746472/9‑98 (MQ=17)
                                                                               ccggataaggcgTTCACGCCGCATCCGGCCT‑‑ACAAAATCGCTTAAATTCAATATATTGCACCATCCTC‑‑GTAGGCCTGATAAGCGTAGCGCATCAGGCaat                      >  1:915809/9‑98 (MQ=21)
                                                                               ccggataaggcgTTCACGCCGCATCCGGCCT‑‑ACAAAAGCATGAAAATTCAATATATTGCAGGAGCTGC‑‑GTAGGCCTGATAAGCGTAGCGCATCAGGCagt                      >  1:3074817/9‑98 (MQ=18)
                                                                               ccggataaggcgTTCACGCCGCATCCGGCCT‑‑ACAAAAGCATGAAAATTCAATATATTGCAGGAGCTGC‑‑GTAGGCCTGATAAGCGTAGCGCATCAGGCagt                      >  1:302215/9‑98 (MQ=17)
                                                                               ccGGATAAGGCGTTCACGCCGCATCCGGCATGAACAAAAACATGCAAATTCAATAGATTGCAGAGATTAT‑‑GTAGGCCTGATAAGCGTAGCGCATCAGGCa                        >  1:201161/1‑100 (MQ=17)
                                                                                 ggataaggcgTTTACGCCGCATCTGGCCT‑‑ACAAATTCTTCCAAATTCGATATATTGCTGAAATCAT‑‑GTAGGCCTGATAAGCGTATCGCATCAGGCAattt                    >  1:148022/7‑100 (MQ=17)
                                                                                 ggataaggcgTTCACGCCGCATCCGGCCT‑‑ACAAAAATCTTGCCAATTCAATATATTGCAGGACCAT‑‑GTAGGCCTGATAAGCGTAGCGCATCAGGCAgttt                    >  2:1516479/7‑100 (MQ=18)
                                                                                 ggataaggcgTTCACGCCGCATCCGGCCT‑‑ACAAAGTCGTATAAATTCAATGTATTGCAGTTTATTG‑‑GTAGGCCTGATAAGCGCAGCGCATCAGGCaatct                    >  1:3203297/7‑96 (MQ=21)
                                                                                 ggataaggcgTTCACGCCGCATCCGGCCT‑‑ACAAAAGCATGAAAATTCAATATATTGCAGGAGCTGC‑‑GTAGGCCTGATAAGCGTAGCGCATCAGGCAgttt                    >  1:212303/7‑100 (MQ=18)
                                                                                 ggATGCGGCGTTAACGCCTTATCCGGCCT‑‑ACAAAATCGTGCAAATTCAATACATTGCAGAATTTTC‑‑GTAGGCCTGATAAGCGGAGCGCATCAGGCAtttt                    <  2:3883383/100‑1 (MQ=255)
                                                                                  gataaggcgTTTACGCCTTATCCGGCCT‑‑ACAAAATCGTGCTAATTCAATATATTGCAGAAACCTT‑‑GTAGGCCTGATAAGCGTAGCGCATCAGGCAgtttt                   >  2:1845336/6‑100 (MQ=18)
                                                                                   aTGCGGCGTTAACGCCTTATCCGGCCT‑‑ACAAAATCGTGCAAATTCAATACATTGCAGAATTTTC‑‑GTAGGCCTGATAAGCGGAGCGCATCAGGCATTtttg                  <  1:1586966/100‑1 (MQ=255)
                                                                                    taaggcgTTTACGCCGCATCCGGCCT‑‑ACAAAAGCGTGCAAATTCAATACATTGCATGGGCCAT‑‑GTAGGCCTGATAAGCGTAGCGCATCAGGCAattttac                 >  2:1429123/4‑98 (MQ=18)
                                                                                    tGCGGCGTTAACGCCTTATCCGGCCT‑‑ACAAACTCGTGTAAATTCAATACATTGCAGAATTTTC‑‑GTAGGCCTGATAAGCGGAGCGCATCAGGCATTtttgc                 <  2:2668498/100‑1 (MQ=255)
                                                                                      cggcgTTAACGCCTTATCCGGCCT‑‑ACAAAATCGTGCAAATTCAATACATTGCAGAATTTTC‑‑GTAGGCCTGATAAGCGGAGCGCATCAGGCATTtttgcct               >  1:36715/1‑100 (MQ=255)
                                                                                      cggcgTTAACGCCTTATCCGGCCT‑‑ACAAAATCGTGCAAATTCAATACATTGCAGAATTTTC‑‑GTAGGCCTGATAAGCGGAGCGCATCAGGCATTtttgcct               >  1:2247401/1‑100 (MQ=255)
                                                                                         cgTTAACGCCTTATCCGGCCT‑‑ACAAAATCGTGCAAATTCAATACATTGCAGAATTTTC‑‑GTAGGCCTGATAAGCGGAGCGCATCAGGCATTtttgcctttg            <  2:3595180/100‑1 (MQ=255)
                                                                                            tAACGCCTTATCCGGCCT‑‑ACAAAATCGTGCAAATTCAATACATTGCAGAATTTTC‑‑GTAGGCCTGATAAGCGGAGCGCATCAGGCATTTTTGCCTTTGCCa         <  1:1116322/100‑1 (MQ=255)
                                                                                             aaCGCCTTATCCGGCCT‑‑ACAAAATCGTGCAAATTCAATACATTGCAGAATTTTC‑‑GTAGGCCTGATAAGCGGAGCGCATCAGGCATTTTTGCCTTTGCCAt        >  2:1267549/1‑100 (MQ=255)
                                                                                              aCGCCTTATCCGGCCT‑‑ACAAAATCGTGCAAATTCAATACATTGCAGAATTTTC‑‑GTAGGCCTGATAAGCGGAGCGCATCAGGCATTTTTGCCTTTGCCATc       <  2:1180630/100‑1 (MQ=255)
                                                                                              aCGCCTTATCCGGCCT‑‑ACAAAATCGTGCAAATTCAATACATTGCAGAATTTTC‑‑GTAGGCCTGATAAGCGGAGCGCATCAGGCATTTTTGCCTTTGCCATc       >  2:1502175/1‑100 (MQ=255)
                                                                                               cGCCTTATCCGGCCT‑‑ACAAAATCGTGCAAATTCAATACATTGCAGAATTTTC‑‑GTAGGCCTGATAAGCGGAGCGCATCAGGCATTTTTGCCTTTGCCATCa      <  1:3749156/100‑1 (MQ=255)
                                                                                                gCCTTATCCGGCCT‑‑ACAAAATCGTGCAAATTCAATACATTGCAGAATTTTC‑‑GTAGGCCTGATAAGCGGAGCGCATCAGGCATTTTTGCCTTTGCCATCAg     >  1:3155332/1‑100 (MQ=255)
                                                                                                gCCTTATCCGGCCT‑‑ACAAAATCGTGCAAATTCAATACATTGCAGAATTTTC‑‑GTAGGCCTGATAAGCGGAGCGCATCAGGCATTTTTGCCTTTGCCATCAg     <  1:3345244/100‑1 (MQ=255)
                                                                                                gCCTTATCCGGCCT‑‑ACAAAATCGTGCAAATTCAATACATTGAAGAATTTTC‑‑GTAGGCCTGATAAGCGGAGCGCATCAGGCATTTTTGCCTTTGCCATCAg     <  2:257960/100‑1 (MQ=255)
                                                                                                  cTTATCCGGCCT‑‑ACAAAATCGTGCAAATTCAATACATTGCAGAATTTTC‑‑GTAGGCCTGATAAGCGGAGCGCATCAGGCATTTTTGCCTTTGCCATCAGTc   <  2:2867253/100‑1 (MQ=255)
                                                                                                   |                                                                                                           
TTGCTGCCGGATGGCGCTGCACGCAGAGGCACAGCGCCACCGGGTATACAGATTGTTGACAAAGTGCGCTTTGTTTATGCCGGATGCGGCGTTAACGCCTTATCCGGCCT‑‑ACAAAATCGTGCAAATTCAATACATTGCAGAATTTTC‑‑GTAGGCCTGATAAGCGGAGCGCATCAGGCATTTTTGCCTTTGCCATCAGTCT  >  NC_000913/3705885‑3706083

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