Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	AE015451	4,741,283	+T	intergenic (+520/+233)	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,283	1	.	T	100.0%	121.3 / NA	29	intergenic (+520/+233)	gltA/yeiW	citrate synthase/putative oxidoreductase
Reads supporting (aligned to +/- strand):  ref base . (0/0);  new base T (15/14);  total (15/14)

CACCGGCCATGGCGTATTGAACTGGGTATAGGTCCAGCTCTGGTCACCTTCGTGCTTGAG‑AA‑CCTCC‑GCCAAC‑TGCAATTCGTGAA‑CCAAGCAAAAGGCCCCGAAAAGGCCCTCC‑GCAGCCCCCCCAACCTGGCCAGCGGCGCCTTGATCACGGTCACACCGCGATAGGCCTTGTACTTGGACCCGGCCACCTCACTGAGGGACACCTTGATCCCGTCCTCGTCCTTGGCCACTTGCCAGTT  >  AE015451/4741169‑4741410                                                                                                                         |                                                                                                                                cACCGGCCATGGCGTATTGAACTGGGTATAGGTCCAGCTCTGGTCACCTTCGTGCTTGAGCAG‑CTTCT‑GCGACT‑TGC‑ATTCGTGAA‑TCCAGGCACAGGCCCCGGACACGTCCTCCTGCAGCGCCTCGACCTTGGCCAGCGGCGCCTTGATc                                                                                              <  1:4032716/151‑1 (MQ=255) cACCGGCCATGGCGTATTGAACTGGGTATAGGTCCAGCTCTGGTCACCTTCGTGCTTGAG‑CA‑GCTTC‑TGCGAC‑TTGCATTCGTGAA‑TCCAGGCACAGGCCCCGGACACGTCCTCCTGCAGCGCCTCGACCTTGGCCAGCGGCGCCTTGATc                                                                                              <  2:4032715/151‑1 (MQ=255)  aCCGGCCATGGCGTATTGAACTGGGTATAGGTCCAGCTCTGGTCACCTTCGTGCTTGAG‑CA‑GCTTC‑TGCGAC‑TTGCATTCGTGAATCCAGGC‑ACAGGCCCCGGACACGTCCTCCTGCAGCGCCTCGACCTTGGCCAGCGGCGCCTTGATc                                                                                              >  2:4032966/1‑150 (MQ=255)   ccGGCCATGGCGTATTGAACTGGGTATAGGTCCAGCTCTGGTCACCTTCGTGCTTGAG‑CAGCTTCT‑GCGACT‑TGC‑ATTCGTGAA‑TCCAGGCACAGGCCCCGGACACGTCCTCCTGCAGCGCCTCGACCTTGGCCAGCGGCGCCTTGATCAc                                                                                            <  1:4032718/151‑1 (MQ=255)     ggCCATGGCGTATTGAACTGGGTATAGGTCCAGCTCTGGTCACCTTCGTGCTTGAGCAG‑CTTCT‑GCGACT‑TGC‑ATTCGTGAATCCAGGC‑ACAGGCCCCGGACACGTCCTCCTGCAGCGCCTCGACCTTGGCCAGCGGCGCCTTGATCACgg                                                                                          >  2:4032970/1‑151 (MQ=255)     ggCCATGGCGTATTGAACTGGGTATAGGTCCAGCTCTGGTCACCTTCGTGCTTGAGCAG‑CTTCT‑GCGACT‑TGC‑ATTCGTGAATCCAGGC‑ACAGGCCCCGGACACGTCCTCCTGCAGCGCCTCGACCTTGGCCAGCGGCGCCTTGATCACgg                                                                                          >  2:4033061/1‑151 (MQ=255)       ccATGGCGTATTGAACTGGGTATAGGTCCAGCTCTGGTCACCTTCGTGCTTGAG‑CA‑GCTTCTGCGACT‑TGC‑ATTCGTGAATCCAGGC‑ACAGGCCCCGGACACGTCCTCCTGCAGCGCCTCGACCTTGGCCAGCGGCGCCTTGATCACGGTc                                                                                        >  2:4032927/1‑151 (MQ=255)        cATGGCGTATTGAACTGGGTATAGGTCCAGCTCTGGTCACCTTCGTGCTTGAG‑CAGCTTCT‑GCGACT‑TGC‑ATTCGTGAA‑TCCAGGCACAGGCCCCGGACACGTCCTCCTGCAGCGCCTCGACCTTGGCCAGCGGCGCCTTGATCACGGTca                                                                                       <  2:4032721/151‑1 (MQ=255)        cATGGCGTATTGAACTGGGTATAGGTCCAGCTCTGGTCACCTTCGTGCTTGAG‑CA‑GCTTCTGCGACT‑TGC‑ATTCGTGAA‑TCCAGGCACAGGCCCCGGACACGTCCTCCTGCAGCGCCTCGACCTTGGCCAGCGGCGCCTTGATCACGGTca                                                                                       <  1:4032722/151‑1 (MQ=255)        cATGGCGTATTGAACTGGGTATAGGTCCAGCTCTGGTCACCTTCGTGCTTGAG‑CA‑GCTTCTGCGACT‑TGC‑ATTCGTGAA‑TCCAGGCACAGGCCCCGGACACGTCCTCCTGCAGCGCCTCGACCTTGGCCAGCGGCGCCTTGATCACGGTca                                                                                       <  1:4032720/151‑1 (MQ=255)        cATGGCGTATTGAACTGGGTATAGGTCCAGCTCTGGTCACCTTCGTGCTTGAG‑CA‑GCTTC‑TGCGAC‑TTGCATTCGTGAATCCAGGC‑ACAGGCCCCGGACACGTCCTCCTGCAGCGCCTCGACCTTGGCCAGCGGCGCCTTGATCACGGTca                                                                                       >  1:4033048/1‑151 (MQ=255)        cATGGCGTATTGAACTGGGTATAGGTCCAGCTCTGGTCACCTTCGTGCTTGAG‑CA‑GCTTC‑TGCGAC‑TTGCATTCGTGAA‑TCCAGGCACAGGCCCCGGACACGTCCTCCTGCAGCGCCTCGACCTTGGCCAGCGGCGCCTTGATCACGGTca                                                                                       >  1:4033001/1‑151 (MQ=255)          tGGCGTATTGAACTGGGTATAGGTCCAGCTCTGGTCACCTTCGTGCTTGAG‑CA‑GCTTC‑TGCGACTTGC‑ATTCGTGAA‑TCCAGGCACAGGCCCCGGACACGTCCTCCTGCAGCGCCTCGACCTTGGCCAGCGGCGCCTTGATCACGGTcaca                                                                                     <  2:4032723/151‑1 (MQ=255)                          ggataGGTCCAGCTCTGGTCACCTTCGTGCTTGAG‑CA‑GCTTC‑TGCGAC‑TTGCATTCGTGAA‑TCCAGGCACAGGCCCCGGACACGTCCTCCTGCAGCGCCTCGACCTTGGCCAgc                                                                                                          >  2:4032929/3‑114 (MQ=255)                                 tCCAGCTCTGGTCACCTTCGTGCTTGAG‑CA‑GCTTC‑TGCGAC‑TTGCATTCGTGAATCCAGGC‑ACAGGCCCCGGACACGTCCTCCTGCAGCGCCTCGACCTTGGCCAGCGGCGCCTTGATCACGGTCACACCGCGATAGGCCTTGTACTTGGa                                                              >  2:4033046/1‑151 (MQ=255)                                  ccAGCTCTGGTCACCTTCGTGCTTGAGCAG‑CTTCT‑GCGACT‑TGC‑ATTCGTGAA‑TCCAGGCACAGGCCCCGGACACGTCCTCCTGCAGCGCCTCGACCTTGGCCAGCGGCGCCTTGATCACGGTCACACCGCGATAGGCCTTGTACTTGGAc                                                             <  1:4032754/151‑1 (MQ=255)                                  ccAGCTCTGGTCACCTTCGTGCTTGAG‑CA‑GCTTCTGCGACT‑TGC‑ATTCGTGAA‑TCCAGGCACAGGCCCCGGACACGTCCTCCTGCAGCGCCTCGACCTTGGCCAGCGGCGCCTTGATCACGGTCACACCGCGATAGGCCTTGTACTTGGa                                                              <  2:4032741/150‑1 (MQ=255)                                  ccAGCTCTGGTCACCTTCGTGCTTGAG‑CA‑GCTTC‑TGCGAC‑TTGCATTCGTGAATCCAGGC‑ACAGGCCCCGGACACGTCCTCCTGCAGCGCCTCGACCTTGGCCAGCGGCGCCTTGATCACGGTCACACCGCGATAGGCCTTGTACTTGGAc                                                             >  2:4033055/1‑151 (MQ=255)                                   cAGCTCTGGTCACCTTCGTGCTTGAG‑CA‑GCTTCTGCGACT‑TGC‑ATTCGTGAA‑TCCAGGCACAGGCCCCGGACACGTCCTCCTGCAGCGCCTCGACCTTGGCCAGCGGCGCCTTGATCACGGTCACACCGCGATAGGCCTTGTACTTGGAcc                                                            <  2:4032760/151‑1 (MQ=255)                                   cAGCTCTGGTCACCTTCGTGCTTGAG‑CA‑GCTTC‑TGCGAC‑TTGCATTCGTGAATCCAGGC‑ACAGGCCCCGGACACGTCCTCCTGCAGCGCCTCGACCTTGGCCAGCGGCGCCTTGATCACGGTCACACCGCGATAGGCCTTGTACTTGGAcc                                                            >  1:4032834/1‑151 (MQ=255)                                    aGCTCTGGTCACCTTCGTGCTTGAG‑CA‑GCTTCTGCGACT‑TGC‑ATTCGTGAA‑TCCAGGCACAGGCCCCGGACACGTCCTCCTGCAGCGCCTCGACCTTGGCCAGCGGCGCCTTGATCACGGTCACACCGCGATAGGCCTTGTACTTGGAccc                                                           >  1:4033111/1‑151 (MQ=255)                                    aGCTCTGGTCACCTTCGTGCTTGAG‑CA‑GCTTC‑TGCGAC‑TTGCATTCGTGAA‑TCCAGGCACAGGCCCCGGACACGTCCTCCTGCAGCGCCTCGACCTTGGCCAGCGGCGCCTTGATCACGGTCACACCGCGATAGGCCTTGTACTTGGAccc                                                           <  1:4032744/151‑1 (MQ=255)                                     gCTCTGGTCACCTTCGTGCTTGAGCAG‑CTTCT‑GCGACT‑TGC‑ATTCGTGAATCCAGGC‑ACAGGCCCCGGACACGTCCTCCTGCAGCGCCTCGACCTTGGCCAGCGGCGCCTTGATCACGGTCACACCGCGATAGGCCTTGTACTTGGACCCg                                                          <  2:4032782/151‑1 (MQ=255)                                     gCTCTGGTCACCTTCGTGCTTGAGCAG‑CTTCT‑GCGACT‑TGC‑ATTCGTGAATCCAGGC‑ACAGGCCCCGGACACGTCCTCCTGCAGCGCCTCGACCTTGGCCAGCGGCGCCTTGATCACGGTCACACCGCGATAGGCCTTGTACTTGGACCCg                                                          <  2:4032771/151‑1 (MQ=255)                                       tctGGTCACCTTCGTGCTTGAGCAG‑CTTCT‑GCGACT‑TGC‑ATTCGTGAATCCAGGC‑ACAGGCCCCGGACACGTCCTCCTGCAGCGCCTCGACCTTGGCCAGCGGCGCCTTGATCACGGTCACACCGCGATAGGCCTTGTACTTGGACCCGGc                                                        >  2:4032863/1‑151 (MQ=255)                                         tGGTCACCTTCGTGCTTGAG‑CA‑GCTTCTGCGACT‑TGC‑ATTCGTGAATCCAGGC‑ACAGGCCCCGGACACGTCCTCCTGCAGCGCCTCGACCTTGGCCAGCGGCGCCTTGATCACGGTCACACCGCGATAGGCCTTGTACTTGGACCCGGcc                                                       >  2:4033010/1‑150 (MQ=255)                                           gTCACCTTCGTGCTTGAG‑CA‑GCTTC‑TGCGAC‑TTGCATTCGTGAATCCAGGC‑ACAGGCCCCGGACACGTCCTCCTGCAGCGCCTCGACCTTGGCCAGCGGCGCCTTGATCACGGTCACACCGCGATAGGCCTTGTACTTGGACCCGGCCAcc                                                    >  2:4033132/1‑151 (MQ=255)                                           gTCACCTTCGTGCTTGAG‑CA‑GCTTC‑TGCGAC‑TTGCATTCGTGAA‑TCCAGGCACAGGCCCCGGACACGTCCTCCTGCAGCGCCTCGACCTTGGCCAGCGGCGCCTTGATCACGGTCACGCCGCGATAGGCCTTGTACTTGGACCCg                                                          >  2:4033254/1‑145 (MQ=255)                                                                                                  cacaGGCCCCGGACACGTCCTCCTGCAGCGCCTCGACCTTGGCCAGCGGCGCCTTGATCACGGTCACACCGCGATAGGCCTTGTACTTGGACCCGGCCACCTCACTGAGGGACACCTTGATCCCGTCCTCGTCCTTGGCCACTTGCCAGtt  <  2:4032795/148‑1 (MQ=255)                                                                                                                         |                                                                                                                                CACCGGCCATGGCGTATTGAACTGGGTATAGGTCCAGCTCTGGTCACCTTCGTGCTTGAG‑AA‑CCTCC‑GCCAAC‑TGCAATTCGTGAA‑CCAAGCAAAAGGCCCCGAAAAGGCCCTCC‑GCAGCCCCCCCAACCTGGCCAGCGGCGCCTTGATCACGGTCACACCGCGATAGGCCTTGTACTTGGACCCGGCCACCTCACTGAGGGACACCTTGATCCCGTCCTCGTCCTTGGCCACTTGCCAGTT  >  AE015451/4741169‑4741410

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