Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 3,705,984 T→G 35.8% intergenic (‑194/+114) 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,9840TG35.8% 346.9 / 290.4 242intergenic (‑194/+114)dppB/dppAdipeptide/heme ABC transporter permease/dipeptide/heme ABC transporter periplasmic binding protein; dipeptide chemotaxis receptor
Reads supporting (aligned to +/- strand):  ref base T (78/77);  new base G (74/12);  total (153/89)
Fisher's exact test for biased strand distribution p-value = 1.57e-08
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 4.69e-02

Only 100 of 246 total aligned reads displayed.

TGCTGCCGGATGGCGCTGCACGCAGAGGCACAGCGCCACCGGGTATACAGATTGTTGACAAAGTGCGCTTTGTTTATGCCGGATGCGGCGTTAACGCCTTATCCGGCCTAC‑AAAATCGTGC‑AAATTCAATACATTGC‑‑AGA‑AT‑TTTC‑‑GTAGGCCTGATAAGCGGAGCGCATCAGGCATTTTTGCCTTTGCCATCAGTCT  >  NC_000913/3705886‑3706083
                                                                                                  |                                                                                                           
tgctgcCGGATGGCGCTGCACGCAGAGGCACAGCGCCACCGGGTATACAGATTGTTGACAAAGTGCGCTTTGTTTATGCCGGATGCGGCGTTAACGCCtt                                                                                                            <  2:2874235/100‑1 (MQ=255)
 gctgcCGGATGGCGCTGCACGCAGAGGCACAGCGCCACCGGGTATACAGATTGTTGACAAAGTGCGCTTTGTTTATGCCGGATGCGGCGTTAACGCCTTa                                                                                                           >  1:721500/1‑100 (MQ=255)
 gctgcCGGATGGCGCTGCACGCAGAGGCACAGCGCCACCGGGTATACAGATTGTTGACAAAGTGCGCTTTGTTTATGCCGGATGCGGCGTTAACGCCTTa                                                                                                           >  1:2170085/1‑100 (MQ=255)
  ctgcCGGATGGCGCTGCACGCAGAGGCACAGCGCCACCGGGTATACAGATTGTTGACAAAGTGCGCTTTGTTTATGCCGGATGCGGCGTTAACGCCTTAt                                                                                                          >  1:1994346/1‑100 (MQ=255)
  ctgcCGGATGGCGCTGCACGCAGAGGCACAGCGCCACCGGGTATACAGATTGTTGACAAAGTGCGCTTTGTTTATGCCGGATGCGGCGTTAACGCCTTAt                                                                                                          >  1:1952097/1‑100 (MQ=255)
    gcCGGATGGCGCTGCACGCAGAGGCACAGCGCCACCGGGTATACAGATTGTTGACAAAGTGCGCTTTGTTTATGCCGGATGCGGCGTTAACGCCTTATcc                                                                                                        <  2:830180/100‑1 (MQ=255)
    gcCGGATGGCGCTGCACGCAGAGGCACAGCGCCACCGGGTATACAGATTGTTGACAAAGTGCGCTTTGTTTATGCCGGATGCGGCGTTAACGCCTTATcc                                                                                                        >  2:590897/1‑100 (MQ=255)
     ccGGATGGCGCTGCACGCAGAGGCACAGCGCCACCGGGTATACAGATTGTTGACAAAGTGCGCTTTGTTTATGCCGGATGCGGCGTTAACGCCTTATCCg                                                                                                       <  1:3371274/100‑1 (MQ=255)
      cGGATGGCGCTGCACGCAGAGGCACAGCGCCACCGGGTATACAGATTGTTGACAAAGTGCGCTTTGTTTATGCCGGATGCGGCGTTAACGCCTTATCCgg                                                                                                      <  2:2024779/100‑1 (MQ=255)
      cGGATGGCGCTGCACGCAGAGGCACAGCGCCACCGGGTATACAGATTGTTGACAAAGTGCGCTTTGTTTATGCCGGATAAGGCGTTAACGCCGCATCCgg                                                                                                      <  2:2199478/100‑1 (MQ=37)
       ggATGGCGCTGCACGCAGAGGCACAGCGCTACCGGGTATACAGATTGTTGACAAAGTGCGCTTTGTTTATGCCGGATAAGGCGTTAACGCCGCATCCGGc                                                                                                     <  1:367745/100‑1 (MQ=25)
       ggATGGCGCTGCACGCAGAGGCACAGCGCCACCGGGTATACAGATTGTTGACAAAGTGCGCTTTGTTTATGCCGGATAAGGCGTTAACGCCGCATCCGGc                                                                                                     <  2:654093/100‑1 (MQ=37)
        gATGGCGCGGCACGCAGAGGCACAGCGCCACCGGGTATACAGATTGTTGACAAAGTGCGCTGTGGTTATGCCGGATGCGGCGTTATCGCCTTATCCGGcc                                                                                                    <  1:1279210/100‑1 (MQ=25)
          tGGCGCTGCACGCAGAGGCACAGCGCCACCGGGTATACAGATTGTTGACAAAGTGCGCTTTGTTTATGCCGGATAAGGCGTTAACGCCGCATCCGGCATa                                                                                                  <  1:1804698/100‑1 (MQ=37)
          tGGCGCTGCACGCAGAGGCACAGCGCCACCGGGTATACAGATTGTTGACAAAGTGCGCTTTGTTTATGCCGGATAAGGCGTTAACGCCGCATCCGGCATa                                                                                                  <  2:3423502/100‑1 (MQ=37)
            gcgcTGCACGCAGAGGCACAGCGTCACCGGGTATACAGATTGTTGACAAAGTGCGCTTTGTTTATGCCGGATGCGGCGTTAACGCCTTATCCGGCCTAC‑a                                                                                               <  1:511058/100‑1 (MQ=255)
            gcgcTGCACGCAGAGGCACAGCGCCACCGGGTATACAGATTGTTGACAAAGTGCGCTTTGTTTATGCCGGATGCGGCGTTAACGCCTTATCCAGCCTAC‑a                                                                                               >  1:363104/1‑100 (MQ=255)
              gcTGCACGCAGAGGCACAGCGCCACCGGGTATACAGATTGTTGACAAAGTGCGCTTTGTTTATGCCGGATGCTGCGTTAACGCCTTATCCGGCCTAC‑aaa                                                                                             >  1:935600/1‑100 (MQ=255)
              gcTGCACGCAGAGGCACAGCGCCACCGGGTATACAGATTGTTGACAAAGTGCGCTTTGTTTATGCCGGATGCGGCGTTAACGCCTTATCCGGCCTAC‑aaa                                                                                             >  2:3026780/1‑100 (MQ=255)
                 gCACGCAGAGGCACGGCGCCACCGGGTATACAGATTGTTGACAAAGTGCGCTTTGTTTATGCCGGATGCGGCGTTAACGCCTTATCCGGCCTAC‑AAAATc                                                                                          >  1:2716316/1‑100 (MQ=255)
                  cACGCAGAGGCACAGCGCCACCGGGTATACAGATTGTTGACAAAGTGCGCTTTGTTTATGCCGGATGCGGCGTTAACGCCTTATCCGGCCTAC‑AAAATCg                                                                                         <  2:2958412/100‑1 (MQ=255)
                  cACGCAGAGGCACAGCGCCACCGGGTATACAGATTGTTGACAAAGTGCGCTTTGTTTATGCCGGATGCGGCGTTAACGCCTTATCCGGCCTAC‑AAAATCg                                                                                         <  2:573253/100‑1 (MQ=255)
                    cGCAGAGGCACAGCGCCACCGGGTATACAGATTGTTGACAAAGTGCGCTTTGTTTATGCCGGATGCGGCGTTAACGCCTTATCCGGCCTAC‑AAAATCGTg                                                                                       <  2:2630196/100‑1 (MQ=255)
                    cGCAGAGGCACAGCGCCACCGGGTATACAGATTGTTGACAAAGTGCGCTTTGTTTATGCCGGATGCGGCGTTAACGCCTTATCCGGCCTAC‑AAAATCGTg                                                                                       <  1:2450053/100‑1 (MQ=255)
                     gCAGAGGCACAGCGCCACCGGGTATACAGATTGTTGACAAAGTGCGCTTTGTTTATGCCGGATGCGGCGTTAACGCCTTATCCGGCCTAC‑AAAATCGTGc                                                                                      <  2:1254493/100‑1 (MQ=255)
                        gagGCACAGCGCCACCGGGTATACAGATTGTTGACAAAGTGCGCTTTGTTTATGCCGGATGCGGCGTTAACGCCTTATCCGGCCTAC‑AAAATCGTGC‑caa                                                                                  >  1:798010/1‑100 (MQ=255)
                            cacaGCGCCACCGGGTATACAGATTGTTGACAAAGTGCGCTTTGTTTATGCCGGATGCGGCGTTAACGCCTTATCCGGCCTAC‑AAAATCGTGC‑AAATTCa                                                                              >  1:2202194/1‑100 (MQ=255)
                              caGCGTCACCGGGTATACAGATTGTTGACAAAGTGCGCTTTGTTTATGCCGGATGCGGCGTTAACGCCTTATCCGGCCTAC‑AAAATCGTGC‑AAATTCAAt                                                                            <  2:1266067/100‑1 (MQ=255)
                               aGCGCCACCGGGTATACAGATTGTTGACAAAGTGCGCTTTGTTTATGCCGGATGCGGCGTTAACGCCTTATCCGGCCTAC‑AAAATCGTGC‑AAATTCAATa                                                                           <  2:1920368/100‑1 (MQ=255)
                                   ccACCGGGTATACAGATTGTTGACAAAGTGCGCTTTGTTTATGCCGGATGCGGCGTTAACGCCTTATCCGGCCTAC‑AAAATCGTGC‑AAATTCAATACAtt                                                                       <  2:2182920/100‑1 (MQ=255)
                                   ccACCGGGTATACAGATTGTTGACAAAGTGCGCTTTGTTTATGACGGATGCGGCGTTAACGCCTTATTCGGCCTAC‑AAAATCGTGC‑AAATTCAATACAtt                                                                       >  1:1894181/1‑100 (MQ=255)
                                    cACCGGGTATACAGATTGTTGACAAAGTTCGTTTTGTTTATGCCGGATGCGGCGTTAACGCCTTACCCGGCCTAC‑TAAATCGTGC‑AAATTCCATACAtgg                                                                      <  1:1386834/100‑3 (MQ=17)
                                     aCCGGGTATACACCTTGTTGACAAAGTCCGCTTTGTTTATCTCGGATGCGTCCTTAACGCCTTATCCGTCCTAC‑AAAATCGTGC‑AAATTCATTACATTGc                                                                     <  1:313699/100‑1 (MQ=17)
                                             taCAGATTGTTGACAAAGTGCGCTTTGTTTATGCCGGATGCGGCGTTAACGCCTTATCCGGCCTAC‑AAAATCGTGC‑AAATTCAATACATTGC‑‑AGA‑At‑ttt                                                         >  1:2480225/1‑100 (MQ=255)
                                              aCAGATTGTTGACAAAGTGCGCTTTGTTTATGCCGGATGCGGCGTTAACGCCTTATCCGGCCTAC‑AAAATCGTGC‑AAATTCAATACATTGC‑‑AGA‑AT‑TTTc                                                        <  1:3467121/100‑1 (MQ=255)
                                               cAGATTGTTGACAAAGTGCGCTTTGTTTATGCCGGATGCGGCGTTAACGCCTTATCCGGCCTAC‑AAAATCGTGC‑AAATTCAATACATTGC‑‑AGA‑AT‑TTTC‑‑g                                                     <  2:1989021/100‑1 (MQ=255)
                                               cAGATTGTTGACAAAGTGCGCTTTGTTTATGCCGGATGCGGCGTTAACGCCTTATCCGGCCTAC‑AAAATCGTGC‑AAATTCAATACATTGC‑‑AGA‑AT‑TTTC‑‑g                                                     <  2:2898213/100‑1 (MQ=255)
                                               cAGATTGTTGACAAAGTGCGCTTTGTTTATGCCGGATGCGGCGTTAACGCCTTATCCGGCCTAC‑AAAATCGTGC‑AAATTCAATACATTGC‑‑AGA‑AT‑TTTC‑‑g                                                     <  1:2800359/100‑1 (MQ=255)
                                               cAGATTGTTGACAAAGTGCGCTTTGTTTATGCCGGATGCGGCGTTAACGCCTTATCCGGCCTAC‑AAAATCGTGC‑AAATTCAATACATTGC‑‑AGA‑AT‑TTTC‑‑g                                                     <  2:492490/100‑1 (MQ=255)
                                                   ttgttgACAAAGTGCGCTTTGTTTATGCCGGATGCGGCGTTAACGCCTTATCCGGCCTAC‑AAAATCGTGC‑AAATTCAATACATTGC‑‑AGA‑AT‑TTTC‑‑GTAgg                                                 <  1:1917468/100‑1 (MQ=255)
                                                   ttgttgACAAAGTGCGCTTTGTTTATGCCGGATGCGGCGTTAACGCCTTATCCGGCCTAC‑AAAATCGTGC‑AAATTCAATACATTGC‑‑AGA‑AT‑TTTC‑‑GTAgg                                                 >  1:3510846/1‑100 (MQ=255)
                                                    tgttgACAAAGTGCGCTTTGTTTATGCCGGATGCGGCGTTAACGCCTTATCCGGCCTAC‑AAAATCGTGC‑AAATTCAATACATTGC‑‑AGA‑AT‑TTTC‑‑GTAGGc                                                <  2:2242086/100‑1 (MQ=255)
                                                     gttgACAAAGTGCGCTTTGTTTATGCCGGATGCGGCGTTAACGCCTTATCCGGCCTAC‑AAAATCGTGC‑AAATTCAATACATTGC‑‑AGA‑AT‑TTTC‑‑GTAGGcc                                               >  2:161064/1‑100 (MQ=255)
                                                         aCAAAGTGCGCTTTGTTTATGCCGGATGCGGCGTTAACGCCTTATCCGGCCTAC‑AAAATCGTGC‑AAATTCAATACATTGC‑‑AGA‑AT‑TTTC‑‑GTAGGCCTGAt                                           <  2:1501042/100‑1 (MQ=255)
                                                          cAAAGTGCGCTTTGTTTATGCCGGATGCGGCGTTAACGCCTTATCCGGCCTAC‑AAAATCGTGC‑AAATTCAATACATTGC‑‑AGA‑AT‑TTTC‑‑GTAGGCCTGATa                                          >  1:3196531/1‑100 (MQ=255)
                                                          cAAAGTGCGCTTTGTTTATGCCGGATGCGGCGTTAACGCCTTATCCGGCCTAC‑AAAATAGTGC‑AAGTTCACTACATTGC‑‑AGA‑CA‑CTTC‑‑GTAGGCCTg                                             >  2:3182428/1‑97 (MQ=14)
                                                             aGTGCGCTTTGTTTATGCCGGATGCGGCGTTAACGCCTTATCCGGCCTAC‑AAAATCGTGC‑AAATTCAATACATTGC‑‑AGA‑AT‑TTTC‑‑GTAGGCCTGATAAGc                                       <  1:3751427/100‑1 (MQ=255)
                                                               tGCGCGTTGTTTATGCCGGATGCGGCGTTAACGCCTTATCCGGCCTAC‑AAAATCGTGC‑AAATTCAATACATTGC‑‑AGA‑AT‑TTTC‑‑GTAGGCCTGATAAGCgg                                     <  1:1307701/100‑1 (MQ=255)
                                                                 cgcTTTGTTTATGCCGGATGCGGCGTTAACGCCTTATCCGGCCTAC‑AAAATCGTGC‑AAATTCAATACATTGC‑‑AGA‑AT‑TTTC‑‑GTAGGCCTGATAAGCGGAg                                   <  1:2709829/100‑1 (MQ=255)
                                                                  gcTTTGTTTATGCCGGATGCGGCGTTAACGCCTTATCCGTCCTAC‑AAAATCGTGC‑AAATTCAATACATTGC‑‑AGA‑AT‑TTTC‑‑GTAGGCCTGATAAGCGGAgc                                  >  1:2611423/1‑100 (MQ=255)
                                                                  gcTTTGTTTATGCCGGATGCGGCGTTAACGCCTTATCCGGCCTAC‑AAAATCGTGCAAAATTCAATACATTGC‑‑AGA‑AT‑TTTC‑‑GTAGGCCTGTTAAGCg                                      >  2:2423031/1‑97 (MQ=255)
                                                                  gcTTTGTTTATGCCGGATGCGGCGTTAACGCCTTATCCGGCCTAC‑AAAATCGTGC‑AAATTCAATACATTGC‑‑AGA‑AT‑TTTC‑‑GTAGGCCTGATAAGCGGa                                    <  1:2325865/98‑1 (MQ=255)
                                                                  gcTTTGTTTATGCCGGATGCGGCGTTAACGCCTTATCCGGCCTAC‑AAAATCGTGC‑AAATTCAATACATTGC‑‑AGA‑AT‑TTTC‑‑GTAGGCCTGATAAGCGGAgc                                  >  2:2063265/1‑100 (MQ=255)
                                                                     ttGTTTATGCCGGATGCGGCGTTAACGCCTTATCCGGCCTAC‑AAAATCGTGC‑AAATTCAATACATTGC‑‑AGA‑AT‑TTTC‑‑GTAGGCCTGATAAGCGGAGCGCa                               <  2:1240099/100‑1 (MQ=255)
                                                                     ttGTTTATGCCGGATGCGGCGTTAACGCCTTATCCGGCCTAC‑AAAATCGTGC‑AAATTCAATACATTGC‑‑AGA‑AT‑TTTC‑‑GTAGGCCTGATAAGCGGAGCGCa                               <  1:838261/100‑1 (MQ=255)
                                                                      tttgtaggCCGGATAAGGCGTTCACGCCGCATCCGGCCTAC‑AAAGTCGTAT‑AAATTCAATGTATTGC‑‑AGT‑TT‑ATTG‑‑GTAGGCCTGATAAGCGCAGCGCAt                              >  2:2698219/8‑100 (MQ=21)
                                                                       ttgtaggCCGGATAAGGCGTTCACGCCGCATCCGGCCTAC‑AAAACTTTGC‑AAATTCAATATATTGC‑‑‑‑A‑AT‑CTCCGTGTAGGCCTGATAAGCGTAGCGCATc                             >  2:1720692/7‑100 (MQ=18)
                                                                         ttATGCCGGATGCGGCGTTAACGCCTTATCCGGCCTAC‑AAAATCGTGC‑AAATTCAATACATTGC‑‑AGA‑AT‑TTTC‑‑GTAGGCCTGATAAGCGGAGCGCATCAg                           <  2:2716330/100‑1 (MQ=255)
                                                                           nggCCGGATAAGGCGTTCACGCCGCATCCGGCCTAC‑‑ATACCGCGC‑AAATTCAATATATTGCGGAGA‑AA‑AT‑‑‑‑GTAGGCCTGATAAGCGTAGCGCGTCAGGct                        >  1:2043341/3‑99 (MQ=18)
                                                                           aggCCGGATAAGGCGTTCACGCCGCATCCGGCCTAC‑AAAATCGCTT‑AAATTCAATATATTGC‑‑ACC‑AT‑CCTC‑‑GTAGGCCTGATAAGCGTAGCGCATCAgg                          >  1:1429532/3‑99 (MQ=18)
                                                                           aTGCCGGATGCGGCGTTAACGCCGTATCCGGCCTAC‑AAAATCGTGC‑AAATTCAATACATTGC‑‑AGA‑AT‑TTTC‑‑GTAGGCCTGATAAGCGGAGCGCATCAgg                          <  1:1470/99‑1 (MQ=255)
                                                                            tGCCGGATGCGGCGTTAACGCCTTATCCGGCCTAC‑AAAATCGTGC‑AAATTCAATACATTGC‑‑AGA‑AT‑TTTC‑‑GTAGGCCTGATAAGCGGAGCGCATCAGGCa                        <  2:406035/100‑1 (MQ=255)
                                                                            tGCCGGATGCGGCGTTAACGCCTTATCCGGCCTAC‑AAAATCGTGC‑AAATTCAATACATTGC‑‑AGA‑AT‑TTTC‑‑GTAGGCCTGATAAGCGGAGCGCATCAGGCa                        >  1:1379054/1‑100 (MQ=255)
                                                                            tGCCGGATGCGGCGTTAACGCCTTATCCGGCCTAC‑AAAATCGTGC‑AAATTCAATACATTGC‑‑AGA‑AT‑TTTC‑‑GTAGGCCTGATAAGCGGAGCGCATCAGGCa                        >  2:1530808/1‑100 (MQ=255)
                                                                            ggCCGGATAAGGCGTTTACGCCGCATCCGGCCTAC‑AAAAGCATGC‑AAATTCAATATATTGC‑‑AGA‑GT‑TCAC‑‑GTAGGCCTGATAAGCGTAGCGCATCAGGCa                        >  1:386876/2‑100 (MQ=14)
                                                                            ggCCGGATAAGGCGTTCACGCCGCATCCGGCCTACAAAAATCTTGC‑CAATTCAATATATTGC‑‑AGG‑AC‑CAT‑‑‑GTAGGCCTGATAAGCGTAGCGCATCAGGCa                        >  2:1389603/2‑100 (MQ=18)
                                                                             gCCGGATGCGGCGTTAACGCCTTATCCGGCCTAC‑AAAATCGTGC‑AAATTCAATACATTGC‑‑AGA‑AT‑TTTC‑‑GTAGGCCTGATAAGCGGAGCGCATCAGGCAt                       >  2:1460892/1‑100 (MQ=255)
                                                                             gCCGGATAAGGCGTTTACGCCGCATCTGGCCTAC‑AAATTCTTCC‑AAATTCAATATATTGC‑‑AGA‑AA‑TCAT‑‑GTAGGCCTGATAAGCGTAGCGCATCAGGCaa                       >  2:2975010/1‑99 (MQ=18)
                                                                             gCCGGATAAGGCGTTTACGCCGCATCCGGCCTAC‑AAAAATGTGC‑AAATTCAATAAATTGC‑‑A‑A‑AT‑CAACATGTAGGCCGGATAAGCGCAGCGCATCAGGCa                        >  1:2571710/1‑100 (MQ=18)
                                                                             gCCGGATAAGGCGTTCACGCCGCATCCGGCCTAC‑AAAACCTTGC‑TAATTCAATATATTGC‑‑AGG‑GA‑CTAT‑‑GTAGGCCTGATAAGCATAGCGCATCAGGCag                       >  1:1217899/1‑99 (MQ=18)
                                                                              ccggataaggcgTTTACGCCGCATCCGGCCTAC‑AAAACCGTGC‑TAATTCAATATATTGC‑‑AGA‑AA‑CCTT‑‑GTAGGCCTGATAAGCGTAGCGCATCAGGCagt                      >  1:2031860/9‑98 (MQ=17)
                                                                              ccggataaggcgTTCACGCCGCATCCGGCCTACAAAAATCTTGC‑CAATTCAATATATTGC‑‑AGG‑AC‑CAT‑‑‑GTAGGCCTGATAAGCGTAGCGCATCAGGCagt                      >  2:3307202/9‑98 (MQ=18)
                                                                              ccggataaggcgTTCACGCCGCATCCGGCCTAC‑AAAATTGTGT‑AAATTCAATAAATTGA‑‑TCG‑AT‑CCAC‑‑GTAGGCCTGATAAGCGTAGCGCATCAGGCagt                      >  2:105779/9‑98 (MQ=18)
                                                                              ccggataaggcgTTCACGCCGCATCCGGCCTAC‑AAAATCGCTT‑AAATTCAATATATTGC‑‑ACC‑AT‑CCTC‑‑GTAGGCCTGATAAGCGTAGCGCATCAGGCaat                      >  1:1152417/9‑98 (MQ=21)
                                                                              ccggataaggcgTTAACGCCGCATCCGGCCTAC‑AAAATCGTGC‑AAATTCAATACATTGC‑‑AGA‑AT‑TTTC‑‑GTAGGCCTGATAAGCGGAGCGCATCAGGCAtt                      >  2:1799584/9‑100 (MQ=18)
                                                                              ccGGATAAGGCGTTCACGCCGCATCCGGCCTAC‑AAAATTGTGT‑AAATTCAATAAATTGA‑‑TCG‑AT‑CCAC‑‑GTAGGCCTGATAAGCGTAGCGCATCAGGCagt                      >  1:3438258/1‑98 (MQ=18)
                                                                               cGGATGCGGCGTTAACGCCTTATCCGGCCTAC‑AAAATCGTGC‑AAATTCAATACATTGC‑‑AGA‑AT‑TTTC‑‑GTAGGCCTGATAAGCGGATCGCATCAGGCAtt                      >  2:1094665/1‑99 (MQ=255)
                                                                                ggataaggcgTTTACGCCGCATCCGGCCTAC‑AAATTC‑TGC‑AAATTCAACAAATTGC‑‑‑GA‑ATCCCTT‑‑GTAGGCCTGATAAGCGTAGCGCATCAGGCAtttgc                   >  1:1032296/7‑98 (MQ=18)
                                                                                ggataaggcgTTCACGCCGCATCCGGCCTACAAAAATCTTGC‑CAATTCAATATATTGC‑‑AGG‑AC‑CAT‑‑‑GTAGGCCTGATAAGCGTAGCGCATCAGGCAgttt                    >  1:2235449/7‑100 (MQ=18)
                                                                                ggataaggcgTTCACGCCGCATCCGGCCTAC‑AAAATCATGA‑AAATTCAATACATTGC‑‑AAG‑AT‑TTTC‑‑GTAGGCCTGATAAGCGTAGCGCATCAGGCacgct                    >  2:732943/7‑96 (MQ=17)
                                                                                ggataaggcgTTCACGCCGCATCCGGCCTAC‑AAAAGCATGA‑AAATTCAATATATTGC‑‑AGG‑AG‑CTGA‑‑GTAGGCCTGATAAGCGTAGCGCATCAGGCAgttt                    >  1:868898/7‑100 (MQ=17)
                                                                                ggataaggcgTTCACGCCGCATCCGGCCTAC‑AAAAACATGC‑AAATTCAATAGATTGC‑‑AGAGAT‑TAT‑‑‑GTAGGCCTGATAAGCGTAGCGCATCAGGCaatag                    >  1:2915465/7‑96 (MQ=17)
                                                                                 gataaggcgTTTACGCCGCATCCGGCCTAC‑AAAAGCATGC‑AAATTCAATATATTGC‑‑AGA‑GT‑TCAC‑‑GTAGGCCTGATAAGCGTAGCGCATCCGGCAgttt                    >  1:1774301/6‑99 (MQ=17)
                                                                                  ataaggcgTTCACGCCGCATCCGGCCTAC‑AAGATCGTGC‑AAATTCAATATATTGC‑‑ATG‑AG‑CCAC‑‑GTAGGCCTGATAAGCGTAGCGCATCAGGCAGTtttg                  <  2:2145583/96‑1 (MQ=17)
                                                                                  aTGCGGCGTTAACGCCTTATCCGGCCTAC‑AAAATCGTGC‑AAATTCAATACATTGC‑‑AGA‑AT‑TTTC‑‑GTAGGCCTGATAAGCGGAGCGCATCAGGCATTtttg                  <  1:2341624/100‑1 (MQ=255)
                                                                                  aTGCGGCGTTAACGCCTTATCCGGCCTAC‑AAAATCGTGC‑AAATTCAATACATTGC‑‑AGA‑AT‑TTTC‑‑GTAGGCCTGATAAGCGGAGCGCATCAGGCATTtttg                  <  1:3025232/100‑1 (MQ=255)
                                                                                     cggcgTTAACGCCTTATCCGGCCTAC‑AAAATCGTGC‑AAATTCAATACATTGC‑‑AGA‑AT‑TTTC‑‑GTAGGCCTGATAAGCGGAGCGCATCAGGCATTtttgcct               <  1:3742334/100‑1 (MQ=255)
                                                                                      ggcgTTCACGCCGCATCCGGCCTAC‑AAAAGCACGC‑AAACTCAATATATTGC‑‑AGAGAT‑CAT‑‑‑GTAGGCCTGATAAGCGTAGCGCATCAGGCAatttagcgtt              >  2:881092/1‑94 (MQ=17)
                                                                                      ggcgTTAACGCCTTATCCGGCCTAC‑AAAATCGTGC‑AAATTCAATACATTGC‑‑AGA‑AT‑TTTC‑‑GTAGGCCTGATAAGCGGAGCGCATCAGGCATTtttgcctt              >  2:3596102/1‑100 (MQ=255)
                                                                                         gTTAACGCCTTATCCGGCCTAC‑AAAATCGTGC‑AAATTCAATACATTGC‑‑AGA‑AT‑TTTC‑‑GTAGGCCTGATAAGCGGAGCGCATCAGGCATTtttgcctttgc           >  2:170255/1‑100 (MQ=255)
                                                                                          ttAACGCCTTATCCGGCCTAC‑AAAATCGTGC‑AAATTCAATACATTGC‑‑AGA‑AT‑TTTC‑‑GTAGGCCTGATAAGCGGAGCGCATCAGGCATTtttgcctttgcn          <  1:281638/100‑2 (MQ=255)
                                                                                            aaCGCCTTATCCGGCCTAC‑AAAATCGTGC‑AAATTCAATACATTGC‑‑AGA‑AT‑TTTC‑‑GTAGGCCTGATAAGCGGAGCGCATCAGGCATTTTTGCCTTTGCCAt        >  1:1838723/1‑100 (MQ=255)
                                                                                             aCGCCTTATCCGGCCTAC‑AAAATCGTGC‑AAATTCAATACATTGC‑‑AGA‑AT‑TTTC‑‑GTAGGCCTGATAAGCGGAGCGCATCAGGCATTTTTGCCTTTGCCATc       <  2:2451783/100‑1 (MQ=255)
                                                                                             aCGCCTTATCCGGACTAC‑AAAATCGTGA‑AAATTAAATACATTGC‑‑AGC‑AT‑TTTA‑‑GTAGGAGTGATAAGAGGAGCGCATCAGGCATTtttgcctttggtatc       >  1:2364465/1‑95 (MQ=18)
                                                                                              ccccTTATCCGGCCTAC‑AAAATCGTGC‑AAATTCAATACATTGC‑‑AGA‑AT‑TTTC‑‑GTAGGCCTGATAAGCGGAGCGCATCAGGCATTTTTGCCTTTGCCCTCa      >  2:2576738/3‑100 (MQ=255)
                                                                                              cGCCTTATCCGGCCTAC‑AAAATCGTGC‑AAATTCAATACATTGC‑‑AGA‑AT‑TTTC‑‑GTAGGCCTGATAAGCGGAGCGCATCAGGCATTTTTGCCTTTGCCATCa      >  1:2086127/1‑100 (MQ=255)
                                                                                              cGCCTTATCCGGCCTAC‑AAAATCGTGC‑AAATTCAATACATTGC‑‑AGA‑AT‑TTTC‑‑GTAGGACTGATAAGCGGAGCGCATCAGGCATTTTTGCCTTTGCCATCa      >  1:3549058/1‑100 (MQ=255)
                                                                                                ccTTATCCGGCCTAC‑AAAATCGTGC‑AAATTCAATACATTGC‑‑AGA‑AT‑TTTC‑‑GTAGGCCTGATAAGCGGAGCGCATCAGGCATTTTTGCCTTTGCCATCAGt    >  2:2224573/1‑100 (MQ=255)
                                                                                                 cTTATCCGGCCTAC‑AAAATCGTGC‑AAATTCAATACATTGC‑‑AGA‑AT‑TTTC‑‑GTAGGCCTGATAAGCGGAGCGCATCAGGCATTTTTGCCTTTGCCATCAGTc   <  2:1773257/100‑1 (MQ=255)
                                                                                                  ttATCCGGCCTAC‑AAAATCGTGC‑AAATTCAATACATTGC‑‑AGA‑AT‑TTTC‑‑GTAGGCCTGATAAGCGGAGCGCATCAGGCATTTTTGCCTTTGCCATCAGTCt  <  2:1311195/100‑1 (MQ=255)
                                                                                                  |                                                                                                           
TGCTGCCGGATGGCGCTGCACGCAGAGGCACAGCGCCACCGGGTATACAGATTGTTGACAAAGTGCGCTTTGTTTATGCCGGATGCGGCGTTAACGCCTTATCCGGCCTAC‑AAAATCGTGC‑AAATTCAATACATTGC‑‑AGA‑AT‑TTTC‑‑GTAGGCCTGATAAGCGGAGCGCATCAGGCATTTTTGCCTTTGCCATCAGTCT  >  NC_000913/3705886‑3706083

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