Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	NC_000913	4,325,000	C→T	G56S (GGC→AGC)	phnC ←	phosphonate/phosphate ABC transporter ATP binding subunit

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	4,325,000	0	C	T	100.0%	224.6 / NA	62	G56S (GGC→AGC)	phnC	phosphonate/phosphate ABC transporter ATP binding subunit
Reads supporting (aligned to +/- strand):  ref base C (0/0);  new base T (44/18);  total (44/18)

GGGTGGAAACGCGCTGATGGGCAAAATGCACCATGCCAACGCGGGTCAGCGCCTGTAACGCGCGCTGTTTCTGCTCGCCGGTGAACCAGCTAAAACAGGTGCGCCAGAACGGCGTGCTGCCGAGCGCGCCAATCAGCACGTTCTCCAGTACGCTCAGGCGGTTCACCAGGTTGAATTGTTGGAATATGTAGCCGGTATGGGCGCGGCTTTTGCGGATATCGCGGGCCAGGCGGCCTTCGCGCTGGACTGTGCGGCCCAGCAGCTCGATATGGCTACCGACAGATTTATCGCCGGTAATCAAACCGCTTAAGTGACGTAAAAGGGTGGATTTTCCGGAACCCGACGGCCCAAGCAGAGCCACCATTTCACCGTGATGAATGTTCAGATCAACCGCATGCAGCGCCTGATGCTGATTGAAGGTTTTGGCGAGCTTCTCGACACGGATAATCGTTTGCATGATGCAGCCTCCCTAAAAAAGTGGCTCCATCGTGGCGGATTATTGTGACGATTTGGTTAATTAAAAGTGACTAACAGATGAAGAGTTAACGGGAATTCGATGACAGTGCGGGAGGGCGGGCCT  >  NC_000913/4324709‑4325288                                                                                                                                                                                                                                                                                                    |                                                                                                                                                                                                                                                                                                 gggTGGAAATGCGCTGATGGGCAAAATGCACCATGCCAACGCGGGTCAGCGCCTGTAACGCGCGCTGTTTCTGCTCGCCGGTGAACCAGCTAAAACAGGTGCGCCAGAACGGCGTGCTGCCGAGCGCGCCAATCAGCACGTTCTCCAGTACGCTCAGGCGGTTCACCAGGTTGAATTGTTGGAATATGTAGCCGGTATGGGCGCGGCTTTTGCGGATATCGCGGGCCAGGCGGCCTTCGCGCTGGACTGTGCGGCCCAGCAGCTCGATATGGCTACCGACAGATTTATCGCTGGTAATCa                                                                                                                                                                                                                                                                                          >  2:31267/1‑300 (MQ=255)                       aaaaTGCACCATGCCAACGCGGGTCAGCGCCTGTAACGCGCGCTGTTTCTGCTCGCCGGTGAACCAGCTAAAACAGGTGCGCCAGAACGGCGTGCTGCCGAGCGCGCCAATCAGCACGTTCTCCAGTACGCTCAGGCGGTTCACCAGGTTGAATTGTTGGAATATGTAGCCGGTATGGGCGCGGCTTTTGCGGATATCGCGGGCCAGGCGGCCTTCGCGCTGGACTGTGCGGCCCAGCAGCTCGATATGGCTACCGACAGATTTATCGCTGGTAATCAAACCGCTTAAGTGACGTAAAAg                                                                                                                                                                                                                                                                    <  1:32007/300‑1 (MQ=255)                                            ggTCAGCGCCTGTAACGCGCGCTGTTTCTGCTCGCCGGTGAACCAGCTAAAACAGGTGCGCCAGAACGGCGTGCTGCCGAGCGCGCCAATCAGCACGTTCTCCAGTACGCTCAGGCGGTTCACCAGGTTGAATTGTTGGAATATGTAGCCGGTATGGGCGCGGCTTTTGCGGATATCGCGGGCCAGGCGGCCTTCGCGCTGGACTGTGCGGCCCAGCAGCTCGATATGGCTACCGACAGATTTATCGCTGGTAATCAAACCGCTTAAGTGACGTAAAAGGGTGGATTTTCCGGAACCCGa                                                                                                                                                                                                                                               >  1:386227/1‑300 (MQ=255)                                                  cgcCTGTAACGCGCGCTGTTTCTGCTCGCCGGTGAACCAGCTAAAACAGGTGCGCCAGAACGGCGTGCTGCCGAGCGCGCCAATCAGCACGTTCTCCAGTACGCTCAGGCGGTTCACCAGGTTGAATTGTTGGAATATGTAGCCGGTATGGGCGCGGCTTTTGCGGATATCGCGGGCCAGGCGGCCTTCGCGCTGGACTGTGCGGCCCAGCAGCTCGATATGGCTACCGACAGATTTATCGCTGGTAATCAAACCGCTTAAGTGACGTAAAAGGGTGGATTTTCCGGAACCCGACGGccc                                                                                                                                                                                                                                         >  2:662992/1‑300 (MQ=255)                                                    ccTGTAACGCGCGCTGTTTCTGCTCGCCGGTGAACCAGCTAAAACAGGTGCGCCAGAACGGCGTGCTGCCGAGCGCGCCAATCAGCACGTTCTCCAGTACGCTCAGGCGGTTCACCAGGTTGAATTGTTGGAATATGTAGCCGGTATGGGCGCGGCTTTTGCGGATATCGCGGGCCAGGCGGCCTTCGCGCTGGACTGTGCGGCCCAGCAGCTCGATATGGCTACCGACAGATTTATCGCTGGTAATCAAACCGCTTAAGTGACGTAAAAGGGTGGATTTTCCGGAACCCGACGGCCCaa                                                                                                                                                                                                                                       >  1:179186/1‑300 (MQ=255)                                                                 cTGTTTCTGCTCGCCGGTGAACCAGCTAAAACAGGTGCGCCAGAACGGCGTGCTGCCGAGCGCGCCAATCAGCACGTTCTCCAGTACGCTCAGGCGGTTCACCAGGTTGAATTGTTGGAATATGTAGCCGGTATGGGCGCGGCTTTTGCGGATATCGCGGGCCAGGCGGCCTTCGCGCTGGACTGTGCGGCCCAGCAGCTCGATATGGCTACCGACAGATTTATCGCTGGTAATCAAACCGCTTAAGTGACGTAAAAGGGTGGATTTTCCGGAACCCGACGGCCCAAGCAGAGcc                                                                                                                                                                                                                               >  1:102888/1‑295 (MQ=255)                                                                   gTTTCTGCTCGCCGGTGAACCAGCTAAAACAGGTGCGCCAGAACGGCGTGCTGCCGAGCGCGCCAATCAGCACGTTCTCCAGTACGCTCAGGCGGTTCACCAGGTTGAATTGTTGGAATATGTAGCCGGTATGGGCGCGGCTTTTGCGGATATCGCGGGCCAGGCGGCCTTCGCGCTGGACTGTGCGGCCCAGCAGCTCGATATGGCTACCGACAGATTTATCGCTGGTAATCAAACCGCTTAAGTGACGTAAAAGGGTGGATTTTCCGGAACCCGACGGCCCAAGCAGAGCCACCAttt                                                                                                                                                                                                                        >  1:467983/1‑300 (MQ=255)                                                                   gTTTCTGCTCGCCGGTGAACCAGCTAAAACAGGTGCGCCAGAACGGCGTGCTGCCGAGCACGCCAATCAGCACGTTCTCCAGTACGCTCAGGCGGTTCACCAGGTTGAATTGTTGGAATATGTAGCCGGTATGGGCGCGGCTTTTGCGGATATCGCGGGCCAGGCGGCCTTCGCGCTGGACTGTGCGGCCCAGCAGCTCGATATGGCTACCGACAGATTTANCGCTGGTAATCAAACCGCTTAAGTGACGTAAAAGGGTGGATTTTCCGGAACCCGACGGCCCAAGCAGAGCCACCAttt                                                                                                                                                                                                                        >  1:487810/1‑300 (MQ=255)                                                                                        aGCTAAAACAGGTGCGCCAGAACGGCGTGCTGCCGAGCGCGCCAATCAGCACGTTCTCCAGTACGCTCAGGCGGTTCACCAGGTTGAATTGTTGGAATATGTAGCCGGTATGGGCGCGGCTTTTGCGGATATCGCGGGCCAGGCGGCCTTCGCGCTGGACTGTGCGGCCCAGCAGCTCGATATGGCTACCGACAGATTTATCGCTGGTAATCAAACCGCTTAAGTGACGTAAAAGGGTGGATTTTCCGGAACCCGACGGCCCAAGCAGAGCCACCATTTCACCGTGATGAATGTTCAGAt                                                                                                                                                                                                   >  2:408469/1‑300 (MQ=255)                                                                                        aGCTAAAACAGGTGCGCCAGAACGGCGTGCTGCCGAGCGCGCCAATCAGCACGTTCTCCAGTACGCTCAGGCGGTTCACCAGGTTGAATTGTTGGAATATGTAGCCGGTATGGGCGCGGCTTTTGCGGATATCGCGGGCCAGGCGGCCTTCGCGCTGGACTGTGCGGCCCAGCAGCTCGATATGGCTACCGACAGATTTATCGCTGGTAATCAAACCGCTTAAGTGACGTAAAAGGGTGGATTTTCCGGAACCCGACGGCCCAAGCAGAGCCACCATTTCAACGTGATGAATGTTCAGAt                                                                                                                                                                                                   >  2:31977/1‑300 (MQ=255)                                                                                         gCTAAAACAGGTGCGCCAGAACGGCGTGCTGCCGAGCGCGCCAATCAGCACGTTCTCCAGTACGCTCAGGCGGTTCACCAGGTTGAATTGTTGGAATATGTAGCCGGTATGGGCGCGGCTTTTGCGGATATCGCGGGCCAGGCGGCCTTCGCGCTGGACTGTGCGGCCCAGCAGCTCGATATGGCTACCGACAGATTTATCGCTGGTAATCAAACCGCTTAAGTGACGTAAAAGGGTGGATTTTCCGGAACCCGACGGCCCAAGCAGAGCCACCATTTCAACGTGATGAATGCTCAgctc                                                                                                                                                                                                  >  2:54109/1‑297 (MQ=255)                                                                                                      cgcCAGAACGGCGTGCTGCCGAGCGCGCCAATCAGCACGTTCTCCAGTACGCTCAGGCGGTTCACCAGGTTGAATTGTTGGAATATGTAGCCGGTATGGGCGCGGCTTTTGCGGATATCGCGGGCCAGGCGGCCTTCGCGCTGGACTGTGCGGCCCAGCAGCTCGATATGGCTACCGACAGATTTATCGCTGGTAATCAAACCGCTTAAGTGACGTAAAAGGGTGGATTTTCCGGAACCCGACGGCCCAAGCAGAGCCACCATTTCACCGTGATGAATGTTCAGATCAACCGCATGCAgc                                                                                                                                                                                     >  2:632100/1‑300 (MQ=255)                                                                                                      cgcCAGAACGGCGTGCTGCCGAGCGCGCCAATCAGCACGTTCTCCAGTACGCTCAGGCGGTTCACCAGGTTGAATTGTTGGAATATGTAGCCGGTATGGGCGCGGCTTTTGCGGATATCGCGGGCCAGGCGGCCTTCGCGCTGGACTGTGCGGCCCAGCAGCTCGATATGGCTACCGACAGATTTATCGCTGGTAATCAAACCGCTTAAGTGACGTAAAAGGGTGGATTTTCCGGAACCCGACGGCCCAAGCAGAGCCACCATTTCACCGTGATGAATGTTCAGATCAACCGCATGCAgc                                                                                                                                                                                     >  1:630098/1‑300 (MQ=255)                                                                                                      cgcCAGAACGGCGTGCTGCCGAGCGCGCCAATCAGCACGTTCTCCAGTACGCTCAGGCGGTTCACCAGGTTGAATTGTTGGAATATGTAGCCGGTATGGGCGCGGCTTTTGCGGATATCGCGGGCCAGGCGGCCTTCGCGCTGGACTGTGCGGCCCAGCAGCTCGATATGGCTACCGACAGATTTATCGCTGGTAATCAAACCGCTTAAGTGACGTAAAAGGGTGGATTTTCCGGAACCCGACGGCCCAAGCAGAGCAACCATTTCACCGTGATGAATGTTCAGATCAACCGCATGCAgc                                                                                                                                                                                     >  1:89921/1‑300 (MQ=255)                                                                                                               ggCGTGCTGCCGAGCGCGCCAATCAGCACGTTCTCCAGTACGCTCAGGCGGTTCACCAGGTTGAATTGTTGGAATATGNAGCCGGTATGGGCGCGGCTTTTGCGGATATCGCGGGCCAGGCGGCCTTCGCGCTGGACTGTGCGGCCCAGCAGCTCGATATGGCTACCGACAGATTTATCGCTGGTAATCAAACCGCTTAAGTGACGTAAAAGGGTGGATTTTCCGGAACCCGACGGCCCAAGCAGAGCCACCATTTCACCGTGATGAATGTTCAGATCAACCGCATGCAGCGCCTGATGc                                                                                                                                                                            <  1:136343/300‑1 (MQ=255)                                                                                                                                   aaTCAGCACGTTCTCCAGTACGCTCAGGCGGTTCACCAGGTTGAATTGTTGGAATATGTAGCCGGTATGGGCGCGGCTTTTGCGGATATCGCGGGCCAGGCGGCCTTCGCGCTGGACTGTGCGGCCCAGCAGCTCGATATGGCTACCGACAGATTTATCGCTGGTAATCAAACCGCTTAAGTGACGTAAAAGGGTGGATTTTCCGGAACCCGACGGCCCAAGCAGAGCCACCATTTCACCGTGATGAATGTTCAGATCAACCGCATGCAGCGCCTGATGCTGATTGAAGGTTTTGGCGa                                                                                                                                                         >  1:544696/1‑299 (MQ=255)                                                                                                                                   aaTCAGCACGTTCTCCAGTACGCTCAGGCGGTTCACCAGGTTGAATTGTTGGAATATGTAGCCGGTATGGGCGCGGCTTTTGCGGATATCGCGGGCCAGGCGGCCTTCGCGCTGGACTGTGCGGCCCAGCAGCTCGATATGGCTACCGACAGATTTATCGCTGGTAATCAAACCGCTTAAGTGACGTAAAAGGGTGGATTTTCCGGAACCCGACGGCCCAAGCAGAGCCAACATTTCACAGTGATGAATGTTCAGATCAACCGCATGCAGCGCCTGATGCTGATTGAAGGTTTTGGCGa                                                                                                                                                         >  1:544962/1‑299 (MQ=255)                                                                                                                                   aaTCAGCACGTTCTCCAGTACGCTCAGGCGGTTCACCAGGTTGAATTGTTGGAATATGTAGCCGGTATGGGCGCGGCTTTTGCGGATATCGCGGGCCAGGCGGCCTTCGCGCTGGACTGTGCGGCCCAGCAGCTCGATATGGCTACCGACAGATTTATCGCTGGTAATCAAACCGCTTAAGTGACGTAAAAGGGTGGATTTTCCGGAACCCGACGGCACAAGCAGAGCCACCATTTCACCGTGATGAATGTTCAGATCAACCGCATGCAGCGCCTGATGCTGCTTGAAGGTTTTGGCGa                                                                                                                                                         >  2:580775/1‑299 (MQ=255)                                                                                                                                   aaTCAGCACGTTCTCCAGTACGCTCAGGCGGTTCACCAGGTTGAATTGTTGGAATATGTAGCAGGTATGGGCGCGGCTTTTGCGGATATCGCGGGCAAGGCGGCCTTCGCGCTGGACTGTGCGGACCAGCAGCTCGATATGGCTACCGACAGATTTATCGCTGGTACTCAAACCGCTTAAGTGACGTCAAAGGGTGGATTTTCCGGAACCCGACGGCCCAAGCAGAGCCACCATTTCAACGTGATGAATTTTCAGATCAAACGCATTCAGCGCCTGATGTTGATTGAGGGTTTTGGCGa                                                                                                                                                         >  2:367969/1‑299 (MQ=255)                                                                                                                                            gTTCTCCAGTACGCTCAGGCGGTTCACCAGGTTGAATTGTTGGAATATGTAGCAGGTATGGGCGCGGCTTTTGCGGATATCGCGGGCCAGGCGGCCTTCGCGCTGGACTGTGCGGCCCAGCAGCTCGATATGGCTAACGACAGATTTATCGCTGGTAATCAAACCGCTTAAGTGACGTAAAAGGGTGGATTTTCCGGAACCCGACGGCCCAAGCAGAGCCACCATTTCACCGTGATGAATGTTCAGATCCACCGCATGCAGCGCCTGATGCTGATTGAAGGTTTTGGCGAGCTTCTCGac                                                                                                                                               >  1:221765/1‑300 (MQ=255)                                                                                                                                              tctcCAGTACGCTCAGGCGGTTCACCAGGTTGAATTGTTGGAATATGTAGCCGGTATGGGCGCGGCTTTTGCGGATATCGCGGGCCAGGCGGCCTTCGCGCTGGACTGTGCGGCCCAGCAGCTCGATATGGCTACCGACAGATTTATCGCTGGTAATCAAACCGCTTAAGTGACGTAAAAGGGTGGATTTTCCGGAACCCGACGGCCCAAGCAGAGCCACCATTTCACCGTGATGAATGTTCAGATCAACCGCATGCAGCGCCTGATGCTGATTGAAGGTTTTGGCGAGCTTCTCGacac                                                                                                                                             <  1:383447/300‑1 (MQ=255)                                                                                                                                               ctcCAGTACGCTCAGGCGGTTCACCAGGTTGAATTGTTGGAATATGTAGCCGGTATGGGCGCGGCTTTTGCGGATATCGCGGGCCAGGCGGCCTTCGCGCTGGACTGTGCGGCCCAGCAGCTCGATATGGCTACCGACAGATTTATCGCTGGTAATCAAACCGCTTAAGTGACGTAAAAGGGTGGATTTTCCGGAACCCGACGGCCCAAGCAGAGCCACCATTTCACCGTGATGAATGTTCAGATCAACCGCATGCAGCGCCTGATGCTGATTGAAGGTTTTGGCGAGCTTCTAGacac                                                                                                                                             >  2:263129/1‑299 (MQ=255)                                                                                                                                                  cAGTACGCTCAGGCGGTTCACCAGGTTGAATTGTTGGAGTATGTAGCCGGTATGGGCGCGGCTTTTGCGGATATCGCGGGCCAGGCGGCCTTCGCGCTGGACTGTGCGGCCCAGCAGCTCGATATGGCTACCGACAGATTTATCGCTGGTACTCAAACCGCTTAAGTGACGTAAAAGGGTGGATTTTCCGGCACCCGACGGCCCAAGCAGAGCCACCATTTCACCGTGATGAATGTTCAGATCAACCGCATGCAGCGCCTGATGCTGATTGAAGGTTTTGGCGAGTTTCTAGACACGGAt                                                                                                                                         >  2:140598/1‑300 (MQ=255)                                                                                                                                                       cGCTCAGGCGGTTCACCAGGTTGAATTGTTGGAATATGTAGCCGGTATGGGCGCGGCTTTTGCGGATATCGCGGGCCAGGCGGCCTTCGCGCTGGACTGTGCGGCCCAGCAGCTCGATATGGCTACCGACAGATTTATCGCTGGTAATCAAACCGCTTAAGTGACGTAAAAGGGTGGATTTTCCGGAACCCGACGGCCCAAGCAGAGCCACCATTTCACCGTGATGAATGTTCAGATCAACCGCATGCAGCGCCTGATGCTGATTGAAGGTTTTGGCGAGCTTCTCGACACGGATAAt                                                                                                                                      >  1:555910/1‑298 (MQ=255)                                                                                                                                                         cTCAGGCGGTTCACCAGGTTGAATTGTTGGAATATGTAGCCGGTATGGGCGCGGCTTTTGCGGATATCGCGGGCCAGGCGGCCTTCGCGCTGGACTGTGCGGCCCAGCAGCTCGATATGGCTACCGACAGATTTATCGCTGGTAATCAAACCGCTTAAGTGACGTAAAAGGGTGGATTTTCCGGAACCCGACGGCCCAAGCAGAGCCACCATTTCACCGTGATGAATGTTCAGATCAACCGCATGCAGCGCCTGATGCTGATTGAAGGTTTTGGCGAGCTTCTCGACACGGATAATCGtt                                                                                                                                  >  1:625959/1‑300 (MQ=255)                                                                                                                                                                 gTTCACCAGGTTGAATTGTTGGAATATGTAGCCGGTATGGGCGCGGCTTTTGCGGATATCGCGGGCCAGGCGGCCTTCGCGCTGGACTGTGCGGCCCAGCAGCTCGATATGGCTACCGACAGATTTATCGCTGGTAATCAAACCGCTTAAGTGACGTAAAAGGGTGGATTTTCCGGAACCCGACGGCCCAAGCAGAGCCACCATTTCACCGTGATGAATGTTCAGATCAACCGCATGCAGCGCCTGATGCTGATTGAAGGTTTTGGCGAGCTTCTCGACACGGATAATCGTTTGCatgat                                                                                                                          >  1:415277/1‑300 (MQ=255)                                                                                                                                                                        aGGTTGAATTTTTGGAATATGTAGCCGGTATGGGCGCTGCTTTTGCGGATATCGCGGGCCAGGCGGCCTTCGCGCTGGACTGTGCGGCCCAGCAGCTCGATATGGCTACCGACAGATTTATCGCTGGTAATCAAACCGCTTAAGTGACGTAAAAGGGTGGATTTTCCGGAACCCGACGGCCCAAGCAGAGCCACCATTTCACCGTGATGAATGTTCAGATCAACCGCATGCAGCGCCTGATGCTGATTGAAGGTTTTGGCGAGCTTCTCGACACGGATAATCGTTTGCATGATGCAGCCt                                                                                                                   <  2:413785/300‑1 (MQ=255)                                                                                                                                                                        aGGTTGAATTGTTGGAATATGTAGCCGGTATGGGCGCGGCTTTTGCGGATATCGCGGGCCAGGCGGCCTTCGCGCTGGACTGTGCGGCCCAGCAGCTCGATATGGCTACCGACAGATTTATCGCTGGTAATCAAACCGCTTAAGTGACGTAAAAGGGTGGATTTTCCGGAACCCGACGGCCCAAGCAGAGCCACCATTTCACCGTGATGAATGTTCAGATCAACCGCATGCAGCGCCTGATGCTGATTGAAGGTTTTGGCGAGCTTCTCGACACGGATAATCGTTTGCATGATGCAGCCt                                                                                                                   <  2:413532/300‑1 (MQ=255)                                                                                                                                                                                 tgttgGAATATGTAGCCGGTATGGGCGCGGCTTTTGCGGATATCGCGGGCCAGGCGGCCTTCGCGCTGGACTGTGCGGCCCAGCAGCTCGATATGGCTACCGACAGATTTATCGCTGGTAATCAAACCGCTTAAGTGACGTAAAAGGGTGGATTTTCCGGAACCCGACGGCCCAAGCAGAGCCACCATTTCACCGTGATGAATGTTCAGATCAACCGCATGCAGCGCCTGATGCTGATTGAAGGTTTTGGCGAGCTTCTCGACACGGATAATCGTTTGCATGATGCAGCCTCCCTaaaaa                                                                                                          <  1:385671/300‑1 (MQ=255)                                                                                                                                                                                        atctGTAGCCGGTATGGGCGCGGCTTTTGCGGATATCGCGGGCCAGGCGGCCTTCGCGCTGGACTGTGCGGCCCAGCAGCTCGATATGGCTACCGACAGATTTATCGCTGGTAATCAAACCGCTTAAGTGACGTAAAAGGGTGGATTTTCCGGAACCCGACGGCCCAAGCAGAGCCACCATTTCACCGTGATGAATGTTCAGATCAACCGCATGCAGCGCCTGATGCTGATTGAAGGTTTTGGCGAGCTTCTCGACACGGATAATCGTTTGCATGATGCAGCCTCCCTAAAAAAGTGGCt                                                                                                   >  1:531931/4‑300 (MQ=255)                                                                                                                                                                                                   gTATGGGCGCGGCTTTTGCGGATATCGCGGGCCAGGCGGCCTTCGCGCTGGACTGTGCGGCCCAGCAGCTCGATATGGCTACCGACAGATTTATCGCTGGTAATCAAACCGCTTAAGTGACGTAAAAGGGTGGATTTTCCGGAACCCGACGGCCCAAGCAGAGCCACCATTTCACCGTGATGAATGTTCAGATCAACCGCATGCAGCGCCTGATGCTGATTGAAGGTTTTGGCGAGCTTCTCGACACGGATAATCGTTTGCATGATGCAGCCTCCCTAAAAAAGTGGCTCAATCGt                                                                                            >  2:607648/1‑296 (MQ=255)                                                                                                                                                                                                   gTATGGGCGCGGCTTTTGCGGATATCGCGGACCAGGCGGCCTTCGCGCTGGACTGTGCGGCCCAGCAGCTCGATATGGCTACCGAAAGATTTATCGCTGGTAATCAAACCGCTTAAGTGACGTAAAAGGGTGGATTTTCCGGAACCCGACGGCCCAAGCAGAGCCACCATTTCACCGTGATGAATATTCAGATCAACCGCATGCAGCGCATGATGCTGATTGAAGGTTTTGGCGAGCTtct                                                                                                                                                   >  2:156949/1‑241 (MQ=255)                                                                                                                                                                                                           gcggcgTTTGCGGATATCGCGGGCCAGGCGGCCTTCGCGCTGGACTGTGCGGCCCAGCAGCTCGATATGGCTACCGACAGATTTATCGCTGGTAATCAAACCGCTTAAGTGACGTAAAAGGGTGGATTTTCCGGAACCCGACGGCCCAAGCAGAGCCACCATTTCACCGTGATGAATGTTCAGACCAACCGCATGCAGCGCCTGATGCTGATTGAAGGTTTTGGCGAGCTTCTCGACACGGATAATCGTTTGCATGATGCAGCCTCCCTAAAAAAGTGGCTCCATCGTGGCGGATTATtg                                                                                <  2:636945/300‑1 (MQ=255)                                                                                                                                                                                                              gCTTTTGCGGATATCGCGGGCCAGGCGGCCTTCGCGCTGGACTGTGCGGCCCAGCAGCTCGATATGGCTACCGACAGATTTATCGCTGGTAATCAAACCGCTTAAGTGACGTAAAAGGGTGGATTTTCCGGAACCCGACGGCCCAAGCAGAGCCACCATTTCACCGTGATGAATGTTCAGATCAACCGCATGCAGCGCCTGATGCTGATTGAAGGTTTTGGNGAGCTTCTCGACACGGATAATCGTTTGCATGATGCAGCCTCCCTAAAAAAGTGGCTCCATCGTGGCGGATTATTGTGa                                                                             >  1:504056/1‑300 (MQ=255)                                                                                                                                                                                                                       ctgcgagatgtgtataagagacaggcgcTTGACTGTGCGGCCCAGCAGCTCGATATGGCTACCGACAGATTTATCGCTGGTAATCACACCGCTTAAGTGACGTAAAAGGGTGGATTTTCCGGAACCCGACGGCCCAAGCAGAGCCACCATTTCCCCGTGATGAATGTTCAGATCAACCGCATGCAGCGCATGATGCTGATTGAAGGTTTTGGCGAGCTTCTCGCCACGGATAATCGTTTGCATGATGCAGCCTACCTAAAAAAGTGGCTCCATCGTGGCGGATTATTGTGACGATTTGGt                                                                    >  1:452709/25‑300 (MQ=255)                                                                                                                                                                                                                           tCGCGGGCCAGGCGGCCTTCGCGCTGGACTGTGCGGCCCAGCAGCTCGATATGGCTACCGACAGATTTATCGCTGGTAATCAAACCGCTTAAGTGACGTAAAAGGGTGGATTTTCCGGAACCCGACGGCCCAAGCAGAGCCACCATTTCACCGTGATGAATGTTCAGATCAACCGCATGCAGCGCCTGATGCTGATTGAAGGTTTTGGCGAGCTTCTCGACACGGACAATCGTTTGCATGATGCAACCCCCCCAAAAAAACGACTCCATCgcgacaaaccaccataacaacccaaccaac                                                                >  1:615797/1‑271 (MQ=255)                                                                                                                                                                                                                                 gCCAGGCGTCCTTCGCGCTGGACTGTGCGGCACAGCAGCTCGATATGGCTACCGACAGATTTATCGCTGGTAATCAAACCGCTTAAGTGACGTAAAAGGGTGGATTTTCCGGAACCCGACGGCCCAAGCAGAGCCACCATTTCACCGTGATGAATGTTCAGATCNACCGCATGCAGCGCCTGATGCTGATTGAAGGTTTTGGCGAGCTTCTCGACACGGATAATCGTTTGCATGATGCAGCCTCCCTAAAAAAGTGGCTCCATCGTGGCGGATTATTGTGACGATTTGGTTAATTAAAAg                                                          <  2:485827/300‑1 (MQ=255)                                                                                                                                                                                                                                           cTTCGCGCTGGACTGTGTGGCCCAGCAGCTCGATATGGCTACCGACAGATTTATCGCTGGTAATCAAACCGCTTAAGTGACGTAAAAGGGTGGATTTTCCGGAACCCGACGGCCCAAGCAGAGCCACCATTTCACCGTGATGAATGTTCAGATCAACCGCATGCAGCGCCTGATGCTGATTGAAGGTTTTGGCGAGCTTCTCGACACGGATAATCGTTTGCATGATGCAGCCTCCATAAAAAAGTGGCTCCATCGTGGCGGATTATTGTGACGATTTGGTTAATTAAAAGTGACTaaaa                                                 >  2:390329/1‑297 (MQ=255)                                                                                                                                                                                                                                            ttCGCGCTGGACTGTGCGGCCCAGCAGCTCGATATGGCTACCGACAGATTTATCGCTGGTAATCAAACCGCTTAAGTGACGTAAAAGGGTGGATTTTCCGGAACCCGACGGCCCAAGCAGAGCCACCATTTCACCGTGATGAATGTTCAGATCAACCGCATGCAGCGCCTGATGCTGATTGAAGGTTTTGGCGAGCTTCTCGACACGGATAATCGTTTGCATGATGCAGCCTCCCTAAAAAAGTGGCTCCATCGTGGCGGATTATTGTGACGATTTGGTTAATTAAAAGTGACTAACAg                                                <  1:649227/299‑1 (MQ=255)                                                                                                                                                                                                                                               gcgcTGGACTGTGCGGNCCAGCAGCTCGATATGGCTACCGACAGATTTATCGCTGGTAATCAAACCGCTTAAGTGACGTAAAAGGGTGGATTTTCCGGAACCCGACGGCCCAAGCAGAGCCACCATTTCACCGTGATGAATGTCCAGATCAACCGCATGCAGCGCCTGATGCTGATTGAAGGTTTTGGCGAGCTTCTCGACACGGATAATCGTTTGCATGAGGCAGCCTCCCTAAAAAAGTGGCTCCATCGTGGCGGATTATTNTGACGATTTGGTTAATTANAAGTGACTAACAGATGa                                            >  1:578192/1‑300 (MQ=255)                                                                                                                                                                                                                                               gcgcTGGACTGTGCGGCCCAGCAGCTGGATATGGCTACCGACAGATTTATCGCTGGTAATCAAACGGCTTAAGTGACGTAAAAGGGTGGATTTTCCGGAACCCGACGGCCCAAGCAGAGCCACCATTTCACCGTGATGAATGTTCAGATCAACCGCATGCAGCGCCTGATGCTGATTGAAGGTTTTGGCGAGCTTCTCGACACGGATAATAGTTTGCATGATGCAGCATACATAAAAAAGTGGCTCCATCGTGGCGGATTGTTGTGACGATTTGGTTAATTAAAAGTGACTAACAGATGa                                            >  2:2115/1‑300 (MQ=255)                                                                                                                                                                                                                                               gcgcTGGACTGTGCGGCCCAGCAGCTCGATATGGCTACCGACAGATTTATCGCTGGTAATCAAACCGCTTAAGTGACGTAAAAGGGTGGATTTTCCGGAACCCGACGGCCCAAGCAGAGCCACCATTTCACCGTGATGAATGTTCAGATCAACCGCATGCAGCGCCTGATGCTGATTGAAGGTTTTGGCGAGCTTCTCGACACGGATAATCGTTTGCATGATGCAGCCTCCCTAAAAAAGTGGCTCCATCGTGGCGGATTATTGTGACGATTTGGTTAATTAAAAGTGACTAACAGATGa                                            >  2:286524/1‑300 (MQ=255)                                                                                                                                                                                                                                               gcgcTGGACTGTGCGGCCCAGCAGCTCGATATGGCTACCGACAGATTTATCGCTGGTAATCAAACCGCTTAAGTGACGTAAAAGGGTGGATTTTCCGGAACCCGACGGCCCAAGCAGAGCCACCATTTCACCGTGATGAATGTTCAGATCAACCGCATGCAGCGCCTGATGCTGATTGAAGGTTTTGGCGAGCTTCTCGACACGGATAATCGTTTGCATGATGCAGCCTCCCTAAAAAAGTGGCTCCATCGTGGCGGATTATTGTGACGATTTGGTTAATTAAAAGTGACTAACAGATGa                                            >  1:87090/1‑300 (MQ=255)                                                                                                                                                                                                                                               gcgcTGGACTGTGCGGCCCAGCAGCTCGATATGGCTACCGACAGATTTATCGCTGGTAATCAAACCGCTTAAGTGACGTAAAAGGGTGGATTTTCCGGAACCCGACGGCCCAAGCAGAGCCACCATTTCACCGTGATGAATGTTCAGATCAACCGCATGCAGCGCCTGATGCTGATTGAAGGTTTTGGCGAGCTTCTCGACACGGATAATCGTTTGCATGATGCAGCCTCCCTAAAAAAGTGGCTCCATCGTGGCGGATTATTGTGAAGATTTGGTTACTTAAAAGTGACTAACAGATGa                                            >  1:647794/1‑300 (MQ=255)                                                                                                                                                                                                                                               gcgcTGGACTGTGCGGCCCAGCAGCTCGATATGGCTACCGACAGATTTATCGCTGGTAATCAAACCGCTTAAGTGACGTAAAAGGGTGGATTTTCCGGAACCCGACGGCCCAAGCAGAGCCACCATTTCACCGTGATGAATGTTCAGATCAACCGCATGCAGCGCCTGATGCTGATTGAAGGTTTTGGCGAGCTTCTCGAAACGGATAATCGTTTGCATGATGCAGCCTCCCTAAAAAAGTGGCTCCATCGTGGCGGATTATTGTGACGATTTGGTTAATTAAAAGTGACTAACAGATGa                                            >  2:549588/1‑300 (MQ=255)                                                                                                                                                                                                                                                   tGGACTGTGCGGCCCAGCAGCTCGATATGGCTACCGACAGATTTATCGCTGGTAATCAAACCGCTTAAGTGACGTAAAAGGGTGGATTTTCCGGAACCCGACGGCCCAAGCAGAGCCACCATTTCACCGTGATGAATGTTCAGATCAACCGCATGCAGCGCCTGATGCTGATTGAAGGTTTTGGCGAGCTTCTCGACACGGATAATCGTTTGCATGATGCAGCCTCCCTAAAAAAGTGGCTCCATCGTGGCGGATTATTGTGACGATTTGGTTAATTAAAAGTGACTAACAGATGAagag                                        <  1:263129/300‑1 (MQ=255)                                                                                                                                                                                                                                                    ggACTGTGCGGCCCAGCAGCTCGATATGGCTACCGACAGATTTATCGCTGGTAATCAAACCGCTTAAGTGACGTAAAAGGGTGGATTTTCCGGAACCCGACGGCCCAAGCAGAGCCACCATTTCACCGTGATGAATGTTCAGATCAACCGCATGCAGCGCCTGATGCTGATTGAAGGTTTTGGCGAGCTTCTCGACACGGATAATCGTTTGCATGATGCAGCCTCCCTAAAAAAGTGGCTCCATCGTGGCGGATTATTGTGACGATTTGGTTAATTAAAAGTGACTAACAGATGAAGAGt                                       >  2:474160/1‑300 (MQ=255)                                                                                                                                                                                                                                                    ggACTGTGCGGCCCAGCAGCTCGATATGGCTACCGACAGATTTATCGCTGGTAATCAAACCGCTTAAGTGACGTAAAAGGGTGGATTTTCCGGAACCCGACGGCCCAAGCAGAGCCACCATTTCACCGTGATGAATGTTCAGATCAACCGCATGCAGCGCCTGATGCTGATTGAAGGTTTTGGCGAGCTTCTCGACACGGATAATCGTTTGCATGATGCAGCCTACCTAAAAAAGTGGCTCCATCGTGGCGGATTATTGTGACGATTTGGTTAATTAAAAGTGACTAACAGATGAAGAGt                                       >  2:474166/1‑300 (MQ=255)                                                                                                                                                                                                                                                    ggACTGTGCGGCCCAGCAGCTCGATATGGCTACCGACAGATTCATCGCTGGTAATCAAACCGCTTAAGTGACGTAAAAGGGTGGATTTTCCGGAACCCGACGGCCCAAGCAGAGCCACCATTTCACCGTGATGCATGTTCAGATCAACCGCATGCAGCGCCTGATGCTGATTAAAGGTTTTGGCAAGCTTCTCGACACGGCTAATCGTTTGCATGATGCAGCCTCCCTAAAAACGTGACTCCATCGTGGCGGATTATGGTGACGATTTGGTTAATTAAAAGTGACTCACATATGCAGAGt                                       >  2:52991/1‑300 (MQ=255)                                                                                                                                                                                                                                                       cTGTGCGGCCCAGCAGCTCGATATGGCTACCGACAGATTTATCGCTGGTAATCAAACCGCTTAAGTGACGTAAAAGGGTGGATTTTCCGGAACCCGACGGCCCAAGCAGAGCCACCATTTCACCGTGATGAATGTTCAGATCAACCGCATGCAGCGCCTGATGCTGATTGAAGGTTTTGGCGAGCTTCTCGACACGGATAATCGTTTGCATGATGCAGCCTCCATAAAAAAGTGGCTCCATCGTGGCGGATTATTGTGACGATTTGGTTAATTAAAAGTGACTAACAGATGAAGAGTTaa                                    >  2:193908/1‑300 (MQ=255)                                                                                                                                                                                                                                                         gtgCGGCCCAGCAGCTCGATATGGCTACCGACAGATTTATCGCTGGTAATCAAACCGCTTAAGTGACGTAAAAGGGTGGATTTTCCGGAACCCGACGGCCCAAGCAGAGCCACCATTTCACCGTGATGAATGTTCAGATTAACCGCATGCAGCGCCTGATGCTGATTGAAGGTTTTGGCGAGCTTCTCGACACGGATAATCGTTTGCATGATGCAGCCTACCTAAAAAAGTGGCTCCATCGTGGCGGATTATTGTGACGATGTGGTTAATTAAAAGTGACTAACAGATGAAGAGTTAAc                                   >  2:336515/1‑299 (MQ=255)                                                                                                                                                                                                                                                         gtgCGGCCCAGCAGCTCGATATGGCTACCGACAGATTTATCGCTGGTAATCAAACCGCTTAAGTGACGTAAAAGGGTGGATTTTCCGGAACCCGACGGCCCAAGCAGAGCCACCATTTCACCGTGATGAATGTTCAGATCAACCGCATGCAGCGCCTGATGCTGATTGAAGGTTTTGGCGAGCTTCTCGACACGGATAATCGTTTGCATGATGCAGCCTCCCTAAAAAAGTGGCTCCATCGTGGCGGATTATTGTGACGATTTGGTTAATTAAAAGTGACTAACAGATGAAGAGTTAAc                                   >  1:203772/1‑299 (MQ=255)                                                                                                                                                                                                                                                         gtgCGGCCCAGCAGCTCGATATGGCTACCGACAGATTTATCGCTGGTAATCAAACCGCTTAAGTGACGTAAAAGGGTGGATTTTCCGGAACCCGACGGCCCAAGCAGAGCCACCATTTCACCGTGATGAATGTTCAGATCAACCGCATGCAGCGCATGATGCTGATTGAAGGTTTTGGCGAGCTTCTCGACACGGATAATCGTTTGCATGATGCAGCCTCCCTAAAAAAGTGGCTCCATCGTGGCGGATTATTGTGACGATTTGGTTAATTAAAAGTGACTAACAGATGAAGAGTTAAc                                   >  1:8921/1‑299 (MQ=255)                                                                                                                                                                                                                                                         gtgCGGCCCAGCAGCTCGATATGGCTACCGACAGATTTATCGCTGGTAATCAAACCGCTTAAGTGACGTAAAAGGGTGGATTTTCCGGAACCCGACGGCACAAGCAGAGCCACCATTTCACCGTGATGAATGTTCAGATCAACCGCATGCAGCGCATGATGCTGATTGAAGGTTTTGGCGAGCTTCTCGACACGGATAATCGTTTGCATGATGCAGCCTACCTAAAAAAGTGGCTCCATCGTGGCGGATTATTGTGACGATTTGGTTAATTAAAAGTGACTAACAGATGAAGAGTTAAc                                   >  2:303960/1‑299 (MQ=255)                                                                                                                                                                                                                                                           gctgCCCAGCAGCTCGATATTGCTACCGACATATTTATCGCTGGTAATCAAACCGCGTAATTGACGTAAAAGGGTGGATTTTCTGGAACCCGACTGCCCAATCAGAGCCACCATTTCACCGTGATGAATGTTCATATCAACTGCATGCAGCGCCTGATGCTGATTGAAGGTTTTGGCGAGCTTCTCGACACGGATAATCGTTTGCATGATGCAGCCTCCCTAAAAAAGTGACTCCACCGTGGCGGATTATTTTGACGATTTGGTTAATTAAAATTGACTAACAGATGAAGAGTTAACggg                                <  2:128523/297‑1 (MQ=255)                                                                                                                                                                                                                                                                  agcagcTCGATATGGCTACCGACAGATTTATCGCTGGTAATCAAACCGCTTAAGTGACGTAAAAGGTTGGATTTTCCGGAACCCGACTGCCCAAGCAGAGCCACCATTTCACCGTGATGAATGTTCAGATCAACCGCATGCAGCGCCTGATGCTGATTGAAGGTTTTGGCGAGCTTCTCGACACGGATAATCGTTTGCATGATGCAGCCTCCCTAAAAAAGTGGCTCCATCGTGGCGGATTATTGTGACGATTTGGTTAATTAAAAGTGACTAACAGATGAAGAGTTAACggg                                <  2:128933/293‑1 (MQ=255)                                                                                                                                                                                                                                                                  agcagcGCGATATGGCTACCGACAGATTTATCGCTGGTAATCAAACCGCTTAAGTGACGTAAAAGTGTGGATTTTCCGGAACCCGACTGCCCAATCAGAGCCACCATTTCACCGTGATGAATGTTCAGATCAACCGCATGCAGCGCCTGATGCTGATTGAAGGTTTTGGCGAGCTTCTCGACACGGATAATCGTTTGCATGATGCAGCCTCCCTAAAAAAGTGGCTCCATCGTGGCGGATTATTGTGACGATTTGGTTAATTAAAAGTGACTAACAGATGAAGAGTTAACggg                                <  2:128921/293‑1 (MQ=255)                                                                                                                                                                                                                                                                      gcTCTATATGGCTACCGACAGATTTATCGCTAGTAATCAAACCGCTTAAGTGACGTAAAATGGTGGACTTTCCGGAACCCGACGGCCCAAGCAGAGCCACCATTTCACCGTGATGAATGTTCAGATCAACCGCATGCAGCGCTTGATGCTGATTGAAGGTTTTGGCGAGCTTCTCGACACGGATAATCGTTTGCATGATGCAGCCTCCCTAAAAAAGTGGCTCCATCGTGGCGAATTATTGTGACGATTTGGTTAATTAAAAGTGACTAACAGGTGGAGAGTTAACGGGAATTCgacgac                     <  2:569495/300‑1 (MQ=255)                                                                                                                                                                                                                                                                                gCTACCGACGGATTTATCGCTGGTAATCAAACCGCTTAAGTGACGTAAAAGGGTGGATTTTCCGGAACCCGACGGCCCAAGCAGAGCCACNATTTCACCGTGATGAATGTTCAGATCAACCGCATGCAGCGCCTGATGCTGATTGAAGGTTTTGGCGAGCTTCTTGACACGGATAATCGTTTGCATGATGCAGCCTCCCTAAAAAAGTGGCTCCATCGTGGCGGATTATTGTGACGATTTGGTTAATTAAAAGTGACTAACAGATGAAGAGTTAACGGGAATTCGATGACAGTGCGGGAg           <  2:538474/300‑1 (MQ=255)                                                                                                                                                                                                                                                                                gCTACCGACAGATTTATCGCTGGTAATCAACCCGCTTAATTGACGTAAAAGTGTGGATTTTCCGGAACCCGACGGCCCAAGCAGAGCCACCATTTCACCTTGATGAATGTTCAGATCAACCGCATGCAGCGCCTGATGCTGATTTAAGGTTTTGGCGAGCTTCTCGACACGGATAATCGTTTGCATGATGCAGCCTCCCTAAAAAAGTGGCTCCATCGTGGCGGATTATTGTGACGATTTGGTTAATTAAAAGGGACGAACAGATGAAGAGTTAACGGGAATTCGATGACAGTGCGGGAg           <  2:561741/300‑1 (MQ=255)                                                                                                                                                                                                                                                                                gCTACCGACAGATTTATCGCTGGTAATCAAACCGCTTAAGTGACGTAAAAGGGTGGATTTTCCGGAACCCGACGGCCCAAGCAGAGCCACCATTTCACCGTGATGAATGTTCAGATCAACCGCATGCAGCGCCTGATGCTGATTGAAGGTTTTGGCGAGCTTCTCGACACGGATAATCGTTTGCATGATGCAGCCTCCCTAAAAAAGTGGCTCCATCGTGGCGGATTATTGTGACGATTTGGTTAATTAAAAGTGACTAACAGATGAAGAGTTAACGGGAATTCGATGACAGTGCGGGAg           <  1:250147/300‑1 (MQ=255)                                                                                                                                                                                                                                                                                         aGATTTATCGCTGGTAATCAAACCGCTTAAGTGACGTAAAAGGGTGGATTTTCCGGAACCCGACGGCCCAAGCAGAGCCACCATTTCACCGTGATGAATGTTCAGATCAACCGCATGCAGCGCCTGATGCTGATTGAAGGTTTTGGCGAGCTTCTCGACACGGATAATCGTTTGCATGATGCAGCCTCCCTAAAAAAGTGGCTCCATCGTGGCGGATTATTGTGACGATTTGGTTAATTAAAAGTGACTAACAGATGAAGAGTTAACGGGAATTCGATGACAGTGCGGGAGGGCGGGCCt  <  1:31267/300‑1 (MQ=255)                                                                                                                                                                                                                                                                                                    |                                                                                                                                                                                                                                                                                                 GGGTGGAAACGCGCTGATGGGCAAAATGCACCATGCCAACGCGGGTCAGCGCCTGTAACGCGCGCTGTTTCTGCTCGCCGGTGAACCAGCTAAAACAGGTGCGCCAGAACGGCGTGCTGCCGAGCGCGCCAATCAGCACGTTCTCCAGTACGCTCAGGCGGTTCACCAGGTTGAATTGTTGGAATATGTAGCCGGTATGGGCGCGGCTTTTGCGGATATCGCGGGCCAGGCGGCCTTCGCGCTGGACTGTGCGGCCCAGCAGCTCGATATGGCTACCGACAGATTTATCGCCGGTAATCAAACCGCTTAAGTGACGTAAAAGGGTGGATTTTCCGGAACCCGACGGCCCAAGCAGAGCCACCATTTCACCGTGATGAATGTTCAGATCAACCGCATGCAGCGCCTGATGCTGATTGAAGGTTTTGGCGAGCTTCTCGACACGGATAATCGTTTGCATGATGCAGCCTCCCTAAAAAAGTGGCTCCATCGTGGCGGATTATTGTGACGATTTGGTTAATTAAAAGTGACTAACAGATGAAGAGTTAACGGGAATTCGATGACAGTGCGGGAGGGCGGGCCT  >  NC_000913/4324709‑4325288

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: ‑