Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	NC_000913	3,542,962	T→C	*79Q (TAA→CAA)	feoC →	ferrous iron transport protein FeoC

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	3,542,962	0	T	C	100.0%	139.6 / NA	40	*79Q (TAA→CAA)	feoC	ferrous iron transport protein FeoC
Reads supporting (aligned to +/- strand):  ref base T (0/0);  new base C (21/19);  total (21/19)

CATGCTGCAACAACTGGAAAGTATGGGCAAAGCCGTGCGGATTCAGGAAGAACCTGACGGCTGCCTCTCTGGCAGTTGTAAAAGCTGCCCGGAAGGAAAAGCCTGTCTGCGCGAGTGGTGGGCGCTGCGTTAACCTTACTCCATCGCCTGTTTTGGAAAGCCCGGTATGCGTCTGCATCCGGGCTTTTTTGCGTGCGGCTTTCCATAAAAATGCAACTCTTGCAGCACGGCGTAAGTTCCTTTGAAAGCATCTCGCAGGGATGAAAA  >  NC_000913/3542832‑3543098                                                                                                                                   |                                                                                                                                         cATGCTGCAACAACTGGAAAGTATGGGCAAAGCCGTGCGGATTCAGGAAGAACCTGACGGCTGCCTCTCTGGCAGTTGTAAAAGCTGCCCGGAAGGAAAAGCCTGTCTGCGCGAGTGGTGGGCGCTGCGTCAACCTTACt                                                                                                                                 >  1:528269/1‑140 (MQ=255)   tgctgcAACAACTGGAAAGTATGGGCAAAGCCGTGCGGATTCAGGAAGAACCTGACGGCTGCCTCTCTGGCAGTTGTAAAAGCTGCCCGGAAGGAAAAGCGTGTGTGCGCGAGTGGTGGGCGCTGCGTCAACCTTACTcc                                                                                                                               >  2:289634/1‑140 (MQ=255)   tgctgcAACAACTGGAAAGTATGGGCAAAGCCGTGCGGATTCAGGAAGAACCTGACGGCTGCCTCTCTGGCAGTTGTAAAAGCTGCCCGGAAGGAAAAGCCTATCTGCGCGAGTGGTGGGCGCTGCGTCAACCTTACTcc                                                                                                                               >  1:490906/1‑140 (MQ=255)        caacaaCTGGAAAGTATGGGCAAAGCCGTGCGGATTCAGGAAGAACCTGACGGCTGCCTCTCTGGCAGTTGTAAAAGCTGCCCGGAAGGAAAAGCCTGTCTGCGCGAGTGGTGGGCGCTGCGTCAACCTTACTCCATCGc                                                                                                                          <  2:906952/140‑1 (MQ=255)        caacaaCTGGAAAGTATGGGCAAAGCCGTGCGGATTCAGGAAGAACCTGACGGCTGCCTCTCTGGCAGTTGTAAAAGCTGCCCGGAAGGAAAAGCCTGTCTGCGCGAGTGGTGGGCGCTGCGTCAACCTTACTCCATCGc                                                                                                                          >  1:523766/1‑140 (MQ=255)        caacaaCTGGAAAGTATGGGCAAAGCCGTGCGGATTCAGGAAGAACCTGACGGCTGCCTCTCTGGCAGTTGTAAAAGCTGCCCGGAAGGAAAAGCCTGTCTGCGCGAGTGGTGGGCGCTGCGTCAACCTTACTCCATCGc                                                                                                                          >  2:234149/1‑140 (MQ=255)        caacaaCTGGAAAGTATGGGCAAAGCCGTGCGGATTCAGGAAGAACCTGACGGCTGCCTCTCTGGCAGTTGTAAAAGCTGCCCGGAAGGAAAAGCCTGTCTGCGCGAGTGGTGGGCGCTGCGTCAACCTTACTCCATCGc                                                                                                                          <  2:696858/140‑1 (MQ=255)        caacaaCTGGAAAGTATGGGCAAAGCCGTGCGGATTCAGGAAGAACCTGACGGCTGCCTCTCTGGCAGTTGTAAAAGCTGCCCGGAAGGAAAAGCCTGTCTGCGCGAGTGGTGGGCGCTGCGTCAACCTTACTCCATCGc                                                                                                                          <  2:195304/140‑1 (MQ=255)                   aaGTATGGGCAAAGCCGTGCGGTTTCAGGAAGAACCTCACGGCTGCCTCTCTGGCAGTTGTAAAAGCTGCCCGGAAGGAAAAGCCTGTCTGCGCGAGTGGTGGGCGCTGCGTCAACCTTACTCCATCGCCTGTTTTGGaa                                                                                                               <  1:373309/140‑1 (MQ=255)                       aTGGGCAAGGCCGTGCCCATTCCGGAAGAATTTTACGGCTGCCTCTCTCGCAGTTGTAAAAGCTGCCCGGAAGGAAAAGCCTGTCTGCGCGAGTGGTGGGCGCTGCGTCAACCTTACTCCATCGCCTGTTTTGGAAAGcc                                                                                                           <  2:490906/140‑1 (MQ=255)                        tGGGCAAAGCCGTGCGGATTCAGGAAGAACCTGACGGCTGCCTCTCTGGCAGTTGTAAAAGCTGCCCGGAAGGAAAAGCCTGTCTGCGCGAGTGGTGGGCGCTGCGTCAACCTTACTCCATCGCCTGTTTTGGAAAGccc                                                                                                          >  2:985288/1‑140 (MQ=255)                        tGGGCAAAGCCGTGCGGATTCAGGAAGAACCTGACGGCTGCCTCTCTGGCAGTTGTAAAAGCTGCCCGGAAGGAAAAGCCTGTCTGCGCGACTGGTGGGCTCTGAGTCGACCTTACTCCATCGCCTGTTTTGGAAAGccc                                                                                                          >  2:340671/1‑140 (MQ=255)                         gggCAAAGCCGTGCGGATTCAGGAAGAACCTGACGGCTGCCTCTCTGGCAGTTGTAAAAGCTGCCCGGAAGGAAAAGCCTGTCTGCGCGAGTGGTGGGCGCTGCGTCAACCTTACTCCATCGc                                                                                                                          >  1:72250/1‑123 (MQ=255)                         gggCAAAGCCGTGCGGATTCAGGAAGAACCTGACGGCTGCCTCTCTGGCAGTTGTAAAAGCTGCCCGGAAGGAAAAGCCTGTCTGCGCGAGTGGTGGGCGCTGCGTCAACCTTACTCCATCGc                                                                                                                          <  2:72250/123‑1 (MQ=255)                           gCAAAGCCGTGCGGATTCAGGAAGAACCTGACGGCTGCCTCTCTGGCAGTTGTAAAAGCTGCCCGGAAGGAAAAGCCTGTCTGCGCGAGTGGTGGGCGCTGCGTCAACCTTACTCCATCGCCTGTTTTGGAAAGCCCGGt                                                                                                       >  2:814225/1‑140 (MQ=255)                             aaaGCCGTGCGGATTCAGGAAGAACCTGACGGCTGCCTCTCTGGCAGTTGTAAAAGCTGCCCGGAAGGAAAAGCCTGTCTGCGCGAGTGGTGGGCGCTGCGTCAACCTTACTCCATCGCCTGTTTTGGAAAGCCCGGTAt                                                                                                     <  2:76054/140‑1 (MQ=255)                                    tGCGGATTCAGGAAGAACCTGACGGCTGCCTCTCTGGCAGTTGTAAAAGCTGCCCGGAAGGAAAAGCCTGTCTGCGCGAGTGGTGGGCGCTGCGTCAACCTTACTCCATCGCCTGTTTTGGAAAGCCCGGTATGCGTCTg                                                                                              >  1:766290/1‑140 (MQ=255)                                    tGCGGATTCAGGAAGAACCTGACGGCTGCCTCTCTGGCAGTTGTAAAAGCTGCCCGGAAGGAAAAGCCTGTCTGCGCGAGTGGTGGGCGCTGCGTCAACCTTACTCCATCGCCTGTTTTGGAAAGCCCGGTATGCGTCTg                                                                                              >  1:620163/1‑140 (MQ=255)                                    tGCGGATTCAGGAAGAACCTGACGGCTGCCTCTCTGGCAGTTGTAAAAGCTGCCCGGAAGGAAAAGCCTGTCTGCGCGAGTGGTGGGCGCTGCGTCAACCTTACTCCATCGCCTGTTTTGGAAAGCCCGGTATGCGTCTg                                                                                              >  2:559683/1‑140 (MQ=255)                                                   aaCCTGACGGCTGCCTCTCTGGCAGTTGTAAAAGCTGCCCGGAAGGAAAAGCCTGTCTGCGCGAGTGGTGGGCGCTGCGTCAACCTTACTCCATCGCCTGTTTTGGAAAGCCCGGTATGCGTCTGCATCCGGGCtttttt                                                                               >  2:727295/1‑140 (MQ=255)                                                             cTGCCTCTCTGGCAGTTGTAAAAGCTGCCCGGAAGGAAAAGCCTGTCTGCGCGAGTGGTGGGCGCTGCGTCAACCTTACTCCATCGCCTGTTTTGGAAAGCCCGGTATGCGTCTGCATCCGGGCTTTTTTGCGTGCGGCt                                                                     >  2:647094/1‑140 (MQ=255)                                                             cTGCCTCTCTGGCAGTTGTAAAAGCTGCCCGGAAGGAAAAGCCTGTCTGCGCGAGTGGTGGGCGCTGCGTCAACCTTACTCCATCGCCTGTTTTGGAAAGCCCGGTATGCGTCTGCATCCGGGCTTTTTTGCGTGCGGCt                                                                     <  1:559683/140‑1 (MQ=255)                                                               gCCTCTCTGGCAGTTGTAAAAGCTGCCCGGAAGGAAAAGCCTGTCTGCGCGAGTGGTGGGCGCTGCGTCAACCTTACTCCATCGCCTGTTTTGGAAAGCCCGGTATGCGTCTGCATCCGGGCTTTTTTGCGTGCGGCttt                                                                   <  2:528269/140‑1 (MQ=255)                                                                ccTCTCTGGCAGTTGTAAAAGCTGCCCGGAAGGAAAAGCCTGTCTGCGCGAGTGGTGGGCGCTGCGTCAACCTTACTCCATCGCCTGTTTTGGAAAGCCCGGTATGCGTCTGCATCCGGGCTTTTTTGCGTGCGGCTTTc                                                                  <  2:523766/140‑1 (MQ=255)                                                                           gTTGTAAAAGCTGCCCGGAAGGAAAAGCCTGTCTGCGCGAGTGGTGGGCGCTGCGTCAACCTTACTCCATCGCCTGTTTTGGAAAGCCCGGTATGCGTCTGCATCCGGGCTTTTTTGCGTGCGGCTTTCCATAAAAATGc                                                       >  2:1009298/1‑140 (MQ=255)                                                                            ttGTAAAAGCTGCCCGGAAGGAAAAGCCTGTCTGCGCGAGTGGTGGGCGCTGCGTCAACCTTACTCCATCGCCTGTTTTGGAAAGCCCGGTATGCGTCTGCATCCGGGCTTTTTTGCGTGCGGCTTTCCATAAAAATGCa                                                      <  1:127386/140‑1 (MQ=255)                                                                               tAAAAGCTGCCCGGAAGGAAAAGCCTGTCTGCGCGAGTGGTGGGCGCTGCGTCAACCTTACTCCATCGCCTGTTTTGGAAAGCCCGGTATGCGTCTGCATCCGGGCTTTTTTGCATGGGGCTTTACATAAAAATGCAAct                                                   >  1:263513/1‑140 (MQ=255)                                                                                      tGCCCGGAAGGAAAAGCCTGTCTGCGCGAGTGGTGGGCGCTGCGTCAACCTTACTCCATCGCCTGTTTTGGAAAGCCCGGTATGCGTCTGCATCCGGGCTTTTTTGCGTGCGGCTTTCCATAAAAATGCAACTCTTgcag                                            >  1:981631/1‑140 (MQ=255)                                                                                       gCCCGGAAGGAAAAGCCTGTCTGCGCGAGTGGTGGGCGCTGCGTCAACCTTACTCCATCGCCTGTTTTGGAAAGCCCGGTATGCGTCTGCATCCGGGCTTTTTTGCGTGCGGCTTTCCATAAAAATGCAACTCTTgcagc                                           <  2:205187/140‑1 (MQ=255)                                                                                         ccGGAAGGAAAAGCCTGTCTGCGCGAGTGGTGGGCGCTGCGTCAACCTTACTCCATCGCCTGTTTTGGAAAGCCCGGTATGCGTCTGCATCCGGGCTTTTTTGCGTGCGGCTTTCCATAAAAATGCAACTCTTGCAGCAc                                         >  1:580756/1‑140 (MQ=255)                                                                                         ccGGAAGGAAAAGCCTGACTGCCCGAGTGGTCGGCCCTGTGTCAACCTCACTCCATCGCCTGTTTTGGAAAGCCCGGTATGCGTCTGCATCCGGGCTTTTTTGCGTGCGGCTTTCCATAAAAATGCAACTCTTGCAGCAc                                         <  2:263513/140‑1 (MQ=255)                                                                                                 aaaaGCCTGTCTGCGCGAGTGGTGGGCGCTGCGTCAACCTTCCTCCATCGCCTGTTTTGGAAAGCCCGGTATGCGTCTGCATCCGGGCTTTTTTGCGTGCGGCTTTCCATAAAAATGCAACTCTTGCAGCACGGCGTAAg                                 <  2:825085/140‑1 (MQ=255)                                                                                                   aaGCCTGTCTGCGCGAGTGGTGGGCGCTGCGTCAACCTTACTCCATCGCCTGTTTTGGAAAGCCCGGTATGCGTCTGCATCCGGGCTTTTTAGCGTGCGGATTTTCATAAAAATGCAACTCTGCCAGCACGACGTAAGtt                               >  2:460740/1‑140 (MQ=255)                                                                                                         gTCTGCGCGAGTGGTGGGCGCTGCGTCAACCTTACTCCATCGCCTGTTTTGGAAAGCCCGGTATGCGTCTGCATCCGGGCTTTTTTGCGTGCGGCTTTCCATAAAAATGCTACTTTTTCAGCACGGCGTAAGTTCCTTTg                         >  1:420733/1‑140 (MQ=255)                                                                                                          tCTGCGCGAGTGGTGGGCGCTGCGTCAACCTTATTCCATCGCCTGTTTTGGAAAGCCCGGTATGCGTCTGCATCCGGGCTTTTTTGCGTGCGGCTTTCCATAAAAATGCAACTCTTGCAGCACGGCGTAAGTTCCTTTGa                        <  1:943124/140‑1 (MQ=255)                                                                                                               gcgAGTGGTGGGCGCTGCGTCATCCTTACTCCATCTCCTTTTTTGGAACGCCCGGTATGCGTCTGCATCCGGGCTTTTTTGCGTGCGGCTTTCCATAAAAATGCAACTCTTGCAGCACGGCGTAAGTTCCTTTGAAAGCa                   <  2:420733/140‑1 (MQ=255)                                                                                                                   gtggtgGGCGCTGCGTCAACCTTACTCCATCGCCTGTTTTGGAAAGCCCGGTATGCGTCTGCATCCGGGCTTTTTTGCGTGCGGCTTTCCATAAAAATGCAACTCTTGCAGCACGGCGTAAGTTCCTTTGAAAGCAtctc               >  2:91585/1‑140 (MQ=255)                                                                                                                    tggtggGCGCTGCGTCAACCTTACTCCATCGCCTGTTTTGGAAAGCCCGGTATGCGTCTGCATCCTGGCTTTTTTGCGTGCGGCTTTCCATAAAAATGCAACTCTTGCAGCACGGCGTAAGTTCCTTTGAAAGCATCTCg              <  2:865844/140‑1 (MQ=255)                                                                                                                           cgcTGCGTCAACCTTACTCCATCGCCTGTTTTGGAAAGCCCGGTATGCGTCTGCATCCGGGCTTTTTTGCGTGCGGCTTTCCATAAAAATGCAACTCTTGCAGCACGGCGTAAGTTCCTTTGAAAGCATCTCGCAGGGAt       <  1:727295/140‑1 (MQ=255)                                                                                                                                cGTCAACCTTACTCCCTCGCCTGTTGTGGAAAGCCCGGTATGCGTCTGCATCCGGGCTTTTTTGCGTGCGGCTTTCCATAAAAATGCAACTCTTGCAGCACGGCGTAAGTTCCTTTGAAAGCATCTCGCAGGGATGaaaa  <  1:412871/140‑1 (MQ=255)                                                                                                                                   |                                                                                                                                         CATGCTGCAACAACTGGAAAGTATGGGCAAAGCCGTGCGGATTCAGGAAGAACCTGACGGCTGCCTCTCTGGCAGTTGTAAAAGCTGCCCGGAAGGAAAAGCCTGTCTGCGCGAGTGGTGGGCGCTGCGTTAACCTTACTCCATCGCCTGTTTTGGAAAGCCCGGTATGCGTCTGCATCCGGGCTTTTTTGCGTGCGGCTTTCCATAAAAATGCAACTCTTGCAGCACGGCGTAAGTTCCTTTGAAAGCATCTCGCAGGGATGAAAA  >  NC_000913/3542832‑3543098

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