Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	NC_000913	4,512,668	G→A	L248F (CTT→TTT)	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,512,668	0	G	A	100.0%	204.3 / NA	58	L248F (CTT→TTT)	fecC	ferric citrate ABC transporter membrane subunit FecC
Reads supporting (aligned to +/- strand):  ref base G (0/0);  new base A (24/34);  total (24/34)

AGCGTGGCCCCCAGCAGCATGCTCACTGGCAGTACGTTGCGCTGATCGAAGCCTGCCCAGAAGCGCGCCAGATGTGGCACCAGCAGACCGATAAACGCCACCGGACCTGCTACGCTGACGCACGCACCAACCAGAAGCAGCACTAACATATTGATGACCAAACGTAGCCTCGTCAGGTTCACTCCCAGCGTATGGGCGGTGCTGTCGCTGAGGTTGAGCAGGTTCAGTTGATTCGCCAGCAGCAACACGACAGGGACTGCAGTGACC  >  NC_000913/4512532‑4512798                                                                                                                                         |                                                                                                                                   aGCGTGGCCCCCAGCAGCATGCTCACTGGCAGTACGTTGCGCTGATCGAAGCCTGCCCAGAAGCGCGCCAGATGTGGCACCAGCAGACCGATAAACGCCACCGGACCTGCTACGCTGACGCACGCACCAACCAGAaacag                                                                                                                                 <  1:433233/140‑1 (MQ=255)   cGTGGCCCCCAGCAGCATGCTCACTGGCAGTACGTTGCGCTGATCGAAGCCTGCCCAGAAGCGCGCCAGATGTGGCACCAGCAGACCGATAAACGCCACCGGACCTGCTACGCTGACGCACGCACCAACCAGAAacagca                                                                                                                               <  2:866872/140‑1 (MQ=255)    gTGGCCCCCAGCAGCATGCTCACTGGCAGTACGTTGCGCTGATCGAAGCCTGCCCAGAAGCGCGCCAGATGTGGCACCAGCAGACCGATAAACGCCACCGGACCTGCTACGCTGACGCACGCACCAACCAGAAACAGCAc                                                                                                                              <  1:1074792/140‑1 (MQ=255)    gTGGCCCCCAGCAGCATGCTCACTGGCAGTACGTTGCGCTGATCGAAGCCTGCCCAGAAGCGCGCCAGATGTGGCACCAGCAGACCGATAAACGCCACCGGACCTGCTACGCTGACGCACGCACCAACCAGAAACAGCAc                                                                                                                              <  2:493001/140‑1 (MQ=255)    gTGGCCCCCAGCAGCATGCTCACTGGCAGTACGTTGCGCTGATCGAAGCCTGCCCAGAAGCGCGCCAGATGTGGCACCAGCAGACCGATAAACGCCACCGGACCTGCTACGCTGACGCACGCACCAACCAGAAACAGCAc                                                                                                                              <  2:903769/140‑1 (MQ=255)    gTGGCCCCCAGCAGCATGCTCACTGGCAGTACGTTGCGCTGATCGAAGCCTGCCCAGAAGCGCGCCAGATGTGGCACCAGCAGACCGATAAACGCCACCGGACCTGCTACGCTGACGCACGCACCAACCAGAAACAGCAc                                                                                                                              <  2:145298/140‑1 (MQ=255)     tGGCCCCCAGCAGCATGCTCACTGGCAGTACGTTGCGCTGATCGAAGCCTGCCCAGAAGCGCGCCAGATGTGGCACCAGCAGACCGATAAACGCCACCGGACCTGCTACGCTGACGCACGCACCAACCAGAAACAGCACt                                                                                                                             <  1:27992/140‑1 (MQ=255)     tGGCCCCCAGCAGCATGCTCACTGGCAGTACGTTGCGCTGATCGAAGCCTGCCCAGAAGCGCGCCAGATGTGGCACCAGCAGACCGATAAACGCCACCGGACCTGCTACGCTGACGCACGCACCAACCAGAAACAGCACt                                                                                                                             <  1:65000/140‑1 (MQ=255)     tGGCCCCCAGCAGCATGCTCACTGGCAGTACGTTGCGCTGATCGAAGCCTGCCCAGAAGCGCGCCAGATGTGGCACCAGCAGACCGATAAACGCCACCGGACCTGCTACGCTGACGCACGCACCAACCAGAAACAGCACt                                                                                                                             <  1:393338/140‑1 (MQ=255)     tGGCCCCCAGCAGCATGCTCACTGGCAGTACGTTGCGCTGATCGAAGCCTGCCCAGAAGCGCGCCAGATGTGGCACCAGCAGACCGATAAACGCCACCGGACCTGCTACGCTGACGCACGCACCAACCAGAAACAGCACt                                                                                                                             >  2:34969/1‑140 (MQ=255)       gCCCCCAGCAGCATGCTCACTGGCAGTACGTTGCGCTGATCGAAGCCTGCCCAGAAGCGCGCCAGATGTGGCACCAGCAGACCGATAAACGCCACCGGACCTGCTACGCTGACGCACGCACCAACCAGAaacag                                                                                                                                 <  1:394941/134‑1 (MQ=255)       gCCCCCAGCAGCATGCTCACTGGCAGTACGTTGCGCTGATCGAAGCCTGCCCAGAAGCGCGCCAGATGTGGCACCAGCAGACCGATAAACGCCACCGGACCTGCTACGCTGACGCACGCACCAACCAGAaacag                                                                                                                                 >  2:394941/1‑134 (MQ=255)        cccccAGCAGCATGCTCACTGGCAGTACGTTGCGCTGATCGAAGCCTGCCCAGAAGCGCGCCAGATGTGGCACCAGCAGACCGATAAACGCCACCGGACCTGCTACGCTGACGCACGCACCAACCAGAAACAGCACTAAc                                                                                                                          <  2:1170181/140‑1 (MQ=255)           ccAGCAGCATGCTCACTGGCAGTACGTTGCGCTGATCGAAGCCTGCCCAGAAGCGCGCCAGATGTGGCACCAGCAGACCGATAAACGCCACCGGACCTGCTACGCTGACGCACGCACCAACCAGAAACAGCACTAACata                                                                                                                       >  1:130743/1‑140 (MQ=255)               cagcaTGCTCACTGGCAGTACGTTGCGCTGATCGAAGCCTGCCCAGAAGCGCGCCAGATGTGGCACCAGCAGACCGATAAACGCCACCGGACCTGCTACGCTGACGCACGCACCAACCAGAAACAGCACTAACATATtga                                                                                                                   <  1:1080249/140‑1 (MQ=255)                  caTGCTCACTGGCAGTACGTTGCGCTGATCGAAGCCTGCCCAGAAGCGCGCCAGATGTGGCACCAGCAGACCGATAAACGCCACCGGACCTGCTACGCTGACGCACGCACCAACCAGAAACAGCACTAACATATtgatga                                                                                                                >  1:344112/1‑140 (MQ=255)                  caTGCTCACTGGCAGTACGTTGCGCTGATCGAAGCCTGCCCAGAAGCGCGCCAGATGTGGCACCAGCAGACCGATAAACGCCACCGGACCTGCTACGCTGACGCACGCACCAACCAGAAACAGCACTAACATATtgatga                                                                                                                >  2:1095507/1‑140 (MQ=255)                   aTGCTCACTGGCAGTACGTTGCGCTGATCGAAGCCTGCCCAGAAGCGCGCCAGATGTGGCACCAGCAGACCGATAAACGCCACCGGACCTGCTACGCTGACGCACGCACCAACCAGAAACAGCACTAACATATTGATGAc                                                                                                               >  2:84003/1‑140 (MQ=255)                    tGCTCACTGGCAGTACGTTGCGCTGATCGAAGCCTGCCCAGAAGCGCGCCAGATGTGGCACCAGCAGACCGATAAACGCCACCGGACCTGCTACGCTGACGCACGCACCAACCAGAAACAGCACTAACATATTGATGAcc                                                                                                              <  2:753152/140‑1 (MQ=255)                    tGCTCACTGGCAGTACGTTGCGCTGATCGAAGCCTGCCCAGAAGCGCGCCAGATGTGGCACCAGCAGACCGATAAACGCCACCGGACCTGCTACGCTGACGCACGCACCAACCAGAAACAGCACTAACATATTGATGAcc                                                                                                              >  1:103950/1‑140 (MQ=255)                           tGGCAGTACGTTGCGCTGATCGAAGCCTGCCCAGAAGCGCGCCAGATGTGGCACCAGCAGACCGATAAACGCCACCGGACCTGCTACGCTGACGCACGCACCAACCAGAAACAGCACTAACATATTGATGACCAAACGTa                                                                                                       <  2:111103/140‑1 (MQ=255)                            ggCAGTACGTTGCGCTGATCGAAGCCTGCCCAGAAGCGCGCCAGATGTGGCACCAGCAGACCGATAAACGCCATCGGACCTGCTACGCTGACGCACGCACCAACCAGAAACAGCACTAACATATCGATGACCAAACGTAg                                                                                                      <  1:1082578/140‑1 (MQ=255)                                 tACGTTGCGCTGATCGAAGCCTGCCCAGAAGCGCGCCAGATGTGGCACCAGCAGACCGATAAACGCCACCGGACCTGCTACGCTGACGCACGCACCAACCAGAAACAGCACTAACATATTGATGACCAAACGTAGCCTCg                                                                                                 >  1:32498/1‑140 (MQ=255)                                       gcgcTGATCGAAGCCTGCCCAGAAGCGCGCCAGATGTGGCACCAGCAGACCGATAAACGCCACCGGACCTGCTACGCTGACGCACGCACCAACCAGAAACAGCACTAACATATTGATGACCAAACGTAGCCTCGTCAGGt                                                                                           >  2:549649/1‑140 (MQ=255)                                        cgcTGATCGAAGCCTGCCCAGAAGCGCGCCAGATGTGGCACCAGCAGACCGATAAACGCCACCGGACCTGCTACGCTGACGCACGCACCAACCAGAAACAGCACTAACATATTGATGACCAAACGTAGCCTCGTCAGGtt                                                                                          >  1:8744/1‑140 (MQ=255)                                             aTCGAAGCCTGCCCAGAAGCGCGCCAGATGTGGCACCAGCAGACCGATAAACGCCACCGGACCTGCTACGCTGACGCACGCACCAACCAGAAACAGCACTAACATATTGATGACCAAACGTAGCCTCGTCAGGTTCACTc                                                                                     >  2:845915/1‑140 (MQ=255)                                             aTCGAAGCCTGCCCAGAAGCGCGCCAGATGTGGCACCAGCAGACCGATAAACGCCACCGGACCTGCTACGCTGACGCACGCACCAACCAGAAACAGCACTAACATATTGATGACCAAACGTAGCCTCGTCAGGTTCACTc                                                                                     >  2:1080347/1‑140 (MQ=255)                                              tCGAATCCTGCCCAGAAGCGCGCCAGCTGTGGCACCAGCAGACCGATAAACGCCACCGGACCTGCTACGCTGACGCACGCACCAACCAGAAACAGCACTAACATATTGATGACCAAACGTAGCCTCGTCAGGTTCACTcc                                                                                    <  2:1076432/140‑1 (MQ=255)                                              tCGAAGCCTGCCCAGAAGCGCGCCAGATGTGGCACCAGCAGACCGATAAACGCCACCGGACCTGCTACGCTGACGCACGCACCAACCAGAAACAGCACTAACATATTGATGACCAAACGTAGCCTCGTCAGGTTCACTcc                                                                                    <  1:845915/140‑1 (MQ=255)                                              tCGAAGCCTGCCCAGAAGCGCGCCAGATGTGGCACCAGCAGACCGATAAACGCCACCGGACCTGCTACGCTGACGCACGCACCAACCAGAAACAGCACTAACATATTGATGACCAAACGTAGCCTCGTCAGGTTCACTcc                                                                                    <  1:85420/140‑1 (MQ=255)                                              tCGAAGCCTGCCCAGAAGCGCGCCAGATGTGGCACCAGCAGACCGATAAACGCCACCGGACCTGCTACGCTGACGCACGCACCAACCAGAAACAGCACTAACATATTGATGACCAAACGTAGCCTCGTCAGGTTCACTcc                                                                                    <  2:147605/140‑1 (MQ=255)                                              tCGAAGCCTGCCCAGAAGCGCGCCAGATGTGGCACCAGCAGACCGATAAACGCCACCGGACCTGCTACGCTGACGCACGCACCAACCAGAAACAGCACTAACATATTGATGACCAAACGTAGCCTCGTCAGGTTCACTcc                                                                                    <  1:246532/140‑1 (MQ=255)                                               cGAAGCCTGCCCAGAAGCGCGCCAGATGTGGCACCAGCAGACCGATAAACGCCACCGGACCTGCTACGCTGACGCACGCACCAACCAGAAACAGCACTAACATATTGATGACCAAACGTAGCCTCGTCAGGTTCACTccc                                                                                   <  1:84003/140‑1 (MQ=255)                                                    ccTGCCCAGAAGCGCGCCAGATGTGGCACCAGCAGACCGATAAACGCCACCGGACCTGCTACGCTGACGCACGCACCAACCAGAAACAGCACTAACATATTGATGACCAAACGTAGCCTCGTCAGGTTCACTCCCAGCGt                                                                              <  2:32498/140‑1 (MQ=255)                                                       gCCCAGAAGCGCGCCAGATGTGGCACCAGCAGACCGATAAACGCCACCGGACCTGCTACGCTGACGCACGCACCAACCAGAAACAGCACTAACATATTGATGACCAAACGTAGCCTCGTCAGGTTCACTCCCAGCGTATg                                                                           >  1:564825/1‑140 (MQ=255)                                                                cgcgcCAGATGTGGCACCAGCAGACCGATAAACGCCACCGGACCTGCTACGCTGACGCACGCACCAACCAGAAACAGCACTAACATATTGATGACCAAACGTAGCCTCGTCAGGTTCACTCCCAGCGTATGGGCGGTGCt                                                                  <  2:103950/140‑1 (MQ=255)                                                                cgcgcCAGATGTGGCACCAGCAGACCGATAAACGCCACCGGACCTGCTACGCTGACGCACGCACCAACCAGAAACAGCACTAACATATTGATGACCAAACGTAGCCTCGTCAGGTTCACTCCCAGCGTATGGGCGGTGCt                                                                  <  1:673749/140‑1 (MQ=255)                                                                  cgcCAGATGTGGCACCAGCAGACCGATAAACGCCACCGGACCTGCTACGCTGACGCACGCACCAACCAGAAACAGCACTAACATATTGATGACCAAACGTAGCCTCGTCAGGTTCACTCCCAGCGTATGGGCGGTGCTGt                                                                >  2:103386/1‑140 (MQ=255)                                                                     cAGATGTGGCACCAGCAGACCGATAAACGCCACCGGACCTGCTACGCTGACGCACGCACCAACCAGAAACAGCACTAACATATTGATGACCAAACGTAGCCTCGTCAGGTTCACTCCCAGCGTATGGGCGGTGCTGTCGc                                                             <  2:8744/140‑1 (MQ=255)                                                                               aCCAGCAGACCGATAAACGCCACCGGACCTGCTACGCTGACGCACGCACCAACCAGAAACAGCACTAACATATTGATGACCAAACGTAGCCTCGTCAGGTTCACTCCCAGCGTATGGGCGGTGCTGTCGCTGAGGTTGAg                                                   >  2:1014180/1‑140 (MQ=255)                                                                                   gcagACCGATAAACGCCACCGGACCTGCTACGCTGACGCACGCACCAACCAGAAACAGCACTAACATATTGATGACCAAACGTAGCCTCGTCAGGTTCACTCCCAGCGTATGGGCGGTGCTGTCGCTGAGGTTGAGCAgg                                               >  1:944340/1‑140 (MQ=255)                                                                                    cagACCGATAAACGCCACCGGACCTGCTACGCTGACGCACGCACCAACCAGAAACAGCACTAACATATTGATGACCAAACGTAGCCTCGTCAGGTTCACTCCCAGCGTATGGGCGGTGCTGTCGCTGAGGTTGAGGAGGt                                              >  2:1118367/1‑140 (MQ=255)                                                                                      gACCGATAAACGCCACCGGACCTGCTACGCTGACGCACGCACCAACCAGAAACAGCACTAACATATTGATGACCAAACGTAGCCTCGTCAGGTTCACTCCCAGCGTATGGGCGGTGCTGTCGCTGAGGTTGAGCAGGTTc                                            <  2:917038/140‑1 (MQ=255)                                                                                                   cACCGGACCTGCTACGCTGACGCACGCACCAACCAGAAACAGCACTAACATATTGATGACCAAACGTAGCCTCGTCAGGTTCACTCCCAGCGTATGGGCGGTGCTGTCGCTGAGGTTGAGCAGGTTCAGTTGATTCGCca                               <  1:1095507/140‑1 (MQ=255)                                                                                                     ccGGACCTGCTACGCTGACGCACGCACCAACCAGAAACAGCACTAACATATTGATGACCAAACGTAGCCTCGTCAGGTTCACTCCCAGCGTATGGGCGGTGCTGTCGCTGAGGTTGAGCAGGTTCAGTTGATTCGCcagc                             <  2:130743/140‑1 (MQ=255)                                                                                                      cGGACCTGCTACGCTGACGCACGCACCAACCAGAAACAGCACTAACATATTGATGACCAAACGTAGCCTCGTCAGGTTCACTCCCAGCGTATGGGCGGTGCTGTCGCTGAGGTTGAGCAGGTTCAGTTGATTCGCcagca                            >  2:377806/1‑140 (MQ=255)                                                                                                         aCCTGCTACGCTGACGCACGCACCAACCAGAAACAGCACTAACATATTGATGACCAAACGTAGCCTCGTCAGGTTCACTCCCAGCGTATGGGCGGTGCTGTCGCTGAGGTTGAGCAGGTTCAGTTGATTCGCCAgcagca                         <  1:1118367/140‑1 (MQ=255)                                                                                                         aCCTGCTACGCTGACGCACGCACCAACCAGAAACAGCACTAACATATTGATGACCAAACGTAGCCTCGTCAGGTTCACTCCCAGCGTATGGGCGGTGCTGTCGCTGAGGTTGAGCAGGTTCAGTTGATTCGCCAgcagca                         <  1:54770/140‑1 (MQ=255)                                                                                                                aCGCTGACGCACGCACCAACCAGAAACAGCACTAACATATTGATGACCAAACGTAGCCTCGTCAGGTTCACTCCCAGCGTATGGGCGGTGCTGTCGCTGAGGTTGAGCAGGTTCAGTTGATTCGCCAGCAGCAACACGAc                  <  2:671210/140‑1 (MQ=255)                                                                                                                     gACGCACGCACCAACCAGAAACAGCACTAACATATTGATGACCAAACGTAGCCTCGTCAGGTTCACTCCCAGCGTATGGGCGGTGCTGTCGCTGAGGTTGAGCAGGTTCAGTTGATTCGCCAGCAGCAACACGACAGGGa             <  1:377806/140‑1 (MQ=255)                                                                                                                     gACGCACGCACCAACCAGAAACAGCACTAACATATTGATGACCAAACGTAGCCTCGTCAGGTTCACTCCCAGCGTATGGGCGGTGCTGTCGCTGAGGTTGAGCAGGTTCAGTTGATTCGCCAGCAGCAACACGACAGGGa             <  1:1014180/140‑1 (MQ=255)                                                                                                                       cgcacgcaCCAACCAGAAACAGCACTAACATATTGATGACCAAACGTAGCCTCGTCAGGTTCACTCCCAGCGTATGGGCGGTGCTGTCGCTGAGGTTGAGCAGGTTCAGTTGATTCGCCAGCAGCAACACGACAGGGACt           <  2:187557/140‑1 (MQ=255)                                                                                                                        gcacgcacCAACCAGAAACAGCACTAACATATTGATGACCAAACGTAGCCTCGTCAGGTTCACTCCCAGCGTATGGGCGGTGCTGTCGCTGAGGTTGAGCAGGTTCAGTTGATTCGCCAGCAGCAACACGACAGGGACTg          >  2:361935/1‑140 (MQ=255)                                                                                                                           cgcacCAACCAGAAACAGCACTAACATATTGATGACCAAACGTAGCCTCGTCAGGTTCACTCCCAGCGTATGGGCGGTGCTGTCGCTGAGGTTGAGCAGGTTCAGTTGATTCGCCAGCAGCAACACGACAGGGACTg          <  2:1126489/137‑1 (MQ=255)                                                                                                                           cgcacCAACCAGAAACAGCACTAACATATTGATGACCAAACGTAGCCTCGTCAGGTTCACTCCCAGCGTATGGGCGGTGCTGTCGCTGAGGTTGAGCAGGTTCAGTTGATTCGCCAGCAGCAACACGACAGGGACTg          >  1:1126489/1‑137 (MQ=255)                                                                                                                           cgcacCAACCAGAAACAGCACTAACATATTGATGACCAAACGTAGCCTCGTCAGGTTCACTCCCAGCGTATGGGCGGTGCTGTCGCTGAGGTTGAGCAGGTTCAGTTGATTCGCCAGCAGCAACACGACAGGGACTGCAg       >  2:3698/1‑140 (MQ=255)                                                                                                                             cacCAACCAGAAACAGCACTAACATATTGATGACCAAACGTAGCCTCGTCAGGTTCACTCCCAGCGTATGGGCGGTGCTGTCGCTGAGGTTGAGCAGGTTCAGTTGATTCGCCAGCAGCAACACGACAGGGACTGCAGTg     >  2:997874/1‑140 (MQ=255)                                                                                                                              accaaccaGAAACAGCACTAACATATTGATGACCAAACGCAGCCTCGTCAGGTTCACTCCCAGCGTATGGGCGGTGCTGTCGCTGAGGTTGAGCAGGTTCAGTTGATTCGCCAGCAGCAACACGACAGGGACTGCAGTGa    >  2:507082/1‑140 (MQ=255)                                                                                                                               ccaaccaGAAACAGCACTAACATATTGATGACCAAACGTAGCCTCGTCAGGTTCACTCCCAGCGTATGGGCGGTGCTGTCGCTGAGGTTGAGCAGGTTCAGTTGATTCGCCAGCAGCAACACGACAGGGACTGCAGTGac   <  1:3698/140‑1 (MQ=255)                                                                                                                                caaccaGAAACAGCACTAACATATTGATGACCAAACGTAGCCTCGTCAGGTTCACTCCCAGCGTATGGGCGGTGCTGTCGCTGAGGTTGAGCAGGTTCAGTTGATTCGCCAGCAGCAACACGAc                  >  2:864305/1‑124 (MQ=255)                                                                                                                                caaccaGAAACAGCACTAACATATTGATGACCAAACGTAGCCTCGTCAGGTTCACTCCCAGCGTATGGGCGGTGCTGTCGCTGAGGTTGAGCAGGTTCAGTTGATTCGCCAGCAGCAACACGAc                  <  1:864305/124‑1 (MQ=255)                                                                                                                                caaccaGAAACAGCACTAACATATTGATGACCAAACGTAGCCTCGTCAGGTTCACTCCCAGCGTATGGGCGGTGCTGTCGCTGAGGTTGAGCAGGTTCAGTTGATTCGCCAGCAGCAACACGACAGGGACTGCAGTGacc  >  1:787327/1‑140 (MQ=255)                                                                                                                                         |                                                                                                                                   AGCGTGGCCCCCAGCAGCATGCTCACTGGCAGTACGTTGCGCTGATCGAAGCCTGCCCAGAAGCGCGCCAGATGTGGCACCAGCAGACCGATAAACGCCACCGGACCTGCTACGCTGACGCACGCACCAACCAGAAGCAGCACTAACATATTGATGACCAAACGTAGCCTCGTCAGGTTCACTCCCAGCGTATGGGCGGTGCTGTCGCTGAGGTTGAGCAGGTTCAGTTGATTCGCCAGCAGCAACACGACAGGGACTGCAGTGACC  >  NC_000913/4512532‑4512798

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