Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	NC_000913	447,020	G→A	A229V (GCT→GTT)	cyoE ←	heme O synthase

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	447,020	0	G	A	100.0%	27.9 / NA	16	A229V (GCT→GTT)	cyoE	heme O synthase
Reads supporting (aligned to +/- strand):  ref base G (0/0);  new base A (14/2);  total (14/2)

GCACCTTTTTAATAGAGAGGTTTTGTTACCACACAGCAGCCAGCAGCGTATGCGAGTCCGGTACCATAAAATCAACGGACATCATCACCGAGAGGGCAGTGATGGCGATGATAGAGAAGCCGAACAGCTTGCGCGCCCAGATTCTGTCATCAGCAACTTTATAACCGCGCAGAGCCATACCTAACCACCAGACGCTAACCGCCGCGGCGACCACCAGATATTTATACCCAGCGTAACCGCCAAGAGAGAGCATCAGCGTGGCAACGGCAAAGGCGATGATATACAGCGTGATGTGATTCTTCGCCACCGAAATGCCTTTTACCACTGGCAATACCGGAATGTTTGCCGCCTGGTAATCCTTAAAGCGG  >  NC_000913/446790‑447157                                                                                                                                                                                                                                       |                                                                                                                                          gCACCTTTTTAATAGAGAGGTTTTGTTACCACACAGCAGCCAGCAGCGTATGCGAGTCCGGTACCATAAAATCAACGGACATCATCACCGAGAGGGCAGTGATGGCGATGATAGAGAAGCCGAACAGCTTGCGCGCCCAGATTCTGTCATCAGCAACTTTATAACCGCGCAGAGCCATACCTAACCACCAGACGCTAACCGCCGCGGCGACCACCAGATATTTATACCCAACGTAACCGCCAAGAGAGAGc                                                                                                                       >  1:415417/1‑251 (MQ=255)                                       gccagcagcGTATGCGAGTCCGGTACCATAAACTCAACGGCCATCATCACCGCGAGGGCCGTGATGGCGATGATCGAGCAGCCGAACAGCTTGCGCGCCCAGATTCTGTCATCAGCAACTTTATAACCGCGCAGAGCCATACCTAACCACCAGACGCTAACCGCCGCGGCGACCACCAGATATTTATACCCAACGTa                                                                                                                                       >  1:179125/1‑197 (MQ=255)                                                      cttcctatcccaaaaaaaCCAACGGACATCATCCCCGAGAGGGCAGTGAGGGCGATGAAAGAGAAGCCAAAAAGCTTGCGCGCCCAGATTCTGTCATCAGCACCTTTATACCCGCGAAGAGCCATCCCAACCCCCCAAACGCAAACCGCCGCGGCGACCACCAGATATTTATACCCAACGTAACCGCCAAGAGAGAGCATCAGCGTGGCAACGGCAAAGGCGATGATATACAGCGTGATGTGAttcttc                                                                    <  1:436069/235‑1 (MQ=255)                                                        gTCCGGTACCATAAAATCAACGGACATCATCACCGAGCGGGCAGTGATGGCGATGATAGAGAAGCCGAACAGCTTGCGCGCCCAGATTCTGTCATCAGCACATTTATAACCGCGCAGAGCCATACCTACCCCCAAAACGCTAACCGCCGCGGCGACCACCCGATATTTATACCCAACGTAACCGCCAAGAGAGAGCATCAGCGTGGCAACGGCAAAGGCGATGATATACAGCGTGATGTGATTCTTCGcc                                                                 >  1:507866/1‑250 (MQ=255)                                                        gTCCGGTACCATAAAATCAACGGACATCATCACCGAGAGGGCAGTGATGGAGATGCTAGAGAAGCCGGACAGCTTGCGAGCCCAGATTCTGTCATCAGCAACTTTATAACCGCGCAGAGCCATACCTAACCACCAGACGCTAACAGCCGCGGCGACCACCAGATATTTATACACAACGTAACCACCAAGAGAGAGCATCa                                                                                                                   >  1:84236/1‑200 (MQ=255)                                                         tCCGGTACCATAAAATCAACGGACATCATCACCGAGAGGGCGGTGATGGCCATGATAGAGAAGCCGAACAGCTTGCGCGCACAGATTCTGCCACCACCAACTCTATAACCGCGCAGAGCCATACCTCACCACCAGACGCTAACCGCCGCGGCGACCACCAGATATTTATACCCAACGTAACAGCCAAGAGAGAGCATCAGCGTGGAAACGGCAAAGGCGATGATATACAGCGTGATGTGATTCTTCGCCa                                                                >  1:236764/1‑250 (MQ=255)                                                                      aaccaccgaaCATCATCCCCGAGAGGGCATTGATGGCGAGGATAGAGAACCCGAACACCTGGCGCGCCCAGTTTCTGCAATCACCAATTTTAAACCCGCGCAGACCCATACCTAACCTCCAACCGCTAACCGCCGCGGCGACCACCATATATGTATACCCAACATAACCGCCAAGAGAGAGCATCAGCGGGGCAACGGCAAAGGCGATGATATACAGCGTGATGTGATTCTTCGCCACCGAAATGCCtttt                                                  <  1:192645/242‑1 (MQ=255)                                                                                                          cGATGATAGAGAAGCCGAACAGCTTGCGCGCCCAGATTCTGTCATCAGCAACTTTATAACCGCGCAGCGCAATAACTCCACACCAGACGCTAACCGCCGCGGAGACCACCAGATATTTATACCCAACGTAACCGCCAAGAGAGAGCATCAGCGTGGCAACGGCCAAGGCGATGATATACAGCGTGATGTGATTCTTCGCCACCGAAATGCCTTTTACCACTGGCAATACCGGAATGTTTGCCGCCTGGTa               >  1:369064/1‑250 (MQ=255)                                                                                                                       nccGAACAGCTTGCGCGCCCAGATTCTGTCATCAGCAACTTTATAACCGCGCAGAGCCATACCTAACCACCAGACGCTAACCGCCGCGGCGACCACCAGATATTTATACCCCACGTAACCGCCAAGAGAGAGCATCAGCGTGGCAACGGCCAAGGTGATGATATACAGCGTGATGTGATTCTTCGCCACCGAAATGCCTTTTACCACTGGCAATACCGGACTGTTTGCCGCCTGGTAATCCTTAAAGCgg  >  1:105984/2‑250 (MQ=255)                                                                                                                         cGAACAGCTTGCGCGCCCAGATTCTGTCATCAGCAACTTTATAACCGCGCAGAGCCATACCTAACCACCAGACGCTAACCGCCGCGGCGACCACCAGATATTTATACCCAACGTAACCGCCAAGAGAGAGCATNAGCNNNNNaann                                                                                                        >  2:305207/1‑144 (MQ=255)                                                                                                                                        cccAGATTCTGTCATCAGCAACTTTATAACCGCGCAGAGCCATACCTAACCACCAGACGCTAACCGCCGCGGCGCCCACCAGATATTTATCCCCAACGCAACCGCCAAGAGAGAGCATCAGCGTGGCAACGGCAAAGGCGATGATATACAGCGTGATGTGATTCTTCGCCACCGAAATGCCTTTTACCACTGGCCATACCGGAATGTTTg                         >  1:431972/1‑210 (MQ=255)                                                                                                                                                       caGCAACTTTATAACCGCGCAGAGCCATACCTAACCACCAGACGCTAACCGCCGCGGCGACCACCAGATATTTATACCCAACGTAACCNNCAagagnn                                                                                                                          >  2:27046/1‑96 (MQ=255)                                                                                                                                                                  taACCGCGCAGAGCCATACCTAACCACCAGACGCTAACCGCCGCGGCGACCACCAGATATTTATACCCAACGTAACCGCCAAGAGAGAGCATCAGCGTGGCAACGGCAAAGGCGATGATATACAGCGTGATGTGATTCTTCGCCACCGAAATGCCTTTTACCACTGGCAATACCGGAATGTTTg                         >  1:285223/1‑184 (MQ=255)                                                                                                                                                                          cAGAGCCATACCTAACCACCAGACGCNAACCGCNNCGNCGCCCACCAGATATTTATANCCNNNGNNNCNGCCAANAGNNAGCATCAGCNTGGCAACGGCAAAGGCGATGATATACAGCGTGATGTGATTCTTCGCCACCGAANTGCCTTTTANNACAGGCAATACCg                                  >  2:120842/1‑167 (MQ=255)                                                                                                                                                                                            ncaGACGCTAACCGCCGCGGCGACCACCAGATATTTATACCCAACGTAACCGCCAAGAGAGAGCATCAGCGTGGCAACGGCCAAGGCGATGATATACAGCGTGATGTGATTCTTCGCCACCGAAATGc                                                       >  1:93336/2‑128 (MQ=255)                                                                                                                                                                                                     aaCCGCCGCGGCGACCACCAGATATTTATACCCAACGTAACCGCCAAGAGAGCGCATACGCGCGGCAACGGCCAAGGCGATGCTATACAGCGAGATGTGATTCTTCGCCACCGCAAt                                                         >  1:245031/1‑117 (MQ=255)                                                                                                                                                                                                                ngACCACCAGATATTTATACCCAACGTAACCGCCAAGAGAGAGCATCAGCGTGGCAACGACAAAGGCGATGATCTACAGCGTGATGTGATTCTTCGCCACCGAAAt                                                         >  1:58052/2‑106 (MQ=255)                                                                                                                                                                                                                                       |                                                                                                                                          GCACCTTTTTAATAGAGAGGTTTTGTTACCACACAGCAGCCAGCAGCGTATGCGAGTCCGGTACCATAAAATCAACGGACATCATCACCGAGAGGGCAGTGATGGCGATGATAGAGAAGCCGAACAGCTTGCGCGCCCAGATTCTGTCATCAGCAACTTTATAACCGCGCAGAGCCATACCTAACCACCAGACGCTAACCGCCGCGGCGACCACCAGATATTTATACCCAGCGTAACCGCCAAGAGAGAGCATCAGCGTGGCAACGGCAAAGGCGATGATATACAGCGTGATGTGATTCTTCGCCACCGAAATGCCTTTTACCACTGGCAATACCGGAATGTTTGCCGCCTGGTAATCCTTAAAGCGG  >  NC_000913/446790‑447157

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