Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	NC_000913	4,318,215	G→A	L184L (CTG→TTG)	phnK ←	carbon‑phosphorus lyase subunit PhnK

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	4,318,215	0	G	A	98.3%	160.0 / NA	59	L184L (CTG→TTG)	phnK	carbon‑phosphorus lyase subunit PhnK
Reads supporting (aligned to +/- strand):  ref base G (0/0);  major base A (43/15);  minor base C (1/0);  total (44/15)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.06e-01

AGGCCGGATAAGGCGTTTACGCCGCATCCGGCACCGGCTCAATTCTGCAAAACCGATGACACCAGCAGCTGTGTATACGGATGATGCGGGTCGTCGAGCACGCGGTCGGTTAACCCACTCTCCACCACTTGCCCCTGCTTCATCACCAGCAAACG‑GTCCGCCAGCAGGCGGGCGACGCCTAAATCATGGGTGACAATCACCACCG‑CGAGGTTCAGCTCCACCACCAGGCCGCGCAGCAGGTCGAGCAGGCGGGCCTGCACCGACACATCCAGCCCGCCGGTCGGTTCATCCATAAACACCAGCTTCGGATGCGTCACCAGGTTGCGGGCAATCTGCAAACGCTGCTGCATACCGCCGGAAAAGGTGGTCGGCAGGTCGTCGATCCGGTTGGCGGGAATCTCCACCTCTTCCAGCCACTTCTGCGCGGTGGCACGAATATCGCCGTAATGACGTGCCCCGGTCGCCATCAGCCGCTCGCCGATATTGCCGC  >  NC_000913/4317968‑4318457                                                                                                                                                                                                                                                          |                                                                                                                                                                                                                                                   aGGCCGGATCAGGAGTTTACGCCGCATCCGGCACCGGCTCAATTCTGCACAACCGATGACACCAGCAGCTGTGTATCCGGATGATGCGGGTCGTCGAGCACGCGGTCGGCTAACCCACTCTCCACCACTTGCCCCTGCTTCATCACCAGCAAACG‑GTCCGCCAGCCGGCGGGCGACGCCTCAATCATGGGTGACAATCACCACCG‑CGAGGTTCCGCTCCACCACCAGGCCGCGCACCAGGTCGAGCAAGCg                                                                                                                                                                                                                                                 >  2:134176/1‑251 (MQ=255)          aaGGCGTTTACGCAGCATCCGGCACCGGCTCAATTCTGCAAAACCGATGACACCAGCAGCTGTGTATACGGATGATGCGGGTCGTCGAGCACGCGGTCGGTTAACCCACTCTCCACCACTTGCCCCTGCTTCATCACCAGCAAACG‑GTCCGCCAGCAGGCGGGCGACGCCTAAATCATGGGGGCCAATCTCCACCG‑CGAGGTTCAGCTCCACCACCAGGCCGCGCAGCAGGTCGAGCAAGCGGGCCTGCAc                                                                                                                                                                                                                                        >  2:46356/1‑251 (MQ=255)                  tACGCCGCATCCGGCACCGGCTCAATTCTGCAAAACCGATGACACCAGCAGCTGTGTATACGGATGATGCGGGTCGTCGAGCACGCGGTCGGTTAACCCACTCTCCACCACTTGCCCCTGCTTCATCACCAGCAAACG‑GTCCGCCAGCAGGCGGGCGACGCCTAAATCATGGGTGCCAATCACCACCG‑CGAGGTTCAGCTCCACCACCAGGCCGCGCCGCAGGTCGAGCAAGCGGGCCTGCACCGAcacac                                                                                                                                                                                                                                >  2:89262/1‑250 (MQ=255)                     gccgcATCCGGCACCGGCTCAATTCTGCACAACCGATGACACCAGCAGCTGTGTATACGGATGATGCGGGTCGTCGAGCACGCGGTCGGTTAACCCACTCTCCACCACTTGCCCCTGCTTCATCACCAGCAAACG‑GTCCGCCAGCAGGCGGGCGACGCCTAAATCATGGGTGACAATCACCACCG‑CGCGGTTCAGCTCCACCACCAGGCCGCGCAGCAGGTCGAGCAAGCGGGCCTGCACCGACACATcc                                                                                                                                                                                                                              >  2:120956/1‑250 (MQ=255)                      ccgcATCCGGCACCGGCTCAATTCTGCAAAACCGATGACACCAGCAGCTGTGTATACGGATGATGCGGGTCGTCGAGCACGCGGTCGGTTAACCCACTCTCCACCACTTGCCCCTGCTTCATCACCAGCAAACG‑GTCCGCCAGCAGGCGGGCGACGCCTAAATCATGGGTGACAATCACCACCG‑CGAGGTTCAGCTCCACCACCAGGCCGCGCAGCAGGTCGCGCAAGCGGGCCTGCACCGACACATcc                                                                                                                                                                                                                              >  2:231807/1‑249 (MQ=255)                         cATCCGGCCCCGGCTCAATTCTGCAAAACCGATGACACCAGCAGCTGTGTATACGGATGATGCGGGTCGTCGAGCACGCGGTCGGTTAACCCACTCTCCACCACTTGCCCCTGCTTCATCACCAGCAAACG‑GTCCGCCAGCAGGCGGGCGACGCCTAAATCATGGGTGACAATCGCCACCG‑CGAGGTTCAGCTCCACCACCAGGCCGCGCAGCAGGTCGAGCAAGCGGGCCTGCACCGACACATCCAGccc                                                                                                                                                                                                                         <  1:257471/251‑1 (MQ=255)                               gCACCGGCTCAATTCTGCAAAACCGATGACACCAGCAGCTGTGTATACGGATGATGCGGGTCGTCGAGCACGCGGTCGGTTAACCCACTCTCCACCACTTGCCCCTGCTTCATCACCAGCAAACG‑GTCCGCCAGCAGGCGGGCGACGCCTAAATCATGGGTGACAATCACCACCG‑CGAGGTTCAGCTCCACCACCAGGCCGCGCAGCAGGTCGAGCAAGCGGGCCTGCACCGACACATCCAGCCCGCcggt                                                                                                                                                                                                                   >  1:159091/1‑251 (MQ=255)                                cACCGGCTCAATTCTGCAAAACCGATGACACCAGCAGCTGTGTATACGGATGATGCGGGTCGTCGAGCACGCGGTCGGTTAACCCACTCTCCACCACTCGCCCCTGCTTCATCACCAGCAAACG‑GTCCGCCAGCAGGCGGGCGACGCCTAAATCATGGGTGACAATCCCCACCG‑CGAGGTTCAGCTCCACCACCAGGCCGCGCAGCAGGTCGAGCAAGCGGGCCTGCACCGACACATCCAGCCCGCcggtc                                                                                                                                                                                                                  <  1:106348/251‑1 (MQ=255)                                 aCCGGCTCAATTCTGCAAACCCGATGACACCAGCAGCTGTGTATACGGATGATGCGGGTCGTCGAGCACGCGGTCGGTTAACCCCCTCTCCACCACTTGCCCCTGCTTCATCACCAGCAAACG‑GTCCGCCAGCAGGCGGGCGACGCCTAAATCATGGGTGACAACCACCACCG‑CGAGGTTCAGCTCCACCACCAGGCCGCGCAGCAGGTCGAGCAAGCGGGCCTGCACCGACACATCCAGCCCGCcggtcg                                                                                                                                                                                                                 <  2:31880/251‑1 (MQ=255)                                     gCTCAATTCTGCAAAACCGATGACACCAGCAGCTGTGTATACGGATGATGCGGGTCGTCGAGCACGCGGTCGGTTAACCCACTCTCCACCACTTGCCCCTGCTTCATCACCAGCAAACG‑GTCCGCCAGCAGGCGGGCGACGCCTAAATCATGGGTGACAATCACCACCG‑CGAGGGTCAGCTCCACCACCAGGCCGCGCAGCAGGTCGCGCAAGCGGGCCTGCACCGACACATCCAGCCCGCCGGTCGGTTc                                                                                                                                                                                                             >  1:46121/1‑251 (MQ=255)                                         aaTTCTGCAAAACCGATGACCCCAGCAGCTTTGTATCCGGATGATGCGGGTCGTCGAGCCCGCGGTCGGTTAACCCTCTCTCCCCCACTTGCCCCTGCTTCATCACCAGCAAACG‑GTCCGCCAGCAGGCGGGCGACGCCTAAATCATGGGTGACAATCCCCACCG‑CGAGGTTCAGCTCCACCACCAGGCCGCGCAGCAGGTCGAGCAAGCGGGCCTGCACCGACACATCCAGCCCGCCGGTCGGTTCATcc                                                                                                                                                                                                         <  1:185006/251‑1 (MQ=255)                                                 aaaaCAGATGACACCAGCAGCTGTGTATACGGATGATGCGGGACGTCGCGCACGCGGTCGGGTAACCCACTCTCCACCACTTGCCCCTGCTTCATCACCAGCAAACG‑GTCCGCCAGCAGGCGGGCGCCGCCTAACTCATGGGCGACAATCACCACCG‑CGAGGTTCAGCTCCACCACCAGGCCGCGCAGCGGGTCGAGCAAGCGGGCCTGCCCCGACACATCCAGCCCGCCGGTCGGTTCATccctcaaccc                                                                                                                                                                                                 >  2:72530/1‑243 (MQ=255)                                                     ccGATGACACCAGCAGCTGTGTATACGGATGATGCGGGTCGTCGAGCACGCGGTCGGTTAACCCACTCTCCACCACTTGCCCCTGCTTCATCACCAGCAAACG‑GTCCGCCAGCAGGCGGGCGACGCCTAAATCATGGGTGACAATCACCACCG‑CGAGGTTCAGCTCCACCACCAGGCCGCGCAGCAGGTCGAGCAAGCGGGCCTGCACCGACACATCCAGCCCGCCGGTCGGTTCATCCATAAACACCAGc                                                                                                                                                                                             >  1:171202/1‑251 (MQ=255)                                                      cGATGACACCAGCAGCTGGGTATCCGGAGGAGGCGGGTCGTCGAGCCCGCGGTCGGTTAACCCACTCTCCACCACTTGCCCCTTCTTCATCACCAGCAAACG‑GTCCGCCAGCAGGCGGGCGACGCCTAAATCATGGGTGACAATCACCACCG‑CGAGGTTCAGCTCCCCCCCCAGGCCGCGCAGCAGGTCGACCAAGCGGGCCTGCACCGACACATCCAGCCCGCCGGTCGGTTCATCCATAAACACCAGCt                                                                                                                                                                                            <  1:46356/251‑1 (MQ=255)                                                             cccAGCAGCTGTGTATGGGGATTATGCGGGTCGGCGCGCACCCGTTCGGTTAGCCCACTCTCCACCACTGGCCCTTGCTTCATCACCATCAACCG‑GTCTGCCGGCAGGCGGGCGACCCCTAAGTCATGGGGGTCAATGTCCACCT‑CGAGGTTCTTCTTCACCACGAGGCCTCGCAGCAGGTCGACTAAGCGGGCCGGCACGGTCTCATCCAGCCCGCCGGTCGGTTCATCCATAAACACCAGTTTCGGATg                                                                                                                                                                                     <  1:255912/250‑1 (MQ=255)                                                                  cagcTGTGTATACGGATGATGCGGGTCGTCGAGCACGCGGTCGGATAACCCACTCTCCACCACTTGCCCCTGCTTCATCACCAGCAAACG‑GTCCGCCAGCAGGCGGGCGACGCCTAAATCATGGGTGACAATCACCACCG‑CGAGGTTCAGCTCCACCACCATGCCGCGCAGCAGGTCGAGCAAGCGGGCCTGCACCGACACATCCAGCCCGCCGGTCGGTTCATCCCTAAACACCAGCTACGGATGCGTcc                                                                                                                                                                                >  1:80482/1‑250 (MQ=255)                                                                    gcTGTGTATACGGATGATGCGGGTCGTCGAGCACGCGGTCGGTTAACCCACTCTCCACCACTTGCCCCTGCTTCATCACCAGCAAACG‑GTCCGCCAGCAGGCGGGCGACGCCTAAATCATGGGTGACAATCACCACCG‑CGAGGTTCAGCTCCACCACCAGGCCGCGCAGCAGGTCGAGCAAGCGGGCCTGCACCGACACATCCAGCCCGCCGGTCGGTTCATCCATAAACACCAGCTTCGGATTCGTCAcc                                                                                                                                                                              >  1:4219/1‑251 (MQ=255)                                                                    gcTGTGTATACGGATGATGCGGGTCGTCGAGCACGCGGTCGGTTAACCCACTCTCCACCACTTGCCCCTGCTTCATCACCAGCAAACG‑GTCCGCCAGCAGGCGGGCGACGCCTAAATCATGGGTGACAATCACCACCG‑CGAGGTTCAGCTCCACCACCAGGCCGCGCAGCAGGTCGAGCAAGCGGGCCTGCACCGACACATCCAGCCCGCCGGTCGGTTCATCCATAAACACCAGCTTCGGATGCGTCAcc                                                                                                                                                                              >  1:4379/1‑251 (MQ=255)                                                                       gtgtATACGGATGATGCGGGTCGTCGAGCACGCGGTCGGTTAACCCACTCTCCACCACTTGCCCCTGCTTCATCACCAGCAAACG‑GTCCGCCAGCAGGCGGGCGACGCCTAAATCATGGGTGACAATCACCACCG‑CGAGGTTCAGCGCCACCACCAGGCCGCGCAGCAGGTCGAGCAAGCGGGCCTGCACCGACCCATCCTGCCCGCCGGTCGGTTCATCCATAAACACCAGCTTCGGATGAGTCACCAgg                                                                                                                                                                           >  2:142709/1‑251 (MQ=255)                                                                           ataCGGATGATGCGGGTCGTCGAGCACGCGGTCGGTTAACCCACTCTCCACCACTTGCCCCTGCTTCATCACCAGCAAACG‑GTCCGCCAGCAGGCGGGCGACGCCTAAATCATGGGTGACAATCACCACCG‑CGAGGTTCAGCTCCACCACCAGGCCGCGCAGCAGGTCGAGCAAGCGGGCCTGCACCGACACATCCAGCCCGCCGGTCGGTTCATCCATAAACACCAGCTTCGGATGCGTCACCAGGTTGc                                                                                                                                                                       >  1:91568/1‑251 (MQ=255)                                                                           agaCGGATGACGCGGGTCGTCGAGCACGCGGGCGGTTAACCCCCCCTCCCCCCCTTGCCCCTGCTTCAACCCCAGCAAACG‑GGCCGCCAGCAGGCGGGGGAGGCCAAAAGCATGGGTGACAATCACCCCCG‑CGGGGTCCAGCTCCCCCACCAGGCCGCTCAGCAGGTCGATCAAGCGGTCCTGCACCGCCACATCCAGCCCGCCGGTCGGTTCATCCATAAACACCAGCGTTGGATGCGGCACCATGTTGc                                                                                                                                                                       <  2:160510/249‑1 (MQ=255)                                                                             aCGGATGATGCGGGTCGTCGAGCACGCGGTCGGTTAACCCACTCTCCACCACTTGCCCCTGCTTCATCACCAGCAAACG‑GTCCGCCAGCAGGCGGGCGACGCCTAAATCATGGGTGACAATCACCACCG‑CGAGGTTCAGCTCCACCACCAGGCCGCGCAGCAGGTCGAGCAAGCGGGCCTGCACCGACACATCCAGCCCGCCGGTCGGTTCATCCATAAACACCAGCTTCGGATGCGTCACCAGGTTGCgg                                                                                                                                                                     >  1:29629/1‑251 (MQ=255)                                                                               ggATGATGCGGGTCGTCGAGCACGCGGTCGGTTAACCCACTCTCCCCCACTTGCCCCTGCTTCATCACCAGCAAACG‑GTCCGCCAGCAGGCGGGCGACGCCTAAATCATGGGTGACAATCACCACCG‑CGAGGTTCAGCTCCACCACCAGGCCGCGCAGCAGGTCGAGCAAGCGGGCCTGCACCGACACATCCAGCCCGCCGGTCGGTTCATCCATAAACACCAGCTTCGGAGGCGTCAc                                                                                                                                                                               <  1:212825/239‑1 (MQ=255)                                                                               ggATGATGCGGGTCGTCGAGCACGCGGTCGGTTAACCCACTCTCCACCACTTGCCCCTGCTTCATCACCAGCAAACG‑GTCCGCCAGCAGGCGGGCGACGCCTAAATCATGGGTGACAATCACCACCG‑CGAGGTTCAGCTCCACCACCAGGCCGCGCAGCAGGTCGAGCAAGCGGGCCTGCACCGACACATCCAGCCCGCCGGTCGGTTCATCCATAAACCCCAGCTTCGGATGCGTCACCAGGTTGCGGGc                                                                                                                                                                   >  2:44376/1‑251 (MQ=255)                                                                               ggATGATGCGGGTCGTCGAGCACGCGGTCGGTTAACCCACTCTCCACCACTTGCCCCTGCTTCATCACCAGCAAACG‑GTCCGCCAGCAGGCGGGCGACGCCTAAATCATGGGTGACAATCACCACCG‑CGAGGTTCAGCTCCACCACCAGGCCGCGCAGCAGGTCGAGCAAGCGGGCCTGCACCGACACATCCAGCCCGCCGGTCGGTTCATCCATAAAAACCAGCTTCGGATGCGTccc                                                                                                                                                                               >  2:212825/1‑237 (MQ=255)                                                                                   gatgCGGGTCGTCGAGCACGCGGTCGGTTAACCCACTCTCCACCACTTGCCCCTGCTTCATCACCAGCAAACG‑GTCCGCCAGCAGGCGGGCGACGCCTAAATCATGGGTGACAATCACCACCG‑CGAGGTTCAGCTCCACCACCAGGCCGCGCAGCAGGTCGAGCAAGCGGGCCTGCACCGACACATCCAGCCCGCCGGTCGGTTCATCCATAAACACCAGCTTCGGATGCGTCACCAGGTTGCGGGCAATc                                                                                                                                                               >  1:103988/1‑251 (MQ=255)                                                                                   gatgCGGGTCGTCGAGCACGCGGTCGGTTAACCCACTCTCCACCACTTGCCCCTGCTTCATCACCAGCAAACG‑GTCCGCCAGCAGGCGGGCGACGCCTAAATCATGGGTGACAATCACCACCG‑CGAGGTTCAGCTCCACCACCAGGCCGCGCAGCAGGTCGAGCAAGCGGGCCTGCACCGACACATCCAGCCCGCCGGTCGGTTCATCCATAAACACCAGCTTCGGATGCGTCACCAGGTTGCGGGCAATc                                                                                                                                                               >  1:145792/1‑251 (MQ=255)                                                                                      gCGGGTCGTCGAGCACGCGGTCGGTTAACCCACTCTCCACCACTTGCCCCTGCTTCATCACCAGCAAACG‑GTCCGCCAGCAGGCGGGCGACGCCTAAATCATGGGTGACAATCACCACCG‑CGAGGTTCAGCTCCACCACCAGGCCGCGCAGCAGGTCGAGCAAGCGGGCCTGCACCGACACATCCAGCCCGCCGGTCGGTTCATCCATAAACACCAGCTTCGGATGCGTCACCAGGTTGCGGGCAATCTGc                                                                                                                                                            <  1:120956/251‑1 (MQ=255)                                                                                        gggTCGTCGAGCACGCGGTCGGTTAACCCCCTCTCCACCACTTGCCCCTGCTTCATCACCAGCAAACG‑GTCCGCCAGCAGGCGGGCGACGCCTAAATCATGGGTGACAATCACCACCG‑CGAGGTTCAGCGCCACCACCAGGCCGCGCAGCAGGGCGAGCAAGCGGGCCTGCACCGACACATCCAGCCCGCCGGTCGGTTCATCCATAAACACCAGCTTCGGATGCGTCACCAGGTTGCGGGCAATCTGCaa                                                                                                                                                          <  1:225104/251‑1 (MQ=255)                                                                                            cgtcgAGCACGCGGTCAAGTGCCCCACTCTCCACCACTTGCCCCTGCTTCATCACCAGCAAACG‑GTCCGCCAGCCGGCGGGCGACGCCTCAATCATGGGTGACAATCACCACCG‑CGAGGTTCAGCTCCACCACCAGGCCGCGCAGCAGGCCGAGCAAGCGGGCCTGCACCGACACATCCAGCCCGCCGGGCGGTTCCTCCATAAACACCCGCTTCGGATGCGTCACCAGGTTGCGGGCAATCTGCAAACgc                                                                                                                                                      >  2:119228/1‑251 (MQ=255)                                                                                                      gcgGTCGGTTAACCCACTCTCCACCACTTGCCCCTGCTTCATCACCAGCAAACG‑GTCCGCCAGCAGGCGGGCGACGCCTAAATCATGGGTGACAATCACCACCG‑CGAGGTTCAGCTCCACCACCAGGCCGCGCAGCAGGTCGAGCAAGCGGGCCTGCACCGACACATCCAGCCCGCCGGTCGGTTCATCCATAAACACCAGCTTCGGATGCGTCACCAGGTTGCGGGCAATCTGCAAACGCTGCTGCATAc                                                                                                                                            >  1:13213/1‑251 (MQ=255)                                                                                                               tAACCCACTCTCCACCACTTGCCCCTGCTTCATCACCAGCAAACG‑GTCCGCCAGCAGGCGGGCGACGCCTAAATCATGGGTGACAATCACCACCG‑CGAGGTTCAGCTCCACCACCAGGCCGCGCAGCAGGTCGAGCAAGCGGGCCTGCACCGACCCATCCAGCCCGCCGGTCGGTTCATCCATAAACACCAGCTTCGGATGCGTCACCAGGTTGCGGGCAATCTGCAAACGCTGCTGCATACCGCCggca                                                                                                                                    >  2:76385/1‑248 (MQ=255)                                                                                                                   ccACTCTCCACCACTTGCCCCTGCTTCATCACCAGCAAACG‑GTCCGCCAGCAGGCGGGCGACGCCTAAATCATGGGCGACACTCACCACCG‑CGAGGTTCAGCTCCACCACCAGGCCGCGCAGCAGGTCGAGCAAGCGGGCCTGCACCGACACATCCAGCCCGCCGGTCGGTTCATCCATAAACACCAGCTTCGGATGCGTCACCAGGTTGCGGGCAATCTGAAAACGCTGATGCATACCGCCGGAAAAggg                                                                                                                               >  1:252022/1‑250 (MQ=255)                                                                                                                    cACTCTCCACCACTTGCCCCTGCTTCATCACCAGCAAACG‑GTCCGCCAGCAGGCGGGCGACGCCTAAATCAAGGGTGACAATCACCACCG‑CGAGGTTCAGCTCCACCACCAGGCCGCGCAGCCGGTCGAGCAAGCGGGCCTGCACCGACACATCCAGCCCGCCGGTCGGTTCATCCATAAACACCAGCTTCGGATGCGTCACCAGGTTGCGGGCAATCTGCAAACGCTGCTGCATACCGCCGGAAAAggtg                                                                                                                              >  1:35235/1‑251 (MQ=255)                                                                                                                      ctctCCACCACTTGCCCCTGCTTCATCAACAGCACACG‑GTCCGCCAGCAGGCGGGCGCCGCCTAAATCATGGGTGACAATCACCACCG‑AGAGGTTCCGCTCCCCCACCAGGCCGCGCAGCAGGTCGAGCCAGCGGGCCTGCACCGACACCTCCCGCCCGCCGGTCGGTTCATCCATAAACACCAGCTTCGGATGCGTCACCAGGGTGCGGGACATCTGCAAACGCTGCTGCATACCGCCGGAAAAggtggt                                                                                                                            >  2:54336/1‑251 (MQ=255)                                                                                                                                ctggCCCCTGCTTCCGCACCAGCAA‑CGTGTCCGCCAGCAGGCGGCCGCTGCTTAAATCGTGGGTGACA‑GCCCCACCGGCGGGGTTCAGCTCCCCCCCCAGGCCGCGCAGCATGTCGGGCAAGCGGGCCTGCACCGACACATCCCGCCCGCCGGTCGGTTCATCCATAAACACCAGCTTCGGATGCGTCACCAGGTTGCGGGCAAGCTGCAAACGCTGCTGCATACCGCCGGAAAAGGTGGTCGGCAGgt                                                                                                                    <  2:57285/246‑1 (MQ=255)                                                                                                                                       cggCTTCATCACAAGCAAACG‑GTCCGCCAGCAGGCGGGCGACGCCTAAATCATGGGTGACAATCACCACCG‑CGAGGTTCAGCTCCACCACCAGGCCGCGCAGCAGGTCGAGCAAGCGGGCCTGCACCGACACATCCAGCCCGCCGGTCGGTTCATCCATAAACACCAGCGTCGGATGCGTCACCAGGTTGCGGGCAATCTGCAAACGCTGCTGCATACCGCCGGAAAAGGTGGTCGGCAGGTCGTCGATcc                                                                                                           <  1:89262/249‑1 (MQ=255)                                                                                                                                          cTTCATCACCAGCAAACG‑GTCCGCCAGCAGGCGGGCGACGCCTAAATCATGGGTGACAATCACCACCG‑CGAGGTTCAGCTCCACCACCAGGCCGCGCAGCCGGTCGAGCAAGCGGGCCTGCACCGACACATCCAGCCCGCCGGTCGGTTCATCCATAAACACCAGCTTCGGATGCGTCACCAGGTTGCGGGCAATCTGCAAACGCTGCTGCATACCGCCGGAAAAGGTGGTCGGCAGGTCGTCGATCCGGt                                                                                                        >  2:137186/1‑251 (MQ=255)                                                                                                                                                    aGCAAACG‑GTCCGCCAGCAGGCGGGCGACGCCTAACTCCTGGGTGACAATCCCCCCCG‑CGAGGATCAGCTCCACCACCAGGCCGCGCAGCAGGTCGAGCAAGCGGGCCTGCACCGACACATCCAGCCCGCCGGTCGGTTCATCCCTAAAACCCAGCTTCGGCTGCGTCACCAGGTTGCGGGCAATCTGCAAACGCTGCTGCATACCGCCGGACAAGGTGGTCGGCAGGTCGTCGATCCGGTTGGCGGGAAt                                                                                              >  2:33511/1‑251 (MQ=255)                                                                                                                                                         ccG‑GTCCGCCAGCAGGCGGGCGACGCCTAAATCATGGGGGACAATCCCCACCG‑CGAGGTTCAGCCCCACCACCAGGCCGCGCAGCAGGTCGAGCAAGCGGGCCTGCACCGACACATCCAGCCCGCCGGTCGGTTCATCCATAAACACCAGCTTCGGATGCGTCACCAGGTTGCGGGCAATCTGCAAACGCTGCTGCATACCGCCGGAAAAGGTGGTCGGCAGGTCGTCGATCCGGTTGGCGGGAATCTCCa                                                                                         <  1:119228/250‑1 (MQ=255)                                                                                                                                                              tCCGCCAGCAGGCGGGCGACGCCTAAATCATGGGTGACAATCACCACCG‑CGAGGTTCAGCTCCACCACCAGGCCGCGCAGCAGGTCGAGCAAGCGGGCCTGCACCGACACATCCAGCCCGCCGGTCGGTTCATCCATAAACACCAGCTTCGGATGCGTCACCAGGTTGCGGGCAATCTGCAAACGCTGCTGCATACCGCCGGAAAAGGTGGTCGGCAGGTCGTCGATCCGGTTGGCGGGAATCTCCACctc                                                                                     >  1:235829/1‑251 (MQ=255)                                                                                                                                                               ccGCCAGCAGGCGGGCGACGCCTAAATCATGGGTGACAATCCCCACCG‑CGAGGTTCAGCTCCACCACCAGGCCGCGCAGCAGGTCGAGCAAGCGGGCCTGCACCGACCCACCCAGCCCGCCGGTCGGTTCATCCATAAACACCAGCTTCGGATGCGTCACCAGGTTGCGGGCAATCTGCAAACGCTGCTGCATACCGCCGGAAAAGGTGGTCGGCAGGTCGTCGATCCGGTTGGCGGGAATCTCCACctct                                                                                    <  2:103988/251‑1 (MQ=255)                                                                                                                                                                 gCCAGCAGGCGGGCGACGCCTCAATCATGGGTGACAATCACCACCG‑CGAGGTTCAGCTCCACCACCAGGCCGCGCAGCAGGTCGAGCAAGCGGGCCTGCACCGACACATCCAGCCCGCCGGTCGGTTCATCCATAAACACCAGCTTCGGATGCGTCACCAGGTTGCGGGCAATCTGCAAACGCTGCTGCATACCGCCGGAAAAGGTGGTCGGCAGGTCGTCGATCCGGTTGGCGGGAATCTCCACCTCTTc                                                                                  >  1:58375/1‑251 (MQ=255)                                                                                                                                                                           gggcgACGCCTAAATCATGAGTGACAATCACCACCG‑CGAGGTTCAGCTCCCCCACCAGGCCGAGCAGCAGGTCGAGCACGCGGGCCTGCACCGACACATCCAGCCCGCCGGTCGGTTCCTCCATAAACCCCAGCTTCGGATGCGTCACCAGGTTGCGGGCCATCTGCAAACGCTGCTGCATACCGCCGGAAAAGGTGGCCGGCAGGTCGTCGATCCGGTTGGCGGGAATCTCCACCGCTTTCAGCCACTTc                                                                        >  2:27261/1‑251 (MQ=255)                                                                                                                                                                                  gCCTAAATCATGGGTGACAATCACCACCG‑CGAGGTTCAGCTCCACCACCAGGCCGCGCAGCAGGTCGAGCAAGCGGGCCTGCACCGACACATCCAGCCCGCCGGTCGGTTCATCCATAAACACCAGCTTCGGATGCGTCACCAGGTTGCGGGCAATCTGCAAACGCTGCTGCATACCGCCGGAAAAGGTGGTCGGCAGGTCGTCGATCCGGTTGGCGGGAATCTCCACCTCTTCCAGCCACTTCTGCGCgg                                                                 >  2:11453/1‑251 (MQ=255)                                                                                                                                                                                    cTAAATCATGGGTGACAATCACCACCG‑CGAGGTTCAGCTCCACCACCAGGCCGCGCAGCAGGTCGAGCAAGCGGGACTGCACCGACACATCCAGCCCGCCGGACGGTTCATCCATAAACACCAGCTTCGGATGCGTCACCAGGTTGCGGGCAATCTGCAAACGCTGCTGCATACCGCCGGAAAAGGTGGTCGGCAGGTCGTCGATCCGGTTGCCGGGAATCTCCACCTCTTCCACCCACTTCTGCGCGGTg                                                               >  2:235706/1‑251 (MQ=255)                                                                                                                                                                                                   cactCACCACCG‑CGAGGTTCAGCTCCACCACCAGGCCGCGCAGCAGGTCGAGCAAGCGGGCCTGCACCGACACATCCAGCCCGCCGGTCGGTTCATCCATAAACACCAGCTTCGGATGCGTCACCAGGTTGCGGGCAATCTGCAAACGCTGCTGCATACCGCCGGAAAAGGTGGTCGGCAGGTCGTCGATCCGGTTGGCGGGAATCTCCACCTCTTCCAGCCACTTCTGCGCGGGGGACCGAATCTCGCCg                                                >  2:81222/4‑251 (MQ=255)                                                                                                                                                                                                    aaTCACCACCG‑CGAGGTTCAGCTCCACCACCAGGCCGCGCAGCAGGTCGAGCAAGCGGGCCTGCACCGACACATCCAGCCCGCCGGTCGGTTCATCCATAAACACCAGCTTCGGATGCGTCACCAGGTTGCGGGCAATCTGCAAACGCTGCTGCATACCGCCGGAAAAGGTGGTCGGCAGGTCGTCGATCCGGTTGGCGGGAATCTCCACCTCTTCCAGCCACTTCTGCGCGGTGGCACGAATATCGCCgg                                               >  2:115251/1‑250 (MQ=255)                                                                                                                                                                                                          caccG‑CGAGGTTCAGCTCCACCACCAGGCCGCGCAGCAGGTCGAGCAAGCGGGCCTGCACCGACACATCCAGCCCGCCGGTCGGTTCATCCATAAACACCAGCTTCGGATGCGTCACCAGGTTGCGGGCAATCTGCAAACGCTGCTGCATACCGCCGGAAAAGGTGGTCGGCAGGTCGTCGATCCGGTTGGCGGGAATCTCCACCTCTTCCAGCCACTTCTGCGCGGTGGCACGAATATCGCCGTAATGAc                                         >  1:92430/1‑251 (MQ=255)                                                                                                                                                                                                           accG‑CGAGGTTCAGCTCCACCACCAGGCCGCGCAGCAGGTCGAGCAAGCGGGCCTGCACCGACACATCCAGCCCGCCGGTCGGTTCATCCATAAACACCAGCTTCGGATGCGTCACCAGGTTGCGGGCAATCTGCAAACGCTGCTGCATACCGCCGGAAAAGGTGGTCGGCAGGTCGTCGATCCGGTTGGCGGGAATCTCCACCTCTTCCAGCCACTTCTGCGCGGTGGCACGAATATCGCCGTAATGACg                                        >  2:216320/1‑251 (MQ=255)                                                                                                                                                                                                                cgAGGTTCAGCTCCACCACCAGGCCGCGCAGCAGGTCGAGCAAGCGGGCCTGCACCGACACATCCAGCCCGCCGGTCGGTTCATCCATAACCACCAGCTTCGGATGCGTCACCAGGTTGCGGGCAATCTGCAAACGCTGCTGCATACCGCCGGAAAAGGTGGTCGGCAGGTCGTCGATCCGGTTGGCGGGACTCTCCACCTCTTCCAGCCCCTTCTGCGCGGTGGCCCGAATATCGCCGTAATgcagtgcc                                    >  2:81182/1‑244 (MQ=255)                                                                                                                                                                                                                   gctccagccccccccccggcccGCGCAGCAGGTCGAGCAAGCGGGCCTGCCCCGACACATCCGGCCCGCCGGTCGGTTCATCCATAAACCCCAGCTTCGGAGGCGTCACCAGGTTGCGGGCAATCTGCAAACGCTGCTGCATACCGCCGGAAAAGGTGGTCGGCAGGTCGTCGATCCGGTTGGCGGGAATCTCCACCTCTTCCAGCCACTTCTGCGCGGTGGCACGAATATCGCCGTAATGACGTGCCCCg                                 <  2:48619/231‑1 (MQ=255)                                                                                                                                                                                                                        ggCTCCACCACCAGGCCGCGCAGAAGGTCGAGCAAGCGGGCCTGCACCGACACATCCAGCCCGCCGGTCGGTTCATCCATAAACACCAGCTTCGGATGCGTCACCAGGTTGCGGGCAATCTGCAAACGCTGCTGCATACCGCCGGAAACGGTGGTCGGCAGGTCGTCGCTCCGGTTGGCGGGAATCCCCCCCTCTTCCAGCCACTTCTGCGCGGTGGCACGCATATCGCCGTAATGACGTGCCCCGGTCGc                            >  2:99562/2‑251 (MQ=255)                                                                                                                                                                                                                        aGCTCCACCACCAGGCCGCGAAGCAGGTCGAGCAAGCGGGCCTGCACCGACACATCCAGCCCGCCGGTCGGTTCATCCATAAACACCAGCTTCGGATGCGTCACCAGGTTGCGGGCAATCTGCAACCGCTGCTGCATACCGCCGGAAAAGGTGGTCGGCAGGTCGTCGATCCGGTTGGCGGGAATCTCCACCTCTTCCAGCCACTTCTGCGCGGTGGCACGAATATCGCCGTAATGACGTGCCCCGGTCGc                            >  2:268596/1‑251 (MQ=255)                                                                                                                                                                                                                                caccaGGCCGCGCAGCAGGTCGAGCAAGCGGGCCTGCACCGACACATCCCGCCCGCCGGTCGGTTCATCCATAAACACCAGATTCGGATGCGTCACCAGGTTGCGGGCAATCTGCAAACGCTGCTGCATACCGCCGGAAAAGGTGGTCGGCAGGTCGTCGATCCGGTTGGCGGGAATCTCCACCTCTTCCAGCCACTTCTGCGCGGTGGCACGAATATCGCCGTAATGACCTGCCCCGGTCGCCATCAGcc                    >  1:222080/1‑251 (MQ=255)                                                                                                                                                                                                                                       ccGCGCAGCAGGTCGAGCAAGCGGGCCTGCACCGACACATCCAGCCCGCCGGTCGGTTCATCCATAAACACCAGCTTCGGATGCGTCACCAGGTTGCGGGCAATATGCAAACGCTGCTGCATACCGCCGGAAAAGGTGGTCGGCAGGTCGTCGATCCGGTTGGCGGGAATCTCCACCTCTTCCAGCCACTTCTGCGCGGTGGCCCGAATATCGCCGTAATGACGTGCCCCGGTCGCCCTCAGCCGCTCGcc             >  2:220209/1‑251 (MQ=255)                                                                                                                                                                                                                                            cagcagGTCGAGCAAGCGGGCCTGCACCGACACATCCAGCCCGCCGGTCGGTTCATCCATAAACACCAGCTTCGGATGCGTCACCAGGTTGCGGGCAATATGCAAACGCTGCTGCATACCGCCGGAACAGGCGGTCGGCAGGTCGTCGATCCGGTTGGCGGGAATCTCCACCTCTTCCAGCCACTTCTGCGCGGTGGCACGCATATCGCCGTAATGACGTGCCCCGGTCGCCTTAAGCCGCTCGCCGatat        >  2:24633/1‑251 (MQ=255)                                                                                                                                                                                                                                                 ggTCGAGCAAGCGGGCCTGCACCGACACATCCAGCCCGCCGGTCGGTTCATCCATAAACACCAGCTTCGGATGCGTCACCAGGTTGCGGGCAATCTGCAAACGCTGCTGCATACCGCCGGAAAAGGTGGTCGGCAGGTCGTCGATCCGGTTGGCGGGAATCTCCACCTCTTCCAGCCACTTCTGCGCGGTGGCACGAATATCGCCGTAATGACGTGCCCCGGTCGCCATCAGCCGCTCGCCGATATTgccg   >  1:262011/1‑251 (MQ=255)                                                                                                                                                                                                                                                  gTCGAGCAAGCGGGCATGCACCGACACCTCAAGCATGACGGTCGGATCATCCATAAACGCCAGCTTCGGATACGTCACCAGGATGCCGGCAATCTGCAAACGATGCTGCATACCGCCGGAAAAGGGGATCGGCAGGACGTCGATCCGGTTGGCGGGAATCTCCACCTCTTCCAGCCACTTCTGCGCGGTAGCACGAATATCGCCGTAATGACGTGCCCCGGTCGCCATCAGCCGCTCTCCGATATTgccgc  >  2:117357/1‑251 (MQ=255)                                                                                                                                                                                                                                                          |                                                                                                                                                                                                                                                   AGGCCGGATAAGGCGTTTACGCCGCATCCGGCACCGGCTCAATTCTGCAAAACCGATGACACCAGCAGCTGTGTATACGGATGATGCGGGTCGTCGAGCACGCGGTCGGTTAACCCACTCTCCACCACTTGCCCCTGCTTCATCACCAGCAAACG‑GTCCGCCAGCAGGCGGGCGACGCCTAAATCATGGGTGACAATCACCACCG‑CGAGGTTCAGCTCCACCACCAGGCCGCGCAGCAGGTCGAGCAGGCGGGCCTGCACCGACACATCCAGCCCGCCGGTCGGTTCATCCATAAACACCAGCTTCGGATGCGTCACCAGGTTGCGGGCAATCTGCAAACGCTGCTGCATACCGCCGGAAAAGGTGGTCGGCAGGTCGTCGATCCGGTTGGCGGGAATCTCCACCTCTTCCAGCCACTTCTGCGCGGTGGCACGAATATCGCCGTAATGACGTGCCCCGGTCGCCATCAGCCGCTCGCCGATATTGCCGC  >  NC_000913/4317968‑4318457

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