Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 3,705,971 C→A 29.6% intergenic (‑181/+127) 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,9710CA29.6% 515.7 / 234.8 261intergenic (‑181/+127)dppB/dppAdipeptide/heme ABC transporter permease/dipeptide/heme ABC transporter periplasmic binding protein; dipeptide chemotaxis receptor
Reads supporting (aligned to +/- strand):  ref base C (97/83);  new base A (71/5);  total (171/90)
Fisher's exact test for biased strand distribution p-value = 7.18e-11
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.85e-02

Only 100 of 275 total aligned reads displayed.

GGGAATGAAGTATTGCTGCCGGATGGCGCTGCACGCAGAGGCACAGCGCCACCGGGTATACAGATTGTTGACAAAGTGCGCTTTGTTTATGCCGGATGCGGCGTTAACGCCTTATCCGGCCTACAAAATCGTGCAAATTCAATACATTGC‑‑AG‑‑‑AATTTTC‑‑GTAGGCCTGATAAGCGGAGCGCATCAGGCATTTTTGCCT  >  NC_000913/3705873‑3706070
                                                                                                  |                                                                                                             
gggAATGAAGTATTGCTGCCGGATGGCGCTGCACGCAGAGGCACAGCGCCACCGGGTATACAGATTGTTGACAAAGTGCGCTTTGTTTATGCCGGATgcg                                                                                                              >  2:209240/1‑100 (MQ=255)
gggAATGAAGTATTGCTGCCGGATGGCGCTGCACGCAGAGGCACAGCGCCACCGGGTATACAGATTGTTGACAAAGTGCGCTTTGTTTATGCCGGATgcg                                                                                                              >  2:3549901/1‑100 (MQ=255)
 ggAATGAAGTATTGCTGCCGGATGGCGCTGCACGCAGAGGCACAGCGCCACCGGGTATACAGATTGTTGACAAAGTGCGCTTTGTTTATGCCGGATgcgg                                                                                                             <  2:2953282/100‑1 (MQ=255)
    aTGAAGTATTGCTGCCGGATGGCGCTGCACGCAGAGGCACAGCGCCACCGGGTATACAGATTGTTGACAAAGTGCGCTTTGTTTATGCCGGATGCGGCGt                                                                                                          >  1:3831039/1‑100 (MQ=255)
    aTGAAGTATTGCTGCCGGATGGCGCTGCACGCAGAGGCACAGCGCCACCGGGTATACAGATTGTTGACAAAGTGCGCTTTGTTTATGCCGGATGCGGCGt                                                                                                          <  1:291740/100‑1 (MQ=255)
      gAAGTATTGCTGCCGGATGGCGCTGCACGCAGAGGCACAGCGCCACCGGGTATACAGATTGTTGACAAAGTGCGCTTTGTTTATGCCGGATGCGGCGTTa                                                                                                        <  1:3546962/100‑1 (MQ=255)
         gTATTGCTGCCGGATGGCGCTGCACGCAGAGGCACAGCGCCACCGGGTATACAGATTGTTGACAAAGTGCGCTTTGTTTATGCCGGATGTGGCGTTAACg                                                                                                     >  1:669645/1‑100 (MQ=255)
          tATTGCTGCCGGATGGCGCTGCACGCAGAGGCACAGCGCCACCGGGTATACAGATTGTTGACAAAGTGCGCTTTGTTTATGCCGGATGCGGCGTTAACGc                                                                                                    <  2:505637/100‑1 (MQ=255)
           aTTGCTGCCGGATAGCGCTGCACGCAGAGGCACAGCGCAACCGGGTATACAGATTGTTGACAAAGTGCGCTTTGTTTATGCCGGATGCGGCGTTAACGcc                                                                                                   >  1:2722073/1‑100 (MQ=255)
           aTTGATGCCTGATGGCGCTGCACGTAGAGGCACAGCGCCACCGGGTATACAGATTGTTGACAAAGTGCGCTTTGTTTATGCCGGATGCGGCGTTAACGcc                                                                                                   <  1:348662/100‑1 (MQ=255)
            ttGCTGCCGGATGGCGCTGCACGCAGAGGCACAGCGCCACCGGGTATACAGATTGTTGACAAAGTGCGCTTTGTTAACGCCGCATccggcataaacaaag                                                                                                  <  2:112005/100‑11 (MQ=37)
               ctgcCGGATGGCGCTGCACGCGGAGGCACAGCGCCACGGGGTATACAGATTGTTGACAAAGTGCGCTTTGTTTATGCCGGATGCGGCGTTAACGCCTTAt                                                                                               >  1:2363313/1‑100 (MQ=255)
               ctgcCGGATGGCGCTGCACGCAGAGGCACAGCGCCACCGGGTATACAGATTGTTGACAAAGTGCGCTTTGTTTATGCCGGATGCGGCGTTAACGCCTTAt                                                                                               >  2:1943922/1‑100 (MQ=255)
               ctgcCGGATGGCGCTGCACGCAGAGGCACAGCGCCACCGGGTATACAGATTGTTGACAAAGTGCGCTTTGTTTATGCCGGATGCGGCGTTAACGCCTTAt                                                                                               >  1:1550010/1‑100 (MQ=255)
               ctgcCGGATGGCGCTGCACGCAGAGGCACAGCGCCACCGGGTATACAGATTGTTGACAAAGTGCGCTTTGTTTATGCCGGATGCGGCGTTAACGCCTTAt                                                                                               >  1:3288729/1‑100 (MQ=255)
                  ccGGATGGCGCTGCACGCAGAGGCACAGCGCCACCGGGTATACAGATTGTTGACAAAGTGCGCTTTGTTTATGCCGGATGCGGCGTTAACGCCTTATCCg                                                                                            >  2:3346639/1‑100 (MQ=255)
                  ccGGATGGCGCTGCACGCAGAGGCACAGCGCCACCGGGTATACAGATTGTTGACAAAGTGCGCTTTGTTTATGCCGGATGCGGCGTTAACGCCTTATCCg                                                                                            >  2:1244450/1‑100 (MQ=255)
                   cGGATGGCGCTGCACGCAGAGGCACAGCGCCACCGGGTATACAGATTGTTGACAAAGTGCGCTTTGTTTATGCCGGATAAGGCGTTAACGCCGCATCCgg                                                                                           <  1:4219811/100‑1 (MQ=37)
                    ggATGGCGCTGCACGCAGAGGCACAGCGCCACCGGGTATACAGATTGTTGACAAAGTGCGCTTTGTTTATGCCGGATGCGGCGTTAACGCCTTATCCGGc                                                                                          >  2:2179141/1‑100 (MQ=255)
                     gATGGCGCTGCACGCAGAGGCACAGCGCCACCGGGTATACAGATTGTTGACAAAGTGCGCTTTGTTTATGCCGGATGCGGCGTTAACGCCTTATCCGGcc                                                                                         <  1:3346619/100‑1 (MQ=255)
                     gATGGCGCTGCACGCAGAGGCACAGCGCCACCGGGTATACAGATTGTTGACAAAGTGCGCTTTGTTTATGCCGGATGCGGCGTTAACGCCTTATCCGGcc                                                                                         <  2:3028986/100‑1 (MQ=255)
                     gATGGCGCTGCACGCAGAGGCACAGCGCCACCGGGTATACAGATTGTTGACAAAGTGCGCTTTGTTTATGCCGGATAAGGCGTTAACGCCGCATCCGGca                                                                                         <  1:2237018/100‑2 (MQ=37)
                     gATGGCGCTGCACGCAGAGGCACAGCGCCACCGGGTATACAGATTGTTGACAAAGTGCGCTTTCTTTATGCCGGATGCGGCGTTAACGCCTTATCCGGcc                                                                                         <  2:3113597/100‑1 (MQ=255)
                      aTGGCGCTGCACGCAGAGGCACAGCGCCACCGGGTATACAGATTGTTGACAAAGTGCGCTTTGTTTATGCCGGATGCGGCGTTAACGCCTTATCCGGCCt                                                                                        <  2:2704199/100‑1 (MQ=255)
                         gcgcTGCACGCAGAGGCACAGCGCCACCGGGTATACAGATTGTTGACAAAGTGCGCTTTGTTTATGCCGTATGCGGCGTTAACGCCTTATCCGGCCTACa                                                                                    >  1:3638596/1‑100 (MQ=255)
                          cgcTGCACGCAGAGGCACAGCGCCACCGGGTATACAGATTGTTGACAAAGTGCGCTTTGTTTATGCCGGATGCGGCGTTAACGCCTTATCCGGCCTACaa                                                                                   >  1:2433928/1‑100 (MQ=255)
                          cgcTGCACGCAGAGGCACAGCGCCACCGGGTATACAGATTGTTGACAAAGTGCGCTTTGTTTATGCCGGATGCGGCGTTAACGCCTTATCCGGCCTACaa                                                                                   >  1:981849/1‑100 (MQ=255)
                           gcTGCACGCAGAGGCACAGCGCCACCGGGTATACAGATTGTTGACAAAGTGCGCTTTGTTTATGCCGGATGCGGCGTTAACGCCTTATCCGGCCTACaaa                                                                                  >  2:836792/1‑100 (MQ=255)
                              gCACGCAGAGGCACAGCGCCACCGGGTATACAGATTGTTGACAAAGTGCGCTTTGTTTATGCCGGATAAGGCGTTAACGCCGCATCCGGcataaacaaag                                                                                <  1:125109/100‑11 (MQ=21)
                               cACGCAGAGGCACAGCGCCACCGGGTATACAGATTGTTGACAAAGTGCGCTTTGTTTATGCCGGATGCGGCGTTAACGCCTTATCCGGCCTACAAAATCg                                                                              >  2:3239664/1‑100 (MQ=255)
                                 cGCAGAGGCACAGCGCCACCGGGTATACAGATTGTTGACAAAGTGCGCTTTGTTTATGCCGGATGCGGCGTTAACGCCTTATCCGGCCTACAAAATCGTg                                                                            >  2:1629355/1‑100 (MQ=255)
                                   cAGAGGCACAGCGCCACCGGGTATACAGATTGTTGACAAAGTGCGCTTTGTTTATGCCGGATGCGGCGTTAACGCCTTATCCGGCCTACAAAATCGTGCa                                                                          >  1:2284197/1‑100 (MQ=255)
                                   cAGAGGCACAGCGCCACCGGGTATACAGATTGTTGACAAAGTGCGCTTTGTTTATGCCGGATGCGGCGTTAACGCCTTATCCGGCCTACAAAATCGTGCa                                                                          >  2:2648129/1‑100 (MQ=255)
                                   cAGAGGCACAGCGCCACCGGGTATACAGATTGTTGACAAAGTGCCCTTTGTTTATGCCGGATGCGGCGTTAACGCCTTATCCGGCCTACAAAATCGTGCa                                                                          <  1:3082814/100‑1 (MQ=255)
                                       ggCACAGCGCCACCGGGTATACAGATTGTTGACAAAGTGCGCTTTGTTTATGCCGGATGTGGCGTTAACGCCTTATCCGGCCTACAAAATCGTGCAAAtt                                                                      >  1:430781/1‑100 (MQ=255)
                                       ggCACAGCGCCACCGGGTATACAGATTGTTGACAAAGTGCGCTTTGTTTATGCCGGATGCGGCGTTAACGCCTTATCCGGCCTACAAAATCGTGCAAAtt                                                                      >  2:1538432/1‑100 (MQ=255)
                                        gCACAGCGCCACCGGGTATACAGATTGTTGACAAAGTGCGCTTTGTTTATGCCGGATGCGGCGTTAACGCCTTATCCGGCCTACAAAATCGTGCAAATTc                                                                     <  2:2374408/100‑1 (MQ=255)
                                        gCACAGCGCCACCGGGTATACAGATTGTTGACAAAGTGCGCTTTGTTTATGCCGGATGCGGCGTTAACGCCTTATCCGGCCTACAAAATCGTGCAAATTc                                                                     >  1:2949552/1‑100 (MQ=255)
                                         cacaGCGCCACCGGGTATACAGATTGTTGACAAAGTGCGCTTTGTTTATGCCGGATGCGGCGTTAACGCCTTATCCGGCCTACAAAATCGTGCAAATTCa                                                                    >  1:1063379/1‑100 (MQ=255)
                                           caGCGCCACCGGGTATACAGATTGTTGACAAAGTGCGCTTTGTTTATGCCGGATGCGGCGTTAACGCCTTATCCGGCCTACAACATCGTGCAAATTCAAt                                                                  >  1:202565/1‑100 (MQ=255)
                                           caGCGCCACCGGGTATACAGATTGTTGACAAAGTGCGCTTTGTTTATGCCGGATGCGGCGTTAACGCCTTATCCGGCCTACAAAATCGTGCAAATTCAAt                                                                  <  2:140986/100‑1 (MQ=255)
                                            aGCGCCACCGGATATACAGATTGTTGACAAAGTGCGCTTTGTTTATGCCGGATGCGGCGTTAACGCCTTATCCGGCCTACAAAATCGTGCAAATTCAAt                                                                  <  1:1886058/99‑1 (MQ=255)
                                                ccACCGGGTATACAGATTGTTGACAAAGTGCGCTTTGTTTATGCCGGATGCGGCGTTAACGCCTTATCCGGCCTACAAAATCGTGCAAATTCAATACAtt                                                             >  2:3938411/1‑100 (MQ=255)
                                                ccACCGGGTATACAGATTGTTGACAAAGTGCGCTTTGTTTATGCCGGATGCGGCGTTAACGCCTTATCCGGCCTACAAAATCGTGCAAATTCAATACAtt                                                             >  2:2841549/1‑100 (MQ=255)
                                                  aCCGGGTATACAGATTGTTGACAAAGTGCGCTTTGTTTATGCCGGATGCGGCGTTAACGCCTTATCCGGCCTACAAAATCGTGCAAATTCAATACATTGc                                                           <  1:2591134/100‑1 (MQ=255)
                                                    cGGGTATACAGATTGTTGACAAAGTGCGCTTTGTTTATGCCGGATGCGGCGTTAACGCCTTATCCGGCCTACAAAATCGTGCAAATTCAATACATTGC‑‑a                                                        <  1:169121/99‑1 (MQ=255)
                                                    cGGGTATACAGATTGTTGACAAAGTGCGCTTTGTTTATGCCGGATGCGGCGTTAACGCCTTATCCGGCCTACAAAATCGTGCAAATTCAATACATTGC‑‑Ag                                                       <  2:378939/100‑1 (MQ=255)
                                                      ggTATACAGATTGTTGACAAAGTGCGCTTTGTTTATGCCGGATGCGGCGTTAACGCCTTATCCGGCCTACAAAATCGTGCAAATTCAATACATTGC‑‑AG‑‑‑aa                                                 >  2:2171039/1‑100 (MQ=255)
                                                        tataCAGATTGTTGACAAAGTGCGCTTTGTTTATGCCCGATTCGGCGTTAACGACTTATCCGGCCTACAAAATCGTGCAGATTCAATACATTGC‑‑AG‑‑‑AAtt                                              >  1:1484584/1‑100 (MQ=18)
                                                            cAGATTGTTGACAAAGTGCGCTTTGTTTATGCCGGATGCGGCGTTAACGCCTTATCCGGCCTACAAAATCGTGCAAATTCAATACATTGC‑‑AG‑‑‑AATTTTC‑‑g                                        >  2:2293739/1‑100 (MQ=255)
                                                            cAGATTGTTGACAAAGTGCGCTTTGTTTATGCCGGATGCGGCGTTAACGCCTTATCCGGCCTACAAAATCGTGCAAATTCAATACATTGC‑‑AG‑‑‑AATTTTC‑‑g                                        <  2:2995149/100‑1 (MQ=255)
                                                             aGATTGTTGACAAAGTGCGCTTTGTTTATGCCGGATGCGGCGTTAACGCCTTATCCGGCCTACAAAATCGTGCAAATTCAATACATTGC‑‑AG‑‑‑AATTTTC‑‑g                                        <  1:3023133/99‑1 (MQ=255)
                                                             aGATTGTTGACAAAGTGCGCTTTGTTTATGCCGGATGCGGCGTTAACGCCTTATCCGACCTACAACATCGTGCAAGTTTAATCCATTGC‑‑AG‑‑‑AATTTTC‑‑Gt                                       >  2:467455/1‑100 (MQ=18)
                                                               ggcataaaCAAAGCGCACTTTGTTTATGCCGGATGCGGCGTTAACGCCTTATCCGGCCTACAAAATCGTGCAAATTCAATACATTGC‑‑AG‑‑‑AATTTTC‑‑GTAg                                     >  2:2017649/8‑100 (MQ=17)
                                                                ttgttgACAAAGTGCGCTTTGTTTATGCCGGATGCGGCGTTAACGCCTTATCCGGCCTACAAAATCGTGCAAATTCAATACATTGC‑‑AG‑‑‑AATTTTC‑‑GTAgg                                    <  1:247495/100‑1 (MQ=255)
                                                                  gttgACAAAGTGCGCTTTGTTTATGCCGGATGCGGCGTTAACGCCTTATCCGGCCTACAAAATCGTGCAAATTCAATACATTGC‑‑AG‑‑‑AATTTTC‑‑GTAGGcc                                  <  2:715622/100‑1 (MQ=255)
                                                                  gttgACAAAGTGCGCTTTGTTTATGCCGGATGCGGCGTTAACGCCTTATCCGGCCTACAAAATCGTGCAAATTCAATACATTGC‑‑AG‑‑‑AATTTTC‑‑GTAGGcc                                  >  2:3589071/1‑100 (MQ=255)
                                                                    tgACAAAGTGCGCTTTGTTTATGCCGGATGCGGCGTTAACGCCTTATCCGGCCTACAAAATCGTGCAAATTCAATACATTGC‑‑AG‑‑‑AATTTTC‑‑GTAGGCCTg                                >  1:938256/1‑100 (MQ=255)
                                                                            tGCGCTTTGTTTATGCCGGATGCGGCGTTAACGCCTTATCCGGCCTACAAAATCGTGCAAATTCAATACATTGC‑‑AG‑‑‑AATTTTC‑‑GTAGGCCTGATAAGCgg                        <  1:2205020/100‑1 (MQ=255)
                                                                             gcgcTTTGTTTATGCCGGATGCGGCGTTAACGCCTTATCCGGCCTACAAAATCGTGCAAATTCAATACATTGC‑‑AG‑‑‑AATTTTC‑‑GTAGGCCTGATAAGCGGa                       <  2:2691949/100‑1 (MQ=255)
                                                                                   tGTTTATGCCGGATGCGGCGTTAACGCCTTATCCGGCCTACAAAATCGTGCAAATTCAATACATTGC‑‑AG‑‑‑AATTTTC‑‑GTAGGCCTGATAAGCGGAGCGCAt                 >  2:3772392/1‑100 (MQ=255)
                                                                                   cttgtaggCCGGATAAGGCGTTCACGCCGCATCCGGCCTACAAGATCGTGCAAATTCAATATATTGC‑‑AT‑‑‑GAGCCAC‑‑GTAGGCCTGATAAGCGTAgcgc                   >  2:945490/8‑98 (MQ=17)
                                                                                    ttgtaggCCGGATAAGGCGTTTACGCCGCATCCGGCCTACAAAATCGTGCAAATTCAACATATTGC‑‑‑‑‑‑‑AATTCTCTTGTAGGCCTGATAAGCGTAGCGCATc                >  2:4343869/7‑100 (MQ=18)
                                                                                    gTTTATGCCGGATGCGGCGTTAACGCCTTATCCGGCCTACAAAATCGTGCAAATTCAATACATTGC‑‑AG‑‑‑AATTTTC‑‑GTAGGCCTGATAAGCGGAGCGCATc                <  2:2062359/100‑1 (MQ=255)
                                                                                    gTTTATGCCGGATGCGGCGTTAACGCCTTATCCGGCCTACAAAATCGTGCAAATTCAATACATTGC‑‑AG‑‑‑AATTTTC‑‑GTAGGCCTGATAAGCGGAGCGCATc                >  1:141880/1‑100 (MQ=255)
                                                                                    gTTTATGCCGGATGCGGCGTTAACGCCTTATCCGGCCTACAAAATCGTGCAAATTCAATACATTGC‑‑AG‑‑‑AATTTTC‑‑GTAGGCCTGATAAGCGGAGCGCATc                <  2:1791456/100‑1 (MQ=255)
                                                                                    gTTTATGCCGGATGCGGCGTTAACGCCTTATCCGGCCTACAAAATCGTGCAAATTCAATACATTGC‑‑AG‑‑‑AATTTTC‑‑GTAGGCCTGATAAGCGGAGCGCATc                >  2:2231739/1‑100 (MQ=255)
                                                                                      gtaggCCGGATAAGGCGTTCACGCTGCATCCGGCCTACAAAATCGCTTAAATTCAATATATTGC‑‑AC‑‑‑CATCCTC‑‑GTAGGCCTGATAAGCGTAGCGCATCAg              >  2:645012/5‑100 (MQ=18)
                                                                                      gtaggCCGGATAAGGCGTTCACGCCGCATCCGGCCTACAAAAATCTTGCCAATTCAATATATTGC‑‑AG‑‑‑GACCAT‑‑GTAGGCCTGATAAGCGTAGCGCATCAg              >  2:1889616/5‑100 (MQ=18)
                                                                                      gtaggCCGGATAAGGCGTTCACGCCGCATCCGGCCTACAAAAACCTGCAAATTCAATATATTGC‑‑AG‑‑‑GATACGT‑‑GTAGGCCTGATAAGCGTAGCGCATCAg              >  2:2582802/5‑100 (MQ=18)
                                                                                       taggCCGGATAAGGCGTTCACGCCGCATCCGGCCTACAAAACTGTGCAAATTCAATATATTGC‑‑AG‑‑‑GAAACAC‑‑GTAGGCCTGATAAGCGAAGCGCATCAgg             >  1:2266878/4‑100 (MQ=18)
                                                                                       taggCCGGATAAGGCGTTCACGCCGCATCCGGCCTACAAAAACATGCAAATTCAATAGATTGC‑‑AG‑‑‑AGATTAT‑‑GTAGGCCTGATAAGCGTAGCGCATCAgg             >  2:2553222/4‑100 (MQ=17)
                                                                                       taggCCGGATAAGGCGTTCACGCCGCATCCGGCCTACAAACCGCGCAAATTCAATATATTGCGGAG‑‑‑AAAAT‑‑‑‑GTAGGCCTGATAAGCGTAGCGCATCAGGc            >  1:3145920/4‑100 (MQ=18)
                                                                                       tATGCCGGATGCGGCGTTAACGCCTTATCCGGCCTACAAAATCGTGCAAATTCAATACATTGC‑‑AG‑‑‑AATTTTC‑‑GTAGGCCTGATAAGCGGAGCGCATCAgg             >  1:3266774/1‑100 (MQ=255)
                                                                                       tATGCCGGATGCGGCGTTAACGCCTTATCCGGCCTACAAAATCGTGCAAATTCAATACATTGC‑‑AG‑‑‑AATTTTC‑‑GTAGGCCTGATAAGCGGAGCGCATCAgg             <  1:3801484/100‑1 (MQ=255)
                                                                                        aggCCGGATAAGGCGTTTACGCCGCATCTGGCCTACAAATTCTTCCAAATTCAATATATTGC‑‑AG‑‑‑AAATCAT‑‑GTAGGCCTGATAAGCGTAGCGCATCAGGc            >  1:3205087/3‑100 (MQ=18)
                                                                                        aggCCGGATAAGGCGTTTACGCCGCATCCGGCCTACAAAAATGTGCAAATTCAATAAATTGC‑‑AATTCAACTT‑‑‑‑GTAGGCCTGATAAGCGCAGCGCATCAgg             >  2:4188952/3‑100 (MQ=21)
                                                                                        aggCCGGATAAGGCGTTTACGCCGCATCCGGCATACAAAAGCGTGCAAATTCAATACATTGC‑‑AT‑‑‑GGGCCAT‑‑GTAGGCCTGATAAGCGTAGCGCATCAGGc            >  1:4033404/3‑100 (MQ=17)
                                                                                        aggCCGGATAAGGCGTTCACGCCGCATCCGGCCTACAAAGTCGTATAAATTCAATGTATTGC‑‑AG‑‑‑TTTATTG‑‑GTAGGCCTGATAAGCGCAGCGCATCAGGc            >  2:2096103/3‑100 (MQ=21)
                                                                                        aggCCGGATAAGGCGTTCACGCCGCATCCGGCCTACAAAAGCATGAAAATTCAATATATTGC‑‑AG‑‑‑GAGCTGC‑‑GTAGGCCTGATAAGCGTAGCGCATCAGGc            >  2:1713708/3‑100 (MQ=18)
                                                                                        aTGCCGGATGCGGCGTTAACGCCTTATCCGGCCTACAAAATCGTGCAAATTCAATACATTGC‑‑AG‑‑‑AATTTTC‑‑GTAGGCCTGATAAGCGGAGCGCATCAGGc            >  2:380383/1‑100 (MQ=255)
                                                                                         tGCCGGATGCGGCGTTAACGCCTTATCCGGCCTACAAAATCGTGCAAATTCAATACATTGC‑‑AG‑‑‑AATTTTC‑‑GTAGGCCTGATAAGCGGAGCGCATCAGGCa           >  2:685377/1‑100 (MQ=255)
                                                                                         ggCCGGATAAGGCGTTTACGCCGCATCCGGCCTACAAATTCTGCAAATTCAACAAATTGC‑‑G‑‑‑AATCCCTT‑‑GTAGGCCTGATAAGCGTAGCGCATCAGGCAt          >  2:4283658/2‑100 (MQ=18)
                                                                                         ggCCGGATAAGGCGTTTACGCCGCATCCGGCCTACAAAATCGTGCTAATTCACTATATTGC‑‑AG‑‑‑AAACCTT‑‑GTAGGCCTGATAAGCGTAGCGCATCAGGCa           >  1:732366/2‑100 (MQ=17)
                                                                                         ggCCGGATAAGGCGTTTACGCCGCATCCGGCCTACAAAATCGTGCAAATTCAACATATTGC‑‑AA‑‑‑TTCTCTT‑‑GTAGGCCTGATAAGCGTAGCGCATCAGGct           >  2:868667/2‑99 (MQ=18)
                                                                                         ggCCGGATAAGGCGTTTACGCCGCATCCGGCCTACAAAAGCATGCAAATTCAATATATTGC‑‑AG‑‑‑AGTTCAC‑‑GTAGGCCTGATAAGCGTAGCGCATCAGGCa           >  1:1609293/2‑100 (MQ=14)
                                                                                         ggCCGGATAAGGCGTTCACTCCGCATCCGGCCTACAAAAGCACGCAAATTCAATATATTGC‑‑AG‑‑‑AGATTGC‑‑GTAGGCCTGATAAGCGTAGCGCATCAGGCa           >  1:4098772/2‑100 (MQ=11)
                                                                                         ggCCGGATAAGGCGTTCACGCCGCATCCGGCCTACAAAATCATGAAAATTCAATACATTGC‑‑AA‑‑‑GATTTTC‑‑GTAGGCCTGATAAGCGTAGCGCATCAGGCa           >  1:1361450/2‑100 (MQ=17)
                                                                                         ggCCGGATAAGGCGTTAACGCCGCATCCGGCCTACAAAATCGTGCAAATTCAATACATTGC‑‑AG‑‑‑AATTTTC‑‑GTAGGCCTGATAAGCGGAGCGCATCAGGCa           >  1:1816571/2‑100 (MQ=18)
                                                                                          gCCGGATAAGGCGTTTACGCCGCATCCGGCCTACAAATTCTGCAAATTCAACAAATTGC‑‑G‑‑‑AATCCCTT‑‑GTAGGCCTGATAAGCGTAGCGCATCAGGCAtt         >  2:2747060/1‑100 (MQ=18)
                                                                                          gCCGGATAAGGCGTTCACGCCGCATCCGGCCTACAAAAATCTTGCCAATTCAATATATTGC‑‑AG‑‑‑GACCAT‑‑GTAGGCCTGATAAGCGTAGCGCATCAGGCag          >  2:1226401/1‑99 (MQ=18)
                                                                                          gCCGGATAAGGCGTTCACGCCGCATCCGGCCTACAAGATCGTGCAAATTCAATATATTGC‑‑AT‑‑‑GAGCCAC‑‑GTAGTCCTGATAAGCGTAGCGCATCAGGCag          >  1:1008699/1‑99 (MQ=17)
                                                                                          gCCGGATAAGGCGTTCACGCCGCATCCGGCCTACAAAAACATGCAAATTCAATAGATTGC‑‑AG‑‑‑AGATTAT‑‑GTAGGCCTGATAAGCGTAGCGCATCAGGCaa          >  2:3361635/1‑99 (MQ=17)
                                                                                          gCCGGATAAGGCGTTCACGCCGCATCCGGCCTACAAATCATGCAAATCCAAAATATTGC‑‑AG‑‑‑GAGTTAT‑‑GTAGGCCTGATAAGCGTAGCGCATCAGGCAtt         >  1:1330197/1‑100 (MQ=21)
                                                                                          gCCGGATAAGGCGTTAACGCCGCATCCGGCCTACAAAATCGTGCAAATTCAATACATTGC‑‑AG‑‑‑AATTTTC‑‑GTAGGCCTGATAAGCGGAGCGCATCAGGCa           >  1:3517417/1‑99 (MQ=17)
                                                                                           ccGGATGCGGCGTTAACGCCTTATCCGGCCTACAAAATCGTGCAAATTCAATACATTGC‑‑AG‑‑‑AATTTTC‑‑GTAGGCCTGATAAGCGGAGCGCATCAGGCAtt         >  2:1357854/1‑100 (MQ=255)
                                                                                              gATGCGGCGTTAACGCCTTATCCGGCCTACAAAATCGTGCAAATTCAATACATTGC‑‑AG‑‑‑AATTTTC‑‑GTAGGCCTGATAAGCGGAGCGCATCAGGCAttttt      <  1:1232349/100‑1 (MQ=255)
                                                                                               aTGCGGCGTTAACGCCTTATCCTGCCTACAAAATCGTGCAAATTCAATACATTGT‑‑AG‑‑‑AATTTTC‑‑GTAGGCCTGATAAGCGGAGCGCATCAGGCATTtttg     <  1:802814/100‑1 (MQ=255)
                                                                                               aTGCGGCGTTAACGCCTTATCCGGCCTACAAAATCGTGCAAATTCAATACATTGC‑‑AG‑‑‑AATTTTC‑‑GTAGGCCTGATAAGCGGAGCGCATCAGGCATTtttg     <  2:2461829/100‑1 (MQ=255)
                                                                                                  cggcgTTAACGCCTTATCCGGCCTACAAAATCGTGCAAATTCAATACATTGC‑‑AG‑‑‑AATTTTC‑‑GTAGGCCTGATAAGCGGAGCGCATCAGGCATTtttgcct  >  2:4160287/1‑100 (MQ=255)
                                                                                                  |                                                                                                             
GGGAATGAAGTATTGCTGCCGGATGGCGCTGCACGCAGAGGCACAGCGCCACCGGGTATACAGATTGTTGACAAAGTGCGCTTTGTTTATGCCGGATGCGGCGTTAACGCCTTATCCGGCCTACAAAATCGTGCAAATTCAATACATTGC‑‑AG‑‑‑AATTTTC‑‑GTAGGCCTGATAAGCGGAGCGCATCAGGCATTTTTGCCT  >  NC_000913/3705873‑3706070

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