Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009133,705,9710CA28.8% 442.2 / 182.2 198intergenic (‑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 (67/74);  new base A (52/5);  total (119/79)
Fisher's exact test for biased strand distribution p-value = 2.65e-09
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 2.29e-01

Only 100 of 206 total aligned reads displayed.

AGGGAATGAAGTATTGCTGCCGGATGGCGCTGCACGCAGAGGCACAGCGCCACCGGGTATACAGATTGTTGACAAAGTGCGCTTTGTTTATGCCGGATGCGGCGTTAACGCCTTATCCGGCCTAC‑AAAATCGTGCAAATTCAATACATTGC‑‑AGA‑AT‑TTTC‑‑GTAGGCCTGATAAGCGGAGCGCATCAGGCATTTTTGC  >  NC_000913/3705872‑3706068
                                                                                                   |                                                                                                        
aGGGAATGAAGTATTGCTGCCGGATGGCGCTGCACGCAGAGGCACAGCGCCACCGGGTATACAGATTGTTGACAAAGTGCGCTTTGTTTATGCCGGATgc                                                                                                          >  1:1108931/1‑100 (MQ=255)
 gggAATGAAGTATTGCTGCCGGATGGCGCTGCACGCAGAGGCACAGCGCCACCGGGTATACAGATTGTTGACAAAGTGCGCTTTGTTTATGCCGGATgcg                                                                                                         >  2:849005/1‑100 (MQ=255)
 gggAATGAAGTATTGCTGCCGGATGGCGCTGCACGCAGAGGCACAGCGCCACCGGGTATACAGATTGTTGACAAAGTGCGCTTTGTTTATGCCGGATgcg                                                                                                         >  2:139990/1‑100 (MQ=255)
  ggAATGAAGTATTGCTGCCGGATGGCGCTGCACGCAGAGGCACAGCGCCACCGGGTATACAGATTGTTGACAAAGTGCGCTTTGTTTATGCCGGATgcgg                                                                                                        <  1:1719637/100‑1 (MQ=255)
    aaTGAAGTATTGCTGCCGGATGGCGCTGCACGCAGAGGCACAGCGCCACCGGGTATACAGATTGTTGACAAAGTGCGCTTTGTTTATGCCGGATgcggcg                                                                                                      <  2:1507491/100‑1 (MQ=255)
     aTGAAGTATTGCTGCCGGATGGCGCTGCACGCAGAGGCACAGCGCCACCGGGTATACAGATTGTTGACAAAGTGCGCTTTGTTTATGCCGGATGCGGCGt                                                                                                     >  2:591130/1‑100 (MQ=255)
     aTGAAGTATTGCTGCCGGATGGCGCTGCACGCAGAGGCACAGCGCCACCGGGTATACAGATTGTTGACAAAGTGCGCTTTGTTTATGCCGGATGCGGCGt                                                                                                     <  2:2278196/100‑1 (MQ=255)
       gAAGTATTGCTGCCGGATGGCGCTGCACGCAGAGGCACAGCGCCACCGGGTATACAGATTGTTGACAAAGTGCGCTTTGTTTATGCCGGATGCGGCGTTa                                                                                                   >  1:2670243/1‑100 (MQ=255)
       gAAGTATTGCTGCCGGATGGCGCTGCACGCAGAGGCACAGCGCCACCGGGTATACAGATTGTTGACAAAGTGCGCTTTGTTTATGCCGGATGCGGCGTTa                                                                                                   >  1:2580983/1‑100 (MQ=255)
         aGTATTGCTGCCGGATGGCGCTGCACGCAGAGGCACAGCGCCACCGGGTATACAGATTGTTGACAAAGTGCGCTTTGTTTATGCCGGATGCGGCGTTAAc                                                                                                 >  1:1294320/1‑100 (MQ=255)
          gTATTGCTGCCGGATGGCGCTGCACGCAGAGGCACAGCGCCACCGGGTATACAGATTGTTGACAAAGTGCGCTTTGTTTATGCCGGATGCGGCGTTAACg                                                                                                <  2:683830/100‑1 (MQ=255)
            aTTGCTGCCGGATGGCGCTGCACGCAGAGGCACAGCGCCACCGGGTATACAGATTGTTGACAAAGTGCGCTTTGTTTATGCCGGATGCGGCGTTAACGcc                                                                                              <  2:2214528/100‑1 (MQ=255)
             ttGCTGCCGGATGGCGCTGCACGCAGAGGCACAGCGCCACCGGGTATACAGATTGTTGACAAAGTGCGCTTTGTTTATGCCGGATGCGGCGTTAACGCCt                                                                                             <  1:2247267/100‑1 (MQ=255)
               gctgcCGGATGGCGCTGCACGCAGAGGCACAGCGCCACCGGGTATACAGATTGTTGACAAAGTGCGCTTTGTTTATGCCGGATGCGGCGTTAACGCCTTa                                                                                           <  1:1249263/100‑1 (MQ=255)
                 tgcCGGATGGCGCTGCACGCAGAGGCACAGCGCCACCGGGTATACAGATTGTTGACAAAGTGCGCTTTGTTTATGCCGGATGCGGCGTTAACGCCTTATc                                                                                         >  1:1148481/1‑100 (MQ=255)
                 tgcCGGATGGCGCTGCACGCAGAGGCACAGCGCCACCGGGTATACAGATTGTTGACAAAGTGCGCTTTGTTTATGCCGGATGCGGCGTTAACGCCTTATc                                                                                         >  2:629211/1‑100 (MQ=255)
                 tgcCGGATGGCGCTGCACGCAGAGGCACAGCGCCACCGGGTATACAGATTGTTGACAAAGTGCGCTTTGTTTATGCCGGATGCGGCGTTAACGCCTTATc                                                                                         >  1:1483994/1‑100 (MQ=255)
                    cGGATGGCGCTGCACGCAGAGGCACAGCGCCACCGGGTATACAGATTGTTGACAAAGTGCGCTTTGTTTATGCCGGATGCGGCGTTAACGCCTTATCCgg                                                                                      >  1:538761/1‑100 (MQ=255)
                    cGGATGGCGCTGCACGCAGAGGCACAGCGCCACCGGGTATACAGATTGTTGACAAAGTGCGCTTTGTTTATGCCGGATGCGGCGTTAACGCCTTATCCgg                                                                                      <  1:1648363/100‑1 (MQ=255)
                     ggATGGCGCTGCACGCAGAGGCACAGCGCCACCGGGTATACAGATTGTTGACAAAGTGTGCTTTGTTTATGCCGGATGCGGCGTTAACGCCTTATCCGGc                                                                                     <  2:2400970/100‑1 (MQ=255)
                     ggATGGCGCTGCACGCAGAGGCACAGCGCCACCGGGTATACAGATTGTTGACAAAGTGCGCTTTGTTTATGCCGGATGCGGCGTTAACGCCTTATCCGGc                                                                                     <  1:958277/100‑1 (MQ=255)
                     ggATGGCGCTGCACGCAGAGGCACAGCGCCACCGGGTATACAGATTGTTGACAAAGTGCGCTTTGTTTATGCCGGATAAGGCGTTAACGCCGCATCCGGc                                                                                     <  1:2954276/100‑1 (MQ=37)
                     ggATGGCGCTGCACGCAGAGGCACAGCGCCACCGGGTATACAGATTGTTGACAAAGTGCGCTTTGTTTATGCCGGATAAGGCGTTAACGCCGCATCCGGc                                                                                     <  2:2357976/100‑1 (MQ=37)
                       aTGGCGCTGCACGCAGAGGCACAGCGCCACCGGGTATACAGATTGTTGACAAAGTGCGCTTTGTTTATGCCGGATAAGGCGTTAACGCCGCATCCGGcat                                                                                   <  1:186130/100‑3 (MQ=37)
                         ggCGCTGCACGCAGAGGCACAGCGCCACCGGGTATACAGATTGTTGACAAAGTGCGCTTTGTTTATGCCGGATGCGGCGTTAACGCCTTATCCGGCCTAc                                                                                 >  1:2136198/1‑100 (MQ=255)
                              tGCACGCAGAGGCACAGCGCCACCGGGTATACAGATTGTTGACAAAGTGCGCTTTGTTTATGCCGGATGCGGCGTTAACGCCTTATCCGGCCTAC‑AAAAt                                                                           >  2:2731062/1‑100 (MQ=255)
                                cACGCAGAGGCACAGCGCCACCGGGTATACAGATTGTTGACAAAGTGCGCTTTGTTTATGCCGGATGCGGCGTTAACGCCTTATCCGGCCTAC‑AAAATCg                                                                         >  1:2499340/1‑100 (MQ=255)
                                 aCGCAGAGGCACAGCGCCACCGGGTATACAGATTGTTGACAAAGTGCGCTTTGTTTATGCCGGATGCGGCGTTAACGCCTTATCCGGCCTAC‑AAAATCGt                                                                        <  2:164655/100‑1 (MQ=255)
                                  cGCAGAGGCACAGCGCCACCGGGTATACAGATTGTTGACAAAGTGCGCTTTGTTTATGCCGGATGCGGCGTTAACGCCTTATCCGGCCTAC‑AAAATCGTg                                                                       <  1:1407678/100‑1 (MQ=255)
                                   gCAGAGGCACAGCGCCACCGGGTATACAGATTGTTGACAAAGTGCGCTTTGTTTATGCCGGATGCGGCGTTAACGCCTTATCCGGCCTAC‑AAAATCGTGc                                                                      <  2:1287182/100‑1 (MQ=255)
                                    cAGAGGCACAGCGCCACCGGGTATACAGATTGTTGACAAAGTGCGCTTTGTTTATGCCGGATGCGGCGTTAACGCCTTATCCGGCCTAC‑AAAATCGTGCa                                                                     >  1:2724700/1‑100 (MQ=255)
                                     agagGCACAGCGCCACCGGGTATACAGATTGTTGACAAAGTGCGCTTTGTTTATGCCGGATGCGGCGTTAACGCCTTATCCGGCCTAC‑AAAATCGTGCaa                                                                    <  1:728278/100‑1 (MQ=255)
                                      gagGCACAGCGCCACCGGGTATACAGATTGTTGACAAAGTGCGCTTTGTTTATGCCGGATGCGGCGTTAACGCCTTATCCGGCCTAC‑AAAATCGTGCaaa                                                                   >  1:1557742/1‑100 (MQ=255)
                                       agGCACAGCGCCACCGGGTATACAGATTGTTGACAAAGTGCGCTTTGTTTATGCCGGATGCGGCGTTAACGCCTTATCCGGCCTAC‑AAAATCGTGCAAAt                                                                  >  2:1060519/1‑100 (MQ=255)
                                         gCACAGCGCCACCGGGTATACAGATTGTTGACAAAGTGCGCTTTGTTTATGCCGGATGCGGCGTTAACGCCTTATCCGGCCTAC‑AAAATCGTGCAAATTc                                                                <  2:2079948/100‑1 (MQ=255)
                                          cacaGCGCCACCGGGTATACAGATTGTTGACAAAGTGCGCTTTGTTTATGCCGGATGCGGCGTTAACGCCTTCTCCGGCCTAC‑AAAATCGTGCAAATTc                                                                >  2:2469047/1‑99 (MQ=255)
                                          cacaGCGCCACCGGGTATACAGATTGTTGACAAAGTGCGCTTTGTTTATGCCGGATGCGGCGTTAACGCCTTATCCGGCCTAC‑AAAATCGTGCAAATTCa                                                               >  2:1637398/1‑100 (MQ=255)
                                            caGCGCCACCGGGTATACAGATTGTTGACAAAGTGCGCTTTGTTTATGCCGGATGCGGCGTTAACGCCTTATCCGGCCTAC‑AAAATCGTGCAAATTCAAt                                                             >  2:2659106/1‑100 (MQ=255)
                                             aGCGCCACCGGGTATACAGATTGTTGACAAAGTGCGCTTTGTTTATGCCGGATGCGGCGTTAACGCCTTATCCGGCCTAC‑AAAATCGTGCAAATTCAATa                                                            <  2:538755/100‑1 (MQ=255)
                                              gcgcCACCGGGTATACAGATTGTTGACAAAGTGCGCTTTGTTTATGCCGGATGCGGCGTTAACGCCTTATCCGGCCTAC‑AAAATCGTGCAAATTCAATAc                                                           <  2:1483440/100‑1 (MQ=255)
                                              gcgcCACCGGGTATACAGATTGTTGACAAAGTGCGCTTTGTTTATGCCGGATGCGGCGTTAACGCCTTATCCGGCCTAC‑AAAATCGTGCAAATTCAATAc                                                           >  1:772337/1‑100 (MQ=255)
                                                 ccACCGGGTATACAGATTGTTGACAAAGTGCGCTTTGTTTATGCCGGATGCGGCGTTAACGCCTTATCCGGCCTAC‑AAAATCGTGCAAATTCAATACAtt                                                        >  1:1303789/1‑100 (MQ=255)
                                                  cACCGGGTATACAGATTGTTGACAAAGTGCGCTTTGTTTATGCCGGATGCGGCGTTAACGCCTTATCCGGCCTAC‑AAAATCGTGCAAATTCAATACATTg                                                       >  2:3006490/1‑100 (MQ=255)
                                                  cACCGGGTATACAGATTGTTGACAAAGTGCGCTTTGTTTATGCCGGATGCGGCGTTAACGCCTTATCCGGCCTAC‑AAAATCGTGCAAATTCAATACATTg                                                       >  2:465919/1‑100 (MQ=255)
                                                   aCCGGGTATACAGATTGTTGACAAAGTGCGCTTTGTTTATGCCGGATGCGGCGTTAACGTCTTATCCGGCCTAC‑AAAATCGTGCTAATTCAATACAtt                                                        <  2:1871805/98‑1 (MQ=255)
                                                     cGGGTATACAGATTGTTGACAAAGTGCTCTTTGTTTATGCCGGATGCGGCGTTAACGCCTTATCCGACCTAC‑AAAATCGTGCAAATTCAATACATTGc                                                      <  2:1137404/98‑1 (MQ=255)
                                                     cGGGTATACAGATTGTTGACAAAGTGCGCTTTGTTTATGCCGGATGCGGCGTTAACGCCTTATCCGGCCTAC‑AAAATCGTGTAAATTCAATACATTGC‑‑Ag                                                  <  1:661341/100‑1 (MQ=255)
                                                      gggTATACAGATTGTTGACAAAGTGCGCTTTGTTTATGCCGGATGCGGCGTTAACGCCTTATCCGGCCTAC‑AAAATCGTGCAAATTCAATACATTGC‑‑AGa                                                 >  1:56044/1‑100 (MQ=255)
                                                      gggTATACAGATTGTTGACAAAGTGCGCTTTGTTTATGCCGGATGCGGCGTTAACGCCTTATCCGGCCTAC‑AAAATCGTGCAAATTCAATACATTGC‑‑AGa                                                 <  2:2580960/100‑1 (MQ=255)
                                                       ggTATACAGATTGTTGACAAAGTGCGCTTTGTTTATGCCGGATGCGGCGTTAACGCCTTATCCGGCCTAC‑AAAATCGTGCAAATTCAATACATTGC‑‑AGa‑a                                               >  1:792665/1‑100 (MQ=255)
                                                        gTATACAGATTGTTGACAAAGTGCGCTTTGTTTATGCCGGATGCGGCGTTAACGCCTTATCCGGCCTAC‑AAAATCGTGCAAATTCAATACATTGC‑‑AGA‑At                                              >  1:1634581/1‑100 (MQ=255)
                                                        gTATACAGATTGTTGACAAAGTGCGCTTTGTTTATGCCGGATGCGGCGTTAACGCCTTATCCGGCCTAC‑AAAATCGTGCAAATTCAATACATTGC‑‑AGA‑At                                              <  2:816131/100‑1 (MQ=255)
                                                          ataCAGATTGTTGACAAAGTGCGCTTTGTTTATGCCGGATGCGGCGTTAACGCCTTATCCGGCCTAC‑AAAATCGTGCAAATTCAATACATTGC‑‑AGA‑At‑tt                                           >  1:582416/1‑100 (MQ=255)
                                                          ataCAGATTGTTGACAAAGTGCGCTTTGTTTATGCCGGATGCGGCGTTAACGCCTTATCCGGCCTAC‑AAAATCGTGCAAATTCAATACATTGC‑‑AGA‑At‑tt                                           <  2:1330559/100‑1 (MQ=255)
                                                            aCAGATTGTTGACAAAGTGCGCTTTGTTTATGCCGGATGCGGCGTTAACGCCTTATCCGGCCTAC‑AAAATCGTGCAAATTCAATACATTGC‑‑AGA‑AT‑TTTc                                         <  2:476873/100‑1 (MQ=255)
                                                             cAGATTGTTGACAAAGTGCGCTTTGTTTATGCCGGATGCGGCGTTAACGCCTTATCCGGCCTAC‑AAAATCGTGCAAATTCAATACATTGC‑‑AGA‑AT‑TTTC‑‑g                                      <  1:139991/100‑1 (MQ=255)
                                                             cAGATTGTTGACAAAGTGCGCTTTGTTTATGCCGGATGCGGCGTTAACGCCTTATCCGGCCTAC‑AAAATCGTGCAAATTCAATACATTGC‑‑AGA‑AT‑TTTC‑‑g                                      <  2:772330/100‑1 (MQ=255)
                                                               gATTGTTGACAAAGTGCGCTTTGTTTATGCCGGATGCGGCGTTAACGCCTTATCCGGCCTAC‑AAAATCGTGCAAATTCAATACATTGC‑‑AGA‑AT‑TTTC‑‑GTa                                    >  2:2918602/1‑100 (MQ=255)
                                                                  tgttgACAAAGTGCGCTTTGTTTATGCCGGATGCGGCGTTAACGCCTTATCCGGCCTAC‑AAAATCGTGCAAATTCAATACATTGC‑‑AGA‑AT‑TTTC‑‑GTAGGc                                 >  2:1267205/1‑100 (MQ=255)
                                                                      gACAAAGTGCGCTTTGTTTATGCCGGATGCGGCGTTAACGCCTTATCCGGCCTAC‑AAAATCGTGCAAATTCAATACATTGC‑‑AGA‑AT‑TTTC‑‑GTAGGCCTGa                             >  1:1542260/1‑100 (MQ=255)
                                                                      gACAAAGTGCGCTTTGTTTATGCCGGATGCGGCGTTAACGCCTTATCCGGCCTAC‑AAAATCGTGCAAATTCAATACATTGC‑‑AGA‑AT‑TTTC‑‑GTAGGCCTGa                             >  2:1433147/1‑100 (MQ=255)
                                                                      gACAAAGTGCGCTTTGTTTATGCCGGATGCGGCGTTAACGCCTTATCCGGCCTAC‑AAAATCGTGCAAATTCAATACATTGC‑‑AGA‑AT‑TTTC‑‑GTAGGCCTGa                             >  1:23376/1‑100 (MQ=255)
                                                                       aCAAAGTGCGCTTTGTTTATGCCGGATGCGGCGTTAACGCCTTATCCGGCCTAC‑AAAATCGTGCAAATTCAATACATTGC‑‑AGA‑AT‑TTTC‑‑GTAGGCCTGAt                            >  2:202671/1‑100 (MQ=255)
                                                                       aCAAAGTGCGCTTTGTTTATGCCGGATGCGGCGTTAACGCCTTATCCGGCCTAC‑AAAATCGTGCAAATTCAATACATTGC‑‑AGA‑AT‑TTTC‑‑GTAGGCCTGAt                            >  1:740102/1‑100 (MQ=255)
                                                                        cAAAGTGCGCTTTGTTTATGCCGGATGCGGCGTTAACGCCTTATCCGGCCTAC‑AAAATCGTGCAAATTCAATACATTGC‑‑AGA‑AT‑TTTC‑‑GTAGGCCTGATa                           <  2:1218333/100‑1 (MQ=255)
                                                                        cAAAGTGCGCTTTGTTTATGCCGGATGCGGCGTTAACGCCTTATCCGGCCTAC‑AAAATCGTGCAAATTCAATACATTGC‑‑AGA‑AT‑TTTC‑‑GTAGGCCTGATa                           <  2:2675614/100‑1 (MQ=255)
                                                                         aaaGTGCGCTTTGTTTATGCCGGATGCGGCGTTAACGCCTTATCCGGCCTAC‑AAAATCGTGCAAATTCAATACATTGC‑‑AGA‑AT‑TTTC‑‑GTAGGCCTGATaa                          >  1:919992/1‑100 (MQ=255)
                                                                           aGTGCGCTTTGTTTATGCCGGATGCGGCGTTAACGCCTTATCCGGCCTAC‑AAAATCGTGCAAATTCAATACATTGC‑‑AGA‑AT‑TTTC‑‑GTAGGCCTGATAAGc                        >  2:1145561/1‑100 (MQ=255)
                                                                              gcgcTTTGTTTATGCCGGATGCGGCGTTAACGCCTTATCCGGCCTAC‑AAAATCGTGCAAATTCAATACATTGC‑‑AGA‑AT‑TTTC‑‑GTAGGCCTGAtgcgctcc                     >  2:2959428/1‑93 (MQ=255)
                                                                              gcgcTTTGTTTATGCCGGATGCGGCGTTAACGCCTTATCCGGCCTAC‑AAAATCGTGCAAATTCAATACATTGC‑‑AGA‑AT‑TTTC‑‑GTAGGCCTGATAAGCGGa                     >  2:1910813/1‑100 (MQ=255)
                                                                                gcTTTGTTTATGCCGGATGCGGCGTTAACGCCTTATCCGGCCTAC‑AAAATCGTGCAAATTCAATACATTGC‑‑AGA‑AT‑TTTC‑‑GTAGGCCTGATAAGCGGAgc                   >  1:804200/1‑100 (MQ=255)
                                                                                 cTTTGTTTATGCCGGATGCGGCGTTAACGCCTTATCCGGCCTAC‑AAAATCGTGCAAATTCAATACATTGC‑‑AGA‑AT‑TTTC‑‑GTAGGCCTGATAAGCGGAgcg                  <  2:578073/100‑1 (MQ=255)
                                                                                 cTTTGTTTATGCCGGATGCGGCGTTAACGCCTTATCCGGCCTAC‑AAAATCGTGCAAATTCAATACATTGC‑‑AGA‑AT‑TTTC‑‑GTAGGCCTGATAAGCGGAgcg                  >  1:1461731/1‑100 (MQ=255)
                                                                                 cTTTGTTTATGCCGGATGCGGCGTTAACGCCTTATCCGGCCTAC‑AAAATCGTGCAAATTCAATACATTGC‑‑AGA‑AT‑TTTC‑‑GTAGGCCTGATAAGCGGAgcg                  <  1:202674/100‑1 (MQ=255)
                                                                                 cTTTGTTTATGCCGGATGCGGCGTTAACGCCTTATCCGGCCTAC‑AAAATCGTGCAAATTCAATACATTGC‑‑AGA‑AT‑TTTC‑‑GTAGGCCTGATAAGCGGAgcg                  >  1:2441741/1‑100 (MQ=255)
                                                                                  tttttgtaggCCGGATAAGGCGTTCACGCCGCATCCGGCCTACAAAAATCTTGCCAATTCAATATATTGC‑‑AGG‑AC‑CAT‑‑‑GTAGGCCTGATAAGCGTAgcgc                 >  1:1794836/10‑100 (MQ=18)
                                                                                   ttttgtaggCCGGATAAGGCGTTCACGCCGCATCCGGCCTAC‑AAAATCGCTTAAATTCAATATATTGC‑‑ACC‑AT‑CCTC‑‑GTAGGCCTGATAAGCGTAGCGCa                >  1:602934/9‑100 (MQ=18)
                                                                                    tttgtaggCCGGATAAGGCGTTCACGCCGCATCCGGCCTAC‑AAAAGCACGCAAACTCAATATATTGC‑‑AGAGAT‑CAT‑‑‑GTAGGCCTGATAAGCGTAGCGCAt               >  1:2353187/8‑100 (MQ=17)
                                                                                     gTTTATGCCGGATGCGGCGTTAACGCCTTATCCGGCCTAC‑AAAATCGTGCAAATTCAATACATTGC‑‑AGA‑AT‑TTTC‑‑GTAGGCCTGATAAGCGGAGCGCATc              >  1:446705/1‑100 (MQ=255)
                                                                                      tgtaggCCGGATAAGGCGTTTACGCCGCATCCGGCCTAC‑AAATTC‑TGCAAATTCAACAAATTGC‑‑‑GA‑ATCCCTT‑‑GTAGGCCTGATAAGCGTAGCGCATCAg            >  1:2423937/6‑100 (MQ=18)
                                                                                      tgtaggCCGGATAAGGCGTTTACGCCGCATCCGGCCTAC‑AAATTC‑TGCAAATTCAACAAATTGC‑‑‑GA‑ATCCCTT‑‑GTAGGCCTGATAAGCGTAGCGCATCAg            >  2:34661/6‑100 (MQ=18)
                                                                                      tgtaggCCGGATAAGGCGTTTACGCCGCATCCGGCCTAC‑AAAATCGTGCTAATTCAATATATTGC‑‑AGA‑AA‑CCTT‑‑GTAGGCCTGATAAGCGTAGCGCATCa             >  2:1101117/6‑100 (MQ=18)
                                                                                      tgtaggCCGGATAAGGCGTTCACGCCGCATCCGGCCTAC‑AAAAGCATGAAAATTCAATATATTGC‑‑AGG‑AG‑CTGC‑‑GTAGGCCTGATAAGCGTAGCGCATCa             >  1:682263/6‑100 (MQ=17)
                                                                                        taggCCGGATAAGGCGTTTACGCCGCATCCGGCCTAC‑AAAACAGTGCAAATTCAATATATTGC‑‑AAT‑CT‑TCTG‑‑GTAGGCCTGATAAGCATCGCGCATCAgg           >  1:14439/4‑100 (MQ=18)
                                                                                        taggCCGGATAAGGCGTTCACGCCGCATCCGGCCTAC‑GAAATCGTGCTAATTCAATATATTGT‑‑GGG‑CT‑CCAC‑‑GTAGGCCTGATAAGCGTAGCGCATCAgg           >  2:2779604/4‑100 (MQ=21)
                                                                                        taggCCGGATAAGGCGTTCACGCCGCATCCGGCCTAC‑AAAACTTTGCAAATTCAATATATTGC‑‑‑‑A‑AT‑CTCCGTGTAGGCCTGATAAGCGTAGCGCATCAgg           >  2:1707380/4‑100 (MQ=18)
                                                                                        taggCCGGATAAGGCGTTCACGCCGCATCCGGCCTAC‑AAAACTTTGCAAATTCAATATATTGC‑‑‑‑A‑AT‑CTCCGTGTAGGCCTGATAAGCGTAGCGCATCAgg           >  2:409985/4‑100 (MQ=18)
                                                                                        taggCCGGATAAGGCGTTCACGCCGCATCCGGCCTAC‑‑AAATCATGCAAATCCAAAATATTGC‑‑AGG‑AG‑TTAT‑‑GTAGGCCTGATAAGCGTAGCGCATCAGGc          >  2:1565302/4‑100 (MQ=21)
                                                                                        taggCCGGATAAGGCGTTCACGCCGCATCCGGCCTAC‑‑AAACCGCGCAAATTCAATATATTGCGGAGA‑AA‑AT‑‑‑‑GTAGGCCTGATAAGCGTAGCGCATCAGGc          >  1:281462/4‑100 (MQ=18)
                                                                                         aggCCGGATAAGGCGTTCACGCCGCATCCGGCCTACAAAAATCTTGCCAATTCAATATATTGC‑‑AGG‑AC‑CAT‑‑‑GTAGGCCTGATAAGCGTAGCGCATCAGGc          >  1:856168/3‑100 (MQ=18)
                                                                                         aggCCGGATAAGGCGTTCACGCCGCATCCGGCCTAC‑AAGTTCGTGCAAATTCAATAAATTGC‑‑A‑A‑TA‑TGAC‑‑GTAGGCCTGATAAGCGTAGCGCATCAGGcg         >  2:2927278/3‑99 (MQ=21)
                                                                                         aggCCGGATAAGGCGTTCACGCCGCATCCGGCCTAC‑AAAATTGTGTAAATTCAATAAATTGA‑‑TCG‑AT‑CCAC‑‑GTAGGCCTGATAAGCGTAGCGCATCAGGc          >  1:37439/3‑100 (MQ=21)
                                                                                         aggCCGGATAAGGCGTTCACGCCGCATCCGGCCTAC‑‑AAACCGCGCAAATTCAATATATTGC‑‑GGA‑GA‑AAAT‑‑GTAGGCCTGATAAGCGTAGCGCATCAGGct         >  1:914955/3‑99 (MQ=18)
                                                                                         aggCCGGATAAGGCGTTAACGCCGCATCCGGCCTAC‑AAAATCGTGCAAATTCAATACATTGC‑‑AGA‑AT‑TTTC‑‑GTAGGCCTGATAAGCGGAGCGCATCCGGc          >  1:1359966/3‑100 (MQ=17)
                                                                                          ggCCGGATAAGGCGTTTACGCCGCATCCGGCCTAC‑AAATTC‑TGCAAATTCAACAAATTGC‑‑‑GA‑ATCCCTT‑‑GTAGGCCTGATAAGCGTAGCGCATCAGGCAt        >  1:2566432/2‑100 (MQ=18)
                                                                                          ggCCGGATAAGGCGTTTACGCCGCATCCGGCCTAC‑AAAATCGTGCTAATTCAATATATTGC‑‑AGA‑AA‑CCTT‑‑GTAGGCCTGATAAGCGTAGCGCATCAGGCa         >  1:2008127/2‑100 (MQ=18)
                                                                                          ggCCGGATAAGGCGTTCACGCCGCATCCGGCCTAC‑AAAACTGTGCAAATTCAATATATTGC‑‑AGG‑AA‑ACAC‑‑GTAGGCCTGATAAGCGAAGCGCATCAGGCa         >  1:493374/2‑100 (MQ=18)
                                                                                          ggCCGGATAAGGCGTTCACGCCGCATCCGGCCTAC‑‑AAATCATGCAAATCCAAAATATTGC‑‑AGG‑AG‑TTAT‑‑GTAGGCCTGATAAGCGTAGCGCATCAGGCAt        >  1:2828617/2‑100 (MQ=21)
                                                                                           gCCGGATAAGGCGTTCACGCCGCATCCGGCCTAC‑AAAGTCGTATAAATTCAATGTATTGC‑‑AGT‑TT‑ATTG‑‑GTAGGCCTGATAAGCGCAGCGCATCAGGCaa        >  1:3034754/1‑99 (MQ=21)
                                                                                                aTGCGGCGTTAACGCCTTATCCGGCCTAC‑AAAATCGTGCAAATTCAATACATTGC‑‑AGA‑AT‑TTTC‑‑GTAGGCCTGATAAGCGGAGCGCATCAGGCATTtttg   <  1:2539529/100‑1 (MQ=255)
                                                                                                   |                                                                                                        
AGGGAATGAAGTATTGCTGCCGGATGGCGCTGCACGCAGAGGCACAGCGCCACCGGGTATACAGATTGTTGACAAAGTGCGCTTTGTTTATGCCGGATGCGGCGTTAACGCCTTATCCGGCCTAC‑AAAATCGTGCAAATTCAATACATTGC‑‑AGA‑AT‑TTTC‑‑GTAGGCCTGATAAGCGGAGCGCATCAGGCATTTTTGC  >  NC_000913/3705872‑3706068

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