Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	NC_000913	773,693	G→A	A218T (GCT→ACT)	cydB →	cytochrome bd‑I subunit 2

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	773,693	0	G	A	100.0%	31.1 / NA	13	A218T (GCT→ACT)	cydB	cytochrome bd‑I subunit 2
Reads supporting (aligned to +/- strand):  ref base G (0/0);  new base A (9/4);  total (9/4)

ACTACACCGGTAACTTCTTCCAGTTGCTTAACCCGTTCGGCCTGCTGGCAGGCGTGGTGAGCGTAGGGATGATCATTACTCAGGGCGCAACCTATCTGCAAATGCGTACCGTGGGC‑GAACTGCACCTGCGTACCCGTGCAACGGCTCAGGTGGCTGCGCTGGTGACACTGGTCTGTTTCGCACTGGCTGGCGTATGGGTGATGTACGGTATCGATGGTTATGTCGTGAAATCGACAATGGACCATTACGCAGCCTCTAACCCACTGAATAAAGAAGTGGTTCGTGAAGCTGGCGCATGGCTGGTTAACTTCAACAACACGCCAATTCTGTGGGCTATTCCGGCACTGGGTGTGGTTCTGCCGCTGCTGACCATCCTGACTGCACGTATGGATAAAGCCGCGTGGGCGTTTGTGTTCTCC  >  NC_000913/773508‑773926                                                                                                                                                                                           |                                                                                                                                                                                                                                          actacAACGGTAACTTCTTCCAGTTGCTTAACCCGTTCGGCCTGCTGGCAGGCGTGGTGAGCGTAGGGATGATCATTACTCAGGGCGCAACCTATCTGCAAATGCGAACCGTGGGC‑GAACTGCACCTGCGTACCCGTGCAACGGCTCAGGTGGCTGCGCTGGTGACACTGGTCTGTTTCGCACTGACTGGCGTATGGGTGATGTACGGTATCGATGGTTATGTCGTGGAATCGACAATGGACCATTACGc                                                                                                                                                                           >  1:58910/1‑250 (MQ=255)          gTAACTTCTTCCAGTTGCTTAACCCGTTCGGCCTGCTGGCAGGCGTGGTGAGCGTAGGGATGATCATTACTCAGGGCGCAACCTATCTGCAAATGCGTACCGTGGGC‑GAACTGCACCTGCGTACCTGTGCAACGGCTCAGGTGGCTGCGCTGGTGGCACTGGTCTGTTTCGCACTGACTGGCGTATGGGTGATGTACgg                                                                                                                                                                                                                     <  1:147648/199‑1 (MQ=255)                               aCCCGTGCGGCCTGATGGCAGGCGTGNGGAGCGNNGANATCACCATTACTCACCGCGNCANNNANNNGNAAATGNGTNNCG‑GGGCAGANCTGCACCTGCGTACCCCTGCAACGGCCAAGGCAGCTGCGCCGGCGACAATGACNTGCTTCGCANNCACTGGCGTATGGGTGATGGANGGTNNcnngngtan                                                                                                                                                                                                         >  2:59265/1‑182 (MQ=255)                                                                                    ggCGCAACCTATCTGCAAATGCGTACCGTGGGC‑GAACTGCACCTGCGTACCCGTGCAACGGCTCAGGTGGCTGCGCTGGTGACACTGGTCTGTTTCGCACTGACTGGCGTATGGGTGATGTACGGTATCGATGGTTATGTCGTGAAATCGACAATGGACCATTACGCAGCCTCTAACCCACTGAATAAAGAAGTGGTTCGTGAAGCTGGCGCATGGCTGGTTAACTTCAacaaca                                                                                                       >  1:5024/1‑235 (MQ=255)                                                                                      cgnAACNTNNCNNCAANTGCGTACCGTGGGC‑GAACTGCNNCTGCGTACCNGTGCAACGGCTCAGGTGGCTGCGCTGGTGACACTGGTCTGTTTCGCACTGACTNGCGTATGGGNNATNTACGTNANNNANNNTTNTGTCGTGAAATCGACAATGNACNNTTACGCNGCCTCTAACCCACTGAATAAAGAAGt                                                                                                                                                <  2:58910/192‑1 (MQ=255)                                                                                                                   gC‑GAACTGCACCTGCGTACCCGTGCAACGGCTCAGGTGGCTGCGCTGGTGACACTGGTCTGTTTCGCACTGACTGGCGAATGGGTGATGTACGGTATCGATGGTTATGTCGTGAAATCGACAATGGACCATTACGCAGCCTCTAACCCACTGAATAAAGAAGTGGTTCGTGAAGCTGGCGCATGGCTGGTTAACTTCAACAACACGCCAATTCTGTGGGCTATTACGGCACTGGGTGTGGTTCTGACgctg                                                        >  1:121074/1‑251 (MQ=255)                                                                                                                                                ggCTCAGGTGGCTGCGCTGGTGACACTGGTCTGTTTCGCactgactgn                                                                                                                                                                                                                                       <  1:155921/48‑2 (MQ=255)                                                                                                                                                ggCTCAGGTGGCTGCGCTGGTGACACTGGTCTGTTTCGCActgactgg                                                                                                                                                                                                                                       >  2:155921/1‑48 (MQ=255)                                                                                                                                                 ncTCAGGTGGCTGCGCTGGTGACACTGGTCTGTTTCGCACTGACTGGCGTATGGGTGATGTACGGTATCGATGGTTATGTCGTGAAATCGACAATGGACCATTACGCAGCCTCTAACCCACTGAATAAAGAAGTGGTTCGTGAAGCt                                                                                                                                   >  1:151712/2‑147 (MQ=255)                                                                                                                                                          ncTGNGCTGGTGACCATGGTCTGTTTCGCACGGACTGGCGTATGGGTGATGTACGGTATCGATGGTTATGTCGTGAAATCGACAATGGACCATTACGCAGCCTCTAACCCACTGAATAAAGAAGTGATTCGTGAAGCt                                                                                                                                   <  2:151712/137‑1 (MQ=255)                                                                                                                                                          gNTGCGCTGGTGACACTGGTCTGTTTCGCACTGACTGGCGTATGGGTGATGTACGGTATCGATGGTTATGTCGTGAAATCGACAATGGACCATTATGCAGCCTCTATCCCACTGAATAACGAAGTGGTTCGTGAAGCTGGCGCATGGCTGGTTAAc                                                                                                                 >  1:198299/1‑156 (MQ=255)                                                                                                                                                          gCTGCGCTGGTGACACTGGTCTGTTTCGCACTGACTGGCGTATGGGTGATGTACGGTATCGATGGTTATGTCGTGAAATCGACAATGGACCATTACGCAGCCTCTACCCCACTGAATAAAGAAGTGGTTCGTGAAGCTGGCGCATGGCTGGTTAACTTCAACAACACGCCAATTCTGTGGGCTATTCCGGCACTGGGTGTGGTTCTGCCGCTGCTGACCATCCTGACTGCACGTATg                                >  1:185836/1‑237 (MQ=255)                                                                                                                                                          gCTGCGCTGGTGACACTGGTCTGTTTCGCACTGACTGGCGTATGGGTGATGTACGGTATCGATGGTTATGTCGTGAAATCGACAATGGACCATTACGCAGCCTCTAACCCACTGAATAAAGAAGTGGTTCGTGAAGCTGGCGCATGGCTGGTTAACTTCAACAACACGCCAATTCTGTGGGCTATTCCGGCACTGGGTGTGGTTCTGCCGCTGCTGACCATCCTGACTGCACGTATg                                >  1:185835/1‑237 (MQ=255)                                                                                                                                                              cgcTGGTGACACTGGTCTGTTTCGCACTGACTGGCGTATGGGTGATGTACGGTATCGATGGTTATGTc                                                                                                                                                                                                     >  1:145116/1‑68 (MQ=255)                                                                                                                                                                                      cACTGACTGGCGTATGGGTGATGTACGGTATCGATGGTTATGTCGTGAAATCGACAATGGACCATTACGCAGCCTCTAACCCACTGAATAAAGAAGTGGTTCGTGAAGCTGGCGCATGGCTGGTTAACTTCAACAACACGCCAATTCTGTGGGCTATTCCGGCACTGGGTGTGGTTCTGCCGCTGCTGACCATCCTGACTGCACGTATGGATAAAGCCGCGTGGGCGTTTGTGTtctcn  <  1:66906/239‑2 (MQ=255)                                                                                                                                                                                           |                                                                                                                                                                                                                                          ACTACACCGGTAACTTCTTCCAGTTGCTTAACCCGTTCGGCCTGCTGGCAGGCGTGGTGAGCGTAGGGATGATCATTACTCAGGGCGCAACCTATCTGCAAATGCGTACCGTGGGC‑GAACTGCACCTGCGTACCCGTGCAACGGCTCAGGTGGCTGCGCTGGTGACACTGGTCTGTTTCGCACTGGCTGGCGTATGGGTGATGTACGGTATCGATGGTTATGTCGTGAAATCGACAATGGACCATTACGCAGCCTCTAACCCACTGAATAAAGAAGTGGTTCGTGAAGCTGGCGCATGGCTGGTTAACTTCAACAACACGCCAATTCTGTGGGCTATTCCGGCACTGGGTGTGGTTCTGCCGCTGCTGACCATCCTGACTGCACGTATGGATAAAGCCGCGTGGGCGTTTGTGTTCTCC  >  NC_000913/773508‑773926

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