Predicted mutation
evidence	seq id	position	mutation	freq	annotation	gene	description
RA	AE015451	6,013,885	T→C	100%	intergenic (+90/+20)	xpt → / ← PP_5266	Xanthine phosphoribosyltransferase/putative acetyl‑CoA hydrolase family protein

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	AE015451	6,013,885	0	T	C	100.0%	64.3 / NA	19	intergenic (+90/+20)	xpt/PP_5266	Xanthine phosphoribosyltransferase/putative acetyl‑CoA hydrolase family protein
Reads supporting (aligned to +/- strand):  ref base T (0/0);  new base C (9/10);  total (9/10)
Rejected as polymorphism: Frequency below/above cutoff threshold.
Rejected as polymorphism: Variant not supported by required number of reads on each strand.

GTGCCACCGTGGCGGGCCTGGGTATCGTGATCGAGAAGTCGTTCCAGGGCGGCCGTGCCGAACTGGACAGCCAGGGTTACCGCGTTGAATCGCTGGCCCGGGTGAAGTCGCTGGAAGGTGGTGTGGTCAGCTTCATCGAATAAACGCGACACCTGCTTCGCGGGCACGCCCGCTCCCACAGGTGCAGCACAGGTCTCGAACGCTGTGCCGTACCTGTGGGAGCGGGCGTGCCTGCGAATGGCCGCAACCCTTTTATCTGGTTGCCGCCAGCCCCGCCAACAACAAGCGCTGGTAAAGTTCTTCTTTAAGCCCCTGCGGCTGCCCCAACTGCATCCGCCCCAGGTGGCTTTCATACCCCTCCGGCCCCGGTGCATCCAACGCTGCCTTGCCCAGCTCCAGCACCTCGCTCAGCTTGAACTTGCTCTTCAACCAGTTCAGCGCC  >  AE015451/6013653‑6014094                                                                                                                                                                                                                                         |                                                                                                                                                                                                                  gTGCCACCGTGGCGGGCCTGGGTATCGTGATCGAGAAGTCGTTCCAGGGCGGCCGTGCCGAACTGGACAGCCAGGGTTACCGCGTTGAATCGCTGGCCCGGGTGAAGTCGCTGGAAGGTGGTGTGGTCAGCTTCATCGAATAAACGCGACACCTGCTTCGCGGGCACGCCCGCTCCCACAGGTGCAGCACAGGTCTCGAACGCTGTGCCGTACCTGTGGGAGCGGGCGTGCCcgcg                                                                                                                                                                                                                <  1:488360/236‑1 (MQ=255)          tGGCGGGCCTGGGTATCGTGATCGAGAAGTCGTTCCAGGGCGGCCGTGCCGAACTGGACAGCCAGGGTTACCGCGTTGAATCGCTGGCCCGGGTGAAGTCGCTGGAAGGTGGTGTGGTCAGCTTCATCGAATAAACGCGACACCTGCTTCGCGGGCACGCCCGCTCCCACAGGTGCAGCACAGGTCTCGAACGCTGTGCCGTACCTGTGGGAGCGGGCGTGCCCGCGAATGGcc                                                                                                                                                                                                         <  2:416987/234‑1 (MQ=255)                            tGATCGAGAAGTCGTTCCAGGGCGGCCGTGCCGAACTGGACAGCCAGGGTTACCGCGTTGAATCGCTGGCCCGGGTGAAGTCGCTGGAAGGTGGTGTGGTCAGCTTCATCGAATAAACGCGACACCTGCTTCGCGGGCACGCCCGCTCCCACAGGTGCAGCACAGGTCTCGAACGCTGTGCCGTACCTGTGGGAGCGGGCGTGCCCGCGAATGGCCGCAACCCTTTTATCTGGTTg                                                                                                                                                                                     >  1:93929/1‑236 (MQ=255)                                               gggCGGCCGTGCCGAACTGGACAGCCAGGGTTACCGCGTTGAATCGCTGGCCCGGGTGAAGTCGCTGGAAGGTGGTGTGGTCAGCTTCATCGAATAAACGCGACACCTGCTTCGCGGGCACGCCCGCTCCCACAGGTGCAGCACAGGTCTCGAACGCTGTGCCGTACCTGTGGGAGCGGGCGTGCCCGCGAATGGCCGCAACCCTTTTATCTGGTTGCCGCCAGCCCCGCcaaca                                                                                                                                                                   <  1:683936/235‑1 (MQ=255)                                                                             ttACCGCGTTGAATCGCTGGCCCGGGTGAAGTCGCTGGAAGGTGGTGTGGTCAGCTTCATCGAATAAACGCGACACCTGCTTCGCGGGCACGCCCGCTCCCACAGGTGCAGCACAGGTCTCGAACGCTGTGCCGTACCTGTGGGAGCGGGCGTGCCCGCGAATGGCCGCAACCCTTTTATCTGGTTGCCGCCAGCCCCGCCAACAACAAGCGCTGGTAAAGATCTTCTTTAAGccc                                                                                                                                    >  1:76954/1‑236 (MQ=255)                                                                                                      gttaaGTCGCTGGAAGGTGGGGTGGTCAGCTTCATCGAATAAACGCGACACCTGCTTCGCGGGCACGCCCGCTCCCACAGGTGCAGCAAAGGTCTAGAACGCTGTGCCGTAACTGTGGGAGCGGGCGTGCCCGCGAATGGcc                                                                                                                                                                                                         <  2:662691/139‑1 (MQ=38)                                                                                                                               tCAGCTTCATCGAATAAACGCGACACCTGCTTCGCGGGCACGCCCGCTCCCACAGGTGCAGCACAGGTCTCGAACGCTGTGCCGTACCTGTGGGAGCGGGCGTGCCCGCGAATGGCCGCAACCCTTTTATCTGGTTGCCGCCAGCCCCGCCAACAACAAGCGCTGGTAAAGTTCTTCTTTAAGCCCCTGCGGCTGCCCCAACTGCATCCGCCCCAGGTGGCTTTCATACCCCTc                                                                                    >  2:691546/1‑234 (MQ=255)                                                                                                                                                   cgACACCTGCTTCGCGGGCACGCCCGCTCCCACAGGTGCAGCACAGGTCTCGAACGCTGTGCCGTACCTGTGGGAGCGGGCGTGCCCGCGAATGGCCGCAACCCTTTTATCTGGTTGCCGCCAGCCCCGCCAACAACAAGCGCTGGTAAAGTTCTTCTTTAAGCCCCTGCGGCTGCCCCAACTGCATCCGCCCCAGGTGGCTTTCATa                                                                                          >  1:249647/1‑208 (MQ=255)                                                                                                                                                      cacCTGCTTCGCGGGCCCGCCCGCTCCCACAGGTGCAGCACAGGTCTCGAACGCTGTGCCGTACCTGTGGGAGCGGGCGTGCCCGCGAATGGCCGCAACCCTTTTATCTGGTTGCCGCCAGCCCCGCCAACAACAAGCGCTGGTAAAGTTcttctt                                                                                                                                           <  2:22123/156‑1 (MQ=255)                                                                                                                                                                                 cacaGGTGCAGCACAGGTCTCGAACGCTGTGCCGTACCTGTGGGAGCGGGCGTGCCCGCGAATGGCCGCAACCCTTTTATCTGGTTGCCGCCAGCCCCGCCAACAACAAGCGCTGGTAAAGTTCTTCTTTAAGCCCCTGCGGCTGCCCCAACTGCATCCGCCCCAGGTGGCTTTCATACCCCTCCGGCCCCGGTGCATCCAACGCTGCCTTGCCCAGCTCCAGCACCTCGCTCAg                                 >  1:558278/1‑235 (MQ=255)                                                                                                                                                                                            cacaGGTCTCGAACGCTGTGCCGTCCCGGTGGGAGCGGGCGTGCCCGCGAATGGCCGCAACCCTTTTATCTGGTTGCCGCCAGCCCCGCCAACAACAAGCGCTGGTAAAGTTCTTCTTTAAGCCCCTGCGGCTGCCCCAACTGCATCCGCCCCAGGTGGCTTTc                                                                                             <  1:251739/164‑1 (MQ=255)                                                                                                                                                                                                 gTCTCGAACGCTGTGCCGTACCTGTGGGAGCGGGCGTGCCCGCGAATGGCCGCAACCCTTTTATCTGGTTGCCGCCAGCCCCGCCAACAACAAGCGCTGGTAAAGTTCTTCGTTAAGCCCCTGCGGCTGCCCCAACTGCATc                                                                                                              >  2:478885/1‑142 (MQ=255)                                                                                                                                                                                                        aCGCTGTGCCGTACCTGTGGGAGCGGGCGTGCCCGCGAATGGCCGCAACCCTTTTATCTGGTTGCCGCCAGCCCCGCCAACAACAAGCGCTGGTAAAGTTCTTCTTTAAGCCCCTGCGGCTGCCCCAACTGCATCCGCCCCAGGTGGCTTTCATACCCCTCCGGCCCCGGTGCATCCAACGCTGCCTTGCCCAGCTCCAGCACCTCGCTCAGCTTGAACTTGCTCTTCAACCAGt          <  1:691548/235‑1 (MQ=255)                                                                                                                                                                                                          ggtgtgCCGGACCTGGGGGAGCGGGCGTGCCCGCGAATGGCCGCAACCCGTTGATCTGGTTGCCGCCAGCCCCGCCAACAACAAGCGCTGGTAAAGTTCTTCTTTAAGc                                                                                                                                      <  2:40590/107‑1 (MQ=38)                                                                                                                                                                                                                       tgtgGGAGCGGGCGTGCCCGCGAATGGCCGCAACCCTTTTATCTGGTTGCCGCCAGCCCCGCCAACAACAAGCGCTGGTAAAGTTCTTCTTTAAGCCCCTGCGGCTGCCCCAACTGCATCCGCCCCAGGTGGCTTTCATACCCCTCCGGCCCCGGTGCATCCAACGCTGCCTTg                                                        >  1:652901/1‑174 (MQ=255)                                                                                                                                                                                                                       tgtgGGAGCGGGCGTGCCCGCGAATGGCCGCAACCCTTTTATCTGGTTGCCGCCAGCCCCGCCAACAACAAGCGCTGGTAAAGTTCTTCTTTAAGCCCCTGCGGCTGCCCCAACTGCATCCGCCCCAGGTGGCTTTCATACCCCTCCGGCCCCGGTGCATCCAACGCTGCCTTg                                                        <  2:652899/174‑1 (MQ=255)                                                                                                                                                                                                                                  gCGTGCCCGCGAATGGCCGCAACCCTTTTATCTGGTTGCCGCCAGCCCCGCCAACAACAAGCGCTGGTAAAGTTCTTCTTTAAGCCCCTGCGGCTGCCCCAACTGCATCCGCCCCAGGTGGCTTTCATACCCCTCCGGCCCCGGTGCATCCAACGCTGCCTTGCCCAGCTCCAGCACCTCGCTCAGCTTGCACTTGCTCTTCAACCAGTTCAGCGcc  >  1:61988/1‑217 (MQ=255)                                                                                                                                                                                                                                   cGTGCCCGCGAATGGCCGCAACCCTTTTATCTGGTTGCCGCCAGCCCCGCCAACAACAAGCGCTGGTAAAGTTCTTCTTTAAGCCCCTGCGGCTGCCCCAACTGCATCCGCCCCAGGTGGCTTTCATAccc                                                                                       <  1:731299/131‑1 (MQ=255)                                                                                                                                                                                                                                   cGTGCCCGCGAATGGCCGCAACCCTTTTATCTGGTTGCCGCCAGCCCCGCCAACAACAAGCGCTGGTAAAGTTCTTCTTTAAGCCCCTGCGGCTGCCCCAACTGCATCCGCCCCAGGTGGCTTTCATAccc                                                                                       >  2:731297/1‑131 (MQ=255)                                                                                                                                                                                                                                    gTGCCCGCGAATGGCCGCACCCCTTTTATCTGGTTGCCGCCAGCCCCGCCAACAACAAGCGCTGGTAAAGTTCTTCTTTAAGCCCCTGCGGCTGCCCCAACTGCATCCGCCCCAGGTGGCTTTCATACCCCTCCGGCCCCGGTGCATCCAACGCTGCCTTGCCCAGCTCCAGCACCTCGCTCAGCTTGAACTTGc                      <  1:209037/195‑1 (MQ=255)                                                                                                                                                                                                                                         |                                                                                                                                                                                                                  GTGCCACCGTGGCGGGCCTGGGTATCGTGATCGAGAAGTCGTTCCAGGGCGGCCGTGCCGAACTGGACAGCCAGGGTTACCGCGTTGAATCGCTGGCCCGGGTGAAGTCGCTGGAAGGTGGTGTGGTCAGCTTCATCGAATAAACGCGACACCTGCTTCGCGGGCACGCCCGCTCCCACAGGTGCAGCACAGGTCTCGAACGCTGTGCCGTACCTGTGGGAGCGGGCGTGCCTGCGAATGGCCGCAACCCTTTTATCTGGTTGCCGCCAGCCCCGCCAACAACAAGCGCTGGTAAAGTTCTTCTTTAAGCCCCTGCGGCTGCCCCAACTGCATCCGCCCCAGGTGGCTTTCATACCCCTCCGGCCCCGGTGCATCCAACGCTGCCTTGCCCAGCTCCAGCACCTCGCTCAGCTTGAACTTGCTCTTCAACCAGTTCAGCGCC  >  AE015451/6013653‑6014094

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