Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	1	1,455,839	2 bp→G	coding (22‑23/168 nt)	unknown ←	hypothetical protein

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	1	1,455,839	0	T	.	100.0%	97.5 / NA	23	coding (23/168 nt)	unknown	hypothetical protein
Reads supporting (aligned to +/- strand):  ref base T (0/0);  new base . (11/12);  total (11/12)
*	1	1,455,840	0	T	G	100.0%	73.8 / NA	23	N8H (AAC→CAC)	unknown	hypothetical protein
Reads supporting (aligned to +/- strand):  ref base T (0/0);  new base G (11/12);  total (11/12)

TAAAAGCTGTTGTTCCGCTTCTTCCATCCAAACTACGGTGTAATCCAATTTCATTTTAAGCCCAATACTTTTTCCAGGTTAGCTAATACTTCTTCATTGGATTTGACTTTAATTTCTCCGTTATCTAGCCGTGCAAAAAATAACGCTTCTTCGGAAAGCGTCACTTTTTGTGGTAGGTAGTCATCTTTACGATAGTTATTC‑TGTGTATTTAAGCTCATAATAATTTCCTCGCTAAAAAGTTAGTCATAATCATAGATAAAACCGACAACAATAGCAACTTATTAACATAAATATGCTAATAAGATTAAACACATAAAAATAAATATGTTTTTATGTGTAAATATTTATCAATATTTATTGATTTATGTTTTTATGTTTTGTATAATCT  >  1/1455644‑1456031                                                                                                                                                                                                    |                                                                                                                                                                                                  tAAAAGCTGCTGCTCCGCTTCTTCCGTCCAAACCACGGTGTAATCCAGTTTCATTTTAAGTCCAATACTTTCTCCAGGTTAGCTAACACTTCTTCATTAGATTTAACTTTAATTTCTCCGTTATCTAACCGAGCAAAAAACAATACTTCTTCGGAAAGTGTCACTTTCTGCGGTAGATAGTCGGCTTTGCGGTAG‑GATTCTTGCATATTTAGGCTCATAATAATTTCCTCGCTAAAAAGTTAGTCATAATc                                                                                                                                           <  1:277438/251‑1 (MQ=255)       ctgctgctCCGCTTCTTCCGTCCAAACCACGGTGTAATCCAGTTTCATTTTAAGTCCAATACTTTCTCCAGGTTAGCTAACACTTCTTCATTAGATTTAACTTTAATTTCTCCGTTATCTAACCGAGCAAAAAACAATACTTCTTCGGAAAGTGTCACTTTCTGCGGTAGATAGTCGTCTTTGCGGTAG‑GATTCTTGCATATTTAGGCTCATAATAATTTCCTCGCTAAAAAGTTAGTCa                                                                                                                                                >  1:473074/8‑240 (MQ=255)              tCCGCTTCTTCCGTCCAAACCACGGTGTAATCCAGTTTCATTTTAAGTCCAATACTTTCTCCAGGTTAGCTAACACTTCTTCATTAGATTTAACTTTAATTTCTCCGTTATCTAACCGAGCAAAAAACAATACTTCTTCGGAAAGTGTCACTTTCTGCGGTAGATAGTCGTCTTTGCGGTAG‑GATTCTTGCATATTTAGGCTCATAATAATTTCCTCGCTAAAAAGTTAGTCATAATCATAGATAAAACCg                                                                                                                              <  1:859939/251‑1 (MQ=255)                cGCTTCTTCCGTCCAAACCACGGTGTAATCCAGTTTCATTTTAAGTCCAATACTTTCTCCAGGTTAGCTAACACTTCTTCATTAGATTTAACTTTAATTTCTCCGTTATCTAACCGAGCAAAAAACAATACTTCTTCGGAAAGTGTCACTTTCTGCGGTAGATAGTCGTCTTTGCGGTAG‑GATTCTTGCATATTTAGGCTCATAATAATTTCCTCGCTAAAAAGTTAGTCATAATc                                                                                                                                           >  1:889129/1‑236 (MQ=255)                        ccgtccaAACCACGGTGTAATCCAGTTTCATTTTAAGTCCAATACTTTCTCCAGGTTAGCTAACACTTCTTCATTAGATTTAACTTTAATTTCTCCGTTATCTAACCGAGCAAAAAACAATACTTCTTCGGAAAGTGTCACTTTCTGCGGTAGATAGTCGTCTTTGCGGTAG‑GATTCTTGCATATTTAGGCTCATAATAATTTCCTCGCTAAAAAGTTAGTCATAATCATAGATAAAACCGACAACAATAg                                                                                                                    <  1:500328/248‑1 (MQ=255)                           tccaAACCACGGTGTAATCCAGTTTCATTTTAAGTCCAATACTTTCTCCAGGTTAGCTAACACTTCTTCATTAGATTTAACTTTAATTTCTCCGTTATCTAACCGAGCAAAAAACAATACTTCTTCGGAAAGTGTCACTTTCTGCGGTAGATAGTCGTCTTTGCGGTAG‑GATTCTTGCATATTTAGGCTCATAATAATTTCCTCGCTAAAAAGTTAGTCATAATCATAGATAAAACCGACAACAATAGCaa                                                                                                                 <  1:716382/251‑1 (MQ=255)                              aaaaaacgaactttttctgtttCATTTTAAGTCCAATACTTTCTCCAGGTTAGCTAACACTTCTTCATTAGATTTAACTTTAATTTCTCCGTTATCTAACCGAGCAAAAAACAATACTTCTTCGGAAAGTGTCACTTTCTGCGGTAGATAGTCGTCTTTGCGGTAG‑GATTCTTGCATATTTAGGCTCATAATAATTTCCTCGCTAAAAAGTTAGTCATAATCATAGATAAAACCGACAACAATAGCAACt                                                                                                               <  1:723304/231‑1 (MQ=255)                                      gtgtAATCCAGTTTCATTTTAAGTCCAATACTTTCTCCAGGTTAGCTAACACTTCTTCATTAGATTTAACTTTAATTTCTCCGTTATCTAACCGAGCAAAAAACAATACTTCTTCGGAAAGTGTCACTTTCTGCGGTAGATAGTCGTCTTTGCGGTAG‑GATTCTTGCATATTTAGGCTCATAATAATTTCCTCGCTAAAAAGTTAGTCATAATCATAGAt                                                                                                                                     >  1:624756/1‑220 (MQ=255)                                      gtgtAATCCAGTTTCATTTTAAGTCCAATACTTTCTCCAGGTTAGCTAACACTTCTTCATTAGATTTAACTTTAATTTCTCCGTTATCTAACCGAGCAAAAAACAATACTTCTTCGGAAAGTGTCACTTTCTGCGGTAGATAGTCGTCTTTGCGGTAG‑GATTCTTGCATATTTAGGCTCATAATAATTTCCTCGCTAAAAAGTTAGTCATAATCATAGATAAAACCGACAACAATAGCAACTTATTAAc                                                                                                        >  1:69249/1‑249 (MQ=255)                                      gtgtAATCCAGTTTCATTTTAAGTCCAATACTTTCTCCAGGTTAGCTAACACTTCTTCATTAGATTTAACTTTAATTTCTCCGTTATCTAACCGAGCAAAAAACAATACTTCTTCGGAAAGTGTCACTTTCTGCGGTAGATAGTCGTCTTTGCGGTAG‑GATTCTTGCATATTTAGGCTCATAATAATTTCCTCGCTAAAAAGTTAGTCATAATCATAGATAAAACCGACAACAATAGCAACTTATTAACAt                                                                                                      >  1:284637/1‑251 (MQ=255)                                           aTCCAGTTTCATTTTAAGTCCAATACTTTCTCCAGGTTAGCTAACACTTCTTCATTAGATTTAACTTTAATTTCTCCGTTATCTAACCGAGCAAAAAACAATACTTCTTCGGAAAGTGTCACTTTCTGCGGTAGATAGTCGTCTTTGCGGTAG‑GATTCTTGCATATTTAGGCTCATAATAATTTCCTCGCTAAAAAGTTAGTCATAATCATAGATAAAACCGACAACAATAGCAACTTATTAACATAAata                                                                                                 <  1:94911/251‑1 (MQ=255)                                                     aTTTTAAGTCCAATACTTTCTCCAGGTTAGCTAACACTTCTTCATTAGATTTAACTTTAATTTCTCCGTTATCTAACCGAGCAAAAAACAATACTTCTTCGGAAAGTGTCACTTTCTGCGGTAGATAGTCGTCTTTGCGGTAG‑GATTCTTGCATATTTAGGCTCATAATAATTTCCTCGCTAAAAAGTTAGTCATAATCATAGATAAAACCGACAACAATAGCAACTTATTAACATAAATATGCTAATAAg                                                                                       <  1:587518/251‑1 (MQ=255)                                                            gtccAATACTTTCTCCAGGTTAGCTAACACTTCTTCATTAGATTTAACTTTAATTTCTCCGTTATCTAACCGAGCAAAAAACAATACTTCTTCGGAAAGTGTCACTTTCTGCGGTAGATAGTCGTCTTTGCGGTAG‑GATTCTTGCATATTTAGGCTCATAATAATTTCCTCGCTAAAAAGTTAGTCATAATCATAGATAAAACCGACAACAATAGCAACTTATTAACATAAATATGCTAATAAGATTaag                                                                                 >  1:752043/3‑249 (MQ=255)                                                                            aGGTTAGCTAACACTTCTTCATTAGATTTAACTTTAATTTCTCCGTTATCTAACCGAGCAAAAAACAATACTTCTTCGGAAAGTGTCACTTTCTGCGGTAGATAGTCGTCTTTGCGGTAG‑GATTCTTGCATATTTAGGCTCATAATAATTTCCTCGCTAAAAAGTTAGTCATAATCATAGATAAAACCGACAACAATAGCAACTTATTAACATAAATATGCTAATAAGATTAAGCATATaaa                                                                         <  1:654542/242‑1 (MQ=255)                                                                                tAGCTAACACTTCTTCATTAGATTTAACTTTAATTTCTCCGTTATCTAACCGAGCAAAAAACAATACTTCTTCGGAAAGTGTCACTTTCTGCGGTAGATAGTCGTCTTTGCGGTAG‑GATTCTTGCATATTTAGGCTCATAATAATTTCCTCGCTAAAAAGTTAGTCATAATCATAGATAAAACCGacaac                                                                                                                         >  1:734837/1‑190 (MQ=255)                                                                                       caCTTCTTCATTAGATTTAACTTTAATTTCTCCGTTATCTAACCGAGCAAAAAACAATACTTCTTCGGAAAGTGTCACTTTCTGCGGTAGATAGTCGTCTTTGCGGTAG‑GATTCTTGCATATTTAGGCTCATAATAATTTCCTCGCTAAAAAGTTAGTCATAATCATAGCTAAAACCGACAACAATAGCAACTTATTAACATAAATATGCTAATAAGATTAAGCATATAAAAACAAATATGTTTTTAtgcg                                                     >  1:441258/2‑249 (MQ=255)                                                                                         cttcttCATTAGATTTAACTTTAATTTCTCCGTTATCTAACCGAGCAAAAAACAATACTTCTTCGGAAAGTGTCACTTTCTGCGGTAGATAGTCGTCTTTGCGGTAG‑GATTCTTGCATATTTAGGCTCATAATAATTTCCTCGCTAAAAAGTTAGTCATAATCATAGATAAAACCGACAACAATAGCaa                                                                                                                 >  1:779827/1‑189 (MQ=255)                                                                                          ttcttcATTAGATTTAACTTTAATTTCTCCGTTATCTAACCGAGCAAAAAACAATACTTCTTCGGAAAGTGTCACTTTCTGCGGTAGATAGTCGTCTTTGCGGTAG‑GATTCTTGCATATTTAGGCTCATAATAATTTCCTCGCTAAAAAGTTAGTCATAATCATAGATAAAACCGACAACAATAGCAACTTATTAAc                                                                                                        <  1:651000/197‑1 (MQ=255)                                                                                              tcATTAGATTTAACTTTAATTTCTCCGTTATCTAACCGAGCAAAAAACAATACTTCTTCGGAAAGTGTCACTTTCTGCGGTAGATAGTCGTCTTTGCGGTAG‑GATTCTTGCATATTTAGGCTCATAATAATTTCCTCGCTAAAAAGTTAGTCATAATCATAGATAAAACCGACAACAATAGCAACTTATTAAc                                                                                                        >  1:298682/1‑193 (MQ=255)                                                                                                          aCTTTAATTTCTCCGTTATCTAACCGAGCAAAAAACAATACTTCTTCGGAAAGTGTCACTTTCTGCGGTAGATAGTCGTCTTTGCGGTAG‑GATTCTTGCATATTTAGGCTCATAATAATTTCCTCGCTAAAAAGTTAGTCATAATCATAGAt                                                                                                                                     <  1:715891/152‑1 (MQ=255)                                                                                                                         ttATCTAACCGAGCAAAAAACAATACTTCTTCGGAAAGTGTCACTTTCTGCGGTAGATAGTCGTCTTTGCGGTAG‑GATTCTTGCATATTTAGGCTCATAATAATTTCCTCGCTAAAAAGTTAGTCATAATc                                                                                                                                           >  1:736513/1‑131 (MQ=255)                                                                                                                                     gCAAAAAACAATACTTCTTCGGAAAGTGTCACTTTCTGCGGTAGATAGTCGTCTTTGCGGTAG‑GATTCTTGCATATTTAGGCTCATAATAATTTCCTCGCTAAAAAGTTAGTCATAATCATAGATAAAACCGACAACAATAGCAACTTATTAACATAAATATGCTAATAAGATTaagca                                                                               <  1:404797/179‑4 (MQ=255)                                                                                                                                          aaacaatacttcttCGGAAAGTGTCACTTTCTGCGGTAGATAGTCGTCTTTGCGGTAG‑GATTCTTGCATATTTAGGCTCATAATAATTTCCTCGCTAAAAAGTTAGTCATAATCATAGATAAAACCGACAACAATAGCAACTTATTAACATAAATATGCTAATAAGATTAAGCATATAAAAACAAATATGTTTTTATGCGTAAATATTTATCAATATTTATTGATTTATGTTTTTATGCTTTGTATAATCt  <  1:737277/243‑1 (MQ=255)                                                                                                                                                                                                    |                                                                                                                                                                                                  TAAAAGCTGTTGTTCCGCTTCTTCCATCCAAACTACGGTGTAATCCAATTTCATTTTAAGCCCAATACTTTTTCCAGGTTAGCTAATACTTCTTCATTGGATTTGACTTTAATTTCTCCGTTATCTAGCCGTGCAAAAAATAACGCTTCTTCGGAAAGCGTCACTTTTTGTGGTAGGTAGTCATCTTTACGATAGTTATTC‑TGTGTATTTAAGCTCATAATAATTTCCTCGCTAAAAAGTTAGTCATAATCATAGATAAAACCGACAACAATAGCAACTTATTAACATAAATATGCTAATAAGATTAAACACATAAAAATAAATATGTTTTTATGTGTAAATATTTATCAATATTTATTGATTTATGTTTTTATGTTTTGTATAATCT  >  1/1455644‑1456031

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