Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009133,705,9850TC43.3% 197.3 / 285.6 217intergenic (‑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 (49/74);  new base C (93/1);  total (142/75)
Fisher's exact test for biased strand distribution p-value = 1.85e-23
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.99e-01

Only 100 of 218 total aligned reads displayed.

GCTGCCGGATGGCGCTGCACGCAGAGGCACAGCGCCACCGGGTATACAGATTGTTGACAAAGTGCGCTTTGTTTATGCCGGATGCGGCGTTAACGCCTTATCCGGCCTAC‑AAAATCGTGCAAATTCAATACATTGCA‑GA‑ATT‑‑‑TTC‑‑GTAGGCCTGATAAGCGGAGCGCATCAGGCATTTTTGCCTTTGCCATCAGTCTT  >  NC_000913/3705887‑3706084
                                                                                                  |                                                                                                           
 ctgcCGGATGGCGCTGCACGCAGAGGCACAGCGCCACCGGGTATACAGATTGTTGACAAAGTGCGCTTTGTTTATGCCGGATGCGGCGTTAACGCCTTAt                                                                                                           <  2:2071471/100‑1 (MQ=255)
  tgcCGGATGGCGCTGCACGCAGAGGCACAGCGCCACCGGGTATACAGATTGTTGACAAAGTGCGCTTTGTTTATGCCGGATGCGGCGTTAACGCCTTATc                                                                                                          >  2:1347484/1‑100 (MQ=255)
  tgcCGGATGGCGCTGCACGCAGAGGCACAGCGCCACCGGGTATACAGATTGTTGACAAAGTGCGCTTTGTTTATGCCGGATGCGGCGTTAACGCCTTATc                                                                                                          >  1:610954/1‑100 (MQ=255)
  tgcCGGATGGCGCTGCACGCAGAGGCACAGCGCCACCGGGTATACAGATTGTTGACAAAGTGCGCTTTGTTTATGCCGGATGCGGCGTTAACGCCTTATc                                                                                                          >  2:2970780/1‑100 (MQ=255)
      ggATGGCGCTGCACGCAGAGGCACAGCGCCACCGGGTATACAGATTGTTGACAAAGTGCGCTTTGTTTATGCCGGATAAGGCGTTAACGCCGCATCCGGc                                                                                                      <  2:2181276/100‑1 (MQ=37)
           gcgcTGCACGCAGAGGCACAGCGCCACCGGGTATACAGATTGTTGACAAAGTGCGCTTTGTTTATGCCGGATGCGGCGTTAACGCCTTATCCGGCCTAC‑a                                                                                                >  1:1464307/1‑100 (MQ=255)
           gcgcTGCACGCAGAGGCACAGCGCCACCGGGTATACAGATTGTTGACAAAGTGCGCTTTGTTTATGCCGGATGCGGCGTTAACGCCTTATCCGGCCTAC‑a                                                                                                <  2:2229045/100‑1 (MQ=255)
           gcgcTGCACGCAGAGGCACAGCGCCACCGGGTATACAGATTGTTGACAAAGTGCGCTTTGTTTATGCCGGATGCGGCGTTAACGCCTTATCCGGCCTAC‑a                                                                                                >  2:2419358/1‑100 (MQ=255)
                 cACGCAGAGGCACAGCGCCACCGGGTATACAGATTGTTGACAAAGTGCGCTTTGTTTATGCCGGATGCGGCGTTAACGCCTTATCCGGCCTAC‑AAAATCg                                                                                          <  2:715507/100‑1 (MQ=255)
                  aCGCAGAGGCACAGCGCCACCGGGTATACAGATTGTTGACAAAGTGCGCTTTGTTTATGCCGGATGCGGCGTTAACGCCTTATCCGGCCTAC‑AAAATCGt                                                                                         <  2:1717997/100‑1 (MQ=255)
                   cGCAGAGGCACAGCGCCACCGGGTATACAGATTGTTGACAAAGTGCGCTTTGTTTATGCCGGATGCGGCGTTAACGCCTTATCCGGCCTAC‑AAAATCGTg                                                                                        >  1:2858014/1‑100 (MQ=255)
                   cGCAGAGGCACAGCGCCACCGGGTATACAGATTGTTGACAAAGTGCGCTTTGTTTATGCCGGATGCGGCGTTAACGCCTTATCCGGCCTAC‑AAAATCGTg                                                                                        <  1:1966050/100‑1 (MQ=255)
                   cGCAGAGGCACAGCGCCACCGGGTATACAGATTGTTGACAAAGTGCGCTTTGTTTATGCCGGATGCGGCGTTAACACCTTATCCGGCCTAC‑AAAATCGTg                                                                                        >  2:1391728/1‑100 (MQ=255)
                    gCAGAGGCACAGCGCCACCGGGTATACAGATTGTTGACAAAGTGCGCTTTGTTTATGCCGGATGCGGCGTTAACGCCTTATCCGGCCTAC‑AAAATCGTGc                                                                                       <  2:904642/100‑1 (MQ=255)
                    gCAGAGGCACAGCGCCACCGGGTATACAGATTGTTGACAAAGTGCGCTTTGTTTATGCCGGATGCGGCGTTAACGCCTTATCCGGCCTAC‑AAAATCGTGc                                                                                       <  1:43043/100‑1 (MQ=255)
                        agGCACAGCGCCACCGGGTATACAGATTGTTGACAAAGTGCGCTTTGTTTATGCCGGATGCGGCGTTAACGCCTTATCCGGCCTAC‑AAAATCGTGCAAAt                                                                                   <  1:2939124/100‑1 (MQ=255)
                        agGCACAGCGCCACCGGGTATACAGATTGTTGACAAAGTGCGCTTTGTTTATGCCGGATGCGGCGTTAACGCCTTATCCGGCCTAC‑AAAATCGTGCAAAt                                                                                   <  2:229758/100‑1 (MQ=255)
                          gCACAGCGCCACCGGGTATACAGATTGTTGACAAAGTGCGCTTTGTTTATGCCGGATGCGGCGTTAACGCCTTATCCGGCCTAC‑AAAATCGTGCAAATTc                                                                                 <  1:2240940/100‑1 (MQ=255)
                           cacaGCGCCACCGGGTATACAGATTGTTGACAAAGTGCGCTTTGTTTATGCCGGATGCGGCGTTAACGCCTTATCCGGCCTAC‑AAAATCGTGCAAATTCa                                                                                >  1:1837613/1‑100 (MQ=255)
                           cacaGCGCCACCGGGTATACAGATTGTTGACAAAGTGCGCTTTGTTTATGCCGGATGCGGCGTTAACGCCTTATCCGGCCTAC‑AAAATCGTGCAAATTCa                                                                                >  2:2132377/1‑100 (MQ=255)
                             caGCGCTACCGGGTATACAGATTGTTGACAAAGTGCGCTTTGTTTATGCCGGATGCGGCGTTAACGCCTTATCCGGCCTAC‑AAAATCGTGCAAATTCAAt                                                                              <  1:1286753/100‑1 (MQ=255)
                             caGCGCCACCGGGTATACAGATTGTTGACAAAGTGCGCTTTGTTTATGCCGGATGCGGCGTTAACGCCTTATCCGGCCTAC‑AAAATCGTGCAAATTCAAt                                                                              <  2:2857986/100‑1 (MQ=255)
                              aGCGCCACCGGGTATACAGATTGTTGACAAAGTGCGCTTTGTTTATGCCGGATGCGGCGTTAACGCCTTATCCGGCCTAC‑AAAATCGTGCAAATTCAATa                                                                             <  1:979603/100‑1 (MQ=255)
                                    aCCGGGTATACAGATTGTTGACAAAGTGCGCTTTGTTTATGCCGGATGCGGCGTTAACGCCTTATTCGGCCTAC‑AAAATCGTGCAAATTCAATACATTg                                                                        >  2:1404725/1‑99 (MQ=255)
                                     ccGGGTATACAGATTGTTGACAAAGTGCGCTTTGTTTATGCCGGATGCGGCGTTAACGCCTTATCCGGCCTAC‑AAAATCGTGCAAATTCAATACATTGCa                                                                      >  1:1375973/1‑100 (MQ=255)
                                            taCAGATTGTTGACAAAGTGCGCTTTGTTTATGCCGGATGCGGCGTTAACGCCTTATCCGGCCTAC‑TAAATCGTGCAAATTCAATACATTGCA‑GA‑Att‑‑‑tt                                                          <  2:2038711/100‑1 (MQ=255)
                                            taCAGATTGTTGACAAAGTGCGCTTTGTTTATGCCGGATGCGGCGTTAACGCCTTATCCGGCCTAC‑AAAATCGTGCAAATTCAATACATTGCA‑GA‑Att‑‑‑tt                                                          >  2:2933668/1‑100 (MQ=255)
                                             aCAGATTGTTGACAAAGTGCGCTTTGTTTATGCCGGATGCGGCGTTAACGCCTTATCCGGCCTAC‑AAAATCGTGCAAATTCAATACATTGCA‑GA‑ATT‑‑‑TTc                                                         <  2:2212386/100‑1 (MQ=255)
                                              cAGATTGTTGACAAAGTGCGCTTTGTTTATGCCGGATGCGGCGTTAACGCCTTATCCGGCCTAC‑AAAATCGTGCAAATTCAATACATTGCA‑GA‑ATT‑‑‑TTC‑‑g                                                      <  2:2478712/100‑1 (MQ=255)
                                              cAGATTGTTGACAAAGTGCGCTTTGTTTATGCCGGATGCGGCGTTAACGCCTTATCCGGCCTAC‑AAAATCGTGCAAATTCAATACATTGCA‑GA‑ATT‑‑‑TTC‑‑g                                                      <  2:2369302/100‑1 (MQ=255)
                                                gATTGTTGACAAAGTGCGCTTTGTTTATGCCGGATGCGGCGTTAACGCCTTATCCGGCCTAC‑AAAATCGTGCAAATTCAATACATTGCA‑GA‑ATT‑‑‑TTC‑‑GTa                                                    <  1:1924279/100‑1 (MQ=255)
                                                 aTTGTTGACAAAGTGCGCTTTGTTTATGCCGGATGCGGCGTTAACGCCTTATCCGGCCTAC‑AAAATCGTGCAAATTCAATACATTGCA‑GA‑ATT‑‑‑TTC‑‑GTAg                                                   <  2:1922372/100‑1 (MQ=255)
                                                  ttgttgACAAAGTGCGCTTTGTTTATGCCGGATGCGGCGTTAACGCCTTATCCGGCCTAC‑AAAATCGTGCAAATTCAATACATTGCA‑GA‑ATT‑‑‑TTC‑‑GTAgg                                                  <  2:426588/100‑1 (MQ=255)
                                                    gttgACAAAGTGCGCTTTGTTTATGCCGGATGCGGCGTTAACGCCTTATCCGGCCTAC‑AAAATCGTGCAAATTCAATACATTGCA‑GA‑ATT‑‑‑TTC‑‑GTAGGcc                                                <  2:390663/100‑1 (MQ=255)
                                                     ttgACAAAGTGCGCTTTGTTTATGCCGGATGCGGCGTTAACGCCTTATCCGGCCTAC‑AAAATCGTGCAAATTCAATACATTGCA‑GA‑ATT‑‑‑TTC‑‑GTAGGCCt                                               <  1:2044402/100‑1 (MQ=255)
                                                     ttgACAAAGTGCGCTTTGTTGATGCCGGATGCGGCGTTAGCGCCTTATCCGGCCTAC‑AAAATCGTGCAAATTCAATACATTGCA‑GA‑ATT‑‑‑TTC‑‑GTAGGCCt                                               <  2:2772865/100‑1 (MQ=255)
                                                       gACAAAGTGCGCTTTGTTTATGCCGGATGCGGCGTTAACGCCTTATCCGGCCTAC‑AAAATCGTGCAAATTCAATACATTGCA‑GA‑ATT‑‑‑TTC‑‑GTAGGCCTGa                                             <  2:2683300/100‑1 (MQ=255)
                                                       aaCAAAGCGCACTTTGTTTATGCCGGATGCGGCGTTAACGCCTTATCCGGCCTAC‑AAAATCGTGCAAATTCAATACATTGCA‑GA‑ATT‑‑‑TTC‑‑GTAGGCCTGa                                             >  1:928465/2‑100 (MQ=255)
                                                              tGCGCTTTGTTTATGCCGGATGCGGCGTTAACGCCTTATCCGGCCTAC‑AAAATCGTGCAAATTCAATACATTGCA‑GA‑ATT‑‑‑TTC‑‑GTAGGCCTGATAAGCgg                                      <  1:1640142/100‑1 (MQ=255)
                                                              tGCGCTTTGTTTATGCCGGATGCGGCGTTAACGCCTTATCCGGCCTAC‑AAAATCGTGCAAATTCAATACATTGCA‑GA‑ATT‑‑‑TTC‑‑GTAGGCCTGATAAGCgg                                      <  2:3061109/100‑1 (MQ=255)
                                                               gcgcTTTGTTTATGCCGGATGCGGCGTTAACGCCTTATCCGGCCTAC‑AAAATCGTGCAAATTCAATACATTGCA‑GA‑ATT‑‑‑TTC‑‑GTAGGCCTGAtgcgctcc                                     >  1:460638/1‑93 (MQ=255)
                                                                  cTTTGTTTATGCCGGATGCGGCGTTAACGCCTTATCCGGCCTAC‑AAAATCGTGCAAATTCAATACATTGCA‑GA‑ATT‑‑‑TTC‑‑GTAGGCCTGATAAGCGGAgcg                                  <  1:1788351/100‑1 (MQ=255)
                                                                  cTTTGTTTATGCCGGATGCGGCGTTAACGCCTTATCCGGCCTAC‑AAAATCGTGCAAATTCAATACATTGCA‑GA‑ATT‑‑‑TTC‑‑GTAGGCCTGATAAGCGGAgcg                                  >  1:1369024/1‑100 (MQ=255)
                                                                  cTTTGTTTATGCCGGATGCGGCGTTAACGCCTTATCCGGCCTAC‑AAAATCGTGCAAATTCAATACATTGCA‑GA‑ATT‑‑‑TTC‑‑GTAGGCCTGATAAGCGGAgcg                                  <  1:1879093/100‑1 (MQ=255)
                                                                    ttttgtaggCCGGATAAGGCGTTCACGCCGCATCCGGCCTAC‑AAAATTGTGTAAATTCAATAAATTGAT‑CG‑ATC‑‑‑CAC‑‑GTAGGCCTGATAAGCGTAGCGCa                                >  2:1769391/9‑100 (MQ=21)
                                                                    ttttgtaggCCGGATAAGGCGTTCACGCCGCATCCGGCCTAC‑AAAAGCACGCAAACTCAATATATTGCA‑GAGATC‑‑‑AT‑‑‑GTAGGCCTGATAAGCGTAGCGCa                                >  1:676191/9‑100 (MQ=17)
                                                                     tGTTTATGCCGGATGCGGCGTTAACGCCTTATCCGGCCTAC‑AAAATCGTGCAAATTCAATACATTGCA‑GA‑ATT‑‑‑TTC‑‑GTAGGCCTGATAAGCGGAGCGCAt                               >  1:1605329/1‑100 (MQ=255)
                                                                      ttgtaggCCGGATAAGGCGTTTACGCCGCATCTGGCCTAC‑AAATTCTTCCAAATTCAATATATTGCA‑GA‑AAT‑‑‑CAT‑‑GTAGGCCTGATAAGCGTAGCGCATc                              >  1:2376834/7‑100 (MQ=18)
                                                                      ttgtaggCCGGATAAGGCGTTCACGCCGCATCCGGCCTACAAAAATCTTGCCAATTCAATATATTGCA‑GG‑ACC‑‑‑AT‑‑‑GTAGGCCTGATAAGCGTAGCGCATc                              >  1:1200845/7‑100 (MQ=18)
                                                                      ttgtaggCCGGATAAGGCGTTCACGCCGCATCCGGCCTAC‑AAGATCGTGCAAATTCAATATATTGCA‑TG‑AGC‑‑‑CAC‑‑GTAGGCCTGATAAGCGTAGCGCATc                              >  1:1046052/7‑100 (MQ=17)
                                                                      ttgtaggCCGGATAAGGCGTTCACGCCGCATCCGGCATAC‑AAAATTGTGTAAATTCAATATATTGAT‑CG‑ATC‑‑‑CAC‑‑GTAGGCCTGATAAGCGTAGCGCATc                              >  2:1183400/7‑100 (MQ=18)
                                                                        gtaggCCGGATAAGGCGTTTACGCCGCATCCGGCATAC‑AAAATCGTGCTAATTCAATATATTGCA‑GA‑AAC‑‑‑CTT‑‑GTAGGCCTGATAAGCGTAGCGCATCAg                            >  2:2590805/5‑100 (MQ=18)
                                                                         taggCCGGATAAGGCGTTTACGCCGCATCCGGCCTAC‑AAAATCGTGCAAATTCAACATATTGCA‑AT‑TCT‑‑‑CTT‑‑GTAGGCCTGATAAGCGTAGCGCATCAgg                           >  1:103512/4‑100 (MQ=18)
                                                                         taggCCGGATAAGGCGTTCACGCCGCATCCGGCCTACAAAAATCTTGCCAATTCAATATATTGCA‑GG‑ACC‑‑‑AT‑‑‑GTAGGCCTGATAAGCGTAGCGCATCAgg                           >  1:2620726/4‑100 (MQ=18)
                                                                         taggCCGGATAAGGCGTTCACGCCGCATCCGGCCTAC‑AAAGTCGTATAAATTCAATGTATTGCA‑GT‑TTA‑‑‑TTG‑‑GTAGGCCTGATAAGCGCAGCGCATCAgg                           >  2:3001535/4‑100 (MQ=21)
                                                                         taggCCGGATAAGGCGTTCACGCCGCATCCGGCCTAC‑AAAACTTTGCAAATTCAATATATTGC‑‑‑A‑ATC‑‑‑TCCGTGTAGGCCTGATAAGCGTAGCGCATCAgg                           >  1:2260622/4‑100 (MQ=18)
                                                                         tATGCCGGATGCGGCGTTAACGCCTTATCCGGCCTAC‑AAAATCGTGCAAATTCAATACATTGCA‑GA‑ATT‑‑‑TTC‑‑GTAGGCCTGATAAGCGGAGCGCATCAgg                           <  2:2636986/100‑1 (MQ=255)
                                                                          aggCCGGATAAGGCGTTTACGCCGCATCCGGCCTAC‑AAAATCGTGCTAATTCAATATATTGCA‑GA‑AAC‑‑‑CTT‑‑GTAGGCCTGATAAGCGTAGCGCATCAGGc                          >  2:78956/3‑100 (MQ=17)
                                                                          aggCCGGATAAGGCGTTCACGCCGCATCCGGCCTAC‑AAAACTTTGCAAATTCAATATATTGC‑‑‑A‑ATC‑‑‑TCCGTGTAGGCCTGATAAGCGTAGCGCATCAGGc                          >  1:1760422/3‑100 (MQ=18)
                                                                          aggCCGGATAAGGCGTTCACGCCGCATCCGGCCTAC‑AAAAACCTGCAAATTCAATATATTGCA‑GG‑ATA‑‑‑CGT‑‑GTAGGCCTGATAAGCGTAGCGCATCAGGc                          >  2:2741117/3‑100 (MQ=18)
                                                                          aggCCGGATAAGGCGTTCACGCCGCATCCGGCCTAC‑AAAAACATGCAAATTCAATAGATTGCA‑GA‑GAT‑‑‑TAT‑‑GTAGGCCTGATAAGCGTAGCGCATCAGGc                          >  2:2109148/3‑100 (MQ=17)
                                                                           ggCCGGATAAGGCGTTTACGCCGCATCCGGCCTAC‑AAAAGCATGCAAATTCAATATATTGCA‑GA‑GTT‑‑‑CAC‑‑GTAGGCCTGATAAGCGTAGCGCATCAGGCa                         >  2:474452/2‑100 (MQ=17)
                                                                           ggCCGGATAAGGCGTTTACGCCGCATCCGGCCTAC‑AAAAGCATGCAAATTCAATATATTGCA‑GA‑GTT‑‑‑CAC‑‑GTAGGCCTGATAAGCGTAGCGCATCAGGCa                         >  2:622228/2‑100 (MQ=17)
                                                                           ggCCGGATAAGGCGTTCACGCCGCATCCGGCCTACAAAAATCTTGCCAATTCAATATATTGCA‑GG‑ACC‑‑‑AT‑‑‑GTAGGCCTGATAAGCGTAGCGCATCAGGCa                         >  1:777102/2‑100 (MQ=18)
                                                                           ggCCGGATAAGGCGTTCACGCCGCATCCGGCCTACAAAAATCTTGCCAATTCAATATATTGCA‑GG‑ACC‑‑‑AT‑‑‑GTAGGCCTGATAAGCGTAGCGCATCAGGCa                         >  1:696460/2‑100 (MQ=18)
                                                                           ggCCGGATAAGGCGTTCACGCCGCATCCGGCCTACAAAAATCTTGCCAATTCAATATATTGCA‑GG‑ACC‑‑‑AT‑‑‑GTAGGCCTGATAAGCGTAGCGCATCAGGCa                         >  1:1835910/2‑100 (MQ=18)
                                                                           ggCCGGATAAGGCGTTCACGCCGCATCCGGCCTAC‑AAAATCGCTTAAATTCAATATATTGCA‑CC‑ATC‑‑‑CTC‑‑GTAGGCCTGATAAGCGTAGCGCATCAGGCa                         >  2:3055622/2‑100 (MQ=21)
                                                                           ggCCGGATAAGGCGTTCACGCCGCATCCGGCCTAC‑AAAACTTTGCAAATTCAATATATTGCA‑AT‑CTC‑‑‑CGT‑‑GTAGGCCTGATAAGCGTAGCGCATCAGGCa                         >  2:2820538/2‑100 (MQ=18)
                                                                            gCCGGATGCGGCGTTAACGCCTTATCCGGCCTAC‑AAAATCGTGCAAATTCAATACATTGCA‑GA‑ATT‑‑‑TTC‑‑GTAGGCCTGATAAGCGGAGCGCATCAGGCAt                        >  1:2808389/1‑100 (MQ=255)
                                                                            gCCGGATAAGGCGTTTACGCCGCATCCGGCCTAC‑AAAATCGTGCTAATTCAATATATTGCA‑GA‑AAC‑‑‑CTT‑‑GTAGGCCTGATAAGCGTAGCGCATCAGGCag                        >  2:1452139/1‑99 (MQ=17)
                                                                            gCCGGATAAGGCGTTTACGCCGCATCCGGCCTAC‑AAAACAGTGCAAATTCAATATATTGCA‑AT‑CTT‑‑‑CTG‑‑GTAGGCCTGATAAGCATCGCGCATCAGGCaa                        >  1:96037/1‑99 (MQ=18)
                                                                            gCCGGATAAGGCGTTCACGCCGCATCCGGCCTAC‑AAAATCGCTTAAATTCAATATATTGCA‑CC‑ATC‑‑‑CTC‑‑GTAGGCCTGATAAGCGTAGCGCATCAGGCaa                        >  2:84063/1‑99 (MQ=21)
                                                                            gCCGGATAAGGCGTTCACGCCGCATCCGGCCTAC‑AAAATCGCTTAAATTCAATATATTGCA‑CC‑ATC‑‑‑CTC‑‑GTAGGCCTGATAAGCGTAGCGCATCAGGCaa                        >  1:1001982/1‑99 (MQ=21)
                                                                             ccggataaggcgTTCACGCCGCATCCGGCCTAC‑AAAATCGCTTAAATTCAATATATTGCA‑CC‑ATC‑‑‑CTC‑‑GTAGGCCTGATAAGCGTAGCGCATCAGGCaat                       >  1:41349/9‑98 (MQ=21)
                                                                             ccGGATGCGGCGTTAACGCCTTATCCGGCCTAC‑AAAATCGTGCAAATTCAATACATTGCA‑GA‑ATT‑‑‑TTC‑‑GTAGGCCTGATAAGCGGAGCGCATCAGGCAtt                       <  2:974133/100‑1 (MQ=255)
                                                                              cggataaggcgTTTACGCCGCATCCGGCCTAC‑AAAAGCATGCAAATTCAATATATTGCA‑GA‑GTT‑‑‑CAC‑‑GTAGGCCTGATAAGCGTAGCGCATCAGGCAgtt                      >  2:1825193/8‑100 (MQ=17)
                                                                              cggataaggcgTTCACGCCTTATCCGACCCAC‑AAAGTCTTGCAAAATCAATAAATTGCAGGA‑ACG‑‑‑AT‑‑‑GTAGGCCTGATAAGCGTAGCGCATCAGGCaact                      >  2:1047157/8‑97 (MQ=25)
                                                                              cggataaggcgTTCACGCCGCATCCGGCCTAC‑AAAACTTTGCAAATTCAATATATTGCA‑AT‑CTC‑‑‑CGT‑‑GTAGGCCTGATAAGCGTAGCGCATCAGGCaatt                      >  1:891575/8‑97 (MQ=18)
                                                                               ggataaggcgTTTACGCCGCATCCGGCCTAC‑AAAAGCGTGCAAATTCAATACATTGCATGG‑GCC‑‑‑AT‑‑‑GTAGGCCTGATAAGCGTAGCGCATCAGGCAattt                     >  2:581699/7‑100 (MQ=18)
                                                                               ggataaggcgTTTACGCCGCATCCGGCCTAC‑AAAAGCATGCAAATTCAATATATTGCA‑GA‑GTT‑‑‑CAC‑‑GTAGGCCTGATAAGCGTAGCGCATCAGGCAgttt                     >  1:2612360/7‑100 (MQ=17)
                                                                               ggataaggcgTTTACGCCGCATCCGGCCTAC‑AAAAGCATGCAAATTCAATATATTGCA‑GA‑GTT‑‑‑CAC‑‑GTAGGCCTGATAAGCGTAGCGCATCAGGCAgttt                     >  1:2107367/7‑100 (MQ=17)
                                                                               ggataaggcgTTTACGCCGCATCCGGCCTAC‑AAAAATGTGCAAATTCAATAAATTGC‑‑‑A‑ATTCAACTT‑‑GTAGGCCTGATAAGCGCAGCGCATCAGGCaatt                      >  1:1127671/7‑97 (MQ=21)
                                                                               ggataaggcgTTCACGCCGCATCCGGCCTAC‑AAGATCGTGCAAATTCAATATATTGCA‑TG‑AGC‑‑‑CAC‑‑GTAGGCCTGATAAGCGTAGCGCATCAGGCAgttt                     >  2:530611/7‑100 (MQ=17)
                                                                               ggataaggcgTTCACGCCGCATCCGGCCTAC‑AAAGTCGTATACATTCAATGTATTGCA‑GT‑TTA‑‑‑TTG‑‑GTAGGCCTGATAAGCGCAGCGCATCACGCaat                       >  1:3074054/7‑96 (MQ=21)
                                                                               ggataaggcgTTCACGCCGCATCCGGCCTAC‑AAAAACATGCAAATTCAATAGATTGCA‑GAGATT‑‑‑AT‑‑‑GTAGGCCTGATAAGCGTAGCGCATCAGGCAattt                     >  1:3046099/7‑100 (MQ=17)
                                                                               ggataaggcgTTCACGCCGCATCCGGCCTAC‑AAAAACATGCAAATTCAATAGATTGCA‑GA‑GAT‑‑‑TAT‑‑GTAGACCTGATAAGCGTAGCGCATCAGGCAattt                     >  2:1859892/7‑100 (MQ=17)
                                                                                gataaggcgTTCACGCCGCATCCGGCCTAC‑AAAATTGTGTAAATTCAATAAATTGAT‑CG‑ATC‑‑‑CAC‑‑GTAGGCCTGATAAGCGTAGCGCATCAGGCAgtttt                    >  1:1487741/6‑100 (MQ=21)
                                                                                gATGCGGCGTTAACGCCTTATCCGGCCTAC‑AAAATCGTGCAAATTCAATACATTGCA‑GA‑ATT‑‑‑TTC‑‑GTAGGCCTGATAAGCGGAGCGCATCAGGCAttttt                    >  1:588275/1‑100 (MQ=255)
                                                                                  taaggcgTTTACGCCGCATCCGGCCTAC‑AAAAGCGTGCAAATTCAATACATTGCATGG‑GCC‑‑‑AT‑‑‑GTAGGCCTGATAAGCGTAGCGCATCAGGCAattttac                  >  2:3040789/4‑98 (MQ=18)
                                                                                  taaggcgTTCACGCCGCATCCGGCCTACAAAAATCTTGCCAATTCAATATATTGCA‑GG‑ACC‑‑‑AT‑‑‑GTAGGCCTGATAAGCGTAGCGCATCAGGCAGTtttgc                  >  1:1687843/4‑100 (MQ=18)
                                                                                    cggcgTTAACGCCTTATCCGGCCTAC‑AAAATCGTGCAAATTCAATACATTGCA‑GA‑ATT‑‑‑TTC‑‑GTAGGCCTGATAAGCGGAGCGCATCAGGCATTtttgcct                >  1:636043/1‑100 (MQ=255)
                                                                                    cggcgTTAACGCCTTATCCGGCCTAC‑AAAATCGTGCAAATTCAATACATTGCA‑GA‑ATT‑‑‑TTC‑‑GTAGGCCTGATAAGCGGAGCGCATCAGGCATTtttgcct                >  2:896604/1‑100 (MQ=255)
                                                                                      gcgTTAACGCCTTATCCGGCCTAC‑AAAATCGTGCAAATTCAATACATTGCA‑GA‑ATT‑‑‑TTC‑‑GTAGGCCTGATAAGCGGAGCGCATCAGGCATTtttgccttt              <  1:2484349/100‑1 (MQ=255)
                                                                                       cgTTAACGCCTTATCCGGCCTAC‑AAAATCGTGCAAATTCCATACATTGCA‑GA‑ATT‑‑‑TTC‑‑GTAGGCCTGATAAGCGGAGCGCATCAGGCATTtttgcctttg             <  1:2933698/100‑1 (MQ=255)
                                                                                         ttAACGCCTTATCCGGCCTAC‑AAAATCGTGCAAATTCAATACATTGCA‑GA‑ATT‑‑‑TTC‑‑GTAGGCCTGATAAGCGGAGCGCATCAGGCATTtttgcctttgcc           <  2:1613808/100‑1 (MQ=255)
                                                                                          tAACGCCTTATCCGGCCTAC‑AAAATCGTGCAAATTCAATACATTGCA‑GA‑ATT‑‑‑TTC‑‑TTTGTTTTTTTAAGCGGAGCGCATCAGGCATTTTTGCCTTTGCCa          >  2:2337558/1‑100 (MQ=21)
                                                                                          tAACGCCTTATCCGGCCTAC‑AAAATCGTGCAAATTCAATACATTGCA‑GA‑ATT‑‑‑TTC‑‑GTAGGCCTGATAAGCGGAGCGCATCAGGCATTTTTGCCTTTGCCa          <  1:2176383/100‑1 (MQ=255)
                                                                                              gCCTTATCCGGCCTAC‑AAAATCGTGCAAATTCAATACATTGCA‑GA‑ATT‑‑‑TTC‑‑GTAGGCCTGATAAGCGGAGCGCATCAGGCATTTTTGCCTTTGCCATCAg      >  1:2297876/1‑100 (MQ=255)
                                                                                               ccTTATCCGGCCTAC‑AAAATCGTGCAAATTCAATACATTGCA‑GA‑ATT‑‑‑TTC‑‑GTAGGCCTGATAAGCGGAGCGCATCAGGCATTTTTGCCTTTGCCATCAGt     >  2:683446/1‑100 (MQ=255)
                                                                                                  tATCCGGCCTAC‑ACAATCGTGCAAATTCAATACATTGCA‑GA‑ATT‑‑‑TTC‑‑GTAGGCCTGATAAGCGGACCGCATCAGGCATTTTTGCCTTTGCCATCAGTCtt  >  2:1910829/1‑100 (MQ=255)
                                                                                                  |                                                                                                           
GCTGCCGGATGGCGCTGCACGCAGAGGCACAGCGCCACCGGGTATACAGATTGTTGACAAAGTGCGCTTTGTTTATGCCGGATGCGGCGTTAACGCCTTATCCGGCCTAC‑AAAATCGTGCAAATTCAATACATTGCA‑GA‑ATT‑‑‑TTC‑‑GTAGGCCTGATAAGCGGAGCGCATCAGGCATTTTTGCCTTTGCCATCAGTCTT  >  NC_000913/3705887‑3706084

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