Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009133,705,9710CA26.8% 414.9 / 148.7 180intergenic (‑181/+127)dppB/dppAdipeptide/heme ABC transporter permease/dipeptide/heme ABC transporter periplasmic binding protein; dipeptide chemotaxis receptor
Reads supporting (aligned to +/- strand):  ref base C (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 = 3.38e-01

Only 100 of 196 total aligned reads displayed.

AGGGAATGAAGTATTGCTGCCGGATGGCGCTGCACGCAGAGGCACAGCGCCACCGGGTATACAGATTGTTGACAAAGTGCGCTTTGTTTATGCCGGATGCGGCGTTAACGCCTTATCCGGCCTAC‑AAAATCGTGCAAA‑TTCAATACATTGCA‑GA‑AT‑TTTCGTAGGCCTGATAAGCGGAGCGCATCAGGCATTTTTGCCT  >  NC_000913/3705872‑3706070
                                                                                                   |                                                                                                        
 gggAATGAAGTATTGCTGCCGGATGGCGCTGCACGCAGAGGCACAGCGCCACCGGGTATACAGATTGTTGACAAAGTGCGCTTTGTTTATGCCGGATgcg                                                                                                         >  2:2478091/1‑100 (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:2590650/1‑100 (MQ=255)
        aaGTATTGCTGCCGGATGGCGCTGCACGCAGAGGCACAGCGCCACCGGGTATACAGATTGTTGACAAAGTGCGCTTTGTTTATGCCGGATGCGGCGTTaa                                                                                                  >  2:1853845/1‑100 (MQ=255)
         aGTATTGCTGCCGGATGGCGCTGCACGCAGAGGCACAGCGCCACCGGGTATACAGATTGTTGACAAAGTGCGCTTTGTTTATGCCGGATGCGGCGTTAAc                                                                                                 <  2:1235513/100‑1 (MQ=255)
         aGTATTGCTGCCGGATGGCGCTGCACGCAGAGGCACAGCGCCACCGGGTATACAGATTGTTGACAAAGTGCGCTTTGTTTATGCCGGATGCGGCGTTAAc                                                                                                 >  1:469123/1‑100 (MQ=255)
          gTATTGCTGCCGGATGGCGCTGCACGCAGAGGCACAGCGCCACCGGGTATACAGATTGTTGACAAAGTGCGCTTTGTTTATGCCGGATGCGGCGTTAACg                                                                                                >  2:624210/1‑100 (MQ=255)
           tATTGCTGCCGGATGGCGCTGCACGCAGAGGCACAGCGCCACCGGGTATACAGATTGTTGACAAAGTGCGCTTTGTTTATGCCGGATGCGGCGTTAACGc                                                                                               <  2:2058736/100‑1 (MQ=255)
            aTTGCTGCCGGATGGCGCTGCACGCAGAGGCACAGCGCCACCGGGTATACAGATTGTTGACAAAGTGCGCTTTGTTTATGCCGGATGCGGCGTTAACGcc                                                                                              <  2:579148/100‑1 (MQ=255)
            aTTGCTGCCGGATGGCGCTGCACGCAGAGGCACAGCGCCACCGGGTATACAGATTGTTGACAAAGTGCGCTTTGTTTATGCCGGATGCGGCGTTAACGcc                                                                                              <  1:1198050/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)
                   ccGGATGGCGCTGCACGCAGAGGCACAGCGCCACCGGGTATACAGATTGTTGACAAAGTGCGCTTTGTTTATGCCGGATGCGGCGTTAACGCCTTATCCg                                                                                       >  2:2534345/1‑100 (MQ=255)
                   ccGGATGGCGCTGCACGCAGAGGCACAGCGCCACCGGGTATACAGATTGTTGACAAAGTGCGCTTTGTTTATGCCGGATGCGGCGTTAACGCCTTATCCg                                                                                       >  1:208606/1‑100 (MQ=255)
                    cGGATGGCGCTGCACGCAGAGGCACAGCGCCACCGGGTATACAGATTGTTGACAAAGTGCGCTTTGTTTATGCCGGATGCGGCGTTAACGCCTTATCCgg                                                                                      <  2:2829031/100‑1 (MQ=255)
                     ggATGGCGCTGCACGCAGAGGCACAGCGCCACCGGGTATACAGATTGTTGACAAAGTGCGCTTTGTTTATGCCGGATAAGGCGTTAACGCCGCATCCGGc                                                                                     <  2:301978/100‑1 (MQ=37)
                      gATGGCGCTGCACGCAGAGGCACAGCGCCACCGGGTATACAGATTGTTGACAAAGTGCGCTTTGTTTATGCCGGATGCGGCGTTAACGCCTTATCCGGcc                                                                                    <  1:2378358/100‑1 (MQ=255)
                      gATGGCGCTGCACGCAGAGGCACAGCGCCACCGGGTATACAGATTGTTGACAAAGTGCGCTTTGTTTATGCCGGATGCGGCGTTAACGCCTTATCCGGcc                                                                                    <  2:33646/100‑1 (MQ=255)
                      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)
                           cgcTGCACGCAGAGGCACAGCGCCACCGGGTATACAGATTGTTGACAAAGTGCGCTTTGTCTATGCCGGATGCGGCGTTAACGCCTTATCCGGCCTAC‑aa                                                                              >  1:2507415/1‑100 (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)
                                cACGAAGAGGCACAGTGCCACCGTGTATACAGATTGTTGACAAAGTGCGCTTTGTTTATGCCGTATTTGGCGTTAACGCCTTATCCGGCCTAC‑AAAATCg                                                                         <  1:610574/100‑1 (MQ=18)
                                 aCGCAGAGGCACAGCGCCACCGGGTATACAGATTGTTGACAAAGTGCGCTTTGTTTATGCCGGATGCGGCGTTAACGCCTTATCCGGCCTAC‑AAAATCGt                                                                        >  2:1716547/1‑100 (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)
                                        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)
                                           acaGCGCCACCGGGTATACAGATTGTTGACAAAGTGCGCTTTGTTTATGCCGGATGCGGCGTTAACGCCTTATCCGGCCTAC‑AAAATCGTGCAAA‑TTCaa                                                             <  2:1760547/100‑1 (MQ=255)
                                            caGCGCCACCGGGTATACAGATTGTTGACAAAGTGCGCTTTGTTTATGCCGGATGCGGCGTTAACGCCTTATCCGGCCTAC‑AAAATCGTGCAAA‑TTCAAt                                                            <  2:2941306/100‑1 (MQ=255)
                                             aGCGCCACCGGGTATACAGATTGTTGACAAAGTGCGCTTTGTTTATGCCGGATGCGGCGTTAACGCCTTATCCGGCCTAC‑AAAATCGTGCAAA‑TTCAATa                                                           <  1:1077300/100‑1 (MQ=255)
                                             aGCGCCACCGGGTATACAGATTGTTGACAAAGTGCGCTTTGTTTATGCCGGATGCGGCGTTAACGCCTTATCCGGCCTAC‑AAAATCGTGCAAA‑TTCAATa                                                           <  1:2460594/100‑1 (MQ=255)
                                             aGCGCCACCGGGTATACAGATTGTTGACAAAGTGCGCTTTGTTTATGCCGGATGCGGCGTTAACGCCTTATCCGGCCTAC‑AAAATCGTGCAAA‑TTCAATa                                                           <  1:363991/100‑1 (MQ=255)
                                                gcCACCGGGTATACAGATTGTTGACAAAGTGCGCTTTGTTTATGCCGGATGCGGCGTTAACGCCTTATCCGGCCTAC‑AAAATCGTGCAAA‑TTCAATACAt                                                        >  2:1426741/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)
                                                     cGGGTATACAGATTGTTGACAAAGTGCGCTTTGTTTATGCCGGATGCGGCGTTAACGCCTTATCCGGCCTAC‑AAAATCGTGCAAA‑TTCAATACATTGCA‑g                                                  <  2:2135164/100‑1 (MQ=255)
                                                       ggTATACAGATTGTTGACAAAGTGCGCTTTGTTTATGCCGGATGCGGCGTTAACGCCTTATCCGGCCTAC‑AAAATCGTGCAAA‑TTCAATACATTGCA‑Ga‑a                                               >  2:1614804/1‑100 (MQ=255)
                                                         tataCAGATTGTTGACAAAGTGCGCTTTGTTTATGCCGGATGCGGCGTTAACGCCTTATCCGGCCTAC‑AAAATCGTGCAAA‑TTCAATACATTGCA‑GA‑At‑t                                            >  2:1294297/1‑100 (MQ=255)
                                                          ataCAGATTGTTGACAAAGTGCGCTTTGTTTATGCCGGATGCGGCGTTAACGCCTTATCCGGCCTAC‑AAAATCGTGCAAA‑TTCAATACATTGCA‑GA‑At‑tt                                           >  1:766061/1‑100 (MQ=255)
                                                            aCAGATTGTTGACAAAGTGCGCTTTGTTTATGCCGGATGCGGCGTTAACGCCTTATCCGGCCTAC‑AAAATCGTGCAAA‑TTCAATACATTGCA‑GA‑AT‑TTTc                                         <  1:1697648/100‑1 (MQ=255)
                                                             cAGATTGTTGACAAAGTGCGCTTTGTTTATGCCGGATGCGGCGTTAACGCCTTATCCGGCCTAC‑AAAATCGTGCAAA‑TTCAATACATTGCA‑GA‑AT‑TTTCg                                        <  2:1088951/100‑1 (MQ=255)
                                                             cAGATTGTTGACAAAGTGCGCTTTGTTTATGCCGGATGCGGCGTTAACGCCTTATCCGGCCTAC‑AAAATCGTGCAAA‑TTCAATACATTGCA‑GA‑AT‑TTTCg                                        <  1:429413/100‑1 (MQ=255)
                                                             cAGATTGCTGACAACGTGCGCGTTGTTCATGCCGGATGCGGCGTGAACGCCTTATCCGGCCTAC‑AAAGTCGTATAAACTGCAATACATTgaatttatacg                                            <  2:709753/100‑11 (MQ=14)
                                                              aGATTGTTGACAAAGTGCGCTTTGTTTATGCCGGATGCGGCGTTAACGCCTTATCCGGCCTAC‑AAAATCGTGCAAA‑TTCAATACATTGCA‑GA‑AT‑TTTCGt                                       <  2:953252/100‑1 (MQ=255)
                                                                 ttgttgACAAAGTGCGCTTTGTTTATGCCGGATGCGGCGTTAACGCCTTATCCGGCCTAC‑AAAATCGTGCAAA‑TTCAATACATTGCA‑GA‑AT‑TTTCGTAgg                                    <  1:701843/100‑1 (MQ=255)
                                                                 gcataaaCAAAGCGCACTTTGTTTATGCCGGATGCGGCGTTAACGCCTTATCCGGCCTAC‑AAAATCGTGCAAA‑TTCAATACATTGCA‑GA‑AT‑TTTCGTAgg                                    >  1:404911/7‑100 (MQ=18)
                                                                   gttgACAAAGTGCGCTTTGTTTATGCCGGATGCGGCGTTAACGCCTTATCCGGCCTAC‑AAAATCGTGCAAA‑TTCAATACATTGCA‑GA‑AT‑TTTCGTAGGcc                                  <  2:571259/100‑1 (MQ=255)
                                                                    ttgACAAAGTGCGCTTTGTTTATGCCGGATGCGGCGTTAACGCCTTATCCGGCCTAC‑AAAATCGTGCAAA‑TTCAATACATTGCA‑GA‑AT‑TTTCGTAGGCCt                                 <  1:2063062/100‑1 (MQ=255)
                                                                    ttgACAAAGTGCGCTTTGTTTATGCCGGATGCGGCGTTAACGCCTTATCCGGCCTAC‑AAAATCGTGCAAA‑TTCAATACATTGCA‑GA‑AT‑TTTCGTAGGCCt                                 >  1:23980/1‑100 (MQ=255)
                                                                     tgACAAAGTGCGCTTTGTTTATGCCGGATGCGGCGTTAACGCCTTATCCGGCCTAC‑AAAATCGTGCAAA‑TTCAATACATTGCA‑GA‑AT‑TTTCGTAGGCCTg                                >  1:2301212/1‑100 (MQ=255)
                                                                     ngACAAAGTGCGCTTTGTTTATGCCGGATGCGGCGTTAACGCCTTATCCGGCCTAC‑AAAATCGTGCAAA‑TTCAATACATTGCA‑GA‑AT‑TTTCGTAGGCCTg                                >  1:1416603/2‑100 (MQ=255)
                                                                      gACAAAGTGCGCTTTGTTTATGCCGGATGCGGCGTTAACGCCTTATCCGGCCTAC‑AAACTCGTGCAAA‑TTCAATACATTGCA‑GA‑AT‑TTTCGTAGGCCTGa                               <  1:1644745/100‑1 (MQ=255)
                                                                            gTGCGCTTTGTTTATGCCGGATGCGGCGTTAACGCCTTATCCGGCCTAC‑AAAATCGTGCAAA‑TTCAATACATTGCA‑GA‑AT‑TTTCGTAGGCCTGATAAGCg                         >  2:1733813/1‑100 (MQ=255)
                                                                              gcgcTTTGTTTATGCCGGATGCGGCGTTAACGCCTTATCCGGCCTAC‑AAAATCGTGCAAA‑TTCAATACATTGCA‑GA‑AT‑TTTCGTAGGCCTGATAAGCGGa                       >  2:1422895/1‑100 (MQ=255)
                                                                               cgcTTTGTTTATGCCGGATGCGGCGTTAACGCCTTATCCGGCCTAC‑AAAATCGTGCAAA‑TTCGATACATTGCA‑GA‑AT‑TTTCGTAGGCCTGATAAGCGGAg                      <  1:659977/100‑1 (MQ=255)
                                                                                gcTTTGTTTATGCCGGATGCGGCGTTAACGCCTTATCCGGCCTAC‑AAAATCGTGCAAA‑TTCAATACATTGCA‑GA‑AT‑TTTCGTAGGCCTGATAAGCGGAgc                     <  1:252980/100‑1 (MQ=255)
                                                                                gcTTTGTTTATGCCGGATGCGGCGTTAACGCCTTATCCGGCCTAC‑AAAATCGTGCAAA‑TTCAATACATTGCA‑GA‑AT‑TTTCGTAGGCCTGATAAGCGGAgc                     <  2:208607/100‑1 (MQ=255)
                                                                                 cTTTGTTTATGCCGGATGCGGCGTTAACGCCTTATCCGGCCTAC‑AAAATCGTGCAAA‑TTCAATACATTGCA‑GA‑AT‑TTTCGTAGGCCTGATAAGCGGAgcg                    >  2:763133/1‑100 (MQ=255)
                                                                                 cTTTGTTTATGCCGGATGCGGCGTTAACGCCTTATCCGGCCTAC‑AAAATCGTGCAAA‑TTCAATACATTGCA‑GA‑AT‑TTTCGTAGGCCTGATAAGCGGAgcg                    <  1:1821906/100‑1 (MQ=255)
                                                                                  cttttgtaggCCGGATAAGGCGTTTACGCCGCATCCGGCCTAC‑AAAAGCATGCAAA‑TTCAATATATTGCA‑GA‑GT‑TCACGTAGGCCTGATAAGCGTAgcgc                   >  2:695212/10‑100 (MQ=17)
                                                                                     ntttATGCCGGATGCGGCGTTAACGCCTTATCCGGCCTAC‑AAAATCGTGCAAA‑TTCAATACATTGCA‑GA‑AT‑TTTCGTAGGCCTGATAAGCGGAGCGCATc                >  2:1638267/2‑100 (MQ=255)
                                                                                      tgtaggCCGGATAAGGCGTTTACGCCGCATCCGGCCTAC‑AAATTC‑TGCAAA‑TTCAACAAATTGC‑‑GA‑ATCCCTTGTAGGCCTGATAAGCGTAGCGCATCAg              >  2:2406988/6‑100 (MQ=18)
                                                                                       gtaggCCGGATAAGGCGTTCACGCCGCATCCGGCCTAC‑AAGTTCGTGCAAA‑TTCAATAAATTGCA‑‑A‑TA‑TGACGTAGGCCTGATAAGCGTAGCGCATCAgg             >  1:1945164/5‑100 (MQ=18)
                                                                                        taggCCGGATAAGGCGTTCACGCCGCATCCGGCCTAC‑AAAATTGTGTAAA‑TTCAATAAATTGAT‑CG‑AT‑CCACGTAGGCCTGATAAGCGTAGCGCATCAgg             >  2:373817/4‑100 (MQ=21)
                                                                                        taggCCGGATAAGGCGTTCACGCCGCATCCGGCCTAC‑AAAACCTTGCTAA‑TTCAATATATTGCAGGG‑AC‑TAT‑GTAGGCCTGATAAGCATAGCGCATCAgg             >  2:2350227/4‑100 (MQ=18)
                                                                                         aggCCGGATAAGGCGTTCACTCCGCATCCGGCCTAC‑AAAAGCACGCAAA‑TTCAATATATTGCA‑GA‑GA‑TTGCGTAGGCCTGATAAGCGTAGCGCATCAGGc            >  1:1820111/3‑100 (MQ=11)
                                                                                         aggCCGGATAAGGCGTTCACGCCGCATCCGGCCTAC‑AAAAGCACGCAAA‑CTCAATATATTGCA‑GAGAT‑CAT‑GTAGGCCTGATAAGCGTAGCGCATCAGGc            >  2:2804013/3‑100 (MQ=17)
                                                                                         aggCCGGATAAGGCGTTAACGCCGCATCCGGCCTAC‑AAAATCGTGCAAA‑TTCAATACATTGCA‑GA‑AT‑TTTCGTAGGCCTGATAAGCGGAGCGCATCAGGc            >  2:2463611/3‑100 (MQ=18)
                                                                                          tGCCGGATGCGGCGTTAACGCCTTATCCGGCCTAC‑AAAATCGTGCAAA‑TTCAATACATTGCA‑GA‑AT‑TTTCGTAGTCCTGATAAGCGTAGCGCATCAGGCa           >  2:438892/1‑100 (MQ=255)
                                                                                          ngCCGGATAAGGCGTTCACGCCGCATCCGGCCTAC‑AAAACTGTGCAAA‑TTCAATATATTGCA‑GG‑AA‑ACACGTAGGCCTGATAAGCGAAGCGCATCAGGCa           >  2:2645807/2‑100 (MQ=18)
                                                                                          ggCCGGATAAGGCGTTTACGCCGCATCCGGCCTAC‑AAAAGCACGCAAA‑CTCAATATATTGCA‑GA‑GA‑TCATGTAGGCCTGATAAGCGTAGCGCATCAGGCa           >  1:2921783/2‑100 (MQ=17)
                                                                                          ggCCGGATAAGGCGTTCACGCCGCATCCGGCCTAC‑AAAGTCGTATAAA‑TTCAATGTATTGCA‑GT‑TT‑ATTGGTAGGCCTGATAAGCGCAGCGCATCAGGCa           >  1:2889630/2‑100 (MQ=21)
                                                                                          ggCCGGATAAGGCGTTCACGCCGCATCCGGCCTAC‑AAAAGCATGAAAA‑TTCAATATATTGCA‑GG‑AG‑CTGCGTAGGCCTGATAAGCGTAGCGCATCAGGCa           >  2:1960166/2‑100 (MQ=18)
                                                                                          ggCCGGATAAGGCGTTCACGCCGCATCCGGCCTAC‑AAAACTGTGCAAA‑TTCAATATATTGCA‑GG‑AA‑ACACGTAGGCCTGATAAGCGAAGCGCATCAGGCa           >  2:545733/2‑100 (MQ=18)
                                                                                           gCCGGATGCGGCGTTAACGCCTTATCCGGCCTAC‑AAAATCGTGCAAA‑TTCAATACATTGGC‑GA‑AT‑TTTCGTAGGCCTGATCAGCAGAGCGCATCAAGcct          >  1:283062/1‑98 (MQ=17)
                                                                                           gCCGGATGCGGCGTTAACGCCTTATCCGGCCTAC‑AAAATCGTGCAAA‑TTCAATACATTGCA‑GA‑AT‑TTTCGTAGGCCTGATAAGCGGAGCGCATCAGGCAt          >  2:285078/1‑100 (MQ=255)
                                                                                           gCCGGATAAGGCGTTTACGCCGCATCCGGCCTAC‑AAAAGCATGCAAA‑TTCAATATATTGCA‑GA‑GT‑TCACGTAGGCCTGATAAGCGTAGCGCATCAGGCag          >  2:1900194/1‑99 (MQ=17)
                                                                                           gCCGGATAAGGCGTTCACGCCGCATCCGGCCTACAAAAATCTTGCCAA‑TTCAATATATTGCA‑GG‑AC‑CAT‑GTAGGCCTGATAAGCGTAGCGCATCAGGCag          >  2:944004/1‑99 (MQ=18)
                                                                                           gCCGGATAAGGCGTTCACGCCGCATCCGGCCTAC‑AAAGTCGTATAAC‑TTCAATGTATTGCA‑GT‑TT‑ATTGGTAGGCCTGATAAGCGCAGCGCATCAGGCaa          >  1:1369201/1‑99 (MQ=18)
                                                                                           gCCGGATAAGGCGTTCACGCCGCATCCGGCCTAC‑AAAGTCGTATAAA‑TTCAATGTATTGCA‑GT‑TT‑ATTGGTAGGCCTGATAAGCGCAGCGCATCAGGCaa          >  2:2669950/1‑99 (MQ=25)
                                                                                           gCCGGATAAGGCGTTCACGCCGCATCCGGCCTAC‑AAAAACCTGCAAA‑TTCAATATATTGCA‑GG‑AT‑ACGTGTAGGCCTGATAAGCGTAGCGCATCAGGCaa          >  2:2385627/1‑99 (MQ=18)
                                                                                           gCCGGATAAGGCGTTCACGCCGCATCCGGCCTAC‑AAAAACATGCAAA‑TTCAATAGATTGCA‑GAGAT‑TAT‑GTAGGCCTGATAAGCGTAGCGCATCAGGCaa          >  1:2236509/1‑99 (MQ=17)
                                                                                           gCCGGATAAGGCGTTCACGCCGCATCCGGCCTAC‑‑AAATCATGCAAA‑TCCAAAATATTGCA‑GG‑AG‑TTATGTAGGCCTGATAAGCGTAGCGCATCAGGCAtt         >  1:990152/1‑100 (MQ=21)
                                                                                                 tGCGGCGTTAACGCCTTATCCGGCCTAC‑AAAATCGTGCAAA‑TTCAATACATTGCA‑GA‑AT‑TTTCGTAGGCCTGATAAGCGGAGCGCATCAGGCATTtttgc    <  2:1510431/100‑1 (MQ=255)
                                                                                                   |                                                                                                        
AGGGAATGAAGTATTGCTGCCGGATGGCGCTGCACGCAGAGGCACAGCGCCACCGGGTATACAGATTGTTGACAAAGTGCGCTTTGTTTATGCCGGATGCGGCGTTAACGCCTTATCCGGCCTAC‑AAAATCGTGCAAA‑TTCAATACATTGCA‑GA‑AT‑TTTCGTAGGCCTGATAAGCGGAGCGCATCAGGCATTTTTGCCT  >  NC_000913/3705872‑3706070

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: ‑