Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	AE015451	4,741,289	C→G	intergenic (+526/+227)	gltA → / ← yeiW	citrate synthase/putative oxidoreductase

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	AE015451	4,741,289	0	C	G	100.0%	49.6 / NA	45	intergenic (+526/+227)	gltA/yeiW	citrate synthase/putative oxidoreductase
Reads supporting (aligned to +/- strand):  ref base C (0/0);  new base G (29/16);  total (29/16)

AACTGGGTATAGGTCCAGCTCTGGTCACCTTCGTGCTTGAG‑AACCTCCGCCAACTGCAATTCGTGAA‑CCAAGCAAAAGGCCCCGAAAAGGCCCTCC‑GCAGCCCCCCCAACCTGGCCAGCGGCGCCTTGATCACGGTCACACCGCGATAGGCCTTGTACTTGGACCCGGCCACCTCACTGAGGGACACCTTGATCCC  >  AE015451/4741188‑4741383                                                                                                         |                                                                                               aaCTGGGTATAGGTCCAGCTCTGGTCACCTTCGTGCTTGAG‑CAGCTTCTGCGACTTGCATTCGTGAA‑TCCAGGCACAGGCCCCGGACACGTCCTCCTGCAGCGCCTCGACCTTGGCCAGCGgcgc                                                                          <  1:4022809/125‑1 (MQ=255)  aCTGGGTATAGGTCCAGCTCTGGTCACCTTCGTGCTTGAGCAGCTTCTGCGACTTGC‑ATTCGTGAA‑TCCAGGCACAGGCCCCGGACACGTCCTCCTGCAGCGCCTCGACCTTGGCCAGCGGCGcc                                                                         <  1:4022811/125‑1 (MQ=255)   cTGGGTATAGGTCCAGCTCTGGTCACCTTCGTGCTTGAGCAGCTTCTGCGACTTGC‑ATTCGTGAA‑TCCAGGCACAGGCCCCGGACACGTCCTCCTGCAGCGCCTCGACCTTGGCCAGCGGCGCCt                                                                        <  1:4022812/125‑1 (MQ=255)    tGGGTATAGGTCCAGCTCTGGTCACCTTCGTGCTTGAG‑CAGCTTCTGCGACTTGCATTCGTGAA‑TCCAGGCACAGGCCCCGGACACGTCCTCCTGCAGCGCCTCGACCTTGGCCAGCGGCGCCtt                                                                       >  2:4022947/1‑125 (MQ=255)     gggTATAGGTCCAGCTCTGGTCACCTTCGTGCTTGAG‑CAGCTTCTGCGACTTGCATTCGTGAA‑TCCAGGCACAGGCCCCGGACACGTCCTCCTGCAGCGCCTCGACCTTGGCCAGCGGCGcc                                                                         <  1:4022814/122‑1 (MQ=255)      ggTATAGGTCCAGCTCTGGTCACCTTCGTGCTTGAG‑CAGCTTCTGCGACTTGCATTCGTGAA‑TCCAGGCACAGGCCCCGGACACGTCCTCCTGCAGCGCCTCGACCTTGGCCAGCGGCGCCTTg                                                                      >  1:4022984/1‑124 (MQ=255)      ggTATAGGTCCAGCTCTGGTCACCTTCGTGCTTGAG‑CAGCTTCTGCGACTTGCATTCGTGAA‑TCCAGGCACAGGCCCCGGACACGTCCTCCTGCAGCGCCTCGACCTTGGCCAGCGGCGCCTTGa                                                                     >  1:4022954/1‑125 (MQ=255)      ggTATAGGTCCAGCTCTGGTCACCTTCGTGCTTGAG‑CAGCTTCTGCGACTTGCATTCGTGAA‑TCCAGGCACAGGCCCCGGACACGTCCTCCTGCAGCGCCTCGACCTTGGCCAGCGGCGCCTTGa                                                                     >  1:4022853/1‑125 (MQ=255)        tataGGTCCAGCTCTGGTCACCTTCGTGCTTGAGCAGCTTCTGCGACTTGC‑ATTCGTGAA‑TCCAGGCACAGGCCCCGGACACGTCCTCCTGCAGCGCCTCGACCTTGGCCAGCGGCGCCTTGATc                                                                   <  1:4022815/125‑1 (MQ=255)         ataGGTCCAGCTCTGGTCACCTTCGTGCTTGAGCAGCTTCTGCGACTTGC‑ATTCGTGAA‑TCCAGGCACAGGCCCCGGACACGTCCTCCTGCAGCGCCTCGACCTTGGCCAGCGGCGCCTTGATc                                                                   >  1:4022978/1‑124 (MQ=255)         ataGGTCCAGCTCTGGTCACCTTCGTGCTTGAGCAGCTTCTGCGACTTGC‑ATTCGTGAA‑TCCAGGCACAGGCCCCGGACACGTCCTCCTGCAGCGCCTCGACCTTGGCCAGCGGCGCCTTGATc                                                                   >  1:4022955/1‑124 (MQ=255)         ataGGTCCAGCTCTGGTCACCTTCGTGCTTGAGCAGCTTCTGCGACTTGC‑ATTCGTGAA‑TCCAGGCACAGGCCCCGGACACGTCCTCCTGCAGCGCCTCGACCTTGGCCAGCGGCGCCTTGATCa                                                                  <  1:4022816/125‑1 (MQ=255)         ataGGTCCAGCTCTGGTCACCTTCGTGCTTGAG‑CAGCTTCTGCGACTTGCATTCGTGAA‑TCCAGGCACAGGCCCCGGACACGTCCTCCTGCAGCGCCTCGACCTTGGCCAGCGGCGCCTTGATc                                                                   >  2:4023027/1‑124 (MQ=255)         ataGGTCCAGCTCTGGTCACCTTCGTGCTTGAG‑CAGCTTCTGCGACTTGCATTCGTGAA‑TCCAGGCACAGGCCCCGGACACATCCTCCTGCAGCGCCTCGACCTTGGCCAGCGGCGCCTTGATCa                                                                  <  1:4022817/125‑1 (MQ=255)           aGGTCCAGCTCTGGTCACCTTCGTGCTTGAGCAGCTTCTGCGACTTGC‑ATTCGTGAA‑TCCAGGCACAGGCCCCGGACACGTCCTCCTGCAGCGCCTCGACCTTGGCCAGCGGCGCCTTGATCACg                                                                >  1:4023097/1‑125 (MQ=255)            ggTCCAGCTCTGGTCACCTTCGTGCTTGAGCAGCTTCTGCGACTTGC‑ATTCGTGAA‑TCCAGGCACAGGCCCCGGACACGTCCTCCTGCAGCGCCTCGACCTTGGCCAGCGGCGCCTTGATCACgg                                                               >  2:4022925/1‑125 (MQ=255)              tCCAGCTCTGGTCACCTTCGTGCTTGAG‑CAGCTTCTGCGACTTGCATTCGTGAA‑TCCAGGCACAGGCCCCGGACACGTCCTCCTGCAGCGCCTCGACCTTGGCCAGCGGCGCCTTGATCACGGTc                                                             >  1:4023041/1‑125 (MQ=255)               ccAGCTCTGGTCACCTTCGTGCTTGAGCAGCTTCTGCGACTTGC‑ATTCGTGAA‑TCCAGGCACAGGCCCCGGACACGTCCTCCTGCAGCGCCTCGACCTTGGCCAGCGGCGCCTTGATCACGGTc                                                             >  2:4022888/1‑124 (MQ=255)                cAGCTCTGGTCACCTTCGTGCTTGAGCAGCTTCTGCGACTTGC‑ATTCGTGAATCCAGGC‑ACAGGCCCCGGACACGTCCTCCTGCAGCGCCTCGACCTTGGCCAGCGGCGCCTTGATCACGGTcac                                                           <  2:4022819/125‑1 (MQ=255)                cAGCTCTGGTCACCTTCGTGCTTGAGCAGCTTCTGCGACTTGC‑ATTCGTGAA‑TCCAGGCACAGGCCCCGGACACGTCCTCCTGCAGCGCCTCGACCTTGGCCAGCGGCGCCTTGATCACGGTcac                                                           >  1:4022914/1‑125 (MQ=255)                 aGCTCTGGTCACCTTCGTGCTTGAGCAGCTTCTGCGACTTGC‑ATTCGTGAATCCAGGC‑ACAGGCCCCGGACACGTCCTCCTGCAGCGCCTCGACCTTGGCCAGCGGCGCCTTGATCACGGTcaca                                                          >  2:4023020/1‑125 (MQ=255)                 aGCTCTGGTCACCTTCGTGCTTGAGCAGCTTCTGCGACTTGC‑ATTCGTGAA‑TCCAGGCACAGGCCCCGGACACGTCCTCCTGCAGCGCCTCGACCTTGGCCAGCGGCGCCTTGATCACGGTcac                                                           >  1:4022991/1‑124 (MQ=255)                  gCTCTGGTCACCTTCGTGCTTGAGCAGCTTCTGCGACTTGC‑ATTCGTGAA‑TCCAGGCACAGGCCCCGGACACGTCCTCCTGCAGCGCCTCGACCTTGGCCAGCGGCGCCTTGATCACGGTCacac                                                         >  2:4023105/1‑125 (MQ=255)                  gCTCTGGTCACCTTCGTGCTTGAGCAGCTTCTGCGACTTGC‑ATTCGTGAA‑TCCAGGCACAGGCCCCGGACACGTCCTCCTGCAGCGCCTCGACCTTGGCCAGCGGCGCCTTGATCACGGTCacac                                                         >  1:4023070/1‑125 (MQ=255)                  gCTCTGGTCACCTTCGTGCTTGAGCAGCTTCTGCGACTTGC‑ATTCGTGAA‑TCCAGGCACAGGCCCCGGACACGTCCTCCTGCAGCGCCTCGACCTTGGCCAGCGGCGCCTTGATCACGGTCacac                                                         >  2:4022957/1‑125 (MQ=255)                  gCTCTGGTCACCTTCGTGCTTGAGCAGCTTCTGCGACTTGC‑ATTCGTGAA‑TCCAGGCACAGGCCCCGGACACGTCCTCCTGCAGCGCCTCGACCTTGGCCAGCGGCGCCTTGATCACGGTCacac                                                         >  1:4023257/1‑125 (MQ=255)                    tctGGTCACCTTCGTGCTTGAGCAGCTTCTGCGACTTGC‑ATTCGTGAATCCAGGC‑ACAGGCCCCGGACACGTCCTCCTGCAGCGCCTCGACCTTGGCCAGCGGCGCCTTGATCACGGTCACACcg                                                       <  2:4022825/125‑1 (MQ=255)                      tGGTCACCTTCGTGCTTGAGCAGCTTCTGCGACTTGC‑ATTCGTGAA‑TCCAGGCACAGGCCCCGGACACGTCCTCCTGCAGCGCCTCGACCTTGGCCAGCGGCGCCTTGATCACGGTCACACcgcg                                                     <  1:4022834/125‑1 (MQ=255)                        gTCACCTTCGTGCTTGAG‑CAGCTTCTGCGACTTGCATTCGTGAA‑TCCAGGCACAGGCCCCGGATACGTCCTCCTGCAGCGCCTCGACCTTGGCCAGCGGCGCCTTGATCACGGTCACACCGCGAt                                                   <  1:4022850/125‑1 (MQ=255)                          cACCTTCGTGCTTGAGCAGCTTCTGCGACTTGC‑ATTCGTGAA‑TCCAGGCACAGGCCCCGGACACGTCCTCCTGCAGCGCCTCGACCTTGGCCAGCGGCGCCTTGATCACGGTCACACCGCGATAg                                                 <  2:4022851/125‑1 (MQ=255)                           aCCTTCGTGCTTGAGCAGCTTCTGCGACTTGC‑ATTCGTGAA‑TCCAGGCACAGGCCCCGGACACGTCCTCCTGCAGCGCCTCGACCTTGGCCAGCGGCGCCTTGATCACGGTCACACCGCGATAgg                                                >  2:4023061/1‑125 (MQ=255)                            ccTTCGTGCTTGAGCAGCTTCTGCGACTTGC‑ATTCGTGAATCCAGGCAC‑AGGCCCCGGACACGTCCTCCTGCAGCGCCTCGACCTTGGCCAGCGGCGCCTTGATCACGGTCACACCGCGATAGGc                                               <  1:4022857/125‑1 (MQ=255)                             cTTCGTGCTTGAGCAGCTTCTGCGACTTGC‑ATTCGTGAA‑TCCAGGCACAGGCCCCGGACACGTCCTCCTGCAGCGCCTCGACCTTGGCCAGCGGCGCCTTGATCACGGTCACACCGCGATAGGcc                                              >  1:4023015/1‑125 (MQ=255)                             cTTCGTGCTTGAGCAGCTTCTGCGACTTGC‑ATTCGTGAA‑TCCAGGCACAGGCCCCGGACACGTCCTCCTGCAGCGCCTCGACCTTGGCCAGCGGCGCCTTGATCACGGTCACACCGCGATAGGcc                                              <  1:4022832/125‑1 (MQ=255)                                                      acttgcaTTCGTGAA‑TCCAGGCACAGGCCCCGGACACGTCCTCCTGCAGCGCCTCGACCTTGGCCAGCGGCGCCTTGATCACGGTCACACCGCGATAGGCCTTGTACTTGGACCCGGCCAcc                         >  2:4023076/7‑122 (MQ=255)                                                      acttgcaTTCGTGAA‑TCCAGGCACAGGCCCCGGACACGTCCTCCTGCAGCGCCTCGACCTTGGCCAGCGGCGCCTTGATCACGGTCACACCGCGATAGGCCTTGTACTTGGACCCGGCCACCTc                       >  1:4023040/7‑124 (MQ=255)                                                      acttgcaTTCGTGAA‑TCCAGGCACAGGCCCCGGACACGTCCTCCTGCAGCGCCTCGACCTTGGCCAGCGGCGCCTTGATCACGGTCACACCGCGATAGGCCTTGTACTTGGACCCGGCCACCTc                       >  1:4023075/7‑124 (MQ=255)                                                      acttgcaTTCGTGAA‑TCCAGGCACAGGCCCCGGACACGTCCTCCTGCAGCGCCTCGACCTTGGCCAGCGGCGCCTTGATCACGGTCACACCGCGATAGGCCTTGTACTTGGACCCGGCCACCTc                       >  2:4022972/7‑124 (MQ=255)                                                      acttgcaTTCGTGAA‑TCCAGGCACAGGCCCCGGACACGTCCTCCTGCAGCGCCTCGACCTTGGCCAGCGGCGCCTTGATCACGGTCACACCGCGATAGGCCTTGTACTTAGACCCGGCCACCTc                       >  2:4023275/7‑124 (MQ=255)                                                       cttgcaTTCGTGAA‑TCCAGGCACAGGCCCCGGACACGTCCTCCTGCAGCGCCTCGACCTTGGCCAGCGGCGCCTTGATCACGGTCACACCGCGATAGGCCTTGTACTTGGACCCGGCCACCTCAc                     >  1:4023053/6‑125 (MQ=255)                                                         tgcaTTCGTGAA‑TCCAGGCACAGGCCCCGGACACGTCCTCCTGCAGCGCCTCGACCTTGGCCAGCGGCGCCTTGATCACGGTCACACCGCGATAGGCCTTGTACTTGGACCCGGCCACCTCACt                    <  1:4022835/121‑1 (MQ=255)                                                           caTTCGTGAATCCAGGC‑ACAGGCCCCGGACACGTCCTCCTGCAGCGCCTCGACCTTGGCCAGCGGCGCCTTGATCACGGTCACACCGCGATAGGCCTTGTACTTGGACCCGGCCACCTCACTGAg                 >  1:4023132/2‑125 (MQ=255)                                                                         aggcacaGGCCCCGGACACGTCCTCCTGCAGCGCCTCGACCTTGGCCAGCGGCGCCTTGATCACGGTCACACCGCGATAGGCCTTGTACTTGGACCCGGCCACCTCACTGAGGGACACCTTGATc    >  1:4023085/7‑125 (MQ=255)                                                                           gcacaGGCCCCGGACACGTCCTCCTGCAGCGCCTCGACCTTGGCCAGCGGCGCCTTGATCACGGTCACACCGCGATAGGCCTTGTACTTGGACCCGGCCACCTCACTGAGGGACACCTTGATccc  <  2:4022865/121‑1 (MQ=255)                                                                             acaGGCCCCGGACACGTCCTCCTGCAGCGCCTCGACCTTGGCCAGCGGCGCCTTGATCACGGTCACACCGCGATAGGCCTTGTACTTGGACCCGGCCACCTCAc                     >  2:4023159/3‑104 (MQ=255)                                                                                                         |                                                                                               AACTGGGTATAGGTCCAGCTCTGGTCACCTTCGTGCTTGAG‑AACCTCCGCCAACTGCAATTCGTGAA‑CCAAGCAAAAGGCCCCGAAAAGGCCCTCC‑GCAGCCCCCCCAACCTGGCCAGCGGCGCCTTGATCACGGTCACACCGCGATAGGCCTTGTACTTGGACCCGGCCACCTCACTGAGGGACACCTTGATCCC  >  AE015451/4741188‑4741383

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