Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	NC_000913	4,151,642	A→G	V486A (GTC→GCC)	ppc ←	phosphoenolpyruvate carboxylase

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	4,151,642	0	A	G	100.0%	141.3 / NA	38	V486A (GTC→GCC)	ppc	phosphoenolpyruvate carboxylase
Reads supporting (aligned to +/- strand):  ref base A (0/0);  new base G (22/16);  total (22/16)

AATCAGGCCACGATACCAGTCAATATTGAGCAGCTGGGTCATGACATCGTTGGCGTTGTTCAGATCATCGAGGGTTTCAAACAGCGGAGCAACCGGCATCGCAAACCCGATACCCGCTTCTTTCAGCAGCAGGTGGACAGCCAGTACGTCGGACGGCGTTTTCGCCATCGAGATCACGTAGGCGGCAATGGAGCCTTGCGGTGCTTCGGCAATCACCTGGCAGGTATCGAGCACTTCGCGCGTTTCGGCGCTTGGTTGCCAGTTGCGC  >  NC_000913/4151506‑4151773                                                                                                                                         |                                                                                                                                    aaTCAGGCCACGATACCAGTCAATATTGAGCAGCTGGGTCATGACATCGTTGGCGTTGTTCAGATCATCGAGGGTTTCAAACAGCGGAGCAACCGGCATCGCAAACCCGATACCCGCTTCTTTCAGCAGCAGGTGGGCa                                                                                                                                   >  1:499478/1‑139 (MQ=255) aaTCAGGCCACGATACCAGTCAATATTGAGCAGCTGGGTCATGACATCGTTGGCGTTGTTCAGATCATCGAGGGTTTCAAACAGCGGAGCAACCGGCATCGCAAACCCGATACCCGCTTCTTTCAGCAGCAGGTGGGCa                                                                                                                                   >  2:804986/1‑139 (MQ=255)       gCCACGATACCAGTCAATATTGAGCAGCTGGGTCATGACATCGTTGGCGTTGTTCAGATCATCGAGGGTTTCAAACAGCGGAGCAACCGGCATCGCAAACCCGATACCCGCTTCTTTCAGCAGCAGGTGGGCAGCCAGt                                                                                                                             >  2:416721/1‑139 (MQ=255)           cGATACCAGTCAATATTGAGCAGCTGGGTCATGACATCGTTGGCGTTGTTCAGATCATCGAGGGTTTCAAACAGCGGAGCAACCGGCATCGCAAACCCGATACCCGCTTCTTTCAGCAGCAGGTGGGCAGCCAGt                                                                                                                             <  2:1023081/135‑1 (MQ=255)           cGATACCAGTCAATATTGAGCAGCTGGGTCATGACATCGTTGGCGTTGTTCAGATCATCGAGGGTTTCAAACAGCGGAGCAACCGGCATCGCAAACCCGATACCCGCTTCTTTCAGCAGCAGGTGGGCAGCCAGt                                                                                                                             >  1:1023081/1‑135 (MQ=255)                     cAATATTGAGCAGCTGGGTCATGACATCGTTGGCGTTGTTCAGATCATCGAGGGTTTCAAACAGCGGAGCAACCGGCATCGCAAACCCGATACCCGCTTCTTTCAGCAGCAGGTGGGCAGCCAGTACGTCGGACGGCGt                                                                                                               >  2:36597/1‑139 (MQ=255)                      aaTATTGAGCAGCTGGGTCATGACATCGTTGGCGTTGTTCAGATCATCGAGGGTTTCAAACAGCGGAGCAACCGGCATCGCAAACCCGATACCCGCTTCTTTCAGCAGCAGGTGGGCAGCCAGTACGTCGGACGGCGtt                                                                                                              >  2:320878/1‑139 (MQ=255)                        tatTGAGCAGCTGGGTCATGACATCGTTGGCGTTGTTCAGATCATCGAGGGTTTCAAACAGCGGAGCAACCGGCATCGCAAACCCGATACCCGCTTCTTTCAGCAGCAGGTGGGCAGCCAGTACGTCGGACGGCGtttt                                                                                                            <  1:320878/139‑1 (MQ=255)                          ttGAGCAGCTGGGTCATGACATCGTTGGCGTTGTTCAGATCATCGAGGGTTTCAAACAGCGGAGCAACCGGCATCGCAAACCCGATACCCGCTTCTTTCAGCAGCAGGTGGGCAGCCAGTACGTCGGACGGCGTTTTCg                                                                                                          <  2:499478/139‑1 (MQ=255)                             agcagcTGGGTCATGACATCGTTGGCGTTGTTCAGATCATCGAGGGTTTCAAACAGCGGAGCAACCGGCATCGCAAACCCGATACCCGCTTCTTTCAGCAGCAGGTGGGCAGCCAGTACGTCGGACGGCGTTTTCGCCa                                                                                                       >  2:1033791/1‑139 (MQ=255)                             agcagcTGGGTCATGACATCGTTGGCGTTGTTCAGATCATCGAGGGTTTCAAACAGCGGAGCAACCGGCATCGCAAACCCGATACCCGCTTCTTTCAGCAGCAGGTGGGCAGCCAGTACGTCGGACGGCGTTTTCGCCa                                                                                                       <  1:1084659/139‑1 (MQ=255)                             agcagcTGGGTCATGACATCGTTGGCGTTGTTCAGATCATCGAGGGTTTCAAACAGCGGAGCAACCGGCATCGCAAACCCGATACCCGCTTCTTTCAGCAGCAGGTGGGCAGCCAGTACGTCGGACGGCGTTTTCGCCa                                                                                                       >  2:1084659/1‑139 (MQ=255)                              gcagcTGGGTCATGACATCGTTGGCGTTGTTCAGATCATCGAGGGTTTCAAACAGCGGAGCAACCGGCATCGCAAACCCGATACCCGCTTCTTTCAGCAGCAGGTGGGCAGCCAGTACGTCGGACGGCGTTTTCGCCAt                                                                                                      >  2:124391/1‑139 (MQ=255)                              gcagcTGGGTCATGACATCGTTGGCGTTGTTCAGATCATCGAGGGGTTCAAACAGCGGAGCAACCGGCATCGCAAACCCGATACCCGCTTCTTTCAGCAGCAGGTGGGCAGCCAGTACGTCGGACGGCGTTTTCGCCAt                                                                                                      >  2:609292/1‑139 (MQ=255)                                       tCATGACATCGTTGGCGTTGTTCAGATCATCGAGGGTTTCAAACAGCGGAGCAACCGGCATCGCAAACCCGATACCCGCTTCTTTCAGCAGCAGGTGGGCAGCCAGTACGTCGGACGGCGTTTTCGCCATCGAGATCAc                                                                                             >  1:618853/1‑139 (MQ=255)                                            aCATCGTTGGCGTTGTTCAGATCATCGAGGGTTTCAAACAGCGGAGCAACCGGCATCGCAAACCCGATACCCGCTTCTTTCAGCAGCAGGTGGGCAGCCAGTACGTCGGACGGCGTTTTCGCCATCGAGATCACGTAgg                                                                                        >  1:633304/1‑139 (MQ=255)                                                 gTTGGCGTTGTTCAGATCATCGAGGGTTTCAAACAGCGGAGCAACCGGCATCGCAAACCCGATACCCGCTTCTTTCAGCAGCAGGTGGGCAGCCAGTACGTCGGACGGCGTTTTCGCCATCGAGATCACGTAGGCGGCa                                                                                   >  2:483502/1‑139 (MQ=255)                                                       gttgttCAGATCATCGAGGGTTTCAAACAGCGGAGCAACCGGCATCGCAAACCCGATACCCGCTTCTTTCAGCAGCAGGTGGGCAGCCAGTACGTCGGACGGCGTTTTCGCCATCGAGATCACGTAGGCGGCAATGGAg                                                                             <  1:987644/139‑1 (MQ=255)                                                       gttgttCAGATCATCGAGGGTTTCAAACAGCGGAGCAACCGGCATCGCAAACCCGATACCCGCTTCTTTCAGCAGCAGGTGGGCAGCCAGTACGTCGGACGGCGTTTTCGCCATCGAGATCACGTAGGCGGCAATGGAg                                                                             <  2:503771/139‑1 (MQ=255)                                                            tCAGATCATCGAGGGTTTCAAACAGCGGAGCAACCGGCATCGCAAACCCGATACCCGCTTCTTTCAGCAGCAGGTGGGCAGCCAGTACGTCGGACGGCGTTTTCGCCATCGAGATCACGTAGGCGGCAATGGAGCCTTg                                                                        <  1:285620/139‑1 (MQ=255)                                                                     cGAGGGTTTCAAACAGCGGAGCAACCGGCATCGCAAACCCGATACCCGCTTCTTTCAGCAGCAGGTGGGCAGCCAGTACGTCGGACg                                                                                                                   <  1:449316/87‑1 (MQ=255)                                                                     cGAGGGTTTCAAACAGCGGAGCAACCGGCATCGCAAACCCGATACCCGCTTCTTTCAGCAGCAGGTGGGCAGCCAGTACGTCGGACg                                                                                                                   >  2:449316/1‑87 (MQ=255)                                                                     cGAGGGTTTCAAACAGCGGAGCAACCGGCATCGCAAACCCGATACCCGCTTCTTTCAGCAGCAGGTGGGCAGCCAGTACGTCGGACGGCGTTTTc                                                                                                           <  1:341410/95‑1 (MQ=255)                                                                     cGAGGGTTTCAAACAGCGGAGCAACCGGCATCGCAAACCCGATACCCGCTTCTTTCAGCAGCAGGTGGGCAGCCAGTACGTCGGACGGCGTTTTc                                                                                                           >  2:341410/1‑95 (MQ=255)                                                                     cGAGGGTTTCAAACAGCGGAGCAACCGGCATCGCAAACCCGATACCCGCTTCTTTCAGCAGCAGGTGGGCAGCCAGTACGTCGGACGGCGTTTTCGCCATCGAGATCACGTAGGCGGCAATGGAGCCTTGCGGTGCTTc                                                               >  2:355769/1‑139 (MQ=255)                                                                          gTTTCAAACAGCGGAGCAACCGGCATCGCAAACCCGATACCCGCTTCTTTCAGCAGCAGGTGGGCAGCCAGTACGTCGGACGGCGTTTTCGCCATCGAGATCACGTAGGCGGCAATGGAGCCTTGCGGTGCTTCGGCaa                                                          >  1:219735/1‑139 (MQ=255)                                                                            ttCAAACAGCGGAGCAACCGGCATCGCAAACCCGATACCCGCTTCTTTCAGCAGCAGGTGGGCAGCCAGTACGTCGGACGGCGTTTTCGCCATCGAGATCACGTAGGCGGCAATGGAGCCTTGCGGTGCTTCGGTAATc                                                        >  2:29149/1‑139 (MQ=255)                                                                              cAAACAGCGGAGCAACCGGCATCGCAAACCCGATACCCGCTTCTTTCAGCAGCAGGTGGGCAGCCAGTACGTCGGACGGCGTTTTCGCCATCGAGATCACGTAGGCGGCAATGGAGCCTTGCGGTGCTTCGGTAATCAc                                                      <  1:29149/139‑1 (MQ=255)                                                                                   aGCGGAGCAACCGGCATCGCAAACCCGATACCCGCTTCTTTCAGCAGCAGGTGGGCAGCCAGTACGTCGGACGGCGTTTTCGCCATCGAGATCACGTAGGCGGCAATGGAGCCTTGCGGTGCTTCGGCAATCACCTGGc                                                 <  1:124391/139‑1 (MQ=255)                                                                                     cGGAGCAACCGGCATCGCAAACCCGATACCCGCTTCTTTCAGCAGCAGGTGGGCAGCCAGTACGTCGGACGGCGTTTTCGCCATCGAGATCACGTAGGCGGCAATGGAGCCTTGCGGTGCTTCGGCAATCACCTGGCAg                                               >  2:77318/1‑139 (MQ=255)                                                                                     cGGAGCAACCGGCATCGCAAACCCGATACCCGCTTCTTTCAGCAGCAGGTGGGCAGCCAGTACGTCGGACGGCGTTTTCGCCATCGAGATCACGTAGGCGGCAATGGAGCCTTGCGGTGCTTCGGCAATCACCTGGCAg                                               <  1:1033791/139‑1 (MQ=255)                                                                                      ggAGCAACCGGCATCGCAAACCCGATACCCGCTTCTTTCAGCAGCAGGTGGGCAGCCAGTACGTCGGACGGCGTTTTCGCCATCGAGATCACGTAGGCGGCAATGGAGCCTTGCGGTGCTTCGGCAATCACCTGGCAgg                                              >  1:1061799/1‑139 (MQ=255)                                                                                               ggCATCGCAAACCCGATACCCGCTTCTTTCAGCAGCAGGTGGGCAGCCAGTACGTCGGACGGCGTTTTCGCCATCGAGATCACGTAGGCGGCAATGGAGCCTTGCGGTGCTTCGGCAATCAc                                                      >  2:886752/1‑122 (MQ=255)                                                                                               ggCATCGCAAACCCGATACCCGCTTCTTTCAGCAGCAGGTGGGCAGCCAGTACGTCGGACGGCGTTTTCGCCATCGAGATCACGTAGGCGGCAATGGAGCCTTGCGGTGCTTCGGCAATCAc                                                      <  1:886752/122‑1 (MQ=255)                                                                                                           ccGATACCCGCTTCTTTCAGCAGCAGGTGGGCAGCCAGTACGTCGGACGGCGTTTTCGCCATCGAGATCACGTAGGCGGCAATGGAGCCTTGCGGTGCTTCGGCAATCACCTGGCAGGTATCGAGCACTTCGCGCGttt                         <  2:633304/139‑1 (MQ=255)                                                                                                               tACCCGCTTCTTTCAGCAGCAGGTGGGCAGCCAGTACGTCGGACGGCGTTTTCGCCATCGAGATCACGTAGGCGGCAATGGAGCCTTGCGGTGCTTCGGCAATCACCTGGCAGGTATCGAGCACTTCGCGCGTTTCGgc                     <  1:77318/139‑1 (MQ=255)                                                                                                                  ccGCTTCTTTCAGCAGCAGGTGGGCAGCCAGTACGTCGGACGGCGTTTTCGCCATCGAGATCACGTAGGCGGCAATGGAGCCTTGCGGTGCTTCGGCAATCACCTGGCAGGTAACGAGCACTTCGCGCGTTTCGGCGCt                  <  2:1061799/139‑1 (MQ=255)                                                                                                                                  cagGTGGGCAGCCAGTACGTCGGACGGCGTTTTCGCCATCGAGATCACGTAGGCGGCAATGGAGCCTTGCGGTGCTTCGGCAATCACCTGGCAGGTATCGAGCACTTCGCGCGTTTCGGCGCTTGGTTGCCAGTTgcgc  >  2:294319/1‑139 (MQ=255)                                                                                                                                         |                                                                                                                                    AATCAGGCCACGATACCAGTCAATATTGAGCAGCTGGGTCATGACATCGTTGGCGTTGTTCAGATCATCGAGGGTTTCAAACAGCGGAGCAACCGGCATCGCAAACCCGATACCCGCTTCTTTCAGCAGCAGGTGGACAGCCAGTACGTCGGACGGCGTTTTCGCCATCGAGATCACGTAGGCGGCAATGGAGCCTTGCGGTGCTTCGGCAATCACCTGGCAGGTATCGAGCACTTCGCGCGTTTCGGCGCTTGGTTGCCAGTTGCGC  >  NC_000913/4151506‑4151773

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