Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	NC_000913	4,513,286	C→T	A42T (GCC→ACC)	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,286	0	C	T	100.0%	119.2 / NA	37	A42T (GCC→ACC)	fecC	ferric citrate ABC transporter membrane subunit FecC
Reads supporting (aligned to +/- strand):  ref base C (0/0);  new base T (19/18);  total (19/18)

TTGTGGGTCAGGGTTTGCAGCAGCGTGCCCGCGAGCGCCAGGCTTGCGCCGATCAGAACGGCGACCAGGCTTCGTGGCAAACGAAGGTTTTGCACCAGCGCTTCTGGTAGCGTTGGCGTGTGTCCAGGCAGCAGGGCGCGGGTTGCATCTGCTCCGGAAACAGGAATGGCCGAGTAGCAAAACAGACTCAGCCAGAAAATAATAATAAGTGCTGCAACGGGAAGCCCCCACAGCAGCACCGGGTGTTTTATCGCGGTCAT  >  NC_000913/4513150‑4513409                                                                                                                                         |                                                                                                                            ttGTGGGTCAGGGTTTGCAGCAGCGTGCCCGCGAGCGCCAGGCTTGCGCCGATCAGAACGGCGACCAGGCTTCGTGGCAAACGAAGGTTTTGCACCAGCGCTTCTGGTAGCGTTGGCGTGTGTCCAGGCAGCAGGGtgcg                                                                                                                          <  2:391817/140‑1 (MQ=255)  tgtgGGTCAGGGTTTGCAGCAGCGTGCCCGCGAGCGCCAGGCTTGCGCCGATCAGAACGGCGACCAGGCTTCGTGGCAAACGAAGGTTTTGCACCAGCGCTTCTGGTAGCGTTGGCGTGTGTCCAGGCAGCAGGGTGCgg                                                                                                                         >  2:186178/1‑140 (MQ=255)         cAGGGTTTGCAGCAGCGTGCCCGCGAGCGCCAGGCTTGCGCCGATCAGAACGGCGACCAGGCTTCGTGGCAAACGAAGGTTTTGCACCAGCGCTTCTGGTAGCGTTGGCGTGTGTCCAGGCAGCAGGGTGCGGGTTGCAt                                                                                                                  >  2:491832/1‑140 (MQ=255)         cAGGGTTTGCAGCAGCGTGCCCGCGAGCGCCAGGCTTGCGCCGATCAGAACGGCGACCAGGCTTCGTGGCAAACGAAGGTTTTGCACCAGCGCTTCTGGTAGCGTTGGCGTGTGTCCAGGCAGCAGGGTGCGGGTTGCAt                                                                                                                  >  1:734682/1‑140 (MQ=255)          aGGGTTTGCAGCAGCGTGCCCGCGAGCGCCAGGCTTGCGCCGATCAGAACGGCGACCAGGCTTCGTGGCAAACGAAGGTTTTGCACCAGCGCTTCTGGTAGCGTTGGCGTGTGTCCAGGCAGCAGGGTGCGGGTTGCATc                                                                                                                 >  2:568814/1‑140 (MQ=255)              tttGCAGCAGCGTGCCCGCGAGCGCCAGGCTTGCGCCGATCAGAACGGCGACCAGGCTTCGTGGCAAACGAAGGTTTTGCACCAGCGCTTCTGGTAGCGTTGGCGTGTGTCCAGGCAGCAGGGTGCGGGTTGCATCTGCt                                                                                                             <  2:356543/140‑1 (MQ=255)                   agcagcGTGCCCGCGAGCGCCAGGCTTGCGCCGATCAGAACGGCGACCAGGCTTCGTGGCAAACGAAGGTTTTGCACCAGCGCTTCTGGTAGCGTTGGCGTGTGTCCAGGCAGCAGGGTGCGGGTTGCATCTGCTCCGGa                                                                                                        >  2:719368/1‑140 (MQ=255)                      agcGTGCCCGCGAGCGCCAGGCTTGCGCCGATCAGAACGGCGACCAGGCTTCGTGGCAAACGAAGGTTTTGCACCAGCGCTTCTGGTAGCGTTGGCGTGTGTCCAGGCAGCAGGGTGCGGGTTGCATCTGCTCCGGAAAc                                                                                                     >  2:65237/1‑140 (MQ=255)                          tGCCCGCGAGCGCCAGGCTTGCGCCGATCAGAACGGCGACCAGGCTTCGTGGCAAACGAAGGTTTTGCACCAGCGCTTCTGGTAGCGTTGGCGTGTGTCCAGGCAGCAGGGTGCGGGTTGCATCTGCTCCGGAAACAGGa                                                                                                 >  2:401461/1‑140 (MQ=255)                                 gAGCGCCAGGCTTGCGCCGATCAGAACGGCGACCAGGCTTCGTGGCAAACGAAGGTTTTGCACCAGCGCTTCTGGTAGCGTTGGCGTGTGTCCAGGCAGCAGGGTGCGGGTTGCATCTGCTCCGGAAACAGGAATGGCCg                                                                                          <  2:991925/140‑1 (MQ=255)                                           cTTGCGCCGATCAGAACGGCGACCAGGCTTCGTGGCAAACGAAGGTTTTGCACCAGCGCTTCTGGTAGCGTTGGCGTGTGTCCAGGCAGCAGGGTGCGGGTTGCATCTGCTCCGGAAACAGGAATGGCCGAGTAGCaaaa                                                                                <  2:986824/140‑1 (MQ=255)                                               cgcCGATCAGAACGGCGACCAGGCTTCGTGGCAAACGAAGGTTTTGCACCAGCGCTTCTGGTAGCGTTGGCGTGTGTCCAGGCAGCAGGGTGCGGGTTGCATCTGCTCCGGAAACAGGAATGGCCGAGTAGCAAAACAGa                                                                            <  1:931887/140‑1 (MQ=255)                                               cgcCGATCAGAACGGCGACCAGGCTTCGTGGCAAACGAAGGTTTTGCACCAGCGCTTCTGGTAGCGTTGGCGTGTGTCCAGGCAGCAGGGTGCGGGTTGCATCTGCTCCGGAAACAGGAATGGCCGAGTAGCAAAACAGa                                                                            >  1:365962/1‑140 (MQ=255)                                                    aTCAGAACGGCGACCAGGCTTCGTGGCAAACGAAGGTTTTGCACCAGCGCTTCTGGTAGCGTTGGCGTGTGTCCAGGCAGCAGGGTGCGGGTTGCATCTGCTCCGGAAACAGGAATGGCCGAGTAGCAAAACAGACTCAg                                                                       >  1:768946/1‑140 (MQ=255)                                                    aTCAGAACGGCGACCAGGCTTCGTGGCAAACGAAGGTTTTGCACCAGCGCTTCTGGTAGCGTTGGCGTGTGTCCAGGCAGCAGGGTGCGGGTTGCATCTGCTCCGGAAACAGGAATGGCCGAGTAGCAAAACAGACTCAg                                                                       >  2:756378/1‑140 (MQ=255)                                                    aTCAGAACGGCGACCAGGCTTCGTGGCAAACGAAGGTTTTGCACCAGCGCTTCTGGTAGCGTTGGCGTGTGTCCAGGCAGCAGGGTGCGGGTTGCATCTGCTCCGGAAACAGGAATGGCCGAGTAGCAAAACAGACTCAg                                                                       >  1:268908/1‑140 (MQ=255)                                                     tCAGAACGGCGACCAGGCTTCGTGGCAAACGAAGGTTTTGCACCAGCGCTTCTGGTAGCGTTGGCGTGTGTCCAGGCAGCAGGGTGCGGGTTGCATCTGCTCCGGAAACAGGAATGGCCGAGTAGCAAAACAGACTCAGc                                                                      <  1:826889/140‑1 (MQ=255)                                                     tCAGAACGGCGACCAGGCTTCGTGGCAAACGAAGGTTTTGCACCAGCGCTTCTGGTAGCGTGGGGGTGTGTCCAGGCAGCTGGGTGCGGGTTGCATCTGCTCCGGAAACAGGAATGGCCGAGTAGCAAAACAGACGCAGc                                                                      <  2:243563/140‑1 (MQ=255)                                                         aaCGGCGACCAGGCTTCGTGGCAAACGAAGGTTTTGCACCAGCGCTTCTGGTAGCGTTGGCGTGTGTCCAGGCAGCAGGGTGCGGGTTGCATCTGCTCCGGAAACAGGAATGGCCGAGTAGCAAAACAGACTCAGCCAGa                                                                  <  2:84225/140‑1 (MQ=255)                                                         aaCGGCGACCAGGCTTCGTGGCAAACGAAGGTTTTGCACCAGCGCTTCTGGTAGCGTTGGCGTGTGTCCAGGCAGCAGGGTGCGGGTTGCATCTGCTCCGGAAACAGGAATGGCCGAGTAGCAAAACAGACTCAGCCAGa                                                                  <  1:568814/140‑1 (MQ=255)                                                         aaCGGCGACCAGGCTTCGTGGCAAACGAAGGTTTTGCACCAGCGCTTCTGGTAGCGTTGGCGTGTGTCCAGGCAGCAGGGTGCGGGTTGCATCTGCTCCGGAAACAGGAATGGCCGAGTAGCAAAACAGACTCAGCCAGa                                                                  <  2:673732/140‑1 (MQ=255)                                                            ggCGACCAGGCTTCGTGGCAAACGAAGGTTTTGCACCAGCGCTTCTGGTAGCGTTGGCGTGTGTCCAGGCAGCAGGGTGCGGGTTGCATCTGCTCCGGAAACAGGAATGGCCGAGTAGCAAAACAGACTCAGCCAGaaaa                                                               <  1:719368/140‑1 (MQ=255)                                                              cGACCAGGCTTCGTGGCAAACGAAGGTTTTGCACCAGCGCTTCTGGTAGCGTTGGCGTGTGTCCAGGCAGCAGGGTGCGGGTTGCATCTGCTCCGGAAACAGGAATGGCCGAGTAGCAAAACAGACTCAGCCAGAAaata                                                             <  2:975757/140‑1 (MQ=255)                                                              cGACCAGGCTTCGTGGCAAACGAAGGTTTTGCACCAGCGCTTCTGGTAGCGTTGGCGTGTGTCCAGGCAGCAGGGTGCGGGTTGCATCTGCTCCGGAAACAGGAATGGCCGAGTAGCAAAACAGACTCAGCCAGAAaata                                                             <  2:858946/140‑1 (MQ=255)                                                                   aGGCTTCGTGGCAAACGAAGGTTTTGCACCAGCGCTTCTGGTAGCGTTGGCGTGTGTCCAGGCAGCAGGGTGCGGGTTGCATCTGCTCCGGAAACAGGAATGGCCGAGTAGCAAAACAGACTCAGCCAGAAaataataat                                                        >  1:285275/1‑140 (MQ=255)                                                                            ggCAAACGAAGGTTTTGCACCAGCGCTTCTGGTAGCGTTGGCGTGTGTCCAGGCAGCAGGGTGCGGGTTGCATCTGCTCCGGAAACAGGAATGGCCGAGTAGCAAAACAGACTCAGCCAGAAAATAATAATAAGtgctgc                                               >  1:48921/1‑140 (MQ=255)                                                                                 aCGAAGGTTTTGCACCAGCGCTTCTGGTAGCGTTGGCGTGTGTCCAGGCAGCAGGGTGCGGGTTGCATCTGCTCCGGAAACAGGAATGGCCGAGTAGCAAAACAGACTCAGCCAGAAAATAATAATAAGTGCTGCAACgg                                          <  2:769277/140‑1 (MQ=255)                                                                                    aaGGTTTTGCACCAGCGCTTCTGGTAGCGTTGGCGTGTGTCCAGGCAGCAGGGTGCGGGTTGCATCTGCTCCGGAAACAGGAATGGCCGAGTAGCAAAACAGACTCAGCCAGAAAATAATAATAAGTGCTGCAACGGGaa                                       <  1:962823/140‑1 (MQ=255)                                                                                       gTTTTGCACCAGCGCTTCTGGTAGCGTTGGCGTGTGTCCAGGCAGCAGGGTGCGGGTTGCATCTGCTCCGGAAACAGGAATGGCCGAGTAGCAAAACAGACTCAGCCAGAAAATAATAATAAGTGCTGCAACGGGAAGcc                                    >  2:805649/1‑140 (MQ=255)                                                                                         tttGCACCAGCGCTTCTGGTAGCGTTGGCGTGTGTCCAGGCAGCAGGGTGCGGGTTGCATCTGCTCCGGAAACAGGAATGGCCGAGTAGCAAAACAGACTCAGCCAGAAAATAATAATAAGTGCTGCAACGGGAAGcccc                                  <  2:952670/140‑1 (MQ=255)                                                                                                       tCTGGTAGCGTTGGCGTGTGTCCAGGCAGCAGGGTGCGGGTTGCATCTGCTCCGGAAACAGGAATGGCCGAGTAGCAAAACAGACTCAGCCAGAAAATAATAATAAGTGCTGCAACGGGAAGCCCCCACAGCAGCACCgg                    <  2:933072/140‑1 (MQ=255)                                                                                                       tCTGGTAGCGTTGGCGTGTGTCCAGGCAGCAGGGTGCGGGTTGCATCTGCTCCGGAAACAGGAATGGCCGAGTAGCAAAACAGACTCAGCCAGAAAATAATAATAAGTGCTGCAACGGGAAGCCCCCACAGCAGCACCgg                    >  2:540418/1‑140 (MQ=255)                                                                                                        cTGGTAGCGTTGGCGTGTGTCCAGGCAGCAGGGTGCGGGTTGCATCTGCTCCGGAAACAGGAATGGCCGAGTAGCAAAACAGACTCAGCCAGAAAATAATAATAAGTGCTGCAACGGGAAGCCCCCACAGCAGCACCggg                   <  1:13720/140‑1 (MQ=255)                                                                                                             aGCGTTGGCGTGTGTCCAGGCAGCAGGGTGCGGGTTGCATCTGCTCCGGAAACAGGAATGGCCGAGTAGCAAAACAGACTCAGCCAGAAAATAATAATAAGTGCTGCAACGGGAAGCCCCCACAGCAGCACCGGGTGttt              >  1:267678/1‑140 (MQ=255)                                                                                                                 ttGGCGTGTGTCCAGGCAGCAGGGTGCGGGTTGCATCTGCTCCGGAAACAGGAATGGCCGAGTAGCAAAACAGACTCAGCCAGAAAATAATAATAAGTGCTGCAACGGGAAGCCCCCACAGCAGCACCGGGCGTTTTATc          >  1:47030/1‑140 (MQ=255)                                                                                                                    gCGTGTGTCCAGGCAGCAGGGTGCGGGTTGCATCTGCTCCGGAAACAGGAATGGCCGAGTAGCAAAACAGACTCAGCCAGAAAATAATAATAAGTGCTGCAACGGGAAGCCCCCACAGCAGCACCGGGTGTTTTATcgcg       <  2:356168/140‑1 (MQ=255)                                                                                                                      gtgtgtCCAGGCAGCAGGGTGCGGGTTGCATCTGCTCCGGAAACAGGAATGGCCGAGTAGCAAAACAGACTCAGCCAGAAAATAATAATAAGTGCTGCAACGGGAAGCCCCCACAGCAGCACCGGGTGTTTTATCGCGGt     <  1:401461/140‑1 (MQ=255)                                                                                                                         tgtCCAGGCAGCAGGGTGCGGGTTGCATCTGCTCCGGAAACAGGAATGGCCGAGTAGCAAAACAGACTCAGCCAGAAAATAATAATAAGTGCTGCAACGGGAAGCCCCCACAGCAGCACCGGGTGTTTTATCGCGGTCAt  >  2:48050/1‑140 (MQ=255)                                                                                                                         tgtCCAGGCAGCAGGGTGCGGGTTGCATCTGCTCCGGAAACAGGAATGGCCGAGTAGCAAAACAGACTCAGCCAGAAAATAATAATAAGTGCTGCAACGGGAAGCCCCCACAGCAGCACCGGGTGTTTTATCGCGGTCAt  >  1:759796/1‑140 (MQ=255)                                                                                                                                         |                                                                                                                            TTGTGGGTCAGGGTTTGCAGCAGCGTGCCCGCGAGCGCCAGGCTTGCGCCGATCAGAACGGCGACCAGGCTTCGTGGCAAACGAAGGTTTTGCACCAGCGCTTCTGGTAGCGTTGGCGTGTGTCCAGGCAGCAGGGCGCGGGTTGCATCTGCTCCGGAAACAGGAATGGCCGAGTAGCAAAACAGACTCAGCCAGAAAATAATAATAAGTGCTGCAACGGGAAGCCCCCACAGCAGCACCGGGTGTTTTATCGCGGTCAT  >  NC_000913/4513150‑4513409

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