Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009133,705,9710CA35.1% 378.4 / 246.7 211intergenic (‑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 (74/63);  new base A (63/11);  total (137/74)
Fisher's exact test for biased strand distribution p-value = 4.44e-06
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.49e-01

Only 100 of 225 total aligned reads displayed.

GGGAATGAAGTATTGCTGCCGGATGGCGCTGCACGCAGAGGCACAGCGCCACCGGGTATACAGATTGTTGACAAAGTGCGCTTTGTTTATGCCGGATGCGGCGTTAACGCCTTATCCGGCCTACAAAATCGTGCAAATTCAATACATTGCAGAATT‑‑‑TTCGTAGGCCTGATAAGCGGAGCGCATCAGGCATTTTTGCCT  >  NC_000913/3705873‑3706070
                                                                                                  |                                                                                                          
gggAATGAAGTATTGCTGCCGGATGGCGCTGCACGCAGAGGCACAGCGCCACCGGGTATACAGATTGTTGACAAAGTGCGCTTTGTTTATGCCGGATgcg                                                                                                           >  1:2654886/1‑100 (MQ=255)
 ggAATGAAGTATTGCTGCCGGATGGCGCTGCACGCAGAGGCACAGCGCCACCGGGTATACAGATTGTTGACAAAGTGCGCTTTGTTTATGCCGGATgcgg                                                                                                          <  1:826429/100‑1 (MQ=255)
      gAAGTATTGCTGCCGGATGGCGCTGCACGCAGAGGCACAGCGCCACCGGGTATACAGATTGTTGACAAAGTGCGCTTTGTTTATGCCGGATGCGGCGTTa                                                                                                     <  1:1580793/100‑1 (MQ=255)
       aaGTATTGCTGCCGGATGGCGCTGCACGCAGAGGCACAGCGCCACCGGGTATACAGATTGTTGACAAAGTGCGCTTTGTTTATGCCGGATGCGGCGTTaa                                                                                                    >  1:203148/1‑100 (MQ=255)
       aaGTATTGCTGCCGGATGGCGCTGCACGCAGAGGCACAGCGCCACCGGGTATACAGATTGTTGACAAAGTGCGCTTTGTTTATGCCGGATGCGGCGTTaa                                                                                                    >  2:48851/1‑100 (MQ=255)
        aGTATTGCTGCCGGATGGCGCTGCACGCAGAGGCACAGCGCCACCGGGTATACAGATTGTTGACAAAGTGCGCTTTGTTTATGCCGGATGCGGCGTTAAc                                                                                                   <  2:2883559/100‑1 (MQ=255)
        aGTATTGCTGCCGGATGGCGCTGCACGCAGAGGCACAGCGCCACCGGGTATACAGATTGTTGACAAAGTGCGCTTTGTTTATGCCGGATGCGGCGTTAAc                                                                                                   >  1:224525/1‑100 (MQ=255)
         gTATTGCTGCCGGATGGCGCTGCACGCAGAGGCACAGCGCCACCGGGTATACAGATTGTTGACAAAGTGCGCTTTGTTTATGCCGGATGCGGCGTTAAc                                                                                                   >  1:1824749/1‑99 (MQ=255)
         gTATTGCTGCCGGATGGCGCTGCACGCAGAGGCACAGCGCCACCGGGTATACAGATTGTTGACAAAGTGCGCTTTGTTTATGCCGGATGCGGCGTTAACg                                                                                                  <  1:1253389/100‑1 (MQ=255)
          tATTGCTGCCGGATGGCGCTGCACGCAGAGGCACAGCGCCACCGGGTATACAGATTGTTGACAAAGTGCGCTTTGTTTATGCCGGATGCGGCGTTAACGc                                                                                                 <  2:1651575/100‑1 (MQ=255)
          tATTGCTGCCGGATGGCGCTGCACGCAGAGGCACAGCGCCACCGGGTATACAGATTGTTGACAAAGTGCGCTTTGTTTATGCCGGATGCGGCGTTAACGc                                                                                                 <  2:10896/100‑1 (MQ=255)
           aTTGCTGCCGGATGGCGCTGCACGCAGAGGCACAGCGCCACCGGGTATACAGATTGTTGACAAAGTGCGCTTTGTTTATGCCGGATGCGGCGTTAACGcc                                                                                                <  2:1659285/100‑1 (MQ=255)
           aTTGCTGCCGGATGGCGCTGCACGCAGAGGCACAGCGCCACCGGGTATACAGATTGTTGACAAAGTGCGCTTTGTTTATGCCGGATGCGGCGTTAACGcc                                                                                                <  2:731255/100‑1 (MQ=255)
            ttGCTGCCGGATGGCGCTGCACGCAGAGGCACAGCGCCACCGGGTATACAGATTGTTGACAAAGTGCGCTTTGTTTATGCCGGATGCGGCGTTAACGCCt                                                                                               <  2:2466725/100‑1 (MQ=255)
            ttGCTGCCGGATGGCGCTGCACGCAGAGGCACAGCGCCACCGGGTATACAGATTGTTGACAAAGTGCGCTTTGTTTATGCCGGATGCGGCGTTAACGCCt                                                                                               >  1:1087709/1‑100 (MQ=255)
                 gcCGGATGGCGCTGCACGCAGAGGCACAGCGCCACCGGGTATACAGATTGTTGACAAAGTGCGCTTTGTTTATGCCGGATGCGGCGTTAACGCCTTATcc                                                                                          >  1:15593/1‑100 (MQ=255)
                   cGGATGGCGCTGCACGCAGAGGCACAGCGCCACCGGGTATACAGATTGTTGACAAAGTGCGCTTTGTTTATGCCGGATGCGGCGTTAACGCCTTATCCgg                                                                                        >  1:101277/1‑100 (MQ=255)
                   cGGATGGCGCTGCACGCAGAGGCACAGCGCCACCGGGTATACAGATTGTTGACAAAGTGCGCTTTGTTTATGCCGGATGCGGCGTTAACGCCTTATCCgg                                                                                        >  2:360503/1‑100 (MQ=255)
                     gATGGCGCTGCACGCAGAGGCACAGCGCCACCGGGTATACAGATTGTTGACAAAGTGCGCTTTGTTTATGCCGGATGCGGCGTTAACGCCTTATCCGGcc                                                                                      <  2:1002249/100‑1 (MQ=255)
                     gATGGCGCTGCACGCAGAGGCACAGCGCCACCGGGTATACAGATTGTTGACAAAGTGCGCTTTGTTTATGCCGGATGCGGCGTTAACGCCTTATCCGGcc                                                                                      <  1:377218/100‑1 (MQ=255)
                     gATGGCGCTGCACGCAGAGGCACAGCGCCACCGGGTATACAGATTGTTGACAAAGTGCGCTTTGTTTATGCCGGATGCGGCGTTAACGCCTTATCCGGcc                                                                                      <  2:1139664/100‑1 (MQ=255)
                     gATGGCGCTGCACGCAGAGGCACAGCGCCACCGGGTATACAGATTGTTGACAAAGTGCGCTTTGTTTATGCCGGATAAGGCGTTAACGCCGCATCCGGca                                                                                      <  1:1807319/100‑2 (MQ=37)
                      aTGGCGCTGCACGCAGAGGCACAGCGCCACCGGGTATACAGATTGTTGACAAAGTGCGCTTTGTTTAGGCCGGATAAGGCGTTAACGCCGCATCCGGcat                                                                                     <  2:3056964/100‑3 (MQ=25)
                        ggCGCTGCACGCAGAGGCACAGCGCCACCGGGTATACAGATTGTTGACAAAGTGCGCTTTGTTTATGCCGGATAAGGCGTTAACGCCGCATCCGGCATaa                                                                                   <  1:1845536/100‑2 (MQ=25)
                         gcgcTGCACGCAGAGGCACAGCGCCACCGGGTATACAGATTGTTGACAAAGTGCGCTTTGTTTATGCCGGATGCGGCGTTAACGCCTTATCCGGCCTACa                                                                                 <  2:210816/100‑1 (MQ=255)
                         gcgcTGCACGCAGAGGCACAGCGCCACCGGGTATACAGATTGTTGACAAAGTGCGCTTTGTTTATGCCGGATAAGGCGTTAACGCCGCATCCGGCATaaa                                                                                  <  1:1461271/100‑3 (MQ=25)
                          cgcTGCACGCAGAGGCACAGCGCCACCGGGTATACAGATTGTTGACAAAGTGCGCTTTGTTTATGCCGGATGCGGCGTTAACGCCTTATCCGGCCTACaa                                                                                >  1:1930134/1‑100 (MQ=255)
                          cgcTGCACGCAGAGGCACAGCGCCACCGGGTATACAGATTGTTGACAAAGTGCGCTTTGTTTATGCCGGATAAGGCGTTAACGCCGCATCCGGcataaac                                                                                 <  2:2290683/100‑7 (MQ=25)
                           gcTGCACGCAGAGGCACAGCGCCACCGTGTATACAGATTGTTGACAAAGTGCGCTTTGTTTATGCCGGATGCGGCGTTAACGCCTTATCCGGCCTACaaa                                                                               >  2:532805/1‑100 (MQ=255)
                             tctaCGCAGAGGCACAGCGCCACCGGGTATACAGATTGTTGACAAAGTGCGCTTTGTTTATGCCGGATGCGGCGTTAACGCCTTATCCGGCCTACAAAAt                                                                             <  1:1624342/97‑1 (MQ=255)
                             tGCACGCAGAGGCACAGCGCCACCGGGTATACAGATTGTTGACAAAGTGCGCTTTGTTTATGCCGGATGCGGCGTTAACGCCTTATCCGGCCTACAAAAt                                                                             >  1:729796/1‑100 (MQ=255)
                                aCGCAGAGGCACAGCGCCACCGGGTATACAGATTGTTGACAAAGTGCGCTTTGTTTATGCCGGATGCGGCGTTAACGCCTTATCCGGCCTACAAAAtaga                                                                          >  2:1624338/1‑97 (MQ=255)
                                aCGCAGAGGCACAGCGCCACCGGGTATACAGATTGTTGACAAAGTGCGCTTTGTTTATGCCGGATGCGGCGTTAACGCCTTATCCGGCCTACAAAATCGt                                                                          >  2:1577149/1‑100 (MQ=255)
                                 cGCAGAGGCACAGCGCCACCGGGTATACAGATTGTTGACAAAGTGCGCTTTGTTTATGCCGGATGCGGCGTTAACGCCTTATCCGGCCTACAAAATCGTg                                                                         <  2:1521408/100‑1 (MQ=255)
                                  gCAGAGGCACAGCGCCACCGGGTATACAGATTGTTGACAAAGTGCGCTTTGTTTATGCCGGATGCGGCGTTAACGCCTTATCCGGCCTACAAAATCGTGc                                                                        <  2:2730701/100‑1 (MQ=255)
                                       ggCACAGCGCCACCGGGTATACAGATTGTTGACAAAGTGCGCTTTGTTTATGCCGGATGCGGCGTTAACGCCTTATCCGGCCTACAAAATCGTGCAAAtt                                                                   >  2:3011433/1‑100 (MQ=255)
                                       ggCACAGCGCCACCGGGTATACAGATTGTTGACAAAGTGCGCTTTGTTTATGCCGGATGCGGCGTTAACGCCTTATCCGGCCTACAAAATCGTGCAAAtt                                                                   <  1:480460/100‑1 (MQ=255)
                                       ggCACAGCGCCACCGGGTATACAGATTGTTGACAAAGTGCGCTTTGTTTATGCCGGATGCGGCGTTAACGCCTTATCCGGCCTACAAAATCGTGCAAAtt                                                                   >  1:1948675/1‑100 (MQ=255)
                                        gCACAGCGCCACCGGGTATACAGATTGTTGACAAAGTGCGCTTTGTTTATGCCGGATGCGGCGTTAACGCCTTATCCGGCCTACAAAATCGTGCAAATTc                                                                  <  2:1427971/100‑1 (MQ=255)
                                        gCACAGCGCCACCGGGTATACAGATTGTTGACAAAGTGCGCTTTGTTTATGCCGGATGCGGCGTTAACGCCTTATCCGGCCTACAAAATCGTGCAAATTc                                                                  <  2:750791/100‑1 (MQ=255)
                                        gCACAGCGCCACCGGGTATACAGATTGTTGACAAAGTGCGCCTTGTTTATGCCGGATGCGGCGTTAACGCCTTATCCGGCCTACAAAATCGTGCAAATTc                                                                  <  1:216976/100‑1 (MQ=255)
                                         cacaGCGCCACCGGGTATACAGATTGTTGACAAAGTGAGCTTTGTTTATGCCGGATGCGGCGTTAACGCCTTATCCGGCCTACAAAATCGTGCAAATTCa                                                                 >  2:1798117/1‑100 (MQ=255)
                                           caGCGCCACCGGGTATACAGATTGTTGACAAAGTGCGCTTTGTTTATGCCGGATGCGGCGTTAACGCCTTATCCGGCCTACAAAATCGTGCAAATTCAAt                                                               <  2:1059568/100‑1 (MQ=255)
                                           caGCGCCACCGGGTATACAGATTGTTGACAAAGTGCGCTTTGTTTATGCCGGATGCGGCGTTAACGCCTTATCCGGCCTACAAAATCGTGCAAATTCAAt                                                               >  1:1917666/1‑100 (MQ=255)
                                               gcCACCGGGTATACAGATTGTTGACAAAGTGCGCTTTGTTTATGCCGGATGCGGCGTTAACGCCTTATCCGGCCTACAAAATCGTGCAAATTCAATACAt                                                           <  2:1494721/100‑1 (MQ=255)
                                                ccACCGGGTATACAGATTGTTGACAAAGTGCGCTTTGTTTATGCCGGATGCGGCGTTAACGCCTTATCCGGCCTACAAAATCGTGCAAATTCAATACAtt                                                          >  2:181037/1‑100 (MQ=255)
                                                ccACCGGGTATACAGATTGTTGACAAAGTGCGCTTTGTTTATGCCGGATGCGGCGTTAACGCCTTATCCGGCCTACAAAATCGTGCAAATTCAATACAtt                                                          <  2:118940/100‑1 (MQ=255)
                                                ccACCGGGTATACAGATTGTTGACAAAGTGCGCTTTGTTTATGCCGGATGCGGCGTTAACGCCTTATCCGGCCTACAAAATCGTGCAAATTCAATACAtt                                                          >  2:3303309/1‑100 (MQ=255)
                                                 cACCGGGTATACAGATTGTTGACAAAGTGCGCTTTGTTTATGCCGGATGCGGCGTTAACGCCTTATCCGGCCTACAAAATCGTGCAAATTCAATACATTg                                                         >  2:1264582/1‑100 (MQ=255)
                                                  aCCGGGTATACAGATTGTTGACAAAGTGCGCTTTGTTTATGCCGGATGCGGCGTTAACGCCTTATCCGGCCTACAAAATCGTGCAAATTCAATACATTGc                                                        <  1:1762730/100‑1 (MQ=255)
                                                  aCCGGGTATACAGATTGTTGACAAAGTGCGCTTTGTTTATGCCGGATGCGGCGTTAACGCCTTATCCGGCCTACAAAATCGTGCAAATTCAATACATTGc                                                        >  1:1597698/1‑100 (MQ=255)
                                                   ccGGGTATACAGATTGTTGACAAAGTGCGCTTTGTTTATGCCGGATGCGGCGTTAACGCCTTATCCGGCCTACAAAATCGTGCAAATTCAATACATTGCa                                                       <  1:2274060/100‑1 (MQ=255)
                                                    cGGGTATACAGATTGTTGACAAAGTGCGCTTTGTTTATGCCGGATGCGGCGTTAACGCCTTATCCGGCCTACAAAATCGTGCAAATTCAATACATTGCAg                                                     <  1:2392299/100‑1 (MQ=255)
                                                         ataCAGATTGTTGACAAAGTGCGCTTTGTTTATGCCGGATGCGGCGTTAACGCCTTATCCGGCCTACAAAATCGTGCAAATTCAATACATTGCAGAAtt‑‑‑t                                           <  1:2149524/100‑1 (MQ=255)
                                                         ataCAGATTGTTGACAAAGTGCGCTTTGTTTATGCCGGATGCGGCGTTAACGCCTTATCCGGCCTACAAAATCGTGCAAATTCAATACATTGCAGAAtt‑‑‑t                                           >  1:3352256/1‑100 (MQ=255)
                                                              gATTGTTGGCAAAGTGCGCTTTGTTTATGCCGGATGCGGCGTTAACGCCTTATCCGGCCTACAAAATCGTGCAAATTCAATACATTGCAGAATT‑‑‑TTCGTa                                      <  2:2649624/100‑1 (MQ=255)
                                                              gATTGTTGACAAAGTGCGCTTTGTTTATGCCGGATGCGGCGTTAACGCCTTATCCGGCCTACAAAATCGTGCAAATTCAATACATTGCAGAATT‑‑‑TTCGt                                       <  2:874430/99‑1 (MQ=255)
                                                              gATTGTTGACAAAGTGCGCTTTGTTTATGCCGGATGCGGCGTTAACGCCTTATCCGGCCTACAAAATCGTGCAAATTCAATACATTGCAGAATT‑‑‑TTCGTa                                      >  1:1043269/1‑100 (MQ=255)
                                                                  gttgACAAAGTGCGCTTTGTTTATGCCGGATGCGGCGTTAACGCCTTATCCGGCCTACAAAATCGTGCAAATTCAATACATTGCAGAATT‑‑‑TTCGTAGGcc                                  <  2:1650893/100‑1 (MQ=255)
                                                                   ttgACAAAGTGCGCTTTGTTTATGCCGGATGCGGCGTTAACGCCTTATCCGGCCTACAAAATCGTGCAAATTCAATACATTGCAGAATT‑‑‑TTCGTAGGCCt                                 >  1:2843728/1‑100 (MQ=255)
                                                                   taaaCAAAGCGCACTTTGTTTATGCCGGATGCGGCGTTAACGCCTTATCCGGCCTACAAAATCGTGCAAATTCAATACATTGCAGAATT‑‑‑TTCGTAGGCCt                                 >  1:1291537/4‑100 (MQ=18)
                                                                     gACAAAGTGCGCTTTGTTTATGCCGGATGCGGCGTTAACGCCTTATCCGGCCTACAAAATCGTGCAAATTCAATACATTGCAGAATT‑‑‑TTCGTAGGCCTGa                               >  1:58313/1‑100 (MQ=255)
                                                                      aCAAAGTGCGCTTTGTTTATGCCGGATGCGGCGTTAACGCCTTATCCGGCCTACAAAATCGTGCAAATTCAATACATTGCAGAATT‑‑‑TTCGTAGGCCTGAt                              <  2:2765929/100‑1 (MQ=255)
                                                                       cAAAGTGCGCTTTGTTTATGCCGGATGCGGCGTTAACGCCTTATCCGGCCTACAAAATCGTGCAAATTCAATACATTGCAGAATT‑‑‑TTCGTAGGCCTGATa                             >  2:282577/1‑100 (MQ=255)
                                                                                 attttgtaggCCGGATAAGGCGTTTACGCCGCATCCGGCCTACAAAATCGTGCTAATTCAATATATTGCAGAAAC‑‑‑CTTGTAGGCCTGATAAGCGTAgcgc                   >  1:1539805/10‑100 (MQ=18)
                                                                                   tttgtaggCCGGATAAGGCGTTCACGCCGCATCCGGCATACAAATCATGCAAATCCAAAATATTGCAGGAGT‑‑‑TATGTAGGCCTGATAAGCGTAGCGCATc                >  1:669997/8‑100 (MQ=18)
                                                                                   tGTTTATGCCGGATGCGGCGTTAACGCCTTATCCGGCCTACAAAATCGTGCAAATTCAATACATTGCAGAATT‑‑‑TTCGTAGGCCTGATAAGCGGAGCGCAt                 >  1:1736916/1‑100 (MQ=255)
                                                                                    gTTTATGCCGGATGCGGCGTTAACGCCTTATCCGGCCTACAAAATCGTGCAAATTCAATACATTGCAGAATT‑‑‑TTCGTAGGCCTGATAAGCGGAGCGCATc                >  2:1326802/1‑100 (MQ=255)
                                                                                     tgtaggCCGGATAAGGCGTTTACGCCGCATCCGGCCTACAAAATCGTGCTAATTCAATATATTGCAGAAAC‑‑‑CTTGTAGGCCTGATAAGCGTAGCGCATCa               >  2:2585220/6‑100 (MQ=17)
                                                                                     tgtaggCCGGATAAGGCGTTCACGCCGCATCCGGCCTACAAAAATCTTGCCAATTCAATATATTGCAGGACC‑‑‑ATGTAGGCCTGATAAGCGTAGCGCATCa               >  2:2406714/6‑100 (MQ=18)
                                                                                     tgtaggCCGGATAAGGCGTTCACGCCGCATCCGGCCTACAAAATCATGAAAATTCAATACATTGCAAGATT‑‑‑TTCGTAGGCCTGATAAGCGTAGCGCATCa               >  2:529276/6‑100 (MQ=17)
                                                                                      ttATGCCGGATGCGGCGTTAACGCCTTATCCGGCCTACAAAATCGTGCAAATTCAATACATTGCAGAATT‑‑‑TTCGTAGGCCTGATAAGCGGAGCGCATCAg              >  1:1336780/1‑100 (MQ=255)
                                                                                      gtaggCCGGATAAGGCGTTTACGCCGCATCCGGCCTACAAAAGCGTGCAAATTCAATACATTGCATGGGCC‑‑‑ATGTAGGCCTGATAAGCGTAGCGCATCAg              >  1:2763613/5‑100 (MQ=17)
                                                                                      gtaggCCGGATAAGGCGTTCACGCCGCATCCGGCCTACAAAAATCTTGCCAATTCAATATATTGCAGGACC‑‑‑ATGTAGGCCTGATAAGCGTAGCGCATCAg              >  2:2228552/5‑100 (MQ=18)
                                                                                      gtaggCCGGATAAGGCGTTCACGCCGCATCCGGCATACAAAAGCATGAAAATTCAATATATTGCAGGAGC‑‑‑TGCGTAGGCCTGATAAGCGTAGCGCATCAg              >  1:2292466/5‑100 (MQ=18)
                                                                                       taggCCGGATAAGGCGTTCACGCCGCATCCGGCCTACAAGATCGTGCAAATTCAATATATTGCATGAGC‑‑‑CACGTAGGCCTGATAAGCGTAGCGCATCAgg             >  2:1840102/4‑100 (MQ=17)
                                                                                       taggCCGGATAAGGCGTTCACGCCGCATCCGGCATACAAAGTCGTATAAATTCAATGTATTGCAGTTTA‑‑‑TTGGTAGGCCTGATAAGCGCAGCGCATCAgg             >  1:3230938/4‑100 (MQ=21)
                                                                                        aggCCGGATAAGGCGTTTACGCCGCATCCGGCCTACAAAATCGTGCTAATTCAATATATTGCAGAAAC‑‑‑CTTGTAGGCCTGATAAGCGTAGCGCATCAGGc            >  1:3119162/3‑100 (MQ=18)
                                                                                        aggCCGGATAAGGCGTTTACGCCGCATCCGGCCTACAAAAATGTGCAAATTCAATAAATTGCAATTCAACTTGTAGGCCTGATAAGCGCAGCGCATCAgg             >  1:2498707/3‑100 (MQ=21)
                                                                                        aggCCGGATAAGGCGTTCACGCCGCATCCGGCCTACAAAAATCTTGCCAATTCAATATATTGCAGGACC‑‑‑ATGTAGGCCTGATAAGCGTAGCGCATCAGGc            >  1:50031/3‑100 (MQ=18)
                                                                                        aggCCGGATAAGGCGTTCACGCCGCATCCGGCCTACAAAAGCATGAAAATTCAATATATTGCAGGAGC‑‑‑TGCGTAGGCCTGATAAGCGTAGCGCATCAGGc            >  2:1592886/3‑100 (MQ=18)
                                                                                         tGCCGGATGCGGCGTTAACGCCTTATCCGGCCTACAAAATCGTGCAAATTCAATACATTGCAGAATT‑‑‑TTCGTAGGCCTGATAAGCGGAGCGCATCAGGCa           >  1:1370537/1‑100 (MQ=255)
                                                                                         ggCCGGATAAGGCGTTTACGCCGCATCCGGCCTACAAAATCGTGCAAATTCAACATATTGCAATTCT‑‑‑CTTGTAGGCCTGATAAGCGTAGCGCATCAGGct           >  1:986839/2‑99 (MQ=18)
                                                                                         ggCCGGATAAGGCGTTCACGCCGCATCCGGCCTACAAGTTCGTGCAAATTCAATAAATTGCAATAT‑‑‑GACGTAGGCCTGATAAGCGTAGCGCATCAGGcga          >  2:643025/2‑98 (MQ=21)
                                                                                         ggCCGGATAAGGCGTTCACGCCGCATCCGGCCTACAAATCCTTGCGAATTCAATATATTGTATGGAAA‑‑‑TTC‑‑AGGCCTGATAAGCGTAGCGCATCAGGctt          >  2:1444908/2‑98 (MQ=21)
                                                                                         ggCCGGATAAGGCGTTCACGCCGCATCCGGCCTACAAAATCGCTTAAATTCAATATATTGCACCATC‑‑‑CTCGTAGGCCTGATAAGCGTAGCGCATCAGGCa           >  1:3014478/2‑100 (MQ=21)
                                                                                         ggCCGGATAAGGCGTTCACGCCGCATCCGGCCTACAAAAGCACGCAAACTCAATATATTGCAGAGATC‑‑‑ATGTAGGCCTGATAAGCGTAGCGCATCAGGCa           >  1:2143376/2‑100 (MQ=17)
                                                                                         ggCCGGATAAGGCGTTAACGCCGCATCCGGCCTACAAAATCGTGCAAATTCAATACATTGCAGAATT‑‑‑TTCGTAGGCCTGATAAGCGGAGCGCATCAGGCa           >  2:596407/2‑100 (MQ=17)
                                                                                          gCCGGATAAGGCGTTTACGCCGCATCCGGCCTACAAATTCTGCAAATTCAACAAATTGCGAATCC‑‑‑CTTGTAGGCCTGATAAGCGTAGCGCATCAGGCAtt         >  2:310594/1‑100 (MQ=18)
                                                                                          gCCGGATAAGGCGTTTACGCCGCATCCGGCCTACAAAATCGTGCTAATTCAATATATTGCAGAAAC‑‑‑CTTGTAGGCCTGATAAGCGTAGCGCATCAGGCag          >  1:2568400/1‑99 (MQ=18)
                                                                                          gCCGGATAAGGCGTTCACGCCGCATCCGGCCTACAAAAATCTTGCCAATTCAATATATTGCAGGACC‑‑‑ATGTAGGCCTGATAAGCGTAGCGCATCAGGCag          >  2:3341620/1‑99 (MQ=18)
                                                                                          gCCGGATAAGGCGTTCACGCCGCATCCGGCCTACGAAATCGTGCTAATTCAATATATTGTGGGCTC‑‑‑CACGTAGGCCTGATAAGCGTAGCGCATCAgagat          >  2:1892466/1‑96 (MQ=21)
                                                                                          gCCGGATAAGGCGTTCACGCCGCATCCGGCCTACGAAATCGTGCTAATTCAATATATTGTGGGCTC‑‑‑CACGTAGGCCTGATAAGCGTAGCGCATCAGGCag          >  1:1417800/1‑99 (MQ=18)
                                                                                          gCCGGATAAGGCGTTCACGCCGCATCCGGCCTACAAAGTCGTATAAATTCAATGTATTGCAGTTTA‑‑‑TTGGTAGGCCTGATAAGCGCAGCGCATCAGGCaa          >  1:846655/1‑99 (MQ=25)
                                                                                          gCCGGATAAGGCGTTCACGCCGCATCCGGCATACAAAATTGTGTAAATTCAATAAATTGATCGATC‑‑‑CACGTAGGCCTGATAAGCGTAGCGCATCAGGCag          >  1:1213618/1‑99 (MQ=21)
                                                                                               aTGCGGCGTTAACGCCTTATCTGGCCTACAAAATCGTGCAAATTCAATACATTGCAGAATT‑‑‑TTCGTAGGCCTGATAAGCGGAGCGCATCAGGCATTtttg     <  2:2056176/100‑1 (MQ=255)
                                                                                                 gcggcgTTAACGCCTTATCCGGCCTACAAAATCGTGCAAATTCAATACATTGCAGAATT‑‑‑TTCGTGGGCCTGATAAGCGGAGCGCATCAGGCATTtttgcc   <  2:3345524/100‑1 (MQ=255)
                                                                                                 gcggcgTTAACGCCTTATCCGGCCTACAAAATCGTGCAAATTCAATACATTGCAGAATT‑‑‑TTCGTAGGCCTGATAAGCGGAGCGCATCAGGCATTtttgcc   <  1:3067158/100‑1 (MQ=255)
                                                                                                 gcggcgTTAACGCCTTATCCGGCCTACAAAATCGTGCAAATTCAATACATTGCAGAATT‑‑‑TTCGTAGGCCTGATAAGCGGAGCGCATCAGGCATTtttgc    >  1:1156189/1‑99 (MQ=255)
                                                                                                  cggcgTTAACGCCTTATCCGGCCTACAAAATCGTGCAAATTCAATACATTGCAGAATT‑‑‑TTCGTAGGCCTGATAAGCGGAGCGCATCAGGCATTtttgcct  >  2:2550205/1‑100 (MQ=255)
                                                                                                  |                                                                                                          
GGGAATGAAGTATTGCTGCCGGATGGCGCTGCACGCAGAGGCACAGCGCCACCGGGTATACAGATTGTTGACAAAGTGCGCTTTGTTTATGCCGGATGCGGCGTTAACGCCTTATCCGGCCTACAAAATCGTGCAAATTCAATACATTGCAGAATT‑‑‑TTCGTAGGCCTGATAAGCGGAGCGCATCAGGCATTTTTGCCT  >  NC_000913/3705873‑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: