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.

CAGGGAATGAAGTATTGCTGCCGGATGGCGCTGCACGCAGAGGCACAGCGCCACCGGGTATACAGATTGTTGACAAAGTGCGCTTTGTTTATGCCGGATGCGGCGTTAACGCCTTATCCGGCCTAC‑AAAATCGTGCAAATTCAATACATTGCAG‑‑‑A‑AT‑TTTC‑‑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)
                            cgcTGCACGCAGAGGCACAGCGCCACCGGGTATACAGATTGTTGACAAAGTGCGCTTTGTTTATGCCGGATGCGGCGTTAACGCCTTATCCGGCCTAC‑aa                                                                                >  1:2433928/1‑100 (MQ=255)
                             gcTGCACGCAGAGGCACAGCGCCACCGGGTATACAGATTGTTGACAAAGTGCGCTTTGTTTATGCCGGATGCGGCGTTAACGCCTTATCCGGCCTAC‑aaa                                                                               >  2:836792/1‑100 (MQ=255)
                              atGCACGCAGAGGCACAGCGCCCCCGGGTATAGAGGTTGTTGACAAAGGGCGCTTTGTTTATGACGGATTCGGCGTTAACGCCTTATCCGGCCTAC‑aaaa                                                                              >  1:3377083/2‑100 (MQ=21)
                                gCACGCAGAGGCACAGCGCCACCGGGTATACAGATTGTTGACAAAGTGCGCTTTGTTTATGCCGGATAAGGCGTTAACGCCGCATCCGGcataaacaaag                                                                             <  1:125109/100‑11 (MQ=21)
                                 cACGCAGAGGCACAGCGCCACCGGGTATACAGATTGTTGACAAAGTGCGCTTTGTTTATGCCGGATGCGGCGTTAACGCCTTATCCGGCCTAC‑AAAATCg                                                                           >  2:3239664/1‑100 (MQ=255)
                                   cgatctggCACAGCGCCACCGGGTATACAGATTGTTGACAAAGTGCGCTTTGTTTATGCCGGATGCGGCGTTAACGCCTTATCCGGCCTAC‑AAAATCGTg                                                                         <  1:556371/94‑1 (MQ=255)
                                    gCAGAGGCACAGCGCCACCGGGTATACAGATTGTTGACAAAGTGCGCTTTGTTTATGCCGGATGCGGCGTTAACGCCTTATCCGGCCTAC‑AAAATCGTGc                                                                        <  2:1562290/100‑1 (MQ=255)
                                    gCAGAGGCACAGCGCCACCGGGTATACAGATTGTTGACAAAGTGCGCTTTGTTTATGCCGGATGCGGCGTTAACGCCTTATCCGGCCTAC‑AAAATCGTGc                                                                        <  2:1161768/100‑1 (MQ=255)
                                     cAGAGGCACAGCGCCACCGGGTATACAGATTGTTGACAAAGTGCGCTTTGTTTATGCCGGATGCGGCGTTAACGCCTTATCCGGCCTAC‑AAAATCGTGCa                                                                       >  1:2284197/1‑100 (MQ=255)
                                     cAGAGGCACAGCGCCACCGGGTATACAGATTGTTGACAAAGTGCGCTTTGTTTATGCCGGATGCGGCGTTAACGCCTTATCCGGCCTAC‑AAAATCGTGCa                                                                       <  2:1373877/100‑1 (MQ=255)
                                     cAGAGGCACAGCGCCACCGGGTATACAGATTGTTGACAAAGTGCGCTTTGTTTATGCCGGATGCGGCGTTAACGCCTTATCCGGCCTAC‑AAAATCGTGCa                                                                       >  2:2648129/1‑100 (MQ=255)
                                     cAGAGGCACAGCGCCACCGGGTATACAGATTGTTGACAAAGTGCGCTTTGTTTATGCCGGATGCGGCGTTAACGCCTTATCCGGCCTAC‑AAAATCGTGCa                                                                       <  1:3809608/100‑1 (MQ=255)
                                          gCACAGCGCCACCGGGTATACAGATTGTTGACAAAGTGCGCTTTGTTTATGCCGGATGCGGCGTTAACGCCTTATCCGGCCTAC‑AAAATCGTGCAAATTc                                                                  <  2:2374408/100‑1 (MQ=255)
                                           cacaGCGCCACCGGGTATACAGATTGTTGACAAAGTGCGCTTTGTTTATGCCGGATGCGGCGTTAACGCCTTATCCGGCCTAC‑AAAATCGTGCAAATTCa                                                                 >  1:1063379/1‑100 (MQ=255)
                                           cacaGCGCCACCGGGTATACAGATTGTTGACAAAGTGCGCTTTGTTTATGCCGGATGCGGCGTTAACGCCTTATCCGGCCTAC‑AAAATCGTGCAAATTCa                                                                 >  1:3290612/1‑100 (MQ=255)
                                             caGCGCCACCGGGTATACAGATTGTTGACAAAGTGCGCTTTGTTTATGCCGGATGCGGCGTTAACGCCTTATCCGGCCTAC‑AACATCGTGCAAATTCAAt                                                               >  1:202565/1‑100 (MQ=255)
                                             caGCGCCACCGGGTATACAGATTGTTGACAAAGTGCGCTTTGTTTATGCCGGATGCGGCGTTAACGCCTTATCCGGCCTAC‑AAAATCGTGCAAATTCAAt                                                               <  2:1412288/100‑1 (MQ=255)
                                              aGCGCCACCGGATATACAGATTGTTGACAAAGTGCGCTTTGTTTATGCCGGATGCGGCGTTAACGCCTTATCCGGCCTAC‑AAAATCGTGCAAATTCAAt                                                               <  1:1886058/99‑1 (MQ=255)
                                                  ccACCGGGTATACAGATTGTTGACAAAGTGCGCTTTGTTTATGCCGGATGCGGCGTTAACGCCTTATCCGGCCTAC‑AAAATCGTGCAAATTCAATACAtt                                                          >  2:2841549/1‑100 (MQ=255)
                                                  ccACCGGGTATACAGATTGTTGACAAAGTGCGCTTTGTTTATGCCGGATGCGGCGTTAACGCCTTATCCGGCCTAC‑AAAATCGTGCAAATTCAATACAtt                                                          >  2:3938411/1‑100 (MQ=255)
                                                    aCCGGGTATACAGATTGTTGACAAAGTGCGCTTTGTTTATGCCGGATGCGGCGTTAACGCCTTATCCGGCCTAC‑AAAATCGTGCAAATTCAATACATTGc                                                        <  1:2591134/100‑1 (MQ=255)
                                                      cGGGTATACAGATTGTTGACAAAGTGCGCTTTGTTTATGCCGGATGCGGCGTTAACGCCTTATCCGGCCTAC‑AAAATCGTGCAAATTCAATACATTGCa                                                       <  1:169121/99‑1 (MQ=255)
                                                      cGGGTATACAGATTGTTGACAAAGTGCGCTTTGTTTATGCCGGATGCGGCGTTAACGCCTTATCCGGCCTAC‑AAAATCGTGCAAATTCAATACATTGCAg                                                      <  2:378939/100‑1 (MQ=255)
                                                       gagTATACAGATTGTTGACAAAGTTCGCTTTGTTTATGCCGGATGCGGCGTTAACGCCTTATCCGGCCTAC‑AAAATCGTGCAAATTCAATACATTGCAG‑‑‑a                                                  <  1:405704/98‑1 (MQ=255)
                                                          tataCAGATTGTTGACAAAGTGCGCTTTGTTTATGCCCGATTCGGCGTTAACGACTTATCCGGCCTAC‑AAAATCGTGCAGATTCAATACATTGCAG‑‑‑A‑At‑t                                             >  1:1484584/1‑100 (MQ=18)
                                                            taCAGATTGTTGACAAAGTGCGCTTTGTTTATGCCGGATGCGGCGTTAACGCCTTATCCGGCCTAC‑AAAATCGTGCAAATTCAATACATTGCAG‑‑‑A‑At‑ttt                                           >  2:1726913/1‑100 (MQ=255)
                                                            taCAGATTGTTGACAAAGAGCGCTGAGTTTATGCCGGATGCGGCGGTAACGCCTTATCCGGCCTAC‑AACATCGTGCAAATTCAATACATTGAAG‑‑‑A‑At‑ttt                                           >  1:3845306/1‑100 (MQ=17)
                                                             aCAGATTGTTGACAAAGTGCGCTTTGTTTATGCCGGATGCGGCGTTAACGCCTTATCCGGCCTAC‑AAAATCGTGCAAATTCAATACATTGCAG‑‑‑A‑AT‑TTTc                                          >  2:1940308/1‑100 (MQ=255)
                                                              cAGATTGTTGACAAAGTGCGCTTTGTTTATGCCGGATGCGGCGTTAACGCCTTATCCGGCCTAC‑AAAATCGTGCAAATTCAATACATTGCAG‑‑‑A‑AT‑TTTC‑‑g                                       <  1:3501177/100‑1 (MQ=255)
                                                              cAGATTGTTGACAAAGTGCGCTTTGTTTATGCCGGATGCGGCGTTAACGCCTTATCCGGCCTAC‑AAAATCGTGCAAATTCAATACATTGCAG‑‑‑A‑AT‑TTTC‑‑g                                       <  2:2995149/100‑1 (MQ=255)
                                                               aGATTGTTGACAAAGTGCGCTTTGTTTATGCCGGATGCGGCGTTAACGCCTTATCCGACCTAC‑AACATCGTGCAAGTTTAATCCATTGCAG‑‑‑A‑AT‑TTTC‑‑Gt                                      >  2:467455/1‑100 (MQ=18)
                                                                gATTGTTGACAAAGTGCGCTTTGTTTATGCCGGATGCGGCGTTAACGCCTTATCCGGCCTAC‑AAAATCGTGCAAATTCAATACATTGCAG‑‑‑A‑AT‑TTTC‑‑GTa                                     <  1:3150194/100‑1 (MQ=255)
                                                                  ttgttgACAAAGTGCGCTTTGTTTATGCCGGATGCGGCGTTAACGCCTTATCCGGCCTAC‑AAAATCGTGCAAATTCAATACATTGCAG‑‑‑A‑AT‑TTTC‑‑GTAgg                                   <  1:247495/100‑1 (MQ=255)
                                                                    gttgACAAAGTGCGCTTTGTTTATGCCGGATGCGGCGTTAACGCCTTATCCGGCCTAG‑AAAATCGTGCAAATTCAATACATTGCAG‑‑‑A‑AT‑TTTC‑‑GTAGGcc                                 >  1:3088648/1‑100 (MQ=255)
                                                                    gttgACAAAGTGCGCTTTGTTTATGCCGGATGCGGCGTTAACGCCTTATCCGGCCTAC‑AAAATCGTGCAAATTCAATACATTGCAG‑‑‑A‑AT‑TTTC‑‑GTAGGcc                                 >  2:3589071/1‑100 (MQ=255)
                                                                    gttgACAAAGTGCGCTTTGTTTATGCCGGATGCGGCGTTAACGCCTTATCCGGCCTAC‑AAAATCGTGCAAATTCAATACATTGCAG‑‑‑A‑AT‑TTTC‑‑GTAGGcc                                 <  2:715622/100‑1 (MQ=255)
                                                                       gACAAAGTGCGCTTTGTTTATGCCGGATGCGGCGTTAACGCCTTATCCGGCCTAC‑AAAATCGTGCAAATTCAATACATTGCAG‑‑‑A‑AT‑TTTC‑‑GTAGGCCTGa                              >  1:3744466/1‑100 (MQ=255)
                                                                              tGCGCTTTGTTTATGCCGGATGCGGCGTTAACGCCTTATCCGGCCTAC‑AAAATCGTGCAAATTCAATACATTGCAG‑‑‑A‑AT‑TTTC‑‑GTAGGCCTGATAAGCgg                       <  1:2205020/100‑1 (MQ=255)
                                                                              ngcgcTTTGTTTATGCCGGATGCGGCGTTAACGCCTTATCCGGCCTAC‑AAAATCGTGCAAATTCAATACATTGCAG‑‑‑A‑AT‑TTTC‑‑GTAGGCCTGATAAGCgg                       >  1:3556827/2‑100 (MQ=255)
                                                                               gcgcTTTGTTTATGCCGGATGCGGCGTTAACGCCTTATCCGGCCTAC‑AAAATCGTGCAAATTCAATACATTGCAG‑‑‑A‑AT‑TTTC‑‑GTAGGCCTGATAAGCGGa                      <  2:2691949/100‑1 (MQ=255)
                                                                                   tttttgtaggCCGGATAAGGCGTTCACGCCGCATCCGGCCTACAAAAATCTTGCCAATTCAATATATTGCAG‑‑‑G‑AC‑CAT‑‑‑GTAGGCCTGATAAGCGTAgcgc                  >  2:209628/10‑100 (MQ=18)
                                                                                    ttGTTTATGCCGGATGCGGCGTTAACGCCTTATCCGGCCTAC‑AAAATCGTGCAAATTCAATACATTGCAG‑‑‑A‑AT‑TTTC‑‑GTAGGCCTGATAAGCGGAGCGCa                 >  1:737346/1‑100 (MQ=255)
                                                                                     tGTTTATGCCGGATGCGGCGTTAACGCCTTATCCGGCCTAC‑AAAATCGTGCAAATTCAATACATTGCAG‑‑‑A‑AT‑TTTC‑‑GTAGGCCTGATAAGCGGAGCGCAt                >  2:3772392/1‑100 (MQ=255)
                                                                                     cttgtaggCCGGATAAGGCGTTCACGCCGCATCCGGCCTAC‑AAGATCGTGCAAATTCAATATATTGCAT‑‑‑G‑AG‑CCAC‑‑GTAGGCCTGATAAGCGTAgcgc                  >  2:945490/8‑98 (MQ=17)
                                                                                      ttgtaggCCGGATAAGGCGTTTACGCCGCATCCGGCCTAC‑AAAATCGTGCAAATTCAACATATTGC‑‑‑‑‑A‑AT‑TCTCTTGTAGGCCTGATAAGCGTAGCGCATc               >  2:4343869/7‑100 (MQ=18)
                                                                                      gTTTATGCCGGATGCGGCGTTAACGCCTTATCCGGCCTAC‑AAAATCGTGCAAATTCAATACATTGCAG‑‑‑A‑AT‑TTTC‑‑GTAGGCCTGATAAGCGGAGCGCATc               >  1:141880/1‑100 (MQ=255)
                                                                                       tgtaggCCGGATAAGGCGTTAACGCCGCATCCGGCCTAC‑AAAATCGTGCAAATTCAATACATTGCAG‑‑‑A‑AT‑TTTC‑‑GTAGGCCTGATAAGCGGAGCGCATc               >  2:2264423/6‑99 (MQ=17)
                                                                                        gtaggCCGGATAAGGCGTTCACGCTGCATCCGGCCTAC‑AAAATCGCTTAAATTCAATATATTGCAC‑‑‑C‑AT‑CCTC‑‑GTAGGCCTGATAAGCGTAGCGCATCAg             >  2:645012/5‑100 (MQ=18)
                                                                                        gtaggCCGGATAAGGCGTTCACGCCGCATCCGGCCTAC‑AAAAACCTGCAAATTCAATATATTGCAG‑‑‑G‑AT‑ACGT‑‑GTAGGCCTGATAAGCGTAGCGCATCAg             >  2:2582802/5‑100 (MQ=18)
                                                                                        gtaggCCGGATAAGGCGTTAACGCCGCATCCGGCCTAC‑AAAATCGTGCAAATGCAATACATTGCAG‑‑‑A‑AG‑TTTC‑‑GTAGACCTGATAAGCGGAGCGCATCAg             >  1:2980673/5‑100 (MQ=18)
                                                                                         taggCCGGATAAGGCGTTCACGCCGCATCCGGCCTAC‑AAAACTTTGCAAATTCAATATATTGC‑‑‑‑‑A‑AT‑CTCCGTGTAGGCCTGATAAGCGTAGCGCATCAgg            >  1:3253062/4‑100 (MQ=18)
                                                                                         taggCCGGATAAGGCGTTCACGCCGCATCCGGCCTAC‑AAAACTGTGCAAATTCAATATATTGCAG‑‑‑G‑AA‑ACAC‑‑GTAGGCCTGATAAGCGAAGCGCATCAgg            >  1:2266878/4‑100 (MQ=18)
                                                                                         taggCCGGATAAGGCGTTCACGCCGCATCCGGCCTAC‑AAAAACATGCAAATTCAATAGATTGCAG‑‑‑A‑GA‑TTAT‑‑GTAGGCCTGATAAGCGTAGCGCATCAgg            >  2:2553222/4‑100 (MQ=17)
                                                                                         tATGCCGGATGCGGCGTTAACGCCTTATCCGGCCTAC‑AAAATCGTGCAAATTCAATACATTGCAG‑‑‑A‑AT‑TTTC‑‑GTAGGCCTGATAAGCGGAGCGCATCAgg            >  1:3266858/1‑100 (MQ=255)
                                                                                          aggCCGGATAAGGCGTTTACGCCGCATCCGGCCTAC‑AAAAATGTGCAAATTCAATAAATTGCAATTCA‑AC‑TT‑‑‑‑GTAGGCCTGATAAGCGCAGCGCATCAgg            >  2:4188952/3‑100 (MQ=21)
                                                                                          aTGCCGGATGCGGCGTTAACGCCTTATCCGGCCTAC‑AAAATCGTGCAAATTCAATACATTGCAG‑‑‑A‑AT‑TTTC‑‑GTAGGCCTGATAAGCGGAGCGCATCAGGc           <  2:1966884/100‑1 (MQ=255)
                                                                                          aTGCCGGATGCGGCGTTAACGCCTTATCCGGCCTAC‑AAAATCGTGCAAATTCAATACATTGCAG‑‑‑A‑AT‑TTTC‑‑GTAGGCCTGATAAGCGGAGCGCATCAGGc           >  2:380383/1‑100 (MQ=255)
                                                                                           tGCCGGATGCGGCGTTAACGCCTTATCCGGCCTAC‑AAAATCGTGCAAATTCAATACATTGCAG‑‑‑A‑AT‑TTTC‑‑GTAGGCCTGATAAGCGGAGCGCATCAGGCa          >  2:685377/1‑100 (MQ=255)
                                                                                           ggCCGGATAAGGCGTTTACGCCGCATCCGGCCTAC‑AAATTC‑TGCAAATTCAACAAATTGC‑G‑‑‑A‑ATCCCTT‑‑GTAGGCCTGATAAGCGTAGCGCATCAGGCAt         >  2:4283658/2‑100 (MQ=18)
                                                                                           ggCCGGATAAGGCGTTTACGCCGCATCCGGCCTAC‑AAAATCGTGCAAATTCAACATATTGCAA‑‑‑T‑TC‑TCTT‑‑GTAGGCCTGATAAGCGTAGCGCATCAGGct          >  2:868667/2‑99 (MQ=18)
                                                                                           ggCCGGATAAGGCGTTTACGCCGCATCCGGCCTAC‑AAAAGCATGCAAATTCAATATATTGCAG‑‑‑A‑GT‑TCAC‑‑GTAGGCCTGATAAGCGTAGCGCATCAGGCa          >  1:1609293/2‑100 (MQ=14)
                                                                                           ggCCGGATAAGGCGTTCACGCCGCATCCGGCCTAC‑GAAATCGTGCTAATTCAATATATTGTGG‑‑‑G‑CT‑CCAC‑‑GTAGGCCTGATAAGCGTAGCGCATCAGGCa          >  2:2075891/2‑100 (MQ=18)
                                                                                           ggCCGGATAAGGCGTTCACGCCGCATCCGGCCTAC‑AAAATCATGAAAATTCAATACATTGCAA‑‑‑G‑AT‑TTTC‑‑GTAGGCCTGATAAGCGTAGCGCATCAGGCa          >  1:1361450/2‑100 (MQ=17)
                                                                                           ggCCGGATAAGGCGTTAACGCCGCATCCGGCCTAC‑AAAATCGTGCAAATTCAATACATTGCAG‑‑‑A‑AT‑TTTC‑‑GTAGGCCTGATAAGCGGAGCGCATCAGGCa          >  1:1816571/2‑100 (MQ=18)
                                                                                            gCCGGATGCGGCGTTAACGCCTTATCCGGCCTAC‑AAAATCGTGCAAATTCAATACATTGCAG‑‑‑A‑AT‑TTTC‑‑GTAGGCCTGATAAGCGGAGCGCATCAGGCAt         >  1:944810/1‑100 (MQ=255)
                                                                                            gCCGGATAAGGCGTTTACGCCGCATCCGGCCTAC‑AAATTC‑TGCAAATTCAACAAATTGC‑G‑‑‑A‑ATCCCTT‑‑GTAGGCCTGATAAGCGTAGCGCATCAGGCAtt        >  2:2747060/1‑100 (MQ=18)
                                                                                            gCCGGATAAGGCGTTTACGCCGCATCCGGCCTAC‑AAAATCGTGCAAATTCAACATATTGC‑‑‑‑‑A‑AT‑TCTCTTGTAGGCCTGATAAGCGTAGCGCATCAGGctg         >  1:4108919/1‑98 (MQ=18)
                                                                                            gCCGGATAAGGCGTTCACGCCGCATCCGGCCTAC‑AAGATCGTGCAAATTCAATATATTGCAT‑‑‑G‑AG‑CCAC‑‑GTAGTCCTGATAAGCGTAGCGCATCAGGCag         >  1:1008699/1‑99 (MQ=17)
                                                                                            gCCGGATAAGGCGTTCACGCCGCATCCGGCCTAC‑AAGATCGTGCAAATTCAATATATTGCAT‑‑‑G‑AG‑CCAC‑‑GTAGGCCTGATAAGCGTAGCGCATCAGGCag         >  2:1820155/1‑99 (MQ=17)
                                                                                            gCCGGATAAGGCGTTCACGCCGCATCCGGCCTAC‑AAGATCGTGCAAATTCAATATATTGCAT‑‑‑G‑AG‑CCAC‑‑GTAGGCCTGATAAGCGTAGCGCATCAGACa          >  2:1010024/1‑99 (MQ=17)
                                                                                            gCCGGATAAGGCGTTCACGCCGCATCCGGCCTAC‑AAAAGCATGAAAATTCAATATATTGCAG‑‑‑G‑AG‑CTGC‑‑GTAGGCCTGATAAGCGTAGCGCATCAGGCag         >  2:2227271/1‑99 (MQ=18)
                                                                                            gCCGGATAAGGCGTTCACGCCGCATCCGGCCTAC‑AAAAACATGCAAATTCAATAGATTGCAG‑‑‑AGAT‑TAT‑‑‑GTAGGCCTGATAAGCGTAGCGCATCAGGCaa         >  2:3361635/1‑99 (MQ=17)
                                                                                            gCCGGATAAGGCGTTCACGCCGCATCCGGCCTAC‑‑AAATCATGCAAATCCAAAATATTGCAG‑‑‑G‑AG‑TTAT‑‑GTAGGCCTGATAAGCGTAGCGCATCAGGCAtt        >  1:1330197/1‑100 (MQ=21)
                                                                                             ccGGATAAGGCGGTCACTCCGCATCCGACCTAC‑‑AAATCATGCAAATCCAAAATATTGCAG‑‑‑G‑AG‑GTAT‑‑GTATGCCTGATAAGCGTAGCGCATCAGGCAtt        >  1:4302543/1‑99 (MQ=18)
                                                                                                gATGCGGCGTTAACGCCTTATCCGGCCTAC‑AAAATCGTGCAAATTCAATACATTGCAG‑‑‑A‑AT‑TTTC‑‑GTAGGCCTGATAAGCGGAGCGCATCAGGCAttttt     <  1:1232349/100‑1 (MQ=255)
                                                                                                 aTGCGGCGTTAACGCCTTATCCGGCCTAC‑AAAATCGTGCAAATTCAATACATTGCAG‑‑‑A‑AT‑TTTC‑‑GTAGGCCTGATAAGCGGAGCGCATCAGGCATTtttg    <  2:2461829/100‑1 (MQ=255)
                                                                                                  tGCGGCGTTAACGCCTTATCCGGCCTAC‑AAAATCGTGCAAATTCAATACATTGCAG‑‑‑A‑AT‑TTTC‑‑GTAGGCCTGATAAGCGGAGCGCATCAGGCATTtttgc   >  2:2022169/1‑100 (MQ=255)
                                                                                                   |                                                                                                           
CAGGGAATGAAGTATTGCTGCCGGATGGCGCTGCACGCAGAGGCACAGCGCCACCGGGTATACAGATTGTTGACAAAGTGCGCTTTGTTTATGCCGGATGCGGCGTTAACGCCTTATCCGGCCTAC‑AAAATCGTGCAAATTCAATACATTGCAG‑‑‑A‑AT‑TTTC‑‑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: ‑