Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	NC_000913	4,513,105	G→A	A102V (GCG→GTG)	fecC ←	ferric citrate ABC transporter membrane subunit FecC

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	4,513,105	0	G	A	100.0%	99.3 / NA	31	A102V (GCG→GTG)	fecC	ferric citrate ABC transporter membrane subunit FecC
Reads supporting (aligned to +/- strand):  ref base G (0/0);  new base A (18/13);  total (18/13)

TATGACGAAATCCGCCTCCTGCGGTCATGACCAGCAGCCAGCTCACGCCGCCCCCGCATGCCGCGATGAATGACAGAGAATAGCCTGCAATCGGCGTCGGACTCAGCGCGCTGGTAAGCGCCATAGCCAGCGCCGCGCCGCTGTTAATGCCGAGCAGTGAAGGAGAGGCCATTGGGTTGTGGGTCAGGGTTTGCAGCAGCGTGCCCGCGAGCGCCAGGCTTGCGCCGATCAGAACGGCGACCAGGCTTCGTGGCAAACGA  >  NC_000913/4512974‑4513233                                                                                                                                    |                                                                                                                                 tATGACGAAATCCGCCTCCTGCGGTCATGACCAGCAGCCAGCTCACGCCGCCCCCGCATGCCGCGATGAATGACAGAGAATAGCCTGCAATCGGCGTCGGACTCAGCGCGCTGGTAAGCGCCATAGCCAGcaccgcgccg                                                                                                                          <  2:664697/140‑1 (MQ=255)    gACGAAATCCGCCTCCTGCGGTCATGACCAGCAGCCAGCTCACGCCGCCCCCGCATGCCGCGATGAATGACAGAGAATAGCCTGCAATCGGCGTCGGACTCAGCGCGCTGGTAAGCGCCATAGCCAGCACCGCGCCGCTg                                                                                                                       >  1:576729/1‑140 (MQ=255)         aaTCCGCCTCCTGCGGTCATGACCAGCAGCCAGCTCACGCCGCCCCCGCATGCCGCGATGAATGACAGAGAATAGCCTGCAATCGGCGTCGGACTCAGCGCGCTGGTAAGGGCCATAACCAGCACCGCGCCGCTGTTAAt                                                                                                                  >  1:270774/1‑140 (MQ=255)         aaTCCGCCTCCTGCGGTCATGACCAGCAGCCAGCTCACGCCGCCCCCGCATGCCGCGATGAATGACAGAGAATAGCCTGCAATCGGCGTCGGACTCAGCGCGCTGGTAAGCGCCATAGCCAGCACCGCGCCGCTGTTAAt                                                                                                                  >  2:397686/1‑140 (MQ=255)          aTCCGCCTCCTGCGGTCATGACCAGCAGCCAGCTCACGCCGCCCCCGCATGCCGCGATGAATGACAGAGAATAGCCTGCAATCGGCGTCGGACTCAGCGCGCTGGTAAGCGCCATAGCCAGCACCGCGCCGCTGTTAATg                                                                                                                 <  2:683833/140‑1 (MQ=255)          aTCCGCCTCCTGCGGTCATGACCAGCAGCCAGCTCACGCCGCCCCCGCATGCCGCGATGAATGACAGAGAATAGCCTGCAATCGGCGTCGGACTCAGCGCCCTGCTAAGCGCCATAGCCAGCACCGCGCCGCTGTTAATg                                                                                                                 >  2:253581/1‑140 (MQ=255)           tCCGCCTCCTGCGGTCATGACCAGCAGCCAGCTCACGCCGCCCCCGCATGCCGCGATGAATGACAGAGAATAGCCTGCAATCGGCGTCGGACTCAGCGCGCTGGTAAGCGCCATAGCCAGCACCGCGCCGCTGTTAATGc                                                                                                                >  1:420430/1‑140 (MQ=255)               cctcctGCGGTCATGACCAGCAGCCAGCTCACGCCGCCCCCGCATGCCGCGATGAATGACAGAGAATAGCCTGCAATCGGCGTCGGACTCAGCGCGCTGGTAAGCGCCATAGCCAGCACCGCGCCGCTGTTAATGCCGAg                                                                                                            <  1:75015/140‑1 (MQ=255)                 tcctGCGGTCATGACCAGCAGCCAGCTCACGCCGCCCCCGCATGCCGCGATGAATGACAGAGAATAGCCTGCAATCGGCGTCGGACTCAGCGCGCTGGTAAGCGCCATAGCCAGCACCGCGCCGCTGTTAATGCCGAg                                                                                                            <  1:149586/138‑1 (MQ=255)                 tcctGCGGTCATGACCAGCAGCCAGCTCACGCCGCCCCCGCATGCCGCGATGAATGACAGAGAATAGCCTGCAATCGGCGTCGGACTCAGCGCGCTGGTAAGCGCCATAGCCAGCACCGCGCCGCTGTTAATGCCGAg                                                                                                            >  2:149586/1‑138 (MQ=255)                 tcctGCGGTCATGACCAGCAGCCAGCTCACGCCGCCCCCGCATGCCGCGATGAATGACAGAGAATAGCCTGCAATCGGCGTCGGACTCAGCGCGCTGGTAAGCGCCATAGCCAGCACCGCGCCGCTGTTAATGCCGAGCa                                                                                                          >  1:376497/1‑140 (MQ=255)                          cATGACCAGCAGCCAGCTCACGCCGCCCCCGCATGCCGCGATGAATGACAGAGAATAGCCTGCAATCGGCGTCGGACTCAGCGCGCTGGTAAGCGCCATAGCCAGCACCGCGCCGCTGTTAATGCCGAGCAGTGAAGgag                                                                                                 >  1:734991/1‑140 (MQ=255)                                   cagccagcTCACGCCGCCCCCGCATGCCGCGATGAATGACAGAGAATAGCCTGCAATCGGCGTCGGACTCAGCGCGCTGGTAAGCGCCATAGCCAGCACCGCGCCGCTGTTAATGCCGAGCAGTGAAGGTGAGGCCATTg                                                                                        >  1:38744/1‑140 (MQ=255)                                   cagccagcTCACGCCGCCCCCGCATGCCGCGATGAATGACAGAGAATAGCCTGCAATCGGCGTCGGACTCAGCGCGCTGGTAAGCGCCATAGCCAGCACCCCGCTGCTGTTAATGCCAAGCATTGAAGGAGAGGCCATTg                                                                                        >  2:290227/1‑140 (MQ=255)                                       cagcTCACGCCGCCCCCGCATGCCGCGATGAATGACAGAGAATAGCCTGCAATCGGCGTCGGACTCAGCGCGCTGGTAAGCGCCATAGCCAGCACCGCGCCGCTGTTAATGCCGAGCAGTGAAGGAGAGGCCATTGGGtt                                                                                    <  2:718716/140‑1 (MQ=255)                                       cagcTCACGCCGCCCCCGCATGCCGCGATGAATGACAGAGAATAGCCTGCAATCGGCGTCGGACTCAGCGCGCTGGTAAGCGCCATAGCCAGCACCGCGCCGCTGTTAATGCCGAGCAGTGAAGGAGAGGCCATTGGGtt                                                                                    <  2:420430/140‑1 (MQ=255)                                       cagcTCACGCCGCCCCCGCATGCCGCGATGAATGACAGAGAATAGCCTGCAATCGGCGTCGGACTCAGCGCGCTGGTAAGCGCCATAGCCAGCACCGCGCCGCTGTTAATGCCGAGCAGTGAAGGAGAGGCCATTGGGtt                                                                                    <  2:496636/140‑1 (MQ=255)                                        agcTCACGCCGCCCCCGCATGCCGCGATGAATGACAGAGAATAGCCTGCAATCGGCGTCGGACTCAGCGCGCTGGTAAGCGCCATAGCCAGCACCGCGCCGCTGTTAATGCCGAGCAGTGAAGGAGAGGCCATTGGGTtg                                                                                   >  1:43247/1‑140 (MQ=255)                                          cTCACGCCGCCCCCGCATGCCGCGATGAATGACAGAGAATAGCCTGCAATCGGCGTCGGACTCAGCGCGCTGGTAAGCGCCATAGCCAGCACCGCGCCGCTGTTAATGCCGAGCAGTGAAGGAGAGGCCATTGGGTtgtg                                                                                 >  2:128371/1‑140 (MQ=255)                                                             ccGCGATGAATGACAGCGAATAGCCTGCAATCGGCGGCGGACTCATCGCGCTGGTAAGCGCCATAACCAGCACCGCGCCGCTGTTAATGCCGAGCAGTGAAGGAGAGGCCATTGGGATGTGGGTCAGGGTTTGCagcagc                                                              >  2:690730/1‑140 (MQ=255)                                                                   tgaatgaCAGAGAATAGCCTGCAATCGGCGTCGGACTCAGCGCGCTGGTAAGCGCCATAGCCAGCACCGCGCCGCTGTTAATGCCGAGCAGTGAAGGAGAGGCCATTGGGTTGTGGGTCAGGGTTTGCAGCAGCGTGccc                                                        >  1:499526/1‑140 (MQ=255)                                                                           agagAATAGCCTGCATTCGGCGCCGGACTCAGCGCGCTGGTAAGCGCCATAGCCAGCACCGCGCCGCTGTTAATGCCGAGCAGTGAAGGAGAGGCCATTGGGTTGTGGGTCAGGGTTTGCAGCAGCGTGCCCGCGAgcgc                                                <  1:290227/140‑1 (MQ=255)                                                                                aTAGCCTGCAATCGGCGTCGGACTCAGCGCGCTGGTAAGCGCCATAGCCAGCACCGCGCCGCTGTTAATGCCGAGCAGTGAAGGAGAGGCCATTGGGTTGTGGGTCAGGGTTTGCAGCAGCGTGCCCGCGAGCGCCAGGc                                           <  2:499526/140‑1 (MQ=255)                                                                                aTAGCCTGCAACCGGCGTCGGACTCAGCGCGATGGTAAGCGCCATAGCTAGCACCGCGCCGCTGTTAAGGCCGAGGAGTGAAGGAGAGGCCATTGGGTTGTGGGTCAGGGTTTGCAGCAGCGTGCCCGCGAGCGCCAGGc                                           <  2:270774/140‑1 (MQ=255)                                                                                    ccTGCAATCGGCGTCGGACTCAGCGCGCTGGTAAGCGCCATAGCCAGCACCGCGCCGCTGTTAATGCCGAGCAGTGAAGGAGAGGCCATTGGGTTGTGGGTCAGGGTTTGCAGCAGCGTGCCCGCGAGCGCCAGGCTTgc                                       >  1:644841/1‑140 (MQ=255)                                                                                     cTGCAATCGGCGTCGGACTCAGCGCGCTGGTAAGCGCCATAGCCAGCACCGCGCCGCTGTTAATGCCGAGCAGTGAAGGAGAGGCCATTGGGTTGTGGGTCAGGGTTTGCAGCAGCGTGCCCGCGAGCGCCAGGCTTgcg                                      <  2:576729/140‑1 (MQ=255)                                                                                       gCAATCGGCGTCGGACTCAGCGCGCTGGTAAGCGCCATAGCCAGCACCGCGCCGCTGTTAATGCCGAGCAGTGAAGGAGAGGCCATTGGGTTGTGGGTCAGGTTTTGCAGTAGCATGACCGCGATCGCCAGGCTTGCGcc                                    >  1:306744/1‑140 (MQ=255)                                                                                         aaTCGGCGTCGGACTCAGCGCGCTGGTAAGCGCCATAGCCAGCACCGCGCCGCTGTTAATGCCGAGCAGTGAAGGAGAGGCCATTGGGTTGTGGGTCAGGGTTTGCAGCAGCGTGCCCGCGAGCGCCAGGCTTgcg                                      <  2:149391/136‑1 (MQ=255)                                                                                         aaTCGGCGTCGGACTCAGCGCGCTGGTAAGCGCCATAGCCAGCACCGCGCCGCTGTTAATGCCGAGCAGTGAAGGAGAGGCCATTGGGTTGTGGGTCAGGGTTTGCAGCAGCGTGCCCGCGAGCGCCAGGCTTgcg                                      >  1:149391/1‑136 (MQ=255)                                                                                         aaTCGGCGTCGGACTCAGCGCGCTGGTAAGCGCCATAGCCAGCACCGCGCCGCTGTTAATGCCGAGCAGTGAAGGAGAGGCCATTGGGTTGTGGGTCAGGGTTTGCAGCAGCGTGCCCGCGAGCGCCAGGCTTGCGCTGa                                  >  1:387677/1‑140 (MQ=255)                                                                                          aTCGGCGTCGGACTCAGCGCGCTGGTAAGCGCCATAGCCAGCACCGCGCCGCTGTTAATGCCGAGCAGTGAAGGAGAGGCCATTGTGTTGTGGGTCAGGGTTTGCAGCAGCGTGCCCGCGAGCGCCAGGCTTGCGCCGAt                                 <  2:376497/140‑1 (MQ=255)                                                                                                gTCGGACTCAGCGCGCTGGTAAGCGCCATAGCCAGCACCGCGCCGCTGTTAATGCCGAGCAGTGAAGGAGAGGCCATTGGGTTGTGGGTCAGGGTTTGCAGCAGCGTGCCCGCGAGCGCCAGGCTTGCGCCGATCAGAAc                           <  2:670991/140‑1 (MQ=255)                                                                                                       tCAGCGCGCTGGTAAGCGCCATAGCCAGCACCGCGCCGCTGTTAATGCCGAGCAGTGAAGGAGAGGCCATTGGGTTGTGGGTCAGGGTTTGGAGCAGCGTGCCCTCTACTGCCAGGCATGCGCTGATCAGAACGGCGAcc                    >  2:340405/1‑140 (MQ=255)                                                                                                                  gTAAGCGCCATAGCCAGCACCGCGCCGCTGTTAATGCCGAGCAGTGAAGGAGAGGCCATTGGGTTGTGGGTCAGGGTTTGCAGCAGCGTGCCCGCGAGCGCCAGGCTTGCGCCGATCAGAACGGCGACCAGGCTTCGTgg         <  2:734991/140‑1 (MQ=255)                                                                                                                         ccATAGCCAGCACCGCGCCGCTGTTAATGCCGAGCAGTGAAGGAGAGGCCATTGGGTTGTGGGTCAGGGTTTGCAGCAGCGTGCCCGCGAGCGCCAGGCTTGCGCCGATCAGAACGGCGACCAGGCTTCGTGGCAAACGa  <  1:690730/140‑1 (MQ=255)                                                                                                                                    |                                                                                                                                 TATGACGAAATCCGCCTCCTGCGGTCATGACCAGCAGCCAGCTCACGCCGCCCCCGCATGCCGCGATGAATGACAGAGAATAGCCTGCAATCGGCGTCGGACTCAGCGCGCTGGTAAGCGCCATAGCCAGCGCCGCGCCGCTGTTAATGCCGAGCAGTGAAGGAGAGGCCATTGGGTTGTGGGTCAGGGTTTGCAGCAGCGTGCCCGCGAGCGCCAGGCTTGCGCCGATCAGAACGGCGACCAGGCTTCGTGGCAAACGA  >  NC_000913/4512974‑4513233

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