Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	NC_000913	1,447,666	T→C	S50P (TCT→CCT)	feaB →	phenylacetaldehyde dehydrogenase

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	1,447,666	0	T	C	100.0%	105.3 / NA	30	S50P (TCT→CCT)	feaB	phenylacetaldehyde dehydrogenase
Reads supporting (aligned to +/- strand):  ref base T (0/0);  new base C (12/18);  total (12/18)

GCATGTAGCAGTATTAAGCCAGGTCCAACAGTTTCTCGATCGTCAACACGGTCTTTATATTGATGGTCGTCCTGGCCCCGCACAAAGTGAAAAACGGTTGGCGATCTTTGATCCGGCCACCGGGCAAGAAATTGCGTCTACTGCTGATGCCAACGAAGCGGATGTAGATAACGCAGTCATGTCTGCCTGGCGGGCCTTTGTCTCGCGTCGCTGGGCCGGGCGATTACCCGCAGAGCGTGAACGTATTCTGCTACGTTTTGCTGATCTGGTGG  >  NC_000913/1447530‑1447801                                                                                                                                         |                                                                                                                                        gCATGTAGCAGTATTAAGCCAGGTCCAACAGTTTCTCGATCGTCAACACGGTCTTTATATTGATGGTCGTCCTGGCCCCGCACAAAGTGAAAAACGGTTGGCGATCTTTGATCCGGCCACCGGGCAAGAAATTGCGCCTa                                                                                                                                      <  2:91856/140‑1 (MQ=255) gCATGTAGCAGTATTAAGCCAGGTCCAACAGTTTCTCGATCGTCAACACGGTCTTTATATTGATGGTCGTCCTGGCCCCGCACAAAGTGAAAAACGGTTGGCGATCTTTGATCCGGCCACCGGGCAAGAAATTGCGCCTa                                                                                                                                      <  2:307733/140‑1 (MQ=255)            tATTAAGCCAGGTCCAACAGTTTCTCGATCGTCAACACGGTCTTTATATTGATGGTCGTCCTGGCCCCGCACAAAGTGAAAAACGGTTGGCGATCTTTGATCCGGCCACCGGGCAAGAAATTGCGCCTACTGCTGATGcc                                                                                                                           >  1:294682/1‑140 (MQ=255)                 aGCCAGGTCCAACAGTTTCTCGATCGTCAACACGGTCTTTATATTGATGGTCGTCCTGGCCCCGCACAAAGTGAAAAACGGTTGGCGGTCTTTGATCCGGCCACCGGGCAAGAAATTGCGCCTACTGCTGATGCCAACGa                                                                                                                      >  2:249926/1‑140 (MQ=255)                 aGCCAGGTCCAACAGTTTCTCGATCGTCAACACGGTCTTTATATTGATGGTCGTCCTGGCCCCGCACAAAGTGAAAAACGGTTGGCGATCTTTGATCCGGCCACCGGGCAAGAAATTGCGCCTACTGCTGATGCCAACGa                                                                                                                      >  2:292475/1‑140 (MQ=255)                   ccAGGTCCAACAGTTTCTCGATCGTCAACACGGTCTTTATATTGATGGTCGTCCTGGCCCCGCACAAAGTGAAAAACGGTTGGCGATCTTTGATCCGGCCACCGGGCAAGAAATTGCGCCTACTGCTGATGCCAACGAAg                                                                                                                    <  2:242865/140‑1 (MQ=255)                   ccAGGTCCAACAGTTTCTCGATCGTCAACACGGTCTTTATATTGATGGTCGTCCTGGCCCCGCACAAAGTGAAAAACGGTTGGCGATCTTTGATCCGGCCACCGGGCAAGAAATTGCGCCTACTGCTGATGCCAACGAAg                                                                                                                    >  2:305222/1‑140 (MQ=255)                      ggTCCAACAGTTTCTCGATCGTCAACACGGTCTTTATATTGATGGTCGTCCTGGCCCCGCACAAAGTGAAAAACGGTTGGCGATCTTTGATCCGGCCACCGGGCAAGAAATTGCGCCTACTGCTGATGCCAACGAAGCgg                                                                                                                 >  2:258064/1‑140 (MQ=255)                        tCCAACAGTTTCTCGATCGTCAACACGGTCTTTATATTGATGGTCGTCCTGGCCCCGCACAAAGTGAAAAACGGTTGGCGATCTTTGATCCGGCCACCGGGCAAGAAATTGCGCCTACTGCTGATGCCAACGAAGCGGAt                                                                                                               <  1:3098/140‑1 (MQ=255)                        tCCAACAGTTTCTCGATCGTCAACACGGTCTTTATATTGATGGTCGTCCTGGCCCCGCACAAAGTGAAAAACGGTTGGCGATCTTTGATCCGGCCACCGGGCAAGAAATTGCGCCTACTGCTGATGCCAACGAAGCGGAt                                                                                                               >  1:168881/1‑140 (MQ=255)                            aCAGTTTCTCGATCGTCAACACGGTCTTTATATTGATGGTCGTCCTGGCCCCGCACAAAGTGAAAAACGGTTGGCGATCTTTGATCCGGCCACCGGGCAAGAAATTGCGCCTACTGCTGATGCCAACGAAGCGGATGTAg                                                                                                           >  2:298524/1‑140 (MQ=255)                             cAGTTTCTCGATCGTCAACACGGTCTTTATATTGATGGTCGTCCTGGCCCCGCACAAAGTGAAAAACGGTTGGCGATCTTTGATCCGGCCACCGGGCAAGAAATTGCGCCTACTGCTGATGCCAACGAAGCGGATGTAGa                                                                                                          <  2:168881/140‑1 (MQ=255)                                      gATCGTCAACACGGTCTTTATATTGATGGTCGTCCTGGCCCCGCACAAAGTGAAAAACGGTTGGCGATCTTTGATCCGGCCACCGGGCAAGAAATTGCGCCTACTGCTGATGCCAACGAAGCGGATGTAGATAACGCAGt                                                                                                 <  2:362420/140‑1 (MQ=255)                                      gATCGTCAACACGGTCTTTATATTGATGGTCGTCCTGGCCCCGCACAAAGTGAAAAACGGTTGGCGATCTTTGATCCGGCCACCGGGCAAGAAATTGCGCCTACTGCTGATGCCAACGAAGCGGATGTAGATAACGCAGt                                                                                                 <  2:516920/140‑1 (MQ=255)                                                    tCTTTATATTGATGGTCGTCCTGGCCCCGCACAAAGTGAAAAACGGTTGGCGATCTCTGATCCGGCCACCGGGCAAGAAATTGCGCCTACTGCTGATGCCAACGAAGCGGATGTAGATAACGCAGTCATGTCTGCCTGGc                                                                                   >  1:402130/1‑140 (MQ=255)                                                     cTTTATATTGATGGTCGTCCTGGCCCCGCACAAAGTGAAAAACGGTTGGCGATCTTTGATCCGGCCACCGGGCAAGAAATTGCGCCTACTGCTGATGCCAACGAAGCGGATGTAGATAACGCAGTCATGTCTGCCTGGCg                                                                                  <  2:294682/140‑1 (MQ=255)                                                             tGATGGTCGTCCTGGCCCCGCACAAAGTGAAAAACGGTTGGCGATCTTTGATCCGGCCACCGGGCAAGAAATTGCGCCTACTGCTGATGCCAACGAAGCGGATGTAGATAACGCAGTCATGTCTGCCTGGCGGGCCTTTg                                                                          <  1:305222/140‑1 (MQ=255)                                                              gATGGTCGTCCTGGCCCCGCACAAAGTGAAAAACGGTTGGCGATCTTTGATCCGGCCACCGGGCAAGAAATTGCGCCTACTGCTGATGCCAACGAAGCGGATGTAGATAACGCAGTCATGTCTGCCTGGCGGGCCTTTGt                                                                         <  1:119110/140‑1 (MQ=255)                                                               aTGGTCGTCCTGGCCCCGCACAAAGTGAAAAACGGTTGGCGATCTTTGATCCGGCCACCGGGCAAGAAATTGCGCCTACTGCTGATGCCAACGAAGCGGATGTAGATAACGCAGTCATGTCTGCCTGGCGGGCCTTTGtc                                                                        >  1:112965/1‑140 (MQ=255)                                                               aTGGTCGTCCTGGCCCCGCACAAAGTGAAAAACGGTTGGCGATCTTTGATCCGGCCACCGGGCAAGAAATTGCGCCTACTGCTGATGCCAACGAAGCGGATGTAGATAACGCAGTCATGTCTGCCTGGCGGGCCTTTGtc                                                                        >  1:457562/1‑140 (MQ=255)                                                                  gtcgtcCTGGCCCCGCACAAAGTGAAAAACGGTTGGCGGTCTTTGATCCGGCCACCGGGCAAGAAATTGCGCCTACTGCTGATGCCAACGAAGCGGATGTAGATAACGCAGTCATGTCTGCCTGGCGGGCCTTTGTCTcg                                                                     <  1:249926/140‑1 (MQ=255)                                                                                           aaaaCGGTTGGCGATCTTTGATCCGGCCACCGGGCAAGAAATTGCGCCTACTGCTGATGCCAACGAAGCGGATGTAGATAACGCAGTCATGTCTGCCTGGCGGGCCTTTGTCTCGCGTCGCTGGGCCGGGCGATTACCCg                                            <  1:298524/140‑1 (MQ=255)                                                                                           aaaaCGGTTGGCGATCTTTGATCCGGCCACCGGGCAAGAAATTGCGCCTACTGCTGATGCCAACGAAGCGGATGTAGATAACGCAGTCATGTCTGCCTGGCGGGCCTTTGTCTCGCGTCGCTGGGCCGGGCGATTACCCg                                            <  1:292475/140‑1 (MQ=255)                                                                                                ggTTGGCGATCTTTGATCCGGCCACCGGGCAAGAAATTGCGCCTACTGCTGATGCCAACGAAGCGGATGTAGATAACGCAGTCATGTCTGCCTGGCGGGCCTTTGTCTCGCGTCGCTGGGCCGGGCGATTACCCGCagag                                       <  2:112965/140‑1 (MQ=255)                                                                                                ggTTGGCGATCTTTGATCCGGCCACCGGGCAAGAAATTGCGCCTACTGCTGATGCCAACGAAGCGGATGTAGATAACGCAGTCATGTCTGCCTGGCGGGCCTTTGTCTCGCGTCGCTGGGCCGGGCGATTACCCGCagag                                       <  1:258064/140‑1 (MQ=255)                                                                                                ggTTGGCGATCTTTGATCCGGCCACCGGGCAAGAAATTGCGCCTACTGCTGATGCCAACGAAGCGGATGTAGATAACGCAGTCATGTCTGCCTGGCGGGCCTTTGTCTCGCGTCGCTGGGCCGGGCGATTACCCGCagag                                       <  2:457562/140‑1 (MQ=255)                                                                                                        aTCTTTGATCCGGCCACCGGGCAAGAAATTGCGCCTACTGCTGATGCCAACGAAGCGGATGTAGATAACGCAGTCATGTCTGCCTGGCGGGCCTTTGTCTCGCGTCGCTGGGCCGGGCGATTACCCGCAGAGCGTGAACg                               <  1:383726/140‑1 (MQ=255)                                                                                                          cTTTGATCCGGCCACCGGGCAAGAAATTGCGCCTACTGCTGATGCCAACGAAGCGGATGTAGATAACGCAGTCATGTCTGCCTGGCGGGCCTTTGTCTCGCGTCGCTGGGCCGGGCGATTACCCGCAGAGCGTGAACGTa                             <  2:136048/140‑1 (MQ=255)                                                                                                                    gCCACCGGGCAAGAAATTGCGCCTACTGCTGATGCCAACGAAGCGGATGTAGATAACGCAGTCATGTCTGCCTGGCGGGCCTTTGTCTCGCGTCGCTGGGCCGGGCGATTACCCGCAGAGCGTGAACGTATTCTGCTACg                   >  1:264857/1‑140 (MQ=255)                                                                                                                                     tGCGCCTACTGCTGATGCCAACGAAGCGGATGTAGATAACGCAGTCATGTCTGCCTGGCGGGCCTTTGTCTCGCGTCGCTGGGCCGGGCGATTACCCGCAGAGCGTGAACGTATTCTGCTACGTTTTGCTGATCtggtgg  >  1:111305/1‑140 (MQ=255)                                                                                                                                         |                                                                                                                                        GCATGTAGCAGTATTAAGCCAGGTCCAACAGTTTCTCGATCGTCAACACGGTCTTTATATTGATGGTCGTCCTGGCCCCGCACAAAGTGAAAAACGGTTGGCGATCTTTGATCCGGCCACCGGGCAAGAAATTGCGTCTACTGCTGATGCCAACGAAGCGGATGTAGATAACGCAGTCATGTCTGCCTGGCGGGCCTTTGTCTCGCGTCGCTGGGCCGGGCGATTACCCGCAGAGCGTGAACGTATTCTGCTACGTTTTGCTGATCTGGTGG  >  NC_000913/1447530‑1447801

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