Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009133,705,9780AC24.3% 430.3 / 156.8 211intergenic (‑188/+120)dppB/dppAdipeptide/heme ABC transporter permease/dipeptide/heme ABC transporter periplasmic binding protein; dipeptide chemotaxis receptor
Reads supporting (aligned to +/- strand):  ref base A (64/82);  new base C (47/0);  total (128/83)
Fisher's exact test for biased strand distribution p-value = 4.85e-14
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 3.08e-01

Only 100 of 215 total aligned reads displayed.

AAGTATTGCTGCCGGATGGCGCTGCACGCAGAGGCACAGCGCCACCGGGTATACAGATTGTTGACAAAGTGCGCTTTGTTTATGCCGGATGCGGCGTTAACGCCTTATCCGGCCTACAAAATCGTGCAAATTCAATACATTGCAGAATT‑‑‑TTC‑‑GTAGGCCTGATAAGCGGAGCGCATCAGGCATTTTTGCCTTTGCCAT  >  NC_000913/3705880‑3706077
                                                                                                  |                                                                                                            
aaGTATTGCTGCCGGATGGCGCTGCACGCAGAGGCACAGCGCCACCGGGTATACAGATTGTTGACAAAGTGCGCTTTGTTTATGCCGGATGCGGCGTTaa                                                                                                             <  2:1908901/100‑1 (MQ=255)
aaGTATTGCTGCCGGATGGCGCTGCACGCAGAGGCACAGCGCCACCGGGTATACAGATTGTTGACAAAGTGCGCTTTGTTTATGCCGGATGCGGCGTTaa                                                                                                             >  2:2590650/1‑100 (MQ=255)
aaGTATTGCTGCCGGATGGCGCTGCACGCAGAGGCACAGCGCCACCGGGTATACAGATTGTTGACAAAGTGCGCTTTGTTTATGCCGGATGCGGCGTTaa                                                                                                             >  2:1853845/1‑100 (MQ=255)
 aGTATTGCTGCCGGATGGCGCTGCACGCAGAGGCACAGCGCCACCGGGTATACAGATTGTTGACAAAGTGCGCTTTGTTTATGCCGGATGCGGCGTTAAc                                                                                                            <  2:1235513/100‑1 (MQ=255)
  gTATTGCTGCCGGATGGCGCTGCACGCAGAGGCACAGCGCCACCGGGTATACAGATTGTTGACAAAGTGCGCTTTGTTTATGCCGGATGCGGCGTTAACg                                                                                                           >  2:624210/1‑100 (MQ=255)
   tATTGCTGCCGGATGGCGCTGCACGCAGAGGCACAGCGCCACCGGGTATACAGATTGTTGACAAAGTGCGCTTTGTTTATGCCGGATGCGGCGTTAACGc                                                                                                          <  2:2058736/100‑1 (MQ=255)
    aTTGCTGCCGGATGGCGCTGCACGCAGAGGCACAGCGCCACCGGGTATACAGATTGTTGACAAAGTGCGCTTTGTTTATGCCGGATGCGGCGTTAACGcc                                                                                                         >  2:862193/1‑100 (MQ=255)
       gctgcCGGATGGCGCTGCACGCAGAGGCACAGCGCCACCGGGTATACAGATTGTTGACAAAGTGCGCTTTGTTTATGCCGGATGCGGCGTTAACGCCTTa                                                                                                      >  1:2187820/1‑100 (MQ=255)
       gctgcCGGATGGCGCTGCACGCAGAGGCACAGCGCCACCGGGTATACAGATTGTTGACAAAGTGCGCTTTGTTTATGCCGGATGCGGCGTTAACGCCTTa                                                                                                      >  2:1877399/1‑100 (MQ=255)
          gcCGGATGGCGCTGCACGCAGAGGCACAGCGCCACCGGGTATACAGATTGTTGACAAAGTGCGCTTTGTTTATGCCGGATGCGGCGTTAACGCCTTATcc                                                                                                   >  1:2124354/1‑100 (MQ=255)
           ccGGATGGCGCTGCACGCAGAGGCACAGCGCCACCGGGTATACAGATTGTTGACAAAGTGCGCTTTGTTTATGCCGGATGCGGCGTTAACGCCTTATCCg                                                                                                  >  2:2534345/1‑100 (MQ=255)
            cGGATGGCGCTGCACGCAGAGGCACAGCGCCACCGGGTATACAGATTGTTGACAAAGTGCGCTTTGTTTATGCCGGATGCGGCGTTAACGCCTTATCCgg                                                                                                 <  2:2829031/100‑1 (MQ=255)
             ggATGGCGCTGCACGCAGAGGCACAGCGCCACCGGGTATACAGATTGTTGACAAAGTGCGCTTTGTTTATGCCGGATAAGGCGTTAACGCCGCATCCGGc                                                                                                <  2:301978/100‑1 (MQ=37)
              gATGGCGCTGCACGCAGAGGCACAGCGCCACCGGGTATACAGATTGTTGACAAAGTGCGCTTTGTTTATGCCGGATGCGGCGTTAACGCCTTATCCGGcc                                                                                               <  2:33646/100‑1 (MQ=255)
                  gcgcTGCACGCAGAGGCACAGCGCCACCGGGTATACAGATTGTTGACAAAGTGCGCTTTGTTTATGCCGGATGCGGCGTTAACGCCTTATCCGGCCTACa                                                                                          >  2:701845/1‑100 (MQ=255)
                  gcgcTGCACGCAGAGGCACAGCGCCACCGGGTATACAGATTGTTGACAAAGTGCGCTTTGTTTATGCCGGATGCGGCGTTAACGCCTTATCCGGCCTACa                                                                                          >  2:1131162/1‑100 (MQ=255)
                   cgcTGCACGCAGAGGCACAGCGCCACCGGGTATACAGATTGTTGACAAAGTGCGCTTTGTTTAGGCCGGATAAGGCGTTAACGCCGCATCCGGcataaac                                                                                          <  2:1203406/100‑7 (MQ=21)
                     cTGCACGCAGAGGCACAGCGCCACCGGGTATACAGATTGTTGACAAAGTGCGCTTTGTTTATGCCGGATGCGGCGTTAACGCCTTATCCGGCCTACaaaa                                                                                       <  1:126618/100‑1 (MQ=255)
                       gCACGCAGAGGCACAGCGCCACCGGGTATACAGATTGTTGACAAAGTGCGCTTTGTTTATGCCGGATGCGGCGTTAACGCCTTATCCGGCCTACAAAAt                                                                                      <  2:129548/99‑1 (MQ=255)
                       gCACGCAGAGGCACAGCGCCACCGGGTATACAGATTGTTGACAAAGTGCGCTTTGTTTATGCCGGATGCGGCGTTAACGCCTTATCCGGCCTACAAAATc                                                                                     >  2:945896/1‑100 (MQ=255)
                        cACGCAGAGGCACAGCGCCACCGGGTATACAGATTGTTGACAAAGTGCGCTTTGTTTATGCCGGATGCGGCGTTAACGCCTTATCCGGCCTACAAAATCg                                                                                    <  2:440234/100‑1 (MQ=255)
                        cACGCAGAGGCACAGCGCCACCGGGTATACAGATTGTTGACAAAGTGCGCTTTGTTTATGCCGGATGCGGCGTTAACGCCTTATCCGGCCTACAAAATCg                                                                                    <  2:2132898/100‑1 (MQ=255)
                        cACGCAGAGGCACAGCGCCACCGGGTATACAGATTGTTGACAAAGTGCGCTTTGTTTATGCCGGATGCGGCGTTAACGACTTATCCGGCCTACAAAATCg                                                                                    <  2:2414023/100‑1 (MQ=255)
                        cACGAAGAGGCACAGTGCCACCGTGTATACAGATTGTTGACAAAGTGCGCTTTGTTTATGCCGTATTTGGCGTTAACGCCTTATCCGGCCTACAAAATCg                                                                                    <  1:610574/100‑1 (MQ=18)
                         aCGCAGAGGCACAGCGCCACCGGGTATACAGATTGTTGACAAAGTGCGCTTTGTTTATGCCGGATGCGGCGTTAACGCCTTATCCGGCCTACAAAATCGt                                                                                   >  2:1716547/1‑100 (MQ=255)
                          cGCAGAGGCACAGCGCCACCGGGTATACAGATTGTTGACAAAGTGCGCTTTGTTTATGCCGGATGCGGCGTTAACGCCTTATCCGGCCTACAAAATCGTg                                                                                  <  2:691862/100‑1 (MQ=255)
                            cAGAGGCACAGCGCCACCGGGTATACAGATTGTTGACAAAGTGCGCTTTGTTTATGCCGGACGCGGCGTTAACGCCTTATCCGGCCTACAAAATCGTGCa                                                                                <  1:1853837/100‑1 (MQ=255)
                                  cacaGCGCCACCGGGTATACAGATTGTTGACAAAGTGCGCTTTGTTTATGCCGGATGCGGCGTTAACGCCTTATCCGGCCTACAAAATCGTGCAAATTCa                                                                          <  1:2534333/100‑1 (MQ=255)
                                  cacaGCGCCACCGGGTATACAGATTGTTGACAAAGTGCGCTTTGTTTATGCCGGATGCGGCGTTAACGCCTTATCCGGCCTACAAAATCGTGCAAATTCa                                                                          >  2:2749211/1‑100 (MQ=255)
                                   acaGCGCCACCGGGTATACAGATTGTTGACAAAGTGCGCTTTGTTTATGCCGGATGCGGCGTTAACGCCTTATCCGGCCTACAAAATCGTGCAAATTCaa                                                                         <  2:1760547/100‑1 (MQ=255)
                                    caGCGCCACCGGGTATACAGATTGTTGACAAAGTGCGCTTTGTTTATGCCGGATGCGGCGTTAACGCCTTATCCGGCCTACAAAATCGTGCAAATTCAAt                                                                        <  2:2941306/100‑1 (MQ=255)
                                     aGCGCCACCGGGTATACAGATTGTTGACAAAGTGCGCTTTGTTTATGCCGGATGCGGCGTTAACGCCTTATCCGGCCTACAAAATCGTGCAAATTCAATa                                                                       <  1:778510/100‑1 (MQ=255)
                                      gcgcCACCGGGTATACAGATTGTTGACAAAGTGCGCTTTGTTTATGCCGGATGCGGCGTTAACGCCTTATCCGGCCTACAAAATCGTGCAAATTCAATAc                                                                      <  1:647090/100‑1 (MQ=255)
                                         ccACCGGGTATACAGATTGTTGACAAAGTGCGCTTTGTTTATGCCGGATGCGGCGTTAACGCCTTATCCGGCCTACAAAATCGTGCAAATTCAATACAtt                                                                   <  2:354389/100‑1 (MQ=255)
                                           aCCGGGTATACAGATTGTTGACAAAGTGCGCTTTGTTTATGCCGGATGCGGCGTTAACGCCTTATCCGGCCTACAAAATCGTGCAAATTCAATACATTGc                                                                 <  1:1047383/100‑1 (MQ=255)
                                               ggTATACAGATTGTTGACAAAGTGCGCTTTGTTTATGCCGGATGCGGCGTTAACGCCTTATCCGGCCTACAAAATCGTGCAAATTCAATACATTGCAGaa                                                           >  2:1614804/1‑100 (MQ=255)
                                                   taCAGATTGTTGACAAAGTGCGCTTTGTTTATGCCGGATGCGGCGTTAACGCCTTATCCGGCCTACAAAATCGTGCAAATTCAATACATTGCAGAAtt‑‑‑tt                                                   >  2:2655802/1‑100 (MQ=255)
                                                   taCAGATTGTTGACAAAGTGCGCTTTGTTTATGCCGGATGCGGCGTTAACGCCTTATCCGGCCTACAAAATCGTGCAAATTCAATACATTGCAGAAtt‑‑‑tt                                                   >  2:806059/1‑100 (MQ=255)
                                                    aCAGATTGTTGACAAAGTGCGCTTTGTTTATGCCGGATGCGGCGTTAACGCCTTATCCGGCCTACAAAATCGTGCAAATTCAATACATTGCAGAATT‑‑‑TTc                                                  <  1:2313453/100‑1 (MQ=255)
                                                     cAGATTGTTGACAAAGTGCGCTTTGTTTATGCCGGATGCGGCGTTAACGCCTTATCCGGCCTACAAAATCGTGCAAATTCAATACATTGCAGAATT‑‑‑TTC‑‑g                                               <  2:1088951/100‑1 (MQ=255)
                                                         ttgttgACAAAGTGCGCTTTGTTTATGCCGGATGCGGCGTTAACGCCTTATCCGGCCTACAAAATCGTGCAAATTCAATACATTGCAGAATT‑‑‑TTC‑‑GTAgg                                           <  1:701843/100‑1 (MQ=255)
                                                          tgttgACAAAGTGCGCTTTGTTTATGCCGGATGCGGCGTTAACGCCTTATCCGGCCTACAAAATCGTGCAAATTCAATACATTGCAGAATT‑‑‑TTC‑‑GTAGGc                                          >  2:1644753/1‑100 (MQ=255)
                                                           gttgACAAAGTGCGCTTTGTTTATGCCGGATGCGGCGTTAACGCCTTATCCGGCCTACAAAATCGTGCAAATTCAATACATTGCAGAATT‑‑‑TTC‑‑GTAGGcc                                         <  2:571258/100‑1 (MQ=255)
                                                           gttgACAAAGTGCGCTTTGTTTATGCCGGATGCGGCGTTAACGCCTTATCCGGCCTACAAAATCGTGCAAATTCAATACATTGCAGAATT‑‑‑TTC‑‑GTAGGcc                                         <  2:1056254/100‑1 (MQ=255)
                                                            ttgACAAAGTGCGCTTTGTTTATGCCGGATGCGGCGTTAACGCCTTATCCGGCCTACAAAATCGTGCAAATTCAATACATTGCAGAATT‑‑‑TTC‑‑GTAGGCCt                                        <  1:2063062/100‑1 (MQ=255)
                                                             tgACAAAGTGCGCTTTGTTTATGCCGGATGCGGCGTTAACGCCTTATCCGGCCTACAAAATCGTGCAAATTCAATACATTGCAGAATT‑‑‑TTC‑‑GTAGGCCTg                                       >  1:419950/1‑100 (MQ=255)
                                                             tgACAAAGTGCGCTTTGTTTATGCCGGATGCGGCGTTAACGCCTTATCCGGCCTACAAAATCGTGCAAATTCAATACATTGCAGAATT‑‑‑TTC‑‑GTAGGCCTg                                       >  1:896462/1‑100 (MQ=255)
                                                             tgACAAAGTGCGCTTTGTTTATGCCGGATGCGGCGTTAACGCCTTATCCGGCCTACAAAATCGTGCAAATTCAATACATTGCAGAATT‑‑‑TTC‑‑GTAGGCCTg                                       <  1:2391242/100‑1 (MQ=255)
                                                             ngACAAAGTGCGCTTTGTTTATGCCGGATGCGGCGTTAACGCCTTATCCGGCCTACAAAATCGTGCAAATTCAATACATTGCAGAATT‑‑‑TTC‑‑GTAGGCCTg                                       >  1:1416603/2‑100 (MQ=255)
                                                              gACAAAGTGCGCTTTGTTTATGCCGGATGCGGCGTTAACGCCTTATCCGGCCTACAAACTCGTGCAAATTCAATACATTGCAGAATT‑‑‑TTC‑‑GTAGGCCTGa                                      <  1:1644745/100‑1 (MQ=255)
                                                                cAAAGTGCGCTTTGTTTATGCCGGATGCGGCGTTAACGCCTTATCCGGCCTACAAAATCGTGCAAATTCAATACATTGCAGAATT‑‑‑TTC‑‑GTAGGCCTGATa                                    >  1:11151/1‑100 (MQ=255)
                                                                    gTGCGCTTTGTTTATGCCGGATGCGGCGTTAACGCCTTATCCGGCCTACAAAATCGTGCAAATTCAATACATTGCAGAATT‑‑‑TTC‑‑GTAGGCCTGATAAGCg                                >  2:1733813/1‑100 (MQ=255)
                                                                     tGCGCTTTGTTTATGCCGGATGCAGCATTAACGCCTTATCCGGCATACAAAATCGTGCAAATTCGAGAAATTGCAGAATT‑‑‑TTC‑‑GTAGGCCTGATAAGCgg                               <  1:1783048/100‑1 (MQ=17)
                                                                      gcgcTTTGTTTATGCCGGATGCGGCGTTAACGCCTTATCCGGCCTACAAAATCGTGCAAATTCAATACATTGCAGAATT‑‑‑TTC‑‑GTAGGCCTGATAAGCGGa                              >  2:1422895/1‑100 (MQ=255)
                                                                       ctcTTTGTTTATGCCGGATGCGGCGTTAACGCCTTATCCGGCCTACAAAATCGTGCAAATTCAATACATTGCAGAATT‑‑‑TTC‑‑GTAGGCCTGATAAGCGGAg                             <  1:972095/98‑1 (MQ=255)
                                                                        gcTTTGTTTATGCCGGATGCGGCGTTAACGCCTTATCCGGCCTACAAAATCGTGCAAATTCAATACATTGCAGAATT‑‑‑TTC‑‑GTAGGCCTGATAAGCGGAgc                            <  1:2398827/100‑1 (MQ=255)
                                                                         cTTTGTTTATGCCGGATGCGGCGTTAACGCCTTATCCGGCCTACAAAATCGTGCAAATTCAATACATTGCAGAATT‑‑‑TTC‑‑GTAGGCCTGATAAGCGGAgcg                           <  1:1821906/100‑1 (MQ=255)
                                                                             ttgtaggCCGGATAAGGCGTTTACGCCGCATCCGGCCTACAAAAATGTGCAAATTCAATAAATTGCAATTCAACTT‑‑GTAGGCCTGATAAGCGCAGCGCAt                        >  1:1993749/7‑100 (MQ=21)
                                                                               ttATGCCGGATGCGGCGTTAACGCCTTATCCGGCCTACAAAATCGTGCAAATTCAATACATTGCAGAATT‑‑‑TTC‑‑GTAGGCCTGATAAGCGGAGCGCATCAg                     <  1:2478080/100‑1 (MQ=255)
                                                                               gtaggCCGGATAAGGCGTTCACGCCGCATCCGGCCTACAAAACTGTGCAAATTCAATATATTGCAGGAAA‑‑‑CAC‑‑GTAGGCCTGATAAGCGAAGCGCATCAg                     >  2:1450379/5‑100 (MQ=18)
                                                                               gtaggCCGGATAAGGCGTTCACGCCGCATCCGGCCTACAAACCGCGCAAATTCAATATATTGCGGAGAA‑‑‑AAT‑‑GTAGGCCTGATAAGCGTAGCGCATCAgg                    >  1:2513595/5‑100 (MQ=18)
                                                                                taggCCGGATAAGGCGTTCACGCCGCATCCGGCCTACAAAAATCTTGCCAATTCAATATATTGCAGGACC‑‑‑AT‑‑GTAGGCCTGATAAGCGTAGCGCATCAgg                    >  2:773755/4‑100 (MQ=18)
                                                                                taggCCGGATAAGGCGTTCACGCCGCATCCGGCCTACAAAATTGTGTAAATTCAATAAATTGATCGATC‑‑‑CAC‑‑GTAGGCCTGATAAGCGTAGCGCATCAgg                    >  2:373817/4‑100 (MQ=21)
                                                                                taggCCGGATAAGGCGTTCACGCCGCATCCGGCCTACAAAACCTTGCTAATTCAATATATTGCAGGGACT‑‑‑AT‑‑GTAGGCCTGATAAGCATAGCGCATCAgg                    >  1:2565217/4‑100 (MQ=18)
                                                                                taggCCGGATAAGGCGTTCACGCCGCATCCGGCCTACAAAACCTTGCTAATTCAATATATTGCAGGGACT‑‑‑AT‑‑GTAGGCCTGATAAGCATAGCGCATCAgg                    >  2:2350227/4‑100 (MQ=18)
                                                                                 aggCCGGATAAGGCGTTCACGCCGCATCCGGCCTACAAAAGCACGCAAACTCAATATATTGCAGAGATC‑‑‑AT‑‑GTAGGCCTGATAAGCGTAGCGCATCAGGc                   >  2:2804013/3‑100 (MQ=17)
                                                                                 aggCCGGATAAGGCGTTAACGCCGCATCCGGCCTACAAAATCGTGCAAATTCAATACATTGCAGAATT‑‑‑TTC‑‑GTAGGCCTGATAAGCGGAGCGCATCAGGc                   >  2:2463611/3‑100 (MQ=18)
                                                                                 aggCCGGATAAGGCGTTAACGCCGCATCCGGCCTACAAAATCGTGCAAATTCAATACATTGCAGAATT‑‑‑TTC‑‑GTAGGCCTGATAAGCGGAGCACATCAGGc                   >  1:1729995/3‑100 (MQ=18)
                                                                                 aTGCCGGATGCGGCGTTAACGCCTTATCCGGCCTACAAAATCGTGCAAATTCAATACATTGCAGAATT‑‑‑TTC‑‑GTAGGCCTGATAAGCGGAGCGCATCAg                     <  1:2073947/98‑1 (MQ=255)
                                                                                  ngCCGGATAAGGCGTTCACGCCGCATCCGGCCTACAAAACTGTGCAAATTCAATATATTGCAGGAAA‑‑‑CAC‑‑GTAGGCCTGATAAGCGAAGCGCATCAGGCa                  >  2:2645807/2‑100 (MQ=18)
                                                                                  ggCCGGATAAGGCGTTTACGCCGCATCCGGCCTACAAAAGCGTGCAAATTCAATACATTGCATGGGCC‑‑‑AT‑‑GTAGGCCTGATAAGCGTAGCGCATCAGGCa                  >  2:957053/2‑100 (MQ=17)
                                                                                  ggCCGGATAAGGCGTTTACGCCGCATCCGGCCTACAAAAGCACGCAAACTCAATATATTGCAGAGAT‑‑‑CAT‑‑GTAGGCCTGATAAGCGTAGCGCATCAGGCa                  >  1:2921783/2‑100 (MQ=17)
                                                                                  ggCCGGATAAGGCGTTCACGCCGCATCCGGCCTACAAAGTCGTATAAATTCAATGTATTGCAGTTTA‑‑‑TTG‑‑GTAGGCCTGATAAGCGCAGCGCATCAGGCa                  >  1:2889630/2‑100 (MQ=21)
                                                                                  ggCCGGATAAGGCGTTCACGCCGCATCCGGCCTACAAAATCGCTTAAATTCAATATATTGCACCATC‑‑‑CTC‑‑GTAGGCCTGATAAGCGTAGCGCATCAGGCa                  >  1:366597/2‑100 (MQ=21)
                                                                                  ggCCGGATAAGGCGTTCACGCCGCATCCGGCCTACAAAACTTTGCAAATTCAATATATTGCAATC‑‑‑TCCGTGTAGGCCTGATAAGCGTAGCGCATCAGGCa                  >  2:576834/2‑100 (MQ=18)
                                                                                   gCCGGATGCGGCGTTAACGCCTTATCCGGCCTACAAAATCGTGCAAATTCAATACATTGGCGAATT‑‑‑TTC‑‑GTAGGCCTGATCAGCAGAGCGCATCAAGcct                 >  1:283062/1‑98 (MQ=17)
                                                                                   gCCGGATAAGGCGTTTACGCCGCATCCGGCCTACAAAACAGTGCAAATTCAATATATTGCAATCTT‑‑‑CTG‑‑GTAGGCCTGATAAGCATCGCGCATCAGGCaa                 >  1:891941/1‑99 (MQ=18)
                                                                                   gCCGGATAAGGCGTTCACGCCGCATCCGGCCTACAAAAATCTTGCCAATTCAATATATTGCAGGACC‑‑‑AT‑‑GTAGGCCTGATAAGCGTAGCGCATCAGGCag                 >  2:936288/1‑99 (MQ=18)
                                                                                   gCCGGATAAGGCGTTCACGCCGCATCCGGCCTACAAAGTCGTATAACTTCAATGTATTGCAGTTTA‑‑‑TTG‑‑GTAGGCCTGATAAGCGCAGCGCATCAGGCaa                 >  1:1369201/1‑99 (MQ=18)
                                                                                   gCCGGATAAGGCGTTCACGCCGCATCCGGCCTACAAAAACCTGCAAATTCAATATATTGCAGGATA‑‑‑CGT‑‑GTAGGCCTGATAAGCGTAGCGCATCAGGCaa                 >  2:2385627/1‑99 (MQ=18)
                                                                                   gCCGGATAAGGCGTTCACGCCGCATCCGGCCTACAAAAACATGCAAATTCAATAGATTGCAGAGATT‑‑‑AT‑‑GTAGGCCTGATAAGCGTAGCGCATCAGGCaa                 >  1:2236509/1‑99 (MQ=17)
                                                                                   gCCGGATAAGGCGTTCACGCCGCATCCGGCCTACAAATCATGCAAATCCAAAATATTGCAGGAGT‑‑‑TAT‑‑GTAGGCCTGATAAGCGTAGCGCATCAGGCAtt                >  1:44038/1‑100 (MQ=21)
                                                                                    ccggataaggcgTTTACGCCGCATCCGGCCTACAAATTCTGCAAATTCAACAAATTGCGAATCC‑‑‑CTT‑‑GTAGGCCTGATAAGCGTAGCGCATCAGGCAttt               >  2:1800455/9‑100 (MQ=18)
                                                                                    ccggataaggcgTTCACGCCGCATCCGGCTTACAAAACTTTGCAAATTCAATATATTGCAATCTC‑‑‑CGT‑‑GTAGGCCTGATAAGCGTAGCGCATCAGGCaat                >  2:2391848/9‑98 (MQ=18)
                                                                                    ccggataaggcgTTCACGCCGCATCCGGCCTACAAAATCGCTTAAATTCAATATATTGCACCATC‑‑‑CTC‑‑GTAGGCCTGATAAGCGTAGCGCATCAGGCaat                >  2:601831/9‑98 (MQ=21)
                                                                                    ccggataaggcgTTCACGCCGCATCCGGCCTACAAAATCATGAAAATTCAATACATTGCAAGATT‑‑‑TTC‑‑GTAGGCCTGATAAGCGTAGCGCATCAGGCacg                >  2:2012178/9‑98 (MQ=17)
                                                                                    ccggataaggcgTTCACGCCGCATCCGGCCTACAAAAACATGCAAATTCAATAGATTGCAGAGATT‑‑‑AT‑‑GTAGGCCTGATAAGCGTAGCGCATCAGGCaat                >  2:719438/9‑98 (MQ=17)
                                                                                     cggataaggcgTTTACGCCGCATCCGGCCTACAAAAGCGTGCAAATTCAATACATTGCATGGGCC‑‑‑AT‑‑GTAGGCCTGATAAGCGTAGCGCATCAGGCaatt               >  2:2476859/8‑97 (MQ=18)
                                                                                     cggataaggcgTTCACGCCTTATCCGACCCACAAAGTCTTGCAAAATCAATAAATTGCAGGAAC‑‑‑GAT‑‑GTAGGCCTGATAAGCGTAGCGCATCAGGCaact               >  1:1586176/8‑97 (MQ=21)
                                                                                      ggataaggcgTTTACGCCGCATCCGGCCTACAAAACAGTGCAAATTCAATATATTGCAATCTT‑‑‑CTG‑‑GTAGGCCTGATAAGCATCGCGCATCAGGCaagct              >  1:730682/7‑96 (MQ=18)
                                                                                      ggataaggcgTTCACGCCGCATCCGGCCTACAAAAATCTTGCCAATTCAATATATTGCAGGACC‑‑‑AT‑‑GTAGGCCTGATAAGCGTAGCGCACCAGGCAgttt              >  2:1600501/7‑100 (MQ=18)
                                                                                      ggataaggcgTTCACGCCGCATCCGGCCTACGAAATCGTGCTAATTCAATATATTGTGGGCTC‑‑‑CAC‑‑GTAGGCCTGATAAGCGTAGCGCATCAGGCAgttt              >  1:1695713/7‑100 (MQ=21)
                                                                                      ggataaggcgTTCACGCCGCATCCGGCCTACAAAATCGCTTAAATTCAATATATTGCACCATC‑‑‑CTC‑‑GTAGACCTGATAAGCGTAGCGCATCAGGCAattc              >  1:1453474/7‑99 (MQ=18)
                                                                                      ggataaggcgTTCACGCCGCATCCGGCCTACAAAATCATGAAAATTCAATACATTGCAAGATT‑‑‑TTC‑‑GTAGGCCTGATAAGCGTAGCGCATCAGGCacgct              >  1:1181880/7‑96 (MQ=17)
                                                                                      ggataaggcgTTCACGCCGCATCCGGCCTACAAAACTGTGCAAATTCAATATATTGCAGGAAA‑‑‑CAC‑‑GTAGGCCTGATAAGCGAAGCGCATCAGGCAgttt              >  1:338876/7‑100 (MQ=18)
                                                                                      ggataaggcgTTAACGCCGCATCCGGCCTACAAAATCGTGCAAATTCAATACATTGCAGAATT‑‑‑TTC‑‑GTAGGCCTGATAAGCGGAGCGCATCAGGCAtttt              >  2:2082386/7‑100 (MQ=18)
                                                                                           cggcgTTAACGCCTTATCCGGCCTACAAAATCGTGCAAATTCAATACATTGCAGAATT‑‑‑TTC‑‑GTAGGCCTGATAAGCGGAGCGCATCAGGCATTtttgcct         >  2:2869997/1‑100 (MQ=255)
                                                                                            ggcgTTAACGCCTTATCCGGCCTACAAAATCGTGCAAATTCAATACATTGCAGAATT‑‑‑TTC‑‑GTAGGCCTGATAAGCGGAGCGCATCAGGCATTtttgcctt        >  2:2099499/1‑100 (MQ=255)
                                                                                             gcgTTAACGCCTTATCCGGCCTACAAAATCGTGCAAATTCAATACATTGCAGAATT‑‑‑TTC‑‑GTAGGCCTGATAAGCGGAGCGCATCAGGCATTtttgccttt       <  1:672352/100‑1 (MQ=255)
                                                                                                  aaCGCCTTATCCGGCCTACAAAATCGTGCAAATTCAATACATTGCAGAATT‑‑‑TTC‑‑GTAGGCCTGATAAGCGGAGCGCATCAGGCATTTTTGCCTTTGCCAt  >  2:1278732/1‑100 (MQ=255)
                                                                                                  |                                                                                                            
AAGTATTGCTGCCGGATGGCGCTGCACGCAGAGGCACAGCGCCACCGGGTATACAGATTGTTGACAAAGTGCGCTTTGTTTATGCCGGATGCGGCGTTAACGCCTTATCCGGCCTACAAAATCGTGCAAATTCAATACATTGCAGAATT‑‑‑TTC‑‑GTAGGCCTGATAAGCGGAGCGCATCAGGCATTTTTGCCTTTGCCAT  >  NC_000913/3705880‑3706077

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