Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009133,706,0170CT32.9% 334.5 / 191.2 231intergenic (‑227/+81)dppB/dppAdipeptide/heme ABC transporter permease/dipeptide/heme ABC transporter periplasmic binding protein; dipeptide chemotaxis receptor
Reads supporting (aligned to +/- strand):  ref base C (79/66);  new base T (69/2);  total (162/69)
Fisher's exact test for biased strand distribution p-value = 3.10e-12
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 8.24e-01

Only 100 of 231 total aligned reads displayed.

CACCGGGTATACAGATTGTTGACAAAGTGCGCTTTGTTTATGCCGGATGCGGCGTTAACGCCTTATCCGGCCTACAAAATCGTGCAAATTCAATACATTGC‑‑AGAATT‑‑‑TTC‑‑GTAGGCCTGATAAGCGGAGCGCATCAGGCATTTTTGCCTTTGCCATCAGTCTTGTATGGCTTTTAATTATTCGATAGAGACGTT  >  NC_000913/3705922‑3706115
                                                                                                |                                                                                                          
cACCGGGTATACAGATTGTTTTTATATTGTGTTTTTTTTATGCCGGATGCGGCGTTAACGCCTTATCTGGCCTACAAAATCGTGCAAATTCAATACAtt                                                                                                         >  2:1357138/1‑99 (MQ=17)
   cGGGTATACAGATTGTTGACAAAGTGCGCTTTGTTTATGCCGGATGCGGCGTTAACGCCTTATCCGGCCTACAAAATCGTGCAAATTCAATACATTGC‑‑Ag                                                                                                  >  1:133567/1‑100 (MQ=255)
    gggTATACAGATTGTTGACAAAGTGCGCTTTGTTTATGCCGGATGCGGCGTTAACGCCTTATCCGGCCTACAAAATCGTGCAAATTCAATACATTGC‑‑AGa                                                                                                 >  1:2703417/1‑100 (MQ=255)
      gTATACAGATTGTTGACAAAGTGCGCTTTGTTTATGCCGGATGCGGCGTTAACGCCTTATCCGGCCTACAAAATCGTGCAAATTCAATACATTGC‑‑AGAAt                                                                                               <  2:2097538/100‑1 (MQ=255)
       tataCAGATTGTTGACAAAGTGCGCTTTGTTTATGCCGGATGCGGCGTTAACGCCTTATCCGGCCTACAAAATCGTGCAAATTCAATACATTGC‑‑AGAAtt                                                                                              >  2:2120147/1‑100 (MQ=255)
        ataCAGATTGTTGACAAAGTGCGCTTTGTTTATGCCGGATGCGGCGTTAACGCCTTATCCGGCCTACAAAATCGTGCAAATTCAATACATTGC‑‑AGAAtt‑‑‑t                                                                                          <  1:1941778/100‑1 (MQ=255)
          aCAGATTGTTGACAAAGTGCGCTTTGTTTATGCCGGATGCGGCGTTAACGCCTTATCCGGCCTACAAAATCGTGCAAATTCAATACATTGC‑‑AGAATT‑‑‑TTc                                                                                        >  1:1169899/1‑100 (MQ=255)
           cAGATTGTTGACAAAGTGCGCTTTGTTTATGCCGGATGCGGCGTTAACGCCTTATCCGGCCTACAAAATCGTGCAAATTCAATACATTGC‑‑AGAATT‑‑‑TTC‑‑g                                                                                     <  2:1381745/100‑1 (MQ=255)
           cAGATTGTTGACAAAGTGCGCTTTGTTTATGCCGGATGCGGCGTTAACGCCTTATCCGGCCTACAAAATCGTGCAAATTCAATACATTGC‑‑AGAATT‑‑‑TTC‑‑g                                                                                     <  1:577068/100‑1 (MQ=255)
           cAGATTGTTGACAAAGTGCGCTTTGTTTATGCCGGATGCGGCGTTAACGCCTTATCCGGCCTACAAAATCGTGCAAATTCAATACATTGC‑‑AGAATT‑‑‑TTC‑‑g                                                                                     <  2:22327/100‑1 (MQ=255)
             gATTGTTGACAAAGTGCGCTTTGTTTATGCCGGATGCGGCGTTAACGCCTTATCCGGCCTACAAAATCGTGCAAATTCAATACATTGC‑‑AGAATT‑‑‑TTC‑‑GTa                                                                                   <  2:1908849/100‑1 (MQ=255)
              aTTGTTGACAAAGTGCGCTTTGTTTATGCCGGATGCGGCGTTAACGCCTTATCCGGCCTACAAAATCGTGCAAATTCAATACATTGC‑‑AGAATT‑‑‑TTC‑‑GTAg                                                                                  >  1:2924171/1‑100 (MQ=255)
               ttgttgACAAAGTGCGCTTTGTTTATGCCGGATGCGGCGTTAACGCCTTATCCGGCCTACAAAATCGTGCAAATTCAATACATTGC‑‑AGAATT‑‑‑TTC‑‑GTAgg                                                                                 >  1:2813786/1‑100 (MQ=255)
               ttgttgACAAAGTGCGCTTTGTTTATGCCGGATGCGGCGTTAACGCCTTATCCGGCCTACAAAATCGTGCAAATTCAATACATTGC‑‑AGAATT‑‑‑TTC‑‑GTAgg                                                                                 <  1:325121/100‑1 (MQ=255)
                 gttgACAAAGTGCGCTTTGTTTATGCCGGATGCGGCGTTAACGCCTTATCCGGCCTACAAAATCGTGCAAATTCAATACATTGC‑‑AGAATT‑‑‑TTC‑‑GTAGGcc                                                                               >  1:715399/1‑100 (MQ=255)
                  ttgACAAAGTGCGCTTTGTTTATGCCGGATGCGGCGTTAACGCCTTATCCGGCCTACAAAATCGTGCAAATTCAATACATTGC‑‑AGAATT‑‑‑TTC‑‑GTAGGCCt                                                                              <  2:1126611/100‑1 (MQ=255)
                  ttgACAAAGTGCGCTTTGTTTATGCCGGATGCGGCGTTAACGCCTTATCCGGCCTACAAAATCGTGCAAATTCAATACATTGC‑‑AGAATT‑‑‑TTC‑‑GTAGGCCt                                                                              >  2:618854/1‑100 (MQ=255)
                  ttgACAAAGTGCGCTTTGTTTATGCCGGATGCGGCGTTAACGCCTTATCCGGCCTACAAAATCGTGCAAATTCAATACATTGC‑‑AGAATT‑‑‑TTC‑‑GTAGGCCt                                                                              >  2:2473709/1‑100 (MQ=255)
                  ntgACAAAGTGCGCTTTGTTTATGCCGGATGCGGCGTTAACGCCTTATCCGGCCTACAAAATCGTGCAAATTCAATACATTGC‑‑AGAATT‑‑‑TTC‑‑GTAGGCCt                                                                              >  1:2180883/2‑100 (MQ=255)
                   tgACAAAGTGCGCTTTGTTTATGCCGGATGCGGCGTTAACGCCTTATCCGGCCTACAAAATCGTGCAAATTCAATACATTGC‑‑AGAATT‑‑‑TTC‑‑GTAGGCCTg                                                                             >  2:177356/1‑100 (MQ=255)
                      cAAAGTGCGCTTTGTTTATGCCGGATGCGGCGTTAACGCCTTATCCGGCCTACAAAATCGTGCAAATTCAATACATTGC‑‑AGAATT‑‑‑TTC‑‑GTAGGCCTGATa                                                                          <  2:1940631/100‑1 (MQ=255)
                         aGTGCGCTTTGTTTATGCCGGATGCGGCGTTAACGCCTTATCCGGCCTACAAAATCGTGCAAATTCAATACATTGC‑‑AGAATT‑‑‑TTC‑‑GTAGGCCTGATAAGc                                                                       >  1:52517/1‑100 (MQ=255)
                          gTGCGCTTTGTTTATGCCGGATGCGGCGTTAACGCCTTATCCGGCCTACAAAATCGTGCAAATTCAATACATTGC‑‑AGAATT‑‑‑TTC‑‑GTAGGCCTGATAAGCg                                                                      >  2:276477/1‑100 (MQ=255)
                           tGCGCTTTGTTTATGCCGGATGCGGCGTTAACGCCTTATCCGGCCTACAAAATCGTGCAAATTCAATACATTGC‑‑AGAATT‑‑‑TTC‑‑GTAGGCCTGATAAGCgg                                                                     >  1:905492/1‑100 (MQ=255)
                                tttGTTTATGCCGGATGCGGCGTTAACGCCTTATCCGGCCTACAAAATCGTGCAAATTCAATACATTGC‑‑AGAATT‑‑‑TTC‑‑GTAGGCCTGATAAGCGGAgcgc                                                                >  2:2034121/1‑100 (MQ=255)
                                 ttGTTTATGCCGGATGCGGCGTTAACGCCTTATCCGGCCTACAAAATCGTGCAAATTCAATACATTGC‑‑AGAATT‑‑‑TTC‑‑GTAGGCCTGATAAGCGGAGCGCa                                                               <  2:589869/100‑1 (MQ=255)
                                  tttgtaggCCGGATAAGGCGTTTACGCCGCATCCGGCCTACAAAAGCGTGCAAATTCAATACATTGC‑‑ATGGGC‑‑‑CAT‑‑GTAGGCCTGATAAGCGTAGCGCAt                                                              >  1:2202428/8‑100 (MQ=17)
                                   ttgtaggCCGGATAAGGCGTTAACGCCGCATCCGGCATACAAAATCATGCAAATTCAATACATTGC‑‑AGAATT‑‑‑TTC‑‑GTAGGCCTGATAAGATGAGCGCATc                                                             >  2:2754846/7‑100 (MQ=14)
                                   gTTTATGCCGGATGCGGCGTTAACGCCTTATCCGGCCTACAAAATCGTGCAAATTCAATACATTGC‑‑AGAATT‑‑‑TTC‑‑GTAGGCCTGATAAGCGGAGCGCATc                                                             <  2:3025709/100‑1 (MQ=255)
                                   atctaggCCGGATAAGGCGTTTACGCCGCATCCGGCCTACAAAATCGTGCTAATTCAATATATTGC‑‑AGAAAC‑‑‑CTT‑‑GTAGGCCTGATAAGCGTAGCGCATc                                                             <  2:1615826/94‑1 (MQ=18)
                                    tgtaggCCGGATAAGGCGTTCACGCCGCATCCGGCCTACAAGATCGTGCAAATTCAATATATTGC‑‑ATGAGC‑‑‑CAC‑‑GTAGGCCTGATAAGCGTAGCGCATCa                                                            >  1:497460/6‑100 (MQ=17)
                                    tgtaggCCGGATAAGGCGTTCACGCCGCATCCGGCCTACAAGATCGTGCAAATTCAATATATTGC‑‑ATGAGC‑‑‑CAC‑‑GTAGGCCTGATAAGCGTAGCGCATCa                                                            >  1:1439528/6‑100 (MQ=17)
                                     ttATGCCGGATGCGGCGTTAACGCCTTATCCGGCCTACAAAATCGTGCAAATTCAATACATTGC‑‑AGAATT‑‑‑TTC‑‑GTAGGCCTGATAAGCGGAGCGCATCAg                                                           <  2:1169304/100‑1 (MQ=255)
                                     gtaggCCGGATAAGGCGTTCACGCCGCATCCGGCCTACAAAAACCTGCAAATTCAATATATTGC‑‑AGGATA‑‑‑CGT‑‑GTAGGCCTGATAAGCGTAGCGCATCAg                                                           >  1:1056087/5‑100 (MQ=17)
                                      taggCCGGATAAGGCGTTTACGCCGCATCCGGCCTACAAAATCGTGCTAATTCAATATATTGC‑‑AGAAAC‑‑‑CTT‑‑GTAGGCCTGATAAGCGTAGCGCATCAgg                                                          >  2:1626141/4‑100 (MQ=18)
                                      taggCCGGATAAGGCGTTCACGCCGCATCCGGCCTACAAAGTCGTATAAATTCAATGTATTGC‑‑AGTTTA‑‑‑TTG‑‑GTAGGCCTGATAAGCGCAGCGCATCAgg                                                          >  1:175202/4‑100 (MQ=21)
                                       aggCCGGATAAGGCGTTTACGCCGCATCTGGCCTACAAATTCTTCCAAATTCAATATATTGC‑‑AGAAAT‑‑‑CAT‑‑GTAGGCCTGATAAGCGTAGCGCATCAGGc                                                         >  2:1534168/3‑100 (MQ=18)
                                       aggCCGGATAAGGCGTTCACGCCGCATCCGGCCTACAAAAATCTTGCCAATTCAATATATTGC‑‑AGGACC‑‑‑AT‑‑GTAGGCCTGATAAGCGTAGCGCATCAGGc                                                         >  2:1845514/3‑100 (MQ=18)
                                       aggCCGGATAAGGCGTTCACGCCGCATCCGGCCTACAAGTTCGTGCAAATTCAATAAATTGC‑‑AATAT‑‑‑GAC‑‑GTAGGCCTGATAAGCGTAGCGCATCAGGcg                                                        >  1:1675086/3‑99 (MQ=21)
                                       aggCCGGATAAGGCGTTCACGCCGCATCCGGCCTACAAAACTGTGCAAATTCAATATATTGC‑‑AGGAAA‑‑‑CAC‑‑GTAGGCCTGATAAGCGAAGCGCATCAGGc                                                         >  1:2272655/3‑100 (MQ=18)
                                       aggCCGGATAAGGCGTTCACGCCGCATCCGGCATACAAAAGCATGAAAATTCAATATATTGC‑‑AGGAGC‑‑‑TGC‑‑GTAGGCCTGATAAGCGTAGCGCATCAGGc                                                         >  1:2800464/3‑100 (MQ=18)
                                       aggCCGGATAAGGCGTTCACGCCGCATCCGGCATACAAACCGCGCAAATTCAATATATTGCGGAGAAAA‑‑‑T‑‑‑‑GTAGGCCTGATAAGCGTAGCGCATCAGGct                                                        >  2:2148907/3‑99 (MQ=18)
                                        tGCCGGATGCGGCGTTAACGCCTTATCCGGCCTACAAAATCGTGCAAATTCAATACATTGC‑‑AGAATT‑‑‑TTC‑‑GTAGGCCTGATAAGCGGAGCGCATCAGGCa                                                        >  2:456903/1‑100 (MQ=255)
                                        ggCCGGATAAGGCGTTTACGCCGCATCCGGCCTACAAAATCGTGCAAATTCAACATATTGC‑‑AATT‑‑‑CTCTTGTAGGCCTGATAAGCGTAGCGCATCAGGct                                                        >  1:1492049/2‑99 (MQ=18)
                                        ggCCGGATAAGGCGTTTACGCCGCATCCGGCCTACAAAAGCGTGCAAATTCAATACATTGC‑‑ATGGGCC‑‑‑AT‑‑GTAGGCCTGATAAGCGTAGCGCATCAGGCa                                                        >  1:873562/2‑100 (MQ=17)
                                        ggCCGGATAAGGCGTTCACGCCGCATCCGGCCTACAAAGTCGTATAAATTCAATGTATTGC‑‑AGTTTA‑‑‑TTG‑‑GTAGGCCTGATAAGCGCAGCGCATCAGGCa                                                        >  1:2752989/2‑100 (MQ=21)
                                        ggCCGGATAAGGCGTTCACGCCGCATCCGGCCTACAAAAGCATGAAAATTCAATATATTGC‑‑AGGAGC‑‑‑TGC‑‑GTAGGCCTGATAAGCGTAGCGCATCAGGCa                                                        >  1:2340239/2‑100 (MQ=18)
                                        ggCCGGATAAGGCGTTCACGCCGCATCCGGCCTACAAAACTTTGCAAATTCAATATATTGC‑‑AATC‑‑‑TCCGTGTAGGCCTGATAAGCGTAGCGCATCAGGCa                                                        >  2:1103927/2‑100 (MQ=18)
                                        ggCCGGATAAGGCGTTCACGCCGCATCCGGCCTACAAAACTTTGCAAATCCAATATATTGC‑‑AATC‑‑‑TCCGTGTAGGCCTGATAAGCGTAGCGCATCAGGCa                                                        >  1:1275579/2‑100 (MQ=18)
                                        ggCCGGATAAGGCGTTCACGCCGCATCCGGCCTACAAATCATGCAAATCCAAAATATTGC‑‑AGGAGT‑‑‑TAT‑‑GTAGGCCTGATAAGCGTAGCGCATCAGGCAt                                                       >  1:1882516/2‑100 (MQ=21)
                                         gCCGGATGCGGCGTTAACGCCTTATCCGGCCTACAAAATCGTGCAAATTCAATACATTGC‑‑AGAATT‑‑‑TTC‑‑GTAGGCCTGATAAGCGGAGCGCATCAGGCAt                                                       >  1:595957/1‑100 (MQ=255)
                                         gCCGGATAAGGCGTTTACGCCGCATCCGGCCTACAAAATCGTGCTAATTCAATATATTGC‑‑AGAAAC‑‑‑CTT‑‑GTAGGCCTGATAAGCGTAGCGCATCAGGCag                                                       >  2:2413821/1‑99 (MQ=18)
                                         gCCGGATAAGGCGTTTACGCCGCATCCGGCCTACAAAAGCATGCAAATTCAATATATTGC‑‑AGAGTT‑‑‑CAC‑‑GTAGGCCTGATAAGCGTAGCGCATCAGGCag                                                       >  1:1824634/1‑99 (MQ=14)
                                         gCCGGATAAGGCGTTCACGCCGCATCCGGCCTACAAAATCGCTTAAATTCAATATATTGC‑‑ACCATC‑‑‑CTC‑‑GTAGGCCTGATAAGCGTAGCGCATCAGGCaa                                                       >  2:2248597/1‑99 (MQ=21)
                                         gCCGGATAAGGCGTTCACGCCGCATCCGGCCTACAAAAGCATGAAAATTCAATATATTGC‑‑AGGAGC‑‑‑TGC‑‑GTAGGCCTGATAAGCGTAGCGCATCAGGCag                                                       >  1:2104094/1‑99 (MQ=17)
                                         gCCGGATAAGGCGTTCACGCCGCATCCGGCCTACAAAAGCATGAAAATTCAATATATTGC‑‑AGGAGC‑‑‑TGC‑‑GTAGGCCTGATAAGCGTAGCGCATCAGGCag                                                       >  2:2967001/1‑99 (MQ=18)
                                         gCCGGATAAGGCGTTCACGCCGCATCCGGCCTACAAAACCTTGCTAATTCAATATATTGC‑‑AGGGAC‑‑‑TAT‑‑GTAGGCCTGATAAGCATAGCGCATCAGGCag                                                       >  1:1545845/1‑99 (MQ=18)
                                          ccggataaggcgTTCACGCCGCATCCGGCCTACAAAACTTTGCAAATTCAATATATTGC‑‑AATC‑‑‑TCCGTGTAGGCCTGATAAGCGTAGCGCATCAGGCaat                                                      >  2:2285866/9‑98 (MQ=18)
                                          ccggataaggcgTTCACGCCGCATCCGGCCTACAAAACTTTGCAAATTCAATATATTGC‑‑AATC‑‑‑TCCGTGTAGGCCTGATAAGCGTAGCGCATCAGGCaat                                                      >  2:2124846/9‑98 (MQ=18)
                                          ccggataaggcgTTCACGCCGCATCCGGCCTACAAAACCTTGCTAATTCAATATATTGC‑‑AGGGACT‑‑‑AT‑‑GTAGGCCTGATAAGCATAGCGCATCAGGCagc                                                      >  1:1232383/9‑98 (MQ=18)
                                            ggataaggcgTTTACGCCGCATCCGGCCTACAAAATCGTGCTAATTCAATATATTGC‑‑AGAAAC‑‑‑CTT‑‑GTAGGCCTGATAAGCGTAGCGCATCAGGCAgttt                                                    >  2:2339514/7‑100 (MQ=18)
                                            ggataaggcgTTTACGCCGCATCCGGCCTACAAAAGCATGCAAATTCAATATATTGC‑‑AGAGTT‑‑‑CAC‑‑GTAGGCCTGATAAGCGTAGCGCATCAGGCAgttt                                                    >  1:159045/7‑100 (MQ=17)
                                            ggataaggcgTTCACGCCGCATCCGGCCTACAAGATCGTGCAAATTCAATATATTGC‑‑ATGAGC‑‑‑CAC‑‑GTAGGCCTGATAAGCGTAGCGCATCAGGCAgttt                                                    >  1:2602746/7‑100 (MQ=17)
                                            ggataaggcgTTCACGCCGCATCCGGCCTACAAAACCTTGCTAATTCAATATATTGC‑‑AGGGACT‑‑‑AT‑‑GTAGGCCTGATAAGCATAGCGCATCAGGCagctt                                                    >  2:570761/7‑96 (MQ=21)
                                            ggataaggcgTTCACGCCGCATCCGGCCTACAAAACCTTGCTAATTCAATATATTGC‑‑AGGGACT‑‑‑AT‑‑GTAGGCCTGATAAGCATAGCGCATCAGGCagctt                                                    >  1:356449/7‑96 (MQ=21)
                                            ggataaggcgTTCACGCCGCATCCGGCCTACAAATCATGCAAATCCAAAATATTGC‑‑AGGAGT‑‑‑TAT‑‑GTAGGCCTGATAAGCGTAGCGCATCAGGCAtttat                                                   >  2:1136070/7‑98 (MQ=21)
                                             gataaggcgTTCACGCCGCATCCGGCCTACAAAATCGCTTAAATTCAATATATTGC‑‑ACCATC‑‑‑GTC‑‑GTAGGCCTGATAAGCGTAGCGCATCAGGCAattct                                                   >  1:1012516/6‑98 (MQ=21)
                                             gataaggcgTTCACGCCGCATCCGGCCTACAAAACCTTGCTAATTCAATATATTGC‑‑AGGGAC‑‑‑TAT‑‑GTAGGCCTGATAAGCATAGCGCATCAGGCagcttt                                                   >  2:169720/6‑95 (MQ=18)
                                                 cggcgTTAACGCCTTATCCGGCCTACAAAATCGTGCAAATTCAATACATTGC‑‑AGAATT‑‑‑TTC‑‑GTAGGCCTGATAAGCGGAGCGCATCAGGCATTtttgcct                                               >  1:630125/1‑100 (MQ=255)
                                                  ggcgTTAACGCCTTATCCGGCCTACAAAATCGTGCAAATTCAATACATTGC‑‑AGAATT‑‑‑TTC‑‑GTAGGCCTGATAAGCGGAGCGCATCAGGCATTtttgcctt                                              >  1:201473/1‑100 (MQ=255)
                                                   gcgTTTACGCCGCATCCGGCCTACAAAAATGTGCAAATTCAATAAATTGC‑‑AATTCAACTT‑‑GTAGGCCTGATAAGCGCAGCGCATCAGGCAatttggcaga                                              >  1:1365975/1‑94 (MQ=21)
                                                      ttAACGCCTTATCCGGCCTACAAAATCGTGCAAATTCAATACATTGC‑‑AGAATT‑‑‑TTC‑‑GTAGGCCTGATAAGCGGAGCGCATCAGGCATTtttgcctttgcc                                          <  1:744428/100‑1 (MQ=255)
                                                       tAACGCCTTATCCGGCCTACAAAATCGTGCAAATTCAATACATTGC‑‑AGAATT‑‑‑TTC‑‑GTAGGCCTGATAAGCGGAGCGCATCAGGCATTTTTGCCTTTGCCa                                         >  2:1334970/1‑100 (MQ=255)
                                                       tAACGCCTTATCCGGCCTACAAAATCGTGCAAATTCAATACATTGC‑‑AGAATT‑‑‑TTC‑‑GTAGGCCTGATAAGCGGAGCGCATCAGGCATTTTTGCCTTTGCCa                                         <  1:1100165/100‑1 (MQ=255)
                                                         aCGCCTTATCCGGCCTACAAAATCGTGCAAATTCAATACATTGC‑‑AGAATT‑‑‑TTC‑‑GTAGGCCTGATAAGCGGAGCGCATCAGGCATTTTTGCCTTTGCCATc                                       <  2:1750927/100‑1 (MQ=255)
                                                          cGCCTTATCCGGCCTACAAAATCGTGCAAATTCAATACATTGC‑‑AGAATT‑‑‑TTC‑‑GTAGGCCTGATAAGCGGAGCGCATCAGGCATTTTTGCCTTTGCCATCa                                      >  1:2167932/1‑100 (MQ=255)
                                                           gCCTTATCCGGCCTACAAAATCGTGCAAATTCAATACATTGC‑‑AGAATT‑‑‑TTC‑‑GTAGGCCTGATAAGCGGAGCGCATCAGGCATTTTTGCCTTTGCCATCAg                                     <  2:2514148/100‑1 (MQ=255)
                                                           gCCTTATCCGGCCTACAAAATCGTGCAAATTCAATACATTGC‑‑AGAATT‑‑‑TTC‑‑GTAGGCCTGATAAGCGGAGCGCATCAGGCATTTTTGCCTTTGCCATCAg                                     >  1:222626/1‑100 (MQ=255)
                                                            ccTTATCCGGCCTACAAAATCGTGCAAATTCAATACATTGC‑‑AGAATT‑‑‑TTC‑‑GTAGGCCTGATAAGCGGAGCGCATCAGGCATTTTTGCCTTTGCCATCAGt                                    >  1:1432319/1‑100 (MQ=255)
                                                                   cGGCCTACAAAATCGTGCAAATTCAATACATTGC‑‑AGAATT‑‑‑TTC‑‑GTAGGCCTGATAAGCGGAGCGCATCAGGCATTTTTGCCTTTGCCATCAGTCTTGTAt                             <  2:71294/100‑1 (MQ=255)
                                                                   cGGCCTACAAAATCGTGCAAATTCAATACATTGC‑‑AGAATT‑‑‑TTC‑‑GTAGGCCTGATAAGCGGAGCGCATCAGGCATTTTTGCCTTTGCCATCAGTCTTGTAt                             <  2:52514/100‑1 (MQ=255)
                                                                        tACAAAATCGTGCAAATTCAATACATTGC‑‑AGAATT‑‑‑TTC‑‑GTAGGCCTGATAAGCGGAGCGCATCAGGCATTTTTGCCTTTGCCATCAGTCTTGTATGGCtt                        <  2:1984560/100‑1 (MQ=255)
                                                                        tACAAAATCGTGCAAATTCAATACATTGC‑‑AGAATT‑‑‑TTC‑‑GTAGGCCTGATAAGCGGAGCGCATCAGGCATTTTTGCCTTTGCCATCAGTCTTGTATGGCtt                        <  2:1415270/100‑1 (MQ=255)
                                                                          cAAAATCGTGCAAATTCAATACATTGC‑‑AGAATT‑‑‑TTC‑‑GTAGGCCTGATAAGCGGAGCGCATCAGGCATTTTTGCCTTTGCCATCAGTCTTGTATGGCtttt                      <  2:847382/100‑1 (MQ=255)
                                                                          cAAAATCGTGCAAATTCAATACATTGC‑‑AGAATT‑‑‑TTC‑‑GTAGGCCTGATAAGCGGAGCGCATCAGGCATTTTTGCCTTTGCCATCAGTCTTGTATGGCtttt                      <  2:1562128/100‑1 (MQ=255)
                                                                            aaaaTCGTGCAAATTCAATACATTGC‑‑AGAATT‑‑‑TTC‑‑GTAGGCCTGATAAGCGGAGCGCATCAGGCATTTTTGCCTTTGCCATCAGTCTTGTATGGCTTTTa                     <  2:1254906/100‑1 (MQ=255)
                                                                              aaTCGTGCAAATTCAATACATTGC‑‑AGAATT‑‑‑TTC‑‑GTAGGCCTGATAAGCGGAGCGCATCAGGCATTTTTGCCTTTGCCATCAGTCTTGTATGGCTTTTAat                   <  2:733000/100‑1 (MQ=255)
                                                                              aaTCGTGCAAATTCAATACATTGC‑‑AGAATT‑‑‑TTC‑‑GTAGGCCTGATAAGCGGAGCGCATCAGGCATTTTTGCCTTTGCCATCAGTCTTCTATGGCTTTTAat                   <  1:2072894/100‑1 (MQ=255)
                                                                                tCGTGCAAATTCGATACATTGC‑‑AGAATT‑‑‑TTC‑‑GTAGGCCTGGTAAGCGGAGCGCATCAGGCATTTTTGCCTTTGCCATCAGTCTTGTATGGCTTTTAatta                 <  1:456910/100‑1 (MQ=255)
                                                                                  gTGCAAATTCAATACATTGC‑‑AGAATT‑‑‑TTC‑‑GTAGGCCTGATAAGCGGAGCGCATCAGGCATTTTTGCCTTTGCCATCAGTCTTGTATGGCTTTTAattatt               >  2:339293/1‑100 (MQ=255)
                                                                                  gTGCAAATTCAATACATTGC‑‑AGAATT‑‑‑TTC‑‑GTAGGCCTGATAAGCGGAGCGCATCAGGCATTTTTGCCTTTGCCATCAGTCTTGTATGGCTTTTAattatt               >  1:2222297/1‑100 (MQ=255)
                                                                                    gCAAATTCAATACATTGC‑‑AGAATT‑‑‑TTC‑‑GTAGGCCTGATAAGCGGAGCGCATCAGGCATTTTTGCCTTTGCCATCAGTCTTGTATGGCTTTTAATTATTCg             <  1:40190/100‑1 (MQ=255)
                                                                                     cAAATTCAATACATTGC‑‑AGAATT‑‑‑TTC‑‑GTAGGCCTGATAAGCGGAGCGCATCAGGCATTTTTTTTTTTTTTTTCAGTCTTGTATGGCTTTTAattatt               >  2:742802/1‑97 (MQ=21)
                                                                                     cAAATTCAATACATTGC‑‑AGAATT‑‑‑TTC‑‑GTAGGCCTGATAAGCGGAGCGCATCAGGCATTTTTGCCTTTGCCATCAGTCTTGTATGGCTTTTAATTATTCGa            >  2:1581655/1‑100 (MQ=255)
                                                                                     cAAATTCAATACATTGC‑‑AGAATT‑‑‑TTC‑‑GTAGGCCTGATAAGCGGAGCGCATCAGGCATTTTTGCCTTTGCCATCAGTCTTGTATGGCTTTTAATTATTCGa            >  2:1231210/1‑100 (MQ=255)
                                                                                      aaaTTCAATACATTGC‑‑AGAATT‑‑‑TTC‑‑GTAGGCCTGATAAGCGGAGCGCATCAGGCATTTTTGCCTTTGCCATCAGTCTTGTATGGCTTTTAATTATTCGAt           <  1:178870/100‑1 (MQ=255)
                                                                                       aaTTCAATACATTGC‑‑AGAATT‑‑‑TTC‑‑GTAGGCCTGATAAGCGGAGCGCATCAGGCATTTTTGCCTTTGCCATCAGTCTTGTATGGCTTTTAATTATTCGATa          >  1:1653870/1‑100 (MQ=255)
                                                                                       aaTTCAATACATTGC‑‑AGAATT‑‑‑TTC‑‑GTAGGCCTGATAAGCGGAGCGCATCAGGCATTTTTGCCTTTGCCATCAGTCTTGTATGGCTTTTAATTATTCGATa          <  2:2571158/100‑1 (MQ=255)
                                                                                        aTTCAATACATTGC‑‑AGAATT‑‑‑TTC‑‑GTAGGCCTGATAAGCGGAGCGCATCAGGCATTTTTGCCTTTGCCATCAGTCTTGTATGGCTTTTAATTATTCGATag         >  2:1657059/1‑100 (MQ=255)
                                                                                           cAATACATTGC‑‑AGAATT‑‑‑TTC‑‑GTAGGCCTGATAAGCGGAGCGCATCAGGCATTTTTGCCTTTGCCATCAGTCTTGTATGGCTTTTAATTATTCGATAgaga      >  2:1902774/1‑100 (MQ=255)
                                                                                                |                                                                                                          
CACCGGGTATACAGATTGTTGACAAAGTGCGCTTTGTTTATGCCGGATGCGGCGTTAACGCCTTATCCGGCCTACAAAATCGTGCAAATTCAATACATTGC‑‑AGAATT‑‑‑TTC‑‑GTAGGCCTGATAAGCGGAGCGCATCAGGCATTTTTGCCTTTGCCATCAGTCTTGTATGGCTTTTAATTATTCGATAGAGACGTT  >  NC_000913/3705922‑3706115

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: