Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009133,705,9700GA26.8% 351.4 / 129.4 180intergenic (‑180/+128)dppB/dppAdipeptide/heme ABC transporter permease/dipeptide/heme ABC transporter periplasmic binding protein; dipeptide chemotaxis receptor
Reads supporting (aligned to +/- strand):  ref base G (57/74);  new base A (43/5);  total (100/80)
Fisher's exact test for biased strand distribution p-value = 9.91e-09
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 4.72e-01

Only 100 of 195 total aligned reads displayed.

CAGGGAATGAAGTATTGCTGCCGGATGGCGCTGCACGCAGAGGCACAGCGCCACCGGGTATACAGATTGTTGACAAAGTGCGCTTTGTTTATGCCGGATGCGGCGTTAACGCCTTATCCGGCCTACAAAATCGTGCAAATTCAATACATTGCAGAATTTTCGTAGGCCTGATAAGCGGAGCGCATCAGGCATTTTTGC  >  NC_000913/3705871‑3706068
                                                                                                   |                                                                                                     
cAGGGAATGAAGTATTGCTGCCGGATGGCGCTGCACGCAGAGGCACAGCGCCACCGGGTATACAGATTGTTGACAAAGTGCGCTTTGTTTATGCCGGATg                                                                                                       >  2:2629751/1‑100 (MQ=255)
cAGGGAATGAAGTATTGCTGCCGGATGGCGCTGCACGCAGAGGCACAGCGCCACCGGGTATACAGATTGTTGACAAAGTGCGCTTTGTTTATGCCGGATg                                                                                                       <  1:527644/100‑1 (MQ=255)
  gggAATGAAGTATTGCTGCCGGATGGCGCTGCACGCAGAGGCACAGCGCCACCGGGTATACAGATTGTTGACAAAGTGCGCTTTGTTTATGCCGGATgcg                                                                                                     <  2:26543/100‑1 (MQ=255)
   ggAATGAAGTATTGCTGCCGGATGGCGCTGCACGCAGAGGCACAGCGCCACCGGGTATACAGATTGTTGACAAAGTGCGCTTTGTTTATGCCGGATgcgg                                                                                                    >  1:2828757/1‑100 (MQ=255)
    gAATGAAGTATTGCTGCCGGATGGCGCTGCACGCAGAGGCACAGCGCCACCGGGTATACAGATTGTTGACAAAGTGCGCTTTGTTTATGCCGGATgcggc                                                                                                   <  1:886222/100‑1 (MQ=255)
     aaTGAAGTATTGCTGCCGGATGGCGCTGCACGCAGAGGCACAGCGCCACCGGGTATACAGATTGTTGACAAAGTGCGCTTTGTTTATGCCGGATgcggcg                                                                                                  <  2:901088/100‑1 (MQ=255)
     aaTGAAGTATTGCTGCCGGATGGCGCTGCACGCAGAGGCACAGCGCCACCGGGTATACAGATTGTTGACAAAGTGCGCTTTGTTTATGCCGGATgcggcg                                                                                                  <  2:866497/100‑1 (MQ=255)
       tGAAGTATTGCTGCCGGATGGCGCTGCACGCAGAGGCACAGCGCCACCGGGTATACAGATTGTTGACAAAGTGCGCTTTGTTTATGCCGGATGCGGCGtt                                                                                                >  1:1130737/1‑100 (MQ=255)
         aaGTATTGCTGCCGGATGGCGCTGCACGCAGAGGCACAGCGCCACCGGGTATACAGATTGTTGACAAAGTGCGCTTTGTTTATGCCGGATGCGGCGTTaa                                                                                              <  2:1908901/100‑1 (MQ=255)
         aaGTATTGCTGCCGGATGGCGCTGCACGCAGAGGCACAGCGCCACCGGGTATACAGATTGTTGACAAAGTGCGCTTTGTTTATGCCGGATGCGGCGTTaa                                                                                              >  2:1853845/1‑100 (MQ=255)
          aGTATTGCTGCCGGATGGCGCTGCACGCAGAGGCACAGCGCCACCGGGTATACAGATTGTTGACAAAGTGCGCTTTGTTTATGCCGGATGCGGCGTTAAc                                                                                             >  2:2391253/1‑100 (MQ=255)
          aGTATTGCTGCCGGATGGCGCTGCACGCAGAGGCACAGCGCCACCGGGTATACAGATTGTTGACAAAGTGCGCTTTGTTTATGCCGGATGCGGCGTTAAc                                                                                             >  1:469123/1‑100 (MQ=255)
          aGTATTGCTGCCGGATGGCGCTGCACGCAGAGGCACAGCGCCACCGGGTATACAGATTGTTGACAAAGTGCGCTTTGTTTATGCCGGATGCGGCGTTAAc                                                                                             <  2:1235513/100‑1 (MQ=255)
           gTATTGCTGCCGGATGGCGCTGCACGCAGAGGCACAGCGCCACCGGGTATACAGATTGTTGACAAAGTGCGCTTTGTTTATGCCGGATGCGGCGTTAACg                                                                                            >  2:624210/1‑100 (MQ=255)
             aTTGCTGCCGGATGGCGCTGCACGCAGAGGCACAGCGCCACCGGGTATACAGATTGTTGACAAAGTGCGCTTTGTTTATGCCGGATGCGGCGTTAACGcc                                                                                          <  2:579148/100‑1 (MQ=255)
             aTTGCTGCCGGATGGCGCTGCACGCAGAGGCACAGCGCCACCGGGTATACAGATTGTTGACAAAGTGCGCTTTGTTTATGCCGGATGCGGCGTTAACGcc                                                                                          >  2:862193/1‑100 (MQ=255)
             aTTGCTGCCGGATGGCGCTGCACGCAGAGGCACAGCGCCACCGGGTATACAGATTGTTGACAAAGTGCGCTTTGTTTATGCCGGATGCGGCGTTAACGcc                                                                                          <  1:1198050/100‑1 (MQ=255)
                gctgcCGGATGGCGCTGCACGCAGAGGCACAGCGCCACCGGGTATACAGATTGTTGACAAAGTGCGCTTTGTTTATGCCGGATGCGGCGTTAACGCCTTa                                                                                       >  2:1877399/1‑100 (MQ=255)
                gctgcCGGATGGCGCTGCACGCAGAGGCACAGCGCCACCGGGTATACAGATTGTTGACAAAGTGCGCTTTGTTTATGCCGGATGCGGCGTTAACGCCTTa                                                                                       >  1:2187820/1‑100 (MQ=255)
                  tgcCGGATGGCGCTGCACGCAGAGGCACAGCGCCACCGGGTATACAGATTGTTGACAAAGTGCGCTTTGTTTATGCCGGATGCGGCGTTAACGCCTTATc                                                                                     <  2:2845829/100‑1 (MQ=255)
                    ccGGATGGCGCTGCACGCAGAGGCACAGCGCCACCGGGTATACAGATTGTTGACAAAGTGCGCTTTGTTTATGCCGGATGCGGCGTTAACGCCTTATCCg                                                                                   >  1:208606/1‑100 (MQ=255)
                     cGGATGGCGCTGCACGCAGAGGCACAGCGCCACCGGGTATACAGATTGTTGACAAAGTGCGCTTTGTTTATGCCGGATGCGGCGTTAACGCCTTATCCgg                                                                                  <  2:2828773/100‑1 (MQ=255)
                      ggATGGCGCTGCACGCAGAGGCACAGCGCCACCGGGTATACAGATTGTTGACAAAGTGCGCTTTGTTTATGCCGGATAAGGCGTTAACGCCGCATCCGGc                                                                                 <  2:301978/100‑1 (MQ=37)
                       gATGGCGCTGCACGCAGAGGCACAGCGCCACCGGGTATACAGATTGTTGACAAAGTGCGCTTTGTTTATGCCGGATAAGGCGTTAACGCCGCATCCGGca                                                                                <  2:788474/100‑2 (MQ=37)
                        ccgatcttGCACGCAGAGGCACAGCGCCACCGGGTATACAGATTGTTGACAAAGTGCGCTTTGTTTATGCCGGATGCGGCGTTAACGCCTTATCCGGCCt                                                                               <  2:959894/93‑1 (MQ=255)
                        aTGGCGCTGCACGCAGAGGCACAGCGCCACCGGGTATACAGATTGTTGACAAAGTGCGCTTTGTTTATGCCGGATAAGGCGTTAACGCCGCATCCGGcat                                                                               <  1:893750/100‑3 (MQ=37)
                           gcgcTGCACGCAGAGGCACAGCGCCACCGGGTATACAGATTGTTGACAAAGTGCGCTTTGTTTATGCCGGATGCGGCGTTAACGCCTTATCCGGCCTACa                                                                           >  2:1131162/1‑100 (MQ=255)
                           gcgcTGCACGCAGAGGCACAGCGCCACCGGGTATACAGATTGTTGACAAAGTGCGCTTTGTTTATGCCCGATGCGGCGTTAACGCCTTATCCGGCCTaaa                                                                            <  2:2017278/100‑3 (MQ=255)
                              cTGCACGCAGAGGCACAGCGCCACCGGGTATACAGATTGTTGACAAAGTGCGCTTTGTTTATGCCGGATGCGGCGTTAACGCCTTATCCGGCCTACaaaa                                                                        <  1:126618/100‑1 (MQ=255)
                               tGCACGCAGAGGCACAGCGCCACCGGGTATACAGATTGTTGACAAAGTGCGCTTTGTTTATGCCGGATGCGGCGTTAACGCCTTATCCGGCCTagatcgg                                                                        >  1:959890/1‑94 (MQ=255)
                                gCACGCAGAGGCACAGCGCCACCGGGTATACAGATTGTTGCCAAAGTGCGCTTTGTTTATGCCGGATGCGGCGTTAACGCCTTATCCGGCCTACAAAATc                                                                      >  1:781395/1‑100 (MQ=255)
                                 cACGCAGAGGCACAGCGCCACCGGGTATACAGATTGTTGACAAAGTGCGCTTTGTTTATGCCGGATGCGGCGTTAACGCCTTATCCGGCCTACAAAATCg                                                                     <  2:2132898/100‑1 (MQ=255)
                                 cACGCAGAGGCACAGCGCCACCGGGTATACAGATTGTTGACAAAGTGCGCTTTGTTTATGCCGGATGCGGCGTTAACGACTTATCCGGCCTACAAAATCg                                                                     <  2:2414023/100‑1 (MQ=255)
                                   cGCAGAGGCACAGCGCCACCGGGTATACAGATTGTTGACAAAGTGCGCTTTGTTTATGCCGGATGCGGCGTTAACGCCTTATCCGGCCTACAAAATCGTg                                                                   <  1:2629738/100‑1 (MQ=255)
                                    gCAGAGGCACAGCGCCACCGGGTATACAGATTGTTGACAAAGTGCGCTTTGTTTATGCCGGATGCGGCGTTAACGCCTTATCCGGCCTACAAAATCGTGc                                                                  >  1:1770296/1‑100 (MQ=255)
                                     cAGAGGCACAGCGCCACCGGGTATACAGATTGTTGACAAAGTGCGCTTTGTTTATGCCGGATGCGGCGTTAACGCCTTATCCGGCCTACAAAATCGTGCa                                                                 >  1:2014835/1‑100 (MQ=255)
                                     cAGAGGCACAGCGCCACCGGGTATACAGATTGTTGACAAAGTGCGCTTTGTTTATGCCGGACGCGGCGTTAACGCCTTATCCGGCCTACAAAATCGTGCa                                                                 <  1:1853837/100‑1 (MQ=255)
                                        agGCACAGCGCCACCGGGTATACAGATTGTTGACAAAGTGCGCTTTGTTTATGCCGGATGCGGCGTTAACGCCTTATCCGGCCTACAAAATCGTGCAAAt                                                             <  2:205921/100‑1 (MQ=255)
                                         ggCACAGCGCCACCGGGTATACAGATTGTTGACAAAGTGCGCTTTGTTTATGCCGGATGCGGCGTTAACGCCTTATCCGGCCTACAAAATCGTGCAAAtt                                                            >  1:1504437/1‑100 (MQ=255)
                                           cacaGCGCCACCGGGTATACAGATTGTTGACAAAGTGCGCTTTGTTTATGCCGGATGCGGCGTTAACGCCTTATCCGGCCTACAAAATCGTGCAAATTCa                                                          <  1:2534333/100‑1 (MQ=255)
                                           cacaGCGCCACCGGGTATACAGATTGTTGACAAAGTGCGCTTTGTTTATGCCGGATGCGGCGTTAACGCCTTATCCGGCCTACAAAATCGTGCAAATTCa                                                          >  2:2796692/1‑100 (MQ=255)
                                             caGCGCCACCGGGTATACAGATTGTTGACAAAGTGCGCTTTGTTTATGCCGGATGCGGCGTTAACGCCTTATCCGGCCTACAAAATCGTGCAAATTCAAt                                                        <  2:2941306/100‑1 (MQ=255)
                                              aGCGCCACCGGGTATACAGATTGTTGACAAAGTGCGCTTTGTTTATGCCGGATGCGGCGTTAACGCCTTATCCGGCCTACAAAATCGTGCAAATTCAATa                                                       <  1:1077300/100‑1 (MQ=255)
                                              aGCGCCACCGGGTATACAGATTGTTGACAAAGTGCGCTTTGTTTATGCCGGATGCGGCGTTAACGCCTTATCCGGCCTACAAAATCGTGCAAATTCAATa                                                       <  1:778510/100‑1 (MQ=255)
                                              aGCGCCACCGGGTATACAGATTGTTGACAAAGTGCGCTTTGTTTATGCCGGATGCGGCGTTAACGCCTTATCCGGCCTACAAAATCGTGCAAATTCAATa                                                       <  1:363991/100‑1 (MQ=255)
                                               gcgcCACCGGGTATACAGATTGTTGACAAAGTGCGCTTTGTTTATGCCGGATGCGGCGTTAACGCCTTATCCGGCCTACAAAATCGTGCAAATTCAATAc                                                      <  1:647090/100‑1 (MQ=255)
                                               gcgcCACCGGGTATACAGATTGTTGACAAAGTGCGCTTTGTTTATGCCGGATGCGGCGTTAACGCCTTATCCGGCCTACAAAATCGTGCAAATTCAATAc                                                      >  2:302547/1‑100 (MQ=255)
                                                  ccACCGGGTATACAGATTGTTGACAAAGTGCGCTTTGTTTATGCCGGATGCGGCGTTAACGCCTTATCCGGCCTACAAAATCGTGCAAATTCAATACAtt                                                   <  2:354389/100‑1 (MQ=255)
                                                   cACCGGGTATACAGATTGTTGACAAAGTGCGCTTTGTTTATGCCGGATGCGGCGTTAACGCCTTATCCGGCCTACAAAATCGTGCAAATTCAATACATTg                                                  >  2:122053/1‑100 (MQ=255)
                                                    aCCGGGTATACAGATTGTTGACAAAGTGCGCTTTGTTTATGCCGGATGCGGCGTTAACGCCTTATCCGGCCTACAAAATCGTGCAAATTCAATACATTGc                                                 <  1:1047383/100‑1 (MQ=255)
                                                     ccGGGTATACAGATTGTTGACAAAGTGCGCTTTGTTTATGCCGGATGCGGCGTTAACGCCTTATCCGGCCTACAAAATCGTGCAAATTCAATACATTGCa                                                <  2:821528/100‑1 (MQ=255)
                                                          tataCAGATTGTTGACAAAGTGCGCTTTGTTTATGCCGGATGCGGCGTTAACGCCTTATCCGGCCTACAAAATCGTGCAAATTCAATACATTGCAGAAtt                                          >  2:1294297/1‑100 (MQ=255)
                                                             aCAGATTGTTGACAAAGTGCGCTTTGTTTATGCCGGATGCGGCGTTAACGCCTTATCCGGCCTACAAAATCGTGCAAATTCAATACATTGCAGAATTTTc                                       <  1:1697648/100‑1 (MQ=255)
                                                             aCAGATTGTTGACAAAGTGCGCTTTGTTTATGCCGGATGCGGCGTTAACGCCTTATCCGGCCTACAAAATCGTGCAAATTCAATACATTGCAGAATTTTc                                       <  1:2313453/100‑1 (MQ=255)
                                                              cAGATTGTTGACAAAGTGCGCTTTGTTTATGCCGGATGCGGCGTTAACGCCTTATCCGGCCTACAAAATCGTGCAAATTCAATACATTGCAGAATTTTCg                                      <  2:1088951/100‑1 (MQ=255)
                                                              cAGATTGTTGACAAAGTGCGCTTTGTTTATGCCGGATGCGGCGTTAACGCCTTATCCGGCCTACAAAATCGTGCAAATTCAATACATTGCAGAATTTTCg                                      <  1:429413/100‑1 (MQ=255)
                                                              cAGATTGCTGACAACGTGCGCGTTGTTCATGCCGGATGCGGCGTGAACGCCTTATCCGGCCTACAAAGTCGTATAAACTGCAATACATTgaatttatacg                                        <  2:709753/100‑11 (MQ=14)
                                                               aGATTGTTGACAAAGTGCGCTTTGTTTATGCCGGATGCGGCGTTAACGCCTTATCCGGCCTACAAAATCGTGCAAATTCAATACATTGCAGAATTTTCGt                                     <  2:953252/100‑1 (MQ=255)
                                                                gATTGTTGACAAAGTGCGCTTTGTTTATGCCGGATGCGGCGTTAACGCCTTATCCGGCCTACAAAATCGTGCAAATTCAATACATTGCAGAATTTTCGTa                                    >  2:1622806/1‑100 (MQ=255)
                                                                 aTTGTTGACAAAGTGCGCTTTGTTTATGCCGGATGCGGCGTTAACGCCTTATCCGGCCTACAAAATCGTGCAAATTCAATACATTGCAGAATTTTCGTAg                                   <  1:714317/100‑1 (MQ=255)
                                                                  gcataaaCAAAGCGCACTTTGTTTATGCCGGATGCGGCGTTAACGCCTTATCCGGCCTACAAAATCGTGCAAATTCAATACATTGCAGAATTTTCGTAgg                                  >  1:404911/7‑100 (MQ=18)
                                                                   tgttgACAAAGTGCGCTTTGTTTATGCCGGATGCGGCGTTAACGCCTTATCCGGCCTACAAAATCGTGCAAATTCAATACATTGCAGAATTTTCGTAGGc                                 >  2:1644753/1‑100 (MQ=255)
                                                                    gttgACAAAGTGCGCTTTGTTTATGCCGGATGCGGCGTTAACGCCTTATCCGGCCTACAAAATCGTGCAAATTCAATACATTGCAGAATTTTCGTAGGcc                                <  2:571259/100‑1 (MQ=255)
                                                                     ttgACAAAGTGCGCTTTGTTTATGCCGGATGCGGCGTTAACGCCTTATCCGGCCTACAAAATCGTGCAAATTCAATACATTGCAGAATTTTCGTAGGCCt                               <  1:2063062/100‑1 (MQ=255)
                                                                      tgACAAAGTGCGCTTTGTTTATGCCGGATGCGGCGTTAACGCCTTATCCGGCCTACAAAATCGTGCAAATTCAATACATTGCAGAATTTTCGTAGGCCTg                              >  1:2301212/1‑100 (MQ=255)
                                                                      tgACAAAGTGCGCTTTGTTTATGCCGGATGCGGCGTTAACGCCTTATCCGGCCTACAAAATCGTGCAAATTCAATACATTGCAGAATTTTCGTAGGCCTg                              <  1:2391242/100‑1 (MQ=255)
                                                                      tgACAAAGTGCGCTTTGTTTATGCCGGATGCGGCGTTAACGCCTTATCCGGCCTACAAAATCGTGCAAATTCAATACATTGCAGAATTTTCGTAGGCCTg                              >  2:2290443/1‑100 (MQ=255)
                                                                      ngACAAAGTGCGCTTTGTTTATGCCGGATGCGGCGTTAACGCCTTATCCGGCCTACAAAATCGTGCAAATTCAATACATTGCAGAATTTTCGTAGGCCTg                              >  1:1416603/2‑100 (MQ=255)
                                                                       gACAAAGTGCGCTTTGTTTATGCCGGATGCGGCGTTAACGCCTTATCCGGCCTACAAACTCGTGCAAATTCAATACATTGCAGAATTTTCGTAGGCCTGa                             <  1:1644745/100‑1 (MQ=255)
                                                                         cAAAGTGCGCTTTGTTTATGCCGGATGCGGCGTTAACGCCTTATCCGGCCTACAAAATCGTGCAAATTCAATACATTGCAGAATTTTCGTAGGCCTGATa                           <  2:1718389/100‑1 (MQ=255)
                                                                               gcgcTTTGTTTATGCCGGATGCGGCGTTAACGCCTTATCCGGCCTACAAAATCGTGCAAATTCAATACATTGCAGAATTTTCGTAGGCCTGATAAGCGGa                     >  2:1422895/1‑100 (MQ=255)
                                                                                 gcTTTGTTTATGCCGGATGCGGCGTTAACGCCTTATCCGGCCTACAAAATCGTGCAAATTCAATACATTGCAGAATTTTCGTAGGCCTGATAAGCGGAgc                   <  1:2398827/100‑1 (MQ=255)
                                                                                  cTTTGTTTATGCCGGATGCGGCGTTAACGCCTTATCCGGCCTACAAAATCGTGCAAATTCAATACATTGCAGAATTTTCGTAGGCCTGATAAGCGGAgcg                  >  2:763133/1‑100 (MQ=255)
                                                                                  cTTTGTTTATGCCGGATGCGGCGTTAACGCCTTATCCGGCCTACAAAATCGTGCAAATTCAATACATTGCAGAATTTTCGTAGGCCTGATAAGCGGAgcg                  <  1:1821906/100‑1 (MQ=255)
                                                                                   cttttgtaggCCGGATAAGGCGTTTACGCCGCATCCGGCCTACAAAAGCATGCAAATTCAATATATTGCAGAGTTCACGTAGGCCTGATAAGCGTAgcgc                 >  2:695212/10‑100 (MQ=17)
                                                                                      ttgtaggCCGGATAAGGCGTTCACGCCGCATCCGGCCTACAAAGTCGTATAAATTCAATGTATTGCAGTTTATTGGTAGGCCTGATAAGCGCAGCGCATc              >  2:2350799/7‑100 (MQ=21)
                                                                                        ttATGCCGGATGCGGCGTTAACGCCTTATCCGGCCTACAAAATCGTGCAAATTCAATACATTGCAGAATTTTCGTAGGCCTGATAAGCGGAGCGCATCAg            <  1:2478080/100‑1 (MQ=255)
                                                                                        gtaggCCGGATAAGGCGTTCACGCCGCATCCGGCCTACAAGTTCGTGCAAATTCAATAAATTGCAATATGACGTAGGCCTGATAAGCGTAGCGCATCAgg           >  1:1945164/5‑100 (MQ=18)
                                                                                        gtaggCCGGATAAGGCGTTCACGCCGCATCCGGCCTACAAAACTGTGCAAATTCAATATATTGCAGGAAACACGTAGGCCTGATAAGCGAAGCGCATCAg            >  2:1450379/5‑100 (MQ=18)
                                                                                        gtaggCCGGATAAGGCGTTCACGCCGCATCCGGCCTACAAATCATGCAAATCCAAAATATTGCAGGAGTTATGTAGGCCTGATAAGCGTAGCGCATCAgg           >  2:253099/5‑100 (MQ=21)
                                                                                        gtaggCCGGATAAGGCGTTCACGCCGCATCCGGCCTACAAACCGCGCAAATTCAATATATTGCGGAGAAAATGTAGGCCTGATAAGCGTAGCGCATCAgg           >  1:2513595/5‑100 (MQ=18)
                                                                                         taggCCGGATAAGGCGTTCACGCCGCATCCGGCCTACAAAATTGTGTAAATTCAATAAATTGATCGATCCACGTAGGCCTGATAAGCGTAGCGCATCAgg           >  2:373817/4‑100 (MQ=21)
                                                                                          aggCCGGATAAGGCGTTCACTCCGCATCCGGCCTACAAAAGCACGCAAATTCAATATATTGCAGAGATTGCGTAGGCCTGATAAGCGTAGCGCATCAGGc          >  1:1820111/3‑100 (MQ=11)
                                                                                          aggCCGGATAAGGCGTTAACGCCGCATCCGGCCTACAAAATCGTGCAAATTCAATACATTGCAGAATTTTCGTAGGCCTGATAAGCGGAGCGCATCAGGc          >  2:2464123/3‑100 (MQ=18)
                                                                                           tGCCGGATGCGGCGTTAACGCCTTATCCGGCCTACAAAATCGTGCAAATTCAATACATTGCAGAATTTTCGTAGTCCTGATAAGCGTAGCGCATCAGGCa         >  2:438892/1‑100 (MQ=255)
                                                                                           tGCCGGATGCGGCGTTAACGCCTTATCCGGCCTACAAAATCGTGCAAATTCAATACATTGCAGAATTTTCGTAGGCCTGATAAGCGGAGCGCATCAGGCa         <  2:161018/100‑1 (MQ=255)
                                                                                           ngCCGGATAAGGCGTTCACGCCGCATCCGGCCTACAAAACTGTGCAAATTCAATATATTGCAGGAAACACGTAGGCCTGATAAGCGAAGCGCATCAGGCa         >  2:2645807/2‑100 (MQ=18)
                                                                                           ggCCGGATAAGGCGTTTACGCCGCATCCGGCCTACAAAAGCACGCAAACTCAATATATTGCAGAGATCATGTAGGCCTGATAAGCGTAGCGCATCAGGCa         >  1:2921783/2‑100 (MQ=17)
                                                                                           ggCCGGATAAGGCGTTCACGCCGCATCCGGCCTACAAAGTCGTATAAATTCAATGTATTGCAGTTTATTGGTAGGCCTGATAAGCGCAGCGCATCAGGCa         >  1:2889630/2‑100 (MQ=21)
                                                                                           ggCCGGATAAGGCGTTCACGCCGCATCCGGCCTACAAAGTCGTATAAATTCAATGTATTGCAGTTTATTGGTAGGCCTGATAAGCGCAGCGCATCAGGCa         >  2:256511/2‑100 (MQ=25)
                                                                                           ggCCGGATAAGGCGTTCACGCCGCATCCGGCCTACAAAACTGTGCAAATTCAATATATTGCAGGAAACACGTAGGCCTGATAAGCGAAGCGCATCAGGCa         >  2:545733/2‑100 (MQ=18)
                                                                                           ggCCGGATAAGGCGTTAACGCCGCATCCGGCCTACAAAATCGTGCAAATTCAATACATTGCAGAATTTTCGTAGGCCTGATAAGCGGAGCGCATCAGGCa         >  2:173911/2‑100 (MQ=18)
                                                                                            gCCGGATGCGGCGTTAACGCCTTATCCGGCCTACAAAATCGTGCAAATTCAATACATTGGCGAATTTTCGTAGGCCTGATCAGCAGAGCGCATCAAGcct        >  1:283062/1‑98 (MQ=17)
                                                                                            gCCGGATAAGGCGTTTACGCCGCATCCGGCCTACAAAAGCATGCAAATTCAATATATTGCAGAGTTCACGTAGGCCTGATAAGCGTAGCGCATCAGGCag        >  2:1900194/1‑99 (MQ=17)
                                                                                            gCCGGATAAGGCGTTTACGCCGCATCCGGCCTACAAAAGCATGCAAATTCAATATATTGCAGAGTTCACGTAGGCCTGATAAGCGTAGCGCATCAGGCag        <  1:695210/100‑2 (MQ=17)
                                                                                            gCCGGATAAGGCGTTCACGCCGCATCCGGCCTACAAAAATCTTGCCAATTCAATATATTGCAGGACCATGTAGGCCTGATAAGCGTAGCGCATCAGGCag        >  2:944004/1‑99 (MQ=18)
                                                                                            gCCGGATAAGGCGTTCACGCCGCATCCGGCCTACAAAGTCGTATAACTTCAATGTATTGCAGTTTATTGGTAGGCCTGATAAGCGCAGCGCATCAGGCaa        >  1:1369201/1‑99 (MQ=18)
                                                                                            gCCGGATAAGGCGTTCACGCCGCATCCGGCCTACAAAGTCGTATAAATTCAATGTATTGCAGTTTATTGGTAGGCCTGATAAGCGCAGCGCATCAGGCaa        >  2:2669950/1‑99 (MQ=25)
                                                                                            gCCGGATAAGGCGTTCACGCCGCATCCGGCCTACAAAAACATGCAAATTCAATAGATTGCAGAGATTATGTAGGCCTGATAAGCGTAGCGCATCAGGCaa        >  1:2236509/1‑99 (MQ=17)
                                                                                            gCCGGATAAGGCGTTCACGCCGCATCCGGCCTACAAATCATGCAAATCCAAAATATTGCAGGAGTTATGTAGGCCTGATAAGCGTAGCGCATCAGGCAtt       >  1:990152/1‑100 (MQ=21)
                                                                                                   |                                                                                                     
CAGGGAATGAAGTATTGCTGCCGGATGGCGCTGCACGCAGAGGCACAGCGCCACCGGGTATACAGATTGTTGACAAAGTGCGCTTTGTTTATGCCGGATGCGGCGTTAACGCCTTATCCGGCCTACAAAATCGTGCAAATTCAATACATTGCAGAATTTTCGTAGGCCTGATAAGCGGAGCGCATCAGGCATTTTTGC  >  NC_000913/3705871‑3706068

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: