Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	NC_000913	4,252,781	T→C	R92R (CGT→CGC)	ubiC →	chorismate lyase

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	4,252,781	0	T	C	100.0%	121.3 / NA	34	R92R (CGT→CGC)	ubiC	chorismate lyase
Reads supporting (aligned to +/- strand):  ref base T (0/0);  new base C (15/19);  total (15/19)

AGCGTGACGATGATCCGCGAAGGGTTTGTCGAGCAGAATGAAATCCCCGAAGAACTGCCGCTGCTGCCGAAAGAGTCTCGTTACTGGTTACGTGAAATTTTGTTATGTGCCGATGGTGAACCGTGGCTTGCCGGTCGTACCGTCGTTCCTGTGTCAACGTTAAGCGGGCCGGAGCTGGCGTTACAAAAATTGGGTAAAACGCCGTTAGGACGCTATCTGTTCACATCATCGACATTAACCCGGGACTTTATTGAGATAGGCCG  >  NC_000913/4252644‑4252906                                                                                                                                          |                                                                                                                              aGCGTGACGATGATCCGCGAAGGGTTTGTCGAGCAGAATGAAATCCCCGAAGAACTGCCGCTGCTGCCGAAAGAGTCTCGTTACTGGTTACGTGAAATTTTGTTATGTGCCGATGGTGAACCGTGGCTTGCCGGTCGCAc                                                                                                                             >  2:456954/1‑140 (MQ=255)      gACGATGATCCGCGAAGGGTTTGTCGAGCAGAATGAAATCCCCGAAGAACTGCCGCTGCTGCCGAAAGAGTCTCGTTACTGGTTACGTGAAATTTTGTTATGTGCCGATGGTGAACCGTGGCTTGCCGGTCGCACcgtcg                                                                                                                        >  1:145447/1‑140 (MQ=255)         gatgatCCGCGAAGGGTTTGTCGAGCAGAATGAAATCCCCGAAGAACTGCCGCTGCTGCCGAAAGAGTCTCGTTACTGGTTACGTGAAATTTTGTTATGTGCCGATGGTGAACCGTGGCTTGCCGGTCGCACCGTCGTTc                                                                                                                     <  1:399174/140‑1 (MQ=255)          atgatCCGCGAAGGGTTTGTCGAGCAGAATGAAATCCCCGAAGAACTGCCGCTGCTGCCGAAAGAGTCTCGTTACTGGTTACGTGAAATTTTGTTATGTGCCGATGGTGAACCGTGGCTTGCCGGTCGCACCGTCGTTcc                                                                                                                    <  1:543101/140‑1 (MQ=255)            gatCCGCGAAGGGTTTGTCGAGCAGAATGAAATCCCCGAAGAACTGCCGCTGCTGCCGAAAGAGTCTCGTTACTGGTTACGTGAAATTTTGTTATGTGCCGATGGTGAACCGTGGCTTGCCGGTCGCACCGTCGTTCCtg                                                                                                                  >  2:324154/1‑140 (MQ=255)             atCCGCGAAGGGTTTGTCGAGCAGAATGAAATCCCCGAAGAACTGCCGCTGCTGCCGAAAGAGTCTCGTTACTGGTTACGTGAAATTTTGTTATGTGCCGATGGTGAACCGTGGCTTGCCGGTCGCACCGTCGTTcc                                                                                                                    >  2:303873/1‑137 (MQ=255)             atCCGCGAAGGGTTTGTCGAGCAGAATGAAATCCCCGAAGAACTGCCGCTGCTGCCGAAAGAGTCTCGTTACTGGTTACGTGAAATTTTGTTATGTGCCGATGGTGAACCGTGGCTTGCCGGTCGCACCGTCGTTcc                                                                                                                    <  1:303873/137‑1 (MQ=255)                  cgAAGGGTTTGTCGAGCAGAATGAAATCCCCGAAGAACTGCCGCTGCTGCCGAAAGAGTCTCGTTACTGGTTACGTGAAATTTTGTTATGTGCCGATGGTGAACCGTGGCTTGCCGGTCGCACCGTCGTTCCTGTGTCaa                                                                                                            >  1:114967/1‑140 (MQ=255)                         tttGTCGAGCAGAATGAAATCCCCGAAGAACTGCCGCTGCTGCCGAAAGAGTCTCGTTACTGGTTACGTGAAATTTTGTTATGTGCCGATGGTGAACCGTGGCTTGCCGGTCGCACCGTCGTTcc                                                                                                                    <  1:165835/125‑1 (MQ=255)                         tttGTCGAGCAGAATGAAATCCCCGAAGAACTGCCGCTGCTGCCGAAAGAGTCTCGTTACTGGTTACGTGAAATTTTGTTATGTGCCGATGGTGAACCGTGGCTTGCCGGTCGCACCGTCGTTcc                                                                                                                    >  2:165835/1‑125 (MQ=255)                         tttGTCGAGCAGAATGAAATCCCCGAAGAACTGCCGCTGCTGCCGAAAGAGTCTCGTTACTGGTTACGTGAAATTTTGTTATGTGCCGATGGTGAACCGTGGCTTGCCGGTCGCACCGTCGTTCCTGTGTCAACGTTAAg                                                                                                     <  1:324154/140‑1 (MQ=255)                          ttGTCGAGCAGAATGAAATCCCCGAAGAACTGCCGCTGCTGCCGAAAGAGTCTCGTTACTGGTTACGTGAAATTTTGTTATGTGCCGATGGTGAACCGTGGCTTGCCGGTCGCACCGTCGTTCCTGTGTCAACGTTAAGc                                                                                                    >  2:497476/1‑140 (MQ=255)                            gTCGAGCAGAGTGAAATCCCCGAAGAATTGCCTCTGCTGCCGAAGGAGTTTCGTTACTGTTGACGTGAAATTTTGTTTTGTGCCGATGGTGAACCGTGGCTTGCCGGCCACACCTGAGTTCCTGGGTTAACGTTaatctg                                                                                                  >  2:133918/1‑136 (MQ=255)                                      aTGAAATCCCCGAAGAACTGCCGCTGCTGCCGAAAGAGTCTCGTTACTGGTTACGTGAAATTTTGTTATGTGCCGATGGTGAACCGTGGCTTGCCGGTCGCACCGTCGTTCCTGTGTCAACGTTAAGCGGGCCGGAGCTg                                                                                        <  2:145447/140‑1 (MQ=255)                                                           cGCTGCTGCCGAAAGAGTCTCGTTACTGGTTACGTGAAATTTTGTTATGTGCCGATGGTGAACCGTGGCTTGCCGGTCGCACCGTCGTTCCTGTGTCAACGTTAAGCGGGCCGGAGCTGGCGTTACAAAAATTGGGTaaa                                                                   <  2:114967/140‑1 (MQ=255)                                                            gctgctGCCGAAAGAGTCTCGTTACTGGTTACGTGAAATTTTGTTATGTGCCGATGGTGAACCGTGGCTTGCCGGTCGCACCGTCGTTCCTGTGTCAACGTTAAGCGGGCCGGAGCTGGCGTTACAAAAATTGGGTaaaa                                                                  <  2:227051/140‑1 (MQ=255)                                                                        agagTCTCGTTACTGGTTACGTGAAATTTTGTTATGTGCCGATGGTGAACCGTGGCTTGCCGGTCGCACCGTCGTTCCTGTGTCAACGTTAAGCGGGCCGGAGCTGGCGTTACAAAAATTGGGTAAAACGCCGTTAGGAc                                                      <  1:133918/140‑1 (MQ=255)                                                                          agTCTCGTTACTGGTTACGTGAAATTTTGTTATGTGCCGATGGTGAACCGTGGCTTGCCGGTCGCACCGTCGTTCCTGTGTCAACGTTAAGCGGGCCGGAGCTGGCGTTACAAAAATTGGGTAAAACGCCGTTAGGACGc                                                    <  2:383694/140‑1 (MQ=255)                                                                                 ttACTGGTTACGTGAAATTTTGTTATGTGCCGATGGTGAACCGTGGCTTGCCGGTCGCACCGTCGTTCCTGTGTCAACGTTAAGCGGGCCGGAGCTGGCGTTACAAAAATTGGGTAAAACGCCGTTAGGACGCTATCTGt                                             >  1:159168/1‑140 (MQ=255)                                                                                   aCTGGTTACGTGAAATTTTGTTATGTGCCGATGGTGAACCGTGGCTTGCCGGTCGCACCGTCGTTCCTGTGTCAACGTTAAGCGGGCCGGAGCTGGCGTTACAAAAATTGGGTAAAACGCCGTTAGGACGCTATCTGTTc                                           >  2:397702/1‑140 (MQ=255)                                                                                     tGGTTACGTGAAATTTTGTTATGTGCCGATGGTGAACCGTGGCTTGCCGGTCGCACCGTCGTTCCTGTGTCAACGTTAAGCGGGCCGGAGCTGGCGTTACAAAAATTGGGTAAAACGCCGTTAGGACGCTATCTGTTcac                                         <  1:497476/140‑1 (MQ=255)                                                                                      ggTTACGTGAAATTTTGTTATGTGCCGATGGTGAACCGTGGCTTGCCGGTCGCACCGTCGTTCCTGTGTCAACGTTAAGCGGGCCGGAGCTGGCGTTACAAAAATTGGGTAAAACGCCGTTAGGACGCTATCTGTTcaca                                        <  1:125436/140‑1 (MQ=255)                                                                                          aCGTGAAATTTTGTTATGTGCCGATGGTGAACCGTGGCTTGCCGGTCGCACCGTCGTTCCTGTGTCAACGTTAAGCGGGCCGGAGCTGGCGTTACAAAAATTGGGTAAAACGCCGTTAGGACGCTATCTGTTCAcatcat                                    <  1:397702/140‑1 (MQ=255)                                                                                            gTGAAATTTTGTTATGTGCCGATGGTGAACCGTGGCTTGCCGGTCGCACCGTCGTTCCTGTGTCAACGTTAAGCGGGCCGGAGCTGGCGTTACAAAAATTGGGTAAAACGCCGTTAGGACGCTATCTGTTCACATCATCg                                  >  1:465117/1‑140 (MQ=255)                                                                                            gTGAAATTTTGTTATGTGCCGATGGTGAACCGTGGCTTGCCGGTCGCACCGTCGTTCCTGTGTCAACGTTAAGCGGGCCGGAGCTGGCGTTACAAAAATTGGGTAAAACGCCGTTAGGACGCTATCTGTTCACATCATCg                                  >  2:344272/1‑140 (MQ=255)                                                                                                aaTTTTGTTATGTGCCGATGGTGAACCGTGGCTTGCCGGTCGCACCGTCGTTCCTGTGTCAACGTTAAGCGGGCCGGAGCTGGCGTTACAAAAATTGGGTAAAACGCCGTTAGGACGCTATCTGTTCACATCATCGACAt                              <  2:137006/140‑1 (MQ=255)                                                                                                aaTTTTGTTATGTGCCGATGGTGAACCGTGGCTTGCCGGTCGCACCGTCGTTCCTGTGTCAACGTTAAGCGGGCCGGAGCTGGCGTTACAAAAATTGGGTAAAACGCCGTTAGGACGCTATCTGTTCACATCATCGACAt                              <  2:465117/140‑1 (MQ=255)                                                                                                 aTTTTGTTATGTGCCGATGGTGAACCGTGGCTTGCCGGTCGCACCGTCGTTCCTGTGTCAACGTTAAGCGGGCCGGAGCTGGCGTTACAAAAATTGGGTAAAACGCCGTTAGGACGCTATCTGTTCACATCATCGACAtt                             >  2:358526/1‑140 (MQ=255)                                                                                                     tGTTATGTGCCGATGGTGAACCGTGGCTTGCCGGTCGCACCGTCGTTCCTGTGTCAACGTTAAGCGGGCCGGAGCTGGCGTTACAAAAATTGGGTAAAACGCCGTTAGGACGCTATCTGTTCACATCATCGACATTAAcc                         <  2:431044/140‑1 (MQ=255)                                                                                                         aTGTGCCGATGGTGAACCGTGGCTTGCCGGTCGCACCGTCGTTCCTGTGTCAACGTTAAGCGGGCCGGAGCTGGCGTTACAAAAATTGGGTAAAACGCCGTTAGGACGCTATCTGTTCACATCATCGACATTAACCCggg                     >  1:239420/1‑140 (MQ=255)                                                                                                          tgtgCCGATGGTGAACCGTGGCTTGCCGGTCGCACCGTCGTTCCTGTGTCAACGTTAAGCGGGCCGGAGCTGGCGTTACAAAAATTGGGTAAAACGCCGTTAGGACGCTATCTGTTCACATCATCGACATTAACCCGGGa                    <  1:355347/140‑1 (MQ=255)                                                                                                           gtgCCGATGGTGAACCGTGGCTTGCCGGTCGCACCGTCGTTCCTGTGTCAACGTTAAGCGGGCCGGAGCTGGCGTTACAAAAATTGGGTAAAACGCCGTTAGGACGCTATCTGTTCACATCATCGACATTAACCCGGGAc                   <  1:344272/140‑1 (MQ=255)                                                                                                                            tGGCTTGCCGGTCGCACCGTCGTTCCTGTGTCAACGTTAAGCGGGCCGGAGCTGGCGTTACAAAAATTGGGTAAAACGCCGTTAGGACGCTATCTGTTCACATCATCGACATTAACCCGGGACTTTATTGAGATAGGCCg  >  2:312680/1‑140 (MQ=255)                                                                                                                              gCTTGCCGGTCGCACCGTCGTTCCTGTGTCAACGTTAAGCGGGCCGGAGCTGGCGTTACAAAAATTGGGTAAAACGCCGTTAGGACGCTATCTGTTCACATCATCGACATTAACCCGGGACTTTATTGAGATAGGCCg  <  1:137271/138‑1 (MQ=255)                                                                                                                              gCTTGCCGGTCGCACCGTCGTTCCTGTGTCAACGTTAAGCGGGCCGGAGCTGGCGTTACAAAAATTGGGTAAAACGCCGTTAGGACGCTATCTGTTCACATCATCGACATTAACCCGGGACTTTATTGAGATAGGCCg  >  2:137271/1‑138 (MQ=255)                                                                                                                                          |                                                                                                                              AGCGTGACGATGATCCGCGAAGGGTTTGTCGAGCAGAATGAAATCCCCGAAGAACTGCCGCTGCTGCCGAAAGAGTCTCGTTACTGGTTACGTGAAATTTTGTTATGTGCCGATGGTGAACCGTGGCTTGCCGGTCGTACCGTCGTTCCTGTGTCAACGTTAAGCGGGCCGGAGCTGGCGTTACAAAAATTGGGTAAAACGCCGTTAGGACGCTATCTGTTCACATCATCGACATTAACCCGGGACTTTATTGAGATAGGCCG  >  NC_000913/4252644‑4252906

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: ‑