Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	NC_000913	47,099	G→A	Q431Q (CAG→CAA)	yaaU →	putative transporter YaaU

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	47,099	0	G	A	98.2%	152.8 / NA	54	Q431Q (CAG→CAA)	yaaU	putative transporter YaaU
Reads supporting (aligned to +/- strand):  ref base G (1/0);  new base A (33/20);  total (34/20)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

GGGCCGGGTAATTTGCAGTGGCTCTATCCTAATGAACTCTTCCCGACAGATATCCGCGCCTCTGCCGTGGGCGTGATTATGTCCTTAAGTCGTATTGGCACCATTGTTTCGACCTGGGCACTACCGATCTTTATCAATAATTACGGTATCAGTAACACGATGCTAATGGGGGCGGGTATCTCGCTGTTTGGCTTGTTGATTTCCGTAGCGTTTGCCCCGGAGACTCGAGGGATGTCACTGGCGCAGACCAGCAATATGACGATCCGCGGGCAGAGAATGGGGTAAATTGTTCAGATTTCTCTCTTTTCTGAATCAATATTATTGACTATAAGCCGCGTGAATATATGACTACACTTTGTGGGAAAACAAAGGCGTAATCACGCGGGCTACCTATGATTCTTATAATTTATGCGCATCCGTATCCGCATCATTCCCATGCGAATAAACGGATGCTTGAACAGGCAAGGACGCTGGAAGGCGTCGAAATTCGCT  >  NC_000913/46854‑47345                                                                                                                                                                                                                                                      |                                                                                                                                                                                                                                                       gggCCGGGTAATTTGCAGTGGCTCTATCCTAATGAACTCTTCCCGACAGATATCCGCGCCTCTGCCGTGGGCGTGATTATGTCCTTAAGCCGTATTGGCACCATTGTTTCGACCTGGGCACTGCCGATCTTTATCAATAATTACGGTATCAGTAACACGATGCTAATGGGGGCGGGTATCTCGCTGTTTGGCTTGTTGATTTCCGTAGCGTTTGCCCCGGAGACTCGAGGGATGTCACTGGCGCAAAcc                                                                                                                                                                                                                                                     >  1:424084/1‑249 (MQ=255)    ccGGGCAATTTGCAGTGGCTCTATCCTAATGACCTCTTCCCGACAGATATCCGCGCCTCTGCCGTGGGCGTGATTATGTCCTTAAGCCGTATTGGCACCATTGTTTCGACCTGGGCACTGCCGATCTTTATCAATAATTACGGTATCAGTAACACGATGCTAATGGGGGCGGGTATCTCGCTGTTTGGCTTGTTGATTTCCGTAGCGTTTGCCCCGGGGACTCGAGGGATGTCACTGGCGCAAACCAGCa                                                                                                                                                                                                                                                 >  2:625085/1‑250 (MQ=255)      gggTAATTTGCAGTGGCTCTATCCTAATGAACTCTTCCCGACAGATATCCGCGCCTCTGCCGTGGGCGTGATTATGTCCTTAAGCCGTATTGGCACCATTGTTTCGACCTGGGCACTGCCGATCTTTATCAATAATTACGGTATCAGTAACACGATGCTAATGGGGGCGGGTATCTCGCTGTTTGGCTTGTTGATTTCCGTAGCGTTTGCCCCGGAGACTCGAGGGATGTCACTGGCGCAAACCAGCAata                                                                                                                                                                                                                                              <  1:134995/251‑1 (MQ=255)         tAATTTGCAGTGGCTCTGTCATAATGACATCTTCCCGACAGATATCCGCGCCTATGCCGTGGGCGTGATTATGTCCTTAAGCCGTATTGGCACCATTGTTTCGACCTGGGCCCTGCCGATCTTTATCAATAATTACGGTATCAGTAACACGATGCTAATGGGGGCGGGTATCTCGCTGTTTGGCTTGTTGATTTCCGTAGCGTTTGCCCCGGAGACTCGAGGGATGTCACTGGCGCAAACCAGCAATATg                                                                                                                                                                                                                                            >  2:319053/1‑250 (MQ=255)               gCAGTGGCTCTATCCTAATGAACTCTTCCCGACAGAGATCCGCGCCTCTGCCGTGGGCGTGATTATGTCCTTAAGCCGTATTGGCACCATTGTTTCGACCTGGGCACTGCCGATCTTTATCAATAATTACGGTATCAGTAACACGATGCTAATGGGGGCGGGTATCTCGCTGTTTGGCTTTTTGATTTCCGTAGCGTTTGCCCCGGAGACTCGAGGGATGTCACTGGCGCAAGCCAGCAATATGACGATcc                                                                                                                                                                                                                                     <  1:301201/251‑1 (MQ=255)                     gCTCTATCCTAATGAACTCTTCCCGACAGATATCCGCGCCTCTGCCGTGGGCGTGATTATGTCCTTAAGCCGTATTGGCACCATTGTTTCGACCTGGGCACTGCCGATCTTTATCAATAATTACGGTATCAGTAACACGATGCTAATGGGGGCGGGTATCTCGCTGTTTGGCTTGTTGATTTCCGTAGCGTTTGCCCCGGAGACTCGAGGGATGTCACTGGCGCAAACCAGCAATATGACGATCCGCGGGc                                                                                                                                                                                                                               <  2:424084/251‑1 (MQ=255)                                 tGAACTCTTCCCGACAGATATCCGCGCCTCTGCCGTGGGCGTGATTATGTCCTTAAGCCGTATTGGCACCATTGTTTCGACCTGGGCACTGCCGATCTTTATCAATAATTACGGTATCAGTAACACGATGCTAATGGGGGCGGGTATCTCGCTGTTTGGCTTGTTGATCTCCGTAGCGTTTGCCCCGGAGACTCGAGGGATGTCACTGGCGCAAACCAGCAATATGACGATCCGCGGGCAGAGAATGGGGt                                                                                                                                                                                                                   >  2:241300/1‑251 (MQ=255)                                    aCTCTTCCCGACAGATATCCGCGCCTCTGCCGTGGGCGTGATTATGTCCTTAAGCCGTATTGGCACCATTGTTTCGACCTGGGCACTGCCGATCTTTATCAATAATTACGGTATCAGTAACACGATGCTAATGGGGGCGGGTATCTCGCTGTTTGGCTTGTTGATTTCCGTAGCGTTTGCCCCGGGGCCTCGAGGGATGTCACTGGCGCAAACCAGCAATATGACGATCCGCGGGCAGGGAATgggggaa                                                                                                                                                                                                                 >  2:123136/1‑247 (MQ=255)                                    aCTCTTCCCGACAGATAGCCGCGCCTCGGCCGTGGGCGTGATTATGTCCTTAAGCCGTATTGGCACCATTGTTTCGACCTGGGCACTGCCGATCTTTATCAATAATTACGGTATCAGTAACACGATGCTAATGGGGGCGGGTATCTCGCTGTTTGGCTTGTTGATTTCCGTAGCGTTTTCCCCGGAGACTCGAGGGATGTCACTGGCGCAAACCAGCAATATGACGATCCGCGGGCATAGAATGGGGTaaa                                                                                                                                                                                                                <  2:315783/251‑1 (MQ=255)                                          cccGAAAGATATCCGCGCCTCTGCCGTGGGCGTGATTATGTCCTTAAGCCGTATTGGCACCATTGTTTCGACCTGGGCACTGCCGATCTTTATCAATAATTACGGTATCAGTAACACGATGCTAATGGGGGCGGGTATCTCGCTGTTTGGCTTGTTGATTTCCGTAGCGTTTGCCCCGGAGACTCGAGGGATGTCACTGGCGCAAACCAGCAATATGACGATCCGCGGGCAGAGAATGGGGTAAATTGTTc                                                                                                                                                                                                          >  2:607837/1‑251 (MQ=255)                                            cGACAGATATCCGCGCCTCTGCCGTGGGCGTGATTATGTCCTGAAGCCGTATGGGCACCATTGTTTCGACCTGGGCACTGCAGATGTTTATCAATAATTGCGGTATCAGGAACACGATGCTAATGGGGGGGGGTATCTCGCTGTTTGGCTTGTTGATGTCCGTAGCGTTTGCCCCGGAGACGCTATGGATGGCACTGGCGCAAACCAGCAATATGACGATCCGCTGGCAGAGAATGGGGTAAATTGTTCAg                                                                                                                                                                                                        <  1:366680/251‑1 (MQ=255)                                                aGATATCCGCGCCTCTGCCGTGGGCGTGATTATGTCCTTAAGCCGTATTGGCACCATTGTTTCGACCTGGGCACTGCCGATCTTTATCAATAATTACGGTATCAGTAACACGATGCTAATGGGTGCGGGTATCTCGCTGTTTGGCTTGTTGATTTCCGGAGCGTTTGCCCCGGAGACTCGAGGTATGTCACTGGCGCAAACCAGCAATATGACGATCCGCGGGCAGAGAATGGGGTAAATTTTTCAGAttt                                                                                                                                                                                                    <  2:49796/251‑1 (MQ=255)                                                         cgcCTCTGCCGTGGGCGTGATTATGTCCTTAAGCCGTATTGGCACCATTGTTTCGACCTGGGCACTGCCGATCTTTATCAATAATTACGGTATCAGTAACACGATGCTAATGGGGGCGGGTATCTCGCTGTTTGGCTTGTTGATTTCCGTAGCGTTTGCCCCGGAGACTCGAGGGATGTCACTGGCGCAAACCAGCAATATGACGATCCGCGGGCAGAGAATGGGGTAAATTGTTCAGATTTCTCTCtttt                                                                                                                                                                                           >  1:151563/1‑251 (MQ=255)                                                             tctGCCGTGGGCGTGATTATGTACTTAAGCCGTCTTGGCAGCATTGCTTCGACCTGGGCACTGCAGATCTTTATCAGTAATTACGGTATCAGTACCACGATGCTACTGGGGGCGGGTATCTCGCTGTTTGGCCTGTTGATTCCCGTAGCGTTTGCCACGGAGACTGGAGGGATGTCACTGGCGCAAACCAGCAATATGACGATCCGCGGGCAGAGAATGGGGAAACCTGTTCAGATTTCACTCTTTTCTg                                                                                                                                                                                        >  2:299733/1‑250 (MQ=255)                                                                gCCGTGGGCGTGATTATGTCCTTAAGCCGTATTGGCACCATTGTTTCGACCTGGGCACTGCCGATCTTTATCAATAATTACGGTATCAGTAACACGATGCTAATGGGGGCGGGTATCTCGCTGTTTGGCTTGTTGATTTCCGTAGCGTTTGCCCCGGAGACTCGAGGGATGTCACTGGCGCAAACCAGCAATATGACGATCCGCGGGCAGAGAATGGGGTAAATTGTTCAGATTTCTCTCTTTTCTGa                                                                                                                                                                                       <  1:9264/248‑1 (MQ=255)                                                                   gTGGGCGTGATTATGTCCTTAAGCCGTATTGGCACCATTGTTTCGACCTGGGCACTGCCGATCTTTATCAATAATTACGGTATCAGTAACACGATGCTAATGGGGGCGGGTATCTCGCTGTTTGGCTTGTTGATTTCCGTAGCGTTTGCCCCGGAGACTCGAGGGATGTCACTGGCGCAAACCAGCAATATGACGATCCGCGGGCAGAGAATGGGGTAAATTGTTCAGATTTCTCTCTTTTCTGAATCAat                                                                                                                                                                                 <  1:16945/251‑1 (MQ=255)                                                                          tGATTATGTCCTTAAGCCGTATTGGCACCATTGTTTCGACCTGGGCACTGCCGATCTTTATCAATAATTACGGTATCAGTAACACGATGCTAATGGGGGCGGGTATCTCGCTGTTTGGCTTGTTGATTTCCGTAGCGTTTGCCCCGGAGACTCGAGGGATGTCACTGGCGCCAACCAGCAATATGACGCTCAGCGGGCAGAGAATGGGGTAAATTGTTCAGATTTCTCTCTTTTCTGAATCAATATTATTg                                                                                                                                                                          >  2:464690/1‑251 (MQ=255)                                                                           gATTATGTCCTTAAGCCGTATTGGCACCATTGTTTCGACCTGGGCACTGCCGCTCTTTATCAATAATTACGGTATCAGTAACACGATGCTAATGGGGGCGGGTATCTCGCTGTTTGGCTTGTTGATTTCCGTAGCGTTTGCCCCGGAGACTCGAGGGATGTCACTGGCGCAAACCAGCAATATGACGATCCGCGGGCAGAGAATGGGGTAAATTGTTCAGATTTCTCTCTTTTCTGAATCAATATTATTg                                                                                                                                                                          >  1:49930/1‑250 (MQ=255)                                                                               aTGTCCTTAAGCCGTATTGGCACCATTGTTTCGACCTGGGCACTGCCGATCTTTATCAATAATTACGGTATCAGTAACACGATGCTAATGGGGGCGGGTATCTCGCTGTTTGGCTTGTTGATTTCCGTAGCGTTTGCCCCGGAGACTCGTGGGATGTCACTGGCGCAAACAAGCAATATGACGATCCGCGGGCAGAGAATGGGGTAAATTGTTCAGATTTCTTTCTTTTCTGAATCAATATTATTGACtat                                                                                                                                                                     >  2:11533/1‑251 (MQ=255)                                                                               aTGTCCTTAAGCCGTATTGGCACCATTGTTTCGAACTGGGCACTGCCGATCTTTATCAATAATTACGGTATCAGTAACACGATGCTAATGGGGGCGGGTATCTCGCTTTTTGGCTTGTTGATTTCCGTAGCGTTTGCCCCGGAGACTCGAGGGATGTCACTGGCGAAAACCAGCAATATGACGATCCGCGGGCAGAGAATGGGGTAAATTGTTCAGATTTCTCTCTTTTCTGAATCAATATTATTGACtat                                                                                                                                                                     <  2:397168/251‑1 (MQ=255)                                                                                      tAAGCCGTATTGGCACCATTGTTTCGACCTGGGCACTGCCGATCTTTATCAATAATTACGGTATCAGTAACACGATGCTAATGGGGGCGGGTATCTCGCTGTTTGGCTTGTTGATTTCCGTAGCGTTTGCCCCGGAGACTCGAGGGATGTCACTGGCGCAAACCAGCAATATGCCGATCCGCGGGCAGAGAATGGGGTAAATTGTTCAGATTTCTCTCTTTTCTGAATCAATATTATTGACTATAAGCcgc                                                                                                                                                              >  2:83614/1‑251 (MQ=255)                                                                                       aagccGTATTGGCACCATTGTTTCGACCTGGGCACTGCCGATCTTTATCAATAATTACGGTATCAGTAACACGATGCTACTGGGGGCGGGTATCTCGCTGTTTGGCTTGTTGATTTCCGTAGCGTTTGCCCCGGAGACTCGAGGGATGTCACTGGCGCAAACCAGCAATATGACGCTCCGCGGGCAGAGAATGGGGTAAATTGTTCAGATTTCTATCTTTTCTGAATCAATATTATTGACTATAAGCcgcg                                                                                                                                                             >  2:57603/5‑251 (MQ=255)                                                                                           cGTATTGGCACCATTGTTTCGACCTGGGCACTGCCGATCTTTATCAATAATTACGGTATCAGTAACACGATGCTAATGGGGGCGGGTATCTCGCTGTTTGGCTTGTTGATTTCCGTAGCGTTTGCCCCGGAGACTCGAGGGATGTCACTGGCGCAAACCAGCAATATGACGATCCGCGGGCAGAGAATGGGGTAAATTGTTCAGATTTCTCTCTTTTCTGAATCAATATTATTGACTATAAGCCGCGTGa                                                                                                                                                          >  1:352475/1‑250 (MQ=255)                                                                                                 ggCACCATTGTTTCGCCCTGGGCACTGCCGATCTTTATCAATAATTACGGTATCAGTAACACGATGCTAATGGGGGCGGGTATCTCGCTGTTTGGCTTGTTGATTTCCGTAGCGTTTGCCCCGGAGACTCGAGGGATGTCACTGGCGCAAACCAGCAATATGCCGATCCGCGGGCCGAGAATGGGGTAAATTGTTCAGATTTCTCTCTTTTCTGAATCAATATTATTGACTATAAGCCGCGTGAATATATg                                                                                                                                                   >  2:526516/1‑251 (MQ=255)                                                                                                     ccATTGTTTCGACCTGGGCACTGCCGATCTTTATCAATAATTACGGTATCAGTAACACGATGCTAATGGGGGCGGGTATCTCGCTGTTTGGCTTGTTGATTTCCGTAGCGTTTGCCCCGGAGACTCGAGGGATGTCACTGGCGCAAACCAGCAATATGACGATCCGCGGGCAGAGAATGGGGTAAATTGTTCAGATTTCTCTCTTTTCTGAATCAATATTATTGACTATAAGCCGCGTGAATATATGACt                                                                                                                                                >  2:235359/1‑250 (MQ=255)                                                                                                      cATTGTTTCGACCTGGGCACTGCCGATCTTTATCAATAATTACGGTATCAGTAACACGATGCTAATGGGGGCGGGTATCTCGCTGTTTGGCTTGTTGATTTCCGTAGCGTTTGCCCCGGAGACTCGAGGGATGTCACTGGCGCAAACCAGCAATATGACGATCCGCGGGCAGAGAATGGGGTAAATTGTTCAGATTTCTCTCTTTTCTGAATCAATATTATTGACTATAAGCCGCGTGAATATATGACTac                                                                                                                                              >  1:122068/1‑251 (MQ=255)                                                                                                          ggttCGACCTGGGCACTGCCGATCTTTATCAATAATTACGGTATCAGTAACACGATGCTAATGGGTGCGGGTATCTCGCTGTTTGGCTTGTTGATTTCCGTAGCGTTTGCCCCGGAGACTCGAGGGATGTCACGGGCGCAATCCAGCAAGATGACGATCCGCGGGCAGATAATGGGGTAAATTGTTCAGATTGCTCTCTGTACTGAATCAATATTATTGACTATAAGCCGCGTGAATATATGACTACACtt                                                                                                                                          <  2:90491/249‑1 (MQ=255)                                                                                                                  ctgggccatACCGATCTTTATCCATAATTCCGGTATCAGTAACACGATGCTAATGGGGGCGGGTATATCGCTGTTTGGCTTGATGATTTCCGTCGCGTTTGCCCCGGAGACTCGCGGGATGTCACTGGCGAAAACCAGCAATATGACGATCCGCGGGCAGAGAATGGGGTAAATTGTTCAGATTTCTCTCTTTTCTGAATCAATATTATTGACTATAAGCCGCGGAAATATATGACTACACTATGTGGGa                                                                                                                                   >  2:96859/9‑250 (MQ=255)                                                                                                                   tGGGCACTGCCGATCTTTATCAATAATTACGGTATCAGTAACACGATGCTAATGGGGGCGGGTATCTCGCTGTTTGGCTTGTTGATTTCCGTAGCGTTTGCCCCGGAGACTCGAGGGATGTCACTGGCGCAAACCAGCAATATGACGATCCGCGGGCAGAGAATGGGGTAAATTGTTCAGATTTCTCTCTTTTCTGAATCAATATTATTGACTATAAGCCGCGTGAATATATGACTACACTTTGTGGGaa                                                                                                                                  >  1:510918/1‑250 (MQ=255)                                                                                                                        actgccGATCTTTATCAATAATTACGGTATCAGTAACACGATGCTAATGGGGGCGGGTATCTCGCTGTTTGGCTTGTTGATTTCCGTAGCGTTTGCCCCGGAGACTCGAGGGATGTCACTGGCGCAAACCAGCAATATGACGATCCGCGGGCAGAGAATGGGGTAAATTGTTCAGATTTCTCTCTTTTCTGAATCAATATTATTGACTATAAGCCGCGTGAATATATGACTACACTTTGTGGGAAAACaaa                                                                                                                            <  2:352475/247‑1 (MQ=255)                                                                                                                                   ttATCAATAATTACTGTATCAGTAACACGATGCTAATGGGGGCGGGTATCTCGCTGTTTGGCTTGTTGATTTCCGTAGCGTTTGCCCCTGAGACTCGAGGGATGTCACTGGCGCAAACCAGCAATATGACGATCCGCGGGCAGAGAATGGGGTAAATTGTTCAGATTTCTCTCTTTTCTGAATCAATATTATTTACTATAAGCCGCGTGAATATATGACTACACTTTGTGGGAAAACAAAGGCGTa                                                                                                                      <  1:42081/246‑1 (MQ=255)                                                                                                                                   ttATCAATAATTACGGTATCAGTAACACGATGCTAATGGGGGCGGGTATCTCGCTGTTTGGCTTGTTGATTTCCGTAGCGTTTGCCCCGGAGACTCGAGGGATGTCACTGGCGCAAACCAGCAATATGACGATCCGCGGGCAGAGAATGGGGTAAATTGTTCAGATTTCTCTCTTTTCTGAATCAATATTATTGACTATAAGCCGCGTGAATATATGACTACACTTTGTGGGAAAACAAAGGCGt                                                                                                                       <  1:130822/245‑1 (MQ=255)                                                                                                                                   ttATCAATAATTACGGTATCAGTAACACGATGCTAATGGGGGCGGGTATCTCGCTGTTTGGCTTGTTGATTTCCGTAGCGTGTGCCCCGGAGACTCGAGGGATGTCACTGGCGCAAACCAGCAATATGACGATCCGCGGGCAGAGAATGGGGTAAATTGTTCAGATTTCTCTCTTTTCTGAATCAATATTATTGACTATAAGCCGCGTGAATATATGACTACACTTTGTGGGAAAACAAAGGCGt                                                                                                                       <  1:542893/245‑1 (MQ=255)                                                                                                                                        aataatTTCGGTATAAGTAACACGATGCTAATGGGGGCGGGTATCTCGCTGTTTGGCTTGTTGATTTCCGAAGCGTTTGCCCCGGCGACTCGTGGGATGTCACTGGCGCAAACCAGAAATATGACGATCCGCGAGCATAGAATGGGGTAAATTGGTCAGATTTCTATCTTTTCTGAATCAATATTATTGAAGATAAGCAGCGTGACTATATGACTACACTTTGTGGGAAAACAACGGCGTAATAACGCggg                                                                                                            >  2:13056/1‑251 (MQ=255)                                                                                                                                              gtctttATCAGTAACACGATGCTAATGGGGGCGGGTACCGCGCTGTGTGGATTGTTGATTTCCGTTACGCTGGCCTCTGAGACTCGAAGGATGTCACTGGCGCAAACTATCAATATGACTAGCCGCGGGCACAGAACGGGGTAAATTGTTCAGATTTCTCTCTTTTCGGAAGCAATATTATTGACGAGAAGCCGCTTGAATATAGGAGCACGCTTTGTGGTTAAACAAAGGCGTAATCATGTGGGCTAcc                                                                                                       <  2:283298/245‑1 (MQ=255)                                                                                                                                                 ggTATCAGTAACACGATGCTAATGGGGGCGGGTATCTCGCTGTTTGGCTTGTTGATTTCCGTAGCGTTTGCCCCGGAGACTCGAGGGATGTCACTGGCGCAAACCAGCAATATGACGATCCGCGGGCAGAGAATGGGGTAAATTGTTCAGATTTCTCTCTTTTCTGAATCAATATTATTGACTATAAGCCGCGTGAATATATGACTACACTTTGTGGGAAAACAAAGGCGTAATCACGCGGGCTACCTATg                                                                                                   >  1:506551/1‑251 (MQ=255)                                                                                                                                                 ggTATCAGTAACACGATGCTAATGGGGGCGGGTATCGCGCTGTTTGGCTTGTTGATTTCCGTAGCGTTTGCCCCGGAGACTCGATGGATGTCACTGGCGCAAGCCAGCAATATGACGATCCGCGGGCAGAGAATGGGGTAAATTGTTCAGATGTCTCTCTTTTCTGAATCAATATTATTGACTATAAGCCGCGTGAATATATGACTTCACTTTGTGGGAAAACAAAGGCGTAATCACGCGGGCTACCTATg                                                                                                   <  2:299963/251‑1 (MQ=255)                                                                                                                                                  gTATCAGTAACACGATGCTAATGGGGGCGGGTATCTCGCTGTTTGGCTTGTTGATTTCCGTAGCGTTTGCCCCGGAGACTCGAGGGATGTCACTGGCGCAAACCAGCAATATGACGATCCGCGGGCAGAGAATGGGGTAAATTGTTCAGATTTCTCTCTTTTCTGAATCAATATTATTGACTATAAGCCGCGTGAATATATGACTACACTTTGTGGGAAAACAAAGGCGTAATCACGCGGGCTACCTATgg                                                                                                  >  1:226039/1‑250 (MQ=255)                                                                                                                                                     tCAGTAACACGATGCTAATGGGGGCGGGTATCTCGCTGTTTGGCTTGTTGATTTCCGTAGCGTTTGCCCCGGAGACTCGAGGGATGTCACTGGCGCAAACCAGCAATATGACGATCCGCGGGCAGAGAATGGGGTAAATTGTTCAGATTTCTCTCTTTTCTGAATCAATATTATTGACTATAAGCCGCGTGAATATATGACTACCCTTTGTGGGAAAACAAAGGCGTAATCACGCGGGCTCCCTATGATTc                                                                                               >  1:145744/1‑251 (MQ=255)                                                                                                                                                     tCAGTAACACGATGCTAATGGGGGCGGGTATCTCGCTGTTTGGCTTGTTGATTTCCGTAACGTTTGCCCCGGAGACTCGAGGGATGTCACTGGCGCAAACCAGCAATATGACGATCCGCGGGCAGAGAATGGGGTAAATTGTTCAGATTTCTCTCTTTTCTGAATCAATATTATTGACTATAAGCCGCGTGAATATATGACTACACTTTGTGGGAAAACAAAGGCGTAATCACGCGGGCTACCTATGATTc                                                                                               >  1:271998/1‑251 (MQ=255)                                                                                                                                                              cGATGCTTATGGGGGCGGGTATCTCGCTGTGTGGCTTGTTGATTTCCATAGCGTTTGCCCCGGCGACTCGAGGGATGTCACTGGCGCAGACCAGCAATATGACGATCCGCGGGCAGAGAATGGGGTCAATTGTTCGGGTTTCTCTCTTTTCTGAATCAATATTATTGACTATAAGCCGCGTGAATATATGACTACACTTTGTGGGAAAACAAAGGCGTAATCACGCGGGCTACCTATGATTCTTATAAttt                                                                                      >  2:60133/1‑251 (MQ=255)                                                                                                                                                                     aaTGGGGGCGGGTATCTCGCTGTTTGGCTTGTTGATTTCCGTAGCGTTTGCCCCGGAGACTCGAGGGATGTCACTGGCGCAAACCAGCAATATGACGATCCGCGGGCAGAGAATGGGGTAAATTGTTCAGATTTCTCTCTTTTCTGAATCAATATTATTGACTATAAGCCGCGTGAATATATGACTACACTTTGTGGGAAAACAAAGGCGTAATCACGCGGGCTACCTATGATTCTTATAATTTATgcgc                                                                                >  1:118345/1‑250 (MQ=255)                                                                                                                                                                        gggggCGGGTATCTCGCTGGTTGGCTTGTTGATTTCCGTAGCGTTTGCCCCGGAGACTCGAGGGATGTCACTGGCGCAAACCAGCAATATGACGATCCGCGGGCAGAGAATGGGGTAAATTGTTCAGATTTCTCTCTTTTCTGAATCAATATTATTGACTATAAGCAGCGTGAATATATGCCTACACTTTGTGGGGAAACACAGGCGTAATCACGCGGGCTCCCTATGATTCTTATAATTTATGCGCATcc                                                                            >  2:189267/1‑251 (MQ=255)                                                                                                                                                                                gTATCTCGCTGTTTGGCTTGTTGATTTCCGTAGCGTTTGCCCCGGAGACTCGAGGGATGTCACTGGCGCAAACCAGCAATATGACGATCCGCGGGCAGAGAATGGGGTAAATTGTTCAGATTTCTCTCTTTTCTGAATCAATATTATTGACTATAAGCCGCGTGAATATATGACTACACTTTGTGGGAAAACAAAGGCGTAATCACGCGGGCTACCTATGATTCTTATAATTTATGCGCATCCGTATCCGc                                                                    >  2:529812/1‑251 (MQ=255)                                                                                                                                                                                    ctcGCTGTTTGGCTTGTTGATTTCCGTAGCGTTTGCCCCGGAGACTCGAGGGATGTCACTGGCGCAAACCAGCAATATGACGATCCGCGGGCAGAGAATGGGGTAAATTGTTCAGATTTCTCTCTTTTCTGAATCAATATTATTGACTATAAGCCGCGTGAATATATGACTACACTTTGTGGGAAAACAAAGGCGTAATCACGCGGGCTACCTATGATTCTTATAATTTATGCGCATCCGTATCCGcatc                                                                 >  2:237923/1‑250 (MQ=255)                                                                                                                                                                                    ctcGCTGTTTGGCTTGTTGATTTCCGTAGCGTTTGCCCCGGAGACTCGAGGGATGTCACTGGCGCAAACCAGCAATATGACGATCCGCGGGCAGAGAATGGGGTAAATTGTTCAGATTTCTCTCTTTTCTGAATCAATATTATTGACTATAAGCCGCGTGAATATATGACTACACTTTGTGGGAAAACAAAGGCGTAATCACGCGGGCTACCTATGATTCTTATAATTTATGCGCATCCGTATCCGcatc                                                                 >  1:327379/1‑250 (MQ=255)                                                                                                                                                                                         tGTTTGGCTTGTTGATTTCCGTAGCGTTTGCCCCGGAGACTCGAGGGATGTCACTGGCGCAAACCAGCAATATGACGATCCGCGGGCAGAGAATGGGGTAAATTGTTCAGATTTCTCTCTTTTCTGAATCAATATTATTGACTATAAGCCGCGTGAATATATGACTACACTTTGTGGGAAAACAAAGGCGTAATCACGCGGGCTACCTATGATTCTTATAATTTATGCGCATCCGTATCCGCATCATTccc                                                           >  2:95452/1‑251 (MQ=255)                                                                                                                                                                                                 ttgttgATTTCCGTAGCGTTTGCCCCGGAGACTCGAGGGATGTAACTGGCGCAAACCAGCAATATGACGATCCGCGGGCAGAGAATGGGGTAAATTGTTCAGATTTCTCTCTTTTCTGAATCAATATTATTGACTATAAGCCGCGTGAATATATGACTACACTTTGTGGGAAAACAAAGGCGTAATCACGCGGGCTACCTATGATTCTTATAATTTATGCGCATCCGTATCCGCATCATTCCCATGCGAAt                                                   >  1:234614/1‑251 (MQ=255)                                                                                                                                                                                                   gttgATTTCCGTAGCGTTTGCCCCGGAGACTCGAGGGATGTCACTGGCGCAAATCAGCAATATGACGATCCGCGGGCAGAGAATGGGGTAAATTGTTCAGATTTCTCTCTTTTCTGAATCAATATTATTGACTATAAGCCGCGTGAATATATGACTACACTTTGTGGGAAAACAAAGGCGTAATCACGCGGGCTACCTATGATTCTTATAATTTATGCGCATCCGTATCCGCATCATTCCCATGCGAATaa                                                 <  1:607837/251‑1 (MQ=255)                                                                                                                                                                                                                    ttGCCCCGGAGACTCGAGGGATGTCACTGGCGCAAACCAGCAATATGACGATCCGCGGGCAGAGAATGGGGTAAATTGTTCAGATTTCTCTCTTTTCTGAATCAATATTATTGACTATAAGCCGCGTGAATATATGACTACACTTTGTGGGAAAACAAAGGCGTAATCACGCGGGCTACCTATGATTCTTATAATTTATGCGCATCCGTATCCGCATCATTCCCATGCGAATAAACGGATGCTTGAACAg                                 <  2:151563/250‑1 (MQ=255)                                                                                                                                                                                                                                        aTGTCACTGGCGCAAACCAGCAATATGACGATCCGCGGGCAGAGACTGGGGCAAATTGTTCAGATTTCTCTCTTTTCTGAATCAATATTATTGACTATAAGCCGCGTGAATATATGACTACACTTTGTGGGAAAACAAAGGCGTAATCACGCGGGCTACCTATGATTCTTCTAATTTATGCGCATCCGTATCCGCATCATTCCCATGCGAATAAACGGATGCTTGAACAGGCAAGGACGCTGGAAGGCGTc            >  1:399547/1‑251 (MQ=255)                                                                                                                                                                                                                                          gTCACTGGCGCAAACCAGCAATATGACGATCCGCGGGCAGAGAATGGGGTAAATTGTTCAGATTTCTCTCTTTTCTGAATCAATATTATTGACTATAAGCCGCGTGAATATATGACTACACTTTGTGGGAAAACAAAGGCGTAATCACGCGGGCTACCTATGATTCTTATAATTTATGCGCATCCGTATCCGCATCATTCCCATGCGAATAAACGGCTGCTTGAACAGGCAAGGACGCTGGAAGGCGTCg           <  1:223585/250‑1 (MQ=255)                                                                                                                                                                                                                                                ggCGCAAACCAGCAATATGACGATCCGCGGGCAGAGAATGGGGTAAATTGTTCAGATTTCTCTCTTTTCTGAATCAATATTATTGACTATAAGCCGCGTGAATATATGACTACACTTTGTGGGAAAACAAAGGCGTAATCACGCGGGCTACCTATGATTCTTATAATTTATGCGCATCCGTATCCGCATCATTCACATGCGAATAAACGGATGCTTGAACAGGCAAGGACGCTGGAAGGCGTCGAAATTc     <  1:83614/250‑1 (MQ=255)                                                                                                                                                                                                                                                  cgcAAACCAGCAATATGACGATCCGCGGGCAGAGAATGGGGTAAATTGTTCAGATTTCTCTCTTTTCTGAATCAATATTATTGACTATAAGCCGCGTGAATATATGACTACACTTTGTGGGAAAACAAAGGCGTAATCACGCGGGCTACCTATGATTCTTATAATTTATGCGCATCCGTATCCGCATCATTCCCATGCGAATAAACGGATGCTTGAACAGGCAAGGACGCTGGAAGGCGTCGAAATTCGct  >  1:527932/1‑251 (MQ=255)                                                                                                                                                                                                                                                      |                                                                                                                                                                                                                                                       GGGCCGGGTAATTTGCAGTGGCTCTATCCTAATGAACTCTTCCCGACAGATATCCGCGCCTCTGCCGTGGGCGTGATTATGTCCTTAAGTCGTATTGGCACCATTGTTTCGACCTGGGCACTACCGATCTTTATCAATAATTACGGTATCAGTAACACGATGCTAATGGGGGCGGGTATCTCGCTGTTTGGCTTGTTGATTTCCGTAGCGTTTGCCCCGGAGACTCGAGGGATGTCACTGGCGCAGACCAGCAATATGACGATCCGCGGGCAGAGAATGGGGTAAATTGTTCAGATTTCTCTCTTTTCTGAATCAATATTATTGACTATAAGCCGCGTGAATATATGACTACACTTTGTGGGAAAACAAAGGCGTAATCACGCGGGCTACCTATGATTCTTATAATTTATGCGCATCCGTATCCGCATCATTCCCATGCGAATAAACGGATGCTTGAACAGGCAAGGACGCTGGAAGGCGTCGAAATTCGCT  >  NC_000913/46854‑47345

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