Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 3,705,985 T→C 32.0% intergenic (‑195/+113) 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,9850TC32.0% 326.7 / 199.3 219intergenic (‑195/+113)dppB/dppAdipeptide/heme ABC transporter permease/dipeptide/heme ABC transporter periplasmic binding protein; dipeptide chemotaxis receptor
Reads supporting (aligned to +/- strand):  ref base T (71/78);  new base C (63/7);  total (134/85)
Fisher's exact test for biased strand distribution p-value = 3.94e-10
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 3.44e-03

Only 100 of 223 total aligned reads displayed.

TGCTGCCGGATGGCGCTGCACGCAGAGGCACAGCGCCACCGGGTATACAGATTGTTGACAAAGTGCGCTTTGTTTATGCCGGATGCGGCGTTAACGCCTTATCCGGCCTAC‑AAAATCGTGCAAATTCAATACATTGC‑‑AGA‑ATTTTC‑‑GTAGGCCTGATAAGCGGAGCGCATCAGGCATTTTTGCCTTTGCCATCAGTCTT  >  NC_000913/3705886‑3706084
                                                                                                   |                                                                                                         
tgctgcCGGATGGCGCTGCACGCAGAGGCACAGCGCCACCGGGTATACAGATTGTTGACAAAGTGCGCTTTGTTTATGCCGGATGCGGCGTTAACGCCtt                                                                                                           >  1:269617/1‑100 (MQ=255)
tgctgcCGGATGGCGCTGCACGCAGAGGCACAGCGCCACCGGGTATACAGATTGTTGACAAAGTGCGCTTTGTTTATGCCGGATGCGGCGTTAACGCCtt                                                                                                           >  1:1364666/1‑100 (MQ=255)
   tgcCGGATGGCGCTGCACGCAGAGGCACAGCGCCACCGGGTATACAGATTGTTGCCAAAGTGCGCTTTGTTTATGCCGGATGCGGCGTTAACGCCTTATc                                                                                                        >  2:792327/1‑100 (MQ=255)
   tgcCGGATGGCGCTGCACGCAGAGGCACAGCGCCACCGGGTATACAGATTGTTGACAAAGTGCGCTTTGTTTATGCCGGATGCGGCGTTAACGCCTTATc                                                                                                        >  2:314627/1‑100 (MQ=255)
    gcCGGATGGCGCTGCACGCAGAGGCACAGCGCCACCGGGTATACAGATTGTTGACAAAGTGCGCTTTGTTTATGCCGGATGCGGCGTTAACGCCTTATcc                                                                                                       >  1:481410/1‑100 (MQ=255)
    gcCGGATGGCGCTGCACGCAGAGGCACAGCGCCACCGGGTATACAGATTGTTGACAAAGTGCGCTTTGTTTATGCCGGATGCGGCGTTAACGCCTTATca                                                                                                       >  2:1258923/1‑99 (MQ=255)
      cGGATGGCGCTGCACGCAGAGGCACAGCGCCACCGGGTATACAGATTGTTGACAAAGTGCGCTTTGTTTATGCCGGATGCGGCGTTAACGCCTTATCCgg                                                                                                     >  2:438106/1‑100 (MQ=255)
      cGGATGGCGCTGCACGCAGAGAAAAAACGCCACCGGGTATACAGATTGTTGACAACATGCGCTTTGTTTATGCAGGATGAGGCGTTAACGCCTTATCCg                                                                                                      >  1:2661478/1‑99 (MQ=21)
       ggATGGCGCTGCACGCAGAGGCACAGCGCCACCGGGTATACAGATTGTTGACAAAGTGCGCTTTGTTTATGCCGGATGCGGCGTTAACGCCTTATCCGGc                                                                                                    <  2:2319675/100‑1 (MQ=255)
        gATGGCGCTGCACGCAGAGGCACAGCGCCACCGGGTATACAGATTGTTGACAAAGTGCGCTTTGTTTATGCCGGATGCGGCGTTAACGCCTTATCCGGcc                                                                                                   <  2:680083/100‑1 (MQ=255)
        gATGGCGCTGCACGCAGAGGCACAGCGCCACCGGGTATACAGATTGTTGACAAAGTGCGCTTTGTTTATGCCGGATAAGGCGTTAACGCCGCATCCGGca                                                                                                   <  1:2060932/100‑2 (MQ=37)
         aTGGCGCTGCACGCAGAGGCACAGCGCCACCGGGTATACAGATTGTTGACAAAGTGCGCTTTGTTTATGCCGGATGCGGCGTTAACGCCTTATCCGGCCt                                                                                                  <  2:1809887/100‑1 (MQ=255)
          tGGCGCTGCACGCAGAGGCACAGCGCCACCGGGTATACAGATTGTTGACAAAGTGCGCTTTGTTTATGCCGGATGCGGCGTTAACGCCTTATCCGGCCTa                                                                                                 >  2:223994/1‑100 (MQ=255)
            gcgcTGCACGCAGAGGCACAGCGCCACCGGGTATACAGATTGTTGACAAAGTGCGCTTTGTTTATGCCGGATGCGGCGTTAACGCCTTATCCGGCCTAC‑a                                                                                              >  1:1311053/1‑100 (MQ=255)
                tGCACGCAGAGGCACAGCGCCACCGGGTATACAGATTGTTGACAAAGTGCGCTTTGTTTATGCCGGATGCGGCGTTAACGCCTTATCCGGCCTAC‑aaagt                                                                                          >  1:1770645/1‑98 (MQ=255)
                tGCACGCACAGGCACAGCGCCACCGGGTATACAGATTGTTGACAAAGTGCGCTTTGTTTATGCCTGATGCGGCGTTAACGCCTTATCCGGCCTAC‑aaa                                                                                            <  1:1911095/98‑1 (MQ=255)
                 gCACGCAGAGGCACAGCGCCACCGGGTATACAGATTGTTGACAAAGTGCGCTTTGTTTATGCCGGATGCGGCGTTAACGCCTTATCCGGCCTAC‑AAAATc                                                                                         >  1:1431514/1‑100 (MQ=255)
                 gCACGCAGAGGCACAGCGCCACCGGGTATACAGATTGTTGACAAAGTGCGCTTTGTTTATGCCGGATGCGGCGTTAACGCCTTATCCGGCCTAC‑AAAATc                                                                                         >  2:1597978/1‑100 (MQ=255)
                  cACGCAGAGGCACAGCGCCACCGGGTATACAGATTGTTGACAAAGTGCGCTTTGTTTATGCCGGATGCGGCGTTAACGCCTTATCCGGCCTAC‑AAAATCg                                                                                        >  2:11003/1‑100 (MQ=255)
                  cACGCAGAGGCACAGCGCCACCGGGTATACAGATTGTTGACAAAGTGCGCTTTGTTTATGCCGGATGCGGCGTTAACGCCTTATCCGGCCTAC‑AAAATCg                                                                                        <  2:1124983/100‑1 (MQ=255)
                  cACGCAGAGGCACAGCGCCACCGGGTATACAGATTGTTGACAAAGTGCGCTTTGTTTATGCCGGATGCGGCGTTAACGCCTTATCCGGCCTAC‑AAAATCg                                                                                        <  2:638664/100‑1 (MQ=255)
                  cACGCAGAGGCACAGCGCCACCGGGTATACAGATTGTTGACAAAGTGCGCTTTGTTTATGCCGGATGCGGCGTTAACGCCTTATCCGGCCTAC‑AAAATCg                                                                                        <  1:2502866/100‑1 (MQ=255)
                    cGCAGAGGCACAGCGCCACCGGGTATACAGATTGTTGACAAAGTGCGCTTTGTTTATGCCGGATGCGGCGTTAACGCCTTATCCGGCCTAC‑AAAATCGTg                                                                                      >  2:2710107/1‑100 (MQ=255)
                    cGCAGAGGCACAGCGCCACCGGGTATACAGATTGTTGACAAAGTGCGCTTTGTTTATGCCGGATGCGGCGTTAACGCCTTATCCGGCCTAC‑AAAATCGTg                                                                                      <  2:1529638/100‑1 (MQ=255)
                     gCAGAGGCACAGCGCCACCGGGTATACAGATTGTTGACAAAGTGCGCTTTGTTTATGCCGGATGCGGCGTTAACGCCTTATCCGGCCTAC‑AAAATCGTGc                                                                                     <  1:1515062/100‑1 (MQ=255)
                       agagGCACAGCGCCACCGGGTATACAGATTGTTGACAAAGTGCGCTTTGTTTATGCCGGATGCGGCGTTAACGCCTTATCCGGCCTAC‑AAAATCGTGCaa                                                                                   <  1:745594/100‑1 (MQ=255)
                          ggCACAGCGCCACCGGGTATACAGATTGTTGACAAAGTGCGCTTTGTTTATGCCGGATGCGGCGTTAACGCCTTATCCGGCCTAC‑AAAATCGTGCAAAtt                                                                                >  2:212578/1‑100 (MQ=255)
                            cacaGCGCCACCGGGTATACAGATTGTTGACAAAGTGCGCTTTGTTTATGCCGGATGCGGCGTTAACGCCTTATCCGGCCTAC‑AAAATCGTGCAAATTCa                                                                              >  1:164047/1‑100 (MQ=255)
                              caGCGCCACCGGGTATACAGATTGTTGACAAAGTGCGCTTTGTTTATGCCGGATGCGGCGTTAACGCCTTATCCGGCCTAC‑AAAATCGTGCAAATTCAAt                                                                            >  2:2791515/1‑100 (MQ=255)
                                  gcCACCGGGTATACAGATTGTTGACAAAGTGCGCTTTGTTTATGCCGGATGCGGCGTTAACGCCTTATCCGGCCTAC‑AAAATCGTGCAAATTCAATAc                                                                          <  1:334449/98‑1 (MQ=255)
                                            ataCAGATTGTTGACAAAGTGCGCTTTGTTTATGCCGGATGCGGCGTTAACGCCTTATCCGGCCTAC‑AAAATCGTGCAAATTCAATACATTGC‑‑AGA‑At                                                             >  1:1897524/1‑98 (MQ=255)
                                             taCAGATTGTTGACAAAGTGCGCTTTGTTTATGCCGGATGCGGCGTTAACGCCTTATCCGGCCTAC‑AAAATCGTGCAAATTCAATACATTGC‑‑AGA‑Atttt                                                          >  1:470312/1‑100 (MQ=255)
                                              aCAGATTGTTGACAAAGTGCGCTTTGTTTATGCCGGATGCGGCGTTAACGCCTTATCCGGCCTAC‑AAAATCGTGCAAATTCAATACATTGC‑‑AGA‑ATTTTc                                                         <  1:2169490/100‑1 (MQ=255)
                                              aCAGATTGTTGACAAAGTGCGCTTTGTTTATGCCGGATGCGGCGTTAACGCCTTATCCGGCCTAC‑AAAATCGTGCAAATTCAATACATTGC‑‑AGA‑ATTTTc                                                         >  1:1246824/1‑100 (MQ=255)
                                               ccgATTGTTGACAAAGTGCGCTTTGTTCATCCCAGAAGCGGCGTTACCGCCTTATCCGGCCTAC‑AAAATCGTGCAAATTCAATACATTGC‑‑AGA‑AATTTC‑‑g                                                      <  1:1740624/98‑1 (MQ=17)
                                               cAGATTGTTGACAAAGTGCGCTTTGTTTATGCCGGATGCGGCGTTAACGCCTTATCCGGCCTAC‑AAAATCGTGCAAATTCAATACATTGC‑‑AGA‑ATTTTC‑‑g                                                      <  2:2439752/100‑1 (MQ=255)
                                                 gATTGTTGACAAAGTGCGCTTTGTTTATGCCGGATGCGGCGTTAACGCCTTATCCGGCCTAC‑AAAATCGTGCAAATTCAATACATTGC‑‑AGA‑ATTTTC‑‑Gt                                                     <  1:223993/99‑1 (MQ=255)
                                                 gATTGTTGACAAAGTGCGCTTTGTTTATGCCGGATGCGGCGTTAACGCCTTATCCGGCCTAC‑AAAATCGTGCAAATTCAATACATTGC‑‑AGA‑ATTTTC‑‑GTa                                                    >  2:814308/1‑100 (MQ=255)
                                                 gATTGTTGACAAAGTGCGCTTTGTTTATGCCGGATGCGGCGTTAACGCCTTATCCGGCCTAC‑AAAATCGTGCAAATTCAATACATTGC‑‑AGA‑ATTTTC‑‑GTa                                                    <  1:1719291/100‑1 (MQ=255)
                                                     gttgACAAAGTGCGCTTTGTTTATGCCGGATGCGGCGTTAACGCCTTATCCGGCCTAC‑AAAATCGTGCAAATTCAATACATTGC‑‑AGA‑ATTTTC‑‑GTAGGcc                                                <  1:105792/100‑1 (MQ=255)
                                                     gttgACAAAGTGCGCTTTGTTTATGCCGGATGCGGCGTTAACGCCTTATCCGGCCTAC‑AAAATCGTGCAAATTCAATACATTGC‑‑AGA‑ATTTTC‑‑GTAGGcc                                                <  1:2887420/100‑1 (MQ=255)
                                                     ataaaCAAAGCGCACTTTGTTTATGCCGGATGCGGCGTTAACGCCTTATCCGGCCTAC‑AAAATCGTGCAAATTCAATACATTGC‑‑AGA‑ATTTTC‑‑GTAGGcc                                                >  1:941227/5‑100 (MQ=18)
                                                         aCAAAGTCCGCTTTGTTTATGCCGGATGCGGCCTTAACGCCTTATCCGGCCTAC‑AAAATCGTTCAAAGTCAATACATTGC‑‑AGA‑ATTTTC‑‑GTAGGCCTGAt                                            <  1:2857788/100‑1 (MQ=18)
                                                           aaaGTGCGCTTTGTTTATGCCGGATGCGGCGTTAACGCCTTATCCGGCCTAC‑AAAATCGTGCAAATTCAATACATTGC‑‑AGA‑ATTTTC‑‑GTAgagatcggaa                                          >  2:1174410/1‑91 (MQ=255)
                                                             aGTGCGCTTTGTTTATGCCGGATGCGGCGTTAACGCCTTATCCGGCCTAC‑AAAATCGTGCAAATTCAATACATTGC‑‑AGA‑ATTTTC‑‑GTAGGCCTGATAAGc                                        >  2:924255/1‑100 (MQ=255)
                                                              gTGCGCTTTGTTTATGCCGGATGCGGCGTTAACGCCTTATCCGGCCTAC‑AAAATCGTGCAAATTCAATACATTGC‑‑AGA‑ATTTTC‑‑GTAGGCCTGATAAGCg                                       <  2:903783/100‑1 (MQ=255)
                                                               tGCGCTTTGTTTATGCCGGATGCGGCGTTAACGCCTTATCCGGCCTAC‑AAAATCGTGCAAATTCAATACATTGC‑‑AGA‑ATTTTC‑‑GTAGGCCTGATAAGCgg                                      <  1:1544549/100‑1 (MQ=255)
                                                               tGCGCTTTGTTTATGACGGATGCGGCGTTAACGCCTTATCCTGCCTAC‑AAAATCGTGCAAATTCTATACATTGC‑‑AGA‑ATTTTC‑‑GTAGGCCTGATAAGCg                                       <  1:835977/99‑1 (MQ=255)
                                                                gcgcTTTGTTTATGCCGGATGCGGCGTTAACGCCTTATCCGGCCTAC‑AAAATCGTGCAAATTCAATACATTGC‑‑AGA‑ATTTTC‑‑GTAGGCCTGATAAGCgg                                      >  2:208679/1‑99 (MQ=255)
                                                                   cTTTGTTTATGCCGGATGCGGCGTTAACGCCTTATCCGGCCTAC‑AAAATCGTGCAAATTCAATACATTGC‑‑AGA‑ATTTTC‑‑GTAGGCCTGATAAGCGGAgcg                                  >  2:868241/1‑100 (MQ=255)
                                                                     ttGTTTATGCCGGATGCGGCGTTAACGCCTTATCCGGCCTAC‑AAAATCGTGCAAATTCAATACATTGC‑‑AGA‑ATTTTC‑‑GTAGGCCTGATAAGCGGAGCGCa                                <  2:2104916/100‑1 (MQ=255)
                                                                       ttgtaggCCGGATAAGGCGTTCACGCCGCATCCGGCCTAC‑AAGTTCGTGCAAATTCAATAAATTGC‑‑A‑A‑TATGAC‑‑GTAGGCCTGATAAGCGTAGCGCATCa                             <  2:2810056/94‑1 (MQ=18)
                                                                       ttgtaggCCGGATAAGGCGTTCACGCCGCATCCGGCCTAC‑AAAATCATGAAAATTCAATACATTGC‑‑AAG‑ATTTTC‑‑GTAGGCCTGATAAGCGTAGCGCATc                              >  2:1850542/7‑100 (MQ=17)
                                                                       ttgtaggCCGGATAAGGCGTTCACGCCGCATCCGGCCTAC‑AAAAACCTGCAAATTCAATATATTGC‑‑AGG‑ATACGT‑‑GTAGGCCTGATAAGCGTAGCGCATc                              >  2:2034384/7‑100 (MQ=18)
                                                                       gTTTATGCCGGATGCGGCGTTAACGCCTTATCCGGCCTAC‑AAAATCGTGCAAATTCAATACATTGC‑‑AGA‑ATTTTC‑‑GTAGGCCTGATAAGCGGAGCGCATc                              <  1:1281196/100‑1 (MQ=255)
                                                                       gTTTATGCCGGATGCGGCGTTAACGCCTTATCCGGCCTAC‑AAAATCGTGCAAATTCAATACATTGC‑‑AGA‑ATTTTC‑‑GTAGGCCTGATAAGCGGAGCGCATc                              >  2:1949783/1‑100 (MQ=255)
                                                                        tgtaggCCGGATAAGGCGTTCACGCCGCATCCGGCATAC‑AAAATTGTGTAAATTCAATAAATTGA‑‑TCG‑ATCCAC‑‑GTAGGCCTGATAAGCGTAGCGCATCa                             >  1:1191559/6‑100 (MQ=18)
                                                                         ttATGCCGGATGCGGCGTTAACGCCTTATCCGGCCTAC‑AAAATCGTGCAAATTCAATACATTGC‑‑AGA‑ATTTTC‑‑GTAGGCCTGATAAGCGGAGCGCATCAg                            >  2:290374/1‑100 (MQ=255)
                                                                         ttATGCCGGATGCGGCGTTAACGCCTTATCCGGCCTAC‑AAAATCGTGCAAATTCAATACATTGC‑‑AGA‑ATTTTC‑‑GTAGGCCTGATAAGCGGAGCGCATCAg                            <  2:689519/100‑1 (MQ=255)
                                                                          taggCCGGATAAGGCGTTCACGCCGCATCCGGCCTAC‑AAAAGCATGAAAATTCAATATATTGC‑‑AGG‑AGCTGC‑‑GTAGGCCTGATAAGCGTAGCGCATCAgg                           >  2:2277299/4‑100 (MQ=17)
                                                                          taggCCGGATAAGGCGTTAACGCCGCATCCGGCCTAC‑AAAATCGTGCAAATTCAATACATTGC‑‑AGA‑ATTTTC‑‑GTAGGCCTGATAAGCGGAGCGCATCAg                            >  1:2308767/4‑99 (MQ=17)
                                                                          tATGCCGGATGCGGCGTTAACGCCTTATCCGGCCTAC‑AAAATCGTGCAAATTCAATACATTGC‑‑AGA‑ATTTTC‑‑GTAGGACTGATAAGCGGAGCGCATCAgg                           >  1:1303748/1‑100 (MQ=255)
                                                                           aggCCGGATAAGGCGTTTACGCCGCATCCGGCCTAC‑AAAACAGTGCAAATTCAATATATTGC‑‑AAT‑CTTGTG‑‑GTAGGCGTGATAAGCATCGCGCATCAGGc                          >  1:1965053/3‑100 (MQ=18)
                                                                           aggCCGGATAAGGCGTTCACGCCGCATCCGGCCTACAAAAATCTTGCCAATTCAATATATTGC‑‑AGG‑ACCAT‑‑‑GTAGGCCTGATAAGCGTAGCGCATCAGGc                          >  2:1385275/3‑100 (MQ=18)
                                                                           aggCCGGATAAGGCGTTCACGCCGCATCCGGCCTAC‑AAAATCGCTTAAATTCAATAGATTGC‑‑ACC‑ATCCTC‑‑GTAGGCCTGATAAGCGTAGCGCATCAGGc                          >  1:1860045/3‑100 (MQ=18)
                                                                           aggCCGGATAAGGCGTTCACGCCGCATCCGGCCTAC‑AAAAGCATGAAAATTCAATATATTGC‑‑AGG‑AGCTGC‑‑GTAGGCCTGATAAGCGTAGCGCATCAGGc                          >  2:2176140/3‑100 (MQ=17)
                                                                           aggCCGGATAAGGCGTTCACGCCGCATCCGGCCTAC‑AAAAACCTGCAAATTCAATATATTGC‑‑AGG‑ATACGT‑‑GTAGGCCTGATAAGCGTAGCGCATCAGGc                          >  1:1079666/3‑100 (MQ=18)
                                                                           aggCCGGATAAGGCGTTCACGCCGCATCCGGCCTAC‑‑AAACCGCGCAAATTCAATATATTGCGGAGA‑AAAT‑‑‑‑GTAGGCCTGATAAGCGTAGCGCATCAGGct                         >  2:2657290/3‑99 (MQ=21)
                                                                           aggCCGGATAAGGCGTTAACGCCGCATCCGGCCTAC‑AAAATCGTGCAAATTCAATACATTGC‑‑AGA‑ATTTTC‑‑GTAGGCCTGATAAGCGGAGCGCATCAGGc                          >  2:1024208/3‑100 (MQ=18)
                                                                           aTGCCGGATGCGGCGTTAACGCCTTATCCGGCCTAC‑AAAATCGTGCAAATTCAATACATTGC‑‑AGA‑ATTTTC‑‑GTAGGCCTGATAAGCGGAGCGCATCAGGc                          >  2:1792401/1‑100 (MQ=255)
                                                                            tGCCGGATGCGGCGTTAACGCCTTATCCGGCCTAC‑AAAATCGTGCAAATTCAATACATTGC‑‑AGA‑ATTTTC‑‑GTAGGCCGGATAAGCGGAGCGCATCAGGCa                         <  1:1687989/100‑1 (MQ=255)
                                                                            ggCCGGATAAGGCGTTCACGCCGCATCCGGCCTACAAAAATCTTGCCAATTCAATATATTGC‑‑AGG‑ACCAT‑‑‑GTAGGCCTGATAAGCGTAGCGCATCAGGCa                         >  2:1554101/2‑100 (MQ=18)
                                                                            ggCCGGATAAGGCGTTCACGCCGCATCCGGCCTAC‑AAAACTTTGCAAATTCAATATATTGC‑‑‑‑A‑ATCTCCGTGTAGGCCTGATAAGCGTAGCGCATCAGGCa                         >  1:2453599/2‑100 (MQ=18)
                                                                            ggCCGGATAAGGCGTTCACGCCGCATCCGGCCTAC‑‑AAATCATGCAAATCCAAAATATTGC‑‑AGG‑AGTTAT‑‑GTAGGCCTGATAAGCGTAGCGCATCAGGCAt                        >  2:2019648/2‑100 (MQ=21)
                                                                             gCCGGATGCGGCGTTAACGCCTTATCCGGCCTAC‑AAAATCGTGCAAATTCAATACATTGC‑‑AGA‑ATTTTC‑‑GTAGGCCTGATAAGCGGAGCGCATCAGGCAt                        >  2:154357/1‑100 (MQ=255)
                                                                             gCCGGATAAGGCGTTCGCGCCGCATCCGGCCTAC‑AACAGTTTGCAAATTCAATAAATTGC‑‑AGA‑ATTCAT‑‑GTAGGCCTGATAAGCGAAGCGCATCAGGCAt                        >  1:2865696/1‑100 (MQ=17)
                                                                             gCCGGATAAGGCGTTCACGCCGCATCCGGCCTACAAAAATCTTGCCAATTCAATATATTGC‑‑AGG‑ACCAT‑‑‑GTAGGCCTGATAAGCGTAGCGCATCAGGCag                        >  2:1312645/1‑99 (MQ=18)
                                                                               cggataaggcgTTCACGCCGCATCCGGCCTAC‑AAAAGCATGAAAATTCAATATATTGC‑‑AGG‑AGCTGC‑‑GTAGGCCTGATAAGCGTAGCGCATCAGGCagtt                      >  1:2731010/8‑97 (MQ=17)
                                                                                ggataaggcgTTTACGCCGCATCCGGCCTAC‑AAAATCGTGCTAATTCAATATATTGC‑‑AGA‑AACCTT‑‑GTCGGCATGATAAGCGTAGCGCCTCAGGCAgttt                     >  2:1875313/7‑100 (MQ=17)
                                                                                ggataaggcgTTTACGCCGCATCCGGCCTAC‑AAAATCGTGCTAATTCAATATATTGC‑‑AGA‑AACCTT‑‑GTAGGCCTGATAAGCGTAGCGCATCAGGCAgttt                     >  1:2284534/7‑100 (MQ=18)
                                                                                ggataaggcgTTTACGCCGCATCCGGCCTAC‑AAAATCGTGCTAATTCAATATATTGC‑‑AGA‑AACCTT‑‑GTAGGCCTGATAAGCGTAGCGCATCAGGCAgttt                     >  1:22064/7‑100 (MQ=18)
                                                                                ggataaggcgTTTACGCCGCATCCGGCCTAC‑AAAAGCATGCAAATTCAATATATTGC‑‑AGA‑GTTCAC‑‑GTAGGCCTGATAAGCGTAGCGCATCAGGCAgtt                      >  2:1618325/7‑99 (MQ=17)
                                                                                ggataaggcgTTCACGCCGCATCCGGCCTACAAAAATCTTGCCAATTCAATATATTGC‑‑AGG‑ACCAT‑‑‑GTAGGCCTGATAAGCGTAGCGCATCAGGCAgttt                     >  2:2382077/7‑100 (MQ=18)
                                                                                ggataaggcgTTCACGCCGCATCCGGCCTAC‑AAAACCTTGCTAATTCAATATATTGC‑‑AGG‑GACTAT‑‑GTAGGCCTGATAAGCATAGCGCATCAGGCagctt                     >  2:1489592/7‑96 (MQ=18)
                                                                                ggataaggcgTTCACGCCGCATCCGGCATAC‑AAAACTGTGCAAATTCAATATATTGC‑‑AGG‑AAACAC‑‑GTAGGCCTGATAAGCGAAGCGCATCAGGCAgttt                     >  1:2247331/7‑100 (MQ=18)
                                                                                 gataaggcgTTCACGCCGCATCCGGCCTAC‑AAAACCTTGCTAATTCAATATATTGC‑‑AGGGACTAT‑‑‑GTAGGCCTGATAAGCATAGCGCATCAGGCagcttt                    >  2:1418324/6‑95 (MQ=21)
                                                                                 gataaggcgTTCACGCCGCATCCGGCCTAC‑AAAAACATGCAAATTCAATAGATTGC‑‑AGA‑GATTAT‑‑GTAGGCCTGATAAGCGTAGCGCATCAGGCAatttt                    >  2:954736/6‑100 (MQ=17)
                                                                                  acgcggcgTTAACGCCTTATCCGGCCTAC‑AAAATCGTGCAAATTCAATACATTGC‑‑AGA‑ATTTTC‑‑GTAGGCCTGATAAGCGGAGCGCATCAGGCATTtttg                   <  2:1864962/98‑1 (MQ=255)
                                                                                     cggcgTTAACGCCTTATCCGGCCTAC‑AAAATCGTGCAAATTCAATACATTGC‑‑AGA‑ATTTTC‑‑GTAGGCCTGATAAGCGGAGCGCATCAGGCATTtttgcct                >  2:1720868/1‑100 (MQ=255)
                                                                                      ggcgTTAACGCCTTATCCGGCCTAC‑AAAATCGTGCAAATTCAATACATTGC‑‑AGA‑ATTTTC‑‑GTAGGCCTGATAAGCGGAGCGCATCAGGCATTtttgcctt               <  1:907132/100‑1 (MQ=255)
                                                                                        cgTTAACGCCTTATCCGGCCTAC‑AAAATCGTGCAAATTCAATACATTGC‑‑AGA‑ATTTTC‑‑GTAGGCCTGATAAGCGGAGCGCATCAGGCATTtttgcctttg             <  2:1107207/100‑1 (MQ=255)
                                                                                        cgTTAACGCCTTATCCGGCCTAC‑AAAATCGTGCAAATTCAATACATTGC‑‑AGA‑ATTTTC‑‑GTAGGCCTGATAAGCGGAGCGCATCAGGCATTtttgcctttg             >  1:2099439/1‑100 (MQ=255)
                                                                                         gTTAACGCCTTATCCGGCCTAC‑AAAATCGTGCAAATTCAATACATTGC‑‑AGA‑ATTTTC‑‑GTAGGCCTGATAAGCGGAGCGCATCAGGCATTtttgcctttgc            >  2:614006/1‑100 (MQ=255)
                                                                                         gTTAACGCCTTATCCGGCCTAC‑AAAATCGTGCAAATTCAATACATTGC‑‑AGA‑ATTTTC‑‑GTAGGCCTGATAAGCGGAGCGCATCAGGCATTtttgcctttgc            <  1:1102315/100‑1 (MQ=255)
                                                                                          ttAACGCCTTATCCGGCCTAC‑AAAATCGTGCAAATTCAATACATTGC‑‑AGA‑ATTTTC‑‑GTAGGCCTGATAAGCGGAGCGCATCAGGCATTtttgcctttgcc           <  1:575528/100‑1 (MQ=255)
                                                                                          ttAACGCCTTATCCGGCCTAC‑AAAATCGTGCAAATTCAATACATTGC‑‑AGA‑ATTTTC‑‑GTAGGCCTGATAAGCGGAGCGCATCAGGCATTtttgcctttgcc           <  1:2082224/100‑1 (MQ=255)
                                                                                            aaCGCCTTATCCGGCCTAC‑AAAATCGTGCAAATTCAATACATTGC‑‑AGA‑ATTTTC‑‑GTAGGCCTGATAAGCGGAGCGCATCAGGCATTTTTGCCTTTGCCAt         >  1:581193/1‑100 (MQ=255)
                                                                                                ccTTATCCGGCCTAC‑AAAATCGTGCAAATTCAATACATTGC‑‑AGA‑ATTTTC‑‑GTAGGCCTGATGAGCGGAGTGCATCAGGCATTTTTGCCTTTGCCATCAg      >  2:379860/1‑99 (MQ=255)
                                                                                                 cTTATCCGGCCTAC‑AAAATCGTGCAAATTCAATACATTGC‑‑AGA‑ATTTTC‑‑GTAGGCCTGATAAGCGGAGCGCATCAGGCATTTTTGCCTTTGCCATCAGTc    >  2:254796/1‑100 (MQ=255)
                                                                                                   tATCCGGCCTAC‑AAAATCGTGCAAATTCAATACATGTC‑‑AGA‑ATTTTC‑‑GTAGGCCTGATAAGCGGAGCGCATCAGGCATTTTTGCCTTTGCCATCAGTCtt  <  1:1574657/100‑1 (MQ=255)
                                                                                                   |                                                                                                         
TGCTGCCGGATGGCGCTGCACGCAGAGGCACAGCGCCACCGGGTATACAGATTGTTGACAAAGTGCGCTTTGTTTATGCCGGATGCGGCGTTAACGCCTTATCCGGCCTAC‑AAAATCGTGCAAATTCAATACATTGC‑‑AGA‑ATTTTC‑‑GTAGGCCTGATAAGCGGAGCGCATCAGGCATTTTTGCCTTTGCCATCAGTCTT  >  NC_000913/3705886‑3706084

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