Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 3,705,970 G→A 21.3% 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,9700GA21.3% 523.7 / 126.1 231intergenic (‑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 (85/95);  new base A (44/5);  total (129/102)
Fisher's exact test for biased strand distribution p-value = 2.92e-08
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.51e-02

Only 100 of 240 total aligned reads displayed.

CAGGGAATGAAGTATTGCTGCCGGATGGCGCTGCACGCAGAGGCACAGCGCCACCGGGTATACAGATTGTTGACAAAGTGCGCTTTGTTTATGCCGGATGCGGCGTTAACGCCTTATCCGGCCTACAAAATCGTGCAAATTCAATACATTGCAGAAT‑‑TTTC‑‑GTAGGCCTGATAAGCGGAGCGCATCAGGCATTTTTGCC  >  NC_000913/3705871‑3706069
                                                                                                   |                                                                                                         
cAGGGAATGAAGTATTGCTGCCGGATGGCGCTGCACGCAGAGGCACAGCGCCACCGGGTATACAGATTGTTGACAAAGTGCGCGTTGTTTATGCCGGATg                                                                                                           >  1:2165297/1‑100 (MQ=255)
  gggAATGAAGTATTGCTGGCGGATGGCGCTGCACGCAGAGGCACAGCGCCACCGGGTATACAGATTGTTGACAAAGTGCGCTTTGTTTATGCCGGATgc                                                                                                          <  1:3421821/99‑1 (MQ=255)
    gAATGAAGTATTGCTGCCGGATGGCGCTGCACGCAGAGGCACAGCGCCACCGGGTATACAGATTGTTGACAAAGTGCGCTTTGTTTATGCCGGATgcggc                                                                                                       <  1:1965877/100‑1 (MQ=255)
      aTGAAGTATTGCTGCCGGATGGCGCTGCACGCAGAGGCACAGCGCCACCGGGTATACAGATTGTTGACAAAGTGCGCTTTGTTTATGCCGGATGCGGCGt                                                                                                     <  2:2345774/100‑1 (MQ=255)
      aTGAAGTATTGCTGCCGGATGGCGCTGCACGCAGAGGCACAGCGCCACCGGGTATACAGATTGTTGACAAAGTGCGCTTTGTTTATGCCGGATGCGGCGt                                                                                                     >  2:2939286/1‑100 (MQ=255)
        gAAGTATTGCTGCCGGATGGCGCTGCACGCAGAGGCACAGCGCCACCGGGTATACAGATTGTTGACAAAGTGCGCTTTGTTTATGCCGGATGCGGCGTTa                                                                                                   >  2:1375669/1‑100 (MQ=255)
        gAAGTATTGCTGCCGGATGGCGCTGCACGCAGAGGCACAGCGCCACCGGGTATACAGATTGTTGACAAAGTGCGCTTTGTTTATGCCGGATGCGGCGTTa                                                                                                   <  1:2097709/100‑1 (MQ=255)
         aaGTATTGCTGCCGGATGGCGCTGCACGCAGAGGCACAGCGCCACCGGGTATACAGATTGTTGACAAAGTGCGCTTTGTTTATGCCGGATGCGGCGTTaa                                                                                                  >  1:1002174/1‑100 (MQ=255)
           gTATTGCTGCCGGATGGCGCTGCACGCAGAGGCACAGCGCCACCGGGTATACAGATTGTTGACAAAGTGCGCTTTGTTTATGCCGTATGCGGCGTTAACg                                                                                                <  1:927692/100‑1 (MQ=255)
           gTATTGCTGCCGGATGGCGCTGCACGCAGAGGCACAGCGCCACCGGGTATACAGATTGTTGACAAAGTGCGCTTTGTTTATGCCGGATGCGGCGTTAACg                                                                                                <  2:1319892/100‑1 (MQ=255)
              ttGCTGCCGGATGGCGCTGCACGCAGAGGCACAGCGCCACCGGGTATACAGATTGTTGACAAAGTGCGCTTTGTTTATGCCGGATGCGGCGTTAACGCCt                                                                                             <  1:2744593/100‑1 (MQ=255)
               tgctgcCGGATGGCGCTGCACGCAGAGGCACAGCGCCACCGGGTATACAGATTGTTGACAAAGTGCGCTTTGTTTATGCCGGATGCGGCGTTAACGCCtt                                                                                            >  2:1664284/1‑100 (MQ=255)
                gctgcCGGATGGCGCTGCACGCAGAGGCACAGCGCCACCGGGTATACAGATTGTTGACAAAGTGCGCTTTGTTTATGCCGGATGCGGCGTTAACGCCTTa                                                                                           >  1:441403/1‑100 (MQ=255)
                gctgcCGGATGGCGCTGCACGCAGAGGCACAGCGCCACCGGGTATACAGATTGTTGACAAAGTGCGCTTTGTTTATGCCGGATGCGGCGTTAACGCCTTa                                                                                           >  2:2561829/1‑100 (MQ=255)
                gctgcCGGATGGCGCTGCACGCAGAGGCACAGCGCCACCGGGTATACAGATTGTTGACAAAGTGCGCTTTGTTTATGCCGGATGCGACGTTAACGCCTTa                                                                                           >  1:605633/1‑100 (MQ=255)
                 ctgcCGGATGGCCCTGCACGCAGAGGCACAGCGCCACCGGGTATACAGATTGTTGACAAAGTGCGCTTTGTTTATGCCGGATGCGGCGTTAACGCCTTAt                                                                                          >  2:2433862/1‑100 (MQ=255)
                   gcCGGATGGCGCTGCACGCAGAGGCACAGCGCCACCGGGTATACAGATTGTTGACAAAGTGCGCTTTGTTTATGCCGGATGCGGCGTTAACGCCTTATcc                                                                                        <  2:2834720/100‑1 (MQ=255)
                      ggATGGCGCTGCACGCAGAGGCACAGCGCCACCGGGTATACAGATTGTTGGCAAAGTGCGCTTTGTTTATGCCGGATGCGGCGTTAACGCCTTATCCGGc                                                                                     <  2:677890/100‑1 (MQ=255)
                      ggATGGCGCTGCACGCAGAGGCACAGCGCCACCGGGTATACAGATTGTTGACAAAGTGCGCTTTGTTTATGCCGGATGCGGCGTTAACGCCTTATCCGGc                                                                                     <  1:2646716/100‑1 (MQ=255)
                       gATGGCGCTGCACGCAGAGGCACAGCGCCACCGGGTATACAGATTGTTGACAAAGTGCGCTTTGTTTATGCCGGATGCGGCGTTAACGCCTTATCCGGcc                                                                                    <  2:3364516/100‑1 (MQ=255)
                          ggCGCTGCACGCAGAGGCACAGCGCCACCGGGTATACAGATTGTTGACAAAGTGCGCTTTGTTTATGCCGGATAAGGCGTTAACGCCGCATCCGGCATaa                                                                                 <  2:1912245/100‑2 (MQ=25)
                           gcgcTGCACGCAGAGGCACAGCGCCACCGGGTATACAGATTGTTGACAAAGTGCGCTTTGTTTATGCCGGATGCGGCGTTAACGCCTTATCCGGCCTACa                                                                                >  2:3253249/1‑100 (MQ=255)
                             gcTGCACGCAGAGGCACAGCGCCACCGGGTATACAGATTGTTGACAAAGTGCGCTTTGTTTATGCCGGATGCGGCGTTAACGCCTTATCCGGCCTACaaa                                                                              >  2:1995120/1‑100 (MQ=255)
                             cctgtaggaagagGCGCAGTGCCACCGGGTATACAGATTGCTGACAAAGTGAGCTCTGATTATGCCCGATGCGTCGTTGACGCCTTATCCGGCCTAc                                                                                 >  1:1775629/10‑97 (MQ=17)
                              cTGCACGCAGAGGCACAGCGCCACCGGGTATACAGATTGTTGACAAAGTGCGCTTTGTTTATGCCGGATGCGGCGTTAACGCCTTATCCGGCCTACaaaa                                                                             >  1:2408557/1‑100 (MQ=255)
                                 cACGCAGAGGCACAGCGCCACCGGGTATACAGATTGTTGACAAAGTGCGCTTTGTTTATGCCGGATGCGGCGTTAACGCCTTATCCGGCCTACAAAATCg                                                                          >  1:1604397/1‑100 (MQ=255)
                                 cACGCAGAGGCACAGCGCCACCGGGTATACAGATTGTTGACAAAGTGCGCTTTGTTTATGCCGGATGCGGCGTTAACGCCTTATCCGGCCTACAAAATCg                                                                          <  1:3311611/100‑1 (MQ=255)
                                   cGCAGAGGCACAGCGCCACCGGGTATACAGATTGTTGACAAAGTGCGCTTTGTTTATGCCGGATGCGGCGTTAACGCCTTATCCGGCCTACAAAATCGTg                                                                        >  2:260996/1‑100 (MQ=255)
                                   cGCAGAGGCACAGCGCCACCGGGTATACAGATTGTTGACAAAGTGCGCTTTGTTTATGCCGGATGCGGCGTTAACGCCTTATCCGGCCTACAAAATCGTg                                                                        >  1:2985076/1‑100 (MQ=255)
                                   cGCAGAGGCACAGCGCCACCGGGTATACAGATTGTTGACAAAGTGCGCTTTGTTTATGCCGGATGCGGCGTTAACGCCTTATCCGGCCTACAAAATCGTg                                                                        <  1:3120527/100‑1 (MQ=255)
                                   cGCAGAGGCACAGCGCCACCGGGTATACAGATTGTTGACAAAGTGCGCTTTGTTTATGCCGGATGCGGCGTTAACGCCTTATCCGGCCTACAAAATCGTg                                                                        <  2:2313373/100‑1 (MQ=255)
                                   cGCAGAGGCACAGCGCCACCGGGTATACAGATTGTTGACAAAGTGCGCTTTGTTTATGCCGGATGCGGCGTTAACGCCTTATCCGGCCTACAAAATCATg                                                                        >  2:1064694/1‑100 (MQ=255)
                                    gCAGAGGCACAGCGCCACCGGGTATACAGATTGTTGACAAAGTGCGCTTTGTTTATGCCGGATGCGGCGTTAACGCCTTATCCGGCCTACAAAATCGTGc                                                                       <  1:299896/100‑1 (MQ=255)
                                     cAGAGGCACAGCGCCACCGGGTATACAGATTGTTGACAAAGTGCGCTTTGTTTATGCCGGATGCGGCGTTAACGCCTTATCCGGCCTACAAAATCGTGCa                                                                      >  2:618931/1‑100 (MQ=255)
                                     cAGAGGCACAGCGCCACCGGGTATACAGATTGTTGACAAAGTGCGCTTTGTTTATGCCGGATGCGGCGTTAACGCCTTATCCGGCCTACAAAATCGTGCa                                                                      >  2:2373106/1‑100 (MQ=255)
                                        agGCACAGCGCCACCGGGTATACAGATTGTTGACAAAGTGCGCTTTGTTTATGCCGGATGCGGCGTTAACGCCTTATCCGGCCTACAAAATCGTGCAAAt                                                                   <  2:1633978/100‑1 (MQ=255)
                                         ggCACAGCGCCACCGGGTATACAGATTGTTGACAAAGTGCGCTTTGTTTATGCCGGATGCGGCGTTAACGCCTTATCCGGCCTACAAAATCGTGCAAAtt                                                                  <  2:3226075/100‑1 (MQ=255)
                                          gCACAGCGCCACCGGGTATACAGATTGTTGACAAAGTGCGCTTTGTTTATGCCGGATGCGGCGTTAACGCCTTATCCGGCCTACAAAATCGTGCAAATTc                                                                 <  2:2658511/100‑1 (MQ=255)
                                          gCACAGCGCCACCGGGTATACAGATTGTTGACAAAGTGCGCTTTGTTTATGCCGGATGCGGCGTTAACGCCTTATCCGGCCTACAAAATCGTGCAAATTc                                                                 <  1:1387367/100‑1 (MQ=255)
                                           cacaGCGCCACCGGGTATACAGATTGTTGACAAAGTGCGCTTTGTTTATGCCGGATGCGGCGTTAACGCCTTATCCGGCCTACAAAATCGTGCAAATTCa                                                                >  2:1215504/1‑100 (MQ=255)
                                           cacaGCGCCACCGGGTATACAGATTGTTGACAAAGTGCGCTTCGTTTATGCCGGATGCGGCGTTAACGCCTTATCCGGCCTACAAAATCGTGCAAATTc                                                                 >  1:3240987/1‑99 (MQ=255)
                                             caGCGCCACCGTGTATACAGATTGTTGACAAAGTGCGCTTTGTTTATGCCGGATGCGGCGTTAACTCCTTATCCGGCCTACAAAATCCTGCAAATTCAAt                                                              <  1:808522/100‑1 (MQ=255)
                                             caGCGCCACCGGGTATACAGATTGTTGACAAAGTGCGCTTTGTTTATGCCGGATGCGGCGTTAACGCCTTATCCGGCCTACAAAATCGTGCAAATTCAAt                                                              <  2:2593043/100‑1 (MQ=255)
                                             caGCGCCACCGGGTATACAGATTGTTGACAAAGTGCGCTTTGTTTATGCCGGATGCGGCGTTAACGCCTTATCCGGCCTACAAAATCGTGCAAATTCAAt                                                              <  2:1107637/100‑1 (MQ=255)
                                                 gcCACCGAGTATGCAGATTGTTGACAAAGTGCTCTTTGTGTATGCCCGATGCGGCGTTAACGCCTTAGCAGGCCTACAAGCTCGTGCAAATTCAAt                                                              <  1:1725959/96‑1 (MQ=17)
                                                  ccACCGGGTATACAGATTGTTGACAAAGTGCGCTTTGTTTATGCCCGATGCGGCGTTAACGCCTTATTCGGCCTACAAAATCGTGCAAATTCAATAcct                                                          >  1:3166936/1‑97 (MQ=255)
                                                   cACCGGGTATACAGATTGTTGACAAAGTGCGCTTTGTTTATGCCGGATGCGGCGTTAACGCCTTATCCGGCCTACAAAATCGTGCAAATTCAATACATTg                                                        >  2:1792365/1‑100 (MQ=255)
                                                    aCCGGGTATACAGATTGTTGACAAAGTGCGCTTTGTTTATGCCGGATGCGGCGTTAACGCCTTATCCGGCCTACAAAATCGTGCAAATTCAATACATTgn                                                       <  1:2920512/100‑2 (MQ=255)
                                                    aCCGGGTATACAGATTGTTGACAAAGTGCGCTTTGTTTATGCCGGATGCGGCGTTAACGCCTTATCCGGCCTACAAAATCGTGCAAATTCAATACATTGc                                                       <  1:2619774/100‑1 (MQ=255)
                                                      cGGGTATACAGATTGTTGACAAAGTGCGCTTTGTTTATGCCGGATGCGGCGTTAACGCCTTATCCGGCCTACAAAATCGTGCAAATTCAATACATTGCAg                                                    <  2:2073775/100‑1 (MQ=255)
                                                       gggTATACAGATTGTTGACAAAGTGCGCTTTGTTTATGCCGGATGCGGCGTTAACGCCTTATCCGGCCTACAAAATCGTGCAAATTCAATACATTGCAGa                                                   <  2:3612945/100‑1 (MQ=255)
                                                       gggTATACAGATTGTTGACAAAGTGCGCTTTGTTTATGCCGGATGCGGCGTTAACGCCTTATCCGCCCTACAAAATCGTGCAAATTCAATACATTGCAgc                                                   >  1:2570598/1‑99 (MQ=255)
                                                           ataCAGATTGTTGACAAAGTGCGCTTTGTTTATGCCGGATGCGGCGTTAACGCCTTATCCGGCCTACAAAATCGTGCAAATTCAATACATTGCAGAAt‑‑tt                                            >  1:2054344/1‑100 (MQ=255)
                                                            taCAGATTGTTGACAAAGTGCGCTTTGTTTATGCCGGATGCGGCGTTAACGCCTTATCCGGCCTACAAAATCGTGCAAATTCAATACATTGCAGAAt‑‑ttt                                           <  2:277542/100‑1 (MQ=255)
                                                              cAGATTGTTGACAAAGTGCGCTTTGTTTATGCCGGATGCGGCGTTAACGCCTTATCCGGCCTACAAAATCGTGCAAATTCAATACATTGCAGAAT‑‑TTTC‑‑g                                       <  2:410958/100‑1 (MQ=255)
                                                              cAGATTGTTGACAAAGTGCGCTTTGTTTATGCCGGATGCGGCGTTAACGCCTTATCCGGCCTACAAAATCGTGCAAATTCAATACATTGCAGAAT‑‑TTTC‑‑g                                       <  1:660342/100‑1 (MQ=255)
                                                               aGATTGTTGACAAAGTGCGCTTTGTTTATGCCGGATGCGGCGTTAACGCCTTATCCGGCCTACAAAATCGTGCAAATTCAATACATTGCAGAAT‑‑TTTC‑‑Gt                                      <  2:3275461/100‑1 (MQ=255)
                                                               aGATTGTTGACAAAGTGCGCTTTGTTTATGCCGGATGCGGCGTTAACGCCTTATCCGGCCTACAAAATCGTGCAAATTCAATACATTGCAGAAT‑‑TTTC‑‑Gt                                      <  2:2096349/100‑1 (MQ=255)
                                                               aGATTGTTGAATAAGTGCGCTTTGTTTATGCCGGATTCGGCGTTAACGCCTTATCCGGACTACAAAATCGTGCAAATTCAATACCTTGCAGAAT‑‑TTTC‑‑Gt                                      <  1:1663792/100‑1 (MQ=17)
                                                                 aTTGTTGACAAAGTGCGCTTTGTTTATGCCGGATGCGGCGTTAACGCCTTATCCGGCCTACAAAATCGTGCAAATTCAATACATTGCAGAAT‑‑TTTC‑‑GTAg                                    <  2:1525862/100‑1 (MQ=255)
                                                                  gcataaaCAAAGCGCACTTTGTTTATGCCGGATGCGGCGTTAACGCCTTATCCGGCCTACAAAATCGTGCAAATTCAATACATTGCAGTAT‑‑TTTC‑‑GTAgg                                   >  1:861013/7‑100 (MQ=17)
                                                                    gttgACAAAGTGCGCTTTGTTTATGCCGGATGCGGCGTTAACGCCTTATCCGGCCTACAAAATCGTGCAAATTCAATACATTGCAGAAT‑‑TTTC‑‑GTAGGcc                                 >  2:1630060/1‑100 (MQ=255)
                                                                      tgACAAAGTGAACTTTGTTTATGCCGGATGCGGCATTAACGCCTTATCCGGCCTACAAAATCGTGAAAATTCAATACATTGCAGAAT‑‑TTTC‑‑GTAgcc                                  >  1:1433280/1‑95 (MQ=17)
                                                                       gACAAAGTGCGCTTTGTTTATGCCGGATGCGGCGTTAACGCCTTATCCGGCCTACAAAATCGTGCAAATTCAATACATTGCAGAAT‑‑TTTC‑‑GTAGAGCTGa                              >  1:3589986/1‑100 (MQ=255)
                                                                         cAAAGTGCGCTTGGTTTATGCCGGATGCGGCGTTCACGTCTTATCCGGCCTACAAAATCGTGCAAATTCAATACATTGCAGAAT‑‑TTTC‑‑GTAGGCCTGATa                            <  2:74374/100‑1 (MQ=255)
                                                                                 gcTTTGTTTATGCCGGATGCGGCGTTAACGCCTTATCCGGCCTACAAAATCGTGCAAATTCAATACATTGCAGAAT‑‑TTTC‑‑GTAGGCCTGATAAGCGGAgc                    >  2:1152879/1‑100 (MQ=255)
                                                                                  cTTTGTTTATGCCGGATGCGGCGTTAACGCCTTATCCGGCCTACAAAATCGTGCAAATTCAATACATTGCAGAAT‑‑TTTC‑‑GTAGGCCTGATAAGCGGAgcg                   <  2:2237953/100‑1 (MQ=255)
                                                                                   tttttgtaggCCGGATAAGGCGTTCACGCCGCATCCGGCATACAAAAACATGCAAATTCAATAGATTGCAGAGAT‑‑TAT‑‑GTAGGCCTGATAAGCGTAgcgc                  >  2:2048626/10‑100 (MQ=17)
                                                                                     tGTTTATGCCGGATGCGGCGTTAACGCCTTATCCGGCCTACAAAATCGTGCAAATTCAATACATTGCAGAAT‑‑TTTC‑‑GTAGGCCTGATAAGCGGAGCGCAt                >  1:3497688/1‑100 (MQ=255)
                                                                                      ttgtaggCCGGATAAGGCGTTTACGCCGCATCCGGCCTACAAAAATGTGCAAATTCAATAAATTGCAATTCAACTT‑‑GTAGGCCTGATAAGCGCAGCGCAt                >  1:3164105/7‑100 (MQ=18)
                                                                                      gTTTATGCCGGATGCGGCGTTAACGCCTTATCCGGCCTACAAAATCGTGCAAATTCAATACATTGCAGAAT‑‑TTTC‑‑GTAGGCCTGATAAGCGGAGCGCATc               >  2:3066254/1‑100 (MQ=255)
                                                                                       tttATGCCGGATGCGGCGTTAACGCCTTATCCGGCCTACAAAATCGTGCAAATTCAATACATTGCAGAAT‑‑TTTC‑‑GTAGGCCTGATAAGCGGAGCGCATCa              >  2:1459496/1‑100 (MQ=255)
                                                                                        ttATGCCGGATGCGGCGTTACCGCCTTATCCGGCCTACAAAATCGTGCAAATTCAATACATTGCAGAAT‑‑TTTC‑‑GTAGGCCTGATAAGCGGAgagatcgga             >  1:2254908/1‑92 (MQ=18)
                                                                                        gtaggCCGGATAAGGCGTTTACGCCGCATCCGGCCTACAAAAATGTGCAAATTCAATAAATTGCAATTCAACTT‑‑GTAGGCCTGATAAGCGCAGCGCATCa              >  2:873352/5‑100 (MQ=18)
                                                                                         tATGCCGGATGCGGCGTTAACGCCTTATCCGGCCTACAAAATCGTGCAAATTCAATACATTGCAGAAT‑‑TTTC‑‑GTAGGCCTGATAAGCGGAGCGCATCAg             <  1:1409824/99‑1 (MQ=255)
                                                                                          aggCCGGATAAGGCGTTTACGCCGCATCCGGCCTACAAAAGCATGCAAATTCAATATATTGCAGAGT‑‑TCAC‑‑GTAGGCCTGATAAGCGTAGCGCATCAGGc           >  1:1558833/3‑100 (MQ=17)
                                                                                          aggCCGGATAAGGCGTTTACGCCGCATCCGGCATACAAAAGCGTGCAAATTCAATACATTGCATGGGC‑‑CAT‑‑GTAGGCCTGATAAGCGTAGCGCATCAGGc           >  1:2681016/3‑100 (MQ=17)
                                                                                          aggCCGGATAAGGCGTTCACGCCGCATCCGGCCTACAAGATCGTGCAAATTCAATATATTGCATGAG‑‑CCAC‑‑GTAGGCCTGATAAGCGTAGCGCATCAgg            >  2:3587946/3‑99 (MQ=17)
                                                                                          aggCCGGATAAGGCGTTCACGCCGCATCCGGCCTACAAAAACATGCAAATTCAATAGATTGCAGAGAT‑‑TAT‑‑GTAGGCCTGATAAGCGTAGCGCATCAGGc           >  1:1026973/3‑100 (MQ=17)
                                                                                          aTGCCGGATGCGGCGTTAACGCCTTATCCGGCCTACAAAATCGTGCAAATTCAATACATTGCAGAAT‑‑TTTC‑‑GTAGGCCTGATAAGCGGAGCGCATCAGGc           >  2:2466547/1‑100 (MQ=255)
                                                                                          aTGCCGGATGCGGCGTTAACGCCTTATCCGGCCTACAAAATCGTGCAAATTCAATACATTGCAGAAT‑‑TTTC‑‑GTAGGCCTGATAAGCGGAGCGCATCAGGc           >  1:750389/1‑100 (MQ=255)
                                                                                           tGCCGGATGCGGCGTTAACGCCTTATCCGGCCTACAAAATCGTGCAAATTCAATACATTGCAGAAT‑‑TTTC‑‑GTAGGCCTGATAAGCGGAGCGCATCAGGCa          >  2:3435657/1‑100 (MQ=255)
                                                                                           ggCCGGATAAGGCGTTTACGCCGCATCCGGCCTACAAATTCTGCAAATTCAACAAATTGCGAATCCCTT‑‑GTAGGCCTGATAAGCGTAGCGCATCAGGCAt         >  1:1062888/2‑100 (MQ=18)
                                                                                           ggCCGGATAAGGCGTTTACGCCGCATCCGGCCTACAAATTCTGCAAATTCAACAAATTGCGAATCCCTT‑‑GTAGGCCTGATAAGCGTAGCGCATCAGGCAt         >  2:2265248/2‑100 (MQ=18)
                                                                                           ggCCGGATAAGGCGTTCACGCCGCATCCGGCCTACAAAATTGTGTAAATTCAATAAATTGATCGAT‑‑CCAC‑‑GTAGGCCTGATAAGCGTAGCGCATCAGGCa          >  1:1284612/2‑100 (MQ=21)
                                                                                           ggCCGGATAAGGCGTTCACGCCGCATCCGGCCTACAAAATCATGAAAATTCAATACATTGCAAGAT‑‑TTTC‑‑GTAGGCCTGATAAGCGTAGCGCATCAGGCa          >  1:2945391/2‑100 (MQ=14)
                                                                                           ggCCGGATAAGGCGTTCACGCCGCATCCGGCCTACAAAACTGTGCAAATTCAATATATTGCAGGAA‑‑ACAC‑‑GTAGGCCTGATAAGCGAAGCGCATCAGGCa          >  1:3572448/2‑100 (MQ=18)
                                                                                           ggCCGGATAAGGCGTTCACGCCGCATCCGGCCTACAAAACCTTGCTAATTCAATATATTGCAGGGA‑‑CTAT‑‑GTAGGCCTGATAAGCATAGCGCATCAGGCa          >  1:1164148/2‑100 (MQ=18)
                                                                                           ggCCGGATAAGGCGTTCACGCCGCATCCGGCCTACAAAAACATGCAAATTCAATAGATTGCAGAGAT‑‑TAT‑‑GTAGGCCTGATAAGCGTAGCGCATCAGGCa          >  2:2084962/2‑100 (MQ=17)
                                                                                           ggCCGGATAAGGCGTTCACGCCGCATCCGGCCTACAAATCATGCAAATCCAAAATATTGCAGGAG‑‑TTAT‑‑GTAGGCCTGATAAGCGTAGCGCATCAGGCa          >  1:3187227/2‑99 (MQ=18)
                                                                                            gCCGGATAAGGCGTTTACGCCGCATCCGGCCTACAAAATCGTGCAAATTCAACATATTGCAAT‑‑TCTCTTGTAGGCCTGATAAGCGTAGCGCATCAGGctg         >  1:2134425/1‑98 (MQ=18)
                                                                                            gCCGGATAAGGCGTTTACGCCGCATCCGGCCTACAAAATCGTGCAAATTCAACATATTGCAAT‑‑TCTCTTGTAGGCCTGATAAGCGTAGCGCATCAGGctg         >  2:350184/1‑98 (MQ=18)
                                                                                            gCCGGATAAGGCGTTTACGCCGCATCCGGCCTACAAAAGCGTGCAAATTCAATACATTGCATGGG‑‑CCAT‑‑GTAGGCCTGATAAGCGTAGCGCATCAGGCaa         >  2:913008/1‑99 (MQ=17)
                                                                                            gCCGGATAAGGCGTTCACGCCGCATCCGGCCTACAAAATCATGAAAATTCAATACATTGCAAGAT‑‑TTTC‑‑GTAGGCCTGATAAGCGTAGCGCATCAGGCac         >  2:1199817/1‑99 (MQ=17)
                                                                                            gCCGGATAAGGCGTTCACGCCGCATCCGGCCTACAAAAGCATGAAAATTCAATATATTGCAGGAG‑‑CTGC‑‑GTAGGCCTGATAAGCGTAGCGCATCAGGCag         >  2:452230/1‑99 (MQ=18)
                                                                                            gCCGGATAAGGCGTTCACGCCGCATCCGGCCTACAAAAGCATGAAAATTCAATATATTGCAGGAG‑‑CTGC‑‑GTAGGCCTGATAAGCGTAGCGCATCAGGCag         >  2:1328885/1‑99 (MQ=17)
                                                                                            gCCGGATAAGGCGTTCACGCCGCATCCGGCCTACAAATCATGCAAATCCAAAATATTGCAGGAG‑‑TTAT‑‑GTAGGCCTGATAAGCGTAGCGCATCAGGCAtt        >  2:2698085/1‑100 (MQ=21)
                                                                                             ccGGATGCGGCGTTAACGCCTTATCCGGCCTACAAAATCGTGCAAATTCAATACATTGCAGAAT‑‑TTTC‑‑GTAGGCCTGATAAGCGGAGCGCATCAGGCAtt        >  2:3079328/1‑100 (MQ=255)
                                                                                                  tGCGGCGTTAACGCCTTATCCGGCCTACAAAATCGTGCAAATTCAATACATTGCAGAAT‑‑TTTC‑‑GTAGGCCTGATAAGCGGAGCGCATCAGGCATTtttgc   <  1:2659298/100‑1 (MQ=255)
                                                                                                   gcggcgTTAACGCCTTATCCGGCCTACAAAATCGTGCAAATTCAATACATTGCAGAAT‑‑TTTC‑‑GTAGGCCTGATAAGCGGAGCGCATCAGGCATTtttgcc  <  1:970168/100‑1 (MQ=255)
                                                                                                   |                                                                                                         
CAGGGAATGAAGTATTGCTGCCGGATGGCGCTGCACGCAGAGGCACAGCGCCACCGGGTATACAGATTGTTGACAAAGTGCGCTTTGTTTATGCCGGATGCGGCGTTAACGCCTTATCCGGCCTACAAAATCGTGCAAATTCAATACATTGCAGAAT‑‑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: