Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	NC_000913	4,602,509	C→T	W117* (TGG→TAG)	yjjP ←	putative succinate exporter YjjP

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	4,602,509	0	C	T	95.7%	172.7 / NA	47	W117* (TGG→TAG)	yjjP	putative succinate exporter YjjP
Reads supporting (aligned to +/- strand):  ref base C (0/0);  major base T (22/23);  minor base G (0/2);  total (22/25)
Fisher's exact test for biased strand distribution p-value = 4.91e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

GCAATGGTGGGGGTATTGCTGAAAGTGGGGAGTTGCAAAAGCAATCCGGAAATGGTGGTGGCGGCGAAAGCGGTAAGGCAAAAGTTGATCTGTGGATGAAGATGACGTTGTGCCAGCAGCTGGCGGATATACATCGCGGTCGTACTGGCAAAGAAGGTGATGACGGCACCATCCCAGCCACCGTTATTGAGTTTACAGAAACAGGCGCAAGAAAGGCCAACCATTAAGGCTACCAGCCATCTTGGGTAACGTAATGGCTGAATTTGGCTAAATCGTTTCTCTACGCCTTTGTAATCCAGCAGATGATGCTCCGCAAGAATCACAATGTGCTGGACTTCAGTCACCACATGCATATTAATGCCGCGATCGTGATTTTTACGTGTCGATGTCAGGCATTGCCCATCTTTAATAGTTGTCAGCACTATGGCGTTCGAAGAGATAGAACTTTCGACGCTGTCCATTCCCAGTGCCCGACCCA  >  NC_000913/4602272‑4602749                                                                                                                                                                                                                                              |                                                                                                                                                                                                                                                 gCAATGGTGGGGGTATTGCTGAAAGTGGGGAGTTGCAAAAGCAATCCGGAAATGGTGGTGGCGGCGAAAGCGGTAAGGCAAAAGTTGATCTGTGGATGAAGATGACGTTGTGCCAGCAGCTGGCGGATATACATCGCGGTCGTACTGGCAAAGAAGGTGATGACGGCACCATCCCAGCCACCGTTATTGAGTTTACAGAAACAGGCGCAAGAAAGGCCAACCATTAAGGCTACCAGCGATCTTGGGTAACg                                                                                                                                                                                                                                     <  1:99750/251‑1 (MQ=255)    aTGGTGGGGGTATTGCTGAAAGTGGGGAGTTGCAAAAGCAATCCGGAAATGGTGGTGGCGGCGAAAGCGGTAAGGCAAAAGTTGATCTGTGGATGAAGATGACGTTGTGCCAGCAGCTGGCGGATATACATCGCGGTCGTACTGGCAAAGAAGGTGATGACGGCACCATCCCAGCCACCGTTATTGAGTTTACAGAAACAGGCGCAAGAAAGGCCAACCATTAAGGCTACCAGCTATCTTGGGTAACGTa                                                                                                                                                                                                                                   >  1:381144/1‑250 (MQ=255)           gggTATTGCTGAAAGTGGGGAGTTGCAACAGCAATCCGGAAATGGTGGTGGCGGCGAAAGCGGTAAGGCAAAAGTTGATCTGTGGATGAAGATGACGTTGTGCCAGCAGCTGGCGGATATACATCGCGGTCGTACTGGCAAAGAAGGTGATGACGGCACCATCCCAGCCACCGTTATTGAGTTTACAGAAACAGGCGCAAGAAAGGCCAACCATTAAGGCTACCAGCTATCTTGGGTAACGTAATGGCTg                                                                                                                                                                                                                            >  2:443623/1‑250 (MQ=255)                 tGCTGAAAGTGGGGAGTTGCAAAAGCAATCCGGAAATGGTGGTGGCGGCGAAAGCGGGAAGGCAAAAGTTGATCTGTGGATGAAGATGACGTTGTGCCAGCAGCTGGCGGATATACATCGCGGTCGTACTGGCAAAGAAGGTGATGACGGCACCATCCCAGCCACCGTTATTGAGTTTACAGAAACAGGCGCAAGAAAGGCCAACCATTAAGGCTACCAGCTATCTTGGGTAACGTAATGGCTGAATTTg                                                                                                                                                                                                                      <  2:422180/250‑1 (MQ=255)                        aGTGGGGGGTTGCAAAAGCAATCCGTAAATGGTGGTGGCGGCGAAAGCGGTAAGGGAAAAGTTGATCTGTGGATGAAGATGACGTTGTGCCAGCAGCGGGCGGATATACATCGCGGTCGTACTTGCAAAGAAGGTGATGACGGCACCATCCCAGCCACCGTTATTGAGTTGACAGAAACAGGCGCAAGAAAGGCCAACCATTAAGGCTACCAGCGATCTTGGGTAACTTAATGGCTGAAGTTGGCTAAATc                                                                                                                                                                                                              <  2:339124/251‑1 (MQ=255)                        aGTGGGGAGTTGCAAAAGCAATCCGGAAATGGTGGTGGCGGCGAAAGCGGTAAGGCAAAAGTTGATCTGTGGATGAAGATGACGTTGTGCCAGCAGCTGGCGGATATACATCGCGGTCGTACTGGCAAAGAAGGTGATGACGGCACCATCCCAGCCACCGTTATTGAGTTTACAGAAACAGGCGCAAGAAAGGCCAACCATTAAGGCTACCAGCTATCTTGGGTAACGTAATGGCTGAATTTGGCTAAATc                                                                                                                                                                                                              >  2:121635/1‑251 (MQ=255)                                  tGCAACAGCAATCCGGAAATGGCGGCGGCCGCGAAAGCGGTAAGGCAAAAGTTGATCTGTGGATGAAGATGACGTTGTGCCAGAAGCTGGCGGATATCCATCGCGGTCGTACTGGCAACGAAGGTGATGACGGCACCATCCCAGCCACCGTTATTGAGTTTCCAGAAACAGGCGCAAGAAAGGCCAACCATTAAGGCTACCAGCTATCTTGGGTAACGTAATGGCTGAATTTGGCTAAATCGTTTCTCTAc                                                                                                                                                                                                    >  2:192929/1‑251 (MQ=255)                                     aaaaGCAATCCGGAAATGGTGGTGGCGGCGAAAGCGGTAAGGCAAAAGTTGATCTGTGGATGAAGATGACGTTGTGCCAGCAGCTGGCGGATATACATCGCGGTCGTACTGGCAAAGAAGGTGATGACGGCACCATCCCAGCCACCGTTATTGAGTTTACAGAAACAGGCGCAAGAAAGGCCAACCATTAAGGCTACCAGCTATCTTGGGTAACGTAATGGCTGAATGTGGCTAAATCGTTTCTCTACGcc                                                                                                                                                                                                 <  2:230254/251‑1 (MQ=255)                                        aGCAATCCGGAAATGGTGGTGGCGGCGAAAGCGGTAAGGCAAAAGTTGATCTGTGGATGAAGATGACGTTGTGCCAGCAGCTGGCGGATATACATCGCGGTCGTACTGGCAAAGAAGGTGATGACGGCACCATCCCAGCCACCGTTATTGAGTTTACAGAAACAGGCGCAAGAAAGGCCAACCATTAAGGCTACCAGCTATCTTGGGTAACGTAATGGCTGAATTTGGCTAAATCGTTTCTCTACGCCttt                                                                                                                                                                                              <  1:382161/251‑1 (MQ=255)                                          cAATCCGGAAATGGTGGTGGCGGCGAAAGCAGTAAGGCAAAAGTTGATCTGGGGATGAAGATGACGTTGTGCCAGCAGCTGGCGGATATACATCGCGGTCGTACTGGCAAAGAAGGTGATGACGGCACCATCCCAGCCACCGTTATTGAGTTTACAGAAACAGGCGCAAGAAAGGCCAACCATTAAGGCTACCAGCTATCTTGGGTAACGTAATGGCTGAATTTGGCTAAATCTTTTCTCTACGCCTTTGt                                                                                                                                                                                            <  2:360082/251‑1 (MQ=255)                                                  aaaTGGTGGTGGCGGCGAAAGCGGTAAGGCAAAAGTTGATCTGTGGATGAAGATGACGTTGTGCCAGCAGCTGGCGGATATACATCGCGGTCGTACTGGCAAAGAAGGTGATGACGGCACCATCCCAGCCACCGTTATTGAGTTTACAGAAACAGGCGCAAGAAAGGCCAACCATTAAGGCTACCAGCTATCTTGGGTAACGTAATGGCTGAATTTGGCTAAATCGTTTCTCTACGCCTTTGTAATCcagc                                                                                                                                                                                    >  1:429744/1‑251 (MQ=255)                                                   aaTGGTGGTGGCGGCGAAAGCGGTAAGGCAAAAGTTGATCTGTGGATGAAGATGACGGTGTGCCAGCAGCTGGCGGATATACATCGCGGTCGTACTGGCAAAGAAGGTGATGACGGCACCAGCCCAGCCACCGTTATTGAGTTTACAGAAACATTCGCAAGAGAGGCCAACCATGAAGTCTTCCAGCTATCTTGGGTAACGTAATGGCTGAATTTGGCTAAATCGTTTCTCTACGCCTTTGTAATCcagca                                                                                                                                                                                   <  2:228093/251‑1 (MQ=255)                                                         ggtggCGGCGAAAGCGGTAAGGCAAAAGTTGATCTGTGGATGAAGATGACGTTGTGCCAGCAGCTGGCGGATATACATCGCGGGCGTACTGGCAAAGAAGGTGATGACGGCACCATCCCAGCCACCGTTATTGAGTTTACAGAAACAGGCGCAATAAAGGCCAACCATTAAGGCTACCAGCTATCTTGGGTAACGTAATGGCTGAATTTGGCTAAATCGTTTCTCTACGCCTTTGTAATCCAGCAgatgat                                                                                                                                                                             <  1:224898/251‑1 (MQ=255)                                                            ggcggcGAAAGCGGTAAGGCAAAAGTTGATCTGTGGATGAAGATGACGTTGTGCCAGCAGCTGGCGGATATACATCGCGGTCGTACTGGCAAAGAAGGTGATGACGGCACCATCCCAGCCACCGTTATTGAGTTTACAGAAACAGGCGCAAGAAAGGCCAACCATTAAGGCTACCAGCTATCTTGGGTAACGTAATGGCTGAATTTGGCTAAATCGTTTCTCTACGCCTTTGTAATCCAGCAGATGATGCt                                                                                                                                                                          >  2:78940/1‑251 (MQ=255)                                                               ggcgAAAGCGGTAAGGCAAAAGTTGATCTGTGGATGAAGATGACGTTGTGCCAGCAGCTGGCGGATATACATCGCGGTCGTACTGGCAAAGAAGGTGATGACGGCACCATCCCAGCCACCGTTATTGAGTTTACAGAAACAGGCGCAAGAAAGGCCAACCATTAAGGCTACCAGCTATCTTGGGTAACGTAATGGCTGAATTTGGCTAAATCGTTTCTCTACGCCTTTGTAATCCAGCAGATGATGCTCCg                                                                                                                                                                       <  1:32975/251‑1 (MQ=255)                                                                  gAAAGCGGTCAGGCAAAAGTTGATCTGTGGATGAAGATGACGTTGTGCCAGCAGCTGGCGGATATACATCGCGGTCGTACTGGCAAAGAAGGTGATGACGGCACCATCCCAGCCACCGTTATTGAGTTTACAGAAACAGGCGCAAGAAAGGCCAACCATTAAGGCTACCAGCTATCTTGGGTAACGTAATGGCTGAATTTGGCTAAATCGTTTCTCTACGCCTTTGTAATCCAGCAGATGATGCTCCGCa                                                                                                                                                                     >  1:95018/1‑250 (MQ=255)                                                                  gAAAGCGGTAAGGCAAAAGTTGATCTGTGGCTGAAGATGACGTTGTGCCAGCAGCTGGCGGATATACATCGCGGTCGTACTGGCAAAGAAGGTGATGACGGCACCATCCCAGCCACCGTTATTGAGTTTACAGAAACAGGCGCAAGAAAGGCCAACCATTAAGGCTACCAGCTATCTTGGGTAACGTAATGGCTGAATTTGGCTAAATCGTTTCTCTACGCCTTTGTAATCCAGCAGATGATGCTCCGCa                                                                                                                                                                     >  1:445113/1‑250 (MQ=255)                                                                  gAAAGCGGTAAGGCAAAAGTTGATCTGTGGATGAAGATGACGTTGTGCCAGCAGCTGGCGGATATACATCGCGGTCGTACTGGCAAAGAAGGTGATGACGGCACCATCCCAGCCACCGTTATTGAGTTTACAGAAACAGGCGCAAGAAAGGCCAACCATTAAGGCTACCAGCTATCTTGGGTAACGTAATGGCTGAATTTGGCTAAATCGTTTCTCTACGCCTTTGTAATCCAGCAGATGATGCTCCGCaa                                                                                                                                                                    <  1:142101/251‑1 (MQ=255)                                                                           aaGGCAAAAGTTGATCTGTGGATGAAGATGACGTTGTGCCAGCAGCTGGCGGATATACATCGCGGTCGTACTGGCAAGGAAGGTGATGACGGCACCATCCCAGCCACCGTTATTGAGTTTACAGAAACAGGCGCAAGAAAGGCCAACCATTAAGGCTACCAGCTATCTTGGGTAACGTAATGGCTGAATTTGGCTAAATCGTTTCTCTACGCCTTTGTAATCCAGCAGATGATGCTCCGCAAGAATcaca                                                                                                                                                            >  2:114218/1‑250 (MQ=255)                                                                            aGGCAAAAGTTGATCTGTGGATGAAGATGACGTTGTGCCAGCAGCTGGCGGATATACATCGCGGTCGTACTGGCAAAGAAGGTGATGACGGCACCATCCCAGCCACCGTTATTGAGTTTACAGAAACAGGCGCATGAAAGGCCAACCATTAAGGCTACCAGCTATCTTGGGTAACGTAATGGCTGAAGTTGGCTAAATCGTTTCTCTACGCCTTTGTAATCCAGCAGATGATGCTCCGCAAGAATCACAAt                                                                                                                                                          <  1:278216/251‑1 (MQ=255)                                                                              gCAAAAGTTGATCTGTGGATGAAGATGACGTTGTGCCAGCAGCTGGCGGATATACATCGCGGTCGTACTGGCAAAGAAGGTGATGACGGCACCATCCCAGCCACCGTTATTGAGTTTACAGAAACAGGCGCAAGAAAGGCCAACCATTAAGGCTACCAGCTATCTTGGGTAACGTAATGGCTGAATTTGGCTAAATCGTTTCTCTACGCCTTTGTAATCCAGCAGATGATGCTCCGCAAGAATCACAAtgt                                                                                                                                                        >  2:434145/1‑251 (MQ=255)                                                                                 aaaGTTGATCTGTGGATGAAGATGACGTTGTGCCAGCAGCTGGCGGATATACATCGCGGTCGTACTGGCAAAGAAGGTGATGACGGCACCATCCCAGCCACCGTTATTGAGTTTACAGAAACAGGCGCAAGAAAGGCCAACCATTAAGGCTACCAGCTATCTTGGGTAAc                                                                                                                                                                                                                                      >  1:249418/1‑170 (MQ=255)                                                                                 aaaGTTGATCTGTGGATGAAGATGACGTTGTGCCAGCAGCTGGCGGATATACATCGCGGTCGTACTGGCAAAGAAGGTGATGACGGCACCATCCCAGCCACCGTTATTGAGTTTACAGAAACAGGCGCAAGAAAGGCCAACCATTAAGGCTACCAGCTATCTTGGGTAAc                                                                                                                                                                                                                                      <  2:249418/170‑1 (MQ=255)                                                                                       gATCTGTGGATGAAGATGACGTTGTGCCAGCAGCTGGCGGATATACATCGCGGTCGTACTGGCAAAGAAGGTGATGACGGCACCATCCCAGCCACCGTTATTGAGTTTACAGAAACAGGCGCAAGAAAGGCCAACCATTAAGGCTACCAGCTATCTTGGGTAACGTAATGGCTGAATTTGGCTAAATCGTTTCTCTACGCCTTTGTAATCCAGCAGATGATGCTCCGCAAGAATCACAATGTGCTGGACtt                                                                                                                                               >  2:307312/1‑251 (MQ=255)                                                                                                  gAAGATGACGTTGTGCCAGCAGCTTGCGGATATACATCGCGGTCGTACTGGCAAAGAAGGTTATGACTGCACCATCCCAGCCACCGTTAGTGATTGTACAGAAACAGGCGCAAGAAAGGCCAACCATTAAGGCTACCAGCTATCGTGGGTAACTTAAGGGCTGAATTGGGCTAAATCGTTGCTCGACGCCTTTGTAATCCAGCAGATGGTGCTCCACAAGAATCACAATGTGCTGGACTTCAGTcaccac                                                                                                                                     <  1:209978/250‑1 (MQ=255)                                                                                                     gATGACGTTGTGCCAGCAGCTGGCGGATATACATCGCGGTCGTACTGGCAAAGAAGGTGATGACGGCACCATCCCAGCCACCGTTATTGAGTTCACAGAAACAGGCGCAAGAAAGGCCAACCATTAAGGCTACCAGCTATCTTGGGTAACGTAATGGCTGAATTTGGCTAAATCGTTTCTCTACGCCTTTGTAATCCAGCAGATGATGCTCCGCAAGAATCACAATGTGCTGGACTTCAGTCACCACATGc                                                                                                                                 >  1:340998/1‑251 (MQ=255)                                                                                                        gACGTTGTGCCAGCAGCTGGCGGATATACATCGCGGTCGTACTGGCAAAGAAGGTGATGACGGCACCATCCCAGCCACCGTTATTGAGTTTACAGAAACAGGCGCAAGAAAGGCCAACCATTAAGACTACCAGCTATCTTGGGTAACGTAATGGCTGAATTTGGCTAAATCGTTTCTCTACGCCTTTGTAATCCAGCAGATGATGCTCCGCAAGAATCACAATGTGCTGGACTTCAGTCACCACATGCat                                                                                                                               <  2:319869/250‑1 (MQ=255)                                                                                                                    gcagcTGGCGGATATACATCGCGGTCGTACTGGCAAAGAAGGTGATGACGGCACCATCCCAGCCACCGTTATTGAGTTTACAGAAACAGGCGCAAGAAAGGCCAACCATTAAGGCTACCAGCTATCTTGGGTAACGTAATGGCTGAATTTGGCTAAATCGTTTCTCTACGCCTTTGTAATCCAGCAGATGATGCTCCGCAAGAATCACAATGTGCTGGACTTCAGTCACCACATGCATATTAATGCcgcg                                                                                                                   >  1:391288/1‑250 (MQ=255)                                                                                                                      agcTGGCGGATATACATCGCGGTCGTACTGGCAAAGAAGGTGATGACGGCACCATCCCAGCCACCGTTATTGAGTTTACAGAAACAGGCGCAAGAAAGGCCAACCATTAAGGCTACCAGCTATCTTGGGTAACGTAATGGCTGAATTTGGCTAAATCGTTTCTCTACGCCTTTGTAATCCAGCAGATGATGCTCCGCAAGAATCACAATGTGCTGGACTTCAGTCACCACATGCATATTAATGCCGCGATc                                                                                                                >  1:445931/1‑251 (MQ=255)                                                                                                                      agcTGGCGGATATACATCGCGGTCGTACTGGCAAAGAAGGTGATGACGGCACCATCCCAGCCACCGTTATTGAGTTTACAGAAACAGGCGCAAGAAAGGCCAACCATTAAGGCTACCAGCTATCTTGGGTAACGTAATGGCTGAATTTGGCTAAATCGTTTCTCTACGCCTTTGTAATCCAGCAGATGATGCTCCGCAAGAATCACAATGTGCTGGACTTCAGTCACCACATGCATATTAATGCCGCGATc                                                                                                                <  2:445113/251‑1 (MQ=255)                                                                                                                         tGGCGGATATACATCGCGGTCGTACTGGCAAAGAAGGTGATGACGGCACCATCCCAGCCACCGTTATTGAGTTTACAGAAACAGGCGCAAGAAAGGCCAACCATTAAGGCTACCAGCTATCTTGGGTAACGTAATGGCTGAATTTGGCTAAATCGTTTCTCTACGCCTTTGTAATCCAGCAGATGATGCTCCGCAAGAATCACAATGTGCTGGACTTCAGTCACCACATGCATATTAATGCCGCGATCGTg                                                                                                             >  2:315907/1‑251 (MQ=255)                                                                                                                              gATATACATCGCGGTCGTACTGGCAAAGAAGGTGATGACGGCACCATCCCAGCCACCGTTATTGAGTTTACAGAAACAGGCGCAAGAAAGGCCAACCATTAAGGCTACCAGCTATCTTGGGTAACGTAATGGCTGAATTTGGCTAAATCGTTTCTCTACGCCTTTGTAATCCAGCAGATGATGCTCCGCAAGAATCACAATGTGCTGGACTTCAGTCACCACATGCATATTAATGCCGCGATCGTGAtttt                                                                                                        >  2:387672/1‑251 (MQ=255)                                                                                                                                                     cAAAGAAGGTGATGACGGCACCATCCCAGCCACCGTTATTGAGGTTACAGAAACAGGCGCAAGAAAGGCCAACCATTAAGGCTACCAGCTATCTTGGGTAACGTAATGGCTGAATTTGGCTAAATCGTTTCTCTACGCCTTTGTAATCCAGCAGATGATGCTCCGCAAGAATCGCAATGTGCTGGTCTTCAGTCACCACATGCATATTAATGCCGCGATCGTGATTTTTACGGGTCGATGTCAGGCATTGc                                                                                 <  2:310617/251‑1 (MQ=255)                                                                                                                                                       aagaagGTGATGACGGCACCATCCCAGCCACCGTTATTGAGTTTACAGAAACAGGCGCAAGAAAGGCCAACCATTAAGGCTACCAGCTATCTTGGGTAACGTAATGGCTGAATTTGGCTAAATCGTTTCTCTACGCCTTTGTAATCCAGCAGATGATGCTCCGCAAGAATCACAATGTGCTGGACTTCAGTCACCACATGCATATTAATGCCGCGATCGTGATTTTTACGTGTCGATGTCAGGCATTGccc                                                                               <  1:391806/251‑1 (MQ=255)                                                                                                                                                        agaagGTGATGACGGCACCATCCCAGCCACCGTTATTGAGTTTACAGAAACAGGCGCAAGAAAGGCCAACCATTAAGGCTACCAGCTATCTTGGGTAACGTAATGGCTGAATTTGGCTAAATCGTTTCTCTACGCCTTTGTAATCCAGCCGATGATGCTCCGCAAGAATCACAATGTGCTGGACTTCAGTCACCACATGCATATTAATGCCGCGATCGTGATTTTTACGTGTCGATGTCAGGCATTGccc                                                                               >  2:455678/1‑250 (MQ=255)                                                                                                                                                               gatgaCGGCACCATCCCAGCCACCGTTATTGAGTTTACAGAAACAGGCGCAAGAAAGGCCAACCATTAAGGCTACCAGCTATCTTGGGTAACGTAATGGCTGAATTTGGCTAAATCGTTTCTCTACGCCTTTGTAATCCAGCAGATGATGCTCCGCAAGAATCACAATGTGCTGGACTTCAGTCACCACATGCATATTAATGCCGCGATCGTGATTTTTACGTGTCGATGTCAGGCATTGCCCATCTTTa                                                                        <  1:93056/250‑1 (MQ=255)                                                                                                                                                               gatgaCGGCACCATCCCAGCCACCGTTATTGAGTTTACAGAAACAGGCGCAAGAAAGGCCAACCATTAAGGCTACCAGCTATCTTGGGTAACGTAATGGCTGAATTTGGCTAAATCGTTTCTCTACGCCTTTGTAATCCAGCAGATGATGCTCCGCAAGAATCACAATGTGCTGGACTTCAGTCACCACATGCATATTAATGCCGCGATCGTGATTTTTACGTGTCGATGTCAGGCATTGCCCATCTTTa                                                                        >  1:254913/1‑250 (MQ=255)                                                                                                                                                                   aCGGCACCATCCCAGCCACCGTTATTGAGTTTACAGAAACAGGCGCAAGAAAGGCCAACCATTAAGGCTACCAGCTATCTTGGGTAACGTAATGGCTGAATTTGGCTAAATCGTTTCTCTACGCCTTTGTAATCCAGCAGATGATGCTCCGCAAGAATCACAATGTGCTGGACTTCAGTCACCACATGCATATTAATGCCGCGATCGTGATTTTTACGTGTCGATGTCAGGCATTGCCCATCTTTAATAGt                                                                   <  2:18659/251‑1 (MQ=255)                                                                                                                                                                        aCCATCCCAGCCCCCGTTATTGAGTTTACAGAAACAGGCGCAAGAAAGGCCAACCATTAAGGCTACCAGCTATCTTGGGTAACGTAATGGCTGAATTTGGCTAAATCGTTTCTCTACGCCTTTGTAATCCAGCAGATGATGCTCCGCAAGAATCACAATGTGCTGGACTTCAGTCACCACATGCATATTAATGCCGCGATCGTGATTTTTACGTGTCGATGTCAGTCATTGCCCATCTTTAATAGTTGTCa                                                              <  1:361680/251‑1 (MQ=255)                                                                                                                                                                              ccAGCCACCGTTATTGAGTTTACAGAAACAGGCGCAAGAAAGGCCAACCATTAAGGCTACCAGCTATCTTGGGAAACGTAATGGCTGAATTTGGCTAAATCGTTTCTCTACGCCTTTGTAATCCAGCAGATGATGCTCCGCAAGAATCACAATGTGCTGGACTTCAGTCACCACATGCATATTAATGCCGCGATCGTGATTTTTACGTGTCGATGTCAGGCATTGCCCATCTTTAATAGTTGTCAGCACt                                                         >  1:300426/1‑250 (MQ=255)                                                                                                                                                                                   cACCGTTATTGAGTTTACAGAAACAGGCGCAAGAAAGGCCAACCATTAAGGCTACCAGCTATCTTGGGTAACGTAATGGCTGAATTTGGCTAAATCGTTTCTCTACGCCTTTGTAATCCAGCAGATGATGCTCCGCAAGAATCACAATGTGCTGGACTTCAGTCACCACATGCATATGAATGCCGCGATCGTGATTTGTACGTGTCGATGTCAGGCATTGCCCATCGTTAATAGGTGTCAGCACTATGGCg                                                   <  2:24170/251‑1 (MQ=255)                                                                                                                                                                                          aTTGAGTTTACAGAAACAGGCGCAAGAAAGGCCAACCATTAAGGCTACCAGCTATCTTGGGTAACGTAATGGCTGAATTTGGCTAAATCGTTTCTCTACGCCTTTGTAATCCAGCAGATGATGCTCCGCAAGAATCACAATGTGCTGGACTTCAGTCACCACATGCATATTAATGCCGCGATCGTGATTTTTACGTGTCGATGTCAGGCATTGCCCATCTTTAATAGTTGTCAGCACTATGGCGTTCGaa                                             <  2:203756/250‑1 (MQ=255)                                                                                                                                                                                             gAGTTTACAGAAACAGGCGCAAGAAAGGCCAACCATTAAGGCTACCAGCTATCTTGGGTAACGTAATGGCTGAATTTGGCTAAATCGTTTCTCTACGCCTTTGTAATCCAGCAGATGATGCTCCGCAAGAATCACAATGTGCTGGACTTCAGTCACCACATGCATATTAATGCCGCGATCGTGATTTTTACGTGTCGATGTCAGGCATTGCCCATCTTTAATAGTTGTCAGCACTATGGCGTTCGAagag                                          <  1:307312/250‑1 (MQ=255)                                                                                                                                                                                               gTTTACAGAAACAGGCGCAAGAAAGGCCAACCATTAAGGCTACCAGCTATCTTGGGTAACGTAATGGCTGAATTTGGCTAAATCGTTTCTCTACGCCTTTGTAATCCAGCAGATGATGCTCCGCAAGAATCACAATGTGCTGGACTTCAGTCACCACATGCATATTAATGCCGCGATCGTGATTTTTACGTGTCGATGTCAGGCATTGCCCATCTTTAATAGTTGTCAGCACTATGGCGTTCGAAGAGAt                                        >  1:41504/1‑250 (MQ=255)                                                                                                                                                                                                              cgcAAGAAAGGCCAACCATTAAGGCTACCAGCTATCTTGGGTAACGTAATGGCTGAATTTGGCTAAATCGTTTCTCTACGCCTTTGTAATCCAGCAGATGATGCTCCGCAAGAATCACAATGTGTTGGACTTCAGTCACCACATGCATATTAATGCCGCGATCGTGATTTTTACGTGTCGATGTCAGGCATTGCCCATCTTTAATAGTTGTCAGCACTATGGCGTTCGAAGAGATAGAACTTTCGACGCTg                        >  2:314496/1‑251 (MQ=255)                                                                                                                                                                                                                                   aGGCTACCAGCTATCTTGGGTAACGTAATGGCTGAATTTGGCTAAATCGTTTCTCTACGCCTTTGTAATCCAGCAGATGATGCTCCGCAAGAATCACAATGTGCTGGACTTCAGTCACCACATGCATATTAATGCCGCGATCGTGATTTTTACGTGTCGATGTCAGGCATTGCCCATCTTTAATAGTTGTCAGCACTATGGCGTTCGAAGAGATAGAACTTTCGACGCTGTCCATTCCCAGTGCCCGAccc   <  1:314496/251‑1 (MQ=255)                                                                                                                                                                                                                                    ggCTACCAGCTATCTTGGGTAACGTAATGGCTGAATTTGGCTAAATCGTTTCTCTACGCCTTTGTAATCCAGCAGATGATGCTCCGCAAGAATCACAATGTGCTGGACTTCAGTCACCACATGCATATTAATGCCGCGATCGTGATTTTTACGTGTCGATGTCAGGCATTGCCCATCTTTAATAGTTGTCAGCACTATGGCGTTCGAAGAGATAGAACTTTCGACGCTGTCCATTCCCAGTGCCCGACCCa  <  2:193902/251‑1 (MQ=255)                                                                                                                                                                                                                                              |                                                                                                                                                                                                                                                 GCAATGGTGGGGGTATTGCTGAAAGTGGGGAGTTGCAAAAGCAATCCGGAAATGGTGGTGGCGGCGAAAGCGGTAAGGCAAAAGTTGATCTGTGGATGAAGATGACGTTGTGCCAGCAGCTGGCGGATATACATCGCGGTCGTACTGGCAAAGAAGGTGATGACGGCACCATCCCAGCCACCGTTATTGAGTTTACAGAAACAGGCGCAAGAAAGGCCAACCATTAAGGCTACCAGCCATCTTGGGTAACGTAATGGCTGAATTTGGCTAAATCGTTTCTCTACGCCTTTGTAATCCAGCAGATGATGCTCCGCAAGAATCACAATGTGCTGGACTTCAGTCACCACATGCATATTAATGCCGCGATCGTGATTTTTACGTGTCGATGTCAGGCATTGCCCATCTTTAATAGTTGTCAGCACTATGGCGTTCGAAGAGATAGAACTTTCGACGCTGTCCATTCCCAGTGCCCGACCCA  >  NC_000913/4602272‑4602749

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