Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009133,705,9850TC25.0% 289.6 / 98.3 140intergenic (‑195/+113)dppB/dppAdipeptide/heme ABC transporter permease/dipeptide/heme ABC transporter periplasmic binding protein; dipeptide chemotaxis receptor
Reads supporting (aligned to +/- strand):  ref base T (50/55);  new base C (29/6);  total (79/61)
Fisher's exact test for biased strand distribution p-value = 3.15e-04
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 3.43e-01

Only 100 of 144 total aligned reads displayed.

TGCTGCCGGATGGCGCTGCACGCAGAGGCACAGCGCCACCGGGTATACAGATTGTTGACAAAGTGCGCTTTGTTTATGCCGGATGCGGCGTTAACGCCTTATCCGGCCTACAAAATCGTGCAAATTCAATACATTGCAGAATT‑‑‑TTC‑‑GTAGGCCTGATAAGCGGAGCGCATCAGGCATTTTTGCCTTTGCCATCAGTCTT  >  NC_000913/3705886‑3706084
                                                                                                   |                                                                                                           
tgctgcCGGATGGCGCTGCACGCAGAGGCACAGCGCCACCGGGTATACAGATTGTTGACAAAGTGCGCTTTGTTTATGCCGGATGCGGCGTTAACGCCtt                                                                                                             <  1:233429/100‑1 (MQ=255)
 gctgcCGGATGGCGCTGCACGCAGAGGCACAGCGCCACCGGGTATACAGATTGTTGACAAAGTGCGCTTTGTTTATGCCGGATGCGGCGTTAACGCCTTa                                                                                                            >  1:102245/1‑100 (MQ=255)
 gctgcCGGATGGCGCTGCACGCAGAGGCACAGCGCCACCGGGTATACAGATTGTTGACAAAGTGCGCTTTGTTTATGCCGGATGCGGCGTTAACGCCTTa                                                                                                            >  2:2099330/1‑100 (MQ=255)
 gctgcCGGATGGCGCTGCACGCAGAGGCACAGCGCCACCGGGTATACAGATTGTTGACAAAGTGCGCTTTGTTTATGCCGGATGCGGCATTAACGCCTTa                                                                                                            <  2:1850961/100‑1 (MQ=255)
   tgcCGGATGGCGCTGCACGCAGAGGCACAGCGCCACCGGGTATACAGATTGTTGACAAAGTGCGCTTTGTTTATGCCGGATGCGGCGTTAACGCCTTATc                                                                                                          <  2:1011647/100‑1 (MQ=255)
    gcCGGATGGCGCTGCACGCAGAGGCACAGCGCCACCGGGTATACAGATTGTTGACAAAGTGCGCTTTGTTTATGCCGGATGCGGCGTTAACGCCTTATcc                                                                                                         >  2:1242004/1‑100 (MQ=255)
     ccGGATGGCGCTGCACGCAGAGGCACAGCGCCACCGGGTATACAGATTGTTGACAAAGTGCGCTTTGTTTATGCCGGATGCGGCGTTAACGCCTTATCCg                                                                                                        >  1:567692/1‑100 (MQ=255)
     ccGGATGGCGCTGCACGCAGAGGCACAGCGCCACCGGGTATACAGATTGTTGACAAAGTGCGCTTTGTTTATGCCGGATGCGGCGTTAACGCCTTATCCg                                                                                                        >  1:2295884/1‑100 (MQ=255)
      cGGATGGCGCTGCACGCAGAGGCACAGCGCCACCGGGTATACAGATTGTTGACAAAGTGCGCTTTGTTTATGCCGGATGCGGCGTTAACGCCTTATCCgg                                                                                                       <  1:2266537/100‑1 (MQ=255)
      cGGATGGCGCTGCACGCAGAGGCACAGCGCCACCGGGTATACAGATTGTTGACAAAGTGCGCTTTGTTTATGCCGGATGCGGCGTTAACGCCTTATACgg                                                                                                       >  2:187018/1‑100 (MQ=255)
      cGGATGGCGCTGCACGCAGAGGCACAGCGCCACCGGGTATACAGATTGTTGACAAAGTGCGCTTTGTTTATGCCGGATGCGGCGTTAACGCCGCATCCgg                                                                                                       <  1:362523/100‑1 (MQ=255)
       ggATGGCGCTGCACGCAGAGGCACAGCGCCACCGGGTATACAGATTGTTGACAAAGTGCGCTTTGTTTATGCCGGATGCGGCGTTAACGCCTTATCCGGc                                                                                                      >  2:910105/1‑100 (MQ=255)
       ggATGGCGCTGCACGCAGAGGCACAGCGCCACCGGGTATACAGATTGTTGACAAAGTGCGCTTTGTTTATGCCGGATGCGGCGTTAACGCCTTATCAGGc                                                                                                      >  1:1763497/1‑100 (MQ=255)
        gATGGCGCTGCACGCAGAGGCACAGCGCCACCGGGTATACAGATTGTTGACAAAGTGCGCTTTGTTTATGCCGGATGCGGCGTTAACGCCTTATCCGGcc                                                                                                     <  1:3898/100‑1 (MQ=255)
        gATGGCGCTGCACGCAGAGGCACAGCGCCACCGGGTATACAGATTGTTGACAAAGTGCGCTTTGTTTATGCCGGATGCGGCGTTAACGCCTTATCCGGcc                                                                                                     <  2:1773164/100‑1 (MQ=255)
        gATGGCGCTGCACGCAGAGGCACAGCGCCACCGGGTATACAGATTGTTGACAAAGTGCGCTTTGTTTATGCCGGATAAGGCGTTAACGCCGCATCCGGca                                                                                                     <  1:1106629/100‑2 (MQ=37)
         aTGGCGCTGCACGCAGAGGCACAGCGCCACCGGGTATACAGATTGTTGACAAAGTGCGCTTTGTTTATGCCGGATGCGGCGTTAACGCCTTATCCGGCCt                                                                                                    <  1:1322538/100‑1 (MQ=255)
         aTGGCGCTGCACGCAGAGGCACAGCGCCACCGGGTATACAGATTGTTGACAAAGTGCGCTTTGTTTATGCCGGATAAGGCGTTAACGCCGCATCCGGcat                                                                                                    <  1:1868756/100‑3 (MQ=37)
           ggCGCTGCACGCAGAGGCACAGCGCCACCGGGTATACAGATTGTTGACAAAGTGCGCTTTGTTTATGCCGGATGCGGCGTTAACGCCTTATCCGGCCTAc                                                                                                  >  2:269211/1‑100 (MQ=255)
                tGCACGCAGAGGCACAGCGCCACCGGGTATACAGATTGTTGACAAAGTGCGCTTTGTTTATGCCGGATGCGGCGTTAACGCCTTATCCGGCCTACAAAAt                                                                                            <  1:805134/100‑1 (MQ=255)
                  cACGCAGAGGCACAGCGCCACCGGGTATACAGATTGTTGACAAAGTGCGCTTTGTTTATGCCGGATGCGGCGTTAACGCCTTATCCGGCCTACAAAATCg                                                                                          <  2:60776/100‑1 (MQ=255)
                     gCAGAGGCACAGCGCCACCGGGTATACAGATTGTTGACAAAGTGCGCTTTGTTTATGCCGGATGCGGCGTTAACGCCTTATCCGGCCTACAAAATCGTGc                                                                                       >  1:1661403/1‑100 (MQ=255)
                     gCAGAGGCACAGCGCCACCGGGTATACAGATTGTTGACAAAGTGCGCTTTGTTTATGCCGGATGCGGCGTTAACGCCTTATCCGGCCTACAAAATCGTGc                                                                                       <  1:1802865/100‑1 (MQ=255)
                     gCAGAGGCACAGCGCCACCGGGTATACAGATTGTTGACAAAGTGCGCTTTGTTTATGCCGGATGCGGCGTTAACGCCTTATCCGGCCTACAAAATCGTGc                                                                                       <  1:215192/100‑1 (MQ=255)
                      cAGAGGCACAGCGCCACCGGGTATACAGATTGTTGACAAAGTGCGCTTTGTTTATGCCGGATGCGGCGTTAACGCCTTATCCGGCCTACAAAATCGTGCa                                                                                      <  2:294214/100‑1 (MQ=255)
                       agagGCACAGCGCCACCGGGTATACAGATTGTTGACAAAGTGCGCTTTGTTTATGCCGGATGCGGCGTTAACGCCTTATCCGGCCTACAAAATCGTGCaa                                                                                     <  2:750395/100‑1 (MQ=255)
                           gCACAGCGCCACCGGGTATACAGATTGTTGACAAAGTGCGCTTTGTTTATGCCGGATGCGGCGTTAACGCCTTATCCGGCCTACAAAATCGTGCAAATTc                                                                                 <  2:173243/100‑1 (MQ=255)
                            nacaGCGCCACCGGGTATACAGATTGTTGACAAAGTGCGCTTTGTTTATGCCGGATGCGGCGTTAACGCCTTATCCGGCCTACAAAATCGTGCAAATTCa                                                                                >  2:131054/2‑100 (MQ=255)
                            cacaGCGCCACCGGGTATACAGATTGTTGACAAAGTGCGCTTTGTTTATGCCGGATGCGGCGTTAACGCCTTATCCGGCCTACAAAATCGTGCAAATTCa                                                                                >  2:769925/1‑100 (MQ=255)
                              caGCGCCACCGGGTATACAGATTGTTGACAAAGTGCGCTTTGTTTATGCCGGATGCGGCGTTAACGCCTTATCCGGCCTACAAAATCGTGCAAATTCAAt                                                                              <  1:1241736/100‑1 (MQ=255)
                              caGCGCCACCGGGTATACAGATTGTTGACAAAGTGCGCTTTGTTTATGCCGGATGCGGCGTTAACGCCTTATCCGGCCTACAAAATCGTGCAAATTCAAt                                                                              >  2:1841885/1‑100 (MQ=255)
                               aGCGCCACCGGGTATACAGATTGTTGACAAAGTGCGCTTTGTTTATGCCGGATGCGGCGTTAACGCCTTATCCGGCCTACAAAATCGTGCAAATTCAATa                                                                             >  2:82293/1‑100 (MQ=255)
                               aGCGCCACCGGGTATACAGATTGTTGACAAAGTGCGCTTTGTTTATGCCGGATGCGGCGTTAACGCCTTATCCGGCCTACAAAATCGTGCAAATTCAATa                                                                             <  1:867530/100‑1 (MQ=255)
                                gcgcCACCGGGTATACAGATTGTTGACAAAGTGCGCTTTGTTTATGCCGGATGCGGCGTTAACGCCTTATCCGGCCTACAAAATCGTGCAAATTCAATAc                                                                            <  1:939965/100‑1 (MQ=255)
                                gcgcCACCGGGTATACAGATTGTTGACAAAGTGCGCTTTGTTTATGCCGGATGCGGCGTTAACGCCTTATCCGGCCTACAAAATCGTGCAAATTCAATAc                                                                            >  1:176020/1‑100 (MQ=255)
                                gcgcCACCGGGTATACAGATTGTTGACAAAGTGCGCTTTGTTTATGCCGGATGCGGCGTTAACGCCTTATCCGGCCTACAAAATCGTGCAAATTCAATAc                                                                            <  1:1595710/100‑1 (MQ=255)
                                  gcCACCGGGTATACAGATTGTTGACAAAGTGCGCTTTGTTTATGCCGGATGCGGCGTTAACGCCTTATCCGGCCTACAAAATCGTGCAAATTCAATACAt                                                                          >  1:301014/1‑100 (MQ=255)
                                  gcCACCGGGTATACAGATTGTTGACAAAGTGCGCTTTGTTTATGCCGGATGCGGCGTTAACGCCTTATCCGGCCTACAAAATCGTGCAAATTCAATACAt                                                                          <  2:1933879/100‑1 (MQ=255)
                                       cGGGTATACAGATTGTTGACAAAGTGCGCTTTGTTTATGCCGGATGCGGCGTTAACGCCTTATCCGGCCTACAAAATCGTGCAAATTCAATACATTGCAg                                                                     <  2:1648356/100‑1 (MQ=255)
                                        gggTATACAGATTGTTGACAAAGTGCGCTTTGTTTATGCCGGATGCGGCGTTAACGCCTTATCCGGCCTACAAAATCGTGCAAATTCAATACATTGCAGa                                                                    <  2:1112582/100‑1 (MQ=255)
                                        gggTATACAGATTGTTGACAAAGTGCGCTTTGTTTATGCCGGATGCGGCGTTAACGCCTTATCCGGCCTACAAAATCGTGCAAATTCAATACATTGCAGa                                                                    <  1:20033/100‑1 (MQ=255)
                                           tataCAGATTGTTGACAAAGTGCGCTTTGTTTATGCCGGATGCGGCGTTAACGCCTTATCCGGCCTACAAAATCGTGCAAATTCAATACATTGCAGAAtt                                                               >  1:890705/1‑100 (MQ=255)
                                             taCAGATTGTTGACAAAGTGCGCTTTGTTTATGCCGGATGCGGCGTTAACGCCTTATCCGGCCTACAAAATCGTGCAAATTCAATACATTGCAGAAtt‑‑‑tt                                                          >  2:1494455/1‑100 (MQ=255)
                                               cAGATTGTTGACAAAGTGCGCTTTGTTTATGCCGGATGCGGCGTTAACGCCTTATCCGGCCTACAAAATCGTGCAAATTCAATACATTGCAGAATT‑‑‑TTC‑‑g                                                      <  1:1507375/100‑1 (MQ=255)
                                                aGATTGTTGACAAAGTGCGCTTTGTTTATGCCGGATGCGGCGTTAACGCCTTATCCGGCCTACAAAATCGTGCAAATTCAATACATTGCAGAATT‑‑‑TTC‑‑Gt                                                     <  2:503316/100‑1 (MQ=255)
                                                 gATTGTTGACAAAGTGCGCTTTGTTTATGCCGGATGCGGCGTTAACGCCTTATCCGGCCTACAAAATCGTGCAAATTCAATACATTGCAGAATT‑‑‑TTC‑‑GTa                                                    <  1:485966/100‑1 (MQ=255)
                                                 gATTGTTGACAAAGTGCGCTTTGTTTATGCCGGATGCGGCGTTAACGCCTTATCCGGCCTACAAAATCGTGCAAATTCAATACATTGCAGAATT‑‑‑TTC‑‑GTa                                                    <  2:301014/100‑1 (MQ=255)
                                                     gttgACAAAGTGCGCTTTGTTTATGCCGGATGCGGCGTTAACGCCTTATCCGGCCTACAAAATCGTGCAAATTCAATACATTGCAGAATT‑‑‑TTC‑‑GTAGGcc                                                >  2:775015/1‑100 (MQ=255)
                                                      ttgACAAAGTGCGCTTTGTTTATGCCGGATGCGGCGTTAACGCCTTATCCGGCCTACAAAATCGTGCAAATTCAATACATTGCAGAATT‑‑‑TTC‑‑GTAGGCCt                                               <  1:1662706/100‑1 (MQ=255)
                                                       tgACAAAGTGCGCTTTGTTTATGCCGGATGCGGCGTTAACGCCTTATCCGGCCTACAAAATCGTGCAAATTCAATACATTGCAGAATT‑‑‑TTC‑‑GTAGGCCTg                                              >  1:1525405/1‑100 (MQ=255)
                                                        gACAAAGTGCGCTTTGTTTATGCCGGATGCGGCGTTAACGCCTTATCCGGCCTACAAAATCGTGCAAATTCAATACATTGCAGAATT‑‑‑TTC‑‑GTAGGCCTGa                                             >  2:354547/1‑100 (MQ=255)
                                                         aCAAAGTGCGCTTTGTTTATGCCGGATGCGGCGTTAACGCCTTATCCGGCCTACAAAATCGTGCAAATTCAATACATTGCAGAATT‑‑‑TTC‑‑GTAGGCCTGAt                                            >  1:71838/1‑100 (MQ=255)
                                                          cAAAGTGCGCTTTGTTTATGCCGGATGCGGCGTTAACGCCTTATCCGGCCTACAAAATCGTGCAAATTCAATACATTGCAGAATT‑‑‑TTC‑‑GTAGGCCTGATa                                           <  1:890441/100‑1 (MQ=255)
                                                          cAAAGTGCGCTTTGTTTATGCCGGATGCGGCGTTAACGCCTTATCCGGCCTACAAAATCGTGCAAATTCAATACATTGCAGAATT‑‑‑TTC‑‑GTAGGCCTGATa                                           <  1:641446/100‑1 (MQ=255)
                                                          cAAAGTGCGCTTTGTTTATGCCGGATGCGGCGTTAACGCCTTATCCGGCCTACAAAATCGTGCAAATTCAATACATTGCAGAATT‑‑‑TTC‑‑GTAGGCCTGATa                                           <  2:1054875/100‑1 (MQ=255)
                                                              gTGCGCTTTGTTTATGCCGGATGCGGCGTTAACGCCTTATCCGGCCTACAAAATCGTGCAAATTCAATACATTGCAGAATT‑‑‑TTC‑‑GTAGGCCTGATAAGCg                                       <  2:2156714/100‑1 (MQ=255)
                                                               tGCGCTTTGTTTATGCCGGATGCGGCGTTAACGCCTTATCCGGCCTACAAAATCGTGCAAATTCAATACATTGTAGAATT‑‑‑TTC‑‑GTAGGCCTGATAAGCgg                                      <  1:807176/100‑1 (MQ=255)
                                                                gcgcTTTGTTTATGCCGGATGCGGCGTTAACGCCTTATCCGGCCTACAAAATCGTGCAAATTCAATACATTGCAGAATT‑‑‑TTC‑‑GTAGGCCTGATAAGCGGa                                     >  1:2379716/1‑100 (MQ=255)
                                                                 cgcTTTGTTTATGCCGGATGCGGCGTTAACGCCTTATCCGGCCTACAAAATCGTGCAAATTCAATACATTGCAGAATT‑‑‑TTC‑‑GTAGGCCTGATAAGCGGAg                                    >  1:1326380/1‑100 (MQ=255)
                                                                 cgcTTTGTTTATGCCGGATGCGGCGTTAACGCCTTATCCGGCCTACAAAATCGTGCAAATTCAATACATTGCAGAATT‑‑‑TTC‑‑GTAGGCCTGATAAGCGGAg                                    <  2:816791/100‑1 (MQ=255)
                                                                 cgcTTTGGTTATGCCGGATGCGGCGTTAACGCTTTATCCGGCCTACAAAATCGTGCACATTCAATACATTGCAGAATT‑‑‑TTC‑‑GTAGGCCTGATAAGCggn                                     <  2:900631/99‑2 (MQ=255)
                                                                  gcTTTGTTTATGCCGGATGCGGCGTTAACGCCTTATCCGGCCTACAAAATCGTGCAAATTCAATACATTGCAGAATT‑‑‑TTC‑‑GTAGGCCTGATAAGCGGAgc                                   >  2:1868082/1‑100 (MQ=255)
                                                                   cTTTGTTTATGCCGGATGCGGCGTTAACGCCTTATCCGGCCTACAAAATCGTGCAAATTCAATACATTGCAGAATT‑‑‑TTC‑‑GTAGGCCTGATAAGCGGAgcg                                  >  1:862490/1‑100 (MQ=255)
                                                                   cTTTGTTTATGCCGGATGCGGCGTTAACGCCTTATCCGGCCTACAAAATCGTGCAAATTCAATACATTGCAGAATT‑‑‑TTC‑‑GTAGGCCTGATAAGCGGAgcg                                  <  1:1215820/100‑1 (MQ=255)
                                                                   cTTTGTTTATGCCGGATGCGGCGTTAACGCCTTATCCGGCCTACAAAATCGTGCAAATTCAATACATTGCAGAATT‑‑‑TTC‑‑GTAGGCCTGATAAGCGGAgcg                                  <  1:468090/100‑1 (MQ=255)
                                                                    tttGTTTATGCCGGATGCGGCGTTAACGCCTTATCCGGCCTACAAAATCGTGCAAATTCAATACATTGCAGAATT‑‑‑TTC‑‑GTAGGCCTGATAAGCGGAgcgc                                 >  1:497555/1‑100 (MQ=255)
                                                                       gTTTATGCCGGATGCGGCGTTAACGCCTTATCCGGCCTACAAAATCGTGCAAATTCAATACATTGCAGAATT‑‑‑TTC‑‑GTAGGCCTGATAAGCGGAGCGCATc                              >  2:1935290/1‑100 (MQ=255)
                                                                           aggCCGGATAAGGCGTTTACGCCGCATCCGGCCTACAAAAATGTGCAAATTCAATAAATTGC‑‑AATTCAACTT‑‑GTAGGCCTGATAAGCGCAGCGCATCAgg                           >  1:1415787/3‑100 (MQ=21)
                                                                           aggCCGGATAAGGCGTTCACGCCGCATCCGGCCTACAAAATCGCTTAAATTCAATATATTGCACCATC‑‑‑CTC‑‑GTAGGCCTGATAAGCGTAGCGCATCAGGc                          >  2:2207506/3‑100 (MQ=18)
                                                                           aggCCGGATAAGGCGTTCACGCCGCATCCGGCCTACAAAAACCTGCAAATTCAATATATTGCAGGATA‑‑‑CGT‑‑GTAGGCCTGATAAGCGTAGCGCATCAGGc                          >  2:1141847/3‑100 (MQ=17)
                                                                            ggCCGGATAAGGCGTTTACGCCGCATCCGGCCTACAAAATCGTGCTAATTCAATATATTGCAGAAAC‑‑‑CTT‑‑GTAGGCCTGATAAGCGTAGCGCATCAGGCa                         >  2:1626817/2‑100 (MQ=18)
                                                                            ggCCGGATAAGGCGTTCACGCCGCATCCGGCCTACAAGTTCGTGCAAATTCAATAAATTGCAATAT‑‑‑GAC‑‑GTAGGCCTGATAAGCGTAGCGCATCAGGcga                        >  1:1441449/2‑98 (MQ=18)
                                                                            ggCCGGATAAGGCGTTCACGCCGCATCCGGCCTACAAGATCGTGCAAATTCAATATATTGCATGAGC‑‑‑CAC‑‑GTAGGCCTGATAAGCGTAGAGCATCa                             >  1:1116599/2‑96 (MQ=17)
                                                                            ggCCGGATAAGGCGTTCACGCCGCATCCGGCCTACAAAACTTTGCAAATTCAATATATTGC‑‑AATC‑‑‑TCCGTGTAGGCCTGATAAGCGTAGCGCATCAGGCa                         >  1:1239455/2‑100 (MQ=18)
                                                                             gCCGGATAAGGCGTTTACGCCGCATCCGGCCTACAAAATCGTGCTAATTCAATATATTGCAGAAAC‑‑‑CTT‑‑GTAGGCCTGATAAGCGTAGCGCATCAGGCag                        >  1:846466/1‑99 (MQ=18)
                                                                             gCCGGATAAGGCGTTCACGCCGCATCCGGCCTACAAGATCGTGCAAATTCAATATATTGCATGAGC‑‑‑CAC‑‑GTAGGCCTGATAAGCGTAGCGCATCAGGCag                        >  2:944556/1‑99 (MQ=17)
                                                                              ccggataaggcgTTCACGCCGCATCCGGCCTACAAAATCGCTTAAATTCAATATATTGCACCATC‑‑‑CTC‑‑GTAGGCCTGATAAGCGTAGCGCATCAGGCaat                       >  1:24749/9‑98 (MQ=21)
                                                                              ccGAATGCGGCGTTAACGCCTTATCCGGCCTACAAAATCGTGCAAATTCAATACATTGCAGAATT‑‑‑TTC‑‑GTAGGCCTGATAAGCGGAGCGCATCAGGCAtt                       >  1:997261/1‑100 (MQ=255)
                                                                                ggataaggcgTTTACGCCGCATCCGGCCTACAAATTCTGCAAATTCAACAAATTGCGAATCC‑‑‑CTT‑‑GTAGGCCTGATAAGCGTAGCGCATCAGGCAtttgc                    >  1:556211/7‑98 (MQ=21)
                                                                                ggataaggcgTTCACGCCGCATCCGGCCTACAAGATCGTGCAAATTCAATATATTGCATGAGC‑‑‑CAC‑‑GTAGGCCTGATAAGCGTAGCGCATCAGGCAgttt                     >  1:1111041/7‑100 (MQ=17)
                                                                                ggataaggcgTTCACGCCGCATCCGGCCTACAAAAACATGCAAATTCAATAGATTGCAGAGATT‑‑‑AT‑‑GTAGGCCTGATAAGCGTAGCGCATCAGGCAattt                     >  1:1017335/7‑100 (MQ=17)
                                                                                ggataaggcgTTCACGCCGCATCCGGCCTACAAAAACATGCAAATTCAATAGATTGCAGAGAT‑‑‑TAT‑‑GTAGGCCTGATAAGCGTAGCGCATCAGGCAattt                     >  2:229942/7‑100 (MQ=17)
                                                                                ggataaggcgTTCACGCCGCATCCAGCCTACAAAATTGTGCAAAATCAATGGATTGCACAGCT‑‑‑AGC‑‑GTAGGCCTGATAAGCGTAGCGCATCAGGCagatt                     >  1:1209252/7‑96 (MQ=37)
                                                                                 gataaggcgTTCACGCCGCATCCGGCCTACAAAAATCTTGCCAATTCAATATATTGCAGGACC‑‑‑AT‑‑GTAGGCCTGATAAGCGTAGCGCATCAGGCAgtttt                    >  1:2253507/6‑100 (MQ=18)
                                                                                 gataaggcgTTCACGCCGCATCCGGCCTACAAAATTGTGTAAATTCAATAAATTGATCGATC‑‑‑CAC‑‑GTAGGCCTGATAAGCGTAGCGCATCAGGCAgtttt                    >  2:942925/6‑100 (MQ=21)
                                                                                  aTGCGGCGTTAACGCCTTATCCGGCCTACAAAATCGTGCAAATTCAATACATTGCAGAATT‑‑‑TTC‑‑GTAGGCCTGATAAGCGGAGCGCATCAGGCATTtttg                   <  2:497555/100‑1 (MQ=255)
                                                                                  aTGCGGCGTTAACGCCTTATCCGGCCTACAAAATCGTGCAAATTCAATACATTGCAGAATT‑‑‑TTC‑‑GTAGGCCTGATAAGCGGAGCGCATCAGGCATTtttg                   <  2:636784/100‑1 (MQ=255)
                                                                                    gcggcgTTAACGCCTTATCCGGCCTACAAAATCGTGCAAATTCAATACATTGCAGAATT‑‑‑TTC‑‑GTAGGCCTGATAAGCGGAGCGCATCAGGCATTtttgcc                 <  2:1586733/100‑1 (MQ=255)
                                                                                    gcggcgTTAACGCCTTATCCGGCCTACAAAATCGTGCAAATTCAATACATTGCAGAATT‑‑‑TTC‑‑GTAGGCCTGATAAGCGGAGCGCATCAGGCATTtttgcc                 <  1:899507/100‑1 (MQ=255)
                                                                                     cggcgTTAACGCCTTATCCGGCCTACAAAATCGTGCAAATTCAATACATTGCAGAATT‑‑‑TTC‑‑GTAGGCCTGATAAGCGGAGCGCATCAGGCATTtttgcct                >  1:2264762/1‑100 (MQ=255)
                                                                                     cggcgTTAACGCCTTATCCGGCCTACAAAATCGTGCAAATTCAATACATTGCAGAATT‑‑‑TTC‑‑GTAGGCCTGATAAGCGGAGCGCATCAGGCATTtttgcct                >  2:1383096/1‑100 (MQ=255)
                                                                                          ttAACGCCTTATCCGGCCTACAAAATCGTGCAAATTCAATACATTGCAGAATT‑‑‑TTC‑‑GTAGGCCTGATATGCGGAGCGCATCAGGCACTtttgcctttgtc           <  1:1617243/100‑3 (MQ=255)
                                                                                          ttAACGCCTTATCCGGCCTACAAAATCGTGCAAATTCAATACATTGCAGAATT‑‑‑TTC‑‑GTAGGCCTGATAAGCGGAGCGCATCAGGCATTtttgcctttgcc           <  2:372409/100‑1 (MQ=255)
                                                                                          ttAACGCCTTATCCGGCCTACAAAATCGTGCAAATTCAATACATTGCAGAATT‑‑‑TTC‑‑GTAGGCCTGATAAGCGGAGCGCATCAGGCATTtttgcctttgcc           <  1:1851719/100‑1 (MQ=255)
                                                                                             aCGCCTTATCCGGCCTACAAAATCGTGCAAATTCAATACATTGCAGAATT‑‑‑TTC‑‑GTAGGCCTGATAAGCGGAGCGCATCAGGCATTTTTGCCTTTGCCATc        <  1:2015400/100‑1 (MQ=255)
                                                                                               gacTTATCCGGCCAACAAAATCGTGCAAATTAAATACATTGCTGAATT‑‑‑TTC‑‑GTAGGCCTGATAAGCGGAGCGCATCAGGCATTTTTGCCTTTGCCATCAg      >  1:1617611/3‑100 (MQ=25)
                                                                                                ccTTATCCGGCCTACAAAATCGTGCAAATTCAATACATTGCAGAATT‑‑‑TTC‑‑GTAGGCCTGATAAGCGGAGCGCATCAGGCATTTTTGCCTTTGCCATCAGt     >  1:1799683/1‑100 (MQ=255)
                                                                                                  ttATCCGGCCTACAAAATCGTGCAAATTCAATACATTGCAGAATT‑‑‑TTC‑‑GTAGGCCTGATAAGCGGAGCGCATCAGGCATTTTTGCCTTTGCCATCAGTCt   <  2:2379713/100‑1 (MQ=255)
                                                                                                   tATCCGGCCTACAAAATCGTGCAAATTCAATACATTGCAGAATT‑‑‑TTC‑‑GTAGGCCTGATAAGCGGAGCGCATCAGGCATTTTTGCCTTTGCCATCAGTCtt  <  1:187019/100‑1 (MQ=255)
                                                                                                   tATCCGGCCTACAAAATCGTGCAAATTCAATACATTGCAGAATT‑‑‑TTC‑‑GTAGGCCTGATAAGCGGAGCGCATCAGGCATTTTTGCCTTTGCCATCAGTCtt  >  1:1275498/1‑100 (MQ=255)
                                                                                                   |                                                                                                           
TGCTGCCGGATGGCGCTGCACGCAGAGGCACAGCGCCACCGGGTATACAGATTGTTGACAAAGTGCGCTTTGTTTATGCCGGATGCGGCGTTAACGCCTTATCCGGCCTACAAAATCGTGCAAATTCAATACATTGCAGAATT‑‑‑TTC‑‑GTAGGCCTGATAAGCGGAGCGCATCAGGCATTTTTGCCTTTGCCATCAGTCTT  >  NC_000913/3705886‑3706084

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: