Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 3,705,970 G→A 29.0% intergenic (‑180/+128) 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,9700GA29.0% 468.1 / 200.3 261intergenic (‑180/+128)dppB/dppAdipeptide/heme ABC transporter permease/dipeptide/heme ABC transporter periplasmic binding protein; dipeptide chemotaxis receptor
Reads supporting (aligned to +/- strand):  ref base G (98/85);  new base A (70/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 = 4.40e-02

Only 100 of 275 total aligned reads displayed.

CAGGGAATGAAGTATTGCTGCCGGATGGCGCTGCACGCAGAGGCACAGCGCCACCGGGTATACAGATTGTTGACAAAGTGCGCTTTGTTTATGCCGGATGCGGCGTTAACGCCTTATCCGGCCTACAAAATCGTGCAAATTCAATACATTGCAG‑‑‑AATTTTC‑‑GTAGGCCTGATAAGCGGAGCGCATCAGGCATTTTTGCC  >  NC_000913/3705871‑3706069
                                                                                                   |                                                                                                           
cAGGGAATGAAGTATTGCTGCCGGATGGCGCTGCACGCAGAGGCACAGCGCCACCGGGTATACAGATTGTTGACAAAGTGCGCTTTGTTTATGCCGGATg                                                                                                             >  2:2318646/1‑100 (MQ=255)
  gggAATGAAGTATTGCTGCCGGATGGCGCTGCACGCAGAGGCACAGCGCCACCGGGTATACAGCTTGTTGACAAAGTGCGCTTTGTTTATGCCGGATgcg                                                                                                           >  1:888111/1‑100 (MQ=255)
  gggAATGAAGTATTGCTGCCGGATGGCGCTGCACGCAGAGGCACAGCGCCACCGGGTATACAGATTGTTGACAAAGTGCGCTTTGTTTATGCCGGATgcg                                                                                                           >  2:209480/1‑100 (MQ=255)
  gggAATGAAGTATTGCTGCCGGATGGCGCTGCACGCAGAGGCACAGCGCCACCGGGTATACAGATTGTTGACAAAGTGCGCTTTGTTTATGCCGGATgcg                                                                                                           >  2:3549901/1‑100 (MQ=255)
   ggAATGAAGTATTGCTGCCGGATGGCGCTGCACGCAGAGGCACAGCGCCACCGGGTATACAGATTGTTGACAAAGTGCGCTTTGTTTATGCCGGATgcgg                                                                                                          <  2:2953282/100‑1 (MQ=255)
      aTGAAGTATTGCTGCCGGATGGCGCTGCACGCAGAGGCACAGCGCCACCGGGTATACAGATTGTTGACAAAGTGCGCTTTGTTTATGCCGGATGCGGCGt                                                                                                       >  1:3028968/1‑100 (MQ=255)
      aTGAAGTATTGCTGCCGGATGGCGCTGCACGCAGAGGCACAGCGCCACCGGGTATACAGATTGTTGACAAAGTGCGCTTTGTTTATGCCGGATGCGGCGt                                                                                                       >  2:1241300/1‑100 (MQ=255)
         aaGTATTGCTGCCGGATGGCGCTGCACGCAGAGGCACAGCGCCACCGGGTATACAGATTGTTGACAAAGTGCGCTTTGTTTATGCCGGATGCGGCGTTaa                                                                                                    <  2:1649013/100‑1 (MQ=255)
            tATTGCTGCCGGATGGCGCTGCACGCAGAGGCACAGCGCCACCGGGTATACAGATTGTTGACAAAGTGCGCTTTGTTTATGCCGGATGCGGCGTTAACGc                                                                                                 <  2:505637/100‑1 (MQ=255)
             aTTGCTGCCGGATGGCGCTGCACGCAGAGGGACAGCGCCACCGGGTATACAGATTGTTGACAAAGTGCGCTTTGTTTATGCCGGATGCGGCGTTAACGcc                                                                                                >  1:701439/1‑100 (MQ=255)
             aTTGCTGCCGGATAGCGCTGCACGCAGAGGCACAGCGCAACCGGGTATACAGATTGTTGACAAAGTGCGCTTTGTTTATGCCGGATGCGGCGTTAACGcc                                                                                                >  1:2722073/1‑100 (MQ=255)
               tgctgcCGGATGGCGCTGCACGCAGAGGCACAGCGCCACCGGGTATACAGATTGTTGACAAAGTGCGCTTTGTTTATGCCGGATGCGGCGTTAACGCCtt                                                                                              <  1:2929442/100‑1 (MQ=255)
                 ctgcCGGATGGCGCTGCACGCGGAGGCACAGCGCCACGGGGTATACAGATTGTTGACAAAGTGCGCTTTGTTTATGCCGGATGCGGCGTTAACGCCTTAt                                                                                            >  1:2363313/1‑100 (MQ=255)
                 ctgcCGGATGGCGCTGCACGCAGAGGCACAGCGCCACCGGGTATACAGATTGTTGACAAAGTGCGCTTTGTTTATGCCGGATGCGGCGTTAACGCCTTAt                                                                                            >  1:1550010/1‑100 (MQ=255)
                 ctgcCGGATGGCGCTGCACGCAGAGGCACAGCGCCACCGGGTATACAGATTGTTGACAAAGTGCGCTTTGTTTATGCCGGATGCGGCGTTAACGCCTTAt                                                                                            >  2:1254826/1‑100 (MQ=255)
                   gcCGGATGGCGCTGCACGCAGAGGCACAGCGCGACCGGGTATACAGATTGTTGACAAAGTGCGCTTTGTTTATGCCGGATGCGGCGTTAACGGCTTATcc                                                                                          >  1:3526511/1‑100 (MQ=255)
                    ccGGATGGCGCTGCACGCAGAGGCACAGCGCCACCGGGTATACAGATTGTTGACAAAGTGCGCTTTGTTTATGCCGGATGCGGCGTTAACGCCTTATCCg                                                                                         >  2:3346639/1‑100 (MQ=255)
                       gATGGCGCTGCACGCAGAGGCACAGCGCCACCGGGTATACAGATTGTTGACAAAGTGCGCTTTGTTTATGCCGGATGCGGCGTTAACGCCTTATCCGGcc                                                                                      <  2:3028986/100‑1 (MQ=255)
                       gATGGCGCTGCACGCAGAGGCACAGCGCCACCGGGTATACAGATTGTTGACAAAGTGCGCTTTGTTTATGCCGGATGCGGCGTTAACGCCTTATCCGGcc                                                                                      <  2:4139052/100‑1 (MQ=255)
                       gATGGCGCTGCACGCAGAGGCACAGCGCCACCGGGTATACAGATTGTTGACAAAGTGCGCTTTGTTTATGCCGGATAAGGCGTTAACGCCGCATCCGGca                                                                                      <  1:2237018/100‑2 (MQ=37)
                       gATGGCGCTGCACGCAGAGGCACAGCGCCACCGGGTATACAGATTGTTGACAAAGTGCGCTTTCTTTATGCCGGATGCGGCGTTAACGCCTTATCCGGcc                                                                                      <  2:3113597/100‑1 (MQ=255)
                        aTGGCGCTGCACGCAGTGGCACAGCGCCACCGGGTATACAGATTGTTGACAAAGTGCGCTTTGTTTATGCCGGATAAGGCGTTAACGCCGCATCCGGcat                                                                                     <  2:2171244/100‑3 (MQ=25)
                        aTGGCGCTGCACGCAGAGGCACAGCGCCACCGGGTATACAGATTGTTGACAAAGTGCGCTTTGTTTATGCCGGATGCGGCGTTAACGCCTTATCCGGCCt                                                                                     <  2:2704199/100‑1 (MQ=255)
                            cgcTGCACGCAGAGGCACAGCGCCACCGGGTATACAGATTGTTGACAAAGTGCGCTTTGTTTATGCCGGATGCGGCGTTAACGCCTTATCCGGCCTACaa                                                                                >  1:2433928/1‑100 (MQ=255)
                             gcTGCACGCAGAGGCACAGCGCCACCGGGTATACAGATTGTTGACAAAGTGCGCTTTGTTTATGCCGGATGCGGCGTTAACGCCTTATCCGGCCTACaaa                                                                               >  2:836792/1‑100 (MQ=255)
                              atGCACGCAGAGGCACAGCGCCCCCGGGTATAGAGGTTGTTGACAAAGGGCGCTTTGTTTATGACGGATTCGGCGTTAACGCCTTATCCGGCCTACaaaa                                                                              >  1:3377083/2‑100 (MQ=21)
                                gCACGCAGAGGCACAGCGCCACCGGGTATACAGATTGTTGACAAAGTGCGCTTTGTTTATGCCGGATAAGGCGTTAACGCCGCATCCGGcataaacaaag                                                                             <  1:125109/100‑11 (MQ=21)
                                 cACGCAGAGGCACAGCGCCACCGGGTATACAGATTGTTGACAAAGTGCGCTTTGTTTATGCCGGATGCGGCGTTAACGCCTTATCCGGCCTACAAAATCg                                                                           >  2:3239664/1‑100 (MQ=255)
                                   cgatctggCACAGCGCCACCGGGTATACAGATTGTTGACAAAGTGCGCTTTGTTTATGCCGGATGCGGCGTTAACGCCTTATCCGGCCTACAAAATCGTg                                                                         <  1:556371/94‑1 (MQ=255)
                                    gCAGAGGCACAGCGCCACCGGGTATACAGATTGTTGACAAAGTGCGCTTTGTTTATGCCGGATGCGGCGTTAACGCCTTATCCGGCCTACAAAATCGTGc                                                                        <  2:1562290/100‑1 (MQ=255)
                                    gCAGAGGCACAGCGCCACCGGGTATACAGATTGTTGACAAAGTGCGCTTTGTTTATGCCGGATGCGGCGTTAACGCCTTATCCGGCCTACAAAATCGTGc                                                                        <  2:1161768/100‑1 (MQ=255)
                                     cAGAGGCACAGCGCCACCGGGTATACAGATTGTTGACAAAGTGCGCTTTGTTTATGCCGGATGCGGCGTTAACGCCTTATCCGGCCTACAAAATCGTGCa                                                                       >  1:2284197/1‑100 (MQ=255)
                                     cAGAGGCACAGCGCCACCGGGTATACAGATTGTTGACAAAGTGCGCTTTGTTTATGCCGGATGCGGCGTTAACGCCTTATCCGGCCTACAAAATCGTGCa                                                                       <  2:1373877/100‑1 (MQ=255)
                                     cAGAGGCACAGCGCCACCGGGTATACAGATTGTTGACAAAGTGCGCTTTGTTTATGCCGGATGCGGCGTTAACGCCTTATCCGGCCTACAAAATCGTGCa                                                                       >  2:2648129/1‑100 (MQ=255)
                                     cAGAGGCACAGCGCCACCGGGTATACAGATTGTTGACAAAGTGCGCTTTGTTTATGCCGGATGCGGCGTTAACGCCTTATCCGGCCTACAAAATCGTGCa                                                                       <  1:3809608/100‑1 (MQ=255)
                                          gCACAGCGCCACCGGGTATACAGATTGTTGACAAAGTGCGCTTTGTTTATGCCGGATGCGGCGTTAACGCCTTATCCGGCCTACAAAATCGTGCAAATTc                                                                  <  2:2374408/100‑1 (MQ=255)
                                           cacaGCGCCACCGGGTATACAGATTGTTGACAAAGTGCGCTTTGTTTATGCCGGATGCGGCGTTAACGCCTTATCCGGCCTACAAAATCGTGCAAATTCa                                                                 >  1:1063379/1‑100 (MQ=255)
                                           cacaGCGCCACCGGGTATACAGATTGTTGACAAAGTGCGCTTTGTTTATGCCGGATGCGGCGTTAACGCCTTATCCGGCCTACAAAATCGTGCAAATTCa                                                                 >  1:3290612/1‑100 (MQ=255)
                                             caGCGCCACCGGGTATACAGATTGTTGACAAAGTGCGCTTTGTTTATGCCGGATGCGGCGTTAACGCCTTATCCGGCCTACAACATCGTGCAAATTCAAt                                                               >  1:202565/1‑100 (MQ=255)
                                             caGCGCCACCGGGTATACAGATTGTTGACAAAGTGCGCTTTGTTTATGCCGGATGCGGCGTTAACGCCTTATCCGGCCTACAAAATCGTGCAAATTCAAt                                                               <  2:1412288/100‑1 (MQ=255)
                                              aGCGCCACCGGATATACAGATTGTTGACAAAGTGCGCTTTGTTTATGCCGGATGCGGCGTTAACGCCTTATCCGGCCTACAAAATCGTGCAAATTCAAt                                                               <  1:1886058/99‑1 (MQ=255)
                                                  ccACCGGGTATACAGATTGTTGACAAAGTGCGCTTTGTTTATGCCGGATGCGGCGTTAACGCCTTATCCGGCCTACAAAATCGTGCAAATTCAATACAtt                                                          >  2:2841549/1‑100 (MQ=255)
                                                  ccACCGGGTATACAGATTGTTGACAAAGTGCGCTTTGTTTATGCCGGATGCGGCGTTAACGCCTTATCCGGCCTACAAAATCGTGCAAATTCAATACAtt                                                          >  2:3938411/1‑100 (MQ=255)
                                                    aCCGGGTATACAGATTGTTGACAAAGTGCGCTTTGTTTATGCCGGATGCGGCGTTAACGCCTTATCCGGCCTACAAAATCGTGCAAATTCAATACATTGc                                                        <  1:2591134/100‑1 (MQ=255)
                                                      cGGGTATACAGATTGTTGACAAAGTGCGCTTTGTTTATGCCGGATGCGGCGTTAACGCCTTATCCGGCCTACAAAATCGTGCAAATTCAATACATTGCa                                                       <  1:169121/99‑1 (MQ=255)
                                                      cGGGTATACAGATTGTTGACAAAGTGCGCTTTGTTTATGCCGGATGCGGCGTTAACGCCTTATCCGGCCTACAAAATCGTGCAAATTCAATACATTGCAg                                                      <  2:378939/100‑1 (MQ=255)
                                                       gagTATACAGATTGTTGACAAAGTTCGCTTTGTTTATGCCGGATGCGGCGTTAACGCCTTATCCGGCCTACAAAATCGTGCAAATTCAATACATTGCAG‑‑‑a                                                  <  1:405704/98‑1 (MQ=255)
                                                          tataCAGATTGTTGACAAAGTGCGCTTTGTTTATGCCCGATTCGGCGTTAACGACTTATCCGGCCTACAAAATCGTGCAGATTCAATACATTGCAG‑‑‑AAtt                                             >  1:1484584/1‑100 (MQ=18)
                                                            taCAGATTGTTGACAAAGTGCGCTTTGTTTATGCCGGATGCGGCGTTAACGCCTTATCCGGCCTACAAAATCGTGCAAATTCAATACATTGCAG‑‑‑AAtttt                                           >  2:1726913/1‑100 (MQ=255)
                                                            taCAGATTGTTGACAAAGAGCGCTGAGTTTATGCCGGATGCGGCGGTAACGCCTTATCCGGCCTACAACATCGTGCAAATTCAATACATTGAAG‑‑‑AAtttt                                           >  1:3845306/1‑100 (MQ=17)
                                                             aCAGATTGTTGACAAAGTGCGCTTTGTTTATGCCGGATGCGGCGTTAACGCCTTATCCGGCCTACAAAATCGTGCAAATTCAATACATTGCAG‑‑‑AATTTTc                                          >  2:1940308/1‑100 (MQ=255)
                                                              cAGATTGTTGACAAAGTGCGCTTTGTTTATGCCGGATGCGGCGTTAACGCCTTATCCGGCCTACAAAATCGTGCAAATTCAATACATTGCAG‑‑‑AATTTTC‑‑g                                       <  1:3501177/100‑1 (MQ=255)
                                                              cAGATTGTTGACAAAGTGCGCTTTGTTTATGCCGGATGCGGCGTTAACGCCTTATCCGGCCTACAAAATCGTGCAAATTCAATACATTGCAG‑‑‑AATTTTC‑‑g                                       <  2:2995149/100‑1 (MQ=255)
                                                               aGATTGTTGACAAAGTGCGCTTTGTTTATGCCGGATGCGGCGTTAACGCCTTATCCGACCTACAACATCGTGCAAGTTTAATCCATTGCAG‑‑‑AATTTTC‑‑Gt                                      >  2:467455/1‑100 (MQ=18)
                                                                gATTGTTGACAAAGTGCGCTTTGTTTATGCCGGATGCGGCGTTAACGCCTTATCCGGCCTACAAAATCGTGCAAATTCAATACATTGCAG‑‑‑AATTTTC‑‑GTa                                     <  1:3150194/100‑1 (MQ=255)
                                                                  ttgttgACAAAGTGCGCTTTGTTTATGCCGGATGCGGCGTTAACGCCTTATCCGGCCTACAAAATCGTGCAAATTCAATACATTGCAG‑‑‑AATTTTC‑‑GTAgg                                   <  1:247495/100‑1 (MQ=255)
                                                                    gttgACAAAGTGCGCTTTGTTTATGCCGGATGCGGCGTTAACGCCTTATCCGGCCTAGAAAATCGTGCAAATTCAATACATTGCAG‑‑‑AATTTTC‑‑GTAGGcc                                 >  1:3088648/1‑100 (MQ=255)
                                                                    gttgACAAAGTGCGCTTTGTTTATGCCGGATGCGGCGTTAACGCCTTATCCGGCCTACAAAATCGTGCAAATTCAATACATTGCAG‑‑‑AATTTTC‑‑GTAGGcc                                 >  2:3589071/1‑100 (MQ=255)
                                                                    gttgACAAAGTGCGCTTTGTTTATGCCGGATGCGGCGTTAACGCCTTATCCGGCCTACAAAATCGTGCAAATTCAATACATTGCAG‑‑‑AATTTTC‑‑GTAGGcc                                 <  2:715622/100‑1 (MQ=255)
                                                                       gACAAAGTGCGCTTTGTTTATGCCGGATGCGGCGTTAACGCCTTATCCGGCCTACAAAATCGTGCAAATTCAATACATTGCAG‑‑‑AATTTTC‑‑GTAGGCCTGa                              >  1:3744466/1‑100 (MQ=255)
                                                                              tGCGCTTTGTTTATGCCGGATGCGGCGTTAACGCCTTATCCGGCCTACAAAATCGTGCAAATTCAATACATTGCAG‑‑‑AATTTTC‑‑GTAGGCCTGATAAGCgg                       <  1:2205020/100‑1 (MQ=255)
                                                                              ngcgcTTTGTTTATGCCGGATGCGGCGTTAACGCCTTATCCGGCCTACAAAATCGTGCAAATTCAATACATTGCAG‑‑‑AATTTTC‑‑GTAGGCCTGATAAGCgg                       >  1:3556827/2‑100 (MQ=255)
                                                                               gcgcTTTGTTTATGCCGGATGCGGCGTTAACGCCTTATCCGGCCTACAAAATCGTGCAAATTCAATACATTGCAG‑‑‑AATTTTC‑‑GTAGGCCTGATAAGCGGa                      <  2:2691949/100‑1 (MQ=255)
                                                                                   tttttgtaggCCGGATAAGGCGTTCACGCCGCATCCGGCCTACAAAAATCTTGCCAATTCAATATATTGCAG‑‑‑GACCAT‑‑GTAGGCCTGATAAGCGTAgcgc                  >  2:209628/10‑100 (MQ=18)
                                                                                    ttGTTTATGCCGGATGCGGCGTTAACGCCTTATCCGGCCTACAAAATCGTGCAAATTCAATACATTGCAG‑‑‑AATTTTC‑‑GTAGGCCTGATAAGCGGAGCGCa                 >  1:737346/1‑100 (MQ=255)
                                                                                     tGTTTATGCCGGATGCGGCGTTAACGCCTTATCCGGCCTACAAAATCGTGCAAATTCAATACATTGCAG‑‑‑AATTTTC‑‑GTAGGCCTGATAAGCGGAGCGCAt                >  2:3772392/1‑100 (MQ=255)
                                                                                     cttgtaggCCGGATAAGGCGTTCACGCCGCATCCGGCCTACAAGATCGTGCAAATTCAATATATTGCAT‑‑‑GAGCCAC‑‑GTAGGCCTGATAAGCGTAgcgc                  >  2:945490/8‑98 (MQ=17)
                                                                                      ttgtaggCCGGATAAGGCGTTTACGCCGCATCCGGCCTACAAAATCGTGCAAATTCAACATATTGC‑‑‑‑‑AATTCTCTTGTAGGCCTGATAAGCGTAGCGCATc               >  2:4343869/7‑100 (MQ=18)
                                                                                      gTTTATGCCGGATGCGGCGTTAACGCCTTATCCGGCCTACAAAATCGTGCAAATTCAATACATTGCAG‑‑‑AATTTTC‑‑GTAGGCCTGATAAGCGGAGCGCATc               >  1:141880/1‑100 (MQ=255)
                                                                                       tgtaggCCGGATAAGGCGTTAACGCCGCATCCGGCCTACAAAATCGTGCAAATTCAATACATTGCAG‑‑‑AATTTTC‑‑GTAGGCCTGATAAGCGGAGCGCATc               >  2:2264423/6‑99 (MQ=17)
                                                                                        gtaggCCGGATAAGGCGTTCACGCTGCATCCGGCCTACAAAATCGCTTAAATTCAATATATTGCAC‑‑‑CATCCTC‑‑GTAGGCCTGATAAGCGTAGCGCATCAg             >  2:645012/5‑100 (MQ=18)
                                                                                        gtaggCCGGATAAGGCGTTCACGCCGCATCCGGCCTACAAAAACCTGCAAATTCAATATATTGCAG‑‑‑GATACGT‑‑GTAGGCCTGATAAGCGTAGCGCATCAg             >  2:2582802/5‑100 (MQ=18)
                                                                                        gtaggCCGGATAAGGCGTTAACGCCGCATCCGGCCTACAAAATCGTGCAAATGCAATACATTGCAG‑‑‑AAGTTTC‑‑GTAGACCTGATAAGCGGAGCGCATCAg             >  1:2980673/5‑100 (MQ=18)
                                                                                         taggCCGGATAAGGCGTTCACGCCGCATCCGGCCTACAAAACTTTGCAAATTCAATATATTGC‑‑‑‑‑AATCTCCGTGTAGGCCTGATAAGCGTAGCGCATCAgg            >  1:3253062/4‑100 (MQ=18)
                                                                                         taggCCGGATAAGGCGTTCACGCCGCATCCGGCCTACAAAACTGTGCAAATTCAATATATTGCAG‑‑‑GAAACAC‑‑GTAGGCCTGATAAGCGAAGCGCATCAgg            >  1:2266878/4‑100 (MQ=18)
                                                                                         taggCCGGATAAGGCGTTCACGCCGCATCCGGCCTACAAAAACATGCAAATTCAATAGATTGCAG‑‑‑AGATTAT‑‑GTAGGCCTGATAAGCGTAGCGCATCAgg            >  2:2553222/4‑100 (MQ=17)
                                                                                         tATGCCGGATGCGGCGTTAACGCCTTATCCGGCCTACAAAATCGTGCAAATTCAATACATTGCAG‑‑‑AATTTTC‑‑GTAGGCCTGATAAGCGGAGCGCATCAgg            >  1:3266858/1‑100 (MQ=255)
                                                                                          aggCCGGATAAGGCGTTTACGCCGCATCCGGCCTACAAAAATGTGCAAATTCAATAAATTGCAATTCAACTT‑‑‑‑GTAGGCCTGATAAGCGCAGCGCATCAgg            >  2:4188952/3‑100 (MQ=21)
                                                                                          aTGCCGGATGCGGCGTTAACGCCTTATCCGGCCTACAAAATCGTGCAAATTCAATACATTGCAG‑‑‑AATTTTC‑‑GTAGGCCTGATAAGCGGAGCGCATCAGGc           <  2:1966884/100‑1 (MQ=255)
                                                                                          aTGCCGGATGCGGCGTTAACGCCTTATCCGGCCTACAAAATCGTGCAAATTCAATACATTGCAG‑‑‑AATTTTC‑‑GTAGGCCTGATAAGCGGAGCGCATCAGGc           >  2:380383/1‑100 (MQ=255)
                                                                                           tGCCGGATGCGGCGTTAACGCCTTATCCGGCCTACAAAATCGTGCAAATTCAATACATTGCAG‑‑‑AATTTTC‑‑GTAGGCCTGATAAGCGGAGCGCATCAGGCa          >  2:685377/1‑100 (MQ=255)
                                                                                           ggCCGGATAAGGCGTTTACGCCGCATCCGGCCTACAAATTCTGCAAATTCAACAAATTGCG‑‑‑AATCCCTT‑‑GTAGGCCTGATAAGCGTAGCGCATCAGGCAt         >  2:4283658/2‑100 (MQ=18)
                                                                                           ggCCGGATAAGGCGTTTACGCCGCATCCGGCCTACAAAATCGTGCAAATTCAACATATTGCAA‑‑‑TTCTCTT‑‑GTAGGCCTGATAAGCGTAGCGCATCAGGct          >  2:868667/2‑99 (MQ=18)
                                                                                           ggCCGGATAAGGCGTTTACGCCGCATCCGGCCTACAAAAGCATGCAAATTCAATATATTGCAG‑‑‑AGTTCAC‑‑GTAGGCCTGATAAGCGTAGCGCATCAGGCa          >  1:1609293/2‑100 (MQ=14)
                                                                                           ggCCGGATAAGGCGTTCACGCCGCATCCGGCCTACGAAATCGTGCTAATTCAATATATTGTGG‑‑‑GCTCCAC‑‑GTAGGCCTGATAAGCGTAGCGCATCAGGCa          >  2:2075891/2‑100 (MQ=18)
                                                                                           ggCCGGATAAGGCGTTCACGCCGCATCCGGCCTACAAAATCATGAAAATTCAATACATTGCAA‑‑‑GATTTTC‑‑GTAGGCCTGATAAGCGTAGCGCATCAGGCa          >  1:1361450/2‑100 (MQ=17)
                                                                                           ggCCGGATAAGGCGTTAACGCCGCATCCGGCCTACAAAATCGTGCAAATTCAATACATTGCAG‑‑‑AATTTTC‑‑GTAGGCCTGATAAGCGGAGCGCATCAGGCa          >  1:1816571/2‑100 (MQ=18)
                                                                                            gCCGGATGCGGCGTTAACGCCTTATCCGGCCTACAAAATCGTGCAAATTCAATACATTGCAG‑‑‑AATTTTC‑‑GTAGGCCTGATAAGCGGAGCGCATCAGGCAt         >  1:944810/1‑100 (MQ=255)
                                                                                            gCCGGATAAGGCGTTTACGCCGCATCCGGCCTACAAATTCTGCAAATTCAACAAATTGCG‑‑‑AATCCCTT‑‑GTAGGCCTGATAAGCGTAGCGCATCAGGCAtt        >  2:2747060/1‑100 (MQ=18)
                                                                                            gCCGGATAAGGCGTTTACGCCGCATCCGGCCTACAAAATCGTGCAAATTCAACATATTGC‑‑‑‑‑AATTCTCTTGTAGGCCTGATAAGCGTAGCGCATCAGGctg         >  1:4108919/1‑98 (MQ=18)
                                                                                            gCCGGATAAGGCGTTCACGCCGCATCCGGCCTACAAGATCGTGCAAATTCAATATATTGCAT‑‑‑GAGCCAC‑‑GTAGTCCTGATAAGCGTAGCGCATCAGGCag         >  1:1008699/1‑99 (MQ=17)
                                                                                            gCCGGATAAGGCGTTCACGCCGCATCCGGCCTACAAGATCGTGCAAATTCAATATATTGCAT‑‑‑GAGCCAC‑‑GTAGGCCTGATAAGCGTAGCGCATCAGGCag         >  2:1820155/1‑99 (MQ=17)
                                                                                            gCCGGATAAGGCGTTCACGCCGCATCCGGCCTACAAGATCGTGCAAATTCAATATATTGCAT‑‑‑GAGCCAC‑‑GTAGGCCTGATAAGCGTAGCGCATCAGACa          >  2:1010024/1‑99 (MQ=17)
                                                                                            gCCGGATAAGGCGTTCACGCCGCATCCGGCCTACAAAAGCATGAAAATTCAATATATTGCAG‑‑‑GAGCTGC‑‑GTAGGCCTGATAAGCGTAGCGCATCAGGCag         >  2:2227271/1‑99 (MQ=18)
                                                                                            gCCGGATAAGGCGTTCACGCCGCATCCGGCCTACAAAAACATGCAAATTCAATAGATTGCAG‑‑‑AGATTAT‑‑GTAGGCCTGATAAGCGTAGCGCATCAGGCaa         >  2:3361635/1‑99 (MQ=17)
                                                                                            gCCGGATAAGGCGTTCACGCCGCATCCGGCCTACAAATCATGCAAATCCAAAATATTGCAG‑‑‑GAGTTAT‑‑GTAGGCCTGATAAGCGTAGCGCATCAGGCAtt        >  1:1330197/1‑100 (MQ=21)
                                                                                             ccGGATAAGGCGGTCACTCCGCATCCGACCTACAAATCATGCAAATCCAAAATATTGCAG‑‑‑GAGGTAT‑‑GTATGCCTGATAAGCGTAGCGCATCAGGCAtt        >  1:4302543/1‑99 (MQ=18)
                                                                                                gATGCGGCGTTAACGCCTTATCCGGCCTACAAAATCGTGCAAATTCAATACATTGCAG‑‑‑AATTTTC‑‑GTAGGCCTGATAAGCGGAGCGCATCAGGCAttttt     <  1:1232349/100‑1 (MQ=255)
                                                                                                 aTGCGGCGTTAACGCCTTATCCGGCCTACAAAATCGTGCAAATTCAATACATTGCAG‑‑‑AATTTTC‑‑GTAGGCCTGATAAGCGGAGCGCATCAGGCATTtttg    <  2:2461829/100‑1 (MQ=255)
                                                                                                  tGCGGCGTTAACGCCTTATCCGGCCTACAAAATCGTGCAAATTCAATACATTGCAG‑‑‑AATTTTC‑‑GTAGGCCTGATAAGCGGAGCGCATCAGGCATTtttgc   >  2:2022169/1‑100 (MQ=255)
                                                                                                   |                                                                                                           
CAGGGAATGAAGTATTGCTGCCGGATGGCGCTGCACGCAGAGGCACAGCGCCACCGGGTATACAGATTGTTGACAAAGTGCGCTTTGTTTATGCCGGATGCGGCGTTAACGCCTTATCCGGCCTACAAAATCGTGCAAATTCAATACATTGCAG‑‑‑AATTTTC‑‑GTAGGCCTGATAAGCGGAGCGCATCAGGCATTTTTGCC  >  NC_000913/3705871‑3706069

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: