Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	NC_000913	3,542,997	T→C	intergenic (+33/‑170)	feoC → / → rpnA	ferrous iron transport protein FeoC/recombination‑promoting nuclease RpnA

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	3,542,997	0	T	C	100.0%	54.8 / NA	17	intergenic (+33/‑170)	feoC/rpnA	ferrous iron transport protein FeoC/recombination‑promoting nuclease RpnA
Reads supporting (aligned to +/- strand):  ref base T (0/0);  new base C (8/9);  total (8/9)

AGCCGTGCGGATTCAGGAAGAACCTGACGGCTGCCTCTCTGGCAGTTGTAAAAGCTGCCCGGAAGGAAAAGCCTGTCTGCGCGAGTGGTGGGCGCTGCGTTAACCTTACTCCATCGCCTGTTTTGGAAAGCCCGGTATGCGTCTGCATCCGGGCTTTTTTGCGTGCGGCTTTCCATAAAAATGCAACTCTTGCAGCACGGCGTAAGTTCCTTTGAAAGCATCTCGCAGGGATGAAAACTCGCTAATACACAGGTGTGGAG  >  NC_000913/3542862‑3543121                                                                                                                                        |                                                                                                                             aGCCGTGCGGATTCAGGAAGAACCTGACGGCTGCCTCTCTGGCAGTTGTAAAAGCTGCCCGGAAGGAAAAGCCTGTCTGCGCGAGTGGTGGGCGCTGCGTTAACCTTACTCCATCGCCTGTTTTGGAAAGCCCGGCATGc                                                                                                                          <  1:186142/140‑1 (MQ=255)  gCCGTGCGGATTCAGGAAGAACCTGACGGCTGCCTCTCTGGCAGTTGTAAAAGCTGCCCGGAAGGAAAAGCCTGTCTGCGCGAGTGGTGGGCGCTGCGTTAACCTTACTCCATCGCCTGTTTTGGAAAGCCCGGCATGCg                                                                                                                         >  1:50875/1‑140 (MQ=255)                                                     aGCTGCCCGGAAGGAAAAGCCTGTCTGCGCGAGTGGTGGGCGCTGCGTTAACCTTACTCCATCGCCTGTTTTGGAAAGCCCGGCATGCTTCTGCATCCGGTCTTTTTTGCGTGCGGCTTTCCATAAAAATGCAACTCTTg                                                                      <  2:264959/140‑1 (MQ=255)                                                     aGCTGCCCGGAAGGAAAAGCCTGTCTGCGCGAGTGGTGGGCGCTGCGTTAACCTTACTCCATCGCCTGTTTTGGAAAGCCCGGCATGCGTCTGCATCCGGGCTTTTTTGCGTGCGGCTTTCCATAAAAATGCAACTCTTg                                                                      <  2:52278/140‑1 (MQ=255)                                                          cccGGAAGGAAAAGCCTGTCTGCGCGAGTGGTGGGCGCTGCGTTAACCTTACTCCATCGCCTGTTTTGGAAAGCCCGGCATGCGTCTGCATCCGGGCTTTTTTGCGTGCGGCTTTCCATAAAAATGCAACTCTTgcagca                                                                 <  2:50875/140‑1 (MQ=255)                                                           ccGGAAGGAAAAGCCTGTCTGCGCGAGTGGTGGGCGCTGCGTTAACCTTACTCCATCGCCTGTTTTGGAAAGCCCGGCATGCGTCTGCATCCGGGCTTTTTTGCGTGCGGCTTTCCATAAAAATGCAACTCTTGCAGCAc                                                                >  2:364345/1‑140 (MQ=255)                                                                                                     tAACCTTACTCCATCGCCTGTTTTGGAAAGCCCGGCATGCGTCTGCATCCGGGCTTTTTTGCGTGCGGCTTTCCATAAAAATGCAACTCTTGCAGCACGGCGTAAGTTCCTTTGAAAGCATCTCGCAGGGATGAAAACTc                      >  1:97481/1‑140 (MQ=255)                                                                                                       aCCTTACTCCATCGCCTGTTTTGGAAAGCCCGGCATGCGTCTGCATCCGGGCTTTTTTGCGTGCGGCTTTCCATAAAAATGCAACTCTTGCAGCACGGCGTAAGTTCCTTTGAAAGCATCTCGCAGGGATGAAAACTCGc                    <  1:241281/140‑1 (MQ=255)                                                                                                            aCTCCATCGCCTGTTTTGGAAAGCCCGGCATGCGTCTGCATCCGGGCTTTTTTGCGTGCGGCTTTCCATAAAAATGCAACTCTTGCAGCACGGCGTAAGTTCCTTTGAAAGCATCTCGCAGGGATGAAAACTCg                     >  2:100278/1‑134 (MQ=255)                                                                                                            aCTCCATCGCCTGTTTTGGAAAGCCCGGCATGCGTCTGCATCCGGGCTTTTTTGCGTGCGGCTTTCCATAAAAATGCAACTCTTGCAGCACGGCGTAAGTTCCTTTGAAAGCATCTCGCAGGGATGAAAACTCg                     <  1:100278/134‑1 (MQ=255)                                                                                                            aCTCCATCGCCTGTTTTGGAAAGCCCGGCATGCGTCTGCATCCGGGCTTTTTTGCGTGCGGCTTTCCATAAAAATGCAACTCTTGCAGCACGGCGTAAGTTCCTTTGAAAGCATCTCGCAGGGATGAAAACTCGCTAATa               >  1:425280/1‑140 (MQ=255)                                                                                                              tCCATCGCCTGTTTTGGAAAGCCCGGCATGCGTCTGCATCCGGGCTTTTTTGCGTGCGGCTTTCCATAAAAATGCAACTCTTGCAGCACGGCGTAAGTTCCTTTGAAAGCATCTCGCAGGGATGAAAACTCGCTAATaca             <  2:224645/140‑1 (MQ=255)                                                                                                                  tCGCCTGTTTTGGAAAGCCCGGCATGCGTCTGCATCCGGGCTTTTTTGCGTGCGGCTTTCCGTAAAAATGCAACTCTTGCAGCACGGCGTAAGTTCCTTTGAAAGCATCTCGCAGGGATGAAAACTCGCTAATaca             >  1:301084/1‑136 (MQ=255)                                                                                                                  tCGCCTGTTTTGGAAAGCCCGGCATGCGTCTGCATCCGGGCTTTTTTGCGTGCGGCTTTCCGTAAAAATGCAACTCTTGCAGCACGGCGTAAGTTCCTTTGAAAGCATCTCGCAGGGATGAAAACTCGCTAATaca             <  2:301084/136‑1 (MQ=255)                                                                                                                  tCGCCTGTTTTGGAAAGCCCGGCATGCGTCTGCATCCGGGCTTTTTTGCGTGCGGCTTTCCATAAAAATGCAACTCTTGCAGCACGGCGTAAGTTCCTTTGAAAGCATCTCGCAGGGATGAAAACTCGCTAATACACAgg         >  1:472529/1‑140 (MQ=255)                                                                                                                  tCGCCTGTTTTGGAAAGCCCGGCATGCGTCTGCATCCGGGCTTTTTTGCGTGCGGCTTTCCATAAAAATGCAACTCTTGCAGCACGGCGTAAGTTCCTTTGAAAGCATCCCGCAGGGATGAAAACTCGCTAATACACAgg         >  1:288297/1‑140 (MQ=255)                                                                                                                         ttttGGAAAGCCCGGCATGCGTCTGCATCCGGGCTTTTTTGCGTGCGGCTTTCCATAAAAATGCAACTCTTGCAGCACGGCGTAAGTTCCTTTGAAAGCATCTCGCAGGGATGAAAACTCGCTAATACACAGGTGTGGAg  <  2:509911/140‑1 (MQ=255)                                                                                                                                        |                                                                                                                             AGCCGTGCGGATTCAGGAAGAACCTGACGGCTGCCTCTCTGGCAGTTGTAAAAGCTGCCCGGAAGGAAAAGCCTGTCTGCGCGAGTGGTGGGCGCTGCGTTAACCTTACTCCATCGCCTGTTTTGGAAAGCCCGGTATGCGTCTGCATCCGGGCTTTTTTGCGTGCGGCTTTCCATAAAAATGCAACTCTTGCAGCACGGCGTAAGTTCCTTTGAAAGCATCTCGCAGGGATGAAAACTCGCTAATACACAGGTGTGGAG  >  NC_000913/3542862‑3543121

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