Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 3,705,970 G→A 22.4% 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,9700GA22.4% 568.3 / 143.9 260intergenic (‑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 (102/98);  new base A (49/9);  total (152/108)
Fisher's exact test for biased strand distribution p-value = 3.65e-06
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.71e-02

Only 100 of 278 total aligned reads displayed.

CAGGGAATGAAGTATTGCTGCCGGATGGCGCTGCACGCAGAGGCACAGCGCCACCGGGTATACAGATTGTTGACAAAGTGCGCTTTGTTTATGCCGGATGCGGCGTTAACGCCTTATCCGGCCTAC‑AAAATCGTGCAAATTCAATACATTG‑CA‑GAAT‑TTTC‑‑GTAGGCCTGATAAGCGGAGCGCATCAGGCATTTTTGCC  >  NC_000913/3705871‑3706069
                                                                                                   |                                                                                                         
cAGGGAATGAAGTATTGCTGCCGGATGGCGCTGCACGCAGAGGCACAGCGCCACCGGGTATACAGATTGTTGACAAAGTGCGCTTTGTTTATGCCGGATg                                                                                                           >  1:1818152/1‑100 (MQ=255)
 aGGGAATGAAGTATTGCTGCCGGATGGCGCTGCACGCAGAGGCACAGCGCCACCGGGTATACAGATTGTTGACAAAGTGCGCTTTGTTTATGCCGGATgc                                                                                                          <  1:40734/100‑1 (MQ=255)
  gggAATGAAGTATTGCTGCCGGATGGCGCTGCACGCAGAGGCACAGCGCCACCGGGTCTACAGCTTGTTGACAAAGTGCGCTTTGTTTATGCCGGATgcg                                                                                                         >  1:3896473/1‑100 (MQ=255)
  gggAATGAAGTATTGCTGCCGGATGGCGCTGCACGCAGAGGCACAGCGCCACCGGGTATACAGATTGTTGACAAAGTGCGCTTTGTTTATGCCGGATgcg                                                                                                         >  1:2992930/1‑100 (MQ=255)
  gggAATGAAGTATTGCTGCCGGATGGCGCTGCACGCAGAGGCACAGCGCCACCGGGTATACAGATTGTTGACAAAGTGCGCTTTGTTTATGCCGGATgcg                                                                                                         <  2:939433/100‑1 (MQ=255)
   ggAATGAAGTATTGCTGCCGGATGGCGCTGCACGCAGAGGCACAGCGCCACCGGGTATACAGATTGTTGACAAAGTGCGCTTTGTTTATGCCGGATgcgg                                                                                                        >  1:4027725/1‑100 (MQ=255)
    gAATGAAGTATTGCTGCCGGATGGCGCTGCACGCAGAGGCACAGCGCCACCGGGTATACAGATTGTTGACAAAGTGCGCTTTGTTTATGCCGGATgcggc                                                                                                       <  2:1722038/100‑1 (MQ=255)
      aTGAAGTATTGCTGCCGGATGGCGCTGCACGCAGAGGCACAGCGCCACCGGGTATACAGATTGTTGACAAAGTGCGCTTTGTTTATGCCGGATgcg                                                                                                         <  1:466135/96‑1 (MQ=255)
        gAAGTATTGCTGCCGGATGGCGCTGCACGCAGAGGCACAGCGCCACCGGGTATACAGATTGTTGACAAAGTGCGCTTTGTTTATGCCGGATGCGGCGTTa                                                                                                   >  2:3772715/1‑100 (MQ=255)
        gAAGTATTACTGCCGGATGGCGCTGCACGCAGAGGCACAGCGCCACCGGGTATACAGATTGTTGACAAAGTGCGCTTTGTTTATGCCGGATGCGGCGTTa                                                                                                   >  2:614698/1‑100 (MQ=255)
         aaGTATTGCTGCCGGATGGCGCTGCACGCAGAGGCACAGCGCCACCGGGTATACAGATTGTTGACAAAGTGCGCTTTGTTTATGCCGGATGCGGCGTTaa                                                                                                  >  1:3458102/1‑100 (MQ=255)
            tATTGCTGCCGGATGGCGCTGCACGCAGAGGCACAGCGCCACCGGGTATACAGATTGTTGACAAAGTGCGCTTTGTTTATGCCGGATGCGGCGTTAACGc                                                                                               <  1:3068014/100‑1 (MQ=255)
             aTTGCTGCCGGATGGCGCTGCATGCAGAGGCACAGCGCCACCGGGTATACAGATTGTTGACAAAGTGCGCTTTGTTTATGCCGGATGCGGCGTTAACGcc                                                                                              <  1:2138372/100‑1 (MQ=255)
               tgctgcCGGATGGCGCTGCACGCAGAGGCACAGCGCCACCGGGTATACAGATTGTTGACAAAGTGCGCTTTGTTTATGCCGGATGCGGCGTTAACGCCtt                                                                                            >  2:1927344/1‑100 (MQ=255)
               tgctgcCGGATGGCGCTGCACGCAGAGGCACAGCGCCACCGGGTATACAGATTGTTGACAAAGTGCGCTTTGTTTATGCCGGATGCGGCGTTAACGCCt                                                                                             <  1:173194/99‑1 (MQ=255)
                gctgcCGGATGGCGCTGCACGCAGAGGCACAGCGCCACCGGGTATACAGATTGTTGACAAAGTGCGCTTTGTTTATGCCGGATGAGGCGTTAACGCCTTa                                                                                           >  2:4205702/1‑100 (MQ=255)
                  tgcCGGATGGCGCTGCACGCAGAGGCACAGCGCCACCGGGTATACAGATTGTTGACAAAGTGCGCTTTGTTTATGCCGGATGCGGCTTTAACGCCTTacc                                                                                         <  2:2898226/100‑3 (MQ=255)
                   gcCGGATGGCGCTGCACGCAGAGGCACAGCGCCACCGGGTATACAGATTGTTGACAAAGTGCGCTTTGTTTATGCCGGATGCGGCGTTAACGCCTTATcc                                                                                        <  1:1128012/100‑1 (MQ=255)
                   gcCGGATGGCGCTGCACGCAGAGGCACAGCGCCACCGGGTATACAGATTGTTGACAAAGTGCGCTTTGTTTATGCCGGATGCGGCGTTAACGCCTTATcc                                                                                        >  2:725844/1‑100 (MQ=255)
                   gcCGGATGGCGCTGCACGCAGAGGCACAGCGCCACCGGGTATACAGATTGTTGACAAAGTGCGCTTTGTTTATGCCGGATGCGGCGTTAACGCCTTATcc                                                                                        >  1:815914/1‑100 (MQ=255)
                   gcCGGATGGCGCTGCACGCAGAGGCACAGCGCCACCGGGTATACAGATTGTTGACAAAGTGCGCTTTGTTTATGCCGGATGCGGCGTTAACGCCTTATcc                                                                                        >  2:250685/1‑100 (MQ=255)
                   gcCGGATGGCGCTGCACGCAGAGGCACAGCGCCACCGGGTATACAGATTGTTGACAAAGTGCGCTTTGTTTATGCCGGATGCGGCGTTAACGCCTTATcc                                                                                        >  2:2697983/1‑100 (MQ=255)
                      ggATGGCGCTGCACGCAGAGGCACAGCGCCACCGGGTATACAGATTGTTGACAAAGTGCGCTTTGTTTATGCCGGATGCGGCGTTAACGCCTTATCCgg                                                                                      <  1:1318253/99‑1 (MQ=255)
                       gATGGCGCTGCACGCAGAGGCACAGCGCCACCGGGTATACAGATTGTTGACAAAGTGCGCTTTGTTTATGCCGGATAAGGCGTTAACGCCGCATCCGGca                                                                                    <  2:3968923/100‑2 (MQ=37)
                       gATGGCGCTGCACGCAGAGGCACAGCGCCACCGGGTATACAGATTGTTGACAAAGTGCGCTTTGTTTATGCCGGATAAGGCGTTAACGCCGCATCCGGca                                                                                    <  2:2985501/100‑2 (MQ=37)
                       gATGGCGCTGCACGCAGAGGCACAGCGCCACCGGGTATACAGATTGTTGACAAAGTGCGCTTTGTTTATGCCGGATAAGGCGTTAACGCCGCATCCGGca                                                                                    <  1:3732519/100‑2 (MQ=37)
                       gATGGCGCTGCACGCAGAGGCACAGCGCCACCGGGTATACAGATTGTTGACAAAGTGCGCTTTGTTCATTCCGGATAAGGCGTTAACGCCGCATCCGGca                                                                                    <  1:2808494/100‑2 (MQ=37)
                       gATGGCGCTGCACGCAGAGGCACAGCGCCACCGGGTATACAGATTGGTGACAAAGTGCGCTTTGTTTATGCCGGATGCGGCGTTAACGCCTTATCCGGcc                                                                                    <  1:3149483/100‑1 (MQ=255)
                        aTGGCGCTGCACGCAGAGGCACAGCGCCACCGGGTATACAGATTGTTGACAAAGTGCGCTTTGTTTATGCCGGATGCGGCGTTAACGCCTTATCCGGCCt                                                                                   <  2:557144/100‑1 (MQ=255)
                          ggCGCTGCACGCAGAGGCACAGCGCCACCGGGTATACAGATTGTTGACAAAGTGCGCTTTGTTTATGCCGGATGCGGCGTTAACGCCTTATCCGGCCTAc                                                                                 <  2:2547563/100‑1 (MQ=255)
                           gcgcTGCACGCAGAGGCACAGCGCCACCGGGTATACAGATTGTTGACAAAGTGCGCTTTGTTTATGCCGGATGCGGCGTTAACGCCTTATCCGGCCTAC‑a                                                                               <  2:3414368/100‑1 (MQ=255)
                           gcgcTGCACGCAGAGGCACAGCGCCACCGGGTATACAGATTGTTGACAAAGTGCGCTTTGTTTATGCCGGATGCGGCGTTAACGCCTTATCCGGCCTAC‑a                                                                               >  1:1868165/1‑100 (MQ=255)
                            cgttGCACGCAGAGGTACAGCGCCCCCGGGTATATAGATTGTTGTCAAAGTGCGCTTTGTTTATGCCGGATGCGGCGTTAACGCCTTATCCGGCCTAC‑a                                                                               <  1:3858042/96‑1 (MQ=21)
                            cgcTGTTCGCAGAGGCACAGCGCTACCGGGTATACAGATTGGTGACAAAGTGCGCTTTGTTTATTCCGGATAAGGCGTTATCGCCGCATCCGGCATaaac                                                                               <  1:1002361/100‑4 (MQ=21)
                               tGCAGGCAGAGGCACAGCGCCACCGGGTATACAGATTGTTGACAAAGTGCGCTTTGTTTATGCCGGATGCGGCGTTAACGCCTTATCCGGCCTAC‑AAAAt                                                                           <  1:2319343/100‑1 (MQ=255)
                               tGCACGCAGAGGCACAGCGCCACCGGGTATACAGATTGTTGACAAAGTGCGCTTTGTTTATGCCGGATGCGGCGTTAACGCCTTATCCGGCCTAC‑AAAAt                                                                           <  2:1124123/100‑1 (MQ=255)
                                  aCGCAGAGGCACAGCGCCACCGGGTATACAGATTGTTGACAAAGTGCGCTTTGTTTATGCCGGATGCGGCGTTAACGCCTTATCCGGCCTAC‑AAAATCGt                                                                        >  1:4011915/1‑100 (MQ=255)
                                   cGCAGAGGCACAGCGCCACCGGGTATACAGATTGTTGACAAAGTGCGCTTTGTTTATGCCGGATGCGGCGTTAACGCCTTATCCGGCCTAC‑AAAATCGTg                                                                       <  2:1050710/100‑1 (MQ=255)
                                   cGCAGAGGCACAGCGCCACCGGGTATACAGATTGTTGACAAAGTGCGCTTTGTTTATGCCGGATGCGGCGTTAACGCCTTATCCGGCCTAC‑AAAATCGTg                                                                       <  1:147500/100‑1 (MQ=255)
                                   cGCAGAGGCACAGCGCCACCGGGTATACAGATTGTTGACAAAGTGCGCTTTGTTTATGCCGGATGCGGCGTTAACGCCTTATCCGGCCTAC‑AAAATCGTg                                                                       >  2:3558195/1‑100 (MQ=255)
                                    gCAGAGGCACAGCGCCACCGGGTATACAGATTGTTGACAAAGTGCGCTTTGTTTATGCCGGATGCGGCGTTAACGCCTTATCCGGCCTAC‑AAAATCGTGc                                                                      <  2:1994244/100‑1 (MQ=255)
                                       gagGCACAGCGCCACCGGGTCTACAGATTGTTGACAAAGTGCGCTTTGTTTATGCCGGATGCGGCGTTAACGCCTTATCCGGCCTAC‑AAAATCGTGCaaa                                                                   >  1:2605970/1‑100 (MQ=255)
                                       gagGCACAGCGCCACCGGGTATACAGATTGTTGACAAAGTGCGCTTTGTTTATGCCGGATGCGGCGTTAACGCCTTATCCGGCCTAC‑AAAATCGTGCaaa                                                                   <  2:2716190/100‑1 (MQ=255)
                                        agGCACAGCGCCACCGGGTATACAGATTGTTGACAAAGTGCGCTTTGTTTATGCCGGATGCGGCGTTAACGCCTTATCCGGCCTAC‑AAAATCGTGCAAAt                                                                  <  2:2293129/100‑1 (MQ=255)
                                             caGCGCCACCGGGTATACAGATTGTTGACAAAGTGCGCTTTGTTTATGCCGGATGCGGCGTTAACGCCTTATCCGGCCTAC‑AAAATCGTGCAAATTCAAt                                                             >  2:318890/1‑100 (MQ=255)
                                                 gcCACCGGGTATACAGATTGTTGACAAAGTGCGCTTTGTTTATGCCGGATGCGGCGTTAACGCCTTATCCGGCCTAC‑AAAATCGTGCAAATTCAATACAt                                                         >  1:224486/1‑100 (MQ=255)
                                                  ccACCGGGTATACAGATTGTTGACAAAGTGCGCTTTGTTTATGCCGGATGCGGCGTTAACGCCTTATCCGGCCTAC‑AAAATCGTGCAAATTCAATACAtt                                                        <  2:2602896/100‑1 (MQ=255)
                                                    aCCGGGTATACAGATTGTTGACAAAGTGCGCTTTGTTTATGCCGGATGCGGCGTTAACGCCTTATCCGGCCTAC‑AAAATCGTGCAAATTCAATACATTG‑c                                                     <  2:2502954/100‑1 (MQ=255)
                                                    aCCGGGTATACAGATTGTTGACAAAGTGCGCTTTGTTTATGCCGGATGCGGCGTTAACGCCTTATCCGGCCTAC‑AAAATCGTGCAAATTCAATACATTG‑c                                                     >  2:3288547/1‑100 (MQ=255)
                                                      cGGGTATACAGATTGTTGACAAAGTGCGCTTTGTTTATGCCGGATGCGGCGTTAACGCCTTATCCGGCCTAC‑AAAATCGTGCAAATTCAATACATTG‑CA‑g                                                  <  1:3330410/100‑1 (MQ=255)
                                                      cGGGTATACAGATTGTTGACAAAGTGCGCTTTGGTTATGCCGGATGCGGCGTTAACGCCTTATCCGGCCTCC‑AAAATCGTGCAAATTCAATACATTG‑CA‑g                                                  >  1:2564638/1‑100 (MQ=255)
                                                          tataCAGATTGTTGACAAAGTGCGCTTTGTTTATGCCGGATGCGGCGTTAACGCCTTATCCGGCCTAC‑AAAATCGTGCAAATTCAATACATTG‑CA‑GAAt‑t                                             >  1:3037125/1‑100 (MQ=255)
                                                             aCAGATTGTTGACAAAGTGCGCTTTGTTTATGCCGGATGCGGCGTTAACGCCTTATCCGGCCTAC‑AAAATCGTGCAAATTCAATACATTG‑CA‑GAAT‑TTTc                                          >  1:577399/1‑100 (MQ=255)
                                                              cAGATTGTTGACAAAGTGCGCTTTGTTTATGCCGGATGCGGCGTTAACGCCTTATCCGGCCTAC‑AAAATCGTGCAAATTCAATACATTG‑CA‑GAAT‑TTTC‑‑g                                       <  2:2109352/100‑1 (MQ=255)
                                                                gATTGTTGACAAAGTGCGCTTTGTTTATGCCGGATGCGGCGTTAACGCCTTATCCGGCCTAC‑AAAATCGTGCAAATTCAATACATTG‑CA‑GAAT‑TTTC‑‑GTa                                     <  2:1379094/100‑1 (MQ=255)
                                                                  gcataaaCAAAGCGCACTTTGTTTATGCCGGATGCGGCGTTAACGCCTTATCCGGCCTAC‑AAAATCGTGCAAATTCAATACATTG‑CA‑GAAT‑TTTC‑‑GTAg                                    >  2:2848951/7‑99 (MQ=17)
                                                                   tgttgACAAAGTGCGCTTTGTTTATGCCGGATGCGGCGTTAACGCCTTATCCGGCCTAC‑AAAATCGTGCAAATTCAATACATTG‑CA‑GAAT‑TTTC‑‑GTAGGc                                  >  2:3213023/1‑100 (MQ=255)
                                                                    gttgACAAAGTGCGCTTTGTTTATGCCGGATGCGGCGTTAACGCCTTATCCGGCCTAC‑AAAATCGTGCAAATTCAATACATTG‑CA‑GAAT‑TTTC‑‑GTAGGcc                                 <  1:644276/100‑1 (MQ=255)
                                                                    gttgACAAAGTGCGCTTTGTTTATGCCGGATGCGGCGTTAACGCCTTATCCGGCCTAC‑AAAATCGTGCAAATTCAATACATTG‑CA‑GAAT‑TTTC‑‑GTAGGcc                                 <  1:2871173/100‑1 (MQ=255)
                                                                    gttgACAAAGTGCGCTTTGTTTATGCCGGATGCGGCGTTAACGCCTTATCCGGCCTAC‑AAAATCGTGCAAATTCAATACATTG‑CA‑GAAT‑TTTC‑‑GTAGGcc                                 >  1:727115/1‑100 (MQ=255)
                                                                      tgACAAAGTGCGCTTTGTTTATGCCGGATGCGGCGTTAACGCCTTATCCGGCCTAC‑AAAATCGTGCAAATTCAATACATTG‑CA‑GAAT‑TTTC‑‑GTAggagat                               >  2:1686246/1‑96 (MQ=255)
                                                                       gACAAAGTGCGCTTTGTTTATGCCGGATGCGGCGTTAACGCCTTATCCGGCCTAC‑AAAATCGTGCAAATTCAATACATTG‑CA‑GAAT‑TTTC‑‑GTAGGCCTGa                              <  1:3536186/100‑1 (MQ=255)
                                                                         cAAAGTGCGCTTTGTTTATGCCGGATGCGGCGTTAACGCCTTATCCGGCCTAC‑AAAATCGTGCAAATTCAATACATTG‑CA‑GAAT‑TTTC‑‑GTAGGCCTGATa                            <  2:67276/100‑1 (MQ=255)
                                                                         cAAAGTGCGCTTTGTTTATGCCGGATGCGGCGTTAACGCCTTATCCGGCCTAC‑AAAATCGTGCAAATTCAATACATTG‑CA‑GAAT‑TTTC‑‑GTAGGCCTGATa                            <  1:1218830/100‑1 (MQ=255)
                                                                             gTGCGCTTTGTTTATGCCGGATGTGGCGTTAACTCCTTATCCGGCCTAC‑AAAATCGTGCAAATTCAATACATTG‑CA‑GAAT‑TTTC‑‑GTAGGCCTGATAAGCg                        <  2:1842329/100‑1 (MQ=255)
                                                                             gTGCGCTTTGTTTATGCCGGATGCGGCGTTAACGCCTTATCCGGCCTAC‑AAAATCGTGCAAATTCAATACATTG‑CA‑GAAT‑TTTC‑‑GTAGGCCTGATAAGCg                        <  1:1978695/100‑1 (MQ=255)
                                                                                 gcTTTGTTTATGCCGGATGCGGCGTTAACGCCTTATCCGGCCTAC‑AAAATCGTGCAAATTCAATACATTG‑CA‑GAAT‑TTTC‑‑GTAGGCCTGATAAGCGGAgc                    <  1:1567754/100‑1 (MQ=255)
                                                                                 gcTTTGTTTATGCCGGATGCGGCGTTAACGCCTTATCCGGCCTAC‑AAAATCGTGCAAATTCAATACATTG‑CA‑GAAT‑TTTC‑‑GTAGGCCTGATAAGCGGAgc                    >  1:2747342/1‑100 (MQ=255)
                                                                                  gattttgtaggCCGGATAAGGCGTTCACGCCGCATCCGGCCTAC‑AAAATCGCTTAAATTCAATATATTG‑CA‑CCAT‑CCTC‑‑GTAGGCCTGATAAGCGTAgcg                   >  2:1553592/11‑100 (MQ=18)
                                                                                  cTTTGTTTATGCCGGATGCGGCGTTAACGCCTTATCCGGCCTAC‑AAAATCGTGCAAATTCAATACATTG‑CA‑GAAT‑TTTC‑‑GTAGGCCTGATAAGCGGAgcg                   <  1:2047092/100‑1 (MQ=255)
                                                                                  cTTTGTTTATGCCGGATGCGGCGTTAACGCCTTATCCGGCCTAC‑AAAATCGTGCAAATTCAATACATTG‑CA‑GAAT‑TTTC‑‑GTAGGCCTGATAAGCGGAgcg                   <  2:1782461/100‑1 (MQ=255)
                                                                                  cTTTGTTTATGCCGGATGCGGCGTTAACGCCTTATCCGGCCTAC‑AAAATCGTGCAAATTCAATACATTG‑CA‑GAAT‑TTTC‑‑GTAGGCCTGATAAGCGGAgcg                   <  1:1531116/100‑1 (MQ=255)
                                                                                    ttGTTTATGCCGGATGCGGCGTTAACGCCTTATCCGGCCTAC‑AAAATCGTGCAAATTCAATACATTG‑CA‑GAAT‑TTTC‑‑GTAGGCCTGATAAGCGGAGCGCa                 <  1:3805391/100‑1 (MQ=255)
                                                                                    ttGTTTATGCCGGATGCGGCGTTAACGCCTTATCCGGCCTAC‑AAAATCGTGCAAATTCAATACATTG‑CA‑GAAT‑TTTC‑‑GTAGGCCTGATAAGCGGAGCGCa                 >  2:3497768/1‑100 (MQ=255)
                                                                                     tttgtaggCCGGATAAGGCGTTTACGCCGCATCCGGCCTAC‑AAAATCGTGCAAATTCAACATATTG‑CA‑ATTC‑TCTT‑‑GTAGGCCTGATAAGCGTAGCGCAt                >  2:27686/8‑100 (MQ=18)
                                                                                     tGTTTATGCCGGATGCGGCGTTAACGCCTTATCCGGCCTAC‑AAAATCGTGCAAATTCAATACATTG‑CA‑GAAT‑TTTC‑‑GTAGGCCTGATAAGCGGAGCGCAt                >  1:1600225/1‑100 (MQ=255)
                                                                                       tgtaggCCGGATAAGGCGTTTACGCCGCATCCGGCCTAC‑AAAATCGTGCTAATTCAATATATTG‑CA‑GAAA‑CCTT‑‑GTAGGCCTGATAAGCGTAGCGCATCa              >  2:2419993/6‑100 (MQ=18)
                                                                                       tgtaggCCGGATAAGGCGCTTACGCCGCATCCGGCCTAC‑AAGAACGTGCAAATTCAACAAATTGTGATGCAT‑CCT‑‑‑GTAGGCCTGATAAGCGTAGCGCATc               >  2:115490/6‑100 (MQ=21)
                                                                                         taggCCGGATAAGGCGTTTACGCCGCATCTGGCCTAC‑AAATTCTTCCAAATTCAATATATTG‑CA‑GAAA‑TCAT‑‑GTAGGCCTGATAAGCGTAGCGCATCAgg            >  2:289032/4‑100 (MQ=18)
                                                                                         taggCCGGATAAGGCGTTTACGCCGCATCCGGCCTAC‑AAAAGCATGCAAATTCAATAAATTG‑CA‑GAGT‑TCAC‑‑GTAGGCCTGATAAGCGTAGCGCATCAgg            >  1:3001215/4‑100 (MQ=17)
                                                                                         taggCCGGATAAGGCGTTCACGCCGCATCCGGCCTAC‑AAAAACATGCAAATTCAATAGATTG‑CA‑GAGA‑TTAT‑‑GTAGGCCTGATAAGCGTAGCGCATCAgg            >  2:3130402/4‑100 (MQ=17)
                                                                                          aggCCGGATAAGGCGTTCACGCCGCATCCGGCCTACAAAAATCTTGCCAATTCAATATATTG‑CA‑GGAC‑CAT‑‑‑GTAGGCCTGATAAGCGTAGCGCATCAGGc           >  2:384718/3‑100 (MQ=18)
                                                                                          aggCCGGATAAGGCGTTCACGCCGCATCCGGCCTAC‑AAAACTGTGCAAATTCAATATATTG‑CA‑GGAA‑ACAC‑‑GTAGGCCTGATAAGCGAAGCGCATCAGGc           >  2:504554/3‑100 (MQ=18)
                                                                                          aggCCGGATAAGGCGTTCACGCCGCATCCGGCCTAC‑AAAACTGTGCAAATTCAATATATTG‑CA‑GGAA‑ACAC‑‑GTAGGCCTGATAAGCGAAGCGCATCAGGc           >  1:955708/3‑100 (MQ=18)
                                                                                          aggCCGGATAAGGCGTTCACGCCGCATCCGGCCTAC‑AAAACCTTGCTAATTCAATATATTG‑CA‑GGGA‑CTAT‑‑GTAGGCCTGATAAGCATAGCGCATCAGGc           >  1:2434933/3‑100 (MQ=18)
                                                                                          aTGCCGGATGCGGCGTTAACGCCTTATCCGGCCTAC‑AAAATCGTGCAAATTCAATACATTG‑CA‑GAAT‑TTTC‑‑GTAGGCCTGATAAGCGGAGCGCATCAGGc           <  1:1432903/100‑1 (MQ=255)
                                                                                           tGCCGGATGCGGCGTTAACGCCTTATCCGGCCTAC‑AAAATCGTGCAAATTCAATACATTG‑CA‑GAAT‑TTTC‑‑GTAGGCCTGATAAGCGGAGCGCATCAGGCa          >  2:1910432/1‑100 (MQ=255)
                                                                                           ggCCGGATAAGGCGTTTACGCCGCATCTGGCCTAC‑AAATTCTTCCAAATTCAATATATTG‑CA‑GAAA‑TCAT‑‑GTAGGCCTGATAAGCGTAGCGCATCAGGCa          >  2:4006677/2‑100 (MQ=18)
                                                                                           ggCCGGATAAGGCGTTCACGCCGCATCCGGCCTAC‑GAAATCGTGCTAATTCAATATATTG‑TG‑GGCT‑CCAC‑‑GTAGGCCTGATAAGCGTAGCGCATCAGGCa          >  1:1959008/2‑100 (MQ=18)
                                                                                           ggCCGGATAAGGCGTTCACGCCGCATCCGGCCTAC‑AAAATCGCTTAAATTCAATATATTG‑CA‑CCAT‑CCTC‑‑GTAGGCCTGATAAGCGTAGCGCATCAGGCa          >  2:2145807/2‑100 (MQ=21)
                                                                                           ggCCGGATAAGGCGTTCACGCCGCATCCGGCCTAC‑AAAACTTTGCAAATTCAATATATTG‑C‑‑‑AAT‑CTCCGTGTAGGCCTGATAAGCGTAGCGCATCAGGCa          >  2:2640927/2‑100 (MQ=18)
                                                                                           ggCCGGATAAGGCGTTCACGCCGCATCCGGCCTAC‑AAAACTGTGCAAATTCAATATATTG‑CA‑GGAA‑ACAC‑‑GTAGGCCTGATAAGCGAAGCGCATCAGGCa          >  2:2077125/2‑100 (MQ=18)
                                                                                            gCCGGATGCGGCGTTAACGCCTTATCCGGCCTAC‑AAAATCGTGCAAATTCAATACATTG‑CA‑GAAT‑TTTC‑‑GTAGGCCTGATAAGCGGAGCGCATCAGGCa          >  1:3434330/1‑99 (MQ=255)
                                                                                            gCCGGATGCGGCGTTAACGCCTTATCCGGCCTAC‑AAAATCGTGCAAATTCAATACATTG‑CA‑GAAT‑TTTC‑‑GTAGGCCTGATAAGCGGAGCGCATCAGGCAt         >  2:4083402/1‑100 (MQ=255)
                                                                                            gCCGGATAAGGCGTTTACGCCGCATCCGGCCTAC‑AAATTC‑TGCAAATTCAACAAATTG‑C‑‑GAATCCCTT‑‑GTAGGCCTGATAAGCGTAGCGCATCAGGCAtt        >  1:1646114/1‑100 (MQ=18)
                                                                                            gCCGGATAAGGCGTTCACGCCGCATCCGGCCTAC‑AAGATCGTGCAAATTCAATATATTG‑CA‑TGAG‑CCAC‑‑GTAGGCCTGATAAGCGTAGCGCATCAGGCag         >  1:2953754/1‑99 (MQ=17)
                                                                                             ccGGATGCGGCGTTAACGCCTTATCCGGCCTAC‑AAAATCGTGCAAATTCAATACATTG‑CA‑GAAT‑TTTC‑‑GTAGGCCTGATAAGCGGAGCGCATCAGGCAtt        >  2:3725697/1‑100 (MQ=255)
                                                                                                 aTGCGGCGTTAACGTCTTATCCGGCCTAC‑AAAATCGGGCAAATTCAATACATTG‑CA‑GAAT‑TTTC‑‑GTAGGCCTGATAAGCGGAGCGCATCAGGCATTtttg    <  1:2836547/100‑1 (MQ=255)
                                                                                                  tGCGGCGTTAACGCCTTATCCGGCCTAC‑AAAATCGTGCAAATTCAATACATTG‑CA‑GAAT‑TTTC‑‑GTAGGCCTGATAGGCGGAGCGCATCAGGCATTtttgc   >  2:1977453/1‑100 (MQ=255)
                                                                                                   gcggcgTTAACGCCTTATCCGGCCTAC‑AAAATCGTGCAAATTCAATACATTG‑CA‑GAAT‑TTTC‑‑GTAGGCCTGATAAGCGGAGCGCATCAGGCATTtttgcc  <  2:224487/100‑1 (MQ=255)
                                                                                                   |                                                                                                         
CAGGGAATGAAGTATTGCTGCCGGATGGCGCTGCACGCAGAGGCACAGCGCCACCGGGTATACAGATTGTTGACAAAGTGCGCTTTGTTTATGCCGGATGCGGCGTTAACGCCTTATCCGGCCTAC‑AAAATCGTGCAAATTCAATACATTG‑CA‑GAAT‑TTTC‑‑GTAGGCCTGATAAGCGGAGCGCATCAGGCATTTTTGCC  >  NC_000913/3705871‑3706069

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: ‑