Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	AM260479	500,512	G→A	W23* (TGG→TAG)	h16_A0473 →	ABC‑type transporter, permease component: PAAT family

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	AM260479	500,512	0	G	A	100.0%	70.2 / NA	19	W23* (TGG→TAG)	h16_A0473	ABC‑type transporter, permease component: PAAT family
Reads supporting (aligned to +/- strand):  ref base G (0/0);  new base A (11/8);  total (11/8)

TGCCGTTGACAGCAGTGTGAAACGGAAGGCATGTCCTTCCGTTTCTTTTTGAGGACGCAACATGAACTACAACTGGCATTGGGCAGTTTTCCTCGAACAGGCCGCCCAGAACGAGACCTACCTGGACTGGATGATCTCGGGCCTGAAGGTAACGCTCGCCCTGGGACTCTCGTCCTGGGTCATCGCCCTGGTCATCGGTTCGGTGCTCGGCGTGCTGCGCACCGTGCCCAACAAGTGGCTGGCGGGCATCGCCGCCA  >  AM260479/500384‑500640                                                                                                                                 |                                                                                                                                 tGCCGTTGACAGCAGTGTGAAACGGAAGGCATGTCCTTCCGTTTCTTTTTGAGGACGCAACATGAACTACAACTGGCATTGGGCAGTTTTCCTCGAACAGGCCGCCCAGAACGAGACCTACCTGGACTAgatg                                                                                                                              >  1:436692/1‑133 (MQ=255) tGCCGTTGACAGCAGTGTGAAACGGAAGGCATGTCCTTCCGTTTCTTTTTGAGGACGCAACATGAACTACAACTGGCATTGGGCAGTTTTCCTCGAACAGGCCGCCCAGAACGAGACCTACCTGGACTAgatg                                                                                                                              <  2:436692/133‑1 (MQ=255)                  tgAAACGG‑AGGCATGTCCTTCCGTTTCTTTTTGAGGACGCAACATGAACTACAACTGGCATTGGGCAGTTTTCCTCGAACAGGCCGCCCAGAACGAGACCTACCTGGACTAGATGATCTCGGGCCTGAAGGTaa                                                                                                           <  1:281813/134‑1 (MQ=255)                      aCGGAAGGCATGTCCTTCCGTTTCTTTTTGAGGACGCAACATGAACTACAACTGGCATTGGGCAGTTTTCCTCGAACAGGCCGCCCAGAACGAGACCTACCTGGACTAGATGATCTCGGGCCTGAAGGTAACGc                                                                                                        >  1:647152/1‑134 (MQ=255)                              cATGTCCTTCCGTTTCTTTTTGAGGACGCAACATGAACTACAACTGGCATTGGGCAGTTTTCCTCGAACAGGCCGCCCAGAACGAGACCTACCTGGACTAGATGATCTCGGGCCTGAAGGTaa                                                                                                           <  2:75570/123‑1 (MQ=255)                              cATGTCCTTCCGTTTCTTTTTGAGGACGCAACATGAACTACAACTGGCATTGGGCAGTTTTCCTCGAACAGGCCGCCCAGAACGAGACCTACCTGGACTAGATGATCTCGGGCCTGAAGGTaa                                                                                                           >  1:75570/1‑123 (MQ=255)                                                tttGAGGACGCAACATGAACTACAACTGGCATTGGGCAGTTTTCCTCGAACAGGCCGCCCAGAACGAGACCTACCTGGACTAGATAATCTCGGGCCTGAAGGTAACGCTCGCCCTGGGACTCTCGTCCTGGGt                                                                               <  2:647152/133‑1 (MQ=255)                                                 ttGAGGACGCAACATGAACTACAACTGGCATTGGGCAGTTTTCCTCGAACAGGCCGCCCAGAACGAGACCTACCTGGACTAGATGATCTCGGGCCTGAAGGTAACGCTCGCCCTGGGACTCTCGTCCTGGGTCa                                                                             >  2:743963/1‑134 (MQ=255)                                                            aCATGAACTACAACTGGCATTGGGCAGTTTTCCTCGAACAGGCCGCCCAGAACGAGACCTACCTGGACTAGATGATCTCGGGCCTGa                                                                                                                 >  2:759747/1‑87 (MQ=255)                                                            aCATGAACTACAACTGGCATTGGGCAGTTTTCCTCGAACAGGCCGCCCAGAACGAGACCTACCTGGACTAGATGATCTCGGGCCTGa                                                                                                                 <  1:759747/87‑1 (MQ=255)                                                                                tGGGCAGTTTTCCTCGAACAGGCCGCCCAGAACGAGACCTACCTGGACTAGATGATCTCGGGCCTGAAGGTAACGCTCGcc                                                                                                   >  1:420846/1‑81 (MQ=255)                                                                                tGGGCAGTTTTCCTCGAACAGGCCGCCCAGAACGAGACCTACCTGGACTAGATGATCTCGGGCCTGAAGGTAACGCTCGcc                                                                                                   <  2:420846/81‑1 (MQ=255)                                                                                tGGGCAGTTTTCCTCGAACAGGCCGCCCAGAACGAGACCTACCTGGACTAGATGATCTCGGGCCTGAAGGTAACGCTCGCCCTGGGACTCTCGTCCTGGGTCATCGCCCTGGTCATCGGTTCGGTGCTCGGCGt                                              >  2:1008007/1‑134 (MQ=255)                                                                                  ggCAGTTTTCCTCGAACAGGCCGCCCAGAACGAGACCTACCTGGACTAGATGATCTCGGGCCTGAAGGTaa                                                                                                           >  2:251362/1‑71 (MQ=255)                                                                                  ggCAGTTTTCCTCGAACAGGCCGCCCAGAACGAGACCTACCTGGACTAGATGATCTCGGGCCTGAAGGTaa                                                                                                           <  1:251362/71‑1 (MQ=255)                                                                                       ttttCCTCGAACAGGCCGCCCAGAACGAGACCTACCTGGACTAGATGATCTCGGGCCTGAAGGTAACGCTCGCCCTGGGACTCTCGTCCTGGGTCATCGCCCTGGTCATCGGTTCGGTGCTCGGCGTGCTgcgc                                       <  1:743963/134‑1 (MQ=255)                                                                                                                       tacctGGACTAGATGATCTCGGGCCTGAAGGTAACGCTCGCCCTGGGACTCTCGTCCTGGGTCATCGCCCTGGTCATCGGTTCGGTGCTCGGCGTGCTGCGCACCGTGCCCAACAAGTGGCTGGCGGGCATcgc       >  1:284702/1‑134 (MQ=255)                                                                                                                        acctGGACTAGATGATCTCGGGCCTGAAGGTAACGCTCGCCCTGGGACTCTCGTCCTGGGTCATCGCCCTGGTCATCGGTTCGGTGCTCGGCGTGCTGCGCACCGTGCCCAACAAGTGGCTGGCGGGCATcgcc      >  2:258651/1‑134 (MQ=255)                                                                                                                        acctGGACTAGATGATCTCGGGCCTGAAGGTAACGCTCGCCCTGGGACTCTCGTCCTGGGTCATCGCCCTGGTCATCGGTTCGGTGCTCGGCGTGCTGCGCACCGTGCCCAACAAGTGGCTGGCGGGCATcgcc      >  2:924962/1‑134 (MQ=255)                                                                                                                            ggactagaTGATCTCGGGCCTGAAGGTAACGCTCGCCCTGGGACTCTCGTCCTGGGTCATCGCCCTGGTCATCGGTTCGGTGCTCGGCGTGCTGCGCACCGTGCCCAACAAGTGGCTGGCGGGCATCGCCGCCa  >  2:646517/1‑134 (MQ=255)                                                                                                                                 |                                                                                                                                 TGCCGTTGACAGCAGTGTGAAACGGAAGGCATGTCCTTCCGTTTCTTTTTGAGGACGCAACATGAACTACAACTGGCATTGGGCAGTTTTCCTCGAACAGGCCGCCCAGAACGAGACCTACCTGGACTGGATGATCTCGGGCCTGAAGGTAACGCTCGCCCTGGGACTCTCGTCCTGGGTCATCGCCCTGGTCATCGGTTCGGTGCTCGGCGTGCTGCGCACCGTGCCCAACAAGTGGCTGGCGGGCATCGCCGCCA  >  AM260479/500384‑500640

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