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‑‑AGAATT‑‑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‑‑Ag                                                                 >  1:1762033/1‑100 (MQ=255)
                                      cGGGTATACAGATTGTTGACAAAGTGCGCTTTGTTTATGCCGGATGCGGCGTTAACGCCTTATCCGGCCTACAAAATCGTGCAAATTCAATACATTGC‑‑Ag                                                                 <  2:541449/100‑1 (MQ=255)
                                       gggTATACAGATTGTTGACAAAGTGCGCTTTGTTTATGCCGGATGCGGCGTTAACGCCTTATCCGGCCTACAAAATCGTGCAAATTCAATACATTGC‑‑AGa                                                                >  1:503843/1‑100 (MQ=255)
                                        ggTATACAGATTGTTGACAAAGTGCGCTTTGTTTATGCCGGATGCGGCGTTAACGCCTTATCCGGCCTACAAAATCGTGCAAATTCAATACATTGC‑‑AGaa                                                               >  2:1191166/1‑100 (MQ=255)
                                         gTATACAGATTGTTGACAAAGTGCGCTTTGTTTATGCCGGATGCGGCGTTAACGCCTTATCCGGCCTACAAAATCGTGCAAATTCAATACATTGC‑‑AGAAt                                                              <  2:170582/100‑1 (MQ=255)
                                          tataCAGATTGTTGACAAAGTGCGCTTTGTTTATGCCGGATGCGGCGTTAACGCCTTATCGGGCCTACAAAATCGTGCAAATTCAATACATTGC‑‑AGAAtt                                                            <  1:620334/100‑1 (MQ=255)
                                          tataCAGATTGTTGACAAAGTGCGCTTTGTTTATGCCGGATGCGGCGTTAACGCCTTATCCGGCCTACAAAATCGTGCAAATTCAATACATTGC‑‑AGAAtt                                                            <  1:595451/100‑1 (MQ=255)
                                            taCAGATTGTTGACAAAGTGCGCTTTGTTTATGCCGGATGCGGCGTTAACGCCTTATCCGGCCTACAAAATCGTGCAAATTCAATACATTGC‑‑AGAAtt‑‑tt                                                        >  1:234890/1‑100 (MQ=255)
                                              cAGATTGTTGACAAAGTGCGCTTTGTTTATGCCGGATGCGGCGTTAACGCCTTATCCGGCCTACAAAATCGTGCAAATTCAATACATTGC‑‑AGAATT‑‑TTCg                                                      <  1:1744610/100‑1 (MQ=255)
                                                gATTGTTGACAAAGTGCGCTTTGTTTATGCCGGATGCGGCGTTAACGCCTTATCCGGCCTACAAAATCGTGCAAATTCAATACATTGC‑‑AGAATT‑‑TTCGTa                                                    >  1:1403868/1‑100 (MQ=255)
                                                   tgttgACAAAGTGCGCTTTGTTTATGCCGGATGCGGCGTTAACGCCTTATCCGGCCTACAAAATCGTGCAAATTCAATACATTGC‑‑AGAATT‑‑TTCGTAGGc                                                 <  2:1740260/100‑1 (MQ=255)
                                                    gttgACAAAGTGCGCTTTGTTTATGCCGGATGCGGCGTTAACGCCTTATCCGGCCTACAAAATCGTGCAAATTCAATACATTGC‑‑AGAATT‑‑TTCGTAGGcc                                                >  1:1169348/1‑100 (MQ=255)
                                                    gttgACAAAGTGCGCTTTGTTTATGCCGGATGCGGCGTTAACGCCTTATCCGGCCTACAAAATCGTGCAAATTCAATACATTGC‑‑AGAATT‑‑TTCGTAGGcc                                                >  1:103363/1‑100 (MQ=255)
                                                       gACAAAGTGCGCTTTGTTTATGCCGGATGCGGCGTTAACGCCTTATCCGGCCTACAAAATCGTGCAAATTCAATACATTGC‑‑AGAATT‑‑TTCGTAGGCCTGa                                             >  2:1340399/1‑100 (MQ=255)
                                                       gACAAAGTGCGCTTTGTTTATGCCGGATGCGGCGTTAACGCCTTATCCGGCCTACAAAATCGTGCAAATTCAATACATTGC‑‑AGAATT‑‑TTCGTAGGCCTGa                                             <  2:1702200/100‑1 (MQ=255)
                                                         cAAAGTGCGCTTTGTTTATGCCGGATGCGGCGTTAACGCCTTATCCGGCCTACAAAATCGTGCAAATTCAATACATTGC‑‑AGAATT‑‑TTCGTAGGCCTGATa                                           >  2:901313/1‑100 (MQ=255)
                                                          aaaGTGCTCTTTGTTTATGCCGGATGCGGCGTTAAAGGCTTATCCGGCCTACAAAATCGTGCAAATTCAATACATTGC‑‑AGAATT‑‑TTCGTAGGCCTGATa                                           >  1:1055188/1‑99 (MQ=255)
                                                          aaaGTGCGCTTTGTTTATGCCGGATGCGGCGTTAACGCCTTATCCGGCCTACAAAATCGTGCAAATTCAATACATTGC‑‑AGAATT‑‑TTCGTAGGCCTGATaa                                          >  1:1055189/1‑100 (MQ=255)
                                                            aGTGCGCTTTGTTTATGCCGGATGCGGCGTTAACGCCTTATCCGGCCTACAAAATCGTGCAAATTCAATACATTGC‑‑AGAATT‑‑TTCGTAGGCCTGATAAGc                                        <  1:520659/100‑1 (MQ=255)
                                                                cgcTTTGTTTATGCCGGATGCGGCGTTAACGCCTTATCCGGCCTACAAAATCGTGCAAATTCAATACATTGC‑‑AGAATT‑‑TTCGTAGGCCTGATAAGCGGAg                                    <  1:679691/100‑1 (MQ=255)
                                                                cgcTTTGTTTATGCCGGATGCGGCGTTAACGCCTTATCCGGCCTACAAAATCGTGCAAATTCAATACATTGC‑‑AGAATT‑‑TTCGTAGGCCTGATAAGCGGAg                                    >  1:1010374/1‑100 (MQ=255)
                                                                cgcTTTGTTTATGCCGGATGCGGCGTTAACGCCTTATCCGGCCTACAAAATCGTGCAAATTCAATACATTGC‑‑AGAATT‑‑TTCGTAGGCCTGATAAGCGGAg                                    <  1:1536712/100‑1 (MQ=255)
                                                                 gcTTTGTTTATGCCGGATGCGGCGTTAACGCCTTATCCGGCCTACAAAATCGTGCAAATTCAATACATTGC‑‑AGAATT‑‑TTCGTAGGCCTGATAAGCGGAgc                                   <  1:363970/100‑1 (MQ=255)
                                                                  tgatttgtaggCCGGATAAGGCGTTCACGCCGCATCCGGCCTACAAATCATGCAAATCCAAAATATTGC‑‑AGGAGT‑‑TATGTAGGCCTGATAAGCGTAgcgc                                 >  1:551126/11‑100 (MQ=18)
                                                                  cTTTGTTTATGCCGGATGCGGCGTTAACGCCTTATCCGGCCTACAAAATCGTGCAAATTCAATACATTGC‑‑AGAATT‑‑TTCGTAGGCCTGATAAGCGGAgcg                                  >  2:1438777/1‑100 (MQ=255)
                                                                  cTTTGTTTATGCCGGATGCGGCGTTAACGCCTTATCCGGCCTACAAAATCGTGCAAATTCAATACATTGC‑‑AGAATT‑‑TTCGTAGGCCTGATAAGCGGAgcg                                  >  2:1548186/1‑100 (MQ=255)
                                                                    ttGTTTATGCCGGATGCGGCGTTAACGCCTTATCCGGCCTACAAAATCGTGCAAATTCAATACATTGC‑‑AGAATT‑‑TTCGTAGGCCTGATAAGCGGAGCGCa                                >  1:10159/1‑100 (MQ=255)
                                                                    ttGTTTATGCCGGATGCGGCGTTAACGCCTTATCCGGCCTACAAAATCGTGCAAATTCAATACATTGC‑‑AGAATT‑‑TTCGTAGGCCTGATAAGCGGAGCGCa                                <  2:1631715/100‑1 (MQ=255)
                                                                    ttGTTTATGCCGGATGCGGCGTTAACGCCTTATCCGGCCTACAAAATCGTGCAAATTCAATACATTGC‑‑AGAATT‑‑TTCGTAGGCCTGATAAGCGGAGCGCa                                <  2:1037093/100‑1 (MQ=255)
                                                                       tttATGCCGGATGCGGCGTTAACGCCTTATCCGGCCTACAAAATCGTGCAAATTCAATACATTGC‑‑AGAATT‑‑TTCGTAGGCCTGATAAGCGGAGCGCATCa                             <  2:1144574/100‑1 (MQ=255)
                                                                       tttATGCCGGATGCGGCGTTAACGCCTTATCCGGCCTACAAAATCGTGCAAATTCAATACATTGC‑‑AGAATT‑‑TTCGTAGGCCTGATAAGCGGAGCGCATCa                             <  1:874024/100‑1 (MQ=255)
                                                                       tttATGCCGGATGCGGCGTTAACGCCTTATCCGGCCTACAAAATCGTGCAAATTCAATACATTGC‑‑AGAATT‑‑TTCGTAGGCCTGATAAGCGGAGCGCATCa                             <  1:59568/100‑1 (MQ=255)
                                                                       tttATGCCGGATGCGGCGTTAACGCCTTATCCGGCCTACAAAATCGTGCAAATTCAATACATTGC‑‑AGAATT‑‑TTCGTAGGCCTGATAAGCGGAGCGCATCa                             <  1:59539/100‑1 (MQ=255)
                                                                        gtaggCCGGATAAGGCGTTCACGCCGCATCCGGCCTACAAAAACATGCAAATTCAATAGATTGC‑‑AGAGAT‑‑TATGTAGGCCTGATAAGCGTAGCGCATCAg                            >  1:311414/5‑100 (MQ=17)
                                                                         taggCCGGATAAGGCGTTAACGCCGCATCCGGCCTACAAAATCGTGCAAATTCAATACATTGC‑‑AGAATT‑‑TTCGTAGGCCTGATAAGCGGAGCGCATCAgg                           >  1:424834/4‑100 (MQ=18)
                                                                          aggCCGGATAAGGCGTTTACGCCGCATCCGGCCTACAAAAGCGTGCAAATTCAATACATTGC‑‑ATGGGCC‑‑ATGTAGGCCTGATAAGCGTAGCGCATCAGGc                          >  1:576835/3‑100 (MQ=17)
                                                                          aggCCGGATAAGGCGTTTACGCCGCATCCGGCCTACAAAAGCATGCAAATTCAATATATTGC‑‑AGAGTT‑‑CACGTAGGCCTGATAAGCGTAGCGCATCAGGc                          >  1:1304631/3‑100 (MQ=17)
                                                                          aggCCGGATAAGGCGTTTACGCCGCATCCGGCCTACAAAAGCATGCAAATTCAATATATTGC‑‑AGAGTT‑‑CACGTAGGCCTGATAAGCGTAGCGCATCAAGc                          >  2:1127876/3‑100 (MQ=14)
                                                                          aggCCGGATAAGGCGTTTACGCCGCATCCGGCCTACAAAACAGTGCAAATTCAATATATTGC‑‑AATCTT‑‑CTGGTAGGCCTGATAAGCATCGCGCATCAGGc                          >  2:1387766/3‑100 (MQ=18)
                                                                          aggCCGGATAAGGCGTTCACGCCGCATCCGGCCTACAAAATTGTGTAAATTCAATAAATTGA‑‑TCGATC‑‑CACGTAGGCCTGATAAGCGTAGCGCATCAgg                           >  2:80290/3‑99 (MQ=18)
                                                                          aggCCGGATAAGGCGTTCACGCCGCATCCGGCCTACAAAAACCTGCAAATTCAATATATTGC‑‑AGGATA‑‑CGTGTAGGCCTGATAAGCGTAGCGCATCAGGc                          >  2:863322/3‑100 (MQ=17)
                                                                          aggCCGGATAAGGCGTTCACGCCGCATCCGGCCTACAAATCATGCAAATCCAAAATATTGC‑‑AGGAGT‑‑TATGTAGGCCTGATAAGCGTAGCGCATCAGGCa                         >  1:1016666/3‑100 (MQ=21)
                                                                          aggCCGGATAAGGCGTTCACGCCGCATCCGGCCTACAAATCATGCAAATCCAAAATATTGC‑‑AGGAGT‑‑TATGTAGGCCTGATAAGCGTAGCGCATCAGGCa                         >  1:1634581/3‑100 (MQ=21)
                                                                          aggCCGGATAAGGCGTTCACGCCGCATCCGGCCTACAAATCATGCAAATCCAAAATATTGC‑‑AGGAGT‑‑TATGTAGGCCTGATAAGCGTAGCGCATCAGGCa                         >  1:371322/3‑100 (MQ=21)
                                                                           tGCCGGATGCGGCGTTAACGCCTTATCCGGCCTACAAAATCGTGCAAATTCAATACATTGC‑‑AGAATT‑‑TTCGTAGGCCTGATAAGCGGAGCGCATCAGGCa                         >  1:1423901/1‑100 (MQ=255)
                                                                           tGCCGGATGCGGCGTTAACGCCTTATCCGGCCTACAAAATCGTGCAAATTCAATACATTGC‑‑AGAATT‑‑TTCGTAGGCCTGATAAGCGGAGCGCATCAGGCa                         >  2:353382/1‑100 (MQ=255)
                                                                           ggCCGGATAAGGCGTTTACGCCGCATCCGGCCTACAAAAATGTGCAAATTCAATAAATTGCAATTCAA‑‑CTTGTAGGCCTGATAAGCGCAGCGCATCAGGc                          >  2:1356032/2‑100 (MQ=18)
                                                                           ggCCGGATAAGGCGTTTACGCCGCATCCGGCCTACAAAAATGTGCAAATTCAATAAATTGC‑‑AATTCAACTAATAGGCCTGATAAGCGCAGCGCATCAGGc                          >  2:1382502/2‑100 (MQ=18)
                                                                           ggCCGGATAAGGCGTTCACGCCGCATCCGGCCTACAAAATCATGAAAATTCAATACATTGC‑‑AAGATT‑‑TTCGTAGGCCTGATAAGAGTAGCGCATCAGGCa                         >  2:1792775/2‑100 (MQ=17)
                                                                           ggCCGGATAAGGCGTTCACGCCGCATCCGGCCTACAAAACTTTGCAAATTCAATATATTGC‑‑AATCTC‑‑CGTGTAGGCCTGATAAGCGTAGCGCATCAGGCa                         >  2:972699/2‑100 (MQ=18)
                                                                            gtcggataaggcgTTCACGCCGCATCCGACCCACAAAGTCTTGCAAAATCAATAAATTGC‑‑AGGAAC‑‑GATGTAGGCCTGATAAGCGTAGCGCATCAGGCaa                        >  2:1383905/10‑99 (MQ=21)
                                                                            gCCGGATAAGGCGTTTACGCCGCATCCGGCCTACAAAATCGTGCTAATTCAATATATTGC‑‑AGAAAC‑‑CTTGTAGGACTGATAAGAGTAGCGCATAAggaa                         >  2:785541/1‑97 (MQ=17)
                                                                            gCCGGATAAGGCGTTCACGCCGCATCCGGCCTACGAAATCGTGCTAATTCAATATATTGT‑‑GGGCTC‑‑CACGTAGGCCTGATAAGCGTAGCGCATCAGGCag                        >  2:1642074/1‑99 (MQ=18)
                                                                            gCCGGATAAGGCGTTCACGCCGCATCCGGCCTACAAAGTCGTATAAATTCAATGTATTGC‑‑AGTTTA‑‑TTGGTAGGCCTGATAAGCGCAGCGCATCAGGCaa                        >  1:466815/1‑99 (MQ=21)
                                                                            gCCGGATAAGGCGTTCACGCCGCATCCGGCCTACAAACCGCGCAAATTCAATATATTGCGGAGAAAA‑‑T‑‑GTAGGCCTGATAAGCGTAGCGCATCAGGctgt                       >  2:1815709/1‑97 (MQ=21)
                                                                             ccggataaggcgTTTACGCCGCATCCGGCCTACAAAAGCGTGCAAATTCAATACATTGC‑‑ATGGGCC‑‑ATGTAGGCCTGATAAGCGTAGCGCATCAGGCaat                       >  1:1461589/9‑98 (MQ=17)
                                                                              cGGATGCGGCGTTAACGCCTTATCCGGCCTACAAAATCGTGCAAATTCAATACATTGC‑‑AGAATT‑‑TTCGTAGGCCTGATAAGCGGAGCGCATCAGGCAttt                      >  2:418232/1‑100 (MQ=255)
                                                                              cGGATGCGGCGTTAACGCCTTATCCGGCCTACAAAATCGTGCAAATTCAATACATTGC‑‑AGAATT‑‑TTCGTAGGCCTGATAAGCGGAGCGCATCAGGCAtt                       >  2:729353/1‑99 (MQ=255)
                                                                               ggataaggcgTTTACGCCGCATCCGGCCTACAAAAGCATGCAAATTCAATATATTGC‑‑AGAGTT‑‑CACGTAGGCCTGATAAGCGTAGCGCATCAGGCAgttt                     >  1:1562793/7‑100 (MQ=17)
                                                                               ggataaggcgTTTACGCCGCAACCGGCCTACAAATTCTGCAAATTCAACAAATTGC‑‑GAATCC‑‑CTTGTAGGCCTGATAAGCGTAGCGCATCAGGCAtttg                     >  1:1809809/7‑98 (MQ=18)
                                                                               ggataaggcgTTCACGCCTTATCCGGCCTACAAAATCGCTTAAATTCAATATATTGC‑‑ACAATA‑‑CTCGTAGGCCTAATAAGCGTAGAGCATCAGGCaat                       >  2:228972/7‑96 (MQ=17)
                                                                               ggataaggcgTTCACGCCGCATCCGGCCTACAAAATCATGAAAATTCAATACATTGC‑‑AAGATT‑‑TTCGTAGGCCTGATAAGCGTAGCGCATCAGGCacgct                     >  2:1450541/7‑96 (MQ=17)
                                                                               ggataaggcgTTCACGCCGCATCCGGCCTACAAAACTGTGCAAATTCAATATATTGC‑‑AGGAAA‑‑CACGTAGGCCTGATAAGCGAAGCGCATCAGGCAgttt                     >  2:128892/7‑100 (MQ=18)
                                                                                     ggcgTTAACGCCTTATCCGGCCTACAAAATCGTGCAAATTCAATACATTGC‑‑AGAATT‑‑TTCGTAGGCCTGATAAGCGGAGCGCATCAGGCATTtttgcctt               >  2:142002/1‑100 (MQ=255)
                                                                                       cgTTAACGCCTTATCCGGCCTACAAAATCGTGCAAATTCAATACATTGC‑‑AGAATT‑‑TTCGTAGGCCTGATAAGCGGAGCGCATCAGGCATTtttgcctttg             <  1:1523404/100‑1 (MQ=255)
                                                                                       cgTTAACGCCTTATCCGGCCTACAAAATCGTGCAAATTCAATACATTGC‑‑AGAATT‑‑TTCGTAGGCCTGATAAGCGGAGCGCATCAGGCATTtttgcctttg             <  1:138339/100‑1 (MQ=255)
                                                                                          tAACGCCTTATCCGACCTACAAAATCGTGCAAATTCAATACATTGC‑‑AGAATT‑‑TTCGTAGGCCTGATAAGCGGAGCGCATCAGGCATTTTTGCCTTTGCCa          >  1:1143213/1‑100 (MQ=255)
                                                                                            aCGCCTTATCCGGCCTACAAAATCGTGCAAATTCAATACATTGC‑‑AGAATT‑‑TTCGTAGGCCTGATAAGCGGAGCGCATCAGGCATTTTTGCCTTTGCCATc        <  2:1449394/100‑1 (MQ=255)
                                                                                             cGCCTTATCCGGCCTACAAAATCGTGCAAATTCAATACATTGC‑‑AGAATT‑‑TTCGTAGGCCTGATAAGCGGAGCGCATCAGGCATTTTTGCCTTTGCCATCa       <  2:37410/100‑1 (MQ=255)
                                                                                                  tATCCGGCCTACAAAATCGTGCAAATTCAATACATTGC‑‑AGAATT‑‑TTCGTAGGCCTGATAAGCGGAGCGCATCAGGCATTTTTGCCTTTGCCATCAGTCtt  >  2:653818/1‑100 (MQ=255)
                                                                                                  tATCCGGCCTACAAAATCGTGCAAATTCAATACATTGC‑‑AGAATT‑‑TTCGTAGGCCTGATAAGCGGAGCGCATCAGGCATTTTTGCCTTTGCCATCAGTAtt  >  2:930055/1‑100 (MQ=255)
                                                                                                  |                                                                                                         
GCTGCCGGATGGCGCTGCACGCAGAGGCACAGCGCCACCGGGTATACAGATTGTTGACAAAGTGCGCTTTGTTTATGCCGGATGCGGCGTTAACGCCTTATCCGGCCTACAAAATCGTGCAAATTCAATACATTGC‑‑AGAATT‑‑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: