Predicted mutation
evidence	seq id	position	mutation	freq	annotation	gene	description
RA	NC_000913	542,105	T→G	21.2%	V73G (GTG→GGG)	glxK →	glycerate kinase II

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	542,105	0	T	G	21.2%	115.3 / 2.5	68	V73G (GTG→GGG)	glxK	glycerate kinase II
Reads supporting (aligned to +/- strand):  ref base T (15/35);  new base G (13/4);  total (29/39)
Fisher's exact test for biased strand distribution p-value = 1.36e-03
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

TGAAAAAGGATTGATAATGAAGATTGTCATTGCGCCAGACTCTTTTAAAGAGAGCTTAAGTGCAGAAAAATGTTGTCAGGCAATTAAAGCCGGGTTTTCGACCCTCTTTCCCGATGCGAACTATATCTGTTTGCCGATAGCGGATGGCGGCGAAGGGACGGTGGATGCGATGGTCGCCGCGACGGGCGGCAACATCGTGACGC‑TTGAAGTCTGCGGGCCGATGGGCGAAAAAGTGAATGCTTTTTATGGCCTTACCGGCGACGGGAAAACGGCGGTGATTGAGATGGCGGCAGCAAGTGGCCTGATGCTGGTCGCGCCTGAAAAGCGTAATCCGTTGCTGGCCTCCAGT‑TTTGGTACGGGGGAGTTAATTCGTCATGCGCTGGATAACGACATTCGCCATATTATTCTCGGCATTGGCGGCAGTGCGACGGTCGACGGCGGTATGGGCATGGCGCAGGCGCT  >  NC_000913/541872‑542333                                                                                                                                                                                                                                           |                                                                                                                                                                                                                                      tgAAAAAGGATTGATAATGAAGATTGTCATTGCGCCAGACTCTTTTAAAGAGAGCTTAAGTGCAGAAAAATGTTGTCAGGCAATTAAAGCCGGGTTTTCGACCCTCTTTCCCGATGCGAACTATATCTGTTTGCCGATAGCGGATGGCGGCGAAGGGACGGTGGATGCGATGGTCGCCGCGACGGGCGGCAACATCGTGACGC‑TTGAAGACTGCGGGCCGATGGGGGAAAAAGGGCATg                                                                                                                                                                                                                                  >  1:707832/1‑239 (MQ=255)  gAAAAAGGATTGATAAGGAAAATTGTCATTGCGCCAGACTCTTTTAAAGAGAGCTTAAGTGCAGAAAAATGTTGTCAGGCAATTAAAGCCGGGTTTTCGACCCTCTTTCCCGATGCGAACTATATCTGTTTGCCGATAGCGGATGGCGGCGAAGGGACGGTGGATGCGATGGTCGCCGCGACGGGCGGCAACATCGTGACGC‑TTGAAGTCTGCGGGCCGATGGGCGAAAAAGTGAATGc                                                                                                                                                                                                                                 <  1:725005/239‑1 (MQ=255)     aaaGGATTGATAATGAAGATTGTCATTGCGCCAGACTCTTTTAAAGAGAGCTTAAGTGCAGAAAAATGTTGTCAGGCAATTAAAGCCGGGTTTTCGACCCTCTTTCCCGATGCGAACTATATCTGTTTGCCGATAGCGGATGGCGGCGAAGGGACGGTGGATGCGATGGTCGCCGCGACGGGCGGCACAATCGGGCGGC‑TTGAAGCTGGCGGGCCGATGGGGGAAAAAGAGAAAGGttt                                                                                                                                                                                                                              >  2:830268/1‑239 (MQ=255)          aTTGATAATGAAGATTGTCATTGCGCCAGACTCTTTTAAAGAGAGCTTAAGTGCAGAAAAATGTTGTCAGGCAATTAAAGCCGGGTTTTCGACCCTCTTTCCCGATGCGAACTATATCTGTTTGCCGATAGCGGATGGCGGCGAAGGGACGGTGGATGCGATGGTCGCCGCGACGGGCGGCAACATCGTGACGC‑TTGAAGTCTGCGGGCCGAGGGGGGAAAAAGGGAATGCTTTTTaag                                                                                                                                                                                                                         >  2:87537/1‑237 (MQ=255)                   gAAGATTGTCATTGCGCCAGACTCTTTTAAAGAGAGCTTAAGTGCAGAAAAATGTTGTCAGGCAATTAAAGCCGGGTTTTCGACCCTCTTTCCCGATGCGAACTATATCTGTTTGCCGATAGCGGATGGCGGCGAAGGGACGGTGGATGCGATGGTCGCCGCGACGGGCGGCAACATCGTGACGC‑TTGAAGTCTGCGGGCCGATGGGCGAAAAAGTGAATGCTTTTTATGGCCTTACCg                                                                                                                                                                                                                >  2:50066/1‑239 (MQ=255)                    aaGATTGTCATTGCGCCAGACTCTTTTAAAGAGAGCTTAAGTGCAGAAAAATGTTGTCAGGCAATTAAAGCCGGGTTTTCAACCCTCTTTCCCGATGCGAACTATATCTGTTTGCCGATAGCGGATGGCGGCGAAGGGACGGGGGATGCGATGGTCGCCGCGACGGGCGGCACCCTCGTGACGC‑TTGGAGTATGCGGGCCGATGGGCACAAAAGGGAAGGCTTTTTATGGCCTCACCgg                                                                                                                                                                                                               >  2:415448/1‑239 (MQ=255)                     aGATTTTCTTTGCGCCAGCTTTTTTTAAAGGGGGGTTAAGTGGAAAAAAATGTTGTCAGGCAATTAAAGCCGGGTTTTCGACCCTCTTTCCCGATGCGAACTATATCTGTTTGCCGATAGCGGATGGCGGCGAAGGGACGGTGGATGCGATGGTCGCCGCGACGGGCGGCAACATCGTGACGC‑TTGAAGTCTGCGGGCCGATGGGCGAAAAAGTGAATGCTTTTTATGGCCTTACCGGc                                                                                                                                                                                                              <  2:1039855/239‑1 (MQ=255)                           tCATTGCGCCAGACTCTTTTAAAGAGAGCTTAAGTGCAGAAAAATGTTGTCAGGCAATTAAAGCCGGGTTTTCGACCCTCTTTCCCGATGCGAACTATATCTGTTTGCCGATAGCGGATGGCGGCGAAGGGACGGTGGATGCGATGGTCGCCGCGACGGGCGGCAACATCGTGACGC‑TTGAAGTCTGCGGGCCGATGGGGGAAAAAGTGAATGCTTTTTTTGGCCTTCCCGGGGGCggg                                                                                                                                                                                                        >  2:721002/1‑239 (MQ=255)                             aTTGCGCCAGACTCTTTTAAAGAGAGCTTAAGTGCAGAAAAATGTTGTCAGGCAATTAAAGCCGGGTTTTCGACCCTCTTTCCCGATGCGAACTATATCTGTTTGCCGATAGCGGATGGCGGCGAAGGGACGGCGGATGCGATGGGTCCCGCGCCGGGGGGCCAACTCGTAACGC‑TTGAAGCTTGCGGGCCGCGAGGCGAGAACGGGAAAGCtctttcggcacttccccggaccgggaa                                                                                                                                                                                                      >  2:871585/1‑217 (MQ=255)                                           ttttAAAGAGAGCTTAAGTGCAGAAAAATGTTGTCAGGCAATTAAAGCCGGGTTTTCGACCCTCTTTCCCGATGCGAACTATATCTGTTTGCCGATAGCGGATGGCGGCGAAGGGACGGTGGATGCGATGGTCGCCGCGACGGGCGGCAACATCGTGACGC‑TTGAAGTCTGCGGGCCGATGGGCGAAAAAGGGAATGCTTTTTATGGCCTTACCGGCGACGGGAAAACGGCGGTGAttt                                                                                                                                                                                        >  1:780243/1‑238 (MQ=255)                                           ttttAAAGAGAGCTTAAGTGCAGAAAAATGTTGTCAGGCAATTAAAGCCGGGTTTTCGACCCTCTTTCCCGATGCGAACTATATCTGTTTGCCGATAGCGGATGGCGGCGAAGGGACGGTGGATGCGATGGTCGCCGCGACGGGCGGCAACATCGTGACGC‑TTGAAGTCTGCGGGCCGATGGGCGAAAAAGGGAAGGCTTTTTCTGGCCTTACCGTCGACGGGaaacagggggggtggg                                                                                                                                                                                        >  2:339733/1‑226 (MQ=255)                                                     agCTTAAGTGCAGAAAAATGTTGTCAGGCAATTAAAGCCGGGTTTTCGACCCTCTTTCCCGATGCGAACTATATCTGTTTGCCGATAGCGGATGGCGGCGAAGGGACGGTGGATGCGATGGTCGCCGCGACGGGCGGGCACATCGTGACGC‑TTGAAGTCTGCGGGCCGGTGGGCGCAAAAGGGAATGCTTTTTATGGCCTTACCGGCGACGGGGAAAcggcgggggtggggtaggtggc                                                                                                                                                                              >  2:323586/1‑223 (MQ=255)                                                     agCTTAAGTGCAGAAAAATGTTGTCAGGCAATTAAAGCCGGGTTTTCGACCCTCTTTCCCGATGCGAACTATATCTGTTTGCCGATAGCGGATGGCGGCGAAGGGACGGTGGATGCGATGGTCGCCGCGACGGGCGGCAACATCGTGCCGC‑TTGAAGTCTGCGGGCCGATGGGCGAAAAAGGGAATGCTTTTTATGGCCTTACCGGCGACGGGAAAACGGCGGGGTATGAGATggcggc                                                                                                                                                                              >  1:743019/1‑239 (MQ=255)                                                     agCTTAAGTGCAGAAAAATGTTGTCAGGCAATTAAAGCCGGGTTTTCGACCCTCTTTCCCGATGCGAACTATATCTGTTTGCCGATAGCGGATGGCGGCGAAGGGACGGTGGATGCGATGGTCGCCGCGACGGGCGGCAACATCGTGACGC‑TTGAAGTCTGGGGGCCGATGGGCGAAAAAGTGAATGCTTTTTATGGCCTTACCGGCGACGGGAAAAAGGGGGGGAATGGGATgggggg                                                                                                                                                                              >  2:343613/1‑238 (MQ=255)                                                     agCTTAAGTGCAGAAAAATGTTGTCAGGCAATTAAAGCCGGGTTTTCGACCCTCTTTCCCGATGCGAACTATATCTGTTTGCCGATAGCGGATGGCGGCGAAGGGACGGTGGATGCGATGGTCGCCGCGACGGGCGGCAACATCGTGACGC‑TTGAAGTCTGCGGGCCGATGGGCGAAAAAGTGAATGCTTTTTATGGCCTTACCGGCGCCGGGAAAACGGCGGGGATTGGGATgggggc                                                                                                                                                                              >  1:164288/1‑239 (MQ=255)                                                     agCTTAAGTGCAGAAAAATGTTGTCAGGCAATTAAAGCCGGGTTTTCGACCCTCTTTCCCGATGCGAACTATATCTGTTTGCCGATAGCGGATGGCGGCGAAGGGACGGTGGATGCGATGGTCGCCGCGACGGGCGGCAACATCGTGACGC‑TTGAAGTCTGCGGGCCGATGGGCGAAAAAGTGAATGCTTTTTATGGCCTTACCGGCGACGGGAAAAcggcg                                                                                                                                                                                               <  1:62428/222‑1 (MQ=255)                                                     agCTTAAGTGCAGAAAAATGTTGTCAGGCAATTAAAGCCGGGTTTTCGACCCTCTTTCCCGATGCGAACTATATCTGTTTGCCGATAGCGGATGGCGGCGAAGGGACGGTGGATGCGATGGTCGCCGCGACGGGCGGCAACATCGTGACGC‑TTGAAGTCTGCGGGCCGATGGGCGAAAAAGTGAATGCTTTTTATGGCCTTACCGGCGACGGGAAAAcggcg                                                                                                                                                                                               >  2:62428/1‑222 (MQ=255)                                                     agCTTAAGTGCAGAAAAATGTTGTCAGGCAATTAAAGCCGGGTTTTCGACCCTCTTTCCCGATGCGAACTATATCTGTTTGCCGATAGCGGATGGCGGCGAAGGGACGGTGGATGCGATGGTCGCCGCGACGGGCGGCAACATCGTGACGC‑TTGAAGTCTGCGGGCCGATGGGCGAAAAAGTGAATGCTTTTTATGGCCTTACCGGCGACGGGAAAACGGCGGAGATTGGGATgtgtgc                                                                                                                                                                              >  1:114999/1‑234 (MQ=255)                                                     agCTTAAGTGCAGAAAAATGTTGTCAGGCAATTAAAGCCGGGTTTTCGACCCTCTTTCCCGATGCGAACTATATCTGTTTGCCGATAGCGGATGGCGGCGAAGGGACGGTGGATGCGATGGTCGCCGCGACGGGCGGCAACATCGTGACGC‑TTGAAGTCTGCGGGCCGATGGGCGAAAAAGTGAATGCTTTTTATGGCCTTACCGGCGACGGGAAAACGGCGGAGATTGAGATggcggc                                                                                                                                                                              >  1:19083/1‑239 (MQ=255)                                                     agCTTAAGTGCAGAAAAATGTTGTCAGGCAATTAAAGCCGGGTTTTCGACCCTCTTTCCCGATGCGAACTATATCTGTTTGCCGATAGCGGATGGCGGCGAAGGGACGGTGGATGCGATGGTCGCCGCGACGGGCGGCAACATCGTGACGC‑TTGAAGTCTGCGGGCCGATGGGCGAAAAAGGGAATGCTTTTTATGGCCTTACCGGCGCCGGGCAAACGGCGGGGATTGAGTTgggggc                                                                                                                                                                              >  1:540061/1‑239 (MQ=255)                                                     agCTTAAGTGCAGAAAAATGTTGTCAGGCAATTAAAGCCGGGTTTTCGACCCTCTTTCCCGATGCGAACTATATCTGTTTGCCGATAGCGGATGGCGGCGAAGGGACGGTGGATGCGATGGTCGCCGCGACGGGCGGCAACATCGTGACGC‑TTGAAGTCTGCGGGCCGATGGGCGAAAAAGGGAATGCTTTTTATGGCCTTACCGGCGACGGGAAAAcggcgggggtttgggatgcgg                                                                                                                                                                               >  1:427303/1‑223 (MQ=255)                                                       cTTAAGTGCAGAAAAATGTTGTCAGGCAATTAAAGCCGGGTTTTCGACCCTCTTTCCCGATGCGAACTATATCTGTTTGCCGATAGCGGATGGCGGCGAAGGGACGGTGGATGCGATGGTCGCCGCGACGGGCGGCAACATCGTGACGC‑TTGAAGTCTGCGGGCCGATGGGCGAAAAAGGGAATGCTTTTTATGGCCTTACCGGCGACGGGAAAACGGCGGGGGTTGAGATGGGGgca                                                                                                                                                                             >  1:762977/1‑238 (MQ=255)                                                          aaGTGCAGAAAAATGTTGTCAGGCAATTAAAGCCGGGTTTTCGACCCTCTTTCCCGATGCGAACTATATCTGTTTGCCGATAGCGGATGGCGGCGAAGGGACGGTGGATGCGATGGTCGCCGCGACGGGCGGCAACGTCGTGAAGC‑TTGAAGTCTGCGGGCCGATGGGCGAAAAAGTGAAGGCTTTTTTTGGCCTTACCGGGGCCGGGAAAAcggcggggggtgaggagggggcgaccg                                                                                                                                                                         >  2:62961/1‑218 (MQ=255)                                                           aGTGCAGAAAAATGTTGTCAGGCAATTAAAGCCGGGTTTTCGACCCTCTTTCCCGATGCGAACTATATCTGTTTGCCGATAGCGGATGGCGGCGAAGGGACGGTGGATGCGATGGTCGCCGCGACGGGCGGCAACATCGTGACGC‑TTGAAGTCTGCGGGCCGATGGGCGAAAAAGTGAATGCTTTTTATGGCCTTACCGGCGACGGGAAAACGGCGGTGATTGAGATGGCGGCAGCAAg                                                                                                                                                                        <  1:59075/239‑1 (MQ=255)                                                                aaaaaatTGGTGTACGGCAATTAAAGCCGGGTTTCCGCCCCTCTTTCCCGGGGCGAACTGTGTTGGTTTGCCGAAAGCGGATGGGGGCGAAGGGACGGTGGAGGCGATGGTCCCCGCGACGGGCGGCAACATCGTGACGC‑TTGAAGTCTGCGGGCCGATGGGCGAAAAAGTGAATGCTTTTTATGGCCTTACCGGCGACGGGAAAACGGCGGTGATTGAGATGGCGGCAGCAAGTGGc                                                                                                                                                                    <  2:660507/236‑1 (MQ=255)                                                                          tggcAGGCAATTAAAGCCGGGTTTGCGACCCTCTTTCCCGATGCGAACTATATCGGTTGGCCGATAGCGGATGGCGGCGAAGGGACGGTGGATGCGATGGTCGACGCGACGGGCGGCAACATCGTGACGC‑TTGAAGTCTGCGGGCCGATGGGCGAAAAAGTGAATGCTTTTTATGGCCTTACCGGCGACGGGAAAACGGCGGTGATTGAGATGGCGGCAGCAAGTGGCCTGATGCTGGt                                                                                                                                                         <  1:204167/236‑1 (MQ=255)                                                                          tgccggggaaTTAAAGCCGGGTTTTGGCACCTTCTTCCGGAGGGGAACTATATCTGGTGGCCGATAGCGGATGGCGGGGAAGGGACGGTGGATGGGATGGTCGCCGCGACGGGCGGCAACATCGTGACGT‑TTGAAGTCTGCGGGCCGATGGGCGAAAAAGTGAATGCTTTTTATGGCCTTACCGGCGACGGGAAAACGGCGGTGATTGAGATGGCGGCAGCAAGTGGCCTGATGCTGGt                                                                                                                                                         <  2:725222/231‑1 (MQ=255)                                                                                    ttAAAGCCGGGTTTTCGACCCTCTTTCCCGATGCGAACTATATCTGTTTGCCGATAGCGGATGGCGGCGAAGGGACGGTGGATGCGATGGTCGCCGCGACGGGCGGCAACATCGTGACGC‑TTGAAGTCTGCGGGCCGATGGGCGAAAAAGTGAATGCTTTTTATGGCCTTACCGGCGACGGGAAAACGGCGGTGATTGAGATGGCGGCAGCAAGTGGCCTGATGCTGGTCGCGCCTGaa                                                                                                                                               <  1:62961/239‑1 (MQ=255)                                                                                    ttAAAGCCGGGTTTTCCACCCTTTTTCCCGATGGGAACTATTTCTGTTTGCCGATAGCGGATGGCGGCGAAGGGGCGGTGGATGCGATGGTCGCCGCGACGGGCGGCAACATCGTGACGC‑TTGAAGTCTGCGGGCCGATGGGCGAAAAAGTGAATGCTTTTTATGGCCTTACCGGCGACGGGAAAACGGCGGTGATTGAGATGGCGGCAGCAAGTGGCCTGATGCTGGTCGCGCCTGaa                                                                                                                                               <  2:1073041/239‑1 (MQ=255)                                                                                                           tttCCCGATGCGAACTATATCTGTTTGGCGATAGCGGATGGCGGCGAAGGGACGGTGGATGCGATGGTCGCCGCGACGGGCGGCAACATCGTGACGC‑TTGAAGTCTGCGGGCCGATGGGCGAAAAAGTGAATGCTTTTTATGGCCTTACCGGCGACGGGAAAACGGCGGTGATTGAGATGGCGGCAGCAAGTGGCCTGATGCTGGTCGCGCCTGAAAAGCGTAATCCGTTGCTGGCCTc                                                                                                                        <  1:345221/239‑1 (MQ=255)                                                                                                           tttCCCGATGCGAACTATATCTGGTTGCCGATAGCGGATGGCGGCGAAGGGACGGTGGATGCGATGGTCGCCGCGACGGGCGGCAACATCGTGACGC‑TTGAAGTCTGCGGGCCGATGGGCGAAAAAGTGAATGCTTTTTATGGCCTTACCGGGGACGGGAAAACGGCGGTGATTGAGATGGCGGCAGCAAGTGGCCTGATGCTGGTCGCGCCTGAAAAGCGTAATCCGTTGCTGGCCTc                                                                                                                        <  1:45116/239‑1 (MQ=255)                                                                                                           gttCCCGAAGCGAACTATATTTGTTTGCCGGTAGCGGATGGCGGTGAAGGGACGGGTGATGCGATGGTCGCCGCGACGGGCGGCAACATCGTGACGC‑TTGAAGTCTGCGGGCCGATGGGCGAAAAAGTGAATGCTTTTTATGGCCTTACCGGCGACGGGAAAACGGCGGTGATTGAGATGGCGGCAGCAAGTGGCCTGATGCTGGTCGCGCCTGAAAAGCGTAATCCGTTGCTGGCCTc                                                                                                                        <  1:557117/238‑1 (MQ=255)                                                                                                                 gATGCGAACTATATCTGTTTGCCGATAGCGGATGGCGGCGAAGGGACGGTGGATGCGATGGTCGCCGCGACGGGCGGCAACATCGTGACGC‑TTGAAGTCTGCGGGCCGATGGGCGAAAAAGTGAATGCTTTTTATGGCCTTACCGGCGACGGGAAAACGGCGGTGATTGAGATGGCGGCAGCAAGTGGCCTGATGCTGGTCGCGCCTGAAAAGCGTAATCCGTTGCTGGCCTCCAGt‑tt                                                                                                                 <  1:50066/239‑1 (MQ=255)                                                                                                                     cGAACTATATCTGTTTGCCGATAGCGGATGGCGGCGAAGGGACGGTGGATGCGATGGTCGCCGCGACGGGCGGCAACATCGTGACGC‑TTGAAGTCTGCGGGCCGATGGGCGAAAAAGTGAATGCTTTTTATGGCCTTACCGGCGACGGGAAAACGGCGGTGATTTGGAATGCGGCAGCAAGTGGCCTGATGCTGGTGCCGCCTGAAAAACGGAATCCGTTGTTGGCCCCCCGG‑TTTTGt                                                                                                             >  1:935500/1‑239 (MQ=255)                                                                                                                       aagtataTCTGTTTGCCGATAGCGGATGGCGGCGAAGGGACGGTGGATGCGATGGTCGCCGCGACGGGCGGCAACATCGTGACGC‑TTGAAGTCTGCGGGCCGATGGGCGAAAAAGTGAATGCTTTTTATGGCCTTACCGGCGACGGGAAAAcg                                                                                                                                                                                                  <  1:595159/150‑1 (MQ=255)                                                                                                                       aaCTATATCTTTTTGCCGATGGGGGGGGGGGGCGAGGGGAGGGTGGAGGCGATGGGCGCCGCGGCGGGGGGCAACATCGTGACGC‑TTGAAGTCTGCGGGCCGATGGGCGAAAAAGTGAATGCTTTTTATGGCCTTACCGGCGACGGGAAAACGGCGGTGATTGAGATGGCGGCAGCAAGTGGCCTGATGCTGGTCGCGCCTGAAAAGCGTAATCCGTTGCTGGCCTCCAGT‑TTTGGTa                                                                                                            <  2:19083/238‑1 (MQ=255)                                                                                                                         cTATATCTGTTTGCCGATAGCGGATGGCGGGGAAGGGACGGTGGATGCGATGGTCGCCGCGACGGGCGGCAACATCGTGACGC‑TTGAAGTCTGCGGGCCGATGGGCGAAAAAGTGAATGCTTTTTATGGCCTTACCGGCGACGGGAAAACGGCGGTGATTGAGATGGCGGCAGCAAGTGGCCTGATGCTGGTCGCGCCTGAAAAGCGTAATCCGTTGCTGGCCTCCAGT‑TTTGGTACgg                                                                                                         <  1:715793/239‑1 (MQ=255)                                                                                                                          tttttcgtgttgcggagggCGGATGGCGGTGAAGGGACGGTGGATGCGTTGGTCGCCGCGACGGGCGGCAAAATCGTGACGA‑TTGAAGTTGGCGGGCCGATGGGCGAAAAAGTGAATGCTTTTTATGGCCTTACCGGCGACGGGAAACCGGCGGTGATTGAGATGGCGGCAGCAAGTGGCCTTATGCTGGTCGCGCTTGAAAAGCTTAATCCGTTGCTGGCCTCCAGT‑TTTGGTGCggg                                                                                                        <  2:707832/221‑1 (MQ=255)                                                                                                                                  tttGCCGATAGCGGATGGCGGCGAAGGGACGGTGGATGCGATGGTCGCCGCGACGGGCGGCAACATCGTGACGC‑TTGAAGTCTGCGGGCCGATGGGCGAAAAAGGGAATGCTTTTTATGGCCTTACCGGCGACGGGAAAACGGCGGTGATTGAGTTGGCGGCAGCAAGTGGCCCGATGCTGGGCGCGCCTGGAAAGCGTAAACCGTTGCTTGCCCCCAGC‑TTTGGGACGGGGGGGTTaa                                                                                                >  1:1054036/1‑239 (MQ=255)                                                                                                                                   ttGCCGATGGGGGGTGGCGGGGAAGGGGGGGGGGATGCGATGGGCGCCGCGACGGGCGGCAAAATCGTGACGC‑TTGAAGTCTGCGGGCCGATGGGCGAAAAAGTGAATGCTTTTTATGGCCTTACCGGCGACGGGAAAACGGCGGTGATTGAGATGGCGGCAGCAAGTGGCCTGATGCTGGTCGCGCCTGAAAAGCGTAATCCGTTGCTGGCCTCCAGT‑TTTGGTACGGGGGAGTTAAt                                                                                               <  2:347816/239‑1 (MQ=255)                                                                                                                                   ttGCCGATGGCGGATGGCGGCGAAGGGACGGTGGATGCGATGGTCGCCGCGACGGGCGGCAACATCGTGACGC‑TTGAAGTCTGCGGGCCGATGGGCGAAAAAGTGAATGCTTTTTATGGCCTTACCGGCGACGGGAAAACGGCGGTGATTGAGATGGCGGCAGCAAGTGGCCTGATGCTGGTCGCGCCTGAAAAGCGTAATCCGTTGCTGGCCTCCAGT‑TTTGGTACGGGGGAGTTAAt                                                                                               <  1:339733/239‑1 (MQ=255)                                                                                                                                    tGCCGATAGCGGATGGCGGCGAAGGGACGGTGGATGCGATGGTCGCCGCGACGGGCGGCAACATCGTGACGC‑TTGAAGTCTGCGGGCCGATGGGCGAAAAAGTGAATGCTTTTTATGGCCTTACCGGCGACGGGAAAACGGCGGTGATTGGGATGGCGGCAGCAAGGGGCCCGGTGCTGGTCGCGCCTGCAAAGCGGAAACCGTTGCTTGCCTCCAGG‑TTTTGTACGGGGGGGTTAAtt                                                                                              >  1:179864/1‑239 (MQ=255)                                                                                                                                    tGCCGATAGCGGATGGCGGCGAAGGGACGGTGGATGCGATGGTCGCCGCGACGGGCGGCAACATCGTGACGC‑TTGAAGTCTGCGGGCCGATGGGCGAAAAAGTGAATGCTTTTTATGGCCTTACCGGCGACGGGAAAACGGCGGTGATTGAGATGGGGGCAGCAAGTGGCCCGAGGCTGGGCGCGCCTGAAAAGCGGAATCCGTTTCTGGGCTCCAGT‑TTTGGTACAGGGGGGGTAAtt                                                                                              >  1:935658/1‑239 (MQ=255)                                                                                                                                    tGCCGATAGCGGATGGCGGCGAAGGGACGGTGGATGCGATGGTCGCCGCGACGGGCGGCAACATCGTGACGC‑TTGAAGTCTGCGGGCCGATGGGCGAAAAAGTGAATGCTTTTTATGGCCTTACCGGCGACGGGAAAACGGCGGGGATTGAGATGGGGGCAGCAAGTGGGCTGATGCTGGTGCGGCCTGAAAAAGGGAAACCGGTTCTGGCCCCCCGG‑TTTTGTCAGGGGGGGGTAAtt                                                                                              >  1:257007/1‑239 (MQ=255)                                                                                                                                                agggcggaGAAGGGACGGTGGATGCGATGGTCGCCGCGACGGGCGGCAACATCGTGACGC‑TTGAAGTCTGCGGGCCGATGGGCGAAAAAGTGAATGCTTTTTATGGCCTTACCGGCGACGGGAAAACGGCGGTGATTGAGATGGCGGCAGCAAGTGGCCTGATGCTGGTCGCGCCTGAAAAGCGTAATCCGTTGCTGGCCTCCAGT‑TTTGGTACGGGGGAGTTAATTCGTCATGCGCt                                                                                   <  1:182150/236‑1 (MQ=255)                                                                                                                                                 tGGCGGCGAAGGGACGGTGGATGCGATGGTCGCCGCGACGGGCGGCAACATCGTGACGC‑TTGAAGTCTGCGGGCCGATGGGCGAAAAAGGGAATGCTTTTTATGGCCTTACCGGCGACGGGAAAACGGCGGGGGTTGGGGGGGCGGCAGAAAGGGGCCGGGGGGTGGGGGCGCCCGAAAAAGGCAAACCGGTGGTGGGCCCCCGTGTTGGGCACgggggggggatctggactaggcggg                                                                                  >  2:145868/1‑219 (MQ=255)                                                                                                                                                             ggcGGTGGATGCGATGGTCGCCGCGACGGGCGGCAACATCGTGACGC‑TTGAAGTCTGCGGGCCGATGGGCGAAAAAGTGAATGCTTTTTATGGCCTTACCGGCGACGGGAAAACGGCGGTGATTGAGATGGCGGCAGCAAGTGGCCTGATGCTGGTCGCGCCTGAAAAGCGTAATCCGTTGCTGGCCTCCAGT‑TTTGGTACGGGGGAGTTAATTCGTCATGCGCTGGATAACGACAtt                                                                      <  1:721002/236‑1 (MQ=255)                                                                                                                                                             gaggggggatgggttgggtcccgcGACGGGCGGCAACATCGTGACGC‑TTGAAGTCTGCGGGCCGATGGGCGAAAAAGTGAATGCTTTTTATGGCCTTACCGGCGACGGGAAAACGGCGGTGATTGAGATGGCGGCAGCAAGTGGCCTGATGCTGGTCGCGCCTGAAAAGCGTAATCCGTTGCTGGCCTCCAGT‑TTTGGTACGGGGGAGTTAATTCGTCATGCGCTGGATAACGACAtt                                                                      <  1:830268/218‑1 (MQ=255)                                                                                                                                                               tgggtgagggggtggtggCCGGGACGGGGGGCAAGA‑CGTGGAGCGTTGAATTTGGCGGGGCGATGGACGAAAAAGGGAATGCTTTTTATGGCCGTACCGGCGACGGGAAAACGGCGGTGATTGAGAGGGCGGCAGCAAGTGGCCGGATGCTGGTCGCGCCTGAAAAGCGTAATCCGTGGCTGGCCTCCAGT‑TTTGGTACGGGGGAGTTAATTCGTCATGCGCGGGATAACGACATTCg                                                                    <  2:46388/226‑1 (MQ=255)                                                                                                                                                               cgtggggtgggggggTCGCCGCGACGGGCGGCAACATCTGGACGC‑TTGAAGTCTGCGGGCCGATGGGCGAAAAAGGGAATGCTTTTTATGGCCTTACCGGCGACGGGAAAACGGCGGTGATTGAGATGGCGGCAGCAAGTGGCCTGATGCTGGTCGCGCCTGAAAAGCGTAATCCGTTGCTGGCCTCCAGT‑TTTGGTACGGGGGAGTTAATTCGTCATGCGCTGGATAACGACATTCg                                                                    <  2:1054036/225‑1 (MQ=255)                                                                                                                                                               cGGTGGATGGGATGGTCGCCGCGACGGGCGGCAACATCGTGACGC‑TTGAAGTCTGCGGGCCGATGGGCGAAAAAGTGAATGCTTTTTATGGCCTTACCGGCGACGGGAAAACGGCGGTGATTGAGATGGCGGCAGCAAGTGGCCTGATGCTGGTCGCGCCTGAAAAGCGTAATCCGTTGCTGGCCTCCAGT‑TTTGGTACGGGGGAGTTAATTCGTCATGCGCTGGATAACGACATTCg                                                                    <  1:145868/238‑1 (MQ=255)                                                                                                                                                               cGGTGGATGCGATGGTCGCCGCGACGGGCGGCAACATCGTGACGC‑TTGAAGTCTGCGGGCCGATGGGCGAAAAAGTGAATGCTTTTTATGGCCTTACCGGCGACGGGAAAACGGCGGTGATTGAGATGGCGGCAGCAAGTGGCCTGATGCTGGTCGCGCCTGAAAAGCGTAATCCGTTGCTGGCCTCCAGT‑TTTGGTACGGGGGAGTTAATTCGTCATGCGCTGGATAACGACATTCg                                                                    <  1:114855/238‑1 (MQ=255)                                                                                                                                                                ggggggtgggggggTCGCCGGGGCGGGCGGCAACATCGTGACGG‑TTGAAGGTTGCGGGCCGATGGGCGAAAAAGTGAATGCTTTTTATGGCCTTACCGGCGACGGGAAAACGGCGGTGATTGAGATGGCGGCAGCAAGTGGCCTGATGCTGGTCGCGCCTGAAAAGCGTAATCCGTTGCTGGCCTCCAGT‑TTTGGTACGGGGGAGTTAATTCGTCATGCGCTGGATAACGACATTCGcc                                                                  <  2:935658/227‑1 (MQ=255)                                                                                                                                                                    gATGCGATGGTCGCCGCGACGGGCGGCAACATCGTGACGC‑TTGAAGTCTGCGGGCCGATGGGCGAAAAAGTGAATGCTTTTTATGGCCTTACCGGCGACGGGAAAACGGCGGTGATTGGGATGGCGGCAGCAAGTGGCCTGGTGCTGGGCGGGCCTGAAAAGGGTAATCCCTTGCTGGGCCCCCGT‑TTTGGGACGGGGGGGTTAAATTGGCATGGGCggggaaaaaacaattgccaaac                                                              >  2:303537/1‑220 (MQ=255)                                                                                                                                                                                       cGGGCGGAAACAGCGGGACGC‑TTGAAGTCTGCGGGACGAGGGGCGAAAAAGTGAATGCTTTTTATGGCCTTACCGGCGACGGGAAAACGGCGGTGATTGAGATGGCGGCAGCAAGTGGCCTGATGCTGGTCGCGCCTGAAAAGCGTAATCCGTTGCTGGCCTCCAGT‑TTTGGTACGGGGGAGTTAATTCGTCATGCGCTGGATAACGACATTCGCCATATTATTCTCGGCATTggcgg                                            <  1:536959/238‑1 (MQ=255)                                                                                                                                                                                         gggggcAACATCGGGACGC‑TTGAAGTCGGGGGGCCGATGGGCGAAAAAGGGAATGTTTTTTATGGCCTTAGCGGGGACGGGAAAACGGGGGTGATTGAGATGGCGGAAGCAAGTGGCCTGATGCTGGTCGGGCCTGAAAAGCGTAATCCGTTGCTGGCCTCCAGT‑TTTGGTACGGGGGAGTTAATTCGTCATGCGCTGGATAACGACATTCGCCATATTATTCTCGGCATTGGCGGCa                                          <  2:762977/235‑1 (MQ=255)                                                                                                                                                                                             gcAACATCTTGGCGC‑TTGAAGTTGGGGGGGCGATGGGCGAAAAAGTGAATGCTTTTTATGGCCTTACCGGGGACGGGAAAACGGCGGTGATTGAGATGGCGGCAGCAAGTGGCCTGATGCTGGTCGCGCCTGAAAAGCGTAATCCGTTGCGGGCCTCCAGT‑TTTGGTACGGGGGAGTTAATTCGTCATGCGCTGGATAACGACATTCGCCATATTATTCTCGGCATTGGCGGCAGTGCg                                     <  2:194984/239‑1 (MQ=255)                                                                                                                                                                                                   atttgggcggttGAGGTCGGGGGGCCGGTGGGCGAAAAAGTGAATGCTTTTTATGGCCTTACCGGGGAGGGGAAAACGGCGGTGATTGGGATGGCGGCAGCAAGTGGCCTGATGCTGGTCGCGCCTGAAAAGCGTAATCCGTTGCTGGCCTCCAGT‑TTTGGTACGGGGGAGTTAATTCGTCATGCGCTGGATAACGACATTCGCCATATTATTCTCGGCATTGGCGGCAGTGCGACGGt                                <  2:427303/229‑1 (MQ=255)                                                                                                                                                                                                       gccGC‑TTGAAGTCTGCGGGCCGATGGGAGAAAAAGTGAATGCTTTTTATGGCCTTACCGGGGACGGGAAAACGGCGGTGATTGAGATGGCGGCAGCAAGTGGCCTGATGCTGGTCGCGCCTGAAAAGCGTAATCCGTTGCTGGCCTCCAGT‑TTTGGTACGGGGGAGTTAATTCGTCATGCGCTGGATAACGACATTCGCCATATTATTCTCGGCATTGGCGGCAGTGCGACGGTCGAcg                           <  1:55087/237‑1 (MQ=255)                                                                                                                                                                                                       gACGC‑TTGAAGTCTGGGGGCCGTTGGGCGAAAAAGTGAATGCTTTTTATGGCCTTACCGGCGACGGGAAAACGGCGGTGATTGAGATGGCGGCAGCAAGTGGCCTGATGCTGGTCGCGCCTGAAAAGCGTAATCCGTTGCTGGCCTCCAGT‑TTTGGTACGGGGGAGTTAATTCGTCATGCGCTGGATAACGACATTCGCCATATTATTCTCGGCATTGGCGGCAGTGCGACGGTCGAcg                           <  1:343613/239‑1 (MQ=255)                                                                                                                                                                                                        gccgtttGAATTTTGGGGGCCGGTGGGGGAAAAAGTGAATGTTTTTTTTGGCCTTACCGGGGAAGGGAAAACGGCGGTGATTGAGATGGCGGCAGCAAGTGGCCTGATGCTGGTCGCGCCTGAAAAGCGTAATCCGTTGCTGGCCTCCAGT‑TTTGGTACGGGGGAGTTAATTCGTCATGCGCTGGATAACGACATTCGCCATATTATTCTCGGCATTGGCGGCAGTGCGACGGTCGAcg                           <  2:935500/234‑1 (MQ=255)                                                                                                                                                                                                          gC‑TTGAAGTCTGCGGGCCGATGGGCGAAAAAGTGAATGCTTTTTATGGCCTTACCGGCGACGGGAAAACGGCGGTGATTGAGATGGCGGCAGCAAGTGGCCTGATGCTGGTCGCGCCTGAAAAGCGTAATCCGTTGCTGGCCTCCAGT‑TTTGGTACGGGGGAGTTAATTCGTCATGCGCTGGATAACGACATTCg                                                                    <  1:49941/195‑1 (MQ=255)                                                                                                                                                                                                          cggttGAAGTCTGCGGGCCGATGGGCGAAAAAGTGAATGCTTTTTATGGCCTTACCGGCGACGGGAAAACGGCGGTGATTGAGATGGCGGCAGCAAGTGGCCTGATGCTGGTCGCGCCTGAAAAGCGTAATCCGTTGCTGGCCTCCAGT‑TTTGGTACGGGGGAGTTAATTCGTCATGCGCTGGATAACGACATTCGCCATATTATTCTCGGCATTGGCGGCAGTGCGACGGTCGAcgg                          <  1:87537/235‑1 (MQ=255)                                                                                                                                                                                                               gAAGTCTGCGGGCCGATGGGCGAAAAAGTGAATGCTTTTTATGGCCTTACCGGCGACGGGCAAACGGCGGTGGTTGAGATGGGGGCAGGAAGTGGCCTGGTGCTGGTGGGCCCTGGAAAGGGGAATTCGGTGCTGGCCTCCAGG‑TTTTGGACGGGGGGGGTAATTCGTCATGCGCTGGGTAACGACAATCGCCACATTGTTTTGGGCATTTGCGGGAGGGCgcggggggaggggggagg                    >  2:691840/1‑222 (MQ=255)                                                                                                                                                                                                                     tgggggCCGGTGGGCGAAAAAGGGAAGGCTTTTTATGGCCTTACCGGCGAAGGGAAAACGGCGGTGATTGAGATGGCGGCAGCAAGTGGCCTGATGCTGGTCGCGCCTGAAAAGCGTAATCCGTTGCTGGCCTCCAGT‑TTTGGTACGGGGGAGTTAATTCGTCATGCGCTGGATAACGACATTCGCCATATTATTCTCGGCATTGGCGGCAGTGCGACGGTCGACGGCGGTATGGGCa               <  1:1060246/235‑1 (MQ=255)                                                                                                                                                                                                                            ccttGGGCGAAAATGTGAATTTTTTTTAGGGCCTTTCCGGGGAGGGGAAAACGGGGGGTATTGGGAGGGGGGCAGCAAGTGGCCTGGTGCTGGTCCAGACTGAAAAGCGTAATCCGTTGCTGGCCTCCAGT‑TTTGGTACGGGGGAGTTAATTCGTCATGCGCTGGATAACGACATTCGCCATATTATTCTCGGCATTGGCGGCAGTGCGACGGTCGACGGCGGTATGGGCATGGCGCa        <  2:531359/235‑1 (MQ=255)                                                                                                                                                                                                                               tGGGCGAAAAAGTGAATGCTTTTTATGGCCTTACCGGCGACGGGAAAACGGCGGTGATTGAGATGGCGGCAGCAAGTGGCCTGATGCTGGTCGCGCCTGAAAAGCGTAATCCGTTGCTGGCCTCCAGT‑TTTGGTACGGGGGAGTTAATTCGTCATGCGCTGGATAACGACCTTCGCCCTATTATTCTCGGCATTGGCGGCAGGTGGGCGGGCGACGGCGGGAAGGGCCTGgcgccagg     >  1:972203/1‑234 (MQ=255)                                                                                                                                                                                                                                  gCGAAAAAGTGAATGCTTTTTATGGCCTTACCGGCGACGGGAAAACGGCGGTGATTGAGATGGCGGCAGCAAGTGGCCTGATGCTGGTCGCGCCTGAAAAGCGTAATCCGTTGCTGGCCTCCAGT‑TTTGGTACGGGGGAGTTAATTCGTCATGCGCTGGATAACGACATTCGCCATACTATTCTCGGCATTGGCGGCAGTGCGACGGTCGACGGCGGTATGGGCATGGCGCAGGCGCt  <  1:303537/238‑1 (MQ=255)                                                                                                                                                                                                                                           |                                                                                                                                                                                                                                      TGAAAAAGGATTGATAATGAAGATTGTCATTGCGCCAGACTCTTTTAAAGAGAGCTTAAGTGCAGAAAAATGTTGTCAGGCAATTAAAGCCGGGTTTTCGACCCTCTTTCCCGATGCGAACTATATCTGTTTGCCGATAGCGGATGGCGGCGAAGGGACGGTGGATGCGATGGTCGCCGCGACGGGCGGCAACATCGTGACGC‑TTGAAGTCTGCGGGCCGATGGGCGAAAAAGTGAATGCTTTTTATGGCCTTACCGGCGACGGGAAAACGGCGGTGATTGAGATGGCGGCAGCAAGTGGCCTGATGCTGGTCGCGCCTGAAAAGCGTAATCCGTTGCTGGCCTCCAGT‑TTTGGTACGGGGGAGTTAATTCGTCATGCGCTGGATAACGACATTCGCCATATTATTCTCGGCATTGGCGGCAGTGCGACGGTCGACGGCGGTATGGGCATGGCGCAGGCGCT  >  NC_000913/541872‑542333

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