Predicted mutation
evidence	seq id	position	mutation	freq	annotation	gene	description
RA	NC_000913	4,070,413	T→G	7.4%	intergenic (‑12/+102)	yihR ← / ← yihS	putative sulphoquinovose mutarotase/sulphoquinovose isomerase

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	4,070,413	0	T	G	7.4%	161.6 / 3.1	57	intergenic (‑12/+102)	yihR/yihS	putative sulphoquinovose mutarotase/sulphoquinovose isomerase
Reads supporting (aligned to +/- strand):  ref base T (35/17);  new base G (2/3);  total (37/20)
Fisher's exact test for biased strand distribution p-value = 3.32e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.87e-01

TGTTCATTAATGGGTAACTGATATTCCTGCCCCTGATAGCGGTAACAGCCGTTGGCGATCCGGTTCGGCCACGGAATAAGCACTTTACCCAGATGCGCCAGCGGCATCTCTTCGGGTTTGTGCGGGATGACTAAATGACAACCCTGAAAGGTTAATTCAGCCAGACCAGCACCGACGGTCACAATGGTCGCGTGATAATCCCCCGCGGCAAGTGAAACGGTCTGTCCACTGCAATGCATATTGGTGATTTGCATACAAGGCACCTTAATTAACGACATACAGCCTCCCGTGAGGCTG‑ATGTCAAACCCAAAAACTGCACCGGACGGTCCCCTCGCCCCTTTGGGGAGAGGGTTAGGGAGAGGGGAAAACCCGCACCCTATTTTCAGATCAGTTATTTCGCATTAATATCCAGCAGACCCGCTGCAACCGCTGGAGCCATGCCGGGCGCTAACGGGATACGCGGGATCACCAGGCAATGCAGCAGGTGATAAATATCCTGTTTGCCGTCCCAGACTTTGGTGGTGACCTTATTGTCCGCATCCAGCTCCTGCCACCAGGAACCATTTTCATAGTCCATCAGGTAC  >  NC_000913/4070124‑4070707                                                                                                                                                                                                                                                                                                  |                                                                                                                                                                                                                                                                                                        tGTTCATTAATGGGTAACTGATATTCCTGCCCCTGATAGCGGTAACAGCCGTTGGCGATCCGGTTCGGCCACGGAATAAGCACTTTACCCAGATGCGCCAGCGGCATCTCTTCGGGTTTGTGCGGGATGACTAAATGACAACCCTGAAAGGTTAATTCAGCCAGACCAGCACCGACGGTCACAATGGTCGCGTGATAATCCCCCGCGGCAAGTGAAACGGTCTGTCCACTGCAATGCATATTGGTGATTTGCATACAAGGCACCTTAATTAACGACATACAGCCTCCCGTGAGGCTG‑At                                                                                                                                                                                                                                                                                               <  1:769231/299‑1 (MQ=255)   ttCATTAATGGGTAACTGATATTCCTGCCCCTGATAGCGGTAACAGCCGTTGGCGATCCGGTTCGGCCACGGAATAAGCACTTTACCCAGATGCGCCAGCGGCATCTCTTCGGGTTTGTGCGGGATGACTAAATGACAACCCTGAAAGGTTAATTCAGCCAGACCAGCACCGACGGTCACAATGGTCGCGTGATAATCCCCCGCGGCAAGTGAAACGGTCTGTCCACTGCAATGCATATTGGTGATTTGCATACAAGGCACCTTAATTAACGACATACAGCCTCCCGTGAGGCTG‑ATg                                                                                                                                                                                                                                                                                              >  1:632601/1‑298 (MQ=255)      aTTAATGGGTAACTGATATTCCTGCCCCTGATAGCGGTAACAGCCGTTGGCGATCCGGTTCGGCCACGGAATAAGCACTTTACCCAGATGCGCCAGCGGCATCTCTTCGGGTTTGTGCGGGATGACTAAATGACAACCCTGAAAGGTTAATTCAGCCAGACCAGCACCGACGGTCACAATGGTCGCGTGATAATCCCCCGCGGCAAGTGAAACGGTCTGTCCACTGCAATGCATATTGGTGATTTGCATACAAGGCACCTTAATTAACGACATACAGCCTCCCGTGAGGCTG‑ATGTCaa                                                                                                                                                                                                                                                                                          >  1:181152/1‑299 (MQ=255)       ttAATGGGTAACTGATATTCCTGCCCCTGATAGCGGTAACAGCCGTTGGCGATCCGGTTCGGCCACGGAATAAGCACTTTACCCAGATGCGCCAGCGGCATCTCTTCGGGTTTGTGCGGGATGACTAAATGACAACCCTGAAAGGTTAATTCAGCCAGACCAGCACCGACGGTCACAATGGTCGCGTGATAATCCCCCGCGGCAAGTGAAACGGTCTGTCCACTGCAATGCATATTGGTGATTTGCATACAAGGCACCTTAATTAACGACATACAGCCTCCCGt                                                                                                                                                                                                                                                                                                         >  2:506625/1‑284 (MQ=255)             ggTAACTGATATTCCTGCCCCTGATAGCGGTAACAGCCGGTGGCGATCCGGTTCGGCCACGGAATAAGCACTTGACCAAGATGCGCCAGCGGCATCTCTTCGGGTGTGTGCGGGATGACTAAAGGACAACCCTGAAAGGTTAATTCAGCCAGACCAGCACCGACGGTCACAAGGGTCGCGTGATAATCCCCCGCGGCAAGTGAAACGGTCTGTCCACTGCAATGCATAGTGGTGATTTGCATACAAGGCACCTTAATTAACGACATACAGCCTCCCGTGAGGCTg                                                                                                                                                                                                                                                                                                  <  2:842122/285‑1 (MQ=255)                aaCTGATCTTCCTGCCCCTGATAGCGGTAACAGCCGTTGGCGATCCGGTTCGGCCATGGAATAAGCACTTTACCCAGATGCGCCAGCGGCATCTCTTCGGGTTTGTGCGGGATGACTAAATGACAACCCTGAAAGGTTAATTCAGCCAGACCAGCACCGAAGGTCACAATGGTCGCGTGATAATCCCCCGCGGCAAGTGAAACGGTCTGTCCACTGCAACGCATCTTGTTGATTTGCATACAAGGCACCTTAATCACCGACATCCAGGCTACCGTTCGGCTG‑a                                                                                                                                                                                                                                                                                                >  2:781687/1‑283 (MQ=255)                   tGATATTCCTGCCCCTGATAGCGGTAACAGCCGTTGGCGATCCGGTTCGGCCACGGAATAAGCACTTTACCCAGATGCGCCAGCGGCATCTCTTCGGGTTTGTGCGGGATGACTAAATGACAACCCTGAAAGGTTAATTCAGCCAGACCAGCACCGACGGTCACAATGGTCGCGTGATAATCCCCCGCGGCAAGTGAAACGGTCTGTCCACTGCAATGCATATTGGTGATTTGCATACAAGGCACCTTAATTAACGACATACAGCCTCCCGTGAGGCTG‑ATGta                                                                                                                                                                                                                                                                                            >  2:190044/1‑283 (MQ=255)                   tGATATTCCTGCCCCTGATAGCGGTAACAGCCGTTGGCGATCCGGTTCGGCCACGGAATAAGCACTTTACCCAGATGCGCCAGCGGCATCTCTTCGGGTTTGTGCGGGATGACTAAATGACAACCCTGAAAGGTTAATTCAGCCAGACCAGCACCGACGGTCACAATGGTCGCGTGATAATCCCCCGCGGCAAGTGAAACGGTCTGTCCACTGCAATGCATATTGGTGATTTGCATACAAGGCACCTTAATTAACGACATACAGCCTCACGGGAGGCTGTATGTcgttaattaaggtgc                                                                                                                                                                                                                                                                              >  1:246043/1‑285 (MQ=255)                    gATATTCCTGCCCCTGATAGCGGTAACAGCCGTTGGCGATCCGGTTCGGCCACGGAATAAGCACTTTACCCAGATGCGCCAGCGGCATCTCTTCGGGTTTGTGCGGGATGACTAAATGACAACCCTGAAAGGTTAATTCAGCCAGACCAGCACCGACGGTCACAATGGTCGCGTGATAATCCCCCGCGGCAAGTGAAACGGTCTGTCAACTGCAATGCATATTGGTTGTTTGCATACAAGGCACATTAATTAACGACATTCAGCCTCCCGTGAGGCTG‑ATGTCa                                                                                                                                                                                                                                                                                           >  2:799218/1‑284 (MQ=255)                             gCCCCTGATAGCGGTAACAGCCGTTGGCGATCCGGTTCGGCCACGGAATAAGCACTTTACCCAGATGCGCCAGCGGCATCTCTTCGGGTTTGTGCGGGATGACTAAATGACAACCCTGAAAGGTTAATTCAGCCAGACCAGCACCGACGGTCACAATGGTCGCGTGATAATCCCCCGCGGCAAGTGAAACGGTCTGTCCACTGCAATGCATATTGGTGATTTGCATACAAGGCACCTTCATTAACGACATACAGCCTCCCGTGAGGCTG‑ATGTCAAACCCaaaa                                                                                                                                                                                                                                                                                  >  2:600476/1‑284 (MQ=255)                             gCCCCTGATAGCGGTAACAGCCGTTGGCGATCCGGTTCGGCCACGGAATAAGCACTTTACCCAGATGCGCCAGCGGCATCTCTTCGGGTTTGTGCGGGATGACTAAATGACAACCCTGAAAGGTTAATTCAGCCAGACCAGCACCGACGGTCACAATGGTCGCGTGATAATCCCCCGCGGAAAGTGAAACGGCCTGTCCACTGCAATGCATATTGGTGATTTGCATACAAGGCACCTTAATTAACGACATACAGCCTCCCGTGAGGCTG‑ATGTCAAACCCAAACACTGCACCGGACgg                                                                                                                                                                                                                                                                    >  1:971836/1‑298 (MQ=255)                                ccTGATAGCGGTAACAGCCGTTGGCGATCCGGTTCGGCCACGGAATAAGCACTTTACCCAGATGCGCCAGCGGCATCTCTTCGGGTTTGTGCGGGATGACTAAATGACAACCCTGAAAGGTTAATTCAGCCAGACCAGCACCGACGGTCACAATGGTCGCGTGATAATCCCCCGCGGCAAGTGAAACGGTCTGTCCACTGCAATGCATATTGGTGATTTGCATACAAGGCACCTTAATTAACGACATACAGCCTCACGGGAGGCTGTATGTcgttaattaaggtg                                                                                                                                                                                                                                                                               <  2:212301/285‑14 (MQ=255)                                       gcggtgaaacgcGTTGGGGGGCCGGTTCGGCCACGGAATAAGCACTTTACCCAGATGCGCCAGCGGCATCTCTTCGGGTTTGTGCGGGATGACTAAATGACAACCCTGAAAGGTTAATTCAGCCAGACCAGCACCGACGGTCACAATGGTCGCGTGATAATCCCCCGCGGCAAGTGAAACGGTCTGTCCACTGCAATGCATATTGGTGATTTGCATACAAGGCACCTTAATTAACGACATACAGCCTCCCGTGAGGCTG‑ATGTCAAACCCAAAAACTGCACCGGa                                                                                                                                                                                                                                                                       <  2:467514/274‑1 (MQ=255)                                        cGGTAACAGCGGTTGGCGATCCGGTTCGGCCACGGAATAAGCACTTTACCCAGATGCGCCAGCGGCATCTCTTCGGGTTTGTGCGGGATGACTAAATGACAACCCTGAAAGGTTAATTCAGCCAGACCAGCACCGACGGTCACAATGGTCGCGTGATAATCCCCCGCGGCAAGTGAAACGGTCTGTCCACTGCAATGCATATTGGTGATTTGCATACAAGGCACCTTAATTAACGACATACAGCCTCCCGTGAGGCTG‑ATGTCAAACCCAAAAACTGCACCGGAc                                                                                                                                                                                                                                                                      <  2:155636/285‑1 (MQ=255)                                              cAGCCGTTGGCGATCCGGTTCGGCCACGGAATAAGCACTTTACCCAGATGCGCCAGCGGCATCTCTTCGGGTTTGTGCGGGATGACTAAATGACAACCCTGAAAGGTTAATTCAGCCCGACCAGCACCGACGGTCACAATGGTCGCGTGATAATCCCCCGCGGCACGTGAAACGGTCTGTCCACTGCAATGCATATTGGTGATTTGCATCCAACGCACCTTAATTAACGACATACAGCCTCCCGTACGGCTG‑ATGTCAAACCCAAACACTGCACCGGACGGTccc                                                                                                                                                                                                                                                                >  2:489023/1‑285 (MQ=255)                                                   gTTGGCGATCCGGTTCGGCCACGGAATAAGCACTTTACCCAGATGCGCCAGCGGCATCTCTTCGGGTTTGTGCGGGTTGACTAAATGACAACCCTGAAAGGTTAATTCAGCCAGACCAGCACCGACGGTCACAATGGTCGCGTGATAATCCCCCGCGGCAAGTGAAACGGTCTGTCCACTGCAATGCATATTGGTGATTTGCATACAAGGCACCTTAATTAACGACATACAGCCTCCCGTGAGGCTG‑ATGTCAAACCCAAAAACTGCACCGGACGGTCCCCTCGCCCCTTTGGGGagag                                                                                                                                                                                                                                             <  1:80089/299‑1 (MQ=255)                                                   gTTGGCGATCCGGTTCGGCCACGGAATAAGCACTTTACCCAGATGCGCCAGCGGCATCTCTTCGGGTTTGTGCGGGATGACTAAATGACAACCCTGAAAGGTTAATTCAGCCAGACCAGCACCGACGGTCACAATGGTCGCGTGATAATCCCCCGCGGCAAGTGAAACGGTCTGTCCACTGCAATGCATATTGGTGATTTGCATACAAGGCACCTTAATTAACGACATACAGCCTCCCGTGAGGCTG‑ATGTCAAACCCAAAAACTGCACCGGACGGTCCCCTCGc                                                                                                                                                                                                                                                           <  2:640491/285‑1 (MQ=255)                                                           tCCGGTTCGGCCACGGAATAAGCACTTTACCCAGATGCGCCAGCGGCATCTCTTCGGGTTTGTGCGGGATGACTAAATGACAACCCTGAAAGGTTAATTCAGCCAGACCAGCACCGACGGTCACAATGGTCGCGTGATAATCCCCCGCGGCAAGTGAAACGGTCTGTCCACTGCAATGCATATTGGTGATTTGCATACAAGGCACCTTAATTAACGACATACAGCCTCCCGTGAGGCTG‑ATGTCAAACCCAAAAACTGCACCGGACGGTCCCCTCGCCCCTTTGGGGAGAGGGTTAgg                                                                                                                                                                                                                                      <  1:413977/298‑1 (MQ=255)                                                            ccGGTTCGGCCACGGAATAAGCACTTTACCCAGATGCGCCAGCGGCATCTCTTCGGGTTTGTGCGGGATGACTAAATGACAACCCTGAAAGGTTAATTCAGCCAGACCAGCACCGACGGTCACAATGGTCGCGTGATAATCCCCCGCGGCAAGTGAAACGGTCTGTCCACTGCAATGCATATTGGTGATTTGCATACAAGGCACCTTAATTAACGACATACAGCCTCCCGTGAGGCTG‑ATGTCAAACCCAAAAACTGCACCGGACGGTCCCCTCGCCCCTTTGGGGAGAGGGTTAGGga                                                                                                                                                                                                                                    <  1:586312/299‑1 (MQ=255)                                                                ttCGGCCACGGAATAAGCACTTTACCCAGATGCGCCAGCGGCATCTCTTCGGGTTTGTGCGGGATGACTAAATGACAACCCTGAAAGGTTAATTCAGCCAGACCAGCACCGACGGTCACAATGGTCGCGTGATAATCCCCCGCGGCAAGTGAAACGGTCTGTCCACTGCAATGCATATTGGTGATTTGCATACAAGGCACCTTAATTAACGACATACAGCCTCCCGTGAGGCTG‑ATGTCAAACCCAAAAACTGCACCGGACGGTCCCCTCGCCCCTTTggg                                                                                                                                                                                                                                                  <  2:545898/281‑1 (MQ=255)                                                                 tCGGCCACGGAATAAGCACTTTACCCAGATGCGCCAGCGGCATCTCTTCGGGTTTGTGCGGGATGACTAAATGACAACCCTGAAAGGTTAATTCAGCCAGACCAGCACCGACGGTCACAATGGTCGCGTGATAATCCCCCGCGGCAAGTGAAACGGTCTGTCCACTGCAATGCATATTGGTGATTTGCATACAAGGCACCTTAATTAACGACATACAGCCTCCCGTGAGGCTG‑ATGTCAAACCCAAAAACTGCACCGGACGGTCCCCTCGCCCCTTTGGGGAGAGGGTTAGGGAGAgg                                                                                                                                                                                                                                <  1:965774/298‑1 (MQ=255)                                                                    gCCACGGAATAAGCACTTTACCCAGATGCGCCAGCGGCATCTCTTCGGGTTTGTGCGGGATGACTAAATGACAACCCTGAAAGGTTAATTCAGCCAGACCAGCACCGACGGTCACAATGGTCGCGTGATAATCCCCCGCGGCAAGTGAAACGGTCTGTCCACTGCAATGCATATTGGTGATTTGCATACAAGGCACCTTAATTAACGACATACAGCCTCCCGTGAGGCTG‑ATGTCCAACCCAAAAACTGCACCGGACGGTCCCCTCGCCCCTTTGGGGAGAGGGt                                                                                                                                                                                                                                          >  2:506941/1‑285 (MQ=255)                                                                                gCACTTTACCCAGATGCGCCAGCGGCATCTCTTCGGGTTTGTGCGGGATGACTAAATGACAACCCTGAAAGGTTAATTCAGCCAGACCAGCACCGACGGTCACAATGGTCGCGTGATAATCCCCCGCGGCAAGTGAAACGGTCTGTCCACTGCAATGCATATTGGTGATTTGCATACAAGGCACCTTAATTAACGACATACAGCCTCCCGTGAGGCTG‑ATGTCAAACCCAAAAACTGCACCGGCCGGTCCCCTCGCCCCTTTGGGGAGAGGGTTAGGGAGAgggg                                                                                                                                                                                                                              >  2:742997/1‑285 (MQ=255)                                                                                                        gCATCTCTTCGGGTTTGTGCGGGATGACTAAATGACAACCCTGAAAGGTTAATTCAGCCAGACCAGCACCGACGGTCACAATGGTCGCGTGATAATCCCCCGCGGCAAGTGAAACGGTCTGTCCACTGCAATGCATATTGGTGATTTGCATACAAGGCACCTTAATTAACGACATACAGCCTCCCGTGAGGCTG‑ATGTCAAACCCAAAAACTGCACCGGACGGTCCCCTCGCCCCTTTGGGGAGAGGGTTAGGGAGAGGCGAAAACCCGCACCCTATTTTCCGATCAGTTATTTCGCAt                                                                                                                                                                                        >  1:325271/1‑299 (MQ=255)                                                                                                                         tgCGGGATGACTAAATGACAACCCTGAAAGGTTAATTCAGCCAGACCAGCACCGACGGTCACAATGGTCGCGTGATAATCCCCCGCGGCAAGTGAAACGGTCTGTCCACTGCAATGCATATTGGTGATTTGCATACAAGGCACCTTAATTAACGACATACAGCCTCCCGTGAGGCTG‑ATGTCAAACCCAAAAACTGCACCGGACGGTCCCCTCGCCCCTTTGGGGAGAGGGTTAGGGAGAGGGGAAAACCCGCACCCTCTTTTCAGATCAGTTCTTTCGCATTa                                                                                                                                                                                      >  2:551393/1‑284 (MQ=255)                                                                                                                              gATGACTAAATGACAACCCTGAAAGGTTAATTCAGCCAGACCAGCACCGACGGTCACAATGGTCGCGTGATAATCCCCCGCGGCAAGTGAAACGGTCTGTCCACTGCAATGCATATTGGTGATTTGCATACAAGGCACCTTAATTAACGACATACAGCCTCCCGTGAGGCTG‑ATGTCAAACCCAAAAACTGCACCGGACGGTCCCCTCGCCCCTTTGGGGAGAGGGTTAGGGAGAGGGGAAAACCCGCACCCTATTTTCAGATCAGTTATTTCGCATTAATATCCAGCAGACCCGCTGc                                                                                                                                                                  >  1:493600/1‑299 (MQ=255)                                                                                                                                 gACTAAATGACAACCCTGAAAGGTTAATTCAGCCAGACCAGCACCGACGGTCACAATGGTCGCGTGATAATCCCCCGCGGCAAGTGAAACGGTCTGTCCACTGCAATGCATATTGGTGATTTGCATACAAGGCACCTTAATTAACGACATACAGCCTCCCGTGAGGCTG‑ATGTCAAACCCAAAAACTGCACCGGACGGTCCCCTCGCCCCTTTGGGGAGAGGGTTAGGGAGAGGGGAACACCCGCACCCTATTTTCAGATCAGTTATTTCGCATTAATATCcag                                                                                                                                                                              >  2:138567/1‑284 (MQ=255)                                                                                                                                              cccTGAAAGGTTAATTCAGCCAGACCAGCACCGACGGACACAATGGTCGCGTGATAATCCCCCGCGGCAAGTGAAACGGTCTGTCCACTGCAATGCATATTGGTGATTTGCATACAAGGCACCTTAATTAACGACATACAGCCTCCCGTGAGGCTG‑AGGTCAAACCCAAAAACTGCACCGGACGGTCCCCTCGCCCCGTTGGGGAAAGGGTTAGGGAGAGGGGAAAACCCGCACCCTATTGTAAGAACAGTTATTTCGCATTAATCTCCAACACACCCGcgcca                                                                                                                                                                 >  2:451487/1‑280 (MQ=255)                                                                                                                                                  gAAAGGTTAATTCAGCCAGACCAGCACCGACGGTCACAATGGTCGCGTGATAATCCCCCGCGGCAAGTGAAACGGTCTGTCCACTGCAATGCATATTGGTGATTTGCATACAAGGCACCTTAATTAACGACATACAGCCTCCCGTGAGGCTG‑ATGTCAAACCCAAAAACTGCACCGGACGGTCCCCTCGCCCCTTTGGGGAGAGGGTTAGGGAGAGGGGAAAACCCGCACCCTATTTTCAGATCAGTTATTTCGCATTAATATCCAGCAGACCCGCTGCAACCGCTGGAGCCACGCCgg                                                                                                                                              >  1:224392/1‑299 (MQ=255)                                                                                                                                                       gTTAATTCAGCCAGACCAGCACCGACGGTCACAATGGTCGCGTGATAATCCCCCGCGGCAAGTGAAACGGTCTGTCCACTGCAATGCATATTGGTGATTTGCATATAAGGCACCTTAATTAACGACATACAGCCTCCCGTGAGGCTG‑ATGTCAAACCCAAAAACTGCACCGGACGGTCCCCTCGCCCCTTTGGGGAGAGGGTTAGGGAGAGGGGAAAACCCGCACCCTATTTTCAGATCAGTTATTTCGCATTAATATCCAGCAGACCCGCTGCAACCGCTGGAGCCATGCCGGGCGCt                                                                                                                                         >  1:325151/1‑299 (MQ=255)                                                                                                                                                            ttCAGCCAGACCAGCACCGACGGTCACAATGGTCGCGTGATAATCCCCCGCGGCAAGTGAAACGGTCTGTCCACTGCAATGCATATTGGTGATTTGCATACAAGGCACCTTAATTAACGACATACAGCCTCCCGTGAGGCTG‑ATGTCAAACCCAAAAACTGCACCGGACGGTCCCCTCGCCCCTTTGGGGAGAGGGTTAGGGAGAGGGGAAAACCCGCACCCTATTTTCAGATCAGTTATTTCGCATTAATATCCAGCAGACCCGCTGCAACCGCTGGAGCCATGCCGGGCGCTAACgg                                                                                                                                    >  1:500501/1‑299 (MQ=255)                                                                                                                                                               aGCCAGACCAGCACCGACGGTCACAATGGTCGCGTGATAATCCCCCGCGGCAAGTGAAACGGTCTGTCCACTGCAATGCATATTGGTGATTTGCATACAAGGCACCTTAATTAACGACATACAGCCTCCCGTGAGGCTG‑ATGTCAAACCCAAAAACTGCACCGGACGGTCCCCTCGCCCCTTTGGGGAGAGGGTTAGGGAGAGGGGGAAACCCGCACCCTATTTTCAGATCAGTTATTTCGCATTAATATCCAGCAGACCCGCTGCAACCGCTGGAGCCATGCCGGGCGCTAACGGGAt                                                                                                                                 >  1:949319/1‑299 (MQ=255)                                                                                                                                                                 ccAGACCAGCACCGACGGTCACAATGGTCGCGTGATAATCCCCCGCGGCAAGTGAAACGGTCTGTCCACTGCAATGCATATTGGTGATTTGCATACAAGGCACCTTAATTAACGACATACAGCCTCCCGTGAGGCTG‑ATGTCAAACCCAAAAACTGCACCGGACGGTCCCCTCGCCCCTTTGGGGAGAGGGTTAGGGAGAGGGGAAAACCCGCACCCTCTTTTCAGATCAGTTATTTCGCATTAATATCCAGCAGACCCGCTGCAACCGCTGGATCCATGCCGGGCGCTAACGGGATa                                                                                                                                >  1:453241/1‑298 (MQ=255)                                                                                                                                                                       cAGCACCGACGGTCACAATGGTCGCGTGATAATCCCCCGCGGCAAGTGAAACGGTCTGTCCACTGCAATGCATATTGGTGATTTGCATACAAGGCACCTTAATTAACGACATACAGCCTCCCGTGAGGCTG‑ATGTCAAACCCAAAAACTGCACCGGACGGTCCCCTCGCCCCTTTGGGGAGAGGGTTAGGGAGAGGGGAAAACCCGCACCCTATTTTCAGATCAGTTATTTCGCATTAATATCCAGCAGACCCGCTGCAACCGCTGGAGCCATGCCGGGCGCTAACGGGATACGCGGGa                                                                                                                         >  1:927944/1‑299 (MQ=255)                                                                                                                                                                       cAGCACCGACGGTCACAATGGTCGCGTGATAATCCCCCGCGGCAAGTGAAACGGTCTGTCCACTGCAATGCATATTGGTGATTTGCATACAAGGCACCTTAATTAACGACATACAGCCTCCCGTGAGGCTG‑ATGTCAAACCCAAAAACTGCACCGGACGGTCCCCTCGCCCCTTTGGGGAGAGGGTTAGGGAGAGGGGAAAACCCGCACCCTATTTTCAGATCAGTTATTTCGCATTAATATCCAGCAGACCCGCTGCAACCGCTGGAGCCATGCCGGGCGCTAACGGGATACGCGGGa                                                                                                                         >  1:569013/1‑299 (MQ=255)                                                                                                                                                                           aCCGACGGTCACAAGGGTCGCGTGATAATCCCCCGCGGCAAGTGAAACGGTCTGTCCACTGCAATGCATATTGGTGATTTGCATACAAGGCACCTTAATTAACGACATACAGCCTCCCGTGAGGCTG‑ATGTCAAACCCAAAACCGGCACCGGACGGTCCCCTCGCCCCTTTGGGGAGAGGGTTAGGGAGAGGGGAAAACCCGCACCCTATTTTCAGATCAGTTATTTCGCATTAATATCCAGCAGACCCGCTGCAACCGCTGGAGCCATGCCGGGCGCTAACGGGATACGCGGGATCa                                                                                                                      <  1:188680/298‑1 (MQ=255)                                                                                                                                                                            ccGACGGTCACAATGGTCGCGTGATAATCCCCCGCGGCAAGTGAAACGGTCTGTCCACTGCAATGCATATTGGTGATTTGCATACAAGGCACCTTAATTAACGACATACAGCCTCCCGTGAGGCTG‑ATGTCAAACCCAAAAACTGCACCGGACGGTCCCCTCGCCCCTTTGGGGAGAGGGTTAGGGAGAGGGGAAAACCCGCACCCTATTTTCAGATCAGTTATTTCGCATTAATATCCAGCAGACCCGCTGCAACCGCTGGAGCCATGCCGGGCGCTAACGGGATACGCGGGATCAcc                                                                                                                    >  1:597576/1‑299 (MQ=255)                                                                                                                                                                             cGACGGTCACAATGGTCGCGTGATAATCCCCCGCGGCAAGTGAAACGGTCTGTCCACTGCAATGCATATTGGTGATTTGCATACAAGGCACCTTAATTAACGACATACAGCCTCCCGTGAGGCTG‑ATGTCAAACCCAAAAACTGCACCGGACGGTCCCCTCGCCCCTTTGGGGAGAGGGTTAGGGAGAGGGGAAAACCCGCACCCTATTTTCAGATCAGTTATTTCGCATTAATATCCAGCAGACCCGCTGCAACCGCTGGAGCCATGCCGGGCGCTAACGGGATACGCGGGATCACCa                                                                                                                   >  1:637245/1‑299 (MQ=255)                                                                                                                                                                               aCGGTCACAATGGTCGCGTGATAATCCCCCGCGGCAAGTGAAACGGTCTGTCCACTGCAATGCATATTGGTGATTTGCATACAAGGCACCTTAATTAACGACATACAGCCTCCCGTGAGGCTG‑ATGTCAAACCCAAAAACTGCACCGGACGGTCCCCTCGCCCCTTTGGGGAGAGGGTTAGGGAGAGGGGAAAACCCGCACCCTCTTTTCAGATCAGTTATTTCGCATTAATATCAAGCAGACCCGCTGCACCCGCTGGCGCCATGCCGAGAGCTAACGGGGTaa                                                                                                                               >  2:468951/1‑284 (MQ=255)                                                                                                                                                                                    cacaATGGTCGGGTGAAAAGCCCCCGCGGAAAGTTAAACGGTCTGGCCAATGCAAGGCATATTGGTGATTTGCATAAAAGGCACCTTAATTAACGACATCCAGCCCCCCGTGAGGCTG‑ATGTCAAACCCAAAAACTGCACCGGACGGTCCCCTCGCCCCTTTGGGGAGAGGGTTAGGGAGAGGGGAAAACCCGCACCCTATTTTCAGATCAGTTATTTCGCATTAATATCCAGCAGACCCGCTGCAACCGCTGGAGCCACGCCGGGCGCTAACGGGATACGCgg                                                                                                                           <  2:192771/284‑1 (MQ=255)                                                                                                                                                                                                   aaaaccccccGCGGCAAGTGAAACGTTCTGTCCACGGCAATGCATATTGGTGATTGGCATACAAGGCCCCTTAATTAACGACAACCCGCCTCCCGTGAGGCTG‑ATGTCAAACCCAAAAACTGCACCGGACGGTCCCCTCGCCCCTTTGGGGAGAGGGTTAGGGAGAGGGGAAAACCCGCACCCTATTTTCAGATCAGTTATTTCGCATTAATATCCAGCAGACCCGCTGCAACCGCTGGAGCCATGCCGGGCGCTAACGGGATACGCGGGATCACCAGGCAATGCAGCAGGTGATAAat                                                                                             <  1:283662/294‑1 (MQ=255)                                                                                                                                                                                                    tAATCCCCCGCGGCAAGTGAAACGGTCTGTCCACTGCAATGCATATTGGTGATTTGCATACAAGGCACCTTAATTAACGACATACAGCCTCCCGTGAGGCTG‑ATGTCAAACCCAAAAACTGCACCGGACGGTCCCCTCGCCCCTTTGGGGAGAGGGTTAGGGAGAGGGGAAAACCCGCACCCTATTTTCAGATCAGTTATTTCGCATTAATATCCAGCAGACCCGCTGCAACCGCTGGAGCCATGCCGGGCGCTAACGGGATACGCGGGATAACCAGGCAATGCAGCAGGTGATAAata                                                                                            >  1:400536/1‑299 (MQ=255)                                                                                                                                                                                                      agcccccGCGGCAAGTGAAACGGTCTGTCCACGGCAATGCATATTGGTGATTTGAATACAAGGCACCTTAATTAACGACATACAGCCTCCCGTGAGGCGG‑ATGTCAAACCCAAAAACGGCACCGGACGGTCCCCTCGCCCCTTTGGGGAGAGGGTTAGGGAGAGGGGAAAACCCGCACCCTATTTTCAGATCAGTTATTTCGCATTAATATCCAGCAGACCCGCTGCAACCGCTGGAGCCATGCCGGGCGCTAACGGGATACGCGGGATCACCAGGCAATGCAGCAGGTGATAAATATc                                                                                          <  1:638920/297‑1 (MQ=255)                                                                                                                                                                                                                                    aCTGCAATGCATATTGGTGATTTGCATACAAGGCACCTTAATTAACGACATACAGCCTCCCGTGAGGCTG‑ATGTCAAACCCAAAAACTGCACCGGACGGTCCCCTCGCCCCTTTGGGGAGAGGGTTAGGGAGAGGGGAAAACCCGCACCCTATTTTCAGATCAGTTATTTCGCATTAATATCCAGCAGACCCGCTGCAACCGCTGGAGCCATGCCGGGCGCTAACGGGATACGCGGGATCACCAGGCAATGCAGCAGGTGATAAATATCCTGTTTGCGGTCCCAGACTTTGGTGGTGa                                                             >  1:647822/1‑298 (MQ=255)                                                                                                                                                                                                                                        cAATGCATATTGGTAATTTGCATACAAGGCACCTTAATTAACGACAAACCGCCTCCCTGGAGGCTG‑ATTTCAAACCCAAAACCTGCACCGTACGGTCCCCTCGCCCCTTTGGGGAGGGGGTTAGGGAGAGGGGAAAACCCGCACCCTATTTTCAGATCAGTTATTTCGCATTAATATCCAGCAGACCCGCTGCAACCGCTGGAGCCATGCCGGGCGCTAACGGGATACGCGGGATCACCAGGCAATGCAGCAGGTGATAAATATCCTGTTTGCCGTCCCAGACTTTGGTGGTGACCTTa                                                        <  1:864981/299‑1 (MQ=255)                                                                                                                                                                                                                                           tGCATATTGGTGATTTGCATACAAGGCACCTTAATTAACGACATACAGCCTCCCGTGAGGCTG‑ATGTCAAACCCAAAAACTGCACCGGACGGTCCCCTCGCCCCTTTGGGGAGAGGGTTAGGGAGAGGGGAAAACCCGCACCCTATTTTCCGTTCAGTTATTTCGCATTAATATCCAGCAGACCCGCTGCAACCGCTGGAGCCATGCCGGGCGCTAACGGGATACGCGGGATCACCAGGCAATGCAGCAGGTGATAAATATCCTGTTTGCCGTCCCAGACTTTGGTGGTGACCTTATTg                                                     >  1:247405/1‑299 (MQ=255)                                                                                                                                                                                                                                                 ttGGTGATTTGCATACAAGGCACCTTAATTAACGACATACAGCCTCCCGTGAGGCTG‑ATGTCAAACCCAAAAACTGCACCGGACGGTCCCCTCGCCCCTTTGGGGAGAGGGTTAGGGAGAGGGGAAAACCCGCACCCTATTTTCAGATCAGTTATTTCGCATTAATATCCAGCAGACCCGCTGCAACCGCTGGAGCCATGCCGGGCGCTAACGGGATACGCGGGATCACCAGGCAATGCAGCAGGTGATAAATATCCTGTTTGCCGTCCCAGACTTTGgtggtg                                                              >  2:805571/1‑284 (MQ=255)                                                                                                                                                                                                                                                               aCAAGGCACCTTAATTAACGACATACAGCCTCCCGTGAGGCTG‑ATGTCAAACCCAAAAACTGCACCGGACGGTCCCCTCGCCCCTTTGGGGAGAGGGTTAGGGAGAGGGGAAAACCCGCACCCTATTTTCAGATCAGTTATTTCGCATTAATATCCAGCAGACCCGCTGCAACCGCTGGAGCCATGCCGGGCGCTAACGGGATACGCGGGATCCCCAGGCAATGCAGCAGGTGATAAATATCCTGTTTGCCGTCCCAGACTTTGGGGGTGACCTTATTGTCCGc                                                >  2:645899/1‑284 (MQ=255)                                                                                                                                                                                                                                                                   ggtaCCTTAATTAACGCCTTACAGCCTCTCGTGAGCTGG‑ATGTCAAACCCAAAACCTGCCCCGGACGTTCCCCCCCCCCCTTGGGGGAGAGGGTTAGGGGGAGGGGAAAACCCGCACCCTATTTTCAGATCAGTTATTTCGCATTAATATCCAGCAGACCCGCTGCAACCGCTGGAGCCATGCCGGGCGCTAACGGGATACGCGGGATCACCAGGCAATGCAGCAGGTGATAAATATCCTGTTTGCCGTCCCAGACTTTGGTGGTGACCTTATTGTCCGCATcc                                            <  2:346888/281‑1 (MQ=255)                                                                                                                                                                                                                                                                   ggCACCTTAATTAACGACATACAGCCTCCCGTGAGGCTG‑ATGTCAAACCCAAAAACTGCACCGGACGGTCCCCTCGCCCCTTTGGGGAGAGGGTTAGGGAGAGGGGAAAACCCGCACCCTATTTTCAGATCAGTTATTTCGCATTAATATCCAGCAGACCCGCTGCAACCGCTGGAGCCATGCCGGGCGCTAACGGGATACGCGGGATCACCAGGCAATGCAGCAGGTGATAAATATCCTGTTTGCAGTCACAGACTTTGGTGGTGACATTATTGTCCGCATCCAGCTCCTGCCACCAg                             >  1:307430/1‑299 (MQ=255)                                                                                                                                                                                                                                                                             ttaaCCACATCCCGCCTCCCGTGGGGCTG‑ATGTCAAACCCAAAAACGGCCCCGGACGGTCCCCTCGCCCCTTTGGGAGGAGGGTTAGGGAGAGGGGAAAACCCGCACCCTATTTTCAGATCAGTTATTTCGCATTAATATCCAGCAGACCCGCTGCAACCGCTGGAGCCATGCCGGGCGCTAACGGGATACGCGGGATCACCAGGCAATGCAGCAGGTGATAAATATCCTGTTTGCCGTCCCAGACTTTGGTGGTGACCTTATTGTCCGCATCCAGCTCCTGCCACCAGGAACCAtttt                   <  1:737993/299‑1 (MQ=255)                                                                                                                                                                                                                                                                                tttgacatcAGCCTCACGGGAGGCTGTATGTCAAACCCAAAAACTGCCCCGGACGGTCCCCTCGCCCCTTTGGGGAGAGGGTTAGGGAGAGGGGAAAACCCGCACCCTATTTTCAGATCAGTTATTTCGCATTAATATCCAGCAGACCCGCTGCCACCCCTGGAGCCCTGCCGGGCGCTAACGGGATACGCGGGATCACCAGGCAATGCAGCAGGTGATAAATATCCTGTTTGCCGTCCCAGACTTTGGTGGTGACCTTATTGTCCGCATCCAACTCCTGccacc                               >  2:541614/9‑285 (MQ=255)                                                                                                                                                                                                                                                                                   acacgcAGCCGCCCTTGAGGCGT‑ATGTCAAACCCAAAACCTGCACCGGACGGCCCCCTCGCCCCTTTGGGGAGAGGGTTAGGTAGAGGGGAAAACCCTCACCCTATTTTCAGATCAGTTATTTCGCATTAATATCCAGCAGACCCGCTGCAACCGCTGGAGCCATGCCGGGCGCTAACGGGATACGCGGGATCACCAGGCAATGCAGCAGGTGATAAATATCCTGTTTGCCGTCCCAGACTTTGGTGGTGACCTTATTGTCCGCATCCAGCTCCTGCCACCAGGa                           <  2:429605/280‑1 (MQ=255)                                                                                                                                                                                                                                                                                    caaccacccTCCCGGGAGGCTG‑ATGTAAAACCCAAAAACTGCCACGGACGGTCCCCCCGCCCCTTTGGGGAGAGGGTTAGGGAGAGGGGAAAACCCGCACCCTATTTTCAGATCAGTTATTTCGCATTAATATCCAGCAGACCCGCTGCAACCGCTGGAGCCATGCCGGGCGCTAACGGGATACGCGGGATCACCAGGCAATGCAGCAGGTGATAAATATCCTGTTTGCCGTCCCAGACTTTGGTGGTGACCTTATTGTCCGCATCCAGCTCCTGCCACCAGGAACCATTTTCATAGt             <  1:799485/291‑1 (MQ=255)                                                                                                                                                                                                                                                                                       aCAGCCTCCCGTGAGGCTG‑ATGTCAAACCCAAAAACTGCACCGGACGGTCCCCTCGCCCCTTTGGGGAGAGGGTTAGGGAGAGGGGAAAACCCGCACCCTATTTTCAGATCAGTTATTTCGCATTAATATCCAGCAGACCCGCTGCAACCGCTGGAGCCATGCCGGGCGCTAACGGGATACGCGGGATCACCAGGCAATGCAGCAGGTGATAAATATCCTGTTTGCCGTCCCAGACTTTGGTGGTGACCTTATTGTCCGCATCCAGCTCCTGCCACCAGGAAcc                        >  2:386059/1‑284 (MQ=255)                                                                                                                                                                                                                                                                                       aCAGCCTCCCGTGAGGCTG‑ATGTCAAACCCAAAAACTGCACCGGACGGTCCCCTCGCCCCTTTGGGGAGAGGGTTAGGGAGAGGGGAAAACCCGCACCCTATTTTCAGATCAGTTATTTCGCATTAATATCCAGCAGACCCGCTGCAACCGCTGGAGCCATGCCGGGCGCTAACGGGATACGCGGGATCACCAGGCAATGCAGCAGGTGATAAATATCCTGTTTGCCGTCCCAGACTTTGGTGGTGACCTTATTGTCCGCATCCAGCTCCTGCCACCAGGAAcc                        >  2:480192/1‑284 (MQ=255)                                                                                                                                                                                                                                                                                            cTCCCGTGAGGCTG‑ATGTCAAACCCAAAAACTGCACCGGACGGTCCCCTCGCCCCTTTGGGGAGAGGGTTAGGGAGAGGGGAAAACCCGCACCCTATTTTCAGATCAGTTATTTCGCATTAATATCCAGCAGACCCGCTGCAACCGCTGGAGCCATGCCGGGCGCTAACGGGATACGCGGGATCACCAGGCAATGCAGCAGGTGATAAATATCCTGTTTGACGTCCCAGACTTTGGTGGTGACCTTATTGTCCGCATCCAGCTCCTGCAACCAGGAACCATTTTCATAGGACATCAGGt    >  1:806776/1‑299 (MQ=255)                                                                                                                                                                                                                                                                                               ccGTGAGGCTG‑ATGTCAAACCCAAAAACTGCACCGGACGGTCCCCTCGCCCCTTTGGGGAGAGGGTTAGGGAGAGGGGAAAACCCGCACCCTATTTTCAGATCAGTTATTTCGCATTAATATCCAGCAGACCCGCTGCAACCGCTGGAGCCATGCCGGGCGCTAACGGGATACGCGGGATCACCAGGCAATGCAGCAGGTGATAAATATCCTGTTTGCCGTCCCAGACTTTGGTGGTGACCTTATTGTCCGCATCCAGCTCCTGCCACCAGGAACCATTTTCATAGtagatcggaaga  >  1:123293/1‑287 (MQ=255)                                                                                                                                                                                                                                                                                                  |                                                                                                                                                                                                                                                                                                        TGTTCATTAATGGGTAACTGATATTCCTGCCCCTGATAGCGGTAACAGCCGTTGGCGATCCGGTTCGGCCACGGAATAAGCACTTTACCCAGATGCGCCAGCGGCATCTCTTCGGGTTTGTGCGGGATGACTAAATGACAACCCTGAAAGGTTAATTCAGCCAGACCAGCACCGACGGTCACAATGGTCGCGTGATAATCCCCCGCGGCAAGTGAAACGGTCTGTCCACTGCAATGCATATTGGTGATTTGCATACAAGGCACCTTAATTAACGACATACAGCCTCCCGTGAGGCTG‑ATGTCAAACCCAAAAACTGCACCGGACGGTCCCCTCGCCCCTTTGGGGAGAGGGTTAGGGAGAGGGGAAAACCCGCACCCTATTTTCAGATCAGTTATTTCGCATTAATATCCAGCAGACCCGCTGCAACCGCTGGAGCCATGCCGGGCGCTAACGGGATACGCGGGATCACCAGGCAATGCAGCAGGTGATAAATATCCTGTTTGCCGTCCCAGACTTTGGTGGTGACCTTATTGTCCGCATCCAGCTCCTGCCACCAGGAACCATTTTCATAGTCCATCAGGTAC  >  NC_000913/4070124‑4070707

Alignment Legend
Aligned base mismatch/match (shaded by quality score): ATCG/ATCG < 3 ≤ ATCG/ATCG < 9 ≤ ATCG/ATCG < 21 ≤ ATCG/ATCG < 34 ≤ ATCG/ATCG < 38 ≤ ATCG/ATCG
Unaligned base: atcg    Masked matching base: atcg    Alignment gap: ‑    Deleted base: ‑