Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	313,543	0	G	A	58.3%	34.0 / 74.1	60	intergenic (‑301/‑1748)	ykgR/insE1	putative membrane protein YkgR/IS3 element protein InsE
Reads supporting (aligned to +/- strand):  ref base G (10/15);  new base A (19/16);  total (29/31)
Fisher's exact test for biased strand distribution p-value = 3.06e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 8.72e-01
Rejected as consensus: Frequency below/above cutoff threshold.

CGAGTTTGGTCAACCAATTTTTTGACATGTATCACAAATTTGAATAGATGTATTACATCAACTATCTTTTATTGCACCAACGTCATTGATATATGTCGCCTGAAGTCAGTTCCGGGAATGAGTCTGATCTCAAGACTGGCCCAGTCCGGGCGTTGATTGGTGCTGAGGAGCATATCGCATCTCATCATAATGTCGTATCTCCTGGGGTGTTATACAAGATATCGTTGTTGGTGACCTGGGAGAGGAATTGAGTTCTATTAAACCGTCAACTATGCCGGATACATACTGGATTACACTGCAGGCACGCCTTATGAGAGAACGTGCCGCAGTGACGGGTTAATTATCTGAAAGAATTTGTGAGGCTGTATCGGTTACTCATTGATTTGATAGTTTTACTCTCGGGAGAATAATAGATATTTAATCCATTAACGGAAACCAGCCAGTTCCTTTCGATGCCTGAATTTGATCCCATAGTTTATCCGGAATGGTT  >  NC_000913/313299‑313788                                                                                                                                                                                                                                                     |                                                                                                                                                                                                                                                      cGAGTTTGGTCAACCAATTTTTTGACATGTATCACAAATTTGAATAGATGTATTACATCAACTATCTTTTATTGCACCAACGTCATTGATATATGTCGCCTGAAGTCAGTTCCGGGAATGAGTCTGATCTCAAGACTGGCCCAGTCCGGGCGTTGATTGGTGCTGAGGAGCATATCGCATCTCATCATAATGTCGTATCTCCTGGGGTGTTATACAAGATATAGTTGTTGGTGACCTGGGAGAGAAATTg                                                                                                                                                                                                                                                  >  1:290516/1‑250 (MQ=255)   aGTTTGGTCAACCAATTTTTTGACATGTATCACAAATTTGAATAGATGTATTACATCAACTATCTTTTATTGCACCAACGTCATTGATATATGTCGCCTGAAGTCAGTTCCGGGAATGAGTCTGATCTCAAGACTGGCCCAGTCCGGGCGTTGATTGGTGCTGAGGAGCATATCGCATCTCATCATAATGTCGTATCTCCTGGGGTGTTATACAAGATATCGTTGTTGGTGACCTGGGAGAGGAATTGAGt                                                                                                                                                                                                                                               >  1:121575/1‑251 (MQ=255)                      ttGACATGTATCACAAATTTGAATAGATGTATTACATCAACTATCTTTTATTGCACCAACGTCATTGATATATGTCGCCTGAAGTCAGTTCCGGGAATGAGTCTGATCTCAAGACTGGCCCAGTCCGGGCGTTGATTGGTGCTGAGGAGCATATCGCATCTCATCATAATGTCGTATCTCCTGGGGTGTTATACAAGATATCGTTGTTGGTGACCTGGGAGAGAAATTGAGTTCTATTAAACCGTCAACt                                                                                                                                                                                                                             <  1:171856/250‑1 (MQ=255)                        gACATGTATCACAAATTTGAATAGCTGTATTACATCAACTATCTTTTATTGCACCAACGTCATTGATATATGTCGCCTGAAGTCAGTTCCGGGAATGAGTCTGATCTCAAGACTGGCCCAGTCCGGGCGTTGATTGGGGCTGAGGAGCATATCGCATCTCATCATAATGTCGTATCTCCTGGGGTGTTATACAAGATATCGTTGTTGGTGACCTGGGAGAGAAATTGAGTTCTATTAAACCGTCAACTATg                                                                                                                                                                                                                          >  1:437544/1‑251 (MQ=255)                             gTATCACAAATTTGAATAGATGTATTACATCAACTATCCTTTATTGCACCAACGTCATTGATATATGTCGCCTGAAGTCAGTTCCGGGAATGAGTCTGATCTCAAGACTGGCCCAGTCCGGGCGTTGATTGGTGCTGAGGAGCATATTGCATCTCATCATAATGTCGTATCTCCTGGGGTGTTATACAAGATATCGTTGTTGGTGACCTGGGAGAGAAATTGAGTTCTa                                                                                                                                                                                                                                           >  2:73583/1‑229 (MQ=255)                             gTATCACAAATTTGAATAGATGTATTACATCAACTATCCTTTATTGCACCAACGTCATTGATATATGTCGCCTGAAGTCAGTTCCGGGAATGAGTCTGATCTCAAGACTGGCCCAGTCCGGGCGTTGATTGGTGCTGAGGAGCATATTGCATCTCATCATAATGTCGTATCTCCTGGGGTGTTATACAAGATATCGTTGTTGGTGACCTGGGAGAGAAATTGAGTTCTa                                                                                                                                                                                                                                           <  1:73583/229‑1 (MQ=255)                                 cacaAATTTGAATAGATGTATTACATCAACTATCTTTTATTGCACCAACGTCATTGATATATGTCGCCTGAAGTCAGTTCCGGGAATGAGTCTGATCTCAAGACTGGCCCAGTCCGGGCGTTGATTGGTGCTGAGGAGCATATCGCATCTCATCATAATGTCGTATCTCCTGGGGTGTTATACAAGATATCGTTGTTGGTGACCTGGGAGAGAAATTGAGTTCTATTAAACCGTCAACTATGCCGGataca                                                                                                                                                                                                                 <  1:457611/251‑1 (MQ=255)                                          gAATAGATTTATTACATCAACTATCTTTTATTGCACCAACGTCATTGATATATGTCGCCTGAAGTCAGTTCCGGGAATGAGTCTGATCTCAAGACTGGCCCAGTCCGGGCGTTGATTGGTGCTGAGGAGCATATCGCATCTCATCATAATGTCGTATCTCCTGGGGTGTTATACAAGATATCGTTGTTGGTGACCTGGGAGAGGAATTGAGTTCTATTAAACCGGCAACTATGCCGGATACATACTGGAtt                                                                                                                                                                                                        <  2:10226/251‑1 (MQ=255)                                           aaTAGATGTATTACATCAACTATCTTTTATTGCACCAACGTCATTGATATATGTCGCCTGACGTCAGTTCCGGGACTGAGTCTGATCTCAAGACTGGCCCAGTCCGGGCGTTGATTGGTGCTGAGGAGCATATCGCATCTCATCATAATGTCGTATCTCCTGGGGTGTTATACAAGATATCGTTGTTGTTGACCTGGGAGAGAAATTGAGTTCTATTAAACCGTCAACTATGCCGGCTACATACTGGAtt                                                                                                                                                                                                        >  2:330698/1‑250 (MQ=255)                                            aTAGATGTATTACATCAACTATCTTTTATTGCACCAACGTCATTGATATATGTCGCCTGAAGTCAGTTCCGGGAATGAGTCTGATCTCAAGACTGGCCCAGTCCGGGCGTTGATTGGTGCTGAGGAGCATATCGCATCTCATCATAATGTCGTATCGCCTGGGGTGTTATACAAGATATCGTTGTTGGTGACCTGGGAGAGGAATTGAGTTCTATTAAACCGTCAACTATGCCGGATACATACTAGATTac                                                                                                                                                                                                      <  1:117873/251‑1 (MQ=255)                                                  gTATTACATCAACTATCTTGTATTGCACCAACGTCATTGATATATGTCGCCTGAAGTCAGTTCCGGGAATGAGTCTGATCTCAAGACTGGCCCAGTCCGGGCGTTGATTGGTGCTGAGGAGCATATCGCATCTCATCATAATGTCGTATCTCCTTGGGTGTTATACAAGATATCGTTGTTGGTGACCTGGGAGAGAAATTGAGTTCTATTAATCCGTCAACTATGCCGGATACATACTGGATGACACTGCa                                                                                                                                                                                                <  1:43205/251‑1 (MQ=255)                                                                    tttATTGCACGAAAGACATTGATATATGTCGCCTGAAGTCAGTTTCAGGAATGAGTCTAATCTCAAGACAGGCCCAGTCCGAGCGTTGATTGGTGCTGAGGAGCATATCGCATCTCATCATAATGTCGTATATCCTGGGGTGTTATACAAGATATCGATGTTGGTGACCTGGGAGAGAAATTGAGTTCTATTAAAACGTCAACTATGCCGGCTACATACTGGATTACACTGAATGCACgacatatgtg                                                                                                                                                                                 >  2:416185/1‑239 (MQ=255)                                                                    tttATTGCACCAACGTCATTGATATATGTCGCCTGAAGTCAGTTCCGGGAATGAGTCTGATCTCAAGACAGGCCCAGTCCGGGCGTTGATTGGTGCTGAGGAGCATCTCGCATCTCATCATAATGTCGTCTCTCATGGGGTGTTATACAAGATATCGTTGTTGGTGACCTGGGAGAGAAATTGAGTTCTATTAAACCGTCAACTATGCCGGATACATACTGGATTACACTGCAGACACGCCTTATgag                                                                                                                                                                                 >  1:440815/1‑248 (MQ=255)                                                                         tGCACCAACGTCATTGATATATGTCGCCTGAAGTCAGTTCCGGGAATGAGTCTGATCTCAAGACTGGCCCAGTCCGGGCGTTGATTGGTGCTGAGGAGCATATCGCATCTCATCATAATGTCGTATCTCCTGGGGTGTTATACAAGATATCGTTGTTGGTGACCTGGGAGAGAAATTg                                                                                                                                                                                                                                                  >  2:46788/1‑178 (MQ=255)                                                                         tGCACCAACGTCATTGATATATGTCGCCTGAAGTCAGTTCCGGGAATGAGTCTGATCTCAAGACTGGCCCAGTCCGGGCGTTGATTGGTGCTGAGGAGCATATCGCATCTCATCATAATGTCGTATCTCCTGGGGTGTTATACAAGATATCGTTGTTGGTGACCTGGGAGAGAAATTg                                                                                                                                                                                                                                                  <  1:46788/178‑1 (MQ=255)                                                                          gcgccatcGTCATGGATATATGTCACCTGAAGTCAGTTCCGGGATTGAGTCTGATCTCAAGACTGGCCCAGTCCGGGCGTTGATTTGTGCTGAGGATCATATCGCATCTCATCATAATGTCGTATCTCCTGGTGTTTTATCCAAGATATCTTTATTGGGGACCTGGGAGAGGAATTGTTTTCTATGAGACCGTCAACTATGCCGCATACATACTGGAGTACACTGCAGCCACGCCttgggagagggcgtt                                                                                                                                                                         <  2:262300/243‑14 (MQ=255)                                                                                aCGTCATTGATATATGTCGCCTGAAGTCAGTTCCGGGAATGAGTCTGATCTCAAGACTGGCCCAGTCCGGGCGTTGATTGGTGCTGAGGAGCATATCGCATCTCATCATAATGTCGTATCTCCTGGGGTGTTATACAAGATATCGTTGTTGGTGACCTGGGAGAGGAATTGAGTTCTATTAAACCGTCAACTATGCCGGATACATACTGGATTACACTGCAGGCACGCCTTATGAGAGAACGTGCCGCAGt                                                                                                                                                                  <  2:401532/251‑1 (MQ=255)                                                                                     aTTGATATATGTCGCCTGAAGTCAGTTCCGGGAATGAGTCTGATCTCCAGACTAGCCCAGTCCGGGCGTTGATTGGTGCTGAGGAGCATATCGCATCTCATCATAATGTCGTATCTCCTGGGGTGTTATACAAGATATCGTTGTTGGTGACCTGGGAGAGGAATTGAGTTCTATTAAACCGTCAACTATGCCGGATACATACTGGATTACACTGCAGGCACGCCTTATGAGAGAACGTGCCGCAGTGACgg                                                                                                                                                             >  1:260429/1‑251 (MQ=255)                                                                                     aTTGATATATGTCGCCTGAAGTCAGTTCCGGGAATGAGTATGATCTCACGACTGGCCCAGTCCGGGCGTTGATTGGTGCTGAGGAGCATATCGCATCTCATCATAATGTCGTATCTCCTGGGGTGTTATACAAGATCTCGTTGTTGGTGACCTGGGAGAGAAATTGAGTTCTATTAAACCGTCAACTATGCCGGATACATACTGGATTACACTGCAGGCACGCCTTATGAGAGAACGTGCCGCAGTGACgg                                                                                                                                                             >  1:446203/1‑251 (MQ=255)                                                                                           atatGTCGCCTGAAGTCAGTTCCGGGAATGAGTCTGATCTCAAGACTGGCCCAGTCCGGGCGTTGATTGGTGCTGAGGAGCATATCGCATCTCATCATAATGTCGTATCTCCTGGGGTGTTATACAAGATATCGTTGTTGGTGACCTGGGAGAGGAATTGAGTTCTATTAAACCGTCAACTATGCCGGATACATACTGGATTACACTGCAGGCACGCCTTATGAGAGAACGTGCCGCAGTGACGGGTTaa                                                                                                                                                        <  2:15968/250‑1 (MQ=255)                                                                                           atatGTCGCCTGAAGTCAGTTCCGGGAATGAGTCTGATCTCAAGACTGGCCCAGTCCGGGCGTTGATTGGTGCTGAGGAGCATATCGCATCTCATCATAATGTCGTATCTCCTGGGGTGTTATACAAGATATCGTTGTTGGTGACCTGGGAGAGGAATTGAGTTCTATTAAACCGTCAACTATGCCGGATACATACTGGATTACACTGCAGGCACGCCTTATGAGAGAACGTGCCGCAGTGACGGGTTAAt                                                                                                                                                       >  2:222955/1‑251 (MQ=255)                                                                                           atatGTCGCCTGAAGTCAGTTCCGGGAATGAGTCTGATCTCAAGACTGGCCCAGTCCGGGCGTTGATTGGTGCTGAGGAGCATATCGCATCTCATCATAATGTCGTATCTCCTGGGGTGTTATACAAGATATCGTTGTTGGTGACCTGGGAGAGGAATTGAGTTCTATTAAACCGTCAACTATGCCGGATACATACTGGATTACACTGCAGGCACGCCTTATGAGAGAACGTGCCGCAGTGACGGGTTAAt                                                                                                                                                       >  1:153882/1‑251 (MQ=255)                                                                                            tatGTCGCCTGAAGTCAGTTCCGGGAATGAGTCTGATCTCAAGACTGGCCCAGTCCGGGCGTTGATTGGTGCTGAGGAGCATATCGCATCTCATCATAATGTCGTATCTCCTGGGGTGTTATACAAGATATCGTTGTTGGTGACCTGGGAGAGAAATTGAGTTCTATTAAACCGTCAACTATGCCGGATACATACTGGATTACACTGCAGGCACGCCTTATGAGAGAACGTGCCGCAGTGACGGGTTAAtt                                                                                                                                                      >  1:377719/1‑251 (MQ=255)                                                                                               gTCGCCTGAAGTCAGTTCCGGGAATGAGTCTGATCTCAAGACTGGCCCAGTCCGGGCGTTGATTGGTGCTGAGGAGCATATCGCATCTCATCATAATGTCGTATCTCCTGGGGTGTTATACAAGATATCGTTGTTGGTGACCTGGGAGAGGAATTGAGTTCTATTAAACCGTCAACTATGCCGGATACATACTGGATTACACTGCAGGCACGCCTTATGAGAGAACGTGCCGCAGTGACGGGTTAATTAt                                                                                                                                                    <  1:311786/250‑1 (MQ=255)                                                                                                tagCCTGAAGTCAGTACCGGGAATGAGTCTGATCTCAAGACTGGCCCAGTCCGGTCGTTGATTGGTGCTGCGGCGCATATCGCATCTCCTCATCATGTCGTATCTCCTGGGGTGTTATACAAGATATCGTTGTTGGTGACCTGTTAGAGGAATTTAGTTCTATTAAACCGTCAACTATGCCGGATACATACTGGTTTACACTGCAGGCACACCTTATGAGAGAACGTGCCGCCGTGACGGGTTAATTATCt                                                                                                                                                  >  2:249187/3‑251 (MQ=255)                                                                                                     ggAAGTCAGTTCCGGGAATGAGTCTGATCTCAAGACTGGCCCAGTCCGGGCGTTGATGGGTGCGGAGGAGCATATCGCATCTCATCATAATGTCGTATCGCCTGGGGTGTTATACAAGATATCGTTGTTGGTGACCTGGGAGAGGAATTGAGTTCTATTAAACCGTCAACTATGCCGGATACATACTGGATTACACTGCAGGCACGCCTTATGAGAGAACGTGCCGCAGTGACGGGTTAATTATCTGaaat                                                                                                                                             <  2:76464/250‑2 (MQ=255)                                                                                                      gAAGTCAGTTCCGGGAATGAGTCTGATCTCAAGACTGGCCCAGTCCGGGCGTTGATTGGTGCTGAGGAGCATATCGCATCTCATCATAATGTCGTATCTCCTGGGGTGTTATACAAGATATCGTTGTTGGTGACCTGGGAGAGAAATTGAGTTc                                                                                                                                                                                                                                             <  1:8724/154‑1 (MQ=255)                                                                                                           cAGTTCCGGGAATGAGTCTGATCACAAGACCTGCCCAGTCCGGTCGTTGATTGGTGCTGAGGAGCATATCGCATCTCATCATAATGTCGTCTCTCCTGGGGTTTTATACAAGATATCGGTGTTGGTGACCTGGGAGAGATATGGATGTCGATTAAACCGTCAACCATGCCGGATACATACTGGATTGCACTGCAGGCACGCCGTAGGAGAGAACGTGCCGCAGTGATGGGTTAATTATCTGAAATAAttt                                                                                                                                        <  2:32536/250‑1 (MQ=255)                                                                                                                   ggAATGAGTCTGATCTCAAGACTGGCCCAGTCCGGGCGTTGATTGGTGCTGAGGAGCATATCGCATCTCATCATAATGTCGTATCTCCTGGGGTGTTATACAAGATATCGTTGTTGGTGACCTGGGAGAGGAATTGAGTTCTATTAAACCGTCAACTATGCCGGATACATACTGGATTACACTGCAGGCACGCCTTATGAGAGAACGTGCCGCAGTGACGGGTTAATTATCTGAAATAATTTGTGAGGCTg                                                                                                                               <  2:276688/251‑1 (MQ=255)                                                                                                                   ggAATGAGGCTGATCTCAAGACTGGCCCAGTCCGGGCGTTGATTTGTGCTGAGTAGCATATCGCATCTCGTCATAATTTCGTATCTCCTGGGGTGTTATACAAGTGATCGTTGTTGGTGACGTGGGAGGGAAATTGGTTTGTATGAATCCGTGAGCTTTGCCGGTTACATACTGGATTACACTGCAGGCACGCCTTATGAGAGAACGTGCCGGAGTGACGGGGTAATTATCTGAAAGAATTTGTGAGGCTg                                                                                                                               <  1:212423/251‑1 (MQ=255)                                                                                                                        gAGTCTGATCTCAAGACTGGCCCAGTCCGGGCGTTGATTGGTGCTGAGGAGCATATCGCATCTCATCATAATGTCGTATCTCCTGGGGTGTTATACAAGATATCGTTGTTGGTGACCTGGGAGAGAAATTGAGTTCTATTAAACCGTCAACTATGCCGGATACATACTGGATTACACTGCAGGCACGCCTTATGAGAGAACGTGCCGCAGTGACGGGTTAATTATCTGAAAGAATTTGTGAGGCTGTATCg                                                                                                                          >  1:29127/1‑251 (MQ=255)                                                                                                                           tCTGATCTCAAGACTGGCCCAGTCCGGGCGTTGATTGGTGCTGAGGAGCATATCGCATCTCATCATAATGTCGTATCTCCTGGGGTGTTATACAAGATATCGTTGTTGGTGACCTGGGAGAGAAATTGAGTTCTATTAAACCGTCAACTATGCCGGATACATACTGGATTACACTGCAGGCACGCCTTATGAGAGAACGTGCCGCAGTGACGGGTTAATTATCTGAAAGAATTTGTGAGGCTg                                                                                                                               <  1:33091/243‑1 (MQ=255)                                                                                                                           tCTGATCTCAAGACTGGCCCAGTCCGGGCGTTGATTGGTGCTGAGGAGCATATCGCATCTCATCATAATGTCGTATCTCCTGGGGTGTTATACAAGATATCGTTGTTGGTGACCTGGGAGAGAAATTGAGTTCTATTAAACCGTCAACTATGCCGGATACATACTGGATTACACTGCAGGCACGCCTTATGAGAGAACGTGCCGCAGTGACGGGTTAATTATCTGAAAGAATTTGTGAGGCTGTATCGGtt                                                                                                                       >  2:37390/1‑251 (MQ=255)                                                                                                                           tCTGATCTCAAGACTGGCCAAGTCCGGGCGTTGATTGGTGCTGCGGAGCATATCGCATCTCATCATAATGTCGTATCTCCTGGGGTGTTATACAAGATATCGTTGTTGGTGACCTGGGAGAGAAATTGAGTTCTATTAAACCGTCAACTATGCCGGATACATACTGGATTACACTGCAGGCACGCCTTATGAGAGAACGTGCCGCAGTGACGGGTTAATTATCTGAAAGAATTTGTGAGGCTg                                                                                                                               >  2:33091/1‑243 (MQ=255)                                                                                                                                       aCTGGCCCAGTCCGGGCGTTGATTGGTGCTGAGGAGCATATCGCATCTCATCATAATGTCGTATCTCCTGGGGTGTTATACAAGATATCGTTGTTGGTGACCTGGGAGAGAAATTGAGTTCTATTAAACCGTCAACTATGCCGGATACATACTGGATTACACTGCAGGCACGCCTTATGAGAGAACGTGCCGCa                                                                                                                                                                    >  2:201728/1‑194 (MQ=255)                                                                                                                                       aCTGGCCCAGTCCGGGCGTTGATTGGTGCTGAGGAGCATATCGCATCTCATCATAATGTCGTATCTCCTGGGGTGTTATACAAGATATCGTTGTTGGTGACCTGGGAGAGAAATTGAGTTCTATTAAACCGTCAACTATGCCGGATACAGACTGGATTACACTGCAGGCACGCCTTATGAGAGAACGTGCCGCa                                                                                                                                                                    <  1:201728/194‑1 (MQ=255)                                                                                                                                           gCCCAGTCCGGGCGTTGATTGGTGCTGAGGAGAATATCGCATCTCATCATAATGTCGTATCTCCTGGGGTGTTATACAAGATATCGTTGTTGGTGACCTGGGAGAGGAATTGAGTTCTATTAAACCGTCAACTATGCCGGATACATACTGGATTACACTGCAGGCACGCCTTATGAGAGAACGTGCCGCAGTGACGGGTTAATTATCTGAAAGAATTTGTGAGGCTGTATCGGTTACTCATTGATTTGATa                                                                                                       <  2:326151/251‑1 (MQ=255)                                                                                                                                                       cGTTGATTGGTGCTGAGGAGCATATCGCATCTCATCATAATGTCGTATCTCCTGGGGTGTTATACAAGATATCGTTGTTGGTGACCTGGGAGAGGAATTGAGTTCTATTAAACCGTCAACTATGCCGGATACATACTGGATTACACTGCAGGCACGCCTTATGAGAGAACGTGCCGCAGTGACGGGTTAATTATCTGAAAGAATTTGTGAGGCTGTATCGGTTACTCATTGATTTGATAGGTTTACTCTCg                                                                                           <  2:45099/251‑1 (MQ=255)                                                                                                                                                       cGTTGATTGGTGCTGAGGAGCATATCGCATCTCATCATAATGTCGTATCTCCTGGGGTGTTATACAAGATATCGTTGTTGGTGACCTGGGAGAGAAATTGAGTTCTATTAAACCGTCAACTATGCCGGATACATACTGGATTACACTGCAGGCACGCCTTATGAGAGAATGTGCCGCAGTGACGGGTTAATTATCTGAAAGAATTTGTGAGGCTGTATCGGTTACTCATTGATTTGATAGTTTTACTCTCg                                                                                           >  1:363815/1‑251 (MQ=255)                                                                                                                                                                gTGCTGAGGAGCCTATCGCATCTCATCATAATGTCGTATCACCTGGGGTGTTATCCACGATATCGTTGTTGGCGACCTGGGAGAGGAATTGAGTTCTATTAAACCGTCAACTATGCCGGATACATACTGGATTACACTGCAGGCACGCCTTATGAGAGAACGTGCCGCAGTGACGGGCTAATTATCTGAAAGAATTTGTGAGGCTGTATCGGTTACTAATTGATTTGATGGTTTTACTCTCGGGAGaata                                                                                   >  2:204603/1‑250 (MQ=255)                                                                                                                                                                  gcggaggagCAGATCGCATCTCATCATAATGTCGTATCTCCTGGGGTGTTATACAAGATATCGTTGTTGGTGACCTGGGAGAGAAATTGAGTTCTATTAAACCGTCAACTATGCCGGATACATACTTGATGACACTGCAGGCACGCCTTATGAGAGAACGTGCCGCAGTGACGGGTTAATTATCTGAAAGAATGTGTGAGGCTGTATCGGTTACTCATTGATTTGATAGTTTTACTCTCGGGAGAataata                                                                                <  2:375613/248‑1 (MQ=255)                                                                                                                                                                     gaggagCATATCGCATCTCATCATAATGTCGTATCTCCTGGGGTGTTATACAAGATATCGTTGTTGGTGACCTGGGAGAGGAATTGAGTTCTATTAAACCGTCAACTATGCCGGATACATACTGGATTACACTGCAGGCACGCCTTATGAGAGAACGTGCCGCAGTGACGGGTTAATTATCTGAAAGAATTTGTGAGGCTGTAGCGGTTACTCATTGATTTGATAGTTTTACTCTCGGGAGAATAATAg                                                                               <  1:135017/249‑1 (MQ=255)                                                                                                                                                                          gCATATCGCATCTCATCATAATGTCGTATCTCCTGGGGTGTTATACAAGATATCGTTGTTGGTGACCTGGGAGAGGAATTGAGTTCTATTAAACCGTCAACTATGCCGGATACATACTGGATTACACTGCAGGCACGCCTTATGAGAGAACGTGCCGCAGTGACGGGTTAATTATCTGAAAGAATTTGTGAGGCTGTATCGGTTACTCATTGATTTGATAGTTTTACTCTCGGGAGAATAATAGATATTTa                                                                        <  1:147987/251‑1 (MQ=255)                                                                                                                                                                                   aTCTCATCATAATGTCGTATCTCCTGGGGTGTTATACAAGATATCGTTGTTGGTGACCTGGGAGAGAAATTGAGTTCTATTAAACCGTCAACTATGCCGGATACATACTGGATTACACTGCAGGCACGCCTTATGAGAGAACGTGCCGCAGTGACGGGTTAATTATCTGAAAGAATTTGTGAGGCTGTATc                                                                                                                           >  1:427865/1‑191 (MQ=255)                                                                                                                                                                                   aTCTCATCATAATGTCGTATCTCCTGGGGTGTTATACAAGATATCGTTGTTGGTGACCTGGGAGAGAAATTGAGTTCTATTAAACCGTCAACTATGCCGGATACATACTGGATTACACTGCAGGCACGCCTTATGAGAGAACGTGCCGCAGTGACGGGTTAATTATCTGAAAGAATTTGTGAGGCTGTATc                                                                                                                           <  2:427865/191‑1 (MQ=255)                                                                                                                                                                                       catcatAATGTCGTATCTCCTGGGGTGTTATACAAGATATCGTTGTTGGTGACCTGGGAGAGGAATTGAGTTCTATTAAACCGTCAACTATGCCGGATACATACTGGATTACACTGCAGGCACGCCTTATGAGAGAACGTGCCGCAGTGACGGGTTAATTATCTGAAAGAATTTGTGAGGCTGTATCGGTTACTCATTGATTTGATAGTTTTa                                                                                                 <  2:331571/213‑1 (MQ=255)                                                                                                                                                                                       catcatAATGTCGTATCTCCTGGGGTGTTATACAAGATATCGTTGTTGGTGACCTGGGAGAGGAATTGAGTTCTATTAAACCGTCAACTATGCCGGATACATACTGGATTACACTGCAGGCACGCCTTATGAGAGAACGTGCCGCAGTGACGGGTTAATTATCTGAAAGAATTTGTGAGGCTGTATCGGTTACTCATTGATTTGATAGTTTTa                                                                                                 >  1:331571/1‑213 (MQ=255)                                                                                                                                                                                               tGTCGTATCTCCTGGGGTGTTATACAAGATATCGTTGTTGGTGACCTGGGAGAGAAATTGAGTTCTATTAAACCGTCAACTATGCCGGATACATACTGGATTACACTGCAGGCACGCCTTATGAGAGAACGTGCCGCAGTGACGGGTTAATTATCTGAAAGAATTTGTGAGGCTGTATCGGTTACTCATTGATTTGATAGTTTTACTCTCGGGAGAATAATAGATATTTAATCCATTAACGGAAAccagcc                                                   >  1:248180/1‑251 (MQ=255)                                                                                                                                                                                                      tctcCTGGGGTGTTATACAAGATATCGTTGTTGGTGACCTGGGGGCGGAATTGAGTTCTATTAAACCGTCAACTATGCCGGATACATACTGGATTACACTGCAGGCACGCCTTATGAGAGAACGTGCCGCAGTGACGGGTTAATTATCTGAAAGAATTTGTGAGGCTGTATCGGTTACTCATTGATTTGATAGTTTTACTCTCGGGAGAATAATAGATATTTAATCCATTAACGGAAACCAGCCAGTTCCt                                            >  2:182131/1‑251 (MQ=255)                                                                                                                                                                                                         ccttgggTGTTATACAAGATATCGTTGGTGGTGACCTGGGAGAGAAATTGAGTTCTATTAAACCGTCAACTATGCCGGATACATACTGGATTACACTGCAGGCACGCCTTATGAGAGAACGTGCCGCAGTGACGGGTTAATTATCTGAAAGAATTTGTGAGGCTGTATCGGTTACTCATTGATTTGATAGTTTTACTCTCGGGAGAATAATAGATATTTAATCCATTAACGGAAACCAGCCAGTTCCttt                                          <  1:247205/246‑1 (MQ=255)                                                                                                                                                                                                            ggggTGTTATACAAGATATCGTTGTTGGTGACCTGGGAGAGGAATTGAGTTCTATTAAACCGTCAACTATGCCGGATACATACTGGATTACACTGCAGGCACGCCTTATGAGAGAACGTGCCGCAGTGACGGGTTAATTATCTGAAAGAATTTGTGAGGCTGTATCGGTTACTCATTGATTTGATAGTTTTACTCTCGGGAGAATAATAGATATTTAATCCATTAACGGAAACCAGCCAGTTCCTTTCGAt                                      <  1:89272/251‑1 (MQ=255)                                                                                                                                                                                                               gtgtTATACAAGATATCGTTGTTGGTGACCTGGGAGAGAAATTGAGTTCTATTAAACCGTCAACTATGCCGGATACATACTGGATTACACTGCAGGCACGCCTTATGAGAGAACGTGCCGCAGTGACGGGTTAATTATCTGAAAGAATTTGTGAGGCTGTATCGGTTACTCATTGATTTGATAGTTTTACTCTCGGGAGAATAATAGATATTGAATCCATTAACGGAACCCAGCCAGTTCCTTTCGATGcc                                   <  2:49100/251‑1 (MQ=255)                                                                                                                                                                                                                  ttATACAAGATATCGTTGTTGGTGACCTGGGAGAGGAATTGAGTTCTATTAAACCGTCAACTATGCCGGATACATACTGGATTACACTGCAGGCACGCCTTATGAGAGAACGTGCCGCAGTGACGGGTTAATTATCTGAAAGAATTTGTGAGGCTGTATCGGTTACTCATTGATTTGATAGTTTTACTCTCGGGAGAATAATAGATATTTAATCCATTAACGGAAACCAGCCAGTTCCTTTCGATGCCTg                                 >  2:57325/1‑250 (MQ=255)                                                                                                                                                                                                                      aCAAGATATCGTTGTTGGTGACCTGGGAGAGAAATTGAGTTCTATTAAACCGTCAACTATGCCGGATACATACTGGATTACACTGCAGGCACGCCTTATGAGAGAACGTGCCGCAGTGACGGGTTAATTATCTGAAAGAATTTGTGAGGCTGTATCGGTTACTCATTGATTTGATAGTTTTACTCTCGGGAGAATAATAGata                                                                            <  2:374054/203‑1 (MQ=255)                                                                                                                                                                                                                      aCAAGATATCGTTGTTGGTGACCTGGGAGAGAAATTGAGTTCTATTAAACCGTCAACTATGCCGGATACATACTGGATTACACTGCAGGCACGCCTTATGAGAGAACGTGCCGCAGTGACGGGTTAATTATCTGAAAGAATTTGTGAGGCTGTATCGGTTACTCATTGATTTGATAGTTTTACTCTCGGGAGAATAATAGata                                                                            >  1:374054/1‑203 (MQ=255)                                                                                                                                                                                                                           atatCGTTGTTGGTGACCTGGGAGAGAAATTGAGTTCTATTAAACCGTCAACTATGCCGGATACATACTGGATTACACTGCAGGCACGCCTTATGAGAGAACGTGCCGCAGTGACGGGTTAATTATCTGAAAGAATTTGTGAGGCTGTATCGGTTACTCATTGATTTGATAGTTTTACTCTCGGGAGAATAATAGATATTTAATCCATTAACGGAAACCAGCCAGTTCCTTTCGATGCCTGAATTTGATcc                       >  2:156505/1‑251 (MQ=255)                                                                                                                                                                                                                             atCGTTGTTGGTGACCTGGGAGAGGAATTGAGTTCTATTAAACCGTCAACTATGCCGGATACATACTGGATTACACTGCAGGCACGCCTTATGAGAGAACGTGCCGCAGTGACGGGTTAATTATCTGAAAGAATTTGTGAGGCTGTATCGGTTACTCATTGATTTGATAGTTTTACTCTCGGGAGAATAATAGATATTTAATCCATTAACGGAAACCAGCCAGTTCCTTTCGATGCCTGAATTTGATcca                      >  1:248585/1‑249 (MQ=255)                                                                                                                                                                                                                                            cTGGGAGAGAAATTGAGTTCTATTAAACCGTCAACTATGCCGGATACATACTGGATTACACt                                                                                                                                                                                                   <  1:224112/62‑1 (MQ=255)                                                                                                                                                                                                                                            cTGGGAGAGAAATTGAGTTCTATTAAACCGTCAACTATGCCGGATACATACTGGATTACACt                                                                                                                                                                                                   >  2:224112/1‑62 (MQ=255)                                                                                                                                                                                                                                                gagaGGAATTGAGTTCTATTAAACCGTCAACTATGCCGGATACATACTGGATTACACTGCAGGCACGCCTTATGAGAGAACGTGCCGCAGTGACGGGTTAATTATCTGAAAGAATTTGTGAGGCTGTATCGGTTACTCATTGATTTGATAGTTTTACTCTCGGGAGAATAATAGATATTTAATCCATTAACGGAAACCAGCCAGTTCCTTTCGATGCCTGAATTTGATCCCATAGTTTATCCGGAATGGtt  <  2:153882/251‑1 (MQ=255)                                                                                                                                                                                                                                                     |                                                                                                                                                                                                                                                      CGAGTTTGGTCAACCAATTTTTTGACATGTATCACAAATTTGAATAGATGTATTACATCAACTATCTTTTATTGCACCAACGTCATTGATATATGTCGCCTGAAGTCAGTTCCGGGAATGAGTCTGATCTCAAGACTGGCCCAGTCCGGGCGTTGATTGGTGCTGAGGAGCATATCGCATCTCATCATAATGTCGTATCTCCTGGGGTGTTATACAAGATATCGTTGTTGGTGACCTGGGAGAGGAATTGAGTTCTATTAAACCGTCAACTATGCCGGATACATACTGGATTACACTGCAGGCACGCCTTATGAGAGAACGTGCCGCAGTGACGGGTTAATTATCTGAAAGAATTTGTGAGGCTGTATCGGTTACTCATTGATTTGATAGTTTTACTCTCGGGAGAATAATAGATATTTAATCCATTAACGGAAACCAGCCAGTTCCTTTCGATGCCTGAATTTGATCCCATAGTTTATCCGGAATGGTT  >  NC_000913/313299‑313788

Alignment Legend
Aligned base mismatch/match (shaded by quality score): ATCG/ATCG < 3 ≤ ATCG/ATCG < 14 ≤ ATCG/ATCG < 17 ≤ ATCG/ATCG < 35 ≤ ATCG/ATCG < 39 ≤ ATCG/ATCG
Unaligned base: atcg    Masked matching base: atcg    Alignment gap: ‑    Deleted base: ‑