Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	NC_000913	4,151,639	T→C	H487R (CAC→CGC)	ppc ←	phosphoenolpyruvate carboxylase

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	4,151,639	0	T	C	100.0%	80.7 / NA	24	H487R (CAC→CGC)	ppc	phosphoenolpyruvate carboxylase
Reads supporting (aligned to +/- strand):  ref base T (0/0);  new base C (15/9);  total (15/9)

CGATACCAGTCAATATTGAGCAGCTGGGTCATGACATCGTTGGCGTTGTTCAGATCATCGAGGGTTTCAAACAGCGGAGCAACCGGCATCGCAAACCCGATACCCGCTTCTTTCAGCAGCAGGTGGACAGCCAGTACGTCGGACGGCGTTTTCGCCATCGAGATCACGTAGGCGGCAATGGAGCCTTGCGGTGCTTCGGCAATCACCTGGCAGGTATCGAGCACTTCGCGCGTTT  >  NC_000913/4151516‑4151750                                                                                                                            |                                                                                                                cGATACCAGTCAATATTGAGCAGCTGGGTCATGACATCGTTGGCGTTGTTCAGATCATCGAGGGTTTCAAACAGCGGAGCAACCGGCATCGCAAACCCGATACCCGCTTCTTTCAGCAGCAGGCGGACAGCCa                                                                                                        >  2:330793/1‑133 (MQ=255)  gATACCAGTCAATATTGAGCAGCTGGGTCATGACATCGTTGGCGTTGTTCAGATCATCGAGGGTTTCAAACAGCGGAGCAACCGGCATCGCAAACCCGATACCCGCTTCTTTCAGCAGCAGGCGGACAGCCAg                                                                                                       >  1:395432/1‑133 (MQ=255)       cAGTCAATATTGAGCAGCTGGGTCATGACATCGTTGGCGTTGTTCAGATCATCGAGGGTTTCAAACAGCGGAGCAACCGGCATCGCAAACCCGATACCCGCTTCTTTCAGCAGCAGGCGGACAGCCAGTACGt                                                                                                  >  2:404752/1‑133 (MQ=255)                 tGAGCAGCTGGGTCATGACATCGTTGGCGTTGTTCAGATCATCGAGGGTTTCAAACAGCGGAGCAACCGGCATCGCAAACCCGATACCCGCTTCTTTCAGCAGCAGGCGGACAGCCAGTACGTCGGACGGCGt                                                                                        >  1:410655/1‑133 (MQ=255)                 tGAGCAGCTGGGTCATGACATCGTTGGCGTTGTTCAGATCATCGAGGGTTTCAAACAGCGGAGCAACCGGCATCGCAAACCCGATACCCGCATCTTTCAGGAGCGGGCTGACAGCCAATACGGCGGACGGCGt                                                                                        >  2:178900/1‑133 (MQ=255)                    gcagcTGGGTCATGACATCGTTGGCGTTGTTCAGATCATCGAGGGTTTCAAACAGCGGAGCAACCGGCATCGCAAACCCGATACCCGCTTCTTTCAGCAGCAGGCGGACAGCCAGTACGTCGGACGGCGtttt                                                                                     <  1:271250/133‑1 (MQ=255)                        cTGGGTCATGACATCGTTGGCGTTGTTCAGATCATCGAGGGTTTCAAACAGCGGAGCAACCGGCATCGCAAACCCGATACCCGCTTCTTTCAGCAGCAGGCGGACAGCCAGTACGTCGGACGGCGTTTTCGcc                                                                                 >  1:466284/1‑133 (MQ=255)                          gggTCATGACATCGTTGGCGTTGTTCAGATCATCGAGGGTTTCAAACAGCGGAGCAACCGGCATCGCAAACCCGATACCCGCTTCTTTCAGCAGCAGGCGGACAGCCAGTACGTCGGACGGCGTTTTCGCCAt                                                                               <  2:466284/133‑1 (MQ=255)                           ggTCATGACATCGTTGGCGTTGTTCAGATCATCGAGGGTTTCAAACAGCGGAGCAACCGGCATCGCAAACCCGATACCCGCTTCTTTCAGCAGCAGGCGGACAGCCAGTACGTCGGACGGCGTTTTCGCCATc                                                                              >  1:114547/1‑133 (MQ=255)                           ggTCATGACATCGTTGGCGTTGTTCAGATCATCGAGGGTTTCAAACAGCGGAGCAACCGGCATCGCAAACCCGATACCCGCTTCTTTCAGCAGCAGGCGGACAGCCAGTACGTCGGACGGCGTTTTCGCCATc                                                                              <  1:315721/133‑1 (MQ=255)                                     tCGTTGGCGTTGTTCAGATCATCGAGGGTTTCAAACAGCGGAGCAACCGGCATCGCAAACCCGATACCCGCTTCTTTCAGCAGCAGGCGGACAGCCAGTACGTCGGACGGCGTTTTCGCCATCGAGATCACGt                                                                    >  2:430472/1‑133 (MQ=255)                                       gTTGGCGTTGTTCAGATCATCGAGGGTTTCAAACAGCGGAGCAACCGGCATCGCAAACCCGATACCCGCTTCTTTCAGCAGCAGGCGGACAGCCAGTACGTCGGACGGCGTTTTCGCCATCGAGATCACGTAg                                                                  >  2:76637/1‑133 (MQ=255)                                         tGGCGTTGTTAAGATTATCGAGGGTTTCAAACCGCGGAGCACCCGGCATCGCAAACCCGATACCCGCTTCTTTCAGCAGCAGGCGGACAGCCAGTACGTCGGACGGCGTTTTCGCCATCGAGATCACGTAggc                                                                <  2:169760/133‑1 (MQ=255)                                                gttCAGATCATCGAGGGTTTCAAACAGCGGAGCAACCGGCATCGCAAACCCGATTCCCGCTTCTTTCAGCAGCAGGCGGACAGCCAGTACGTCGGACGGCGTTTTCGCCATCGAGATCACGTAGGCGGCAATg                                                         <  2:395432/133‑1 (MQ=255)                                                 ttCAGATCATCGAGGGTTTCAAACAGCGGAGCAACCGGCATCGCAAACCCGATACCCGCTTCTTTCAGCAGCAGGCGGACAGCCAGTACGTTGGACGGCGGTTTAGCCATCGAGATCACGTACGCGGCAATgg                                                        >  2:189561/1‑133 (MQ=255)                                                     gATCATCGCGGGTTTCAATCAGTGGAGCAACCGGCATCGCAAACCCGATACCCGCTTCTTTCAGCAGCAGGCGGACAGCCAGTACGTCGGACGGCGTTTTCGCCATCGAGATCACGTAGGCGGCAATGGAGcc                                                    <  1:178900/133‑1 (MQ=255)                                                          tcGAGGGTTTTAAACAGCGGAGCAGCCGTCATCGCAAACCCGATACCCGCTTCTTTCAGCAGCAGGCGGACAGCCAGTACGTCGGACGGCGTTTTCGCCATCGAGATCACGTAGGCGGCAATGGAGCCTTGCg                                               <  2:206756/133‑1 (MQ=255)                                                           cGAGGGTTTCAAACAGCGGAGCAACCGGCATCGCAAACCCGATACCCGCTTCTTTCAGCAGCAGGCGGACAGCCAGTACGTCGGACGGCGTTTTCGCCATCGAGATCACGTAGGCGGCAATGGAGCCTTGCgg                                              >  1:411244/1‑133 (MQ=255)                                                               ggTTTCAAACAGCGGAGCAACCGGCATCGCAAACCCGATACCCGCTTCTTTCAGCAGCAGGCGGACAGCCAGTACGTCGGACGGCGTTTTCGCCATCGAGATCACGTAGGCGGCAATGGAGCCt                                                   <  1:64614/124‑1 (MQ=255)                                                               ggTTTCAAACAGCGGAGCAACCGGCATCGCAAACCCGATACCCGCTTCTTTCAGCAGCAGGCGGACAGCCAGTACGTCGGACGGCGTTTTCGCCATCGAGATCACGTAGGCGGCAATGGAGCCt                                                   >  2:64614/1‑124 (MQ=255)                                                                   tCAAACAGCGGAGCAACCGGCATCGCAAACCCGATACCCGCTTCTTTCAGCAGCAGGCGGACAGCCAGTACGTCGGACGGCGTTTTCGCCATCGAGATCAAGTAAGCGACAATGGAGACTTGCGGTGCTTCgg                                      >  2:220681/1‑133 (MQ=255)                                                                    cAAACAGCGGAGCAACCGGCATCGCAAACCCGATACCCGCTTCTTTCAGCAGCAGGCGGACAGCCAGTACGTCGGACGGCGTTTTCGCCATCGAGATCACGTAGGCGGCAATGGAGCCTTGCGGTGCTTCGGc                                     >  2:78123/1‑133 (MQ=255)                                                                                 aaCCGGCATCGCAAACCCGATACCCGCTTCTTTCAGCAGCAGGCGGACAGCCAGTACGTCGGACGGCGTTTTCGCCATCGAGATCACGTAGGCGGCAATGGAGCCTTGCGGTGCTTCGGCAATCACCTGGCAg                        <  1:78123/133‑1 (MQ=255)                                                                                                       cccGCTTCTTTCAGCAGCAGGCGGACAGCCAGTACGTCGGACGGCGTTTTCGCCATCGAGATCACGTAGGCGGCAATGGAGCCTTGCGGTGCTTCGGCAATCACCTGGCAGGTATCGAGCACTTCGCGCGttt  >  1:65128/1‑133 (MQ=255)                                                                                                                            |                                                                                                                CGATACCAGTCAATATTGAGCAGCTGGGTCATGACATCGTTGGCGTTGTTCAGATCATCGAGGGTTTCAAACAGCGGAGCAACCGGCATCGCAAACCCGATACCCGCTTCTTTCAGCAGCAGGTGGACAGCCAGTACGTCGGACGGCGTTTTCGCCATCGAGATCACGTAGGCGGCAATGGAGCCTTGCGGTGCTTCGGCAATCACCTGGCAGGTATCGAGCACTTCGCGCGTTT  >  NC_000913/4151516‑4151750

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