Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009133,705,9840TG43.2% 254.8 / inf 236intergenic (‑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 (65/69);  new base G (96/6);  total (161/75)
Fisher's exact test for biased strand distribution p-value = 3.33e-15
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 2.15e-01

Only 100 of 239 total aligned reads displayed.

TTGCTGCCGGATGGCGCTGCACGCAGAGGCACAGCGCCACCGGGTATACAGATTGTTGACAAAGTGCGCTTTGTTTATGCCGGATGCGGCGTTAACGCCTTATCCGGCCTACAAAATCGTGCAAATTCAATACATTGC‑‑AGAATTTTC‑‑GTAGGCCTGATAAGCGGAGCGCATCAGGCATTTTTGCCTTTGCCATCAGTCT  >  NC_000913/3705885‑3706083
                                                                                                   |                                                                                                           
  gctgcCGGATGGCGCTGCACGCAGAGGCACAGCGCCACCGGGTATACAGATTGTTGACAAAGTGCGCTTTGTTTATGCCGGATGCGGCGTTAACGCCTTa                                                                                                           >  1:2184057/1‑100 (MQ=255)
   ctgcCGGATGGCGCTGCACGCAGAGGCACAGCGCCACCGGGTATACAGATTGTTGACAAAGTGCGCTTTGTTTATGCCGGATGCGGCGTTAACGCCTTAt                                                                                                          <  2:944258/100‑1 (MQ=255)
     gcCGGATGGCGCTGCACGCAGAGGCACAGCGCCACCGGGTATACAGATTGTTGACAAAGTGCGCTTTGTTTATGCCGGATGCGGCGTTAACGCCTTATcc                                                                                                        >  1:1875217/1‑100 (MQ=255)
      cctgATGGCGCTGCACGCAGAGGCACAGCGCCACCGGGTATACAGATTGTTGACAAAGTGCGCTTTGTTTATGCCGGATGCGGCGTTAACGCCTTATCCg                                                                                                       <  1:57057/97‑1 (MQ=255)
      ccGGATGGCGCTGCACGCAGAGGCACAGCGCCACCGGGTATACAGATTGTTGACAAAGTGCGCTTTGTTTATGCCGGATGCGGCGTTAACGCCTTATCCg                                                                                                       >  1:377671/1‑100 (MQ=255)
        ggATGGCGCTGCACGCAGAGGCACAGCGCCACCGGGTATACAGATTGTTGACAAAGTGCGCTTTGTTTATGCCGGATGCGGCGTTAACGCCTTATCCGGc                                                                                                     <  2:2049852/100‑1 (MQ=255)
         gATGGCGCTGCACGCAGAGGCACAGCGCCACCGGGTATACAGATTGTTGACAAAGTGCGCTTTGTTTATGCCGGATGCGGCGTTAACGCCTTATCCGGcc                                                                                                    <  2:689082/100‑1 (MQ=255)
          aTGGCGCTGCACGCAGAGGCACAGCGCCACCGGGTATACAGATTGTTGACAAAGTGCGCTTTGTTTATGCCGGATGCGGCGTTAACGCCTTATCCGGCCt                                                                                                   <  1:2802245/100‑1 (MQ=255)
           tGGCGCTGCACGCAGAGGCACAGCGCCACCGGGTATACAGATTGTTGACAAAGTGCGCTTTGTTTATGCCGGATGCGGCGTTAACGCCTTATCCGGCCTa                                                                                                  >  2:1322059/1‑100 (MQ=255)
             gcgcTGCACGCAGAGGCACAGCGCCACCGGGTATACAGATTGTTGACAAAGTGCGCTTTGTTTATGCCGGATGCGGCGTTAACGCCTTATCCGGCCTACa                                                                                               <  1:1420361/100‑1 (MQ=255)
              cgcTGCACGCAGAGGCACAGCGCCACCGGGTATACAGATTGTTGACAAAGTGCGCTTTGTTTATGCCGGATGCGGCGTTAACGCCTTATCCGGCCTACaa                                                                                              <  1:2191731/100‑1 (MQ=255)
               gcTGCACGCAGAGGCACAGCGCCACCGGGTATACAGATTGTTGACAAAGTGCGCTTTGTTTATGCCGGATGCGGCGTTAACGCCTTATCCGGCCTACaaa                                                                                             >  1:1808665/1‑100 (MQ=255)
                 tGCACGCAGAGGCACAGCGCCACCGGGTATACAGATTGTTGACAAAGTGCGCTTTGTTTATGCCGGATGCGGCGTTAACGCCTTATCCGGCCTACAAAAt                                                                                           >  1:1007048/1‑100 (MQ=255)
                   cACGCAGAGGCACAGCGCCACCGGGTATACAGATTGTTGACAAAGTGCGCTTTGTTTATGCCGGATGCGGCGTTAACGCCTTATCCGGCCTACAAAATCg                                                                                         <  1:270829/100‑1 (MQ=255)
                   cACGCAGAGGCACAGCGCCACCGGGTATACAGATTGTTGACAAAGTGCGCTTTGTTTATGCCGGATGCGGCGTTAACGCCTTATCCGGCCTACAAAATCg                                                                                         <  2:354577/100‑1 (MQ=255)
                   cACGCAGAGGCACAGCGCCACCGGGTATACAGATTGTTGACAAAGTGCGCTTTGTTTATGCCGGATGCGGCGTTAACGCCTTATCCGGCCTACAAAATCg                                                                                         >  1:2057103/1‑100 (MQ=255)
                       cAGAGGCACAGCGCCACCGGGTATACAGATTGTTGACAAAGTGCGCTTTGTTTATGCCGGATGCGGCGTTAACGCCTTATCCGGCCTACAAAATCGTGCa                                                                                     <  1:1726175/100‑1 (MQ=255)
                       cAGAGGCACAGCGCCACCGGGTATACAGATTGTTGACAAAGTGCGCTTTGTTTATGCCGGATGCGGCGTTAACGCCTTATCCGGCCTACAAAATCGTGCa                                                                                     >  2:1065124/1‑100 (MQ=255)
                       cAGAGGCACAGCGCCACCGGGTATACAGATTGTTGACAAAGTGCGCTTTGTTTATGCCGGATGCGGCGTTAACGCCTTATCCGGCCTACAAAATCGTGCa                                                                                     <  2:522884/100‑1 (MQ=255)
                             nacaGCGCCACCGGGTATACAGATTGTTGACAAAGTGCGCTTTGTTTATGCCGGATGCGGCGTTAACGCCTTATCCGGCCTACAAAATCGTGCAAATTCa                                                                               >  1:397039/2‑100 (MQ=255)
                             cacaGCGCCACCGGGTATACAGATTGTTGACAAAGTGCGCTTTGTTTATGCCGGATGCGGCGTTAACGCCTTATCCGGCCTACAAAATCGTGCAAATTCa                                                                               >  1:181582/1‑100 (MQ=255)
                               caGCGCCACCGGGTATACAGATTGTTGACAAAGTGCGCTTTGTTTATGCCGGATGCGGCGTTAACGCCTTATCCGGCCTACAAAATCGTGCAAATTCAAt                                                                             <  1:2364772/100‑1 (MQ=255)
                               caGCGCCACCGGGTATACAGATTGTTGACAAAGTGCGCTTTGTTTATGCCGGATGCGGCGTTAACGCCTTATCCGGCCTACAAAATCGTGCAAATTCAAt                                                                             <  2:2709499/100‑1 (MQ=255)
                                    ccACCGGGTATACAGATTGTTGACAAAGTGCGCTTTGTTTATGCCGGATGCGGCGTTAACGCCTTATCCGGCCTACAAAATCGTGCAAATTCAATACAtt                                                                        <  1:2166533/100‑1 (MQ=255)
                                     cACCGGGTATACAGATTGTTGACAAAGTGCGCTTTGTTTATGCCGGATGCGGCGTTAACGCCTTATCCGGCCTACAAAATCGTGCAAATTCAATACATTg                                                                       >  1:1263404/1‑100 (MQ=255)
                                      aCCGGGTATACAGATTGTTGACAAAGTGCGCTTTGTTTATGCCGGATGCGGCGTTAACGCCTTATCCGGCCTACAAAATCGTGCAAATTCAATACATTGc                                                                      <  1:66680/100‑1 (MQ=255)
                                        cGGGTATACAGATTGTTGACAAAGTGCGCTTTGTTTATGCCGGATGCGGCGTTAACGCCTTATCCGGCCTACAAAATCGTGCAAATTCAATACATTGC‑‑Ag                                                                 <  1:2334974/100‑1 (MQ=255)
                                          ggTATACAGATTGTTGACAAAGTGCGCTTTGTTTATGCCGGATGCGGCGTTAACGCCTTATCCGGCCTACAAAATCGTGCAAATTCAATACATTGC‑‑AGaa                                                              >  2:2846732/1‑100 (MQ=255)
                                           gTATACAGATTGTTGACAAAGTGCGCTTTGTTTATGCCGGATGCGGCGTTAACGCCTTATCCGGCCTACAAAATCGTGCAAATTCAATACATTGC‑‑AGAAt                                                             <  2:233905/100‑1 (MQ=255)
                                               aCAGATTGTTGACAAAGTGCGCTTTGTTTATGCCGGATGCGGCGTTAACGCCTTATCCGGCCTACAAAATCGTGCAAATTCAATACATTGC‑‑AGAATTTTc                                                        >  2:1407566/1‑100 (MQ=255)
                                                cAGATTGTTGACAAAGTGCGCTTTGTTTATGCCGGATGCGGCGTTAACGCCTTATCCGGCCTACAAAATCGTGCAAATTCAATACATTGC‑‑AGAATTTTC‑‑g                                                     <  1:2569771/100‑1 (MQ=255)
                                                    gcataaaCAAAGCGCACTTTGTTTATGCCGGATGCGGCGTTAACGCCTTATCCGGCCTACAAAATCGTGCAAATTCAATACATTGC‑‑AGAATTTTC‑‑GTAgg                                                 >  2:97915/7‑100 (MQ=18)
                                                         gACAAAGTGCGCTTTGTTTATGCCGGATGCGGCGTTAACGCCTTATCCGGCCTACAAAATCGTGCAAATTCAATACATTGC‑‑AGAATTTTC‑‑GTAGGCCTGa                                            >  1:1091185/1‑100 (MQ=255)
                                                           cAAAGTGCGCTTTGTTTATGCCGGATGCGGCGTTAACGCCTTATCCGGCCTACAAAATCGTGCAAATTCAATACATTGC‑‑AGAATTTTC‑‑GTAGGCCTGATa                                          <  2:477514/100‑1 (MQ=255)
                                                               gTGCGCTTTGTTTATGCCGGATGCGGCGTTAACGCCTTATCCGGCCTACAAAATCGTGCAAATTCAATACATTGC‑‑AGAATTTTC‑‑GTAGGCCTGATAAGcn                                      <  2:397038/100‑2 (MQ=255)
                                                                tGCGCTTTGTTTATGCCGGATGCGGCGTTAACGCCTTATCCGGCCTACAAAATCGTGCAAATTCAATACATTGC‑‑AGAATTTTC‑‑GTAGGCCTGATAAGCgg                                     >  2:1836634/1‑100 (MQ=255)
                                                                tGCGCTTTGTTTATGCCGGATGCGGCGTTAACGCCTTATCCGGCCTACAAAATCGTGCAAATTCAATACATTGC‑‑AGAATTTTC‑‑GTAGGCCTGATAAGCgg                                     <  1:820293/100‑1 (MQ=255)
                                                                 gcgcTTTGTTTATGCCGGATGCGGCGTTAACGCCTTATCCGGCCTACAAAATCGTGCAAATTCAATACATTGC‑‑AGAATTTTC‑‑GTAGGCCTGATAAGCGGa                                    <  2:1482778/100‑1 (MQ=255)
                                                                  cgcTTTGTTTATGCCGGATGCGGCGTTAACGCCTTATCCGGCCTACAAAATCGTGCAAATTCAATACATTGC‑‑AGAATTTTC‑‑GTAGGCCTGATAAGCGGAg                                   <  1:1779697/100‑1 (MQ=255)
                                                                   gcTTTGTTTATGCCGGATGCGGCGTTAACGCCTTATCCGGCCTACAAAATCGTGCAAATTCAATACATTGC‑‑AGAATTTTC‑‑GTAGGCCTGATAAGCGGAgc                                  <  1:216799/100‑1 (MQ=255)
                                                                    gcttttgtaggCCGGATAAGGCGTTTACGCCGCATCCGGCCTACAAAAGCATGCAAATTCAATATATTGC‑‑AGAGTTCAC‑‑GTAGGCCTGATAAGCGTAgcg                                 >  2:2197671/11‑100 (MQ=17)
                                                                    cTTTGTTTATGCCGGATGCGGCGTTAACGCCTTATCCGGCCTACAAAATCGTGCAAATTCAATACATTGC‑‑AGAATTTTC‑‑GTAGGCCTGATAAGCGGAgcg                                 >  1:433220/1‑100 (MQ=255)
                                                                      ttGTTTATGCCGGATGCGGCGTTAACGCCTTATCCGGCCTACAAAATCGTGCAAATTCAATACATTGC‑‑AGAATTTTC‑‑GTAGGCCTGATAAGCGGAGCGCa                               <  2:2571594/100‑1 (MQ=255)
                                                                       tGTTTATGCCGGATGCGGCGTTAACGCCTTATCCGGCCTACAAAATCGTGCAAATTCACTACATTGC‑‑AGAATTTTC‑‑GTAGGCCTGATAAGCGGAGCGCAt                              >  2:2594305/1‑100 (MQ=255)
                                                                       tGTTTATGCCGGATGCGGCGTTAACGCCTTATCCGGCCTACAAAATCGTGCAAATTCAATACATTGC‑‑AGAATTTTC‑‑GTAGGCCTGATAAGCGGAGCGCAt                              >  2:757035/1‑100 (MQ=255)
                                                                       tGTTTATGCCGGATGCGGCGTTAACGCCTTATCCGGCCTACAAAATCGTGCAAATTCAATACATTGC‑‑AGAATTTTC‑‑GTAGGCCTGATAAGCGGAGCGCAt                              >  1:1977426/1‑100 (MQ=255)
                                                                        ttgtaggCCGGATAAGGCGTTCGCGCCGCATCCGGCATACAAAAGTTTGCAAATTCAATAAATTGC‑‑AGAATGCAT‑‑GTAGGCCTGATAAGCGAAGCACATc                             >  2:2452989/7‑100 (MQ=12)
                                                                        gTTTATGCCGGATGCGGCGTTAACGCCTTATCCGGCCTACAAAATCGTGCAAATTCACTACATTGC‑‑AGCATTTTC‑‑GTAGGCCTGATAAGCGGAGCGCagc                             <  1:1951477/100‑3 (MQ=255)
                                                                         tgtaggCCGGATAAGGCGTTCACGCCGCATCCGGCATACAAGTTCGTGCAAATTCAATAAATTGC‑‑AATATGAC‑‑GTAGGCCTGATAAGCGTAGCGCATCAg                           >  1:747543/6‑100 (MQ=18)
                                                                          ttATGCCGGATGCGGCGTTAACGCCGTATCCGGCCTACAAAATCGTGCAAATTCAATACATTGC‑‑AGAATTTTC‑‑GTAGGCCTGATAAGCGGAGCGCATCAg                           <  2:500621/100‑1 (MQ=255)
                                                                          gtaggCCGGATAAGGCGTTTACGCCGCATCCGGCCTACAAAAGCGTGCAAATTCAATACATTGC‑‑ATGGGCCAT‑‑GTAGGCCTGATAAGCGTAGCGCATCAg                           >  2:1245512/5‑100 (MQ=17)
                                                                          gtaggCCGGATAAGGCGTTCACGCCGCATCCGGCCTACAAAATCGCTTAAATTCAATATATTGC‑‑ACCATCCTC‑‑GTAGGCCTGATAAGCGTAGCGCATCAg                           >  2:1591017/5‑100 (MQ=18)
                                                                          gtaggCCGGATAAGGCGTTCACGCCGCATCCGGCCTACAAAATCGCTTAAATTCAATATATTGC‑‑ACCATCCTC‑‑GTAGGCCTGATAAGCGTAGCGCATCAg                           >  2:2422358/5‑100 (MQ=21)
                                                                           taggCCGGATAAGGCGTTCACGCCGCATCCGGCCTACAAATCATGCAAATCCAAAATATTGC‑‑AGGAGTTAT‑‑GTAGGCCTGATAAGCGTAGCGCATCAGGc                         >  1:356265/4‑100 (MQ=21)
                                                                            nggCCGGATAAGGCGTTTACGCCGCATCCGGCCTACAAAATCGTGCAAATTCAACATATTGC‑‑AATTCTCTTGTAGGCCTGATAAGCGTAGCGCATCAGGc                         >  2:1712541/3‑100 (MQ=18)
                                                                            aggCCGGATAAGGCGTTTACGCCGCATCCGGCCTACAAAATCGTGCTAATTCAATATATTGC‑‑AGAAACCTT‑‑GTAGGCCTGATAAGCGTAGCGCATCAGGc                         >  1:2728560/3‑100 (MQ=17)
                                                                            aggCCGGATAAGGCGTTCACGCCGCATCCGGCCTACAAATCATGCAAATCCAAAATATTGC‑‑AGGAGTTAT‑‑GTAGGCCTGATAAGCGTAGCGCATCAGGCa                        >  1:2127215/3‑100 (MQ=21)
                                                                            aggCCGGATAAGGCGTTAACGCCGCATCCGGCCTACAAAATCGTGCAAATTCAATACATTGC‑‑AGAATTTTC‑‑GTAGGCCTGATAAGCGGAGCGCATCAGGc                         >  2:2019591/3‑100 (MQ=18)
                                                                             ggCCGGATAAGGCGTTTACGCCGCATCCGGCCTACAAAATCGTGCAAATTCAACATATTGC‑‑AATTCTCTT‑‑GTAGGCCTGATAAGCGTAGCGCATCAGGct                        >  1:1626519/2‑99 (MQ=18)
                                                                             ggCCGGATAAGGCGTTCACGCCGCATCCGGCCTACAAAAATCTTGCCAATTCAATATATTGC‑‑AGGACCAT‑‑GTAGGCCTGATAAGCGTAGCGCATCAGGCa                        >  1:2647303/2‑100 (MQ=18)
                                                                             ggCCGGATAAGGCGTTCACGCCGCATCCGGCCTACAAGATCGTGCAAATTCAATATATTGC‑‑ATGAGCCAC‑‑GTAGGCCTGATAAGCGTAGCGCATCAGGCa                        >  1:1055559/2‑100 (MQ=17)
                                                                             ggCCGGATAAGGCGTTCACGCCGCATCCGGCCTACAAAGTCGTATAAATTCAATGTATTGC‑‑AGTTTATTG‑‑GTAGGCCTGATAAGCGCAGCGCATCAGGCa                        >  1:2689107/2‑100 (MQ=25)
                                                                             ggCCGGATAAGGCGTTCACGCCGCATCCGGCCTACAAAAGCACGCAAACTCAATATATTGC‑‑AGAGATCAT‑‑GTAGGCCTGATAAGCGTAGCGCATCAGGCa                        >  1:1217599/2‑100 (MQ=17)
                                                                             ggCCGGATAAGGCGTTCACGCCGCATCCGGCCTACAAAAGCACGCAAACTCAATATATTGC‑‑AGAGATCAT‑‑GTAGGCCTGATAAGCGTAGCGCATCAGGCa                        >  2:2614083/2‑100 (MQ=17)
                                                                             ggCCGGATAAGGCGTTCACGCCGCATCCGGCCTACAAAACCTTGCTAATTCAATATATTGC‑‑AGGGACTAT‑‑GTAGGCCTGATAAGCATAGCGCATCAGGCa                        >  2:2094932/2‑100 (MQ=18)
                                                                             ggCCGGATAAGGCGTTCACGCCGCATCCGGCCTACAAATCATGCAAATCCAAAATATTGC‑‑AGGAGTTAT‑‑GTAGGCCTGATAAGCGTAGCGCATCAGGCAt                       >  1:2828885/2‑100 (MQ=21)
                                                                             ggCCGGATAAGGCGTTAACGCCGCATCCGGCCTACAAAATCGTGCAAATTCAATACTTTGC‑‑AGAATTTTC‑‑GTAGGCCTGATAAGCGGAGAGCATCAGGc                         >  2:2493784/2‑99 (MQ=17)
                                                                              gCCGGATGCGGCGTTAACGCCTTATCCGGCCTACAAAATCGTGCAAATTCAATACATTGC‑‑AGAATTTTC‑‑GTAGGCCTGATAAGCGGAGCGCATCAGGCAt                       >  1:310663/1‑100 (MQ=255)
                                                                              gCCGGATAAGGCGTTTACGCCGCATCCGGCCTACAAATTCTGCAAATTCAACAAATTGC‑‑GAATCCCTT‑‑GTAGGCCTGATAAGCGTAGCGCATCAGGCAtt                      >  1:2079449/1‑100 (MQ=18)
                                                                              gCCGGATAAGGCGTTTACGCCGCATCCGGCCTACAAAAGCATGCAAATTCAATATATTGC‑‑AGAGTTCAC‑‑GTAGGCCTGATAAGCGTAGCGCATCAGGCag                       >  1:806694/1‑99 (MQ=17)
                                                                              gCCGGATAAGGCGTTTACGCCGCATCCGGCCTACAAAAGCATGCAAATTCAATATATTGC‑‑AGAGTTCAC‑‑GTAGGCCTGATAAGCGTAGCGCATCAGGCag                       >  2:1200681/1‑99 (MQ=17)
                                                                              gCCGGATAAGGCGTTCACGCCGCATCCGGCCTACGAAATCGTGCTAATTCAATATATTGT‑‑GGGCTCCAC‑‑GTAGGCCTGATAAGCGTAGCGCATCAGGCag                       >  2:638621/1‑99 (MQ=21)
                                                                              gCCGGATAAGGCGTTCACGCCGCATCCGGCCTACAAAATCGCTTAAATTCAATATATTGC‑‑ACCATCCTC‑‑GTAGGCCTGATAAGCGTAGCGCATCAGGCaa                       >  1:1572285/1‑99 (MQ=21)
                                                                              gCCGGATAAGGCGTTCACGCCGCATCCGGCCTACAAAACTGTGCAAATTCAATATATTGC‑‑AGGAAACAC‑‑GTAGGCCTGATAAGCGAAGCGCATCAGGCag                       >  1:879250/1‑99 (MQ=18)
                                                                              gCCGGATAAGGCGTTCACGCCGCATCCGGCCTACAAACCGCGCAAATTCAATATATTGCGGAGAAAAT‑‑‑‑GTAGGCCTGATAAGCGTAGCGCATCAGGctgt                      >  2:2659404/1‑97 (MQ=21)
                                                                               ccggataaggcgTTCACGCCGCATCCGGCCTACAAAACTGTGCAAATTCAATATATTGC‑‑AGGAAACAC‑‑GTAGGCCTGATAAGCGAAGCGCATCAGGCagt                      >  1:990819/9‑98 (MQ=18)
                                                                               ccGGATGCGGCGTTAACGCCTTATCCGGCCTACGAAATCGTGCAAATTCAATACATTGC‑‑AGAATTTTC‑‑GTAGGCCTGATAAGCGGAGCGGATCAGGCAtt                      >  2:1894587/1‑100 (MQ=255)
                                                                               ccGGATAAGGCGTTTACGCCGCATCCGGCCTACAAAATCGTGCTAATTCAATATATTGC‑‑AGAAACCTT‑‑GTAGGCCTGATAAGCGTAGCGCATCAGGCag                       >  1:2224785/1‑98 (MQ=17)
                                                                                cggataaggcgTTTACGCCGCATCCGGCATACAAAAGCATGCAAATTCAATATATTGC‑‑AGAGTTCAC‑‑GTAGGCCTGATAAGCGTAGCGCATCAGGCagtt                     >  2:1429748/8‑97 (MQ=17)
                                                                                cggataaggcgTTCACGCCGCATCCGGCCTACAAAATCGCTTAAATTCAATATATTGC‑‑ACCATCCTC‑‑GTAGGCCTGATAAGCGTAGCGCATCAGGCaatt                     >  2:720531/8‑97 (MQ=21)
                                                                                cgatctaggcgTTCACGCCGCATCCGACCCACAAAGTCTTGCAAAATCAATAAATTGC‑‑AGGAACGAT‑‑GTAGGCCTGATAAGCGTAGCGCATCAGGCaact                     <  1:1146744/93‑4 (MQ=21)
                                                                                cgatctaggcgTTCACGCCGCATCCGACCCACAAAGTCTTGCAAAATCAATAAATTGC‑‑AGGAACGAT‑‑GTAGGCCTGATAAGCGTAGCGCATCAGGCaact                     <  1:2597716/93‑4 (MQ=25)
                                                                                 ggataaggcgTTTACGCCGCATCCGGCCTACAAAATCGTGCTAATTCAATATATTGC‑‑AGAAACCTT‑‑GTAGGCCTGATAAGCGTAGCGCATCAGGCAgttt                    >  1:2001385/7‑100 (MQ=18)
                                                                                 ggataaggcgTTTACGCCGCATCCGGCCTACAAAAGCATGCAAATTCAATATATTGC‑‑AGAGTTCAC‑‑GTAGGCCTGATAAGCGTAGCGCATCAGGCAgttt                    >  1:2457023/7‑100 (MQ=17)
                                                                                 ggataaggcgTTCACGCCGCATCCGGCCTACAAAACCTTGCTAATTCAATATATTGC‑‑AGGGACTAT‑‑GTAGGCCTGATAAGCATAGCGCATCAGGCagctt                    >  1:1891837/7‑96 (MQ=18)
                                                                                 ggataaggcgTTCACGCCGCATCCGGCCTACAAATCATGCAAATCCAAAATATTGC‑‑AGGAGTTAT‑‑GTAGGCCTGATAAGCGTAGCGCATCAGGCAtttat                   >  2:2271105/7‑98 (MQ=21)
                                                                                 ggataaggcgTTCACGCCGCATCCGGCCTACAAATCATGCAAATCCAAAATATTGC‑‑AGGAGTCAT‑‑GTACGCCTGATACGCGTAGCGCATCAGGCAtttat                   >  2:1952617/7‑98 (MQ=21)
                                                                                  gataaggcgTTTACGCCTTATCCGGCCTACAAATTCTTCCAAATTCAATATATTGC‑‑AGAAATCAT‑‑GTAGGCCTGATAAGCGTAGCGCATCAGGCAatttt                   >  2:1231712/6‑100 (MQ=18)
                                                                                  gataaggcgTTTACGCCGCATCCGGCCTACAAAACAGTGCAAATTCAATATATTGC‑‑AATCTTCTG‑‑GTAGGCCTGATAAGCATCGCGCATCAGGCaagctt                   >  1:2242303/6‑95 (MQ=21)
                                                                                  gataaggcgTTCACGCCGCATCCGGCCTACAAAATCGCTTAAATTCAATATATTGC‑‑ACCATCCTC‑‑GTAGGCCTGATAAGCGTAGCGCATCAGGCaattct                   >  1:1641059/6‑95 (MQ=21)
                                                                                   aTGCGGCGTTAACGCCTTATCCGGCCTACAAAATCGTGCAAATTCAATACATTGC‑‑AGAATTTTC‑‑GTAGGCCTGATAAGCGGAGCGCATCAGGCATTtttg                  <  2:2833220/100‑1 (MQ=255)
                                                                                    taaggcgTTCACGCCGCATCCGGCCTACAAAATCGCTTAAATTCAATATATTGC‑‑ACCATCCTC‑‑GTAGGCCTGATAAGCGTAGCGCATCAGGCAattctgt                 >  2:1276733/4‑96 (MQ=21)
                                                                                     aaggcgTTCACGCCGCATCCGGCCTACAAAGTCGTATAAATTCAATGTATTGC‑‑AGTTTATTG‑‑GTAGGCCTGATAAGCGCAGCGCATCAGGCaatcttacg                <  1:1327297/98‑9 (MQ=37)
                                                                                         cgTTAACGCCTTATCCGGCCTACAAAATCGTGCAAATTCAATACATTGC‑‑AGAATTTTC‑‑GTAGGCCTGATAAGCGGAGCGCATCAGGCATTtttgcctttg            <  2:1091182/100‑1 (MQ=255)
                                                                                         cgTTAACGCCTTATCCGGCCTACAAAATCGTGCAAATTCAATACATTGC‑‑AGAATTTTC‑‑GTAGGCCTGATAAGCGGAGCGCATCAGGCATTtttgcctttg            <  2:1630585/100‑1 (MQ=255)
                                                                                           ttAACGCCTTATCCGGCCTACAAAATCGTGCAAATTCAATACATTGC‑‑AGAATTTTC‑‑GTAGGCCTGATAAGCGGAGCGCATCAGGCATTtttgcctttgcc          >  1:134289/1‑100 (MQ=255)
                                                                                                 ccTTATCCGGCCTACAAAATCGTGCAAATTCAATACATTGC‑‑AGAATTTTC‑‑GTAGGCCTGATAAGCGGAGCGCATCAGGCATTTTTGCCTTTGCCATc       >  2:1759215/1‑97 (MQ=255)
                                                                                                  cTTATCCGGCCTACAAAATCGTGCAAATTCAATACATTGC‑‑AGAATTTTC‑‑GTAGGCCTGATAAGCGGAGCGCATCAGGCATTTTTGCCTTTGCCATCAGTc   <  2:1127792/100‑1 (MQ=255)
                                                                                                  cTTATCCGGCCTACAAAATCGTGCAAATTCAATACATTGC‑‑AGAATTTTC‑‑GTAGGCCTGATAAGCGGAGCGCATCAGGCATTTTTGCCTTTGCCATCAGTc   <  2:769853/100‑1 (MQ=255)
                                                                                                   ttATCCGGCCTACAAGATCGTGCAAATTCAATACATTGC‑‑AGAATTTTC‑‑GTAGGCCTGATAAGCGGAGCGCATCAGGCATTTTTGCCTTTGCCATCAGTCt  <  1:2846753/100‑1 (MQ=255)
                                                                                                   |                                                                                                           
TTGCTGCCGGATGGCGCTGCACGCAGAGGCACAGCGCCACCGGGTATACAGATTGTTGACAAAGTGCGCTTTGTTTATGCCGGATGCGGCGTTAACGCCTTATCCGGCCTACAAAATCGTGCAAATTCAATACATTGC‑‑AGAATTTTC‑‑GTAGGCCTGATAAGCGGAGCGCATCAGGCATTTTTGCCTTTGCCATCAGTCT  >  NC_000913/3705885‑3706083

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