Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	NC_000913	3,873,698	G→A	Y259Y (TAC→TAT)	dgoK ←	2‑dehydro‑3‑deoxygalactonokinase

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	3,873,698	0	G	A	100.0%	171.7 / NA	51	Y259Y (TAC→TAT)	dgoK	2‑dehydro‑3‑deoxygalactonokinase
Reads supporting (aligned to +/- strand):  ref base G (0/0);  new base A (29/22);  total (29/22)

CAGGTTTCAGTACCGTACCTGCGCCAATCAACGCCTTGTCGCCGTACGCATCAACGATGGCGGGAATGCTTTGCTCCCATTGTGGGGAATTCAGCGGGATTTCAACCGCGTCGAACCCGGCGTCAATCACCGCGCCAACATGCGCCAGCGCCTCGTCGGGCGTAATACCGCGCAAAATGGCGATCAGCGGGAGTTTAGTTTGCCACTGCATGAGCGATGCTCCTTATACCAGCCTGAAATGCCGTGTCGCCCGCCACCGCCGTCACGTCGCAACCCATCGCCTGAAAGGCTTGCTGGTAGCGCGCGGTCAGCGATGTTCCGGCGACAAGGGTGATGGCGTGTTGATGGGCCACATAGTCGCGCATACTGGCGACCTCTGCGCCAATCAACAAACCAGAGAGAAATTCGCTGACCTGTTCGCGGGGAAGTGTTCCCAGCACATGCGAGGCGCGAACTTCAAAAAGCTGCGGCAATATGGCGGGCGTATTAAGACCACGCTCAAGGCCAGCTGTGAAGGCATCGGCAGAGTTTTCCTGCGGCGGCAAACCTGCGCCAATCAATGAGTGATTTAACAGTAAATGATGTAATTCA  >  NC_000913/3873402‑3873992                                                                                                                                                                                                                                                                                                         |                                                                                                                                                                                                                                                                                                       cAGGTTTCAGTACCGTACCTGCGCCAATCAACGCCTTGTCGCCGTACGCATCAACGATGGCGGGAATGCTTTGCTCCCATTGTGGGGAATTCAGCGGGATTTCAACCGCGTCGAACCCGGCGTCAATCACCGCGCCAACATGCGCCAGCGCCTCGTCGGGCGTAATACCGCGCAAAATGGCGATCAGCGGGAGTTTAGTTTGCCACTGCATGAACGATGCTCCTTATACCAGCCTGAAATGCCGTGTCGCCCGCCACCGCCGTCACGTCGCAACCCATCGCCTGAAAGGCTTGCTGATAg                                                                                                                                                                                                                                                                                                     <  1:104533/300‑1 (MQ=255)     tttCAGTACCGTACCTGCGCCAATCAACGCCTTGTCGCCGTACGCATCAACGATGGCGGGAATGCTTTGCTCCCATTGTGGGGAATTCAGCGGGATTTCAACCGCGTCGAACCCGGCGTCAATCACCGCGCCAACATGCGCCAGCGCCTCGTCGGGCGTAATACCGCGCAAAATGGCGATCAGCGGGAGTTTAGTTTGCCACTGCATGAACGATGCTCCTTATACCAGCCTGAAATGCCGTGTCGCCCGCCACCGCCGTCACGTCGCAACCCATCGCCTGAAAGGCTTGCTGATAgcgcg                                                                                                                                                                                                                                                                                                 <  2:279952/300‑1 (MQ=255)                              aaCGCCTTGTCATCGTACGCATCGACGATCGCGGGAATGCTTTGCTCCCATTGTGGGGAATTCAGCTGGATGTCAACCGCGTCGAACCCGTCGTTAATCGCCGCGCCAACATGCGCCAGCGCCTCGTCGTTCGTAATACCGCGCAAAACTACAATCAGCGGGAGTTTAGTTTGCCACTGCATGAACGATGCTCTTTATACCAGACTGAGATGCCGTGTCGCCCGCCACCGCCGTTACGTCGCAACCCATCTCCTGAAAGACTTGCTGATAACGCGCGGTCAGCGATGTTCCGTCGAcggg                                                                                                                                                                                                                                                                        <  2:234684/300‑4 (MQ=255)                                  ccTTGTCGCCGTACGCATCAACGATGGCGGGAATGCTTTGCTCCCATTGTGGGGAATTCAGCGGGATTTCAACCGCGTCGAACCCGGCGTCAATCACCGCGCCAACATGCGCCAGCGCCTCGTCGGGCGTAATACCGCGCAAAATGGCGATCAGCGGGAGTTTAGTTTGCCACTGCATGAACGATGCTCCTTATACCAGCCTGAAATGCCGTGTCGCCCGCCACCGCCGTCACGTCGCAACCCATCGCCTGAAAGGCTTGCTGATAGCGCGCGGTCAGCGATGTTCCGGCGAAAAGGGTg                                                                                                                                                                                                                                                                    >  2:543018/1‑300 (MQ=255)                                                gCATCAACGATGGCGGGAATGCTTTGCTCCCATTGTGGGGAATTCAGCGGGATTTCAACCGCGTCGAACCCGGCGTCAATCACCGCGCCAACATGCGCCAGCGCCTCGTCGGGCGTAATACCGCGCAAAATGGCGATCAGCGGGAGTTTAGTTTGCCACTGCATGAACGATGCTCCTTATACCAGCCTGAAATGCCGTGTCGCCCGCCACCGCCGTCACGTCGCAACCCATCGCCTGAAAGGCTTGCTGATAGCGCGCGGTCAGCGATGTTCCGGCGACAAGGGTGATGGCGTGTTGATg                                                                                                                                                                                                                                                      <  1:438428/300‑1 (MQ=255)                                                 cATCAACGATGGCGGGAATGCTTTGCTCCCATTGTGGGGAATTCAGCGGGATTTCAACCGCGTCGAACCCGGCGTCAATCACCGCGCCAACATGCGCCAGCGCCTCGTCGGGCGTAATACCGCGCAAAATGGCGATCAGCGGGAGTTTAGTTTGCCACTGCATGAACGATGCTCCTTATACCAGCCTGAAATGCCGTGTCGCCTGCCACCGCCGTCACGTCGCAACCCATCGCCTGAAAGGCTTGCTGATAGCGCGCGGTCAGCGATGTTCCGGCGACAAGGGTGATGGCGTGTTGATgg                                                                                                                                                                                                                                                     >  1:842242/1‑300 (MQ=255)                                                          tGGCGGGAATGCTTTGCTCCCATTGTGGGGAATTCAGCGGGATTTCAACCGCGTCGAACCCGGCGTCAATCACCGCGCCAACATGCGCCAGCGCCTCGTCGGGCGTAATACCGCGCAAAATGGCGATCAGCGGGAGTTTAGTTTGCCACTGCATGAACGATGCTCCTTATACCAGCCTGAAATGCCGTGTCGCCCGCCACCGCCGTCACGTCGCAACCCATCGCCTGAAAGGCTTGCTGATAGCGCGCGGTCAGCGATGTTCCGGCGACAAGGGTGATGGCGTGTTGATGGGCCACATAg                                                                                                                                                                                                                                            <  2:538349/300‑1 (MQ=255)                                                                     cTTTGCTCCCATTGTGGGGAATTCAGCGGGATTTCAACCGCGTCGAACCCGGCGTCAATCACCGCGCCAACATGCGCCAGCGCCTCGTCGGGCGTAATACCGCGCAAAATGGCGATCAGCGGGAGTTTAGTTTGCCACTGCATGAACGATGCTCCTTATACCAGCCTGAAATGCCGTGTCGCCCGCCACCGCCGTCACGTCGCAACCCATCGCCTGAAAGGCTTGCTGATAGCGCGCGGTCAGCGATGTTCCGGCGACAAGGGTGATGGCGTGTTGATGGGCCACATAGTCGCGCATACt                                                                                                                                                                                                                                 <  1:90361/300‑1 (MQ=255)                                                                                   tgGGGAATTCATCGGGATTTCAACCGCTTCGAACCCGGCGTCAATCACCGCGCCAACATGCGCCAGCGCCTCGTCGGGCGTAATACCGCGCAAAATGGCGATCAGCGGGAGTTTAGTTTGCCACTGCATGAACGATGCTCCTTATACCAGCCTGAAATGCCGTGTCGCCCGCCACCGCCGTCACGTCGCAACCCATCGCCTGAAATGCTTGCTGATAGCGCGCGGTCAGCGATGTTCCGGCGACAAGGGTGATGGCGTGTTGATGGGCCACATAGTCGCGCATACTGGCGACCTCTgcgc                                                                                                                                                                                                                   <  2:337052/300‑1 (MQ=255)                                                                                   tgGGGAATTCAGCGGGATTTCAACCGCGTCGAACCCGGCGTCAATCACCGCGTCAACATGCGCCAGCGCCTCGTCGGGCGTAATACCGCGCAAAATTGCGATCAGCGGTAGTTTAGTTTGCCACTGCATGAACGATGCTCCTTATACCAGCCTGAAATGCCGTGTCGCCCGCCACCGCCGTCACGTCGCAACCCATCGCTTGAAAGGCTTGCTGATAGCGCGCGGTCAGCGAGGTTCCGGCGACAAGGGTGATGGCGTGTTGATGGGCCACATAGTCGCGCATACTGGCGACCGCTgcgc                                                                                                                                                                                                                   <  1:786765/300‑1 (MQ=255)                                                                                   tgGGGAATTCAGCGGGATTTCAACCGCGTCGAACCCGGCGTCAATCACCGCGCCAACATGCGCCAGCGCCTCGTCGGGCGTAATACCGCGCAAAATGGCGATCAGCGGGAGTTTAGTTTGCCACTGCATGAACGATGCTCCTTATACCAGCCTGAAATGCCGTGTCGCCCGCCACCGCCGTCACGTCGCAACCCATCGCCTGAAAGGCTTGCTGATAGCGCGCGGTCAGCGATGTTCCGGCGACAAGGGTGATGGCGTGTTGATGGGCCACATAGTCGCGCATACTGGCGACCTCTgcgc                                                                                                                                                                                                                   <  1:50587/300‑1 (MQ=255)                                                                                   tgGGGAATTCAGCGGGATTTCAACCGCGTCGAACCCGGCGTCAATCACCGCGCCAACATGCGCCAGCGCCTCGTCGGGCGTAATACCGCGCAAAATGGCGATCAGCGGGAGTTTAGTTTGCCACTGCATGAACGATGCTCCTTATACCAGCCTGAAATGCCGTGTCGCCCGCCACCGCCGTCACGTCGCAACCCATCGCCTGAAAGGCTTGCTGATAGCGCGCGGTCAGCGATGTTCCGGCGACAAGGGTGATGGCGTGTTGATGGGCCACATAGTCGCGCATACTGGCGACCTCTgcgc                                                                                                                                                                                                                   <  2:185203/300‑1 (MQ=255)                                                                                       gAATTCAGCGGGATTTCAACCGCGTCGAACCCGGCGTCAATCACCGCGCCAACATGCGCCAGCGCCTCGTCGGGCGTAATACCGCGCAAAATGGCGATCAGCGGGAGTTTAGTTTGCCACTGCATGAACGATGCTCCTTATACCAGCCTGAAATGCCGTGTCGCCCGCCACCGCCGTCACGTCGCAACCCATCGCCTGAAAGGCTTGCTGATAGCGCGCGGTCAGCGATGTTCCGGCGACAAGGGTGATGGCGTGTTGATGGGCCACATAGTCGCGCATACTGGCGACCTCTGCGCCAAt                                                                                                                                                                                                               >  1:452560/1‑300 (MQ=255)                                                                                                            gcgTCGAACCCGGCGTCAATCACCGCGCCAACATGCGCCAGCGCCTCGTCGGGCGTAATACCGCGCAAAATGGCGATCAGCGGGAGTTTAGTTTGCCACTGCATGAACGATGCTCCTTATACCAGCCTGAAATGCCGTGTCGCCCGCCACCGCCGTCACGTCGCAACCCATCGACTGAAAGGCTTGCTGATAGCGCGCGGTCAGCGATGTTCCGGCGACAAGGGTGATGGCGTGTTGATGGGCCACATAGTCGCGCATACTGGCGACCTCTGCGCCAATCAACAAACCAGAGAGAAATTc                                                                                                                                                                                          >  1:637029/1‑300 (MQ=255)                                                                                                                              aTCACCGCGCCAACATGCGCCAGCGCCTCGTCGGGCGTAATACCGCGCAAAATGGCGATCAGCGGGAGTTTAGTTTGCCACTGCATGAACGATGCTCCTTATACCAGCCTGAAATGCCGTGTCGCCCGCCACCGCCGTCACGTCGCAACCCATCGCCTGAAAGGCTTGCTGATAGCGCGCGGTCAGCGATGTTCCGGCGACAAGGGTGATGGCGTGTTGATGGGCAACATAGTCGCGCATACTGGCGACCTCTGCGCCAATCAACAAACCAGAGAGAAATTCGCTGACCTGTTCGCgggg                                                                                                                                                                        >  1:426598/1‑300 (MQ=255)                                                                                                                                 aCCGCGCCAACATGCGCCAGCGCCTCGTCGGGCGTAATACCGCGCAAAATGGCGATCAGCGGGAGTTTAGTTTGCCACTGCATGAACGATGCTCCTTATACCAGCCTGAAATGCCGTGTCGCCCGCCACCGCCGTCACGTCGCAACCCATCGCCTGAAAGGCTTGCTGATAGCGCGCGGTCAGCGATGTTCCGGCGACAAGGGTGATGGCGTGTTGATGGGCCACATAGTCGCGCATACTGGCGACCTCTGCGCCAATCAACAAACCAGAGAGAAATTCGCTGACCTGTTCGCGGGGAAg                                                                                                                                                                     <  2:827226/300‑1 (MQ=255)                                                                                                                                    gcgcCAACATGCGCCAGCGCCTCGTCGGGCGTAATACCGCGCAAAATGGCGATCAGCGGGAGTTTAGTTTGCCACTGCATGAACGATGCTCCTTATACCAGCCTGAAATGCCGTGTCGCCCGCCACCGCCGTCACGTCGCAACCCATCGCCTGAAAGGCTTGCTGATAGCGCGCGGTCAGCGATGTTCCGGCGACAAGGGTGATGGCGTGTTGATGGGCCACATAGTCGCGCATACTGGCGACCTCTGCGCCAATCAACAAACCAGAGAGAAATTCGCTGACCTGTTCGCGGGGAAgtgt                                                                                                                                                                  >  1:521136/1‑300 (MQ=255)                                                                                                                                             tGCGCCAGCGCCTCTTCGGGCGTAATACCGCGCAAAATGGCGATCAGCGGGAGTTTAGTTTGCCACTGCATGAACGATGCTCCTTATACCAGCCTGAAATGCCGTGTCGCCCGCCACCGCCGTCACGTCGCAACCCATCGCCTGAAAGGCTTGCTGATAGCGCGCGGTCAGCGATGTTCCGGCGACAAGGGTGATGGCGTGTTGATGGGCCACATAGTCGCGCATACTGGCGACCTCTGCGCCAATCAACAAACCAGAGAGAAATTCGCTGACCTGTTCGCGGGGAAGTGTTCCCAGcac                                                                                                                                                         <  1:101210/300‑1 (MQ=255)                                                                                                                                                   aGCGCCTCGTCGGGCGTAATACCGCGCAAAATGGCGATCAGCGGGAGTTTAGTTTGCCACTGCATGAACGATGCTCCTTATACCAGCCTGAAATGCCGTGTCGCCCGCCACCGCCGTCACGTCGCAACCCATCGCCTGAAAGTCTTGCCGATAGCGCGCGGTCATCGATGTTCCGGCGACAAGGTCGATGGCATGTTGATGGGCCACATAGTCGCGCATACTGGCGACCTCTGCGCCAATCAACAAACCAGAGAGAAATTCGCCGACCTGTTCGCGTGGAAGTGTTCCCAGCACATGCGa                                                                                                                                                   <  2:411589/300‑1 (MQ=255)                                                                                                                                                                ggggaaTACCGCGCAAAATGGCGATCAGCGGGAGTTTAGTTTGCCACTGCATGAACGATGCTCCTTATACCAGCCTGAAATGCCGTGTCGCCCGCCGCCGCCGTCACGTAGCAACCCATCGCCTGAAAGGCTTGCTGATAGCGCGCGGTCAGCGGTGTTCCGGCGAAAAGGGTGATGGCGTGTTGATGGGCCACATAGTCGCGCATACTGGCGACCTCTGCGCCAATCAACAAACCCGAGAGAAATTCGCTGACCTGTTCGCGGGGAAGTATTCCCAGCACATGCGAGGTGCGAACTTCa                                                                                                                                      >  2:521363/5‑300 (MQ=255)                                                                                                                                                                gCGTAATCCCGCGCAAAATGGCGATCAGCGGGAGTTTAGTTTGCCACTGCATGAACGATGCTCCTTATACCAGCCTGAAATGCCGTGTCGCCCGCCACCGCCGTCACGTCGCCACCCATCGCCTGAAAGGCTTGCTGATAGCGCGCGGTCAGCGATGTTCCGGCGACAAGGGTGATGGCGTGTTGATGGGCCACATAGTCGCGCATACTGGCGACCTCTGCGCCAATCAACAAACCAGAGAGAAATTCGCTGACCTGTTCGCGGGGAAGTGTTCCCAGCACATGCGAGGCGCGAACTTcc                                                                                                                                      >  2:354101/1‑299 (MQ=255)                                                                                                                                                                gCGTAATACCGCGCAAAATGGCGATCAGCGGGAGTTTAGTTTGCCACTGCATGAACGATGCTCCTTATACCAGCCTGAAATGCCGTGTCGCCCGCCACCGCCGTCACGTCGCAACCCATCGCCTGAAAGGCTTGCTGATAGCGCGCGGTCAGCGATGTTCCGGCGACAAGGGTGATGGCGTGTTGATGGGCCACATAGTCGCGCATACTGGCGACCTCTGCGCCAATCAACACACCAGAGAGAAATTCGCTGACCTGTTCGCGGGGGAGTGTTCACAGCACATGCGAGGCGCGAACTTCa                                                                                                                                      >  2:382504/1‑300 (MQ=255)                                                                                                                                                                gCGTAATACCGCGCAAAATGGCGATCAGCGGGAGTTTAGTTTGCCACTGCATGAACGATGCTCCTTATACCAGCCTGAAATGCCGTGTCGCCCGCCACCGCCGTCACGTCGCAACCCATCGCCTGAAAGGCTTGCTGATAGCGCGCGGTCAGCGATGTTCCGGCGACAAGGGTGATGGCGTGTTGATGGGCCACATAGTCGCGCATACTGGCGACCTCTGCGCCAATCAACAAACCAGAGAGAAATTCGCTGACCTGTTCGCGGGGAAGTGTTCCCAGCACATGCGAGGCGCGAACTTCa                                                                                                                                      >  1:64512/1‑300 (MQ=255)                                                                                                                                                                gCGTAATACCGCGCAAAATGGCGATCAGCGGGAGTTTAGTTTGCCACTGCATGAACGATGCTCCTTATACCAGCCTGAAATGCCGTGTCGCCCGCCACCGCCGTCACGTCGCAACCCATCGCCTGAAAGGCTTGCTGATAGCGCGCGGTCAGCGATGTTCCGGCGACAAGGGTGATGGCGTGTTGATGGGCCACATAGTCGCGCATACTGGCGACCTCTGCGCCAATCAACAAACCAGAGAGAAATTCGCTGACCTGTTCGCGGGGAAGTGTTCCCAGCACATGCGAGGCGCGAACTTCa                                                                                                                                      >  2:441702/1‑300 (MQ=255)                                                                                                                                                                gCGTAATACCGCGCAAAATGGCGATCAGCGGGAGTTTAGTTTGCCACTGCATGAACGATGCTCCTTATACCAGCCTGAAATGCCGTGTCGCCCGCCACCGCCGTCACGTCGCAACCCATCGCCTGAAAGGCTTGCTGATAGCGCGCGGTCAGCGATGTTCCGGCGACAAGGGTGATGGCGTGTTGATGGGCCACATAGTCGCGCATACTGGCGACATCTGCGCCAATCAACAAACCAGAGAGAAANTCGCTGACCTGTTCGCGGGGAAGTGTTCCCAGCACATGCGAGGAGCGNACTTCa                                                                                                                                      >  2:710627/1‑300 (MQ=255)                                                                                                                                                                gCGTAATACCGCGCAAAATGGCGATCAGCGGGAGTTTAGTTTGCCACTGCATGAACGATGCTCCTTATACCAGCATGAAATGCCGTGTCGCCCGCCACCGCCGTCACGTCGCAACCCATCGCCTGAAAGGCTTGCTGATAGCGCGCGGTCAGCGATGTTCCGGCGACAAGGGTGATGGCGTGTTGATGGGCCACATAGTCGCGCATACTGGCGACCTCTGCGCCAATCAACAAACCAGAGAGAAATTCGCTGACCTGTTCGCGGGGAAGTGTTCCCAGCACATGCGAGGCGCGAACTTCa                                                                                                                                      >  2:521452/1‑300 (MQ=255)                                                                                                                                                                  gTAATACCGCGCAAAATGGCGATCAGCGGGAGTTTAGTTTGCCACTGCATGAACGATGCTCCTTATACCAGCCTGAAATGCCGTGTCGCCCGCCACCGCCGTCACGTCGCAACCCATCGCCTGAAAGGCTTGCTGATAGCGCGCGGTCAGCGATGTTCCGGCGACAAGGGTGATGGCGTGTTGATGGGCCACATAGTCGCGCATACTGGCGACCTCTGCGCCAATCAACAAACCAGAGAGAAATTCGCTGACCTGTTCGCGGGGAAGTGTTCCCAGCACATGCGAGGCGCGAACTTCaaa                                                                                                                                    >  1:844727/1‑300 (MQ=255)                                                                                                                                                                    aaTACCGCGCAAAATTGCGATCAGCGTTAGTTTAGTTTGCCACTGCATGAACGATGCTCCTTATACCAGCCTGAAATGCCGTGTCGCCCGCCACCGCCGTCACGTCGCAACCCATCGCCTGAAAGGCTTGCTGATAGCGCGCGGTCAGCGATGTTCCGGCGACAAGGGTGATGGCGTGTTGATGGGCCACATAGTCGCGCATACTGGCGACCTCTGCGCCAATCAACAAACCAGATAGAAATTCGCTGACCTGTTCGCGGGGAAGTGTTCCCAGCACATGCGAGGCGCGAACTTCaaaaa                                                                                                                                  <  2:163870/300‑1 (MQ=255)                                                                                                                                                                                         cAGCGGGAGTTTAGTTTGCCACTGCATGAACGATGCTCCTTATACCAGCCTGAAATGCCGTGTCGCCCGCCACCGCCGTCACGTCGCAACCCATCGCCTGAAAGGCTTGCTGATAGCGCGCGGTCAGCGATGTTCCGGCGACAAGGGTGATGGCGTGTTGATGGGCCACATAGTCGCGCATACTGGCGACCTCTGCGCCAATCAACAAACCAGAGAGAAATTCGCTGACCTGTTCGCGGGGAAGTGTTCCCAGCACATGCGAGGCGCGAACTTCAAAAAGCTGCGGCAATATggcgggcg                                                                                                             <  2:771911/300‑1 (MQ=255)                                                                                                                                                                                              tgAGTTTAGTTTGCCACTGCATGAACGATGCTCCTTATACCAGCCTGAAATGCCGTGTCGCCCGCCACCGCCGTCACGTCGCAACCCATCGCCTGAAAGGCTTGCTGATAGCGCGCGGTCAGCGATGTTCCTGCGACAAGGGTGATGGCGTGTTGATGGGCCACATAGTCGCGCATACTGGCGACCTCTGCGCCAATCAACAAACCAGAGAGAAATTCGCTGACCTGTTCGCGGGGAAGTGTTCCCAGCACATGCGAGGCGCGAACTTCAAAAAGCTGCGGCAATATGGCGGGCGTATTa                                                                                                        <  1:50252/299‑1 (MQ=255)                                                                                                                                                                                                 gtttagtttGCCACTGCATGAACGATGCTCCTTATACCAGCCTGAAATGCCGTGTCGCCCGCCACCGCCGTCACGTCGCAACCCATCGCATGAAAGGCTTGCTGATAGCGCGCGGTCAGCGATGTTCCGGCGACAAGGGTGATGGCGTGTTGATGGGCCACATAGTCGCGCATACTGGCGACCTCTGCGCCAATCAACAAACCAGAGAGAAATTCGCTGACCTGTTCGCGGGGAAGTGTTCCCAGNACATGCGAGGCGCGAACTTCAAAAAGCTGCGGCAATATGGCGGGCGTNTTaaaa                                                                                                     >  2:706948/1‑298 (MQ=255)                                                                                                                                                                                                      gtttGCCACTGCATGAACGATGCTCCTTATACCAGCCTGAAATGCCGTGTCGCCCGCCACCGCCGTCACGTCGCAACCCATCGCCTGAAAGGCTTGCTGATAGCGCGCGGTCAGCGATGTTCCGGCGACAAGGGTGATGGCGTGTTGATGGGCCACATAGTCGCGCATACTGGCGACCTCTGCGCCAATCAACAAACAAGAGAGAAATTCGCTGACCTGTTCGCGGGGAAGTGTTCACAGCACATNCGAGGCGCGAACTTCAAAAAGATGCGGCAATATGGCGGGCGTAGTAAAACCACg                                                                                                >  2:636297/1‑300 (MQ=255)                                                                                                                                                                                                      gtttGCCACTGCATGAACGATGCTCCTTATACCAGCCTGAAATGCCGTGTCGCCCGCCACCGCCGTCACGTCGCAACCCATCGCCTGAAAGGCTTGCTGATAGCGCGCGGTCAGCGATGTTCCGGCGACAAGGGTGATGGCGTGTTGATGGGCCACATAGTCGCGCANACTGGCGACCTCTGCGCCAATCAACAAACCAGAGAGAAATTCGCTGACCTGTTCGCGGGGAAGTGTTCCCAGCACATGCGAGGCGCGAACTTCAAAAAGCTGCGGCAATATGGCGGGCGTATTAAAACCACg                                                                                                >  2:191600/1‑300 (MQ=255)                                                                                                                                                                                                                    gaacGATGCTCCTTATACCAGCCTGAAATGCCGTGTCGCCCGCCACCGCCGTCANGTCGCAACCCATCGCCTGAAAGGCTTGCTGATAGCGCGCGGTCAGCGATGTTCCTGCGACAGGGGTGATGGCGTGTTGATGGGCCACATAGTCGCGCATACTGGCGACCTCTGCGCCAATCAACAAACCAGAGAGAAATTCGCTGACCTGTTCGCGGGGAAGTGTTCCCAGCACATGCGAGGCGCGAANTTCAAAAAGCTGCGGCAATATGGCGGGCGTATTAAAACCACGCTCAAGGCCAGcgg                                                                                  <  2:795558/297‑3 (MQ=255)                                                                                                                                                                                                                         aTGCTCCTTATACCAGCCTGAAATGCCGTGTCGCCCGCCACCGCCGTCACGTCGCAACCCATCGCCTGAAAGGCTTGCTGATAGCGCGCGGTCAGCGATGTTCCGGCGACAAGGTTGATGGCGTGTTGATGGGCCACATAGTCGCGCATACTGGCGACCTCTGCGCCAATCAACAAACCAGAGAGAAATTCGCTGACCTGTTCGCGGGGAAGTGTTCCCAGCACATGCGAGGCGCGAACTTCAAAAAGCTGCGGCAATATGGCGGGCGTATTAAAACCACGCTCAAGGCCAGCTGTGAAg                                                                             <  2:649261/300‑1 (MQ=255)                                                                                                                                                                                                                         aTGCTCCTTATACCAGCCTGAAATGCCGTGTCGCCCGCCACCGCCGTCACGTCGCAACCCATCGCCTGAAAGGCTTGCTGATAGCGCGCGGTCAGCGATGTTCCGGCGACAAGGGTGATGGCGTGTTGATGGGCCACATAGTCGCGCATACTGGCGACCTCTGCGCCAATCAACAAACCAGAGAGAAATTCGCTGACCTGTTCGCGGGGAAGTGTTCCCAGCACATGCGAGGCGCGAACTTCAAANAGCTGCGGCAATATGGCGGGCGTATTAAAACCACGCTCAAGGCCAGCTGTGAAg                                                                             >  2:641849/1‑300 (MQ=255)                                                                                                                                                                                                                         aTGCTCCTTATACCAGCCTGAAATGCCGTGTCACCCGCCACCGCCGTCACGTCGCAACCCATCGCCTGAAATGCTTGCTGATAGCGCGCGGTCAGCGATGTTCCGGCGACAAGGGTGATGGCGTGTTGATGTGCCACATAGTCGCGCATACTGGCGACCTCTGCGCCAATCAACAAACCATAGAGAAATTCGCTGACCTGTTCGCGGGGAAGTGTTCCCAGCACATGCGAGGCGCGAACTTCAAAAAGCTGCGGCAATATGGCGGGCGTATTAAAACCACGCTCAAGGCCAGCTGTGAAg                                                                             <  2:60620/300‑1 (MQ=255)                                                                                                                                                                                                                          tGCTCCTTATACCAGCCTGAAATGCCGTGTCGCCCGCCACCGCCGTCACGTCGCAACCCATCGCCTGAAAGGCTTGCTGATAGCGCGCGGTCAGCGATGTTCCGGCGACAAGGGTGATGGCGTGTTGATGGGCCACATAGTCGCGCATACTGGCGACCTCTGCGCCAATCAACAAACCAGAGAGAAATTCGCTGACCTGTTCGCGGGGAAGTGTTCCCAGCACATGCGAGGCGCGAACTTCAAAAAGCTGCGGCAATATGGCGGGCGTATTAAAACCACGCTCAAGGCCAGCTGTGAAgg                                                                            <  1:483139/300‑1 (MQ=255)                                                                                                                                                                                                                              ccTTATACCAGCCTGAAATGCCGTGTCGCCCGCCACCGCCGTCACGTCGCAACCCATCGCCTGAAAGGCTTGCTGATAGCGCGCGGTCAGCGATGTTCCGGCGACAAGGGTGATGGCGTGTTGATGGGCCACATAGTCGCGCATACTGGCGACCTCTGCGCCAATCAACAAACCAGAGAGAAATTCGCTGACCTGTTCGCGGGGAAGTGTTCCCAGCACATGCGAGGCGCGAACTTCAAAAAGCTGCGGCAATATGGCGGGCGTATTAAAACCACGCTCAAGGCCAGCTGTGAAGGCATc                                                                        >  1:17847/1‑300 (MQ=255)                                                                                                                                                                                                                                  ataCCAGCCTGAAATGCCGTGTCGCCCGCCACCGCCGTCACGTCGCAACCCATCGCCTGAAAGGCTTGCTGATAGCGCGCGGTCAGCGATGTTCCGGCGACAAGGGTGATGGCGTGTTGATGGGCCACATAGTCGCGCATACTGGCGACCTCTGCGCCAATCAACAAACCAGAGAGAAATTCGCTGACCTGTTCGCGGGGAAGTGTTCCCAGCACATGCGAGGCGCGAACTTCAAAAAGCTGCGGCAATATGGCGGGCGTATTAAAACCACGCTCAAGGCCAGCTGTGAAGGCATCGGCa                                                                    >  1:805560/1‑300 (MQ=255)                                                                                                                                                                                                                                       aGCCTGAAATGCCGTGTCGCCCGCCACCGCCGTCACGTCGCAACCCATCGCCTGAAAGGCTTGCTGATAGCGCGCGGTCAGCGATGTTCCGGCGACAAGGGTGATGGCGTGTTGATGGGCCACATAGTCGCGCATACTGGCGACCTCTGCGCCAATCAACAAACCAGAGAGAAATTCGCTGACCTGTTCGCGGGGAAGTGTTCCCAGCACATGCGAGGCGCGAACTTCAAAAAGCTGCGGCAATATGGCGGGCGTATTAAAACCACGCTCAAGGCCAGCTGTGAAGGCATCGGCAGAGtt                                                               >  1:675599/1‑300 (MQ=255)                                                                                                                                                                                                                                             aaaTGCCGTGTCGCCCGCCACCGCCGTCACGTCGCAACCCATCGCCTTAAAGGCNTGCTGATAGCGCGCGGTCAGCGATGTTCCGGCGACAAGGGTGATGGCGTGTTGATGGGCCACATAGTCGCGCATACTGGCGACCTCTGCGCCAATCAACAAACCAGAGAGAAATTCGCTGACCTGTTCGCGGGGAAGTGTTCCCAGCACATGCGAGGCGCGAACTTCAAAAAGCTGCGGCAATATGGCNGGCGTATTAAAACCACGCTCAAGGCCAGCTGTGAAGGCATCGGCAGAGTTTTCCTg                                                         <  2:790545/300‑1 (MQ=255)                                                                                                                                                                                                                                                         gcccgccACCGCCGTCACGTCGCAACCCATCGCCTGAAAGGCTTGCTGATAGCGCGCGGTCAGCGATGTTCCGGCGACAAGGGTGATGGCGTGTTGATGGGCCACATAGTCGCGCATACTGGCGACCTCTGCGCCAATCAACAAACCAGAGAGAAATTCGCTGACCTGTTCGCGGGGAAGTGTTCCCAGCACATGCGAGGCGCGAACTTCAAAAAGCTGCGGCAATATGGCGGGCGTATTAAAACCACGCTCAAGGCCAGCTGTGAAGGCATCGGCAGAGTTTTCCTGCGGCGGCAAAcc                                             >  1:809435/1‑300 (MQ=255)                                                                                                                                                                                                                                                            cgccACCGCCGTCACGTCGCAACCCATCGCCTGAAAGGCTTGCTGATAGCGCGCGGTCAGCGATGTTCCGGCGAAAAGGGTGATGGCGTGTTGATGGGCCACATAGTCGCGCATACTGGCGACCTCTGCGCCAATCAACAAACCAGAGAGAAATTCGCTGACCTGTTCGCGGGGAAGTGTTCCCAGCACATGCGAGGCGCGAACTTCAAAAAGCTGCGGCAATATGGCGGGCGTATTAAAACCACGCTCAAGGCCAGCTGTGAAGGCATCGGCAGAGTTTTCCTGCGGCGGCAAACCTgc                                          >  1:482535/1‑300 (MQ=255)                                                                                                                                                                                                                                                                          cGTCGCAACCCATCGCCTGAAAGGCTTGCTGATAGCGCGCGGTCAGCGATGTTCCGGCGACAAGGGTGATGGCGTGTTGATGGGCCACATAGTCGCGCATACTGGCGACCTCTGCGCCAATCAANAAACCAGAGAGAAATTCGCTGACCTGTTCGCGGGGAAGTGTTCCCAGCACATGCGAGGCGCNAACTTCAAAAAGCTGCGGCAATATGGCGGGCGTATTAAAAACACGCTCAAGGCCAGCTGTGAAGGCATCGGCAGAGTTTTCCTGCGGCGGCAAACCTGCGCCAATCAATGAGt                            >  2:675075/1‑300 (MQ=255)                                                                                                                                                                                                                                                                          cGTCGCAACCCATCGCCTGAAAGGCTTGCTGATAGCGCGCGGTCAGCGATGTTCCGGCGACAAGGGTGATGGCGTGTTGATGGGCCACATAGTCGCGCATACTGGCGACCTCTGCGCCAATCAACAAACCAGAGAGAAATTCGCTGACCTGTTCGCGGGGAAGTGTTCCCAGCACATGCGAGGCGCGAACTTCAAAAAGCTGCGGCAATATGGCGGGCGTATTAAAACCACGCTCAAGGCCAGCTGTGAAGGCATCGGCAGAGTTTTCCTGCGGCGGCAAACCTGCGCCAATCAATGAGt                            >  2:446835/1‑300 (MQ=255)                                                                                                                                                                                                                                                                            tCGCAACCCATCGCCTGAAAGGCTTGCTGATAGCGCGCGGTCAGCGATGTTCCGGCGACAAGGGTGATGGCGTGTTGATGGGCCACATAGTCGCGCATACTGGCGACCTCTGCGCCAATCAACAAACCAGAGAGAAATTCGCTGACCTGTTCGCGGGGAAGTGTTCCCAGCACATGCGAGGCGCGAACTTCAAAAAGCTGCGGCAATATGGCGGGCGTATTAAAACCACGCTCAAGGCCAGCTGTGAAGGCATCGGCAGAGTTTTCCTGCGGCGGCAAACCTGCGCCAATCAATGAGTGa                          >  2:500384/1‑300 (MQ=255)                                                                                                                                                                                                                                                                                         ccTGAAAGGCTTGCTGATAGCGCGCGGTCAGCGATGTTCCGGCGACAAGGGTGATGGCGTGTTGATGGGCCACATAGTCGCGCATACTGGCGACCTCTGCGCCAATCAACAAACCAGAGAGAAATTCGCTGACCTGTTCGCGGGGAAGTGTTCCCAGCACATGCGAGGCGCGAACTTCAAAAAGCTGCGGCAATATGGCGGGCGTATTAAAACCACGCTCAAGGCCAGCTGTGAAGGCATCGGCAGAGTTTTCCTGCGGCGGCAAACCTGCGCCAATCAATGAGTGATTTAACAGTAAat             >  1:190903/1‑300 (MQ=255)                                                                                                                                                                                                                                                                                            gAAAGGCTTGCTGATAGCGCGCGGTCAGCGATGTTCCGGCGAAAAGGGTGATGGCGTGTTGATGGGCCACATAGTCGCGCATACTGGCGACCTCTGCGCCAATCAACAAACCAGAGAGAAATTCGCTGACCTGTTCGCGGGGAAGTGTTCCCAGCACATGCGAGGCGCGAACTTCAAAAAGCTGCGGCAATATGGCGGGCGTATTAAAACCACGCTCAAGGCCAGCTGTGAAGGCATCGGCAGAGTTTTCCTGCGGCGGCAAAACTGCGCCAATCAATGAGTGATTTAACAGTAAatgat          >  2:450082/1‑300 (MQ=255)                                                                                                                                                                                                                                                                                                   ttGCTGATAGCACGCGGTCAGCGATGTTCCGGCGACAAGGGTGATGGCGTGTTGATGGGCCACATAGTCGCGCATACTGGCGACCTCTGCGCCAATCAACAAACCAGAGAGAAATTCGCTGACCTGTTCGCGTGGAAGTGTTCCCAGCACATGCGAGGCGCGAACTTCAAAAAGCTGCGGCAATATGGCGGGCGTATTAAAACCACGCTCAATGCCAGCTGTGAAGGCATCGGCAGAGTTTTCCTGCGGCGGCAAACCTGCGCCAATCAATGAGTGATTTAACAGTAAATGATGTAAtgc   <  2:355646/300‑3 (MQ=255)                                                                                                                                                                                                                                                                                                    tGCTGATAGCGCGCGGTCAGCGATGTTCCGGCGACAAGGGTGATGGCGTGTTGATGGGCCACATAGTCGCGCATACTGGCGACCTCTGCGCCAATCAAAAAACCAGAGAGAAATTCGCTGACCTGTTCGCGGGGAAGTGTTCCCAGCACATGCGAGGCGCGAACTTCAAAAAGCTGCGGCAATATGGCGGGCGTATTAAAACCACGCTCAAGGCCAGCTGTGAAGGCATCGGCAGAGTTTTCCTGCGGCGGCAAACCTGCGCCAATCAATGAGTGATTTAACAGTAAATGATGTAATTCa  >  2:417960/1‑300 (MQ=255)                                                                                                                                                                                                                                                                                                         |                                                                                                                                                                                                                                                                                                       CAGGTTTCAGTACCGTACCTGCGCCAATCAACGCCTTGTCGCCGTACGCATCAACGATGGCGGGAATGCTTTGCTCCCATTGTGGGGAATTCAGCGGGATTTCAACCGCGTCGAACCCGGCGTCAATCACCGCGCCAACATGCGCCAGCGCCTCGTCGGGCGTAATACCGCGCAAAATGGCGATCAGCGGGAGTTTAGTTTGCCACTGCATGAGCGATGCTCCTTATACCAGCCTGAAATGCCGTGTCGCCCGCCACCGCCGTCACGTCGCAACCCATCGCCTGAAAGGCTTGCTGGTAGCGCGCGGTCAGCGATGTTCCGGCGACAAGGGTGATGGCGTGTTGATGGGCCACATAGTCGCGCATACTGGCGACCTCTGCGCCAATCAACAAACCAGAGAGAAATTCGCTGACCTGTTCGCGGGGAAGTGTTCCCAGCACATGCGAGGCGCGAACTTCAAAAAGCTGCGGCAATATGGCGGGCGTATTAAGACCACGCTCAAGGCCAGCTGTGAAGGCATCGGCAGAGTTTTCCTGCGGCGGCAAACCTGCGCCAATCAATGAGTGATTTAACAGTAAATGATGTAATTCA  >  NC_000913/3873402‑3873992

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