Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	NC_000913	4,298,401	G→A	L322L (CTC→CTT)	fdhF ←	formate dehydrogenase H

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	4,298,401	0	G	A	98.6%	217.9 / NA	69	L322L (CTC→CTT)	fdhF	formate dehydrogenase H
Reads supporting (aligned to +/- strand):  ref base G (0/0);  major base A (32/36);  minor base C (1/0);  total (33/36)
Fisher's exact test for biased strand distribution p-value = 4.78e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

TCTTCGCCCATAATGTACGCGGCACGCACTTCGCCATGCGCTGCGCGGTGCGGCAGCTCGCTGATGCGATAGCCGGTATGCGCTGGCAGGCTTTCCACGCCCCAGGCTTTGGCGAATTTCTCGCGGTTAGCCGGATCTTTCACGTACTGGTATCCCGGATACGTATCCGGCAGCGCGCCCATATCGCAGGCACCCTGAACGTTGTTCTGACCACGAACCGGGTTAACACCCGCATGCGGCTTACCGAGGTTACCGGTCAGCATCGCGAGGCTGGTCAGAGAACGCACGGTTTCCACGCCCTGGTAGAACTGGGTTACACCCATGCCCCACAGGATGGCGGCGCTTTTCGCCTGGGCATACATCCGTGCCGCCTGACGAATCTCACTGGCGCTGACGCCGGTGATATCTTCAACCGACTCCGGCGTGTAGCCTTCAACGATTTTACGATACTCTTCAAAGCCTTCTGTACGTGAAGCGACGAACGCTT  >  NC_000913/4298156‑4298642                                                                                                                                                                                                                                                      |                                                                                                                                                                                                                                                  tCTTGGCCCATAAGGTGCGCGGCACGCCCTTCGGCAGGCGCTGCGCGGTGCGGCAGCGCGCTGATGCGATGTCCGGTATGCGCTGGCAGGCTTTCCACGCCCCAGGCTTTAGCGAATTTCTCGCGGTTAGCCGGATCTTTCACGTACTGGTATCCCGGATACGTATCCGGCAGCGCGCCCATATCGCAGGCACCCTGAACGTTGTTCTGACCACGAACCGGGTTAACACCCGCATGCGGCTTGCCAAGGt                                                                                                                                                                                                                                               <  1:270091/250‑1 (MQ=255)    tCGCCCATAAAGTACGCGGCACGCACGTCGCCATGCGCTGCGCTGTGCGGCAGCTCGCTGATGCGATAGCCGGTAGGCGCTGGCAGGCTTTCCACGCCCCAGGCTTTGGCGAATTTCTCGCGGTTAGCCGGATCTTTCACGTACTGGTATCCCGGATACGTATCCGGCAGCGCGCCCATATCGCCGGCACCCTGAACGTTGTTCTGACCACGAACCGGGTTACCACCCGCATGCGGCTTACCAAGGTTAc                                                                                                                                                                                                                                            <  2:335004/250‑1 (MQ=255)       cccATAATGTACGCGGCACGCACTTCGCCATGCGCTGCGCGGTGCGGCAGCTCGCTGATGCGATAGCCGGTATGCGCTGGCAGGCTTTCCACGCCCCAGGCTTTGGCGAATTTCTCGCGGTTAGCCGGATCTTTCACGTACTGGTATCCCGGATACGTATCCGGCAGCGCGCCCATATCGCAGGCACCCTGAACGTTGTTCTGACCACGAACCGGGTTAACACCCGCATGCGGCTTACCAAGGTTACCGGt                                                                                                                                                                                                                                        >  1:400609/1‑251 (MQ=255)             atggaCGCGGCACGCACTTCGCCATGCGCTGCGCGGTGCGGCAGCTCGCTGATGCGATAGCCGGTATGCGCTGGCAGGCTTTCCACGCCCCAGGCTTTGGCGAATTTCTCGCGGTTAGCCGGATCTTTCACGTACTGGTATCCCGGACACGTATCCGGCAGCGCGCCCCTATCGCAGGCACCCTGAACGTTGTTCTGACCACGAACCGGGTTAACACCCGCATGCGGCTTACCCAGGTTTCCGGTCAGCa                                                                                                                                                                                                                                   >  2:16397/5‑250 (MQ=255)               gTACGCGGCACGCACTTCGCCATGCGCTGCGCGGTGCGGCAGCTCGCTGATGCGATAGCCGGTATGCGCTGGCAGGCTTTCCACGCCCCAGGCTTTGGCGAATTTCTCGCGGTTAGCCGGATCTTTCACGTACTGGTATCCCGGATACGTATCCGGCAGCGCGCCCATATCGCAGGCACCCTGAACGTTGTTCTGACCACGAACCGGGTTAACACCCGCATGCGGCTTACCACGGTTACCGGTCAGCATcc                                                                                                                                                                                                                                >  2:233025/1‑250 (MQ=255)               gTACGCGGCACGCACTTCGCCATGCGCTGCGCGGTGCGGCAGCTCGCTGATGCGATAGCCGGTATGCGCTGGCAGGCTTTCCACGCCCCAGGCTTTGGCGAATTTCTCGCGGTTAGCCGGATCTTTCACGTACTGGTATCCCGGATACGTATCCGGCAGCGCGCCCATATCGCAGGCACCCTGAACGTTGTTCTGACCACGAACCGGGTTAACACCCGCATGCGGCTTACCAAGGTTACCGGTCAGCATc                                                                                                                                                                                                                                 <  1:184934/250‑1 (MQ=255)                tACGCGGCACGCACTTCGCCATGCGCTGCGCGGTGCGGCAGCTCGCTGATGCGATAGCCGGTATGCGCTGGCAGGCTTTCCACGCCCCAGGCTTTGGCGAATTTCTCGCGGTTAGCCGGATCTTTCACGTACTGGTATCCCGGATACGTATCCGGCAGCGCGCCCATATCGCAGGCACCCTGAACGTTGTTCTGACCACGAACCGGGTTAACACCCGCATGCGGCTTACCAAGGTTACCGGTCAGCATcgc                                                                                                                                                                                                                               >  1:370175/1‑251 (MQ=255)                             cTTCGCCATGCGCTGCGCGGTGCGGCAGCTCGCTGATGCGATAGCCGGTATGCGCTGGCAGGCTTTCCACGCCCCAGGCTTTGGCGAATTTCTCGCGGTTAGCCGGATCTTTCACGTACTGGTATCCCGGATACGTATCCGGCAGCGCGCCCATATCGCAGGCACCCTGAACGTTGTTCTGACCACGAACCGGGTTAACACCCGCATGCGGCTTACCAAGGTTACCGGTCAGCATCGCGAGGCTGGTc                                                                                                                                                                                                                     >  1:287891/1‑248 (MQ=255)                                      gcgctgcgcGGTGCGGCAGCTCGCTGATGCGATAGCCGGTATGCGCTGGCAGGCTTTCCACGCCCCAGGCTTTGGCGAATTTCTCGCGGTTAGCCGGATCTTTCACGTACTGGTATCCCGGATACGTATCCGGCAGCGCGCCCATATCGCAGGCACCCTGAACGTTGTTCTGACCACGAACCGGGTTAACACCCGCATGCGGCTTACCAAGGTTACCGGTCAGCATCGCGAGGCTGGTCAGAGAACGCAc                                                                                                                                                                                                          <  1:418649/250‑1 (MQ=255)                                        gcggcgcGGTGCGGCAGCTCGCTGATGCGATAGCCGGTATGCGCTGGCAGTCTTTCCACGCCCCAGGCTGTGGCGAATTTCTCGCGGTTAGCCGGATCTCTCACGTACTGGTATCCCGGATACGTAGCCGGCAGCGCGCCCATATCGCAGGCACCCTGAACGTTGTTCTGACCACGAACCGGGTTAACACCCGCATGCGGCTTACCAAGGTTACCGGTCAGCATCGCGAGGCTGGTCAGAGAACGCACGGt                                                                                                                                                                                                       <  2:370175/248‑1 (MQ=255)                                               ggTGCGGCAGCTCGCTGATGCGATAGCCGGTATGCGCTGGCAGGCTTTCCACGCCCCAGGCTTTGGCGAATTTCTCGCGGTTAGCCGGATCTTTCACGTACTGGTATCCCGGATACGTATCCGGCAGCGCGCCCATATCGCAGGCACCCTGAACGTTGTTCTGACCACGAACCGGGTTAACACCCGCATGCGGCTTACCAAGGTTACCGGTCAGCATCGCGAGGCTGGTCAGAGAACGCACGGTTTCCACg                                                                                                                                                                                                <  2:153987/251‑1 (MQ=255)                                                                  gCGATAGCCGGTATGCGCTGGCAGGCTTTCCACGCCCCAGGCTTTGGCGAATTTCTCGCGGTTAGCCGGATCTTTCACGTACTGGTATCCCGGATACGTATCCGGCAGCGCGCCCATATCGCAGGCACCCTGAACGTTGTTCTGACCACGAACCGGGTTAACACCCGCATGCGGCTTACCAAGGTTACCGGTCAGCATCGCGAGGCTGGTCAGAGAACGCACGGTTTCCACGCCCTGGTAGAACTGGGTTa                                                                                                                                                                             <  1:341129/251‑1 (MQ=255)                                                                       agccggttgtcgcTGGCAGGCTTTCCACGCCCCAGGCTTTGGCGAATTTCTCGCGGTTAGCCGGATCTTTCACGTACTGGTATCCCGGATACGTATCCGGCAGCGCGCCCATATCGCAGGCACCCTGAACGTTGTTCTGACCACGAACCGGGTTAACACCCGCATGCGGCTTACCAAGGTTACCGGTCAGCATCGCGAGGCTGGTCAGAGAACGCACGGTTTCCACGCCCTGGTAGAACTGGGTTACAccc                                                                                                                                                                        <  2:140335/241‑1 (MQ=255)                                                                            gTATGCGCTGGCAGGCTTTCCACGCCCCAGGCTTTGGCGAATTTCTCGCGGTTAGCCGGATCTTTCACGTACTGGTATCCCGGATACGTATCCGGCAGCGCGCCCATATCGCAGGCACCCTGAACGTTGTTCTGACCACGAACCGGGTTAACACCCGCATGCGGCTTACCAAGGTTACCGGTCAGCATCGCGAGGCTGGTCAGAGAACGCACGGTTTCCACGCCCTGGTAGAACTGGGTTACACCCATGcc                                                                                                                                                                   <  2:6243/251‑1 (MQ=255)                                                                                gggcgtgccggCTTTCCACGCCCCGGGCTTGTCCGAATTCCCCGCGGGTAGCCGGACCTTGCACGTACTGGTATCCCGGATTCGTATCCGGCAGCGCGCCCATATCGCAGGCACCCTGAACGGTGTTCTGACCTCGAACCGGGTTAACACCCGCATGCGGCCTACCAAGGGTACCGGTCAGCAGCGCGAGGCTGGTCAGAGAACGCACGGTTTCCACGCCCTGGTAGAACTGGGTTACACCCATGCCCca                                                                                                                                                                <  2:55164/241‑1 (MQ=255)                                                                                   cTGGCAGGCTTTCCACGCCCCAGGCTTTGGCGAATTTCTCGCGGTTAGCCGGATCTTTCACGTACTGGTATCCCGGATACGTATCCGGCAGCGCGCCCATATCGCAGGTACCCTGAACGTTGTTCTGACCACGAACCGGGTTAACACCCGCATGCGGCTTACCAAGGTTACCGGTCAGCATCGCGAGGCTGGTCAGAGAACGCACGGTTTCCACGCCCTGGTAGAACTGGGTTACACCCATGCCCCACAgg                                                                                                                                                            >  1:427432/1‑251 (MQ=255)                                                                                   cTGGCAGGCTTTCCACGCCCCAGGCTTTGGCGAATTTCTCGCGGTTAGCCGGATCTTTCACGTACTGGTATCCCGGATACGTATCCGGCAGCGCGCCCATATCGCAGGCACCCTGAACGTTGTTCTGACCACGAACCGGGTTAACACCCGCATGCGGCTTACCAAGGTTACCGGTCAGCATCGCGAGGCTGGTCAGAGAACGCACGGTTTCCACGCCCTGGTAGAACTGGGTTACACCCATGCCCCACAgg                                                                                                                                                            >  1:461270/1‑251 (MQ=255)                                                                                     ggccggccttCCACGCCCCCGGCTTTGGCGAATTTCTCGCGGTTAGCCGTATCGTTCACGTCCTGGTCTCCCGGATTCGGATCCGGCAGCGCGCCCATAGCGCAGGCACCCGGAACGTTGTTCTGACCACGAACCGGGTTAACACCCGCATGCGGCTTACCAAGGTTACCGGTCAGCATCGCGAGGCTGGTCAGAGAACGCACGGTTTCCACGCCCTGGTAGAACTGGGTTACACCCATGCCCCACAgg                                                                                                                                                            <  2:168429/241‑1 (MQ=255)                                                                                        aGGCTTTCCACGCCCCAGGCTTTGGCGAATTTCTCGCGGTTAGCCGGATCTTTCACGTACTGGTATCCCGGATACGTATCCGGCAGCGCGCCCATATCGCAGGCACCCTGAACGTTGTTCTGACCACGAACCGGGTTAACACCCGCATGCGGCTTACCAAGGTTACCGGTCAGCATCGCGAGGCTGGTCAGAGAACGCACGGTTTCCACGCCCTGGTAGAACTGGGTTACACCCATGCCCCACAGGATggc                                                                                                                                                       >  1:44284/1‑251 (MQ=255)                                                                                               ccACGCCCCAGGCTTTGGCGAATTTCTCGCGGTTAGCCGGATCTTTCACGTACTGGTATCCCGGATACGTATCCGGCAGCGCGCCCATATCGCAGGCACCCTGAACGTTGTTCTGACCACGAACCGGGTTAACACCCGCATGCGGCTTACCAAGGTTACCGGTCAGCATCGCGAGGCTGGTCAGAGAACGCACGGTTTCCACGCCCTGGTAGAACTGGGTTacac                                                                                                                                                                          <  2:196448/225‑1 (MQ=255)                                                                                               ccACGCCCCAGGCTTTGGCGAATTTCTCGCGGTTAGCCGGATCTTTCACGTACTGGTATCCCGGATACGTATCCGGCAGCGCGCCCATATCGCAGGCACCCTGAACGTTGTTCTGACCACGAACCGGGTTAACACCCGCATGCGGCTTACCAAGGTTACCGGTCAGCATCGCGAGGCTGGTCAGAGAACGCACGGTTTCCACGCCCTGGTAGAACTGGGTTacac                                                                                                                                                                          >  1:196448/1‑225 (MQ=255)                                                                                               ccACGCCCCAGGCTTTGGCGAATTTCTCGCGGTTAGCCGGATCTTTCACGTACTGGTATCCCGGATACGTATCCGGCAGCGCGCCCATATCGCAGGCACCCTGAACGTTGTTCTGACCACGAACCGGGTTAACACCCGCATGCGGCTTACCAAGGTTACCGGTCAGCATCGCGAGGCTGGTCAGAGAACGCACGGTTTCCACGCCCTGGTAGAACTGGGTTACACCCATGCCCCACAGGATGGCGGCGCtt                                                                                                                                                <  1:52260/251‑1 (MQ=255)                                                                                               ccACGCCCCAGGCTTTGGCGAATTTCTCGCGGTTAGCCGGATCTTTCACGTACGGGTAGCCCGGATACGTATCCGGCAGCGCGCCCATAGCGCAGGCACCCTGAACGTTGTTCTGACCACGAACCGGGTTAACACCCGCATGCGGCTTACCAAGGTTACCGGTCAGCATCGCGAGGCTGGTCAGAGAACGCACGGTTTCCACGCCCTGGTAGAACTGGGTTACACCCATGCCCCACAGGATGGCGGCGCtt                                                                                                                                                <  2:483717/251‑1 (MQ=255)                                                                                                      ccAGGCTTTGGCGAATTTCTCGCGGTTAGCCGGATCTTTCACGTACTGGTATCCCGGATACGTATCCGGCAGCGCGCCCATATCGCAGGCACCCTGAACGTTGTTCTGACCACGAACCGGGTTAACACCCGCATGCGGCTTACCAAGGTTACCGGTCAGCATCGCGAGGCTGGTCAGAGAACGCACGGTTTCCACGCCCTGGTAGAACTGGGTTACACCCATGCCCCACAGGATGGCGGCGCTTTTCGCCt                                                                                                                                         <  1:46836/251‑1 (MQ=255)                                                                                                        aGGCTATGGCGAATTTCTCGCGGTTAGACGGATCTTTCACGAACTGGTATCCCGGATAAGTATCCGGCAGCGCGCCCATATCGCAGGCACCCTGAACGTTGTTCTGACCACGAACCGGGTTAACACCCGCATGAGGCTGACCAAGGTTACCGGTCAGCATCGCGAGGCTGGTCAGAGAACGCACGGTTTCCACGCCCTGGTAGAACTGAGTTAAACCCATGCCCCACAGGATGGCGGCGCTTTTAGCCTgg                                                                                                                                       >  2:328743/1‑251 (MQ=255)                                                                                                         ggCTTTGGCGAATTTCTCGCGGTTAGCCGGATCTTTCACGTACTGGTATCCCGGATACGTATCCGGCAGCGCGCCCATATCGCAGGCACCCTGAACGTTGTTCTGACCACGAACCGGGTTAACACCCGCATGCGGCTTACCAAGGTTACCGGTCAGCATCGCGAGGCTGGTCAGAGAACGCACGGTTTCCACGCCCTGGTAGAACTGGGTTACACCCATGCCCCACAGGATGGCGGCGCTTTTCGCCTggg                                                                                                                                      <  1:214095/251‑1 (MQ=255)                                                                                                             ttGGCGAATTTCTCGCGGTTAGCCGGATCTTTCACGTACTGGTATCCCGGATACGTATCCGGCAGCGCGCCCATATCGCAGGCACCCTGAACGTTGTTCTGACCACGAACCGGGTTAACACCCGCATGCGGCTTACCAAGGTTACCGGTCAGCATCGCGAGGCTGGTCAGAGAACGCACGGTTTCCACGCCCTGGTAGAACTGGGTTACACCCATGCCCCACAGGATGGCGGCGCTTTTCGCCTGGGCAt                                                                                                                                   <  2:189475/250‑1 (MQ=255)                                                                                                                    attgatccgggTTAGCCGGATCTTTCACGTACTGGTAGCCCGGATACGTATCCGGCAGCGCGCCCATATCGCAGGCGCCCTGAACGTTGTTCTGACCACGAACCGGGTTAGCACCCGCAGGCGGCGTCCCAAGGTTTCCGGTCAGCATCGCGAGGCTGGTCAGAGAACGCGCGGTTTCCACGCCCTGGTAGAACTGGGTTACACCCATGCCCCACAGGATGGCGGCGCTTTGCGCCTGGGCATACATCCGt                                                                                                                           <  2:359723/242‑1 (MQ=255)                                                                                                                       tcttgcgGTTAGCCGGATCTTTCCTGTACTGGTATCCCGGATACGTATCCGGCAGCGCGCCCATATCGCAGGCACCCTGAACGTTGTTCTGACCACGAACCGGGTTAACACCCGCATGCGGCTTACCAAGGGTACCGGTCAGCATCGCGAGGCTGGTCAGAGAACGCAAGGTTTCCACGCCCTGGTAGAACTGGGTTACACCCATGCCCCACAGGATGGCCGCGCTTTTCGCCTGGGCATACATCCGTgc                                                                                                                         <  2:252453/246‑1 (MQ=255)                                                                                                                                     ggATCTTTCACGTACTGGTATCCCGGATACGTATCCGGCAGCGCGCCCATATCGCAGGCACCCTGAACGTTGTTCTGACCACGAACCGGGTTAACACCCGCATGCGGCTTACCAAGGTTACCGGTCa                                                                                                                                                                                                                                      >  1:199685/1‑127 (MQ=255)                                                                                                                                     ggATCTTTCACGTACTGGTATCCCGGATACGTATCCGGCAGCGCGCCCATATCGCAGGCACCCTGAACGTTGTTCTGACCACGAACCGGGTTAACACCCGCATGCGGCTTACCAAGGTTACCGGTCAGCATCGCGAGGCTGGTCAGAGAACGCACGGTTTCCACGCCCTGGTAGAACTGGGTTACACCCATGCCCCACAGGATGGCGGCGCTTTTCGCCTGGGCATACATCCGTACCGCCTGACGAAtctc                                                                                                          >  2:92339/1‑251 (MQ=255)                                                                                                                                     ggATCTGTCACGTACTGGTATCCCTGATACGTATCCGGCAGCGCGCCCATATCGCAGGCACCCTGAACGTTGTTCTGACCACGAACCGGGTTAACACCCGCATGCGGCTTACCAAGGTTACCGGTCa                                                                                                                                                                                                                                      <  2:199685/127‑1 (MQ=255)                                                                                                                                         cTTTCACGTACTGGTATCCCGGATACGAATCCGGCAGCGCGCCCATATCGCAGGCACACTGAACGTTGTTCTGACCACGAACCGGGTTAACACCCGCATGCGGCTTACCAAGGTTACCGGTCAGCATCGCGAGGCTGGTCAGAGAACGCACGGTTTCCACGCCCTGGTAGAACTGGGTTAAACCCATGCCCCACAGGATGCCGGCGCTTTTCGCCTAGGCATACATCCGTGCCGCCTGACGAATCTCACTg                                                                                                      >  1:144662/1‑251 (MQ=255)                                                                                                                                           ttCACGTACTGGTATCCCGGATACGTATCCGGCAGCGCGCCCATATCGCAGGCACCCTGAACGTTGTTCTGACCACGAACCGGGTTAACACCCGCATGCGGCTTACCAAGGTTACCGGTCAGCATCGCGAGGCTGGTCAGAGAACGCACGGTTTCCACGCCCTGGTAGAACTGGGTTACACCCATGCCCCACAGGATGGCGGCGCTTTTCGCCTGGGCATACATCCGTGCCGCCTGACGAATCTCACTg                                                                                                      <  1:148389/249‑1 (MQ=255)                                                                                                                                                gTACTGGTATCCCGGATACGTATCCGGCAGCGCGCCCATATCGCAGGCACCCTGAACGTTGTTCTGCCCACGAACCGGGTTAACACCCGCATGCGGCTTACCAAGGTTACCGGTCAGCATCGCGAGGCTGGTCAGAGAACGCACGGTTTCCACGCCCTGGTAGAACTGGGTTACACCCATGCCCCACAGGATGGCGGCGCTTTTCGCCTGGGCATACATCCGTGCCGCCTGACGAATCTCCCTGACGCTgc                                                                                               >  2:105815/1‑250 (MQ=255)                                                                                                                                                          cccGGATACGTATCCGGCAGCGCGCCCATATCGCAGGCACCCTGAACGTTGTTCTGACCACGAACCGGGTTAACACCCGCATGCGGCTTACCAAGGTTACCGGTCAGCATCGCGAGGCTGGTCAGAGAACGCACGGTTTCCACGCCCTGGTAGAACTGGGTTACACCCATGCCCCACAGGATGGCGGCGCTTTTCGCCTGGGCATACATCCGTGCCGCCTGACGAATCTCACTGGCGCTGCCGCCGGGGat                                                                                     >  2:97530/1‑251 (MQ=255)                                                                                                                                                               aTACGTATCCGGCAGCGCGCCCCTATCGCAGGCACCCTGAACGTTGTTCTGACCACGAACCGGGTTAACACCCGCATGCGGCTTACCAAGGTTACCGGTCAGCATCGCGAGGCTGGTCAGAGAACGCACGGTTTCCACGCCCTGGTAGAACTGGGTTACACCCATGCCCCACAGGATGGCGGCGCTTTTCGCCTGGGCCTACATCCGTGCCGCCTGACGAATCTACCTGGCGCTGACGCCGGTGatatatt                                                                                >  2:76855/1‑248 (MQ=255)                                                                                                                                                                 aCGTATCCGGCAGCGCGCCCATATCGCAGGCACCCTGAACGTTGTTCTGACCACGAACCGGGTTAACACCCGCATGCGGCTTACCAAGGTTCCCGGTCAGCATCGCGAGGCTGGTCAGAGAACGCACGGTTTCCACGCCCTGGTAGAACTGGGTTACACCCATGCCCCACAGGATGGCGGCGCTTTTCGCCTGGGCATACATCCGTGCCGCCTGACGAATCTCACTGGCGCTGACGCCGGTGATATCTTc                                                                               >  2:134805/1‑250 (MQ=255)                                                                                                                                                                 aCGTATCCGGCAGCGCGCCCATATCGCAGGCACCCTGAACGTTGTTCTGACCACGAACCGGGTTAACACCCGCATGCGGCTTACCAAGGTTACCGGTCAGCATCGCGAGGCTGGTCAGAGAACGCACGGTTTCCACGCCCTGGTAGAACTGGGTTACACCCATGCCCCACAGGATGGCGGCGCTTTTCGCCTGGGCATACATCCGTGCCGCCTGACGAATCTCACTGGCGCTGACGCCGGTGATATCTTc                                                                               >  1:4493/1‑250 (MQ=255)                                                                                                                                                                     aTCCGGCAGCGCGCCCCTATCGCCGGCACCCTGAACGTTGTTCTGACCACGAACCGGGTTAACACCCGCATGCGGCTTACCAAGGTTACCGGTCAGCATCGCGAGGCTGGTCAGAGCACGCACGGTTTCCACGCCCTGGAAGAACTGGGTTACACCCATGCCCCACAGGATGGCGGCGCTTTTCGCCTGGGCATACATCCGTGCCGCCTGACGAATCTCACTGGCGATGACGCCGGTGATATCTTCAACCg                                                                          >  2:438271/1‑251 (MQ=255)                                                                                                                                                                       ccGGCAGCGCGCCCATATCGCAGGCACCCTGAACGTTGTTCTGACCACGAACCGGGTTAACACCCGCATGCGGCTTACCAAGGTTACCGGTCAGCATCGCGAGGCTGGTCAGAGAACGCACGGTTTCCACGCCCTGGTAGAACTGGGTTACACCCATGCCCCACAGGATGGCGGCGCTTTTCGCCTGGGCATACATCCGTGCCGCCTGACGAATCTCACTGGCGCTGACGCCGGTGATATCTTCAACCGAc                                                                        <  1:72466/251‑1 (MQ=255)                                                                                                                                                                        cGGCAGCGCGCCCATATCGCAGGCACCCTGAACGTTGTTCTGACCACGAACCGGGGTAACTCCCGCATGCGGCTTACCAAGGTTGCCGGTCAGCATCGCTAGGCTGGTCAGAGAACGCACGGTTTCCACGCCCTGGTAGAACTGGGTTACACCCATGCCCCACAGGATGGCGGCGCTTTTCGCCTGGGCATACATCCGTGCCGCCTGACGAATCTCACTGGCGCTGACGCCGGTGATATCTTCAACCGACt                                                                       <  2:464872/251‑1 (MQ=255)                                                                                                                                                                              cgcgcCCATATCGCAGGCACCCTGAACGTTGTTCTGACCACGAACCGGGTTAACACCCGCATGCGGCTTACCAAGGTTACCGGTCAGCATCGCGAGGCTGGTCAGAGAACGCACGGTTTCCACGCCCTGGTAGAACTGGGTTACACCCATGCCCCACAGGATGGCGGCGCTTTTCGCCTGGGCATACATCCGTGCCGCCTGACGAATCTCACTGGCGCTGACGCCGGTGATATCTTCAACCGACTCCGGCg                                                                 >  2:189094/1‑251 (MQ=255)                                                                                                                                                                                 gcCCATATCGCAGTCACCCTGAACGTTGTTCTGACCCCGAACCGGGTTAACACCCGCATGCGGCTTACCAAGGTTACCGGTCAGCATCGCGAGACTGGTCAGAGAACGCACGGTTTCCACGCCCTGGTAGAACGGGGTTACACCCATGCCCCACAGGATGGCGGCGCTTTTCGCCTGGGCATACATCCGTGCCGCCTGACGAATCTCACTGGCGCTGACGCCGGCGATATCTTCAAACGACTCCGGCgtgt                                                              >  2:483369/1‑251 (MQ=255)                                                                                                                                                                                 gcCCATATCGCAGGCACCCTGAACGTTGTTCTGACCACGAACCGGGTTAACACCCGCATGCGGCTTACCAAGGTTACCGGTCAGCATCGCGAGGCTGGTCAGAGAACGCACGGTTTCCACGCCCTGGTAGAACTGGGTTACACCCATGCCCCACAGGATGGCGGCGCTTTTCGCCTGGGCATACATCCGTGCCGCCTGACGAATCTCACTGGCGCTGCCGCCGGTGATATCTTCAACCGACTCCGGCgtgg                                                              >  2:208856/1‑250 (MQ=255)                                                                                                                                                                                         cGCAGGCACCCTGAACGTTGTTCTGACCACGAACCGGGTTAACACCCGCATGCGGCTTACCAAGGTTACCGGTCAGCATCGCGAGGCTGGTCAGAGAACGCACGGTTTCCACGCCCTGGTAGAACTGGGTTACACCCATGCCCCACAGGATGGCGGCGCTTTTCGCCTGGGCATACATCCGTGCCGCCTGACGAATCTCACTGGCGCTGACGCCGGTGATATCTTCAACCGACTCCGGCGTGTAGCCTTCa                                                      <  2:44284/251‑1 (MQ=255)                                                                                                                                                                                            aGGCACCCTGACCGTTGTTCTGACCACGAACCGGGTTAACACCCGCATGCGGCTTACCAAGGTTACCGGTCAGCATCGCGAGGCTGGTCAGAGAACGCACGGTTTCCACGCCCTGGTAGAACTGGGTTACACCCATGCCCCACAGGATGGCGGCGCTTTTCGCCTGGGCATACATCCGTGCTGCCTGACGAATCTCACTGGCGCTGACGCCGGTGATATCTTCAACCGACTCCGGCGTGTAGCCTTCAACg                                                   <  1:301838/251‑1 (MQ=255)                                                                                                                                                                                            aGGCACCCTGAACGTTGTTCTGACCACGAACCGGGTTAACACCCGCATGCGGCTTACCAAGGTTACCGGTCAGCATCGCGAGGCTGGTCAGAGAACGCACGGTTTCCACGCCCTGGTAGAACTGGGTTACACCCATGCCCCACAGGATGGCGGCGCTTTTCGCCTGGGCATACATCCGTGCCGCCTGACGAATCTCACTGGCGCTGACGCCGGTGATATCTTCAACCGACTCCGGCGTGTAGCCTTCAACg                                                   <  2:217709/251‑1 (MQ=255)                                                                                                                                                                                                    tGAACGTTGTTCTGACCACGAACCGGGTTAACACCCGCATGCGGCTTACCAAGGTTACCGGTCAGCATCGCGAGGCTGGTCAGAGAACGCACGGTTTCCACGCCCTGGTAGAACTGGGTTACACCCATGCCCCACAGGATGGCGGCGCTTTTCGCCTGGGCATACATCCGTGCCGCCTGACGAATCTCACTGGCGCTGACGCCGGTGATATCTTCAACCGCCTCCGGCGTGTAGCCTTCAACGATTTTACg                                           >  2:489268/1‑251 (MQ=255)                                                                                                                                                                                                        cGTTGTTCTGACCACGAACCGGGTTAACACCCGCATGCGGCTTACCAAGGTTACCGGTCAGCATCGCGAGGCTGGTCAGAGAACGCACGGTTTCCACGCCCTGGTAGAACTGGGTTACACCCATGCCCCACAGGATGGCGGCGCTTTTCGCCTGGGCATACATCCGTGCCGCCTGACGAATCTCACTGGCGCTGACGCCGGTGATATCtt                                                                                >  2:123570/1‑210 (MQ=255)                                                                                                                                                                                                        cGTTGTTCTGACCACGAACCGGGTTAACACCCGCATGCGGCTTACCAAGGTTACCGGTCAGCATCGCGAGGCTGGTCAGAGAACGCACGGTTTCCACGCCCTGGTAGAACTGGGTTACACCCATGCCCCACAGGATGGCGGCGCTTTTCGCCTGGGCATACATCCGTGCCGCCTGACGAATCTCACTGGCGCTGACGCCGGTGATATCtt                                                                                <  1:123570/210‑1 (MQ=255)                                                                                                                                                                                                         gttgttCTGACCACGAACCGGGTTAACACCCGCATGCGGCTTACCAAGGTTACCGGTCAGCATCGCGAGGCTGGTCAGAGAACGCACGGTTTCCACGCCCTGGTAGAACTGGGTTACACCCATGCCCCACAGGATGGCGGCGCTTTTCGCCTGGGCATACATCCGTGCCGCCTGACGAATCTCACTGGCGCTGACGCCGGTGATATCTTCAACCGACTCCGGCGTGTAGCCTTCaa                                                     <  1:398854/236‑1 (MQ=255)                                                                                                                                                                                                         gttgttCTGACCACGAACCGGGTTAACACCCGCATGCGGCTTACCAAGGTTACCGGTCAGCATCGCGAGGCTGGTCAGAGAACGCACGGTTTCCACGCCCTGGTAGAACTGGGTTACACCCATGCCCCACAGGATGGCGGCGCTTTTCGCCTGGGCATACATCCGTGCCGCCTGACGAATCTCACTGGCGCTGACGCCGGTGATATCTTCAACCGACTCCGGCGTGTAGCCTTCaa                                                     >  2:398854/1‑236 (MQ=255)                                                                                                                                                                                                         gttgttCTGACCACGAACCGGGTTAACACCCGCATGCGGCTTACCAAGGTTACCGGTCAGCATCGCGAGGCTGGTCAGAGAACGCACGGTTTCCACGCCCTGGTAGAACTGGGTTACACCCATGCCCCACAGGATGGCGGCGCTTTTCGCCTGGGCATACATCCGTGCCGCCTGACGAATCTCACTGGCGCTGACGCCGGTGATATCTTCAACCGACTCCGGCGTGTAGCCTTCAACGATTTTACGATAct                                      >  1:476135/1‑251 (MQ=255)                                                                                                                                                                                                         gttgttCTGACCACGAACCAGGTTAACTCCCGCATGCGGCTTACCAAGGTTACCGGTCAGCATCGCGAGGCTGGTCAGAGAACGCACGGTTTCCACGCCCTGGTAGAACTGGGTTACACCCATGCCCCACAGGATGGCGGCGCTTTTCGCCTGGGCATACATCCGTGCCGCCTGACGAATCTCACTGGCGCTGACGCCGGTGATATCTTCAACCGACACCGGCGTGTAGCCTTCAACGATTTTACGATAct                                      <  2:461270/251‑1 (MQ=255)                                                                                                                                                                                                            gttCTGACCACGAACCGGGTTAACACCCGCATGCGGCTTACCAAGGTTACCGGTCAGCATCGCGAGGCTGGTCAGAGAACGCACGGTTTCCACGCCCTGGTAGAACTGGGTTACACCCATGCCCCACAGGATGGCGGCGCTTTTCGCCTGGGCATACATCCGTGCCGCCTGACGAATCTCACTGGCGCTGACGCCGGTGATATCTTCAACCGACTCCGGCGTGTAGCCTTCAACGATTTTACGATACTCtt                                   <  1:74124/251‑1 (MQ=255)                                                                                                                                                                                                                tGACCACGAACCGGGTTAACACCCGCATGCGGCTTACCAAGGTTACCGGTCAGCATCGCGAGGCTGGTCAGAGAACGCACGGTTTCCACGCCCTGGAAGAACTGGGTAACACCCATGCCCCACAGGATGGCGGCGCTTTTCGCCTGGGCATACATCCGTGCCGCCTGACGAATCTCACTGGCGCTGACGCCGGTGATATCTTAAACCGACTCCGGCGTGTAGCCTTCAACGAATTTCAGATCCTCTTCaa                                >  2:381594/1‑250 (MQ=255)                                                                                                                                                                                                                      cGAACCGGGTTAACACCCGCATGCGGCTTACCAAGGTTACCGGTCAGCATCGCGAGGCTGGTCAGAGAACGCACGGTTTCCACGCCCTGGTAGAACTGGGTTACACCCATGCCCCACAGGATGGCGGCGCTTTTCGCCTGGGCATACATCCGTGCCGCCTGACGAATCTCACTGGCGCTGACGCCGGTGATATCTTCAACCGACTCCGGCGTGTAGCCTTCAACGATTTTACGATACTCTTCAAAGCCTTc                         <  1:89310/251‑1 (MQ=255)                                                                                                                                                                                                                            gggTTAACACCCGCATGCGGCTTACCAAGGTTACCGGTCAGCATCGCGAGGCTGGTCAGAGAACGCACGGTTTCCACGCCCTGGTAGAACTGGGTTACACCCATGCCCCACAGGATGGCGGCGCTTTTCGCCTGGGCATACATCCGTGCCGCCTGACGAATCTCACTGGCGCTGACGCCGGTGATATCTTCAACCGACTCCGGCGTGTAGCCTTCAACGATTTTACGATACTCTTCAAAGCCTTCTGTACg                   >  1:136237/1‑251 (MQ=255)                                                                                                                                                                                                                                  acacCCGCATGCGGCTTACCAAGGTTACCGGTCAGCATCGCGAGGCTGGTCAGAGAACGCACGGATTCCACGCCCTGGTAGAACTGGGTTACACCCATGCCCCACAGGATGGCGGCGCTTTTCGCCTGGGCATACATCCGTGCCGCCTGACGAATCTCACTGGCGCTGACGCCGGTGATATCTTCAACCGACTCCGGCGTGGAGCCTTCAACGATTTTACGATACTCTTCAAAGCCTTCTGTACGTGAAGc             >  1:183940/1‑251 (MQ=255)                                                                                                                                                                                                                                   cacCCGCATGCGGCTTACCAAGGTTACCGGTCAGCATCGCGAGGCTGGCCAGAGAACGCACGGTTTCCACGCCCTGGAAGAACTGGGTTACACCCATGCCCCACAGGATGGCGGCGCTTTTCGCCTGGGCATACATCCGTGCCGCCTGACGAATCTCACTGGCGCTGACGCCGGTGATATCTTCAACCGACCCCGGCGTGTAGCCTTCAACGATTTTCCGATACTCTTCCACGCCTTCTGAACGTGAAGcg            >  2:273674/1‑251 (MQ=255)                                                                                                                                                                                                                                     cccGCATGCGGCTTACCAAGGTTACCGGTCAGCATCGCGTGGCTGGTCAGAGAACGCACGGTTTCCACGCCCTGGTAGAACTGGGTTAAACCCATGCCCCACAGGATGGCGGCGCTTTTCGCCTGGGCATACATCCGTGCCGCCTGACGAATCTCACTGGCGCTGACGCCGGTGATATCTTCAACCGACTCCGGCGTGTAGGCTTc                                                       >  1:376956/1‑206 (MQ=255)                                                                                                                                                                                                                                     cccGCATGCGGCTTACCAAGGTAACCGGTCAGCATCGCGAGGCGGGTCAGAGAACGCACGGTTTCCACTCCCTGGTAGAACTGGGTTACACCCAAGCCCCACAGGATTGCGGCGCTTTTCGCCTGGGCATACATCCGTGCCGCCTGACGAATCTCACTGGCGCTGACGCCGGTGATATCTTCAACCGACTCCGGCGTGTAGCCTTc                                                       <  2:376956/206‑1 (MQ=255)                                                                                                                                                                                                                                      ccGCATGCGGCTTACCAAGGTTACCGGTCAGCATCGCGAGGCTGGTCAGAGAACGCACGGTTTCCACGCCCTGGTAGAACTGGGTTACACCCATGCCCCACAGGATGGCGGCGCTTTTCGCCTGGGCATACATCCGTGCCGCCTGACGAATCTCACTGGCGCTGACGCCGGTGATATCTTCAACCGACTCCGGCGTGTAGCCTTCAACGATTTTACGATACTCTTCAAAGCCTTCTGTACGTGAAGcgacg         <  1:285991/251‑1 (MQ=255)                                                                                                                                                                                                                                            gcgccTTACCAAGGTTACCGGTCAGCATCGCGAGGCTGGTCAGAGAACGCACGGTTTCCACGCCCTGGTAGAACTGGGTTTCACCCATGCCCAACAGGAGGGCGGCGTTTTTCGCCTGGGCATACATCCGTGCCGCCTGACGAATCGTACTGGCGTTTACGCCGGTGATATCTTCAACCGACTCCGGCGTGTAGCCTGCAACGATGTTACGATTCTCTTCAAAGCCTTCTGTACGTGAAGCGACGAACGCt   <  1:71587/247‑1 (MQ=255)                                                                                                                                                                                                                                            gCGGCTTACCAAGGTTCCCGGTCAGCATCGCGAGGCTGGTCAGAGAACGCACGGTTTCCACGCCCTGGTAGAACTGGGTTACACCCATGCCCCACAGGATGGCGGCGCTTTTCGCCTGGGCATACATCCGTGCCGCCTGACGAATCTCACTGGCGCTGACGCCGGTGATATCTTCAACCGACTCCGGCGTGTAGCCTTCAACGATTTTACGATACTCTTCAAAGCCTTCTGTACGTGAAGCGACGAACGCt   <  1:362340/251‑1 (MQ=255)                                                                                                                                                                                                                                            gCGGCTTACCAAGGTTACCGGTCAGCATCGCGAGGCTGGTCAGAGAACGCACGGTTTCCACGCCCTGGTAGAACTGGGTTACACCCATGCCCCACAGGATGGCGGCGCTTTTCGCCTGGGCATACATCCGTGCCGCCTGACGAATCTCACTGGCGCTGACGCCGGTGATATCTTCAACCGACTCCGGCGTGTAGCCTTCACCGCTTTTACGATACTCTTCAAAGCCTTCTGTACGTGAAGCGGCGAACGCt   >  2:231702/1‑251 (MQ=255)                                                                                                                                                                                                                                            gCGGCTTACCAAGGTTACCGGTCAGCATCGCGAGGCTGGTCAGAGAACGCACGGTTTCCACGCCCTGGAAGAAATGGGTTACACCCATGCCCCACAGGATGGCGGCGCTTTTCGCCTGGGCATACATCCGTGCCGCCTGACGAATCTCACTGGCGCTGACGCCGGTGATATCTTCAACCGACTCCGGCGTGTAGCCTTCAACGATTTTACGATACTCTTCACAGCCTTCTGTACGTGAAGCGACGAACGCt   >  1:289205/1‑251 (MQ=255)                                                                                                                                                                                                                                             cGGCTTACCAAGGTTACCGGTCAGCATCGCGAGGCTGGTCAGAGAACGCACGGTTTCCACGCCCTGGTAGAACTGGGTTACACCCATGCCCCACAGGATGGCGGCGCTTTTCGCCTGGGCATACATCCGTGCCGCCTGACGAATCTCACTGGCGCTGACGCCGGTGATATCTTCAACCGACTCCGGCGTGTAGCCTTCAACGATTTTACGATACTCTTCAAAGCCTTCTGTACGTGAAGCGACGAACGCtt  <  1:231702/251‑1 (MQ=255)                                                                                                                                                                                                                                                      |                                                                                                                                                                                                                                                  TCTTCGCCCATAATGTACGCGGCACGCACTTCGCCATGCGCTGCGCGGTGCGGCAGCTCGCTGATGCGATAGCCGGTATGCGCTGGCAGGCTTTCCACGCCCCAGGCTTTGGCGAATTTCTCGCGGTTAGCCGGATCTTTCACGTACTGGTATCCCGGATACGTATCCGGCAGCGCGCCCATATCGCAGGCACCCTGAACGTTGTTCTGACCACGAACCGGGTTAACACCCGCATGCGGCTTACCGAGGTTACCGGTCAGCATCGCGAGGCTGGTCAGAGAACGCACGGTTTCCACGCCCTGGTAGAACTGGGTTACACCCATGCCCCACAGGATGGCGGCGCTTTTCGCCTGGGCATACATCCGTGCCGCCTGACGAATCTCACTGGCGCTGACGCCGGTGATATCTTCAACCGACTCCGGCGTGTAGCCTTCAACGATTTTACGATACTCTTCAAAGCCTTCTGTACGTGAAGCGACGAACGCTT  >  NC_000913/4298156‑4298642

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