Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	3,705,971	0	C	A	22.1%	249.1 / 61.5	95	intergenic (‑181/+127)	dppB/dppA	dipeptide/heme ABC transporter permease/dipeptide/heme ABC transporter periplasmic binding protein; dipeptide chemotaxis receptor
Reads supporting (aligned to +/- strand):  ref base C (38/36);  new base A (20/1);  total (58/37)
Fisher's exact test for biased strand distribution p-value = 2.26e-04
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 2.39e-02

AGGGAATGAAGTATTGCTGCCGGATGGCGCTGCACGCAGAGGCACAGCGCCACCGGGTATACAGATTGTTGACAAAGTGCGCTTTGTTTATGCCGGATGCGGCGTTAACGCCTTATCCGGCCTACAAAATCGTGCAAATTCAATACATTGC‑‑A‑GAATT‑‑TTCGTAGGCCTGATAAGCGGAGCGCATCAGGCATTT  >  NC_000913/3705872‑3706064                                                                                                    |                                                                                                   aGGGAATGAAGTATTGCTGCCGGATGGCACTGCACGCAGAGGCACAGCGCCACCGGGTATACAGATTGTTGACAAAGTGCGCTTTGTTTATGCCGGATgc                                                                                                    <  1:1105132/100‑1 (MQ=255)  gggAATGAAGTATTGCTGCCGGATGGCGCTGCACGCAGAGGCACAGCGCCACCGGGTATACAGATTGTTGACAAAGTGCGCTTTGTTTATGCCGGATgcg                                                                                                   >  1:1663205/1‑100 (MQ=255)     aaTGAAGTATTGCTGCCGGATGGCGCTGCACGCAGAGGCACAGCGCCACCGGGTATACAGATTGTTGACAAAGTGCGCTTTGTTTATGCCGGATgcggcg                                                                                                >  2:1217027/1‑100 (MQ=255)      aTGAAGTATTGCTGCCGGATGGCGCTGCACGCAGAGGCACAGCGCCACCGGGTATACAGATTGTTGACAAAGTGCGCTTTGTTTATGCCGGATGCGGCGt                                                                                               <  2:1813762/100‑1 (MQ=255)       tGAAGTATTGCTGCCGGATGGCGCTGCACGCAGAGGCACAGCGCCACCGGGTATACAGATTGTTGACAAAGTGCGCTTTGTTTATGCCGGATGCGGCGtt                                                                                              <  1:109338/100‑1 (MQ=255)        gAAGTATTGCTGCCGGATGGCGCTGCACGCAGAGGCACAGCGCCACCGGGTATACAGATTGTTGACAAAGTGCGCTTTGTTTATGCCGGATGCGGCGTTa                                                                                             >  1:1685349/1‑100 (MQ=255)           gTATTGCTGCCGGATGGCGCTGCACGCAGAGGCACAGCGCCACCGGGTATACAGATTGTTGACAAAGTGCGCTTTGTTTATGCCGGATGCGGCGTTAACg                                                                                          <  2:1430295/100‑1 (MQ=255)                gctgcCGGATGGCGCTGCACGCAGAGGCACAGCGCCACCGGGTATACAGATTGTTGACAAAGTGCGCTTTGTTTATGCCGGATGCGGCGTTAACGCCTTa                                                                                     >  2:297678/1‑100 (MQ=255)                gctgcCGGATGGCGCTGCACGCAGAGGCACAGCGCCACCGGGTATACAGATTGTTGACAAAGTGCGCTTTGTTTATGCCGGATGCGGCGTTAACGCCTTa                                                                                     <  2:1681662/100‑1 (MQ=255)                    ccGGATGGCGCTGCACGCAGAGGCACAGCGCCACCGGGTATACAGATTGTTGACAAAGTGCGCTTTGTTTATGCCGGATGCGGCGTTAACGCCTTATCCg                                                                                 >  1:821922/1‑100 (MQ=255)                     cGGATGGCGCTGCACGCAGAGGCACAGCGCCACCGGGTATACAGATTGTTGACAAAGTGCGCTTTGTTTATGCCGGATGCGGCGTTAACGCCTTATCCgg                                                                                >  1:1236090/1‑100 (MQ=255)                      ggATGGCGCTGCACGCAGAGGCACAGCGCCACCGGGTATACAGATTGTTGACAAAGTGCGCTTTGTTTATGCCGGATGCGGCGTTAACGCCTTATCCGGc                                                                               >  1:1702187/1‑100 (MQ=255)                      ggATGGCGCTGCACGCAGAGGCACAGCGCCACCGGGTATACAGATTGTTGACAAAGTGCGCTTTGTTTATGCCGGATGCGGCGTTAACGCCTTATCCGGc                                                                               <  1:264075/100‑1 (MQ=255)                       gATGGCGCTGCACGCAGAGGCACAGTGCCACCGGGTATACAGATTGTTGACAAAGTGCGCTTTGTTTATGCCGGATGCGGCGTTAACGCCTTATCCGGcc                                                                              <  2:677213/100‑1 (MQ=255)                             gcTGCACGCAGAGGCACAGCGCCACCGGGTATACAGATTGTTGACAAAGTGCGCTTTGTTTATGCCGGATGCGGCGTTAACGCCTTATCCGGCCTACaaa                                                                        >  2:1256370/1‑100 (MQ=255)                             gcTGCACGCAGAGGCACAGCGCCACCGGGTATACAGATTGTTGACAAAGTGCGCTTTGTTTATGCCGGATGCGGCGTTAACGCCTTATCCGGCCTACaaa                                                                        >  2:1068255/1‑100 (MQ=255)                              cTGCACGCAGAGGCACAGCGCCACCGGGTATACAGATTGTTGACAAAGTGCGCTTTGTTTATGCCGGATAAGGCGTTAACGCCGCATCCGGcataaacaa                                                                       <  1:1763549/100‑9 (MQ=25)                              cTGCACGCAGAGGCACAGCGCCACCGGGTATACAGATTGTTGACAAAGTGCGCCTGGTTTATGCCGGATGCGGCGTTAACGCCTTATCCGGCCTACaaa                                                                        >  1:58529/1‑99 (MQ=255)                                 cACGCAGAGGCACAGCGCCACCGGGTATACAGATTGTTGACAAAGTGCGCTTTGTTTATGCCGGATGCGGCGTTAACGCCTTATCCGGCCTACAAAATCg                                                                    <  2:1149932/100‑1 (MQ=255)                                 cACGCAGAGGCACAGCGCCACCGGGTATACAGATTGTTGACAAAGTGCGCTTTGTTTATGCCGGATGCGGCGTTAACGCCTTATCCGGCCTACAAAATCg                                                                    <  1:1100084/100‑1 (MQ=255)                                 cACGCAGAGGCACAGCGCCACCGGGTATACAGATTGTTGACAAAGTGCGCTTTGTTTATGCCGAATCCGGCGTTAACACCTTATCCGGCCTACAAAATCg                                                                    >  2:592340/1‑100 (MQ=255)                                   cGCAGAGGCACAGCGCCACCGGGTATACAGATTGTTGACAAAGTGCGCTTTGTTTATGCCGGATGCGGCGTTAACGCCTTATCCGGCCTACAAAATCGTg                                                                  >  2:991970/1‑100 (MQ=255)                                   cGCAGAGGCACAGCGCCACCGGGTATACAGATTGTTGACAAAGTGCGCTTTGTTTATGCCGGATGCGGCGTTAACGCCTTATCCGGCCTACAAAATCGTg                                                                  <  2:734117/100‑1 (MQ=255)                                    gCAGAGGCACAGCGCCACCGGGTATACAGATTGTTGACAAAGTGCGCTTTGTTTATGCCGGATGCGGCGTTAACGCCTTATCCGGCCTACAAAATCGTGc                                                                 <  2:706941/100‑1 (MQ=255)                                        agGCACAGCGCCACCGGGTATACAGATTGTTGACAAAGTGCGCTTTGTTTATGCCGGATGCGGCGTTAACGCCTTATCCGGCCTACAAAATCGTGCAAAt                                                             <  1:1606729/100‑1 (MQ=255)                                         ggCACAGCGCCACCGGGTATACAGATTGTTGACAAAGTGCGCTTTGTTTATGCCGGATGCGGCGTTAACGCCTTATCCGGCCTACAAAATCGTGCAAAtt                                                            >  2:1523415/1‑100 (MQ=255)                                         ggCACAGCGCCACCGGGTATACAGATTGTTGACAAAGTGCGCTTTGTTTATGCCGGATGCGGCGTTAACGCCTTATCCGGCCTACAAAATCGTGCAAAtt                                                            >  2:1532130/1‑100 (MQ=255)                                           cacaGCGCCACCGGGTATACAGATTGTTGACAAAGTGCGCTTTGTTTATGCCGGATGCGGCGTTAACGCCTTATCCGGCCTACAAAATCGTGCAAATTCa                                                          >  2:1680175/1‑100 (MQ=255)                                           cacaGCGCCACCGGGTATACAGATTGTTGACAAAGTGCGCTTTGTTTATGCCGGATGCGGCGTTAACGCCTTATCCGGCCTACAAAATCGTGCAAATTCa                                                          >  2:1309767/1‑100 (MQ=255)                                           cacaGCGCCACCGGGTATACAGATTGTTGACAAAGTGCGCTTTGTTTATGCCGGATGCGGCGTTAACGCCTTATCCGGCCTACAAAATCGTGCAAATTCa                                                          >  2:488539/1‑100 (MQ=255)                                            acaGCGCCACCGGGTATACAGATTGTTGACAAAGTGCGCTTTGTTTATGCCGGATGCGGCGTTAACGCCTTATCCGGCCTACAAAATCGTGCAAATTCa                                                          <  2:1476944/99‑1 (MQ=255)                                             caGCGCCACCGGGTATACAGATTGTTGACAAAGTGCGCTTTGTTTATGCCGGATGCGGCGTTAACGCCTTATCCGGCCTACAAAATCGTGCAAATTCAAt                                                        >  2:67366/1‑100 (MQ=255)                                             caGCGCCACCGGGTATACAGATTGTTGACAAAGTGCGCTTTGTTTATGCCGGATGCGGCGTTAACGCCTTATCCGGCCTACAAAATCGTGCAAATTCAAt                                                        <  2:338379/100‑1 (MQ=255)                                               gcgcCACCGGGTATACAGATTGTTGACAAAGTGCGCTTTGTTTATGCCGGATGCGGCGTTAACGCCTTATCCGGCCTACAAAATCGTGCAAATTCAATAc                                                      <  1:1032586/100‑1 (MQ=255)                                               gcgcCACCGGGTATACAGATTGTTGACAAAGTGCGCTTTGTTTATGCCGGATGCGGCGTTAACGCCTTATCCGGCCTACAAAATCGTGCAAATTCAATAc                                                      <  1:1033347/100‑1 (MQ=255)                                                cgcCACCGGGTATACAGATTGTTGACAAAGTGCGCTTTGTTTATGCCGGATGCGGCGTTAACGCCTTATCCGGCCTACAAAATCGTGCAAATTCAATAc                                                      <  2:1047255/99‑1 (MQ=255)                                                cgcCACCGGGTATACAGATTGTTGACAAAGTGCGCTTTGTTTATGCCGGATGCGGCGTTAACGCCTTATCCGGCCTACAAAATCGTGCAAATTCAATACa                                                     >  2:911581/1‑100 (MQ=255)                                                    aCCGGGTATACAGATTGTTGACAAAGTGCGCTTTGTTTATGCCGGATGCGGCGTTAACGCCTTATCCGGCCTACAAAATCGTGCAAATTCAATACAtt                                                   <  2:1293695/98‑1 (MQ=255)                                                    aCCGGGGATAGAGATTGTTGACAAAGTGCGCTTTGTTTATGCCGGATGCGGCGTTAACGCCTTATCCGGCCTACAAAATCGTGCAAATTCAATACATTGc                                                 <  1:143326/100‑1 (MQ=255)                                                      cGGGTATACAGATTGTTGACAAAGTGCGCTTTGTTTATGCCGGATGCGGCGTTAACGCCTTATCCGGCCTACAAAATCGTGCAAATTCAATACATTGC‑‑A‑g                                            <  2:541449/100‑1 (MQ=255)                                                      cGGGTATACAGATTGTTGACAAAGTGCGCTTTGTTTATGCCGGATGCGGCGTTAACGCCTTATCCGGCCTACAAAATCGTGCAAATTCAATACATTGC‑‑A‑g                                            >  1:1762033/1‑100 (MQ=255)                                                       gggTATACAGATTGTTGACAAAGTGCGCTTTGTTTATGCCGGATGCGGCGTTAACGCCTTATCCGGCCTACAAAATCGTGCAAATTCAATACATTGC‑‑A‑Ga                                           >  1:503843/1‑100 (MQ=255)                                                        ggTATACAGATTGTTGACAAAGTGCGCTTTGTTTATGCCGGATGCGGCGTTAACGCCTTATCCGGCCTACAAAATCGTGCAAATTCAATACATTGC‑‑A‑Gaa                                          >  2:1191166/1‑100 (MQ=255)                                                         gTATACAGATTGTTGACAAAGTGCGCTTTGTTTATGCCGGATGCGGCGTTAACGCCTTATCCGGCCTACAAAATCGTGCAAATTCAATACATTGC‑‑A‑GAAt                                         <  2:170582/100‑1 (MQ=255)                                                          tataCAGATTGTTGACAAAGTGCGCTTTGTTTATGCCGGATGCGGCGTTAACGCCTTATCGGGCCTACAAAATCGTGCAAATTCAATACATTGC‑‑A‑GAAtt                                        <  1:620334/100‑1 (MQ=255)                                                          tataCAGATTGTTGACAAAGTGCGCTTTGTTTATGCCGGATGCGGCGTTAACGCCTTATCCGGCCTACAAAATCGTGCAAATTCAATACATTGC‑‑A‑GAAtt                                        <  1:595451/100‑1 (MQ=255)                                                            taCAGATTGTTGACAAAGTGCGCTTTGTTTATGCCGGATGCGGCGTTAACGCCTTATCCGGCCTACAAAATCGTGCAAATTCAATACATTGC‑‑A‑GAAtt‑‑tt                                    >  1:234890/1‑100 (MQ=255)                                                              cAGATTGTTGACAAAGTGCGCTTTGTTTATGCCGGATGCGGCGTTAACGCCTTATCCGGCCTACAAAATCGTGCAAATTCAATACATTGC‑‑A‑GAATT‑‑TTCg                                  <  1:1744610/100‑1 (MQ=255)                                                                gATTGTTGACAAAGTGCGCTTTGTTTATGCCGGATGCGGCGTTAACGCCTTATCCGGCCTACAAAATCGTGCAAATTCAATACATTGC‑‑A‑GAATT‑‑TTCGTa                                >  1:1403868/1‑100 (MQ=255)                                                                   tgttgACAAAGTGCGCTTTGTTTATGCCGGATGCGGCGTTAACGCCTTATCCGGCCTACAAAATCGTGCAAATTCAATACATTGC‑‑A‑GAATT‑‑TTCGTAGGc                             <  2:1740260/100‑1 (MQ=255)                                                                    gttgACAAAGTGCGCTTTGTTTATGCCGGATGCGGCGTTAACGCCTTATCCGGCCTACAAAATCGTGCAAATTCAATACATTGC‑‑A‑GAATT‑‑TTCGTAGGcc                            >  1:103363/1‑100 (MQ=255)                                                                    gttgACAAAGTGCGCTTTGTTTATGCCGGATGCGGCGTTAACGCCTTATCCGGCCTACAAAATCGTGCAAATTCAATACATTGC‑‑A‑GAATT‑‑TTCGTAGGcc                            >  1:1169348/1‑100 (MQ=255)                                                                       gACAAAGTGCGCTTTGTTTATGCCGGATGCGGCGTTAACGCCTTATCCGGCCTACAAAATCGTGCAAATTCAATACATTGC‑‑A‑GAATT‑‑TTCGTAGGCCTGa                         <  2:1702200/100‑1 (MQ=255)                                                                       gACAAAGTGCGCTTTGTTTATGCCGGATGCGGCGTTAACGCCTTATCCGGCCTACAAAATCGTGCAAATTCAATACATTGC‑‑A‑GAATT‑‑TTCGTAGGCCTGa                         >  2:1340399/1‑100 (MQ=255)                                                                       gACAAAGTGCGCTTTGTTTATGCCGGATGCGGCGTTAACGCCTTATCCGGCCTACAAAATCGTGCAAATTCAATACATTGC‑‑A‑GAATT‑‑TTCGTAGGCCTGa                         >  1:1683413/1‑100 (MQ=255)                                                                         cAAAGTGCGCTTTGTTTATGCCGGATGCGGCGTTAACGCCTTATCCGGCCTACAAAATCGTGCAAATTCAATACATTGC‑‑A‑GAATT‑‑TTCGTAGGCCTGATa                       >  2:901313/1‑100 (MQ=255)                                                                          aaaGTGCTCTTTGTTTATGCCGGATGCGGCGTTAAAGGCTTATCCGGCCTACAAAATCGTGCAAATTCAATACATTGC‑‑A‑GAATT‑‑TTCGTAGGCCTGATa                       >  1:1055188/1‑99 (MQ=255)                                                                          aaaGTGCGCTTTGTTTATGCCGGATGCGGCGTTAACGCCTTATCCGGCCTACAAAATCGTGCAAATTCAATACATTGC‑‑A‑GAATT‑‑TTCGTAGGCCTGATaa                      >  1:1055189/1‑100 (MQ=255)                                                                            aGTGCGCTTTGTTTATGCCGGATGCGGCGTTAACGCCTTATCCGGCCTACAAAATCGTGCAAATTCAATACATTGC‑‑A‑GAATT‑‑TTCGTAGGCCTGATAAGc                    >  1:12707/1‑100 (MQ=255)                                                                            aGTGCGCTTTGTTTATGCCGGATGCGGCGTTAACGCCTTATCCGGCCTACAAAATCGTGCAAATTCAATACATTGC‑‑A‑GAATT‑‑TTCGTAGGCCTGATAAGc                    <  1:520659/100‑1 (MQ=255)                                                                                cgcTTTGTTTATGCCGGATGCGGCGTTAACGCCTTATCCGGCCTACAAAATCGTGCAAATTCAATACATTGC‑‑A‑GAATT‑‑TTCGTAGGCCTGATAAGCGGAg                >  1:1010374/1‑100 (MQ=255)                                                                                cgcTTTGTTTATGCCGGATGCGGCGTTAACGCCTTATCCGGCCTACAAAATCGTGCAAATTCAATACATTGC‑‑A‑GAATT‑‑TTCGTAGGCCTGATAAGCGGAg                <  1:679691/100‑1 (MQ=255)                                                                                cgcTTTGTTTATGCCGGATGCGGCGTTAACGCCTTATCCGGCCTACAAAATCGTGCAAATTCAATACATTGC‑‑A‑GAATT‑‑TTCGTAGGCCTGATAAGCGGAg                <  1:1536712/100‑1 (MQ=255)                                                                                 gcTTTGTTTATGCCGGATGCGGCGTTAACGCCTTATCCGGCCTACAAAATCGTGCAAATTCAATACATTGC‑‑A‑GAATT‑‑TTCGTAGGCCTGATAAGCGGAgc               <  1:363970/100‑1 (MQ=255)                                                                                  tgatttgtaggCCGGATAAGGCGTTCACGCCGCATCCGGCCTAC‑AAATCATGCAAATCCAAAATATTGC‑‑A‑GGAGT‑‑TATGTAGGCCTGATAAGCGTAgcgc             >  1:551126/11‑100 (MQ=18)                                                                                  cTTTGTTTATGCCGGATGCGGCGTTAACGCCTTATCCGGCCTACAAAATCGTGCAAATTCAATACATTGC‑‑A‑GAATT‑‑TTCGTAGGCCTGATAAGCGGAgcg              >  2:1438777/1‑100 (MQ=255)                                                                                  cTTTGTTTATGCCGGATGCGGCGTTAACGCCTTATCCGGCCTACAAAATCGTGCAAATTCAATACATTGC‑‑A‑GAATT‑‑TTCGTAGGCCTGATAAGCGGAgcg              >  2:1548186/1‑100 (MQ=255)                                                                                    ttGTTTATGCCGGATGCGGCGTTAACGCCTTATCCGGCCTACAAAATCGTGCAAATTCAATACATTGC‑‑A‑GAATT‑‑TTCGTAGGCCTGATAAGCGGAGCGCa            <  2:1631715/100‑1 (MQ=255)                                                                                    ttGTTTATGCCGGATGCGGCGTTAACGCCTTATCCGGCCTACAAAATCGTGCAAATTCAATACATTGC‑‑A‑GAATT‑‑TTCGTAGGCCTGATAAGCGGAGCGCa            >  1:10159/1‑100 (MQ=255)                                                                                    ttGTTTATGCCGGATGCGGCGTTAACGCCTTATCCGGCCTACAAAATCGTGCAAATTCAATACATTGC‑‑A‑GAATT‑‑TTCGTAGGCCTGATAAGCGGAGCGCa            <  2:1037093/100‑1 (MQ=255)                                                                                      gTTTATGCCGGATGCGGCGTTAACGCCTTATCCGGCTTACAAAATCGTGCAAATTCAATACATTGC‑‑A‑GAATT‑‑TTCGTAGGCCTGATAAGCGGAGCGCATc          <  2:1502798/100‑1 (MQ=255)                                                                                       tttATGCCGGATGCGGCGTTAACGCCTTATCCGGCCTACAAAATCGTGCAAATTCAATACATTGC‑‑A‑GAATT‑‑TTCGTAGGCCTGATAAGCGGAGCGCATCa         <  1:874024/100‑1 (MQ=255)                                                                                       tttATGCCGGATGCGGCGTTAACGCCTTATCCGGCCTACAAAATCGTGCAAATTCAATACATTGC‑‑A‑GAATT‑‑TTCGTAGGCCTGATAAGCGGAGCGCATCa         <  2:1144574/100‑1 (MQ=255)                                                                                       tttATGCCGGATGCGGCGTTAACGCCTTATCCGGCCTACAAAATCGTGCAAATTCAATACATTGC‑‑A‑GAATT‑‑TTCGTAGGCCTGATAAGCGGAGCGCATCa         <  1:59568/100‑1 (MQ=255)                                                                                       tttATGCCGGATGCGGCGTTAACGCCTTATCCGGCCTACAAAATCGTGCAAATTCAATACATTGC‑‑A‑GAATT‑‑TTCGTAGGCCTGATAAGCGGAGCGCATCa         <  1:59539/100‑1 (MQ=255)                                                                                        gtaggCCGGATAAGGCGTTCACGCCGCATCCGGCCTACAAAAACATGCAAATTCAATAGATTGC‑‑A‑GAGAT‑‑TATGTAGGCCTGATAAGCGTAGCGCATCAg        >  1:311414/5‑100 (MQ=17)                                                                                         taggCCGGATAAGGCGTTAACGCCGCATCCGGCCTACAAAATCGTGCAAATTCAATACATTGC‑‑A‑GAATT‑‑TTCGTAGGCCTGATAAGCGGAGCGCATCAgg       >  1:424834/4‑100 (MQ=18)                                                                                          aggCCGGATAAGGCGTTTACGCCGCATCCGGCCTACAAAAGCGTGCAAATTCAATACATTGC‑‑ATGGGCC‑‑AT‑GTAGGCCTGATAAGCGTAGCGCATCAGGc      >  1:576835/3‑100 (MQ=17)                                                                                          aggCCGGATAAGGCGTTTACGCCGCATCCGGCCTACAAAAGCATGCAAATTCAATATATTGC‑‑A‑GAGTT‑‑CACGTAGGCCTGATAAGCGTAGCGCATCAGGc      >  1:1304631/3‑100 (MQ=17)                                                                                          aggCCGGATAAGGCGTTTACGCCGCATCCGGCCTACAAAAGCATGCAAATTCAATATATTGC‑‑A‑GAGTT‑‑CACGTAGGCCTGATAAGCGTAGCGCATCAAGc      >  2:1127876/3‑100 (MQ=14)                                                                                          aggCCGGATAAGGCGTTTACGCCGCATCCGGCCTACAAAACAGTGCAAATTCAATATATTGC‑‑A‑ATCTT‑‑CTGGTAGGCCTGATAAGCATCGCGCATCAGGc      >  2:1387766/3‑100 (MQ=18)                                                                                          aggCCGGATAAGGCGTTCACGCCGCATCCGGCCTACAAAATTGTGTAAATTCAATAAATTGA‑‑T‑CGATC‑‑CACGTAGGCCTGATAAGCGTAGCGCATCAgg       >  2:80290/3‑99 (MQ=18)                                                                                          aggCCGGATAAGGCGTTCACGCCGCATCCGGCCTACAAAAACCTGCAAATTCAATATATTGC‑‑A‑GGATA‑‑CGTGTAGGCCTGATAAGCGTAGCGCATCAGGc      >  2:863322/3‑100 (MQ=17)                                                                                          aggCCGGATAAGGCGTTCACGCCGCATCCGGCCTAC‑AAATCATGCAAATCCAAAATATTGC‑‑A‑GGAGT‑‑TATGTAGGCCTGATAAGCGTAGCGCATCAGGCa     >  1:1016666/3‑100 (MQ=21)                                                                                          aggCCGGATAAGGCGTTCACGCCGCATCCGGCCTAC‑AAATCATGCAAATCCAAAATATTGC‑‑A‑GGAGT‑‑TATGTAGGCCTGATAAGCGTAGCGCATCAGGCa     >  1:1634581/3‑100 (MQ=21)                                                                                          aggCCGGATAAGGCGTTCACGCCGCATCCGGCCTAC‑AAATCATGCAAATCCAAAATATTGC‑‑A‑GGAGT‑‑TATGTAGGCCTGATAAGCGTAGCGCATCAGGCa     >  1:371322/3‑100 (MQ=21)                                                                                           tGCCGGATGCGGCGTTAACGCCTTATCCGGCCTACAAAATCGTGCAAATTCAATACATTGC‑‑A‑GAATT‑‑TTCGTAGGCCTGATAAGCGGAGCGCATCAGGCa     >  2:353382/1‑100 (MQ=255)                                                                                           tGCCGGATGCGGCGTTAACGCCTTATCCGGCCTACAAAATCGTGCAAATTCAATACATTGC‑‑A‑GAATT‑‑TTCGTAGGCCTGATAAGCGGAGCGCATCAGGCa     >  1:1423901/1‑100 (MQ=255)                                                                                           ggCCGGATAAGGCGTTTACGCCGCATCCGGCCTACAAAAATGTGCAAATTCAATAAATTGCA‑A‑TTCAA‑‑CTTGTAGGCCTGATAAGCGCAGCGCATCAGGc      >  2:1356032/2‑100 (MQ=18)                                                                                           ggCCGGATAAGGCGTTTACGCCGCATCCGGCCTACAAAAATGTGCAAATTCAATAAATTGC‑‑A‑‑ATTCAACTAATAGGCCTGATAAGCGCAGCGCATCAGGc      >  2:1382502/2‑100 (MQ=18)                                                                                           ggCCGGATAAGGCGTTCACGCCGCATCCGGCCTACAAAATCATGAAAATTCAATACATTGC‑‑A‑AGATT‑‑TTCGTAGGCCTGATAAGAGTAGCGCATCAGGCa     >  2:1792775/2‑100 (MQ=17)                                                                                           ggCCGGATAAGGCGTTCACGCCGCATCCGGCCTACAAAACTTTGCAAATTCAATATATTGC‑‑A‑ATCTC‑‑CGTGTAGGCCTGATAAGCGTAGCGCATCAGGCa     >  2:972699/2‑100 (MQ=18)                                                                                            gCCGGATAAGGCGTTTACGCCGCATCCGGCCTACAAAATCGTGCTAATTCAATATATTGC‑‑A‑GAAAC‑‑CTTGTAGGACTGATAAGAGTAGCGCATAAggaa     >  2:785541/1‑97 (MQ=17)                                                                                            gCCGGATAAGGCGTTCACGCCGCATCCGGCCTACGAAATCGTGCTAATTCAATATATTGT‑‑G‑GGCTC‑‑CACGTAGGCCTGATAAGCGTAGCGCATCAGGCag    >  2:1642074/1‑99 (MQ=18)                                                                                            gCCGGATAAGGCGTTCACGCCGCATCCGGCCTACAAAGTCGTATAAATTCAATGTATTGC‑‑A‑GTTTA‑‑TTGGTAGGCCTGATAAGCGCAGCGCATCAGGCaa    >  1:466815/1‑99 (MQ=21)                                                                                            gCCGGATAAGGCGTTCACGCCGCATCCGGCCTAC‑AAACCGCGCAAATTCAATATATTGCGGA‑GAAAA‑‑T‑‑GTAGGCCTGATAAGCGTAGCGCATCAGGctgt   >  2:1815709/1‑97 (MQ=21)                                                                                              cGGATGCGGCGTTAACGCCTTATCCGGCCTACAAAATCGTGCAAATTCAATACATTGC‑‑A‑GAATT‑‑TTCGTAGGCCTGATAAGCGGAGCGCATCAGGCAttt  >  2:418232/1‑100 (MQ=255)                                                                                              cGGATGCGGCGTTAACGCCTTATCCGGCCTACAAAATCGTGCAAATTCAATACATTGC‑‑A‑GAATT‑‑TTCGTAGGCCTGATAAGCGGAGCGCATCAGGCAtt   >  2:729353/1‑99 (MQ=255)                                                                                                    |                                                                                                   AGGGAATGAAGTATTGCTGCCGGATGGCGCTGCACGCAGAGGCACAGCGCCACCGGGTATACAGATTGTTGACAAAGTGCGCTTTGTTTATGCCGGATGCGGCGTTAACGCCTTATCCGGCCTACAAAATCGTGCAAATTCAATACATTGC‑‑A‑GAATT‑‑TTCGTAGGCCTGATAAGCGGAGCGCATCAGGCATTT  >  NC_000913/3705872‑3706064

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