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.

CAGGGAATGAAGTATTGCTGCCGGATGGCGCTGCACGCAGAGGCACAGCGCCACCGGGTATACAGATTGTTGACAAAGTGCGCTTTGTTTATGCCGGATGCGGCGTTAACGCCTTATCCGGCCTAC‑AAAATCGTGCAAA‑TTCAATACATTGCAGA‑ATTTTCGTAGGCCTGATAAGCGGAGCGCATCAGGCATTTTTGC  >  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)
                           gcgcTGCACGCAGAGGCACAGCGCCACCGGGTATACAGATTGTTGACAAAGTGCGCTTTGTTTATGCCGGATGCGGCGTTAACGCCTTATCCGGCCTAC‑a                                                                           >  2:1131162/1‑100 (MQ=255)
                           gcgcTGCACGCAGAGGCACAGCGCCACCGGGTATACAGATTGTTGACAAAGTGCGCTTTGTTTATGCCCGATGCGGCGTTAACGCCTTATCCGGCCTaaa                                                                            <  2:2017278/100‑3 (MQ=255)
                              cTGCACGCAGAGGCACAGCGCCACCGGGTATACAGATTGTTGACAAAGTGCGCTTTGTTTATGCCGGATGCGGCGTTAACGCCTTATCCGGCCTAC‑aaaa                                                                        <  1:126618/100‑1 (MQ=255)
                               tGCACGCAGAGGCACAGCGCCACCGGGTATACAGATTGTTGACAAAGTGCGCTTTGTTTATGCCGGATGCGGCGTTAACGCCTTATCCGGCCTagatcgg                                                                        >  1:959890/1‑94 (MQ=255)
                                gCACGCAGAGGCACAGCGCCACCGGGTATACAGATTGTTGCCAAAGTGCGCTTTGTTTATGCCGGATGCGGCGTTAACGCCTTATCCGGCCTAC‑AAAATc                                                                      >  1:781395/1‑100 (MQ=255)
                                 cACGCAGAGGCACAGCGCCACCGGGTATACAGATTGTTGACAAAGTGCGCTTTGTTTATGCCGGATGCGGCGTTAACGCCTTATCCGGCCTAC‑AAAATCg                                                                     <  2:2132898/100‑1 (MQ=255)
                                 cACGCAGAGGCACAGCGCCACCGGGTATACAGATTGTTGACAAAGTGCGCTTTGTTTATGCCGGATGCGGCGTTAACGACTTATCCGGCCTAC‑AAAATCg                                                                     <  2:2414023/100‑1 (MQ=255)
                                   cGCAGAGGCACAGCGCCACCGGGTATACAGATTGTTGACAAAGTGCGCTTTGTTTATGCCGGATGCGGCGTTAACGCCTTATCCGGCCTAC‑AAAATCGTg                                                                   <  1:2629738/100‑1 (MQ=255)
                                    gCAGAGGCACAGCGCCACCGGGTATACAGATTGTTGACAAAGTGCGCTTTGTTTATGCCGGATGCGGCGTTAACGCCTTATCCGGCCTAC‑AAAATCGTGc                                                                  >  1:1770296/1‑100 (MQ=255)
                                     cAGAGGCACAGCGCCACCGGGTATACAGATTGTTGACAAAGTGCGCTTTGTTTATGCCGGATGCGGCGTTAACGCCTTATCCGGCCTAC‑AAAATCGTGCa                                                                 >  1:2014835/1‑100 (MQ=255)
                                     cAGAGGCACAGCGCCACCGGGTATACAGATTGTTGACAAAGTGCGCTTTGTTTATGCCGGACGCGGCGTTAACGCCTTATCCGGCCTAC‑AAAATCGTGCa                                                                 <  1:1853837/100‑1 (MQ=255)
                                        agGCACAGCGCCACCGGGTATACAGATTGTTGACAAAGTGCGCTTTGTTTATGCCGGATGCGGCGTTAACGCCTTATCCGGCCTAC‑AAAATCGTGCAAA‑t                                                             <  2:205921/100‑1 (MQ=255)
                                         ggCACAGCGCCACCGGGTATACAGATTGTTGACAAAGTGCGCTTTGTTTATGCCGGATGCGGCGTTAACGCCTTATCCGGCCTAC‑AAAATCGTGCAAA‑tt                                                            >  1:1504437/1‑100 (MQ=255)
                                           cacaGCGCCACCGGGTATACAGATTGTTGACAAAGTGCGCTTTGTTTATGCCGGATGCGGCGTTAACGCCTTATCCGGCCTAC‑AAAATCGTGCAAA‑TTCa                                                          <  1:2534333/100‑1 (MQ=255)
                                           cacaGCGCCACCGGGTATACAGATTGTTGACAAAGTGCGCTTTGTTTATGCCGGATGCGGCGTTAACGCCTTATCCGGCCTAC‑AAAATCGTGCAAA‑TTCa                                                          >  2:2796692/1‑100 (MQ=255)
                                             caGCGCCACCGGGTATACAGATTGTTGACAAAGTGCGCTTTGTTTATGCCGGATGCGGCGTTAACGCCTTATCCGGCCTAC‑AAAATCGTGCAAA‑TTCAAt                                                        <  2:2941306/100‑1 (MQ=255)
                                              aGCGCCACCGGGTATACAGATTGTTGACAAAGTGCGCTTTGTTTATGCCGGATGCGGCGTTAACGCCTTATCCGGCCTAC‑AAAATCGTGCAAA‑TTCAATa                                                       <  1:1077300/100‑1 (MQ=255)
                                              aGCGCCACCGGGTATACAGATTGTTGACAAAGTGCGCTTTGTTTATGCCGGATGCGGCGTTAACGCCTTATCCGGCCTAC‑AAAATCGTGCAAA‑TTCAATa                                                       <  1:778510/100‑1 (MQ=255)
                                              aGCGCCACCGGGTATACAGATTGTTGACAAAGTGCGCTTTGTTTATGCCGGATGCGGCGTTAACGCCTTATCCGGCCTAC‑AAAATCGTGCAAA‑TTCAATa                                                       <  1:363991/100‑1 (MQ=255)
                                               gcgcCACCGGGTATACAGATTGTTGACAAAGTGCGCTTTGTTTATGCCGGATGCGGCGTTAACGCCTTATCCGGCCTAC‑AAAATCGTGCAAA‑TTCAATAc                                                      <  1:647090/100‑1 (MQ=255)
                                               gcgcCACCGGGTATACAGATTGTTGACAAAGTGCGCTTTGTTTATGCCGGATGCGGCGTTAACGCCTTATCCGGCCTAC‑AAAATCGTGCAAA‑TTCAATAc                                                      >  2:302547/1‑100 (MQ=255)
                                                  ccACCGGGTATACAGATTGTTGACAAAGTGCGCTTTGTTTATGCCGGATGCGGCGTTAACGCCTTATCCGGCCTAC‑AAAATCGTGCAAA‑TTCAATACAtt                                                   <  2:354389/100‑1 (MQ=255)
                                                   cACCGGGTATACAGATTGTTGACAAAGTGCGCTTTGTTTATGCCGGATGCGGCGTTAACGCCTTATCCGGCCTAC‑AAAATCGTGCAAA‑TTCAATACATTg                                                  >  2:122053/1‑100 (MQ=255)
                                                    aCCGGGTATACAGATTGTTGACAAAGTGCGCTTTGTTTATGCCGGATGCGGCGTTAACGCCTTATCCGGCCTAC‑AAAATCGTGCAAA‑TTCAATACATTGc                                                 <  1:1047383/100‑1 (MQ=255)
                                                     ccGGGTATACAGATTGTTGACAAAGTGCGCTTTGTTTATGCCGGATGCGGCGTTAACGCCTTATCCGGCCTAC‑AAAATCGTGCAAA‑TTCAATACATTGCa                                                <  2:821528/100‑1 (MQ=255)
                                                          tataCAGATTGTTGACAAAGTGCGCTTTGTTTATGCCGGATGCGGCGTTAACGCCTTATCCGGCCTAC‑AAAATCGTGCAAA‑TTCAATACATTGCAGA‑Att                                          >  2:1294297/1‑100 (MQ=255)
                                                             aCAGATTGTTGACAAAGTGCGCTTTGTTTATGCCGGATGCGGCGTTAACGCCTTATCCGGCCTAC‑AAAATCGTGCAAA‑TTCAATACATTGCAGA‑ATTTTc                                       <  1:1697648/100‑1 (MQ=255)
                                                             aCAGATTGTTGACAAAGTGCGCTTTGTTTATGCCGGATGCGGCGTTAACGCCTTATCCGGCCTAC‑AAAATCGTGCAAA‑TTCAATACATTGCAGA‑ATTTTc                                       <  1:2313453/100‑1 (MQ=255)
                                                              cAGATTGTTGACAAAGTGCGCTTTGTTTATGCCGGATGCGGCGTTAACGCCTTATCCGGCCTAC‑AAAATCGTGCAAA‑TTCAATACATTGCAGA‑ATTTTCg                                      <  2:1088951/100‑1 (MQ=255)
                                                              cAGATTGTTGACAAAGTGCGCTTTGTTTATGCCGGATGCGGCGTTAACGCCTTATCCGGCCTAC‑AAAATCGTGCAAA‑TTCAATACATTGCAGA‑ATTTTCg                                      <  1:429413/100‑1 (MQ=255)
                                                              cAGATTGCTGACAACGTGCGCGTTGTTCATGCCGGATGCGGCGTGAACGCCTTATCCGGCCTAC‑AAAGTCGTATAAACTGCAATACATTgaatttatacg                                        <  2:709753/100‑11 (MQ=14)
                                                               aGATTGTTGACAAAGTGCGCTTTGTTTATGCCGGATGCGGCGTTAACGCCTTATCCGGCCTAC‑AAAATCGTGCAAA‑TTCAATACATTGCAGA‑ATTTTCGt                                     <  2:953252/100‑1 (MQ=255)
                                                                gATTGTTGACAAAGTGCGCTTTGTTTATGCCGGATGCGGCGTTAACGCCTTATCCGGCCTAC‑AAAATCGTGCAAA‑TTCAATACATTGCAGA‑ATTTTCGTa                                    >  2:1622806/1‑100 (MQ=255)
                                                                 aTTGTTGACAAAGTGCGCTTTGTTTATGCCGGATGCGGCGTTAACGCCTTATCCGGCCTAC‑AAAATCGTGCAAA‑TTCAATACATTGCAGA‑ATTTTCGTAg                                   <  1:714317/100‑1 (MQ=255)
                                                                  gcataaaCAAAGCGCACTTTGTTTATGCCGGATGCGGCGTTAACGCCTTATCCGGCCTAC‑AAAATCGTGCAAA‑TTCAATACATTGCAGA‑ATTTTCGTAgg                                  >  1:404911/7‑100 (MQ=18)
                                                                   tgttgACAAAGTGCGCTTTGTTTATGCCGGATGCGGCGTTAACGCCTTATCCGGCCTAC‑AAAATCGTGCAAA‑TTCAATACATTGCAGA‑ATTTTCGTAGGc                                 >  2:1644753/1‑100 (MQ=255)
                                                                    gttgACAAAGTGCGCTTTGTTTATGCCGGATGCGGCGTTAACGCCTTATCCGGCCTAC‑AAAATCGTGCAAA‑TTCAATACATTGCAGA‑ATTTTCGTAGGcc                                <  2:571259/100‑1 (MQ=255)
                                                                     ttgACAAAGTGCGCTTTGTTTATGCCGGATGCGGCGTTAACGCCTTATCCGGCCTAC‑AAAATCGTGCAAA‑TTCAATACATTGCAGA‑ATTTTCGTAGGCCt                               <  1:2063062/100‑1 (MQ=255)
                                                                      tgACAAAGTGCGCTTTGTTTATGCCGGATGCGGCGTTAACGCCTTATCCGGCCTAC‑AAAATCGTGCAAA‑TTCAATACATTGCAGA‑ATTTTCGTAGGCCTg                              >  1:2301212/1‑100 (MQ=255)
                                                                      tgACAAAGTGCGCTTTGTTTATGCCGGATGCGGCGTTAACGCCTTATCCGGCCTAC‑AAAATCGTGCAAA‑TTCAATACATTGCAGA‑ATTTTCGTAGGCCTg                              <  1:2391242/100‑1 (MQ=255)
                                                                      tgACAAAGTGCGCTTTGTTTATGCCGGATGCGGCGTTAACGCCTTATCCGGCCTAC‑AAAATCGTGCAAA‑TTCAATACATTGCAGA‑ATTTTCGTAGGCCTg                              >  2:2290443/1‑100 (MQ=255)
                                                                      ngACAAAGTGCGCTTTGTTTATGCCGGATGCGGCGTTAACGCCTTATCCGGCCTAC‑AAAATCGTGCAAA‑TTCAATACATTGCAGA‑ATTTTCGTAGGCCTg                              >  1:1416603/2‑100 (MQ=255)
                                                                       gACAAAGTGCGCTTTGTTTATGCCGGATGCGGCGTTAACGCCTTATCCGGCCTAC‑AAACTCGTGCAAA‑TTCAATACATTGCAGA‑ATTTTCGTAGGCCTGa                             <  1:1644745/100‑1 (MQ=255)
                                                                         cAAAGTGCGCTTTGTTTATGCCGGATGCGGCGTTAACGCCTTATCCGGCCTAC‑AAAATCGTGCAAA‑TTCAATACATTGCAGA‑ATTTTCGTAGGCCTGATa                           <  2:1718389/100‑1 (MQ=255)
                                                                               gcgcTTTGTTTATGCCGGATGCGGCGTTAACGCCTTATCCGGCCTAC‑AAAATCGTGCAAA‑TTCAATACATTGCAGA‑ATTTTCGTAGGCCTGATAAGCGGa                     >  2:1422895/1‑100 (MQ=255)
                                                                                 gcTTTGTTTATGCCGGATGCGGCGTTAACGCCTTATCCGGCCTAC‑AAAATCGTGCAAA‑TTCAATACATTGCAGA‑ATTTTCGTAGGCCTGATAAGCGGAgc                   <  1:2398827/100‑1 (MQ=255)
                                                                                  cTTTGTTTATGCCGGATGCGGCGTTAACGCCTTATCCGGCCTAC‑AAAATCGTGCAAA‑TTCAATACATTGCAGA‑ATTTTCGTAGGCCTGATAAGCGGAgcg                  >  2:763133/1‑100 (MQ=255)
                                                                                  cTTTGTTTATGCCGGATGCGGCGTTAACGCCTTATCCGGCCTAC‑AAAATCGTGCAAA‑TTCAATACATTGCAGA‑ATTTTCGTAGGCCTGATAAGCGGAgcg                  <  1:1821906/100‑1 (MQ=255)
                                                                                   cttttgtaggCCGGATAAGGCGTTTACGCCGCATCCGGCCTAC‑AAAAGCATGCAAA‑TTCAATATATTGCAGA‑GTTCACGTAGGCCTGATAAGCGTAgcgc                 >  2:695212/10‑100 (MQ=17)
                                                                                      ttgtaggCCGGATAAGGCGTTCACGCCGCATCCGGCCTAC‑AAAGTCGTATAAA‑TTCAATGTATTGCAGT‑TTATTGGTAGGCCTGATAAGCGCAGCGCATc              >  2:2350799/7‑100 (MQ=21)
                                                                                        ttATGCCGGATGCGGCGTTAACGCCTTATCCGGCCTAC‑AAAATCGTGCAAA‑TTCAATACATTGCAGA‑ATTTTCGTAGGCCTGATAAGCGGAGCGCATCAg            <  1:2478080/100‑1 (MQ=255)
                                                                                        gtaggCCGGATAAGGCGTTCACGCCGCATCCGGCCTAC‑AAGTTCGTGCAAA‑TTCAATAAATTGCA‑A‑TATGACGTAGGCCTGATAAGCGTAGCGCATCAgg           >  1:1945164/5‑100 (MQ=18)
                                                                                        gtaggCCGGATAAGGCGTTCACGCCGCATCCGGCCTAC‑AAAACTGTGCAAA‑TTCAATATATTGCAGG‑AAACACGTAGGCCTGATAAGCGAAGCGCATCAg            >  2:1450379/5‑100 (MQ=18)
                                                                                        gtaggCCGGATAAGGCGTTCACGCCGCATCCGGCCTAC‑‑AAATCATGCAAA‑TCCAAAATATTGCAGG‑AGTTATGTAGGCCTGATAAGCGTAGCGCATCAgg           >  2:253099/5‑100 (MQ=21)
                                                                                        gtaggCCGGATAAGGCGTTCACGCCGCATCCGGCCTAC‑‑AAACCGCGCAAA‑TTCAATATATTGCGGA‑GAAAATGTAGGCCTGATAAGCGTAGCGCATCAgg           >  1:2513595/5‑100 (MQ=18)
                                                                                         taggCCGGATAAGGCGTTCACGCCGCATCCGGCCTAC‑AAAATTGTGTAAA‑TTCAATAAATTGATCG‑ATCCACGTAGGCCTGATAAGCGTAGCGCATCAgg           >  2:373817/4‑100 (MQ=21)
                                                                                          aggCCGGATAAGGCGTTCACTCCGCATCCGGCCTAC‑AAAAGCACGCAAA‑TTCAATATATTGCAGA‑GATTGCGTAGGCCTGATAAGCGTAGCGCATCAGGc          >  1:1820111/3‑100 (MQ=11)
                                                                                          aggCCGGATAAGGCGTTAACGCCGCATCCGGCCTAC‑AAAATCGTGCAAA‑TTCAATACATTGCAGA‑ATTTTCGTAGGCCTGATAAGCGGAGCGCATCAGGc          >  2:2464123/3‑100 (MQ=18)
                                                                                           tGCCGGATGCGGCGTTAACGCCTTATCCGGCCTAC‑AAAATCGTGCAAA‑TTCAATACATTGCAGA‑ATTTTCGTAGTCCTGATAAGCGTAGCGCATCAGGCa         >  2:438892/1‑100 (MQ=255)
                                                                                           tGCCGGATGCGGCGTTAACGCCTTATCCGGCCTAC‑AAAATCGTGCAAA‑TTCAATACATTGCAGA‑ATTTTCGTAGGCCTGATAAGCGGAGCGCATCAGGCa         <  2:161018/100‑1 (MQ=255)
                                                                                           ngCCGGATAAGGCGTTCACGCCGCATCCGGCCTAC‑AAAACTGTGCAAA‑TTCAATATATTGCAGG‑AAACACGTAGGCCTGATAAGCGAAGCGCATCAGGCa         >  2:2645807/2‑100 (MQ=18)
                                                                                           ggCCGGATAAGGCGTTTACGCCGCATCCGGCCTAC‑AAAAGCACGCAAA‑CTCAATATATTGCAGA‑GATCATGTAGGCCTGATAAGCGTAGCGCATCAGGCa         >  1:2921783/2‑100 (MQ=17)
                                                                                           ggCCGGATAAGGCGTTCACGCCGCATCCGGCCTAC‑AAAGTCGTATAAA‑TTCAATGTATTGCAGT‑TTATTGGTAGGCCTGATAAGCGCAGCGCATCAGGCa         >  1:2889630/2‑100 (MQ=21)
                                                                                           ggCCGGATAAGGCGTTCACGCCGCATCCGGCCTAC‑AAAGTCGTATAAA‑TTCAATGTATTGCAGT‑TTATTGGTAGGCCTGATAAGCGCAGCGCATCAGGCa         >  2:256511/2‑100 (MQ=25)
                                                                                           ggCCGGATAAGGCGTTCACGCCGCATCCGGCCTAC‑AAAACTGTGCAAA‑TTCAATATATTGCAGG‑AAACACGTAGGCCTGATAAGCGAAGCGCATCAGGCa         >  2:545733/2‑100 (MQ=18)
                                                                                           ggCCGGATAAGGCGTTAACGCCGCATCCGGCCTAC‑AAAATCGTGCAAA‑TTCAATACATTGCAGA‑ATTTTCGTAGGCCTGATAAGCGGAGCGCATCAGGCa         >  2:173911/2‑100 (MQ=18)
                                                                                            gCCGGATGCGGCGTTAACGCCTTATCCGGCCTAC‑AAAATCGTGCAAA‑TTCAATACATTGGCGA‑ATTTTCGTAGGCCTGATCAGCAGAGCGCATCAAGcct        >  1:283062/1‑98 (MQ=17)
                                                                                            gCCGGATAAGGCGTTTACGCCGCATCCGGCCTAC‑AAAAGCATGCAAA‑TTCAATATATTGCAGA‑GTTCACGTAGGCCTGATAAGCGTAGCGCATCAGGCag        >  2:1900194/1‑99 (MQ=17)
                                                                                            gCCGGATAAGGCGTTTACGCCGCATCCGGCCTAC‑AAAAGCATGCAAA‑TTCAATATATTGCAGA‑GTTCACGTAGGCCTGATAAGCGTAGCGCATCAGGCag        <  1:695210/100‑2 (MQ=17)
                                                                                            gCCGGATAAGGCGTTCACGCCGCATCCGGCCTACAAAAATCTTGCCAA‑TTCAATATATTGCAGG‑ACCAT‑GTAGGCCTGATAAGCGTAGCGCATCAGGCag        >  2:944004/1‑99 (MQ=18)
                                                                                            gCCGGATAAGGCGTTCACGCCGCATCCGGCCTAC‑AAAGTCGTATAAC‑TTCAATGTATTGCAGT‑TTATTGGTAGGCCTGATAAGCGCAGCGCATCAGGCaa        >  1:1369201/1‑99 (MQ=18)
                                                                                            gCCGGATAAGGCGTTCACGCCGCATCCGGCCTAC‑AAAGTCGTATAAA‑TTCAATGTATTGCAGT‑TTATTGGTAGGCCTGATAAGCGCAGCGCATCAGGCaa        >  2:2669950/1‑99 (MQ=25)
                                                                                            gCCGGATAAGGCGTTCACGCCGCATCCGGCCTAC‑AAAAACATGCAAA‑TTCAATAGATTGCAGAGATTAT‑GTAGGCCTGATAAGCGTAGCGCATCAGGCaa        >  1:2236509/1‑99 (MQ=17)
                                                                                            gCCGGATAAGGCGTTCACGCCGCATCCGGCCTAC‑‑AAATCATGCAAA‑TCCAAAATATTGCAGG‑AGTTATGTAGGCCTGATAAGCGTAGCGCATCAGGCAtt       >  1:990152/1‑100 (MQ=21)
                                                                                                   |                                                                                                     
CAGGGAATGAAGTATTGCTGCCGGATGGCGCTGCACGCAGAGGCACAGCGCCACCGGGTATACAGATTGTTGACAAAGTGCGCTTTGTTTATGCCGGATGCGGCGTTAACGCCTTATCCGGCCTAC‑AAAATCGTGCAAA‑TTCAATACATTGCAGA‑ATTTTCGTAGGCCTGATAAGCGGAGCGCATCAGGCATTTTTGC  >  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: ‑