Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	NC_000913	4,512,341	A→G	V25A (GTG→GCG)	fecD ←	ferric citrate ABC transporter membrane subunit FecD

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	4,512,341	0	A	G	100.0%	151.0 / NA	43	V25A (GTG→GCG)	fecD	ferric citrate ABC transporter membrane subunit FecD
Reads supporting (aligned to +/- strand):  ref base A (0/0);  new base G (24/19);  total (24/19)

ACGCCCGCCACGGCGAGGGCTGCACCGACAAACAGTGCCAGCAGCAAGCGCGGCAGTCGGTACTCCATCAATACATAATAATGCTCGTGTCCGGCCTGCCAGTCGGTCAGCAGCGCGCGCCACGGCACGGGGATCACTCCCATATGGAGTGATAACAGCGCACAGCCCGCCAGGGCAAGGGTGATGAAAATAACCAGCGCAATTTTCATCCTCGCCTCCTCACAAGCCAGACAAAGCAAGGGCTGCCAATCAGCGCCAGCACTGCGCCTGC  >  NC_000913/4512206‑4512476                                                                                                                                        |                                                                                                                                        aCGCCCGCCACGGCGAGGGCTGCACCGACAAACAGTGCCAGCAGCAAGCGCGGCAGTCGGTACTCCATCAATACATAATAATGCTCGTGTCCGGCCTGCCAGTCGGTCAGCAGCGCGCGCCACGGCACGGGGATCGCTcc                                                                                                                                     <  2:188615/140‑1 (MQ=255)    cccgccACGGCGAGGGCTGCACCGACAAACAGTGCCAGCAGCAAGCGCGGCAGTCGGTACTCCATCAATACATAATAATGCTCGTGTCCGGCCTGCCAGTCGGTCAGCAGCGCGCGCCACGGCACGGGGATCGCt                                                                                                                                       >  1:95665/1‑135 (MQ=255)    cccgccACGGCGAGGGCTGCACCGACAAACAGTGCCAGCAGCAAGCGCGGCAGTCGGTACTCCATCAATACATAATAATGCTCGTGTCCGGCCTGCCAGTCGGTCAGCAGCGCGCGCCACGGCACGGGGATCGCTCCCat                                                                                                                                  >  1:189287/1‑140 (MQ=255)    cccgccACGGCGAGGGCTGCACCGACAAACAGTGCCAGCAGCAAGCGCGGCAGTCGGTACTCCATCAATACATAATAATGCTCGTGTCCGGCCTGCCAGTCGGTCAGCAGCGCGCGCCACGGCACGGGGATCGCTCCCat                                                                                                                                  >  2:196426/1‑140 (MQ=255)    cccgccACGGCGAGGGCTGCACCGACAAACAGTGCCAGCAGCAAGCGCGGCAGTCGGTACTCCATCAATACATAATAATGCTCGTGTCCGGCCTGCCAGTCGGTCAGCAGCGCGCGCCACGGCACGGGGATCGCTCCCat                                                                                                                                  >  1:317101/1‑140 (MQ=255)      cgccACGGCGAGGGCTGCACCGACAAACAGTGCCAGCAGCAAGCGCGGCAGTCGGTACTCCATCAATACATAATAATGCTCGTGTCCGGCCTGCCAGTCGGTCAGCAGCGCGCGCCACGGCACGGGGATCGCTCCCatat                                                                                                                                <  1:196426/140‑1 (MQ=255)             gCGAGGGCTGCACCGACAAACAGTGCCAGCAGCAAGCGCGGCAGTCGGTACTCCATCAATACATAATAATGCTCGTGTCCGGCCTGCCAGTCGGTCAGCAGCGCGCGCCACGGCACGGGGATCGCTCCCATATGGAGTGa                                                                                                                         <  2:363509/140‑1 (MQ=255)                     tGCACCGACAAACAGTGCCAGCAGCAAGCGCGGCAGTCGGTACTCCATCAATACATAATAATGCTCGTGTCCGGCCTGCCAGTCGGTCAGCAGCGCGCGCCACGGCACGGGGATCGCTCCCATATGGAGTGATAACAgcg                                                                                                                 >  1:982226/1‑140 (MQ=255)                           gACAAACAGTGCCAGCAGCAAGCGCGGCAGTCGGTACTCCATCAATACATAATAATGCTCGTGTCCGGCCTGCCAGTCGGTCAGCAGCGCGCGCCACGGCACGGGGATCGCTCCCATATGGAGTGATAACAGCGCACAGc                                                                                                           <  2:311527/140‑1 (MQ=255)                               aaCAGTGCCAGCAGCAAGCGCGGCAGTCGGTACTCCATCAATACATAATAATGCTCGTGTCCGGCCTGCCAGTCGGTCAGCAGCGCGCGCCACGGCACGGGGATCGCTCCCATATGGAGTGATAACAGCGCACAGCCCGc                                                                                                       <  2:189287/140‑1 (MQ=255)                               aaCAGTGCCAGCAGCAAGCGCGGCAGTCGGTACTCCATCAATACATAATAATGCTCGTGTCCGGCCTGCCAGTCGGTCAGCAGCGCGCGCCACGGCACGGGGATCGCTCCCATATGGAGTGATAACAGCGCACAGCCCGc                                                                                                       <  1:728364/140‑1 (MQ=255)                                 cAGTGCCAGCAGCAAGCGCGGCAGTCGGTACTCCATCAATACATAATAATGCTCGTGTCCGGCCTGCCAGTCGGTCAGCAGCGCGCGCCACGGCACGGGGATCGCTCCCATATGGAGTGATAACAGCGCACAGCCCGCCa                                                                                                     >  2:93588/1‑140 (MQ=255)                                    tGCCAGCAGCAAGCGCGGCAGTCGGTACTCCATCAATACATAATAATGCTCGTGTCCGGCCTGCCAGTCGGTCAGCAGCGCGCGCCACGGCACGGGGATCGCTCCCATATGGAGTGATAACAGCGCACAGCCCGc                                                                                                       >  1:856880/1‑135 (MQ=255)                                    tGCCAGCAGCAAGCGCGGCAGTCGGTACTCCATCAATACATAATAATGCTCGTGTCCGGCCTGCCAGTCGGTCAGCAGCGCGCGCCACGGCACGGGGATCGCTCCCATATGGAGTGATAACAGCGCACAGCCCGc                                                                                                       <  2:856880/135‑1 (MQ=255)                                      ccAGCAGCAAGCGCGGCAGTCGGTACTCCATCAATACATAATAATGCTCGTGTCCGGCCTGCCAGTCGGTCAGCAGCGCGCGCCACGGCACGGGGATCGCTCCCATATGGAGTGATAACAGCGCACAGCCCGCCAGGGCa                                                                                                >  1:201414/1‑140 (MQ=255)                                      ccAGCAGCAAGCGCGGCAGTCGGTACTCCATCAATACATAATAATGCTCGTGTCCGGCCTGCCAGTCGGTCAGCAGCGCGCGCCACGGCACGGGGATCGCTCCCATATGGAGTGATAACAGCGCACAGCCCGCCAGGGCa                                                                                                >  2:82577/1‑140 (MQ=255)                                             caAGCGCGGCAGTCGGTACTCCATCAATACATAATAATGCTCGTGTCCGGCCTGCCAGTCGGTCAGCAGCGCGCGCCACGGCACGGGGATCGCTCCCATATGGAGTGATAACAGCGCACAGCCCGCCAGGGCAAGGGtga                                                                                         <  2:317101/140‑1 (MQ=255)                                             caAGCGCGGCAGTCGGTACTCCATCAATACATAATAATGCTCGTGTCCGGCCTGCCAGTCGGTCAGCAGCGCGCGCCACGGCACGGGGATCGCTCCCATATGGAGTGATAACAGCGCACAGCCCGCCAGGGCAAGGGtga                                                                                         <  2:141468/140‑1 (MQ=255)                                              aaGCGCGGCAGTCGGTACTCCATCAATACATAATAATACTCGTGTCCGGCCTGCCAGTCGGTCAGCAGCGCGCGCCACGGCACGGGGATCGCTCCCATATGGAGTGATAACAGCGCACAGCCCGCCAGGGCAAGGGtga                                                                                         <  1:204547/139‑1 (MQ=255)                                              aaGCGCGGCAGTCGGTACTCCATCAATACATAATAATACTCGTGTCCGGCCTGCCAGTCGGTCAGCAGCGCGCGCCACGGCACGGGGATCGCTCCCATATGGAGTGATAACAGCGCACAGCCCGCCAGGGCAAGGGtga                                                                                         >  2:204547/1‑139 (MQ=255)                                                                             gataatGCTCGTGTCCGGCCTGCCAGTCGGTCAGCAGCGCGCGCCACGGCACGGGGATCGCTCCCATATGGAGTGATAACAGCGCACAGCCCGCCAGGGCAAGGGTGATGAAAATAACCAGCGCAATTTTCATCCTCGcc                                                         >  1:267802/2‑140 (MQ=255)                                                                                  tGCTCGTGTCCGGCCTGCCAGTCGGTCAGCAGCGCGCGCCACGGCACGGGGATCGCTCCCATATGGAGTGATAACAGCGTACAGCCCGCCAGGGCAAGGGTGATGAAAATAACCAGCGCAATTTTCATCCTCGCctcctc                                                    >  2:894459/1‑140 (MQ=255)                                                                                  tGCTCGTGTCCGGCCTGCCAGTCGGTCAGCAGCGCGCGCCACGGCACGGGGATCGCTCCCATATGGAGTGATAACAGCGCACAGCCCGCCAGGGCAAGGGTGATGAAAATAACCAGCGCAATTTTCATCCt                                                             <  2:628508/131‑1 (MQ=255)                                                                                  tGCTCGTGTCCGGCCTGCCAGTCGGTCAGCAGCGCGCGCCACGGCACGGGGATCGCTCCCATATGGAGTGATAACAGCGCACAGCCCGCCAGGGCAAGGGTGATGAAAATAACCAGCGCAATTTTCATCCt                                                             >  1:628508/1‑131 (MQ=255)                                                                                  tGCTCGTGTCCGGCCTGCCAGTCGGTCAGCAGCGCGCGCCACGGCACGGGGATCGCTCCCATATGGAGTGATAACAGCGCACAGCCCGCCAGGGCAAGGGTGATGAAAATAACCAGCGCAATTTTCATCCTCGCctcctc                                                    >  2:802625/1‑140 (MQ=255)                                                                                  tGCTCGTGTCCGGCCTGCCAGTCGGTCAGCAGCGCGCGCCACGGCACGGGGATCGCTCCCATATGGAGTGATAACAGCGCACAGCCCGCCAGGGCAAGGGTGATGAAAATAACCAGCGCAATTTTCATCCTCAcct                                                        <  1:152157/136‑1 (MQ=255)                                                                                  tGCTCGTGTCCGGCCTGCCAGTCGGTCAGCAGCGCGCGCCACGGCACGGGGATCGCTCCCATATGGAGTGATAACAGCGCACAGCCCGCCAGGGCAAGGGTGATGAAAATAACCAGCGCAATTTTCATCCTCAcct                                                        >  2:152157/1‑136 (MQ=255)                                                                                   gCTCGTGTCCGGCCTGCCAGTCGGTCAGCAGCGCGCGCCACGGCACGGGGATCGCTCCCATATGGAGTGATAACAGCGCACAGCCCGCCAGGGCAAGGGTGATGAAAATAACCAGCGCAATTTTCATCCTCGCCTCCTca                                                   <  2:30609/140‑1 (MQ=255)                                                                                    cTCGTGTCCGGCCTGCCAGTCGGTCAGCAGCGCGCGCCACGGCACGGGGATCGCTCCCATATGGAGTGATAACAGCGCACAGCCCGCCAGGGCAAGGGTGATGAAAATAACCAGCGCAATTTTCATCCTCGCCTCCTcac                                                  >  2:749924/1‑140 (MQ=255)                                                                                    cTCGTGTCCGGCCTGCCAGTCGGTCAGCAGCGCGCGCCACGGCACGGGGATCGCTCCCATATGGAGTGATAACAGCGCACAGCCCGCCAGGGCAAGGGTGATGAAAATAACCAGCGCAATTTTCATCCTCGCCTCCTcac                                                  >  2:461704/1‑140 (MQ=255)                                                                                    cTCGTGTCCGGCCTGCCAGTCGGTCAGCAGCGCGCGCCACGGCACGGGGATCGCTCCCATATGGAGTGATAACAGCGCACAGCCCGCCAGGGCAAGGGTGATGAAAATAACCAGCGCAATTTTCATCCTCGCCTCCTcac                                                  >  2:566205/1‑140 (MQ=255)                                                                                         gtCCGGCCTGCCAGTCGGTCAGCAGCGCGCGCCACGGCACGGGGATCGCTCCCATATGGAGTGATAACAGCGCACAGCCCGCCAGGGCAAGGGTGATGAAAATAACCAGCGCAATTTTCATCCTCGcct                                                        <  1:28803/129‑1 (MQ=255)                                                                                         gtCCGGCCTGCCAGTCGGTCAGCAGCGCGCGCCACGGCACGGGGATCGCTCCCATATGGAGTGATAACAGCGCACAGCCCGCCAGGGCAAGGGTGATGAAAATAACCAGCGCAATTTTCATCCTCGcct                                                        >  2:28803/1‑129 (MQ=255)                                                                                         gtCCGGCCTGCCAGTCGGTCAGCAGCGCGCGCCACGGCACGGGGATCGCTCCCATATGGAGTGATAACAGCGCACAGCCCGCCAGGGCAAGGGTGATGAAAATAACCAGCGCAATTTTCATCCTCGCCTCCTCACAAGcc                                             >  2:266580/1‑140 (MQ=255)                                                                                            cGGCCTGCCAGTCGGTCAGCAGCGCGCGCCACGGCACGGGGATCGCTCCCATATGGAGTGATAACAGCGCACAGCCCGCCAGGGCAAGGGTGATGAAAATAACCAGCGCAATTTTCATCCTCGCCTCCTCACAAGCCAGa                                          <  1:93588/140‑1 (MQ=255)                                                                                              gCCTGCCAGTCGGTCAGCAGCGCGCGCCACGGCACGGGGATCGCTCCCATATGGAGTGATAACAGCGCACAGCCCGCCAGGGCAAGGGTGATGAAAATAACCAGCGCAATTTTCATCCTCGCCTCCTCACAAGCCAGACa                                        >  1:806771/1‑140 (MQ=255)                                                                                                     aGTCGGTCAGCAGCGCGCGCCACGGCACGGGGATCGCTCCCATATGGAGTGATAACAGCGCACAGCCCGCCAGGGCAAGGGTGATGAAAATAACCAGCGCAATTTTCATCCTCGCCTCCTCACAAGCCAGACAAAGCAAg                                 <  2:267802/140‑1 (MQ=255)                                                                                                     aGTCGGTCAGCAGCGCGCGCCACGGCACGGGGATCGCTCCCATATGGAGTGATAACAGCGCACAGCCCGCCAGGGCAAGGGTGATGAAAATAACCAGCGCAATTTTCATCCTCGCCTCCTCACAAGCCAGACAAAGCAAg                                 <  2:853069/140‑1 (MQ=255)                                                                                                        cGGTCAGCAGCGCGCGCCACGGCACGGGGATCGCTCCCATATGGAGTGATAACAGCGCACAGCCCGCCAGGGCAAGGGTGATGAAAATAACCAGCGCAATTTTCATCCTCGCCTCCTCACAAGCCAGACAAAGCAAGGGc                              >  1:939837/1‑140 (MQ=255)                                                                                                                agcGCGCGCCACGGCACGGGGATCGCTCCCATATGGAGTGATAACAGCGCACAGCCCGCCAGGGCAAGGGTGATGAAAATAACCAGCGCAATTTTCATCCTCGCCTCCTCACAAGCCAGACAAAGCAAGGGCTGCCAATc                      <  1:82577/140‑1 (MQ=255)                                                                                                                  cgcgcgcCACGGCACGGGGATCGCTCCCATATGGAGTGATAACAGCGCACAGCCCGCCAGGGCAAGGGTGATGAAAATAACCAGCGCAATTTTCATCCTCGCCTCCTCACAAGCCAGACAAAGCAAGGGCTGCCAATCAg                    >  1:420620/1‑140 (MQ=255)                                                                                                                    cgcgcCACGGCACGGGGATCGCTCCCATATGGAGTGATAACAGCGCACAGCCCGCCAGGGCAAGGGTGATGAAAATAACCAGCGCAATTTTCATCCTCGCCTCCTCACAAGCCAGACAAAGCAAGGGCTg                            <  2:297561/130‑1 (MQ=255)                                                                                                                    cgcgcCACGGCACGGGGATCGCTCCCATATGGAGTGATAACAGCGCACAGCCCGCCAGGGCAAGGGTGATGAAAATAACCAGCGCAATTTTCATCCTCGCCTCCTCACAAGCCAGACAAAGCAAGGGCTg                            >  1:297561/1‑130 (MQ=255)                                                                                                                                    gATCGCTCCCATATGGAGTGATAACAGCGCACAGCCCGCCAGGGCAAGGGTGATGAAAATAACCAGCGCAATTTTCATCCTCGCCTCCTCACAAGCCAGACAAAGCAAGGGCTGCCAATCAGCGCCAGCACTGCGCCTGc  <  2:982226/140‑1 (MQ=255)                                                                                                                                        |                                                                                                                                        ACGCCCGCCACGGCGAGGGCTGCACCGACAAACAGTGCCAGCAGCAAGCGCGGCAGTCGGTACTCCATCAATACATAATAATGCTCGTGTCCGGCCTGCCAGTCGGTCAGCAGCGCGCGCCACGGCACGGGGATCACTCCCATATGGAGTGATAACAGCGCACAGCCCGCCAGGGCAAGGGTGATGAAAATAACCAGCGCAATTTTCATCCTCGCCTCCTCACAAGCCAGACAAAGCAAGGGCTGCCAATCAGCGCCAGCACTGCGCCTGC  >  NC_000913/4512206‑4512476

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