Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009133,705,9850TC26.2% 207.9 / 73.1 107intergenic (‑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 (41/38);  new base C (26/2);  total (67/40)
Fisher's exact test for biased strand distribution p-value = 7.35e-05
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 8.63e-01

Only 100 of 109 total aligned reads displayed.

GCTGCCGGATGGCGCTGCACGCAGAGGCACAGCGCCACCGGGTATACAGATTGTTGACAAAGTGCGCTTTGTTTATGCCGGATGCGGCGTTAACGCCTTATCCGGCCTACAAAATCGTGCAAATTCAATACATTGC‑‑A‑GAAT‑T‑‑TTCGTAGGCCTGATAAGCGGAGCGCATCAGGCATTTTTGCCTTTGCCATCAGTCTT  >  NC_000913/3705887‑3706084
                                                                                                  |                                                                                                         
gctgcCGGATGGCGCTGCACGCAGAGGCACAGCGCCACCGGGTATACAGATTGTTGACAAAGTGCGCTTTGTTTATGCCGGATGCGGCGTTAACGCCTTa                                                                                                          <  2:1681662/100‑1 (MQ=255)
gctgcCGGATGGCGCTGCACGCAGAGGCACAGCGCCACCGGGTATACAGATTGTTGACAAAGTGCGCTTTGTTTATGCCGGATGCGGCGTTAACGCCTTa                                                                                                          >  2:297678/1‑100 (MQ=255)
    ccGGATGGCGCTGCACGCAGAGGCACAGCGCCACCGGGTATACAGATTGTTGACAAAGTGCGCTTTGTTTATGCCGGATGCGGCGTTAACGCCTTATCCg                                                                                                      >  1:821922/1‑100 (MQ=255)
     cGGATGGCGCTGCACGCAGAGGCACAGCGCCACCGGGTATACAGATTGTTGACAAAGTGCGCTTTGTTTATGCCGGATGCGGCGTTAACGCCTTATCCgg                                                                                                     >  1:1236090/1‑100 (MQ=255)
      ggATGGCGCTGCACGCAGAGGCACAGCGCCACCGGGTATACAGATTGTTGACAAAGTGCGCTTTGTTTATGCCGGATGCGGCGTTAACGCCTTATCCGGc                                                                                                    <  1:264075/100‑1 (MQ=255)
      ggATGGCGCTGCACGCAGAGGCACAGCGCCACCGGGTATACAGATTGTTGACAAAGTGCGCTTTGTTTATGCCGGATGCGGCGTTAACGCCTTATCCGGc                                                                                                    >  1:1702187/1‑100 (MQ=255)
       gATGGCGCTGCACGCAGAGGCACAGTGCCACCGGGTATACAGATTGTTGACAAAGTGCGCTTTGTTTATGCCGGATGCGGCGTTAACGCCTTATCCGGcc                                                                                                   <  2:677213/100‑1 (MQ=255)
             gcTGCACGCAGAGGCACAGCGCCACCGGGTATACAGATTGTTGACAAAGTGCGCTTTGTTTATGCCGGATGCGGCGTTAACGCCTTATCCGGCCTACaaa                                                                                             >  2:1256370/1‑100 (MQ=255)
             gcTGCACGCAGAGGCACAGCGCCACCGGGTATACAGATTGTTGACAAAGTGCGCTTTGTTTATGCCGGATGCGGCGTTAACGCCTTATCCGGCCTACaaa                                                                                             >  2:1068255/1‑100 (MQ=255)
              cTGCACGCAGAGGCACAGCGCCACCGGGTATACAGATTGTTGACAAAGTGCGCTTTGTTTATGCCGGATAAGGCGTTAACGCCGCATCCGGcataaacaa                                                                                            <  1:1763549/100‑9 (MQ=25)
              cTGCACGCAGAGGCACAGCGCCACCGGGTATACAGATTGTTGACAAAGTGCGCCTGGTTTATGCCGGATGCGGCGTTAACGCCTTATCCGGCCTACaaa                                                                                             >  1:58529/1‑99 (MQ=255)
                 cACGCAGAGGCACAGCGCCACCGGGTATACAGATTGTTGACAAAGTGCGCTTTGTTTATGCCGGATGCGGCGTTAACGCCTTATCCGGCCTACAAAATCg                                                                                         <  1:1100084/100‑1 (MQ=255)
                 cACGCAGAGGCACAGCGCCACCGGGTATACAGATTGTTGACAAAGTGCGCTTTGTTTATGCCGGATGCGGCGTTAACGCCTTATCCGGCCTACAAAATCg                                                                                         <  2:1149932/100‑1 (MQ=255)
                 cACGCAGAGGCACAGCGCCACCGGGTATACAGATTGTTGACAAAGTGCGCTTTGTTTATGCCGAATCCGGCGTTAACACCTTATCCGGCCTACAAAATCg                                                                                         >  2:592340/1‑100 (MQ=255)
                   cGCAGAGGCACAGCGCCACCGGGTATACAGATTGTTGACAAAGTGCGCTTTGTTTATGCCGGATGCGGCGTTAACGCCTTATCCGGCCTACAAAATCGTg                                                                                       <  2:734117/100‑1 (MQ=255)
                    gCAGAGGCACAGCGCCACCGGGTATACAGATTGTTGACAAAGTGCGCTTTGTTTATGCCGGATGCGGCGTTAACGCCTTATCCGGCCTACAAAATCGTGc                                                                                      <  2:706941/100‑1 (MQ=255)
                        agGCACAGCGCCACCGGGTATACAGATTGTTGACAAAGTGCGCTTTGTTTATGCCGGATGCGGCGTTAACGCCTTATCCGGCCTACAAAATCGTGCAAAt                                                                                  <  1:1606729/100‑1 (MQ=255)
                         ggCACAGCGCCACCGGGTATACAGATTGTTGACAAAGTGCGCTTTGTTTATGCCGGATGCGGCGTTAACGCCTTATCCGGCCTACAAAATCGTGCAAAtt                                                                                 >  2:1523415/1‑100 (MQ=255)
                         ggCACAGCGCCACCGGGTATACAGATTGTTGACAAAGTGCGCTTTGTTTATGCCGGATGCGGCGTTAACGCCTTATCCGGCCTACAAAATCGTGCAAAtt                                                                                 >  2:1532130/1‑100 (MQ=255)
                           cacaGCGCCACCGGGTATACAGATTGTTGACAAAGTGCGCTTTGTTTATGCCGGATGCGGCGTTAACGCCTTATCCGGCCTACAAAATCGTGCAAATTCa                                                                               >  2:488539/1‑100 (MQ=255)
                           cacaGCGCCACCGGGTATACAGATTGTTGACAAAGTGCGCTTTGTTTATGCCGGATGCGGCGTTAACGCCTTATCCGGCCTACAAAATCGTGCAAATTCa                                                                               >  2:1680175/1‑100 (MQ=255)
                           cacaGCGCCACCGGGTATACAGATTGTTGACAAAGTGCGCTTTGTTTATGCCGGATGCGGCGTTAACGCCTTATCCGGCCTACAAAATCGTGCAAATTCa                                                                               >  2:1309767/1‑100 (MQ=255)
                            acaGCGCCACCGGGTATACAGATTGTTGACAAAGTGCGCTTTGTTTATGCCGGATGCGGCGTTAACGCCTTATCCGGCCTACAAAATCGTGCAAATTCa                                                                               <  2:1476944/99‑1 (MQ=255)
                             caGCGCCACCGGGTATACAGATTGTTGACAAAGTGCGCTTTGTTTATGCCGGATGCGGCGTTAACGCCTTATCCGGCCTACAAAATCGTGCAAATTCAAt                                                                             <  2:338379/100‑1 (MQ=255)
                               gcgcCACCGGGTATACAGATTGTTGACAAAGTGCGCTTTGTTTATGCCGGATGCGGCGTTAACGCCTTATCCGGCCTACAAAATCGTGCAAATTCAATAc                                                                           <  1:1032586/100‑1 (MQ=255)
                               gcgcCACCGGGTATACAGATTGTTGACAAAGTGCGCTTTGTTTATGCCGGATGCGGCGTTAACGCCTTATCCGGCCTACAAAATCGTGCAAATTCAATAc                                                                           <  1:1033347/100‑1 (MQ=255)
                                cgcCACCGGGTATACAGATTGTTGACAAAGTGCGCTTTGTTTATGCCGGATGCGGCGTTAACGCCTTATCCGGCCTACAAAATCGTGCAAATTCAATAc                                                                           <  2:1047255/99‑1 (MQ=255)
                                cgcCACCGGGTATACAGATTGTTGACAAAGTGCGCTTTGTTTATGCCGGATGCGGCGTTAACGCCTTATCCGGCCTACAAAATCGTGCAAATTCAATACa                                                                          >  2:911581/1‑100 (MQ=255)
                                    aCCGGGGATAGAGATTGTTGACAAAGTGCGCTTTGTTTATGCCGGATGCGGCGTTAACGCCTTATCCGGCCTACAAAATCGTGCAAATTCAATACATTGc                                                                      <  1:143326/100‑1 (MQ=255)
                                      cGGGTATACAGATTGTTGACAAAGTGCGCTTTGTTTATGCCGGATGCGGCGTTAACGCCTTATCCGGCCTACAAAATCGTGCAAATTCAATACATTGC‑‑A‑g                                                                 >  1:1762033/1‑100 (MQ=255)
                                      cGGGTATACAGATTGTTGACAAAGTGCGCTTTGTTTATGCCGGATGCGGCGTTAACGCCTTATCCGGCCTACAAAATCGTGCAAATTCAATACATTGC‑‑A‑g                                                                 <  2:541449/100‑1 (MQ=255)
                                       gggTATACAGATTGTTGACAAAGTGCGCTTTGTTTATGCCGGATGCGGCGTTAACGCCTTATCCGGCCTACAAAATCGTGCAAATTCAATACATTGC‑‑A‑Ga                                                                >  1:503843/1‑100 (MQ=255)
                                        ggTATACAGATTGTTGACAAAGTGCGCTTTGTTTATGCCGGATGCGGCGTTAACGCCTTATCCGGCCTACAAAATCGTGCAAATTCAATACATTGC‑‑A‑Gaa                                                               >  2:1191166/1‑100 (MQ=255)
                                         gTATACAGATTGTTGACAAAGTGCGCTTTGTTTATGCCGGATGCGGCGTTAACGCCTTATCCGGCCTACAAAATCGTGCAAATTCAATACATTGC‑‑A‑GAAt                                                              <  2:170582/100‑1 (MQ=255)
                                          tataCAGATTGTTGACAAAGTGCGCTTTGTTTATGCCGGATGCGGCGTTAACGCCTTATCGGGCCTACAAAATCGTGCAAATTCAATACATTGC‑‑A‑GAAt‑t                                                            <  1:620334/100‑1 (MQ=255)
                                          tataCAGATTGTTGACAAAGTGCGCTTTGTTTATGCCGGATGCGGCGTTAACGCCTTATCCGGCCTACAAAATCGTGCAAATTCAATACATTGC‑‑A‑GAAt‑t                                                            <  1:595451/100‑1 (MQ=255)
                                            taCAGATTGTTGACAAAGTGCGCTTTGTTTATGCCGGATGCGGCGTTAACGCCTTATCCGGCCTACAAAATCGTGCAAATTCAATACATTGC‑‑A‑GAAt‑t‑‑tt                                                        >  1:234890/1‑100 (MQ=255)
                                              cAGATTGTTGACAAAGTGCGCTTTGTTTATGCCGGATGCGGCGTTAACGCCTTATCCGGCCTACAAAATCGTGCAAATTCAATACATTGC‑‑A‑GAAT‑T‑‑TTCg                                                      <  1:1744610/100‑1 (MQ=255)
                                                gATTGTTGACAAAGTGCGCTTTGTTTATGCCGGATGCGGCGTTAACGCCTTATCCGGCCTACAAAATCGTGCAAATTCAATACATTGC‑‑A‑GAAT‑T‑‑TTCGTa                                                    >  1:1403868/1‑100 (MQ=255)
                                                   tgttgACAAAGTGCGCTTTGTTTATGCCGGATGCGGCGTTAACGCCTTATCCGGCCTACAAAATCGTGCAAATTCAATACATTGC‑‑A‑GAAT‑T‑‑TTCGTAGGc                                                 <  2:1740260/100‑1 (MQ=255)
                                                    gttgACAAAGTGCGCTTTGTTTATGCCGGATGCGGCGTTAACGCCTTATCCGGCCTACAAAATCGTGCAAATTCAATACATTGC‑‑A‑GAAT‑T‑‑TTCGTAGGcc                                                >  1:1169348/1‑100 (MQ=255)
                                                    gttgACAAAGTGCGCTTTGTTTATGCCGGATGCGGCGTTAACGCCTTATCCGGCCTACAAAATCGTGCAAATTCAATACATTGC‑‑A‑GAAT‑T‑‑TTCGTAGGcc                                                >  1:103363/1‑100 (MQ=255)
                                                       gACAAAGTGCGCTTTGTTTATGCCGGATGCGGCGTTAACGCCTTATCCGGCCTACAAAATCGTGCAAATTCAATACATTGC‑‑A‑GAAT‑T‑‑TTCGTAGGCCTGa                                             >  2:1340399/1‑100 (MQ=255)
                                                       gACAAAGTGCGCTTTGTTTATGCCGGATGCGGCGTTAACGCCTTATCCGGCCTACAAAATCGTGCAAATTCAATACATTGC‑‑A‑GAAT‑T‑‑TTCGTAGGCCTGa                                             <  2:1702200/100‑1 (MQ=255)
                                                         cAAAGTGCGCTTTGTTTATGCCGGATGCGGCGTTAACGCCTTATCCGGCCTACAAAATCGTGCAAATTCAATACATTGC‑‑A‑GAAT‑T‑‑TTCGTAGGCCTGATa                                           >  2:901313/1‑100 (MQ=255)
                                                          aaaGTGCTCTTTGTTTATGCCGGATGCGGCGTTAAAGGCTTATCCGGCCTACAAAATCGTGCAAATTCAATACATTGC‑‑A‑GAAT‑T‑‑TTCGTAGGCCTGATa                                           >  1:1055188/1‑99 (MQ=255)
                                                          aaaGTGCGCTTTGTTTATGCCGGATGCGGCGTTAACGCCTTATCCGGCCTACAAAATCGTGCAAATTCAATACATTGC‑‑A‑GAAT‑T‑‑TTCGTAGGCCTGATaa                                          >  1:1055189/1‑100 (MQ=255)
                                                            aGTGCGCTTTGTTTATGCCGGATGCGGCGTTAACGCCTTATCCGGCCTACAAAATCGTGCAAATTCAATACATTGC‑‑A‑GAAT‑T‑‑TTCGTAGGCCTGATAAGc                                        <  1:520659/100‑1 (MQ=255)
                                                                cgcTTTGTTTATGCCGGATGCGGCGTTAACGCCTTATCCGGCCTACAAAATCGTGCAAATTCAATACATTGC‑‑A‑GAAT‑T‑‑TTCGTAGGCCTGATAAGCGGAg                                    <  1:679691/100‑1 (MQ=255)
                                                                cgcTTTGTTTATGCCGGATGCGGCGTTAACGCCTTATCCGGCCTACAAAATCGTGCAAATTCAATACATTGC‑‑A‑GAAT‑T‑‑TTCGTAGGCCTGATAAGCGGAg                                    >  1:1010374/1‑100 (MQ=255)
                                                                cgcTTTGTTTATGCCGGATGCGGCGTTAACGCCTTATCCGGCCTACAAAATCGTGCAAATTCAATACATTGC‑‑A‑GAAT‑T‑‑TTCGTAGGCCTGATAAGCGGAg                                    <  1:1536712/100‑1 (MQ=255)
                                                                 gcTTTGTTTATGCCGGATGCGGCGTTAACGCCTTATCCGGCCTACAAAATCGTGCAAATTCAATACATTGC‑‑A‑GAAT‑T‑‑TTCGTAGGCCTGATAAGCGGAgc                                   <  1:363970/100‑1 (MQ=255)
                                                                  tgatttgtaggCCGGATAAGGCGTTCACGCCGCATCCGGCCTAC‑AAATCATGCAAATCCAAAATATTGC‑‑A‑GGAG‑T‑‑TATGTAGGCCTGATAAGCGTAgcgc                                 >  1:551126/11‑100 (MQ=18)
                                                                  cTTTGTTTATGCCGGATGCGGCGTTAACGCCTTATCCGGCCTACAAAATCGTGCAAATTCAATACATTGC‑‑A‑GAAT‑T‑‑TTCGTAGGCCTGATAAGCGGAgcg                                  >  2:1438777/1‑100 (MQ=255)
                                                                  cTTTGTTTATGCCGGATGCGGCGTTAACGCCTTATCCGGCCTACAAAATCGTGCAAATTCAATACATTGC‑‑A‑GAAT‑T‑‑TTCGTAGGCCTGATAAGCGGAgcg                                  >  2:1548186/1‑100 (MQ=255)
                                                                    ttGTTTATGCCGGATGCGGCGTTAACGCCTTATCCGGCCTACAAAATCGTGCAAATTCAATACATTGC‑‑A‑GAAT‑T‑‑TTCGTAGGCCTGATAAGCGGAGCGCa                                >  1:10159/1‑100 (MQ=255)
                                                                    ttGTTTATGCCGGATGCGGCGTTAACGCCTTATCCGGCCTACAAAATCGTGCAAATTCAATACATTGC‑‑A‑GAAT‑T‑‑TTCGTAGGCCTGATAAGCGGAGCGCa                                <  2:1631715/100‑1 (MQ=255)
                                                                    ttGTTTATGCCGGATGCGGCGTTAACGCCTTATCCGGCCTACAAAATCGTGCAAATTCAATACATTGC‑‑A‑GAAT‑T‑‑TTCGTAGGCCTGATAAGCGGAGCGCa                                <  2:1037093/100‑1 (MQ=255)
                                                                       tttATGCCGGATGCGGCGTTAACGCCTTATCCGGCCTACAAAATCGTGCAAATTCAATACATTGC‑‑A‑GAAT‑T‑‑TTCGTAGGCCTGATAAGCGGAGCGCATCa                             <  2:1144574/100‑1 (MQ=255)
                                                                       tttATGCCGGATGCGGCGTTAACGCCTTATCCGGCCTACAAAATCGTGCAAATTCAATACATTGC‑‑A‑GAAT‑T‑‑TTCGTAGGCCTGATAAGCGGAGCGCATCa                             <  1:874024/100‑1 (MQ=255)
                                                                       tttATGCCGGATGCGGCGTTAACGCCTTATCCGGCCTACAAAATCGTGCAAATTCAATACATTGC‑‑A‑GAAT‑T‑‑TTCGTAGGCCTGATAAGCGGAGCGCATCa                             <  1:59568/100‑1 (MQ=255)
                                                                       tttATGCCGGATGCGGCGTTAACGCCTTATCCGGCCTACAAAATCGTGCAAATTCAATACATTGC‑‑A‑GAAT‑T‑‑TTCGTAGGCCTGATAAGCGGAGCGCATCa                             <  1:59539/100‑1 (MQ=255)
                                                                        gtaggCCGGATAAGGCGTTCACGCCGCATCCGGCCTACAAAAACATGCAAATTCAATAGATTGC‑‑A‑GAGA‑T‑‑TATGTAGGCCTGATAAGCGTAGCGCATCAg                            >  1:311414/5‑100 (MQ=17)
                                                                         taggCCGGATAAGGCGTTAACGCCGCATCCGGCCTACAAAATCGTGCAAATTCAATACATTGC‑‑A‑GAAT‑T‑‑TTCGTAGGCCTGATAAGCGGAGCGCATCAgg                           >  1:424834/4‑100 (MQ=18)
                                                                          aggCCGGATAAGGCGTTTACGCCGCATCCGGCCTACAAAAGCGTGCAAATTCAATACATTGC‑‑ATGGGC‑C‑‑AT‑GTAGGCCTGATAAGCGTAGCGCATCAGGc                          >  1:576835/3‑100 (MQ=17)
                                                                          aggCCGGATAAGGCGTTTACGCCGCATCCGGCCTACAAAAGCATGCAAATTCAATATATTGC‑‑A‑GAGT‑T‑‑CACGTAGGCCTGATAAGCGTAGCGCATCAGGc                          >  1:1304631/3‑100 (MQ=17)
                                                                          aggCCGGATAAGGCGTTTACGCCGCATCCGGCCTACAAAAGCATGCAAATTCAATATATTGC‑‑A‑GAGT‑T‑‑CACGTAGGCCTGATAAGCGTAGCGCATCAAGc                          >  2:1127876/3‑100 (MQ=14)
                                                                          aggCCGGATAAGGCGTTTACGCCGCATCCGGCCTACAAAACAGTGCAAATTCAATATATTGC‑‑A‑ATCT‑T‑‑CTGGTAGGCCTGATAAGCATCGCGCATCAGGc                          >  2:1387766/3‑100 (MQ=18)
                                                                          aggCCGGATAAGGCGTTCACGCCGCATCCGGCCTACAAAATTGTGTAAATTCAATAAATTGA‑‑T‑CGAT‑C‑‑CACGTAGGCCTGATAAGCGTAGCGCATCAgg                           >  2:80290/3‑99 (MQ=18)
                                                                          aggCCGGATAAGGCGTTCACGCCGCATCCGGCCTACAAAAACCTGCAAATTCAATATATTGC‑‑A‑GGAT‑A‑‑CGTGTAGGCCTGATAAGCGTAGCGCATCAGGc                          >  2:863322/3‑100 (MQ=17)
                                                                          aggCCGGATAAGGCGTTCACGCCGCATCCGGCCTAC‑AAATCATGCAAATCCAAAATATTGC‑‑A‑GGAG‑T‑‑TATGTAGGCCTGATAAGCGTAGCGCATCAGGCa                         >  1:1016666/3‑100 (MQ=21)
                                                                          aggCCGGATAAGGCGTTCACGCCGCATCCGGCCTAC‑AAATCATGCAAATCCAAAATATTGC‑‑A‑GGAG‑T‑‑TATGTAGGCCTGATAAGCGTAGCGCATCAGGCa                         >  1:1634581/3‑100 (MQ=21)
                                                                          aggCCGGATAAGGCGTTCACGCCGCATCCGGCCTAC‑AAATCATGCAAATCCAAAATATTGC‑‑A‑GGAG‑T‑‑TATGTAGGCCTGATAAGCGTAGCGCATCAGGCa                         >  1:371322/3‑100 (MQ=21)
                                                                           tGCCGGATGCGGCGTTAACGCCTTATCCGGCCTACAAAATCGTGCAAATTCAATACATTGC‑‑A‑GAAT‑T‑‑TTCGTAGGCCTGATAAGCGGAGCGCATCAGGCa                         >  1:1423901/1‑100 (MQ=255)
                                                                           tGCCGGATGCGGCGTTAACGCCTTATCCGGCCTACAAAATCGTGCAAATTCAATACATTGC‑‑A‑GAAT‑T‑‑TTCGTAGGCCTGATAAGCGGAGCGCATCAGGCa                         >  2:353382/1‑100 (MQ=255)
                                                                           ggCCGGATAAGGCGTTTACGCCGCATCCGGCCTACAAAAATGTGCAAATTCAATAAATTGCA‑A‑TTCA‑A‑‑CTTGTAGGCCTGATAAGCGCAGCGCATCAGGc                          >  2:1356032/2‑100 (MQ=18)
                                                                           ggCCGGATAAGGCGTTTACGCCGCATCCGGCCTACAAAAATGTGCAAATTCAATAAATTGC‑‑A‑‑ATT‑CAACTAATAGGCCTGATAAGCGCAGCGCATCAGGc                          >  2:1382502/2‑100 (MQ=18)
                                                                           ggCCGGATAAGGCGTTCACGCCGCATCCGGCCTACAAAATCATGAAAATTCAATACATTGC‑‑A‑AGAT‑T‑‑TTCGTAGGCCTGATAAGAGTAGCGCATCAGGCa                         >  2:1792775/2‑100 (MQ=17)
                                                                           ggCCGGATAAGGCGTTCACGCCGCATCCGGCCTACAAAACTTTGCAAATTCAATATATTGC‑‑A‑ATCT‑C‑‑CGTGTAGGCCTGATAAGCGTAGCGCATCAGGCa                         >  2:972699/2‑100 (MQ=18)
                                                                            gtcggataaggcgTTCACGCCGCATCCGACCCACAAAGTCTTGCAAAATCAATAAATTGC‑‑A‑GGAA‑C‑‑GATGTAGGCCTGATAAGCGTAGCGCATCAGGCaa                        >  2:1383905/10‑99 (MQ=21)
                                                                            gCCGGATAAGGCGTTTACGCCGCATCCGGCCTACAAAATCGTGCTAATTCAATATATTGC‑‑A‑GAAA‑C‑‑CTTGTAGGACTGATAAGAGTAGCGCATAAggaa                         >  2:785541/1‑97 (MQ=17)
                                                                            gCCGGATAAGGCGTTCACGCCGCATCCGGCCTACGAAATCGTGCTAATTCAATATATTGT‑‑G‑GGCT‑C‑‑CACGTAGGCCTGATAAGCGTAGCGCATCAGGCag                        >  2:1642074/1‑99 (MQ=18)
                                                                            gCCGGATAAGGCGTTCACGCCGCATCCGGCCTACAAAGTCGTATAAATTCAATGTATTGC‑‑A‑GTTT‑A‑‑TTGGTAGGCCTGATAAGCGCAGCGCATCAGGCaa                        >  1:466815/1‑99 (MQ=21)
                                                                            gCCGGATAAGGCGTTCACGCCGCATCCGGCCTAC‑AAACCGCGCAAATTCAATATATTGCGGA‑GAAA‑A‑‑T‑‑GTAGGCCTGATAAGCGTAGCGCATCAGGctgt                       >  2:1815709/1‑97 (MQ=21)
                                                                             ccggataaggcgTTTACGCCGCATCCGGCCTACAAAAGCGTGCAAATTCAATACATTGC‑‑ATGGGC‑C‑‑AT‑GTAGGCCTGATAAGCGTAGCGCATCAGGCaat                       >  1:1461589/9‑98 (MQ=17)
                                                                              cGGATGCGGCGTTAACGCCTTATCCGGCCTACAAAATCGTGCAAATTCAATACATTGC‑‑A‑GAAT‑T‑‑TTCGTAGGCCTGATAAGCGGAGCGCATCAGGCAttt                      >  2:418232/1‑100 (MQ=255)
                                                                              cGGATGCGGCGTTAACGCCTTATCCGGCCTACAAAATCGTGCAAATTCAATACATTGC‑‑A‑GAAT‑T‑‑TTCGTAGGCCTGATAAGCGGAGCGCATCAGGCAtt                       >  2:729353/1‑99 (MQ=255)
                                                                               ggataaggcgTTTACGCCGCATCCGGCCTACAAAAGCATGCAAATTCAATATATTGC‑‑A‑GAGT‑T‑‑CACGTAGGCCTGATAAGCGTAGCGCATCAGGCAgttt                     >  1:1562793/7‑100 (MQ=17)
                                                                               ggataaggcgTTTACGCCGCAACCGGCCTACAAATTC‑TGCAAATTCAACAAATTGC‑‑‑‑GAATCC‑‑CTTGTAGGCCTGATAAGCGTAGCGCATCAGGCAtttg                     >  1:1809809/7‑98 (MQ=18)
                                                                               ggataaggcgTTCACGCCTTATCCGGCCTACAAAATCGCTTAAATTCAATATATTGC‑‑A‑CAAT‑A‑‑CTCGTAGGCCTAATAAGCGTAGAGCATCAGGCaat                       >  2:228972/7‑96 (MQ=17)
                                                                               ggataaggcgTTCACGCCGCATCCGGCCTACAAAATCATGAAAATTCAATACATTGC‑‑A‑AGAT‑T‑‑TTCGTAGGCCTGATAAGCGTAGCGCATCAGGCacgct                     >  2:1450541/7‑96 (MQ=17)
                                                                               ggataaggcgTTCACGCCGCATCCGGCCTACAAAACTGTGCAAATTCAATATATTGC‑‑A‑GGAA‑A‑‑CACGTAGGCCTGATAAGCGAAGCGCATCAGGCAgttt                     >  2:128892/7‑100 (MQ=18)
                                                                                     ggcgTTAACGCCTTATCCGGCCTACAAAATCGTGCAAATTCAATACATTGC‑‑A‑GAAT‑T‑‑TTCGTAGGCCTGATAAGCGGAGCGCATCAGGCATTtttgcctt               >  2:142002/1‑100 (MQ=255)
                                                                                       cgTTAACGCCTTATCCGGCCTACAAAATCGTGCAAATTCAATACATTGC‑‑A‑GAAT‑T‑‑TTCGTAGGCCTGATAAGCGGAGCGCATCAGGCATTtttgcctttg             <  1:1523404/100‑1 (MQ=255)
                                                                                       cgTTAACGCCTTATCCGGCCTACAAAATCGTGCAAATTCAATACATTGC‑‑A‑GAAT‑T‑‑TTCGTAGGCCTGATAAGCGGAGCGCATCAGGCATTtttgcctttg             <  1:138339/100‑1 (MQ=255)
                                                                                          tAACGCCTTATCCGACCTACAAAATCGTGCAAATTCAATACATTGC‑‑A‑GAAT‑T‑‑TTCGTAGGCCTGATAAGCGGAGCGCATCAGGCATTTTTGCCTTTGCCa          >  1:1143213/1‑100 (MQ=255)
                                                                                            aCGCCTTATCCGGCCTACAAAATCGTGCAAATTCAATACATTGC‑‑A‑GAAT‑T‑‑TTCGTAGGCCTGATAAGCGGAGCGCATCAGGCATTTTTGCCTTTGCCATc        <  2:1449394/100‑1 (MQ=255)
                                                                                             cGCCTTATCCGGCCTACAAAATCGTGCAAATTCAATACATTGC‑‑A‑GAAT‑T‑‑TTCGTAGGCCTGATAAGCGGAGCGCATCAGGCATTTTTGCCTTTGCCATCa       <  2:37410/100‑1 (MQ=255)
                                                                                                  tATCCGGCCTACAAAATCGTGCAAATTCAATACATTGC‑‑A‑GAAT‑T‑‑TTCGTAGGCCTGATAAGCGGAGCGCATCAGGCATTTTTGCCTTTGCCATCAGTCtt  >  2:653818/1‑100 (MQ=255)
                                                                                                  tATCCGGCCTACAAAATCGTGCAAATTCAATACATTGC‑‑A‑GAAT‑T‑‑TTCGTAGGCCTGATAAGCGGAGCGCATCAGGCATTTTTGCCTTTGCCATCAGTAtt  >  2:930055/1‑100 (MQ=255)
                                                                                                  |                                                                                                         
GCTGCCGGATGGCGCTGCACGCAGAGGCACAGCGCCACCGGGTATACAGATTGTTGACAAAGTGCGCTTTGTTTATGCCGGATGCGGCGTTAACGCCTTATCCGGCCTACAAAATCGTGCAAATTCAATACATTGC‑‑A‑GAAT‑T‑‑TTCGTAGGCCTGATAAGCGGAGCGCATCAGGCATTTTTGCCTTTGCCATCAGTCTT  >  NC_000913/3705887‑3706084

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