Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	NC_000913	3,810,304	G→A	*160* (TAG→TAA)	coaD →	pantetheine‑phosphate adenylyltransferase

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	3,810,304	0	G	A	100.0%	80.7 / NA	26	*160* (TAG→TAA)	coaD	pantetheine‑phosphate adenylyltransferase
Reads supporting (aligned to +/- strand):  ref base G (0/0);  new base A (12/14);  total (12/14)

TGTGTTTCTGATGCCGTCGAAAGAGTGGTCGTTTATCTCTTCATCGTTGGTGAAAGAGGTGGCGCGCCATCAGGGCGATGTCACCCATTTCCTGCCGGAGAATGTCCATCAGGCGCTGATGGCGAAGTTAGCGTAGCGTTTATGCCGGATGGTATGCCATCCGGCGCGCATGAATTACTTCTGGCACTGCCGACAATAAAACGTTGCCCGCTGCGCATGTTTAGTCGCCACAATCGGCGTACCGCACACCCGACACGGCTCACCT  >  NC_000913/3810169‑3810433                                                                                                                                        |                                                                                                                                  tgtgtTTCTGATGCCGTCGAAAGAGTGGTCGTTTATCTCTTCATCGTTGGTGAAAGAGGTGGCGCGCCATCAGGGCGATGTCACCCATTTCCTGCCGGAGAATGTCCATCAGGCGCTGATGGCGAAGTTagcgtaacgt                                                                                                                                >  2:513017/1‑139 (MQ=255)        cTGATGCCGTCGAAAGAGTGGTCGTTTATCTCTTCATCGTTGGTGAAAGAGGTGGCGCGCCATCAGGGCGATGTCACCCATTTCCTGCCGGAGAATGTCCATCAGGCGCTGATGGCGAAGTTAGCGTAACGTTTATGcc                                                                                                                         <  2:146009/139‑1 (MQ=255)        cTGATGCCGTCGAAAGAGTGGTCGTTTATCTCTTCATCGTTGGTGAAAGAGGTGGCGCGCCATCAGGGCGATGTCACCCATTTCCTGCCGGAGAATGTCCATCAGGCGCTGATGGCGAAGTTAGCGTAACGTTTATGcc                                                                                                                         <  2:1026480/139‑1 (MQ=255)                        agTGGTCGTTTATCTCTTCATCGTTGGTGAAAGAGGTGGCGCGCCATCAGGGCGATGTCACCCATTTCCTGCCGGAGAATGTCCATCAGGCGCTGATGGCGAAGTTAGCGTAACGTTTATGCCGGATGGTATGCCATcc                                                                                                         >  1:1066446/1‑139 (MQ=255)                               gTTTATCTCTTCATCGTTGGTGAAAGAGGTGGCGCGCCATCAGGGCGATGTCACCCATTTCCTGCCGGAGAATGTCCATCAGGCGCTGATGGCGAAGTTAGCGTAACGTTTATGCCGGATGGTATGCCATCCGgcgcgc                                                                                                  <  2:1066446/139‑1 (MQ=255)                                      tctTCATCGTTGGTGAAAGAGGTGGCGCGCCATCAGGGCGATGTCACCCATTTCCTGCCGGAGAATGTCCATCAGGCGCTGATGGCGAAGTTAGCGTAACGTTTATGCCGGATGGTATGCCATCCGGCGCGCATGAAtt                                                                                           <  2:688630/139‑1 (MQ=255)                                        ttCATCGTTGGTGAAAGAGGTGGCGCGCCATCAGGGCGATGTCACCCATTTCCTGCCGGAGAATGTCCATCAGGCGCTGATGGCGAAGTTAGCGTAACGTTTATGCCg                                                                                                                        >  1:345018/1‑108 (MQ=255)                                        ttCATCGTTGGTGAAAGAGGTGGCGCGCCATCAGGGCGATGTCACCCATTTCCTGCCGGAGAATGTCCATCAGGCGCTGATGGCGAAGTTAGCGTAACGTTTATGCCg                                                                                                                        <  2:345018/108‑1 (MQ=255)                                           aTCGTTGGTGAAAGAGGTGGCGCGCCATCAGGGCGATGTCACCCATTTCCTGCCGGAGAATGTCCATCAGGCGCTGATGGCGAAGTTAGCGTAACGTTTATGCCGGATGGTATGCCATCCGGCGCGCa                                                                                                 <  1:205691/128‑1 (MQ=255)                                           aTCGTTGGTGAAAGAGGTGGCGCGCCATCAGGGCGATGTCACCCATTTCCTGCCGGAGAATGTCCATCAGGCGCTGATGGCGAAGTTAGCGTAACGTTTATGCCGGATGGTATGCCATCCGGCGCGCa                                                                                                 >  2:205691/1‑128 (MQ=255)                                               ttGGTGAAAGAGGTGGCGCGCCATCAGGGCGATGTCACCCATTTCCTGCCGGAGAATGTCCATCAGGCGCTGATGGCGAAGTTAGCGTAACGTTTATGCCGGATGGTATGCCATCCGGCGCGCATGAATTACTTCTGGc                                                                                  >  2:13376/1‑139 (MQ=255)                                                     aaaGAGGTGGCGCGCCATCAGGGCGATGTCACCCATTTCCTGCCGGAGAATGTCCATCAGGCGCTGATGGCGAAGTTAGCGTAACGTTTATGCCGGATGGTATGCCATCCGGCGCGCATGAATTACTTCTGGCACTGcc                                                                            <  1:513017/139‑1 (MQ=255)                                                            tGGCGCGCCATCAGGGCGATGTCACCCATTTCCTGCCGGAGAATGTCCATCAGGCGCTGATGGCGAAGTTAGCGTAACGTTTATGCCGGATGGTATGCCATCCGGCGCGCATGAATTACTTCTGGCACTGCCGACAATa                                                                     >  1:252765/1‑139 (MQ=255)                                                            tGGCGCGCCATCAGGGCGATGTCACCCATTTCCTGCCGGAGAATGTCCATCAGGCGCTGATGGCGAAGTTAGCGTAACGTTTATGCCGGATGGTATGCCATCCGGCGCGCATGAATTACTTCTGGCACTGCCGACAATa                                                                     >  1:737324/1‑139 (MQ=255)                                                                                              gCCGGAGAATGTCCATCAGGCGCTGATGGCGAAGTTAGCGTAACGTTTATGCCg                                                                                                                        <  2:1006659/54‑1 (MQ=255)                                                                                              gCCGGAGAATGTCCATCAGGCGCTGATGGCGAAGTTAGCGTAACGTTTATGCCg                                                                                                                        >  1:1006659/1‑54 (MQ=255)                                                                                                     aaTGTCCATCAGGCGCTGATGGCGAAGTGAGCGTAACGTTTATGCCGGATGGTATGCCATCCGGCGCGCATGAATTACTTCTGGCACTGCCGACAATAAAACGTTGCCCGCTGCGCATGTTTAGTCGCCACAATCGGCg                            <  1:643502/139‑1 (MQ=255)                                                                                                         tCCATCAGGCGCTGATGGCGAAGTTAGCGTAACGTTTATGCCGGATg                                                                                                                    <  2:304526/47‑1 (MQ=255)                                                                                                         tCCATCAGGCGCTGATGGCGAAGTTAGCGTAACGTTTATGCCGGATg                                                                                                                    >  1:304526/1‑47 (MQ=255)                                                                                                         tCCATCAGGCGCTGATGGCGAAGTTAGCGTAACGTTTATGCCGGATGGTATGCCATCCGGCGCGCa                                                                                                 <  1:481599/66‑1 (MQ=255)                                                                                                         tCCATCAGGCGCTGATGGCGAAGTTAGCGTAACGTTTATGCCGGATGGTATGCCATCCGGCGCGCa                                                                                                 >  2:481599/1‑66 (MQ=255)                                                                                                                 gcgcTGATGGCGAAGTTAGCGTAACGTTTATGCCGGa                                                                                                                      >  1:919838/1‑37 (MQ=255)                                                                                                                 gcgcTGATGGCGAAGTTAGCGTAACGTTTATGCCGGa                                                                                                                      <  2:919838/37‑1 (MQ=255)                                                                                                                           cGAAGTTAGCGTAACGTTTATGCCGGATGGTATGCCATCCGGCGCGCa                                                                                                 <  1:734416/48‑1 (MQ=255)                                                                                                                           cGAAGTTAGCGTAACGTTTATGCCGGATGGTATGCCATCCGGCGCGCa                                                                                                 >  2:734416/1‑48 (MQ=255)                                                                                                                           cGAAGTTAGCGTAACGTTTATGCCGGATGGTATGCCATCCGGCGCGCATGAATTACTTCTGGCACTGCCGACAATAAAACGTTGCCCGCTGCGCATGTTTAGTCGCCACAATCGGCGTACCGCACACCCGACACGGCTc      <  2:737324/139‑1 (MQ=255)                                                                                                                               gTTAGCGTAACGTTTATGCCGGATGGTATGCCATCCGGCGCGCATGAATTACTTCTGGCACTGCCGACAATAAAACGTTGCCCGCTGCGCATGTTTAGTCGCCACAATCGGCGTACCGCACACCCGACACGGCTCACCt  >  1:465146/1‑139 (MQ=255)                                                                                                                                        |                                                                                                                                  TGTGTTTCTGATGCCGTCGAAAGAGTGGTCGTTTATCTCTTCATCGTTGGTGAAAGAGGTGGCGCGCCATCAGGGCGATGTCACCCATTTCCTGCCGGAGAATGTCCATCAGGCGCTGATGGCGAAGTTAGCGTAGCGTTTATGCCGGATGGTATGCCATCCGGCGCGCATGAATTACTTCTGGCACTGCCGACAATAAAACGTTGCCCGCTGCGCATGTTTAGTCGCCACAATCGGCGTACCGCACACCCGACACGGCTCACCT  >  NC_000913/3810169‑3810433

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