Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	NC_000913	3,216,464	G→A	R10C (CGC→TGC)	yqjH ←	putative siderophore interacting protein

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	3,216,464	0	G	A	100.0%	142.9 / NA	41	R10C (CGC→TGC)	yqjH	putative siderophore interacting protein
Reads supporting (aligned to +/- strand):  ref base G (0/0);  new base A (16/25);  total (16/25)

GCTGTGATCGTCAAAGCCACGCGATGTAAAACCGTCCAGCGCCTCGCCGCCGAGGACAATGCGCTGAAAACCGGCGCTGATGCGCTCAACGCGTAACACAGTCAGTTCACGGAAGCGCAGATCATTGCGAACGCGCTGCGGGTAGCGGGGGGTGTTATTCATTGTTATCGCCTTCGTGATGGTAATCAGATATATCTAAATAAAACTCGCAAATGATAATGATTGTTAATCATGATAAATGCAAGCGATTTGTAGAACTGATA  >  NC_000913/3216330‑3216592                                                                                                                                       |                                                                                                                                 gCTGTGATCGTCAAAGCCACGCGATGTAAAACCGTCCAGCGCCTCGCCGCCGAGGACAATGCGCTGAAAACCGGCGCTGATGCGCTCAACGCGTAACACAGTCAGTTCACGGAAGCGCAGATCATTGCGAACGCACTGc                                                                                                                              >  2:682630/1‑139 (MQ=255)          gTCAAAGCCACGCGATGTAAAACCGTCCAGCGCCTCGCCGCCGAGGACAATGCGCTGAAAACCGGCGCTGATGCGCTCAACGCGTAACACAGTCAGTTCACGGAAGCGCAGATCATTGCGAACGCACTGCgg                                                                                                                            >  1:100954/1‑132 (MQ=255)          gTCAAAGCCACGCGATGTAAAACCGTCCAGCGCCTCGCCGCCGAGGACAATGCGCTGAAAACCGGCGCTGATGCGCTCAACGCGTAACACAGTCAGTTCACGGAAGCGCAGATCATTGCGAACGCACTGCgg                                                                                                                            <  2:100954/132‑1 (MQ=255)           tCAAAGCCACGCGATGTAAAACCGTCCAGCGCCTCGCCGCCGAGGACAATGCGCTGAAAACCGGCGCTGATGCGCTCAACGCGTAACACAGTCAGTTCACGGAAGCGCAGATCATTGCGAACGCACTGCGGGTAGCggg                                                                                                                    <  1:124619/139‑1 (MQ=255)               aGCCACGCGATGTAAAACCGTCCAGCGCCTCGCCGCCGAGGACAATGCGCTGAAAACCGGCGCTGATGCGCTCAACGCGTAACACAGTCAGTTCACGGAAGCGCAGATCATTGCGAACGCACTGCGGGTAGCGGGGGgt                                                                                                                <  2:236421/139‑1 (MQ=255)                    cgcgATGTAAAACCGTCCAGCGCCTCGCCGCCGAGGACAATGCGCTGAAAACCGGCGCTGATGCGCTCAACGCGTAACACAGTCAGTTCACGGAAGCGCAGATCATTGCGAACGCACTGCGGGTAGCGGGGGGTGTTAt                                                                                                           >  2:609923/1‑139 (MQ=255)                       gATGTAAAACCGTCCAGCGCCTCGCCGCCGAGGACAATGCGCTGAAAACCGGCGCTGATGCGCTCAACGCGTAACACAGTCAGTTCACGGAAGCGCAGATCATTGCGAACGCACTGCGGGTAGCGGGGGGTGTTATTCa                                                                                                        <  1:533441/139‑1 (MQ=255)                       gATGTAAAACCGTCCAGCGCCTCGCCGCCGAGGACAATGCGCTGAAAACCGGCGCTGATGCGCTCAACGCGTAACACAGTCAGTTCACGGAAGCGCAGATCATTGCGAACGCACTGCGGGTAGCGGGGGGTGTTATTCa                                                                                                        <  1:682630/139‑1 (MQ=255)                          gTAAAACCGTCCAGCGCCTCGCCGCCGAGGACAATGCGCTGAAAACCGGCGCTGATGCGCTCAACGCGTAACACAGTCAGTTCACGGAAGCGCAGATCATTGCGAACGCACTGCGGGTAGCGGGGGGTGTTATTCa                                                                                                        <  2:618236/136‑1 (MQ=255)                          gTAAAACCGTCCAGCGCCTCGCCGCCGAGGACAATGCGCTGAAAACCGGCGCTGATGCGCTCAACGCGTAACACAGTCAGTTCACGGAAGCGCAGATCATTGCGAACGCACTGCGGGTAGCGGGGGGTGTTATTCa                                                                                                        >  1:618236/1‑136 (MQ=255)                              aaCCGTCCAGCGCCTCGCCGCCGAGGACAATGCGCTGAAAACCGGCGCTGATGCGCTCAACGCGTAACACAGTCAGTTCACGGAAGCGCAGATCATTGCGAACGCACTGCGGGTAGCGGGGGGTGTTATTCATTGTTAt                                                                                                 <  1:191683/139‑1 (MQ=255)                              aaCCGTCCAGCGCCTCGCCGCCGAGGACAATGCGCTGAAAACCGGCGCTGATGCGCTCAACGCGTAACACAGTCAGTTCACGGAAGCGCAGATCATTGCGAACGCACTGCGGGTAGCGGGGGGTGTTATTCATTGTTAt                                                                                                 <  2:30895/139‑1 (MQ=255)                                   tCCAGCGCCTCGCCGCCGAGGACAATGCGCTGAAAACCGGCGCTGATGCGCTCAACGCGTAACACAGTCAGTTCACGGAAGCGCAGATCATTGCGAACGCACTGc                                                                                                                              >  1:368464/1‑105 (MQ=255)                                   tCCAGCGCCTCGCCGCCGAGGACAATGCGCTGAAAACCGGCGCTGATGCGCTCAACGCGTAACACAGTCAGTTCACGGAAGCGCAGATCATTGCGAACGCACTGc                                                                                                                              <  2:368464/105‑1 (MQ=255)                                   tCCAGCGCCTCGCCGCCGAGGACAATGCGCTGAAAACCGGCGCTGATGCGCTCAACGCGTAACACAGTCAGTTCACGGAAGCGCAGATCATTGCGAACGCACTGCGGGTAGCggg                                                                                                                    >  1:452444/1‑115 (MQ=255)                                   tCCAGCGCCTCGCCGCCGAGGACAATGCGCTGAAAACCGGCGCTGATGCGCTCAACGCGTAACACAGTCAGTTCACGGAAGCGCAGATCATTGCGAACGCACTGCGGGTAGCggg                                                                                                                    <  2:452444/115‑1 (MQ=255)                                   tCCAGCGCCTCGCCGCCGAGGACAATGCGCTGAAAACCGGCGCTGATGCGCTCAACGCGTAACACAGTCAGTTCACGGAAGCGCAGATCATTGCGAACGCACTGCGGGTAGCGGGGGGTGTTATTCATTGTTATCGCCt                                                                                            >  1:274924/1‑139 (MQ=255)                                   tCCAGCGCCTCGCCGCCGAGGACAATGCGCTGAAAACCGGCGCTGATGCGCTCAACGCGTAACACAGTCAGTTCACGGAAGCGCAGATCATTGCGAACGCACTGCGGGTAGCGGGGGGTGTTATTCATTGTTATCGCCt                                                                                            >  2:33962/1‑139 (MQ=255)                                   tCCAGCGCCTCGCCGCCGAGGACAATGCGCTGAAAACCGGCGCTGATGCGCTCAACGCGTAACACAGTCAGTTCACGGAAGCGCAGATCATTGCGAACGCACTGCGGGTAGCGGGGGGTGTTATTCATTGTTATCGCCt                                                                                            >  1:740267/1‑139 (MQ=255)                                          ccTCGCCGCCGAGGACAATGCGCTGAAAACCGGCGCTGATGCGCTCAACGCGTAACACAGTCAGTTCACGGAAGCGCAGATCATTGCGAACGCACTGCGGGTAGCGGGGGGTGTTATTCATTGTTATCGCCTTCGTGAt                                                                                     <  1:609923/139‑1 (MQ=255)                                           cTCGCCGCCGAGGACAATGCGCTGAAAACCGGCGCTGATGCGCTCAACGCGTAACACAGTCAGTTCACGGAAGCGCAGATCATTGCGAACGCACTGCGGGTAGCGGGGGGTGTTATTCATTGTTATCGCCTTCGTGATg                                                                                    <  2:274924/139‑1 (MQ=255)                                           cTCGCCGCCGAGGACAATGCGCTGAAAACCGGCGCTGATGCGCTCAACGCGTAACACAGTCAGTTCACGGAAGCGCAGATCATTGCGAACGCACTGCGGGTAGCGGGGGGTGTTATTCATTGTTATCGCCTTCGTGATg                                                                                    <  1:163197/139‑1 (MQ=255)                                                    gAGGACAATGCGCTGAAAACCGGCGCTGATGCGCTCAACGCGTAACACAGTCAGTTCACGGAAGCGCAGATCATTGCGAACGCACTGCGGGTAGCGGGGGGTGTTATTCATTGTTATCGCCTTCGTGATGGTAATCAGa                                                                           <  1:749701/139‑1 (MQ=255)                                                         cAATGCGCTGAAAACCGGCGCTGATGCGCTCAACGCGTAACACAGTCAGTTCACGGAAGCGCAGATCATTGCGAACGCACTGCGGGTAGCGGGGGGTGTTATTCATTGTTATCGCCTTCGTGATGGTAATCAGatatat                                                                      <  2:251535/139‑1 (MQ=255)                                                                 tGAAAACCGGCGCTGATGCGCTCAACGCGTAACACAGTCAGTTCACGGAAGCGCAGATCATTGCGAACGCACTGCGGGTAGCGGGGGGTGTTATTCATTGTTATCGCCTTCGTGATGGTAATCAGATATATCtaaataa                                                              >  1:66953/1‑139 (MQ=255)                                                                   aaaaCCGGCGCTGATGCGCTCAACGCGTAACACAGTCAGTTCACGGAAGCGCAGATCATTGCGAACGCACTGCGGGTAGCGGGGGGTGTTATTCATTGTTATCGCCTTCGTGATGGTAATCAGATATATCTAAATaaaa                                                            >  1:297173/1‑139 (MQ=255)                                                                        cGGCGCTGATGCGCTCAACGCGTAACACAGTCAGTTCACGGAAGCGCAGATCATTGCGAACGCACTGCGGGTAGCGGGGGGTGTTATTCATTGTTATCGCCTTCGTGATGGTAATCAGATATATCTAAATAAAACTCGc                                                       >  1:760038/1‑139 (MQ=255)                                                                               gATGCGCTCAACGCGTAACACAGTCAGTTCACGGAAGCGCAGATCATTGCGAACGCACTGCGGGTAGCGGGGGGTGTTATTCATTGTTATCGCCTTCGTGATGGTAATCAGATATATCTAAATAAAACTCGCAaatgat                                                <  2:297173/139‑1 (MQ=255)                                                                                  gcgcTCAACGCGTAACACAGTCAGTTCACGGAAGCGCAGATCATTGCGAACGCACTGCGGGTAGCGGGGGGTGTTATTCATTGTTATCGCCTTCGTGATGGTAATCAGATATATCTAAATAAAACTCGCAaatgataat                                             <  1:411472/139‑1 (MQ=255)                                                                                                   cagtcagtTCACGGAAGCGCAGATCATTGCGAACGCACTGCGGGTAGCGGGGGGTGTTATTCATTGTTATCGCCTTCGTGATGGTAATCAGATATATCTAAATAAAACTCGCAAATGATAATGATTGTTAATCATGATa                            <  1:96308/139‑1 (MQ=255)                                                                                                    agtcagtTCACGGAAGCGCAGATCATTGCGAACGCACTGCGGGTAGCGGGGGGTGTTATTCATTGTTATCGCCTTCGTGATGGTAATCAGatata                                                                       >  2:53440/1‑95 (MQ=255)                                                                                                    agtcagtTCACGGAAGCGCAGATCATTGCGAACGCACTGCGGGTAGCGGGGGGTGTTATTCATTGTTATCGCCTTCGTGATGGTAATCAGatata                                                                       <  1:53440/95‑1 (MQ=255)                                                                                                             aCGGAAGCGCAGATCATTGCGAACGCACTGCGGGTAGCGGGGGGTGTTATTCATTGTTATCGCCTTCGTGATGGTAATCAGATATATCTAAATAAAACTCGCAAATGATAATGATTGTTAATCATGATAAATGCAAGCg                  <  1:84093/139‑1 (MQ=255)                                                                                                               ggAAGCGCAGATCATTGCGAACGCACTGCGGGTAGCGGGGGGTGTTATTCATTGTTATCGCCTTCGTGATGGTAATCAGATATATCTAAATAAAACTCGCAAATGATAATGATTGTTAATCATGATAAATGCAAGCGAt                <  2:760038/139‑1 (MQ=255)                                                                                                                    cgcAGATCATTGCGAACGCACTGCGGGTAGCGGGGGGTGTTATTCATTGTTATCGCCTTCGTGATg                                                                                    >  2:355732/1‑66 (MQ=255)                                                                                                                    cgcAGATCATTGCGAACGCACTGCGGGTAGCGGGGGGTGTTATTCATTGTTATCGCCTTCGTGATg                                                                                    <  1:355732/66‑1 (MQ=255)                                                                                                                    cgcAGATCATTGCGAACGCACTGCGGGTAGCGGGGGGTGTTATTCATTGTTATCGCCTTCGTGATGGTAATCAGATATATCTAAATAAAACTCGCAAATGATAATGATTGTTAATCATGATAAATGCAAGCGATTTGTa           >  1:74285/1‑139 (MQ=255)                                                                                                                            aTTGCGAACGCACTGCGGGTAGCGGGGGGTGTTATTCATTGTTATCGCCTTCGTGATGGTAATCAGATATATCTAAATAAAACTCGCAAATGATAATGATTGTTAATCATGATAAATGCAAGCGATTTGTAGAACTGat   <  2:66953/139‑1 (MQ=255)                                                                                                                             ttGCGAACGCACTGCGGGTAGCGGGGGGTGTTATTCATTGTTATCGCCTTCGTGATGGTAATCAGATATATCTAAATAAAACTCGCAAATGATAATGATTGTTAATCATGATAAATGCAAGCGATTTGTAGAACTGata  <  2:176543/139‑1 (MQ=255)                                                                                                                              tGCGAACGCACTGCGGGTAGCGGGGGGTGTTATTCATTGTTATCGCCTTCg                                                                                         >  1:534779/1‑51 (MQ=255)                                                                                                                              tGCGAACGCACTGCGGGTAGCGGGGGGTGTTATTCATTGTTATCGCCTTCg                                                                                         <  2:534779/51‑1 (MQ=255)                                                                                                                                       |                                                                                                                                 GCTGTGATCGTCAAAGCCACGCGATGTAAAACCGTCCAGCGCCTCGCCGCCGAGGACAATGCGCTGAAAACCGGCGCTGATGCGCTCAACGCGTAACACAGTCAGTTCACGGAAGCGCAGATCATTGCGAACGCGCTGCGGGTAGCGGGGGGTGTTATTCATTGTTATCGCCTTCGTGATGGTAATCAGATATATCTAAATAAAACTCGCAAATGATAATGATTGTTAATCATGATAAATGCAAGCGATTTGTAGAACTGATA  >  NC_000913/3216330‑3216592

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