Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	NC_000913	4,513,886	A→G	A141A (GCT→GCC)	fecB ←	ferric citrate ABC transporter periplasmic binding protein

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	4,513,886	0	A	G	100.0%	175.1 / NA	48	A141A (GCT→GCC)	fecB	ferric citrate ABC transporter periplasmic binding protein
Reads supporting (aligned to +/- strand):  ref base A (0/0);  new base G (26/22);  total (26/22)

CTGTTCGCGTGATGTGCCAAAGGCCACGCGTGTCCCTTTGGGAAGCTGGCTGGCCCACTGCGCCATCCTCTCTTTATGTTGTTCCAGACGTGCCTGCATCTCTCGCTTTTTACCCACCATTTCGCCGATGATAGCCGCAGATTGCAAATTTTCAGCGTAGGTTTCGTTGCGGGACTTAAGCAGCAGTACCGGCGCGATTTGCTGCAAGGCGATGTAAACCCCCGCATGGCGACTGCTGTCGGCAATGATCAGGTCTGGTTTCAGAG  >  NC_000913/4513754‑4514019                                                                                                                                     |                                                                                                                                      cTGTTCGCGTGATGTGCCAAAGGCCACGCGTGTCCCTTTGGGAAGCTGGCTGGCCCACTGCGCCATCCTCTCTTTATGTTGTTCCAGACGTGCCTGCATCTCTCGCTTTTTACCCACCATTTCGCCGATGATGGCCGCAg                                                                                                                                <  1:156266/140‑1 (MQ=255)  tGTTCGCGTGATGTGCCAAAGGCCACGCGTGTCCCTTTGGGAAGCTGGCTGGCCCACTGCGCCATCCTCTCTTTATGTTGTTCCAGACGTGCCTGCATCTCTCGCTTTTTACCCACCATTTCGCCGATGATGGCCGCAGa                                                                                                                               >  1:24589/1‑140 (MQ=255)     tCGCGTGATGTGCCAAAGGCCACGCGTGTCCCTTTGGGAAGCTGGCTGGCCCACTGCGCCATCCTCTCTTTATGTTGTTCCAGACGTGCCTGCATCTCTCGCTTTTTACCCACCATTTCGCCGATGATGGCCGCAGATTg                                                                                                                            >  1:160876/1‑140 (MQ=255)     tCGCGTGATGTGCCAAAGGCCACGCGTGTCCCTTTGGGAAGCTGGCTGGCCCACTGCGCCATCCTCTCTTTATGTTGTTCCAGACGTGCCTGCATCTCTCGCTTTTTACCCACCATTTCGCCGATGATGGCCGCAGATTg                                                                                                                            >  1:921631/1‑140 (MQ=255)      cgcgTGATGTGCCAAAGGCCACGCGTGTCCCTTTGGGAAGCTGGCTGGCCCACTGCGCCATCCTCTCTTTATGTTGTTCCAGACGTGCCTGCATCTCTCGCTTTTTACCCACCATTTCGCCGATGATGGCCGCAGATTGc                                                                                                                           >  2:900740/1‑140 (MQ=255)         gTGATGTGCCAAAGGCCACGCGTGTCCCTTTGGGAAGCTGGCTGGCCCACTGCGCCATCCTCTCTTTATGTTGTTCCAGACGTGCCTGCATCTCTCGCTTTTTACCCACCATTTCGCCGATGATGGCCGc                                                                                                                                  <  1:159634/130‑1 (MQ=255)         gTGATGTGCCAAAGGCCACGCGTGTCCCTTTGGGAAGCTGGCTGGCCCACTGCGCCATCCTCTCTTTATGTTGTTCCAGACGTGCCTGCATCTCTCGCTTTTTACCCACCATTTCGCCGATGATGGCCGc                                                                                                                                  >  2:159634/1‑130 (MQ=255)            aTGTGCCAAAGGCCACGCGTGTCCCTTTGGGAAGCTGGCTGGCCCACTGCGCCATCCTCTCTTTATGTTGTTCCAGACGTGCCTGCATCTCTCGCTTTTTACCCACCATTTCGCCGATGATGGCCGCAGATTGCAAAttt                                                                                                                     >  2:702822/1‑140 (MQ=255)             tgtgCCAAAGGCCACGCGTGTCCCTTTGGGAAGCTGGCTGGCCCACTGCGCCATCCTCTCTTTATGTTGTTCCAGACGTGCCTGCATCTCTCGCTTTTTACCCACCATTTCGCCGATGATGGCCGCAGATTGCAAAtttt                                                                                                                    <  2:921631/140‑1 (MQ=255)             tgtgCCAAAGGCCACGCGTGTCCCTTTGGGAAGCTGGCTGGCCCACTGCGCCATCCTCTCTTTATGTTGTTCCAGACGTGCCTGCATCTCTCGCTTTTTACCCACCATTTCGCCGATGATGGCCGCAGATTGCAAAtttt                                                                                                                    <  2:160876/140‑1 (MQ=255)                         cACGCGTGTCCCTTTGGGAAGCTGGCTGGCCCACTGCGCCATCCTCTCTTTATGTTGTTCCAGACGTGCCTGCATCTCTCGCTTTTTACCCACCATTTCGCCGATGATGGCCGCAGATTGCAAATTTTCAGCGTAGGttt                                                                                                        >  1:136853/1‑140 (MQ=255)                           cgcgTGTCCCTTTGGGAAGCTGGCTGGCCCACTGCGCCATCCTCTCTTTATGTTGTTCCAGACGTGCCTGCATCTCTCGCTTTTTACCCACCATTTCGCCGATGATGGCCGCAGATTGCAAATTTTCAGCGTAGGTTTCg                                                                                                      >  1:354544/1‑140 (MQ=255)                           cgcgTGTCCCTTTGGGAAGCTGGCTGGCCCACTGCGCCATCCTCTCTTTATGTTGTTCCAGACGTGCCTGCATCTCTCGCTTTTTACCCACCATTTCGCCGATGATGGCCGCAGATTGCAAATTTTCAGCGTAGGTTTCg                                                                                                      >  1:576563/1‑140 (MQ=255)                               tgtCCCTTTGGGAAGCTGGCTGGCCCACTGCGCCATCCTCTCTTTATGTTGTTCCAGACGTGCCTGCATCTCTCGCTTTTTACCCACCATTTCGCCGATGATGGCCGCAGATTGCAAATTTTCAGCGTAGGTTTCGTTGc                                                                                                  <  1:837484/140‑1 (MQ=255)                                 tCCCTTTGGGAAGCTGGCTGGCCCACTGCGCCATCCTCTCTTTATGTTGTTCCAGACGTGCCTGCATCTCTCGCTTTTTACCCACCATTTCGCCGATGATGGCCGCAGATTGCAAATTTTCAGCGTAGGTTTCGTTGCgg                                                                                                <  2:89186/140‑1 (MQ=255)                                      ttGGGAAGCTGGCTGGCCCACTGCGCCATCCTCTCTTTATGTTGTTCCAGACGTGCCTGCATCTCTCGCTTTTTACCCACCATTTCGCCGATGATGGCCGCAGATTGCAATTTTTCAGCGTAGGTTTCGTTGCGGGACtt                                                                                           >  2:373723/1‑140 (MQ=255)                                       tGGGAAGCTGGCTGGCCCACTGCGCCATCCTCTCTTTATGTTGTTCCAGACGTGCCTGCATCTCTCGCTTTTTACCCACCATTTCGCCGATGATGGCCGCAGATTGCAAATTTTCAGCGTAGGTTTCGTTGCGGGACTTa                                                                                          <  2:136853/140‑1 (MQ=255)                                              ctggctggCCCACTGCGCCATCCTCTCTTTATGTTGTTCCAGACGTGCCTGCATCTCTCGCTTTTTACCCACCATTTCGCCGATGATGGCCGCAGATTGCAAATTTTCAGCGTAGGTTTCGTTGCGGGACTTAAgcagca                                                                                   >  2:301116/1‑140 (MQ=255)                                                    ggcCCACTGCGCCATCCTCTCTTTATGTTGTTCCAGACGTGCCTGCATCTCTCGCTTTTTACCCACCATTTCGCCGATGATGGCCGCAGATTGCAAATTTTCAGCGTAGGTTTCGTTGCGGGACTTAAGCAGCAGTACCg                                                                             <  2:124235/140‑1 (MQ=255)                                                    ggcCCACTGCGCCATCCTCTCTTTATGTTGTTCCAGACGTGCCTGCATCTCTCGCTTTTTACCCACCATTTCGCCGATGATGGCCGCAGATTGCAAATTTTCAGCGTAGGTTTCGTTGCGGGACTTAAGCAGCAGTACCg                                                                             >  2:771437/1‑140 (MQ=255)                                                       ccACTGCGCCATCCTCTCTTTATGTTGTTCCAGACGTGCCTGCATCTCTCGCTTTTTACCCACCATTTCGCCGATGATGGCCGCAGATTGCAAATTTTCAGCGTAGGTTTCGTTGTGGGACTTAAGCAGCAGTACCGgcg                                                                          <  2:354544/140‑1 (MQ=255)                                                                   ccTCTCTTTATGTTGTTCCAGACGTGCCTGCATCTCTCGCTTTTTACCCACCATTTCGCCGATGATGGCCGCAGATTGCAAATTTTCAGCGTAGGTTTCGTTGCGGGACTTAAGCAGCAGTACCGGCGCGATTTGCTGCa                                                              <  1:771437/140‑1 (MQ=255)                                                                      ctctTTATGTTGTTCCAGACGTGCCTGCATCTCTCGCTTTTTACCCACCATTTCGCCGATGATGGCCGCAGATTGCAAATTTTCAGCGTAGGTTTCGTTGCGGGACTTAAGCAGCAGTACCGGCGCGATTTGCTGCAAgg                                                           <  1:175248/140‑1 (MQ=255)                                                                                tgttCCAGACGTGCCTGCATCTCTCGCTTTTTACCCACCATTTCGCCGATGATGGCCGCAGATTGCAAATTTTCAGCGTAGGTTTCGTTGCGGGACTTAAGCAGCAGTACCGGCGCGATTTGCTGCAAGGCGATGTAAAc                                                 <  1:702822/140‑1 (MQ=255)                                                                                 gttCCAGACGTGCCTGCATCTCTCGCTTTTTACCCACCATTTCTCCGATGATGGCCGCAGATTGCAAATTTTCAGCGTAGGTTTCGTTGCGGGACTTAAGCAGCAGTACCGGCGCGATTTGCTGCAAGGCGATGTaaa                                                  <  2:132133/138‑1 (MQ=255)                                                                                 gttCCAGACGTGCCTGCATCTCTCGCTTTTTACCCACCATTTCTCCGATGATGGCCGCAGATTGCAAATTTTCAGCGTAGGTTTCGTTGCGGGACTTAAGCAGCAGTACCGGCGCGATTTGCTGCAAGGCGATGTaaa                                                  >  1:132133/1‑138 (MQ=255)                                                                                 gttCCAGACGTGCCTGCATCTCTCGCTTTTTACCCACCATTTCGCCGATGATGGCCGCAGATTGCAAATTTTCAGCGTAGGTTTCGTTGCGGGACTTAAGCAGCAGTACCGGCGCGATTTGCTGCAAGGCGATGTa                                                    >  2:391351/1‑136 (MQ=255)                                                                                 gttCCAGACGTGCCTGCATCTCTCGCTTTTTACCCACCATTTCGCCGATGATGGCCGCAGATTGCAAATTTTCAGCGTAGGTTTCGTTGCGGGACTTAAGCAGCAGTACCGGCGCGATTTGCTGCAAGGCGATGTa                                                    <  1:391351/136‑1 (MQ=255)                                                                                 gttCCAGACGTGCCTGCATCTCTCGCTTTTTACCCACCATTTCGCCGATGATGGCCGCAGATTGCAAATTTTCAGCGTAGGTTTCGTTGCGGGACTTAAGCAGCAGTACCGGCGCGATTTGCTGCAAGGCGATGTAAAcc                                                >  2:459407/1‑140 (MQ=255)                                                                                 gttCCAGACGTGCCTGCATCTCTCGCTTTTTACCCACCATTTCGCCGATGATGGCCGCAGATTGCAAATTTTCAGCGTAGGTTTCGTTGCGGGACTTAAGCAGCAGTACCGGCGCGATTTGCTGCAAGGCGATGTAAAcc                                                >  2:103160/1‑140 (MQ=255)                                                                                            gCCTGCATCTCTCGCTTTTTACCCACCATTTCGCCGATGATGGCCGCAGATTGCAAATTTTCAGCGTAGGTTTCGTTGCGGGACTTAAGCAGCAGTACCGGCGCGATTTGCTGCAAGGCGATGTAAACCCCCGCATGGCg                                     >  2:774756/1‑140 (MQ=255)                                                                                            gCCTGCATCTCTCGCTTTTTACCCACCATTTCGCCGATGATGGCCGCAGATTGCAAATTTTCAGCGTAGGTTTCGTTGCGGGACTTAAGCAGCAGTACCGGCGCGATTTGCTGCAAGGCGATGTAAACCCCCGCATGGCg                                     >  1:613994/1‑140 (MQ=255)                                                                                            gCCTGCATCTCTCGCTTTTTACCCACCATTTCGCCGATGATGGCCGCAGATTGCAAATTTTCAGCGTAGGTTTCGTTGCGGGACTTAAGCAGCAGTACCGGCGCGATTTGCTGCAAGGCGATGTAAACCCCCGCATGGCg                                     >  2:315289/1‑140 (MQ=255)                                                                                             ccTGCATCTCTCGCTTTTTACCCACCATTTCGCCGATGATGGCCGCAGATTGCAAATTTTCAGCGTAGGTTTCGTTGCGGGACTTAAGCAGCAGTACCGGCGCGATTTGCTGCAAGGCGATGTAAACCCCCGCATGGCGa                                    >  1:135454/1‑140 (MQ=255)                                                                                                gCATCTCTCGCTTTTTACCCACCATTTCGCCGATGATGGCCGCAGATTGCAAATTTTCAGCGTAGGTTTCGTTGCGGGACTTAAGCAGCAGTACCGGCGCGATTTGCTGCAAGGCGATGTAAACCCCCGCATGGCGActg                                 <  1:103160/140‑1 (MQ=255)                                                                                                    ctctcGCTTTTTACCCACCATTTCGCCGATGATGGCCGCAGATTGCAAATTTTCAGCGTAGGTTTCGTTGCGGGACTTAAGCAGCAGTACCGGCGCGATTTGCTGCAAGGCGATGTAAACCCCCGCAt                                         >  1:776206/1‑128 (MQ=255)                                                                                                    ctctcGCTTTTTACCCACCATTTCGCCGATGATGGCCGCAGATTGCAAATTTTCAGCGTAGGTTTCGTTGCGGGACTTAAGCAGCAGTACCGGCGCGATTTGCTGCAAGGCGATGTAAACCCCCGCAt                                         <  2:776206/128‑1 (MQ=255)                                                                                                    ctctcGCTTTTTACCCACCATTTCGCCGATGATGGCCGCAGATTGCAAATTTTCAGCGTAGGTTTCGTTGCGGGACTTAAGCAGCAGTACCGGCGCGATTTGCTGCAAGGCGATGTAAACCCCCGCATGGCGACTGCTGt                             <  1:301116/140‑1 (MQ=255)                                                                                                    ctctcGCTTTTTACCCACCATTTCGCCGATGATGGCCGCAGATTGCAAATTTTCAGCGTAGGTTTCGTTGCGGGACTTAAGCAGCAGTACCGGCGCGATTTGCTGCAAGGCGATGTAAACCCCCGCATGGCGACTGCTGt                             <  1:315289/140‑1 (MQ=255)                                                                                                      ctcGCTTTTTACCCACCATTTCGCCGATGATGGCCGCAGATTGCAAATTTTCAGCGTAGGTTTCGTTGCGGGACTTAAGCAGCAGTACCGGCGCGATTTGCTGCAAGGCGATGTAAACCCCCGCATGGCGACTGCTGTCg                           >  1:969041/1‑140 (MQ=255)                                                                                                      ctcGCTTTTTACCCACCATTTCGCCGATGATGGCCGCAGATTGCAAATTTTCAGCGTAGGTTTCGTTGCGGGACTTAAGCAGCAGTACCGGCGCGATTTGCTGCAAGGCGATGTAAACCCCCGCATGGCGACTGCTGTCg                           >  1:816033/1‑140 (MQ=255)                                                                                                      ctcGCTTTTTACCCACCATTTCGCCGATGATGGCCGCAGATTGCAAATTTTCAGCGTAGGTTTCGTTGCGGGACTTAAGCAGCAGTACCGGCGCGATTTGCTGCAAGGCGATGTAAACCCCCGCATGGCGACTGCTGTCg                           >  1:67694/1‑140 (MQ=255)                                                                                                         gCTTTTTACCCACCATTTCGCCGATGATGGCCGCAGATTGCAAATTTTCAGCGTAGGTTTCGTTGCGGGACTTAAGCAGCAGTACCGGCGCGATTTGCTGCAAGGCGATGTAAACCCCCGCATGGCGACTGCTGTCGGCa                        <  2:67694/140‑1 (MQ=255)                                                                                                            ttttACCCACCATTTCGCCGATGATGGCCGCAGATTGCAAATTTTCAGCGTAGGTTTCGTTGCGGGACTTAAGCAGCAGTACCGGCGCGATTTGCTGCAAGGCGATGTAAACCCCCGCATGGCGACTGCTGTCGGCAATg                     >  1:342687/1‑140 (MQ=255)                                                                                                                   caccaTTTCGCCGATGATGGCCGCAGATTGCAAATTTTCAGCGTAGGTTTCGTTGCGGGACTTAAGCAGCAGTACCGGCGCGATTTGCTGCAAGGCGATGTAAACCCCCGCATGGCGACTGCTGTCGGCAATGATCAGGt              <  1:774756/140‑1 (MQ=255)                                                                                                                         ttCGCCGATGATGGCCGCAGATTGCAAATTTTCAGCGTAGGTTTCGTTGCGGGACTTAAGCAGCAGTACCGGCGCGATTTGCTGCAAGGCGATGTAAACCCCCGCATGGCGACTGCTGTCGGCa                        <  2:643816/124‑1 (MQ=255)                                                                                                                         ttCGCCGATGATGGCCGCAGATTGCAAATTTTCAGCGTAGGTTTCGTTGCGGGACTTAAGCAGCAGTACCGGCGCGATTTGCTGCAAGGCGATGTAAACCCCCGCATGGCGACTGCTGTCGGCa                        >  1:643816/1‑124 (MQ=255)                                                                                                                               gatgatGGCCGCAGATTGCAAATTTTCAGCGTAGGTTTCGTTGCGGGACTTAAGCAGCAGTACCGGCGCGATTTGCTGCAAGGCGATGTAAACCCCCGCATGGCGACTGCTGTCGGCAATGATCAGGTCTGGTTTCagag  <  2:24589/140‑1 (MQ=255)                                                                                                                                     |                                                                                                                                      CTGTTCGCGTGATGTGCCAAAGGCCACGCGTGTCCCTTTGGGAAGCTGGCTGGCCCACTGCGCCATCCTCTCTTTATGTTGTTCCAGACGTGCCTGCATCTCTCGCTTTTTACCCACCATTTCGCCGATGATAGCCGCAGATTGCAAATTTTCAGCGTAGGTTTCGTTGCGGGACTTAAGCAGCAGTACCGGCGCGATTTGCTGCAAGGCGATGTAAACCCCCGCATGGCGACTGCTGTCGGCAATGATCAGGTCTGGTTTCAGAG  >  NC_000913/4513754‑4514019

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