Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	NC_000913	4,512,124	(C)5→6	coding (291/957 nt)	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,119	1	.	C	100.0%	187.5 / NA	51	A99G (GCT→GGT)	fecD	ferric citrate ABC transporter membrane subunit FecD
Reads supporting (aligned to +/- strand):  ref base . (0/0);  new base C (31/20);  total (31/20)

GCCACGCCGGTGAGCGCCAGCTTCATCGGCTGGTGGGTCTTTGCCAGCATCTTCAGCAATATCAACCCCGCCATGCCGCCCGCAAAGGCCAGCAGCGGCAGCACCATCACGGGCAGTGACGGCATAAGAAGTAGAG‑CCCCCACAGAGGCCAGGCTGGCGGCATGGTTAACACCGAGAATATCCGGTGATGCCAGAGGGTTGCGCACAATCCCCTGTATCAGCACGCCCGCCACGGCGAGGGCTGCACCGACAAACAGTGCCAGCAGC  >  NC_000913/4511984‑4512250                                                                                                                                         |                                                                                                                                    cccgcGCCGGTGAGCGCCAGCTTCATCGGCGGGTGGGTTTTTGCCAGCATCTTCAGCAATACCACCCCCGCCATGCCGCCCGCAAAGGCCAGCAGCGGCAGCACCATCACGGGCAGTGACGGCATAAGAAGTAGAG‑cccc                                                                                                                                 <  2:136559/136‑1 (MQ=255) cccACGCGGGTGAGCCCCAGTTTCAGCGGCTGGTGGTTTTTTGCCAGCATCTTCAGCAATTTCAACCCCGCCATGCCGCCCGCAAAGGCCAGCAGCGGCAGCACCATCACGGGCAGTGACGGCATAAGAAGTAGAG‑cccc                                                                                                                                 <  2:27639/139‑1 (MQ=255)   cACGCCGGTGAGCGCCAGCTTCATCGGCTGGTGGGTCTTTGCCAGCATCTTCAGCAATATCAACCCCGCCATGCCGCCCGCAAAGGCCAGCAGCGGCAGCACCATCACGGGCAGTGACGGCATAAGAAGTAGAG‑cccccc                                                                                                                               >  1:453118/1‑139 (MQ=255)    aCGCCGGTGAGCGCCAGCTTCATCGGCTGGTGGGTCTTTGCCAGCATCTTCAGCAATATCAACCCCGCCATGCCGCCCGCAAAGGCCAGCAGCGGCAGCACCATCACGGGCAGTGACGGCATAAGAAGTAGAG‑cccccca                                                                                                                              <  2:429482/140‑3 (MQ=255)     cGCCGGTGAGCGCCAGCTTCATCGGCTGGTGGGTCTTTGCCAGCATCTTCAGCAATATCAACCCCGCCATGCCGCCCGCAAAGGCCAGCAGCGGCAGCACCATCACGGGCAGTGACGGCATAAGAAGTAGAG‑ccccccac                                                                                                                             >  1:889669/1‑137 (MQ=255)               cgcCAGCTTCATCGGCTGGTGGGTCTTTGCCAGCATCTTCAGCAATATCACCCCCGCCATGCCGCCCGCAAAGGCCAGCAGCGGCAGCACCATCACGGGCAGTGACGGCATAAGAAGTAGAGCCCCCCACAGAGGCCAgg                                                                                                                    <  1:45351/140‑1 (MQ=255)               cgcCAGCTTCATCGGCTGGTGGGTCTTTGCCAGCATCTTCAGCAATATCAACCCCGCCATGCCGCCCGCAAAGGCCAGCAGCGGCAGCACCATCACGGGCAGTGACGGCATAAGAAGTAGAGCCCCCCACAGAGGCCAgg                                                                                                                    >  1:398714/1‑140 (MQ=255)                  cAGCTTCATCGGCTGGTGGGTCTTTGCCAGCATCTTCAGCAATATCAACCCCGCCATGCCGCCCGCAAAGGCCAGCAGCGGCAGCACCATCACGGGCAGTGACGGCATAAGAAGTAGAGCCCCCCACAGAGGCCAGGCTg                                                                                                                 >  1:256430/1‑140 (MQ=255)                    gCTTCATCGGCTGGTGGGTCTTTGCCAGCATCTTCAGCAATATCAACCCCGCCATGCCGCCCGCAAAGGCCAGCAGCGGCAGCACCATCACGGGCAGTGACGGCATAAGAAGTAGAGCCCCCCACAGAGGCCAGGCTggc                                                                                                               >  2:673311/1‑140 (MQ=255)                       tccccGGCTGGTGGGTTTTTGCAAGCACCTTCAGCAATAACAACCCCGCCATGCCGCCCGCAAAGGCCAGCAGCGGCAGCACCATCACGGGCAGTGACGGCATAAGAAGTAGAGCCCCCCACAGAGGCCAGGCTggcggc                                                                                                            <  2:129697/136‑1 (MQ=255)                       tCATCGGCTGGTGGGTCTTTGCCAGCATCTTCAGCAATATCAACCCCGCCATGCCGCCCGCAAAGGCCAGCAGCGGCAGCACCATCACGGGCAGTGACGGCATAAGAAGTAGAGCCCCCCACAGAGGCCAGGCTggcggc                                                                                                            >  2:149752/1‑140 (MQ=255)                            ggggggtggGTCTTTGCCAGCATCTTCAGCAATACCACCCCCGCCATGCCGCCCGCAAAGGCCAGCAGCGGCAGCACCATCACGGGCAGTGACGGCATAAGAAGTAGAGCCCCCCACAGAGGCCAGGCTGGCGGCATGGt                                                                                                       <  2:839789/136‑1 (MQ=255)                            ggCTGGTGGGTCTTTGCCAGCATCTTCAGCAATATCAACCCCGCCATGCCGCCCGCAAAGGCCAGCAGCGGCAGCACCATCACGGGCAGTGACGGCATAAGAAGTAGAGCCCCCCACAGAGGCCAGGCTGGCGGCATGGt                                                                                                       >  2:12000/1‑140 (MQ=255)                             gCTGGTGGGTCTTTGCCAGCATCTTCAGCAATATCAACCCCGCCATGCCGCCCGCAAAGGCCAGCAGCGGCAGCACCATCACGGGCAGTGACGGCATAAGAAGTAGAGCCCCCCACAGAGGCCAGGCTGGCGGCATGGtt                                                                                                      <  1:149752/140‑1 (MQ=255)                              cTGGTGGGTCTTTGCCAGCATCTTCAGCAATATCAACCCCGCCATGCCGCCCGCAAAGGCCAGCAGCGGCAGCACCATCACGGGCAGTGACGGCATAAGAAGTAGAGCCCCCCACAGAGGCCAGGCTGGCGGCATGGTTa                                                                                                     <  2:256430/140‑1 (MQ=255)                               tggtggGTCTTTGCCAGCATCTTCAGCAATATCAACCCCGCCATGCCGCCCGCAAAGGCCAGCAGCGGCAGCACCATCACGGGCAGTGACGGCATAAGAAGTAGAGCCCCCCACAGAGGCCAGGCTGGCGGCATGGTTaa                                                                                                    >  2:566619/1‑140 (MQ=255)                                  ggggTTTTTGCCAGCATCTTCAGCAATACCACCCCCGCCATGCCGCCCGCAAAGGCCAGCAGCGGCAGCACCATCACGGGCAGTGACGGCATAAGAAGTAGAGCCCCCCACAGAGGCCAGGCTGGCGGCATGGTTAacac                                                                                                 <  2:306274/139‑1 (MQ=255)                                        tttGCCAGCATCTTCAGCAATATCAACCCCGCCATGCCGCCCGCAAAGGCCAGCAGCGGCAGCACCATCACGGGCAGTGACGGCATAAGAAGTAGAGCCCCCCACAGAGGCCAGGCTGGCGGCATGGTTAACACCGAGaa                                                                                           >  1:415851/1‑140 (MQ=255)                                                cATCTTCAGCAATATCAACCCCGCCATGCCGCCCGCAAAGGCCAGCAGCGGCAGCACCATCACGGGCAGTGACGGCATAAGAAGTAGAGCCCCCCACAGAGGCCAGGCTGGCGGCATGGTTAACACCGAGAATATCCgg                                                                                    <  2:62271/139‑1 (MQ=255)                                                cATCTTCAGCAATATCAACCCCGCCATGCCGCCCGCAAAGGCCAGCAGCGGCAGCACCATCACGGGCAGTGACGGCATAAGAAGTAGAGCCCCCCACAGAGGCCAGGCTGGCGGCATGGTTAACACCGAGAATATCCgg                                                                                    >  1:62271/1‑139 (MQ=255)                                                cATCTTCAGCAATATCAACCCCGCCATGCCGCCCGCAAAGGCCAGCAGCGGCAGCACCATCACGGGCAGTGACGGCATAAGAAGTAGAGCCCCCCACAGAGGCCAGGCTGGCGGCATGGTTAACACCGAGAATATCCGGt                                                                                   <  2:173563/140‑1 (MQ=255)                                                cATCTTCAGCAATATCAACCCCGCCATGCCGCCCGCAAAGGCCAGCAGCGGCAGCACCATCACGGGCAGTGACGGCATAAGAAGTAGAGCCCCCCACAGAGGCCAGGCTGGCGGCATGGTTAACACCGAGAATATCCGGt                                                                                   <  1:673311/140‑1 (MQ=255)                                                cATCTTCAGCAATATCAACCCCGCCATGCCGCCCGCAAAGGCCAGCAGCGGCAGCACCATCACGGGCAGTGACGGCATAAGAAGTAGAGCCCCCCACAGAGGCCAGGCTGGCGGCATGGTTAACACCGAGAATATCCGGt                                                                                   >  1:261249/1‑140 (MQ=255)                                                  tCTTCAGCAATATCAACCCCGCCATGCCGCCCGCAAAGGCCAGCAGCGGCAGCACCATCACGGGCAGTGACGGCATAAGAAGTAGAGCCCCCCACAGAGGCCAGGCTGGCGGCATGGTTAACACCGAGAATATCCGGTGa                                                                                 <  1:555457/140‑1 (MQ=255)                                                  tCTTCAGCAATATCAACCCCGCCATGCCGCCCGCAAAGGCCAGCAGCGGCAGCACCATCACGGGCAGTGACGGCATAAGAAGTAGAGCCCCCCACAGAGGCCAGGCTGGCGGCATGGTTAACACCGAGAATATCCGGTGa                                                                                 >  1:728992/1‑140 (MQ=255)                                                      cAGCAATATCAACCCCGCCATGCCGCCCGCAAAGGCCAGCAGCGGCAGCACCATCACGGGCAGTGACGGCATAAGAAGTAGAGCCCCCCACAGAGGCCAGGCTGGCGGCATGGTTAACACCGAGAATATCCGGTg                                                                                  <  2:469279/135‑1 (MQ=255)                                                      cAGCAATATCAACCCCGCCATGCCGCCCGCAAAGGCCAGCAGCGGCAGCACCATCACGGGCAGTGACGGCATAAGAAGTAGAGCCCCCCACAGAGGCCAGGCTGGCGGCATGGTTAACACCGAGAATATCCGGTg                                                                                  >  1:469279/1‑135 (MQ=255)                                                            tatCAACCCCGCCATGCCGCCCGCAAAGGCCAGCAGCGGCAGCACCATCACGGGCAGTGACGGCATAAGAAGTAGAGCCCCCCACAGAGGCCAGGCTGGCGGCATGGTTAACACCGAGAATATCCGGTGATGCCAGAggg                                                                       >  1:486016/1‑140 (MQ=255)                                                            tatCAACCCCGCCATGCCGCCCGCAAAGGCCAGCAGCGGCAGCACCATCACGGGCAGTGACGGCATAAGAAGTAGAGCCCCCCACAGAGGCCAGGCTGGCGGCATGGTTAACACCGAGAATATCCGGTGATGCCAGAggg                                                                       >  2:754212/1‑140 (MQ=255)                                                            tatCAACCCCGCCATGCCGCCCGCAAAGGCCAGCAGCGGCAGCACCATCACGGGCAGTGACGGCATAAGAAGTAGAGCCCCCCACAGAGGCCAGGCTGGCGGCATGGTTAACACCGAGAATATCCGGTGATGCCAGAggg                                                                       >  1:36281/1‑140 (MQ=255)                                                            tatCAACCCCGCCATGCCGCCCGCAAAGGCCAGCAGCGGCAGCACCATCACGGGCAGTGACGGCATAAGAAGTAGAGCCCCCCACAGAGGCCAGGCTGGCGGCATGGTTAACACCGAGAATATCCGGTGATGCCAGAggg                                                                       >  1:130660/1‑140 (MQ=255)                                                                 aCCCCGCCATGCCGCACGCAAAGGCCAGCAGCGGCAGCACCATCACGGGCAGTGACGGCATAAGAAGTAGAGCCCCCCACAGAGGCCAGGCTGGCGGCATGGTTAACACCGAGAATATCCGGTGATGCCAGAGGGTTgcg                                                                  >  2:740544/1‑140 (MQ=255)                                                                    ccGCCATGCCGCCCGCAAAGGCCAGCAGCGGCAGCACCATCACGGGCAGTGACGGCATAAGAAGTAGAGCCCCCCACAGAGGCCAGGCTGGCGGCATGGTTAACACCGAGAATATCCGGTGa                                                                                 >  1:65379/1‑122 (MQ=255)                                                                    ccGCCATGCCGCCCGCAAAGGCCAGCAGCGGCAGCACCATCACGGGCAGTGACGGCATAAGAAGTAGAGCCCCCCACAGAGGCCAGGCTGGCGGCATGGTTAACACCGAGAATATCCGGTGa                                                                                 <  2:65379/122‑1 (MQ=255)                                                                      gCCATGCCGCCCGCAAAGGCCAGCAGCGGCAGCACCATCACGGGCAGTGACGGCATAAGAAGTAGAGCCCCCCACAGAGGCCAGGCTGGCGGCATGGTTAACACCGAGAATATCCGGTGATGCCAGAGGGTTGCGCACaa                                                             <  2:758843/140‑1 (MQ=255)                                                                      gCCATGCCGCCCGCAAAGGCCAGCAGCGGCAGCACCATCACGGGCAGTGACGGCATAAGAAGTAGAGCCCCCCACAGAGGCCAGGCTGGCGGCATGGTTAACACCGAGAATATCCGGTGATGCCAGAGGGTTGCGCACaa                                                             >  1:190045/1‑140 (MQ=255)                                                                      gCCATGCCGCCCGCAAAGGCCAGCAGCGGCAGCACCATCACGGGCAGTGACGGCATAAGAAGTAGAGCCCCCCACAGAGGCCAGGCTGGCGGCATGGTTAACACCGAGAATATCCGGTGATGCCAGAGGGTTGCGCACaa                                                             <  2:191047/140‑1 (MQ=255)                                                                       ccATGCCGCCCGCAAAGGCCAGCAGCGGCAGCACCATCACGGGCAGTGACGGCATAAGAAGTAGAGCCCCCCACAGAGGCCAGGCTGGCGGCATGGTTAACACCGAGAATATCCGGTGATGCCAGAGGGTTGCGCACAAt                                                            >  1:213422/1‑140 (MQ=255)                                                                          tGCCGCCCGCAAAGGCCAGCAGCGGCAGCACCATCACGGGCAGTGACGGCATAAGAAGTAGAGCCCCCCACAGAGGCCAGGCTGGCGGCATGGTTAACACCGAGAATATCCGGTGATGCCAGAGGGTTGCGCACAATccc                                                         <  2:213422/140‑1 (MQ=255)                                                                          tGCCGCCCGCAAAGGCCAGCAGCGGCAGCACCATCACGGGCAGTGACGGCATAAGAAGTAGAGCCCCCCACAGAGGCCAGGCTGGCGGCATGGTTAACACCGAGAATATCCGGTGATGCCAGAGGGTTGCGCACAATccc                                                         >  2:312122/1‑140 (MQ=255)                                                                          tGCCGCCCGCAAAGGCCAGCAGCGGCAGCACCATCACGGGCAGTGACGGCATAAGAAGTAGAGCCCCCCACAGAGGCCAGGCTGGCGGCATGGTTAACACCGAGAATATCCGGTGATGCCAGAGGGTTGCGCACAATccc                                                         >  2:934279/1‑140 (MQ=255)                                                                                       ggCCAGCAGCGGCAGCACCATCACGGGCAGTGACGGCATAAGAAGTAGAGCCCCCCACAGAGGCCAGGCTGGCGGCATGGTTAACACCGAGAATATCCGGTGATGCCAGAGGGTTGCGCACAACCCCCTGTATCAGCAcg                                            <  2:130660/140‑1 (MQ=255)                                                                                          cagcagCGGCAGCACCATCACGGGCAGTGACGGCATAAGAAGTAGAGCCCCCCACAGAGGCCAGGCTGGCGGCATGGTTAACACCGAGAATATCCGGTGATGCCAGAGGGTTGCGCACAATCCCCTGTATCAGCAcgccc                                         >  1:912408/1‑140 (MQ=255)                                                                                             cagcGGCAGCACCATCACGGGCAGTGACGGCATAAGAAGTAGAGCCCCCCACAGAGGCCAGGCTGGCGGCATGGTTAACACCGAGAATATCCGGTGATGCCAGAGGGTTGCGCACAATCCCCTGTATCAGCAcgcccgcc                                      <  1:12000/140‑1 (MQ=255)                                                                                                  gcagcaCCATCACGGGCAGTGACGGCATAAGAAGTAGAGCCCCCCACAGAGGCCAGGCTGGCGGCATGGTTAACACCGAGAATATCCGGTGATGCCAGAGGGTTGCGCACAATCCCCTGTATCAGCACGCCCGCCACGGc                                 <  2:398714/140‑1 (MQ=255)                                                                                                  gcagcaCCATCACGGGCAGTGACGGCATAAGAAGTAGAGCCCCCCACAGAGGCCAGGCTGGCGGCATGGTTAACACCGAGAATATCCGGTGATGCCAGAGGGTTGCGCACAATCCCCTGTATCAGCACGCCCGCCACGGc                                 >  2:227629/1‑140 (MQ=255)                                                                                                  gcagcaCCATCACGGGCAGTGACGGCATAAGAAGTAGAGCCCCCCACAGAGGCCAGGCTGGCGGCATGGTTAACACCGAGAATATCCGGTGATGCCAGAGGGTTGCGCACAATCCCCTGTATCAGCACGCCCGCCACGGc                                 >  1:1137927/1‑140 (MQ=255)                                                                                                  gcagcaCCATCACGGGCAGTGACGGCATAAGAAGTAGAGCCCCCCACAGAGGCCAGGCTGGCGGCATGGTTAACACCGAGAATATCCGGTGATGCCAGAGGGTTGCGCACAATCCCCTGTATCAGCACGCCCGCCACGGc                                 >  1:270297/1‑140 (MQ=255)                                                                                                          aTCACGGGCAGTGACGGCATAAGAAGTAGAGCCCCCCACAGAGGCCAGGCTGGCGGCATGGTTAACACCGAGAATATCCGGTGATGCCAGAGGGTTGCGCACAATCCCCTGTATCAGCACGCCCGCCACGGCGAGGGCTg                         <  2:453118/140‑1 (MQ=255)                                                                                                          aTCACGGGCAGTGACGGCATAAGAAGTAGAGCCCCCCACAGAGGCCAGGCTGGCGGCATGGTTAACACCGAGAATATCCGGTGATGCCAGAGGGTTGCGCACAATCCCCTGTATCAGCACGCCCGCCACGGCGAGGGCTg                         >  1:246759/1‑140 (MQ=255)                                                                                                                   aGTGACGGCATAAGAAGTAGAGCCCCCCACAGAGGCCAGGCTGGCGGCATGGTTAACACCGAGAATATCCGGTGATGCCAGAGGGTTGCGCACAATCCCCTGTATCAGCACGCCCGCCACGGCGAGGGCTGCACCGACaa                >  1:747181/1‑140 (MQ=255)                                                                                                                   aGTGACGGCATAAGAAGTAGAGCCCCCCACAGAGGCCAGGCTGGCGGCATGGTTAACACCGAGAATATCCGGTGATGCCAGAGGGTTGCGCACAATCCCCTGTATCAGCACGCCCGCCACGGCGAGGACTGCACCGACaa                >  2:440100/1‑140 (MQ=255)                                                                                                                           cATAAGAAGTAGAGCCCCCCACAGAGGCCAGGCTGGCGGCATGGTTAACACCGAGAATATCCGGTGATGCCAGAGGGTTGCGCACAATCCCCTGTATCAGCACGCCCGCCACGGCGAGGGCTGCACCGACAAACAGTGcc        >  2:151725/1‑140 (MQ=255)                                                                                                                           cATAAGAAGTAGAGCCCCCCACAGAGGCCAGGCTGGCGGCATGGTTAACACCGAGAATATCCGGTGATGCCAGAGGGTTGCGCACAATCCCCTGTATCAGCACGCCCGCCACGGCGAGGGCTGCACCGACAAACAGTGcc        >  2:752801/1‑140 (MQ=255)                                                                                                                           cATAAGAAGTAGAGCCCCCCACAGAGGCCAGGCTGGCGGCATGGTTAACACCGAGAATATCCGGTGATGCCAGAGGGTTGCGCACAATCCCCTGTATCAGCACGCCCGCCACGGCGAGGGCTGCACCGACAAACAGTGcc        >  1:421968/1‑140 (MQ=255)                                                                                                                                 aagTAGAGCCCCCCACAGAGGCCAGGCTGGCGGCATGGTTAACACCGAGAATATCCGGTGATGCCAGAGGGTTGCGCACAATCCCCTGTATCAGCACGCCCGCCACGGCGAGGGCTGCACCGACAAACAGTGCCagcagc  <  2:728992/140‑1 (MQ=255)                                                                                                                                         |                                                                                                                                    GCCACGCCGGTGAGCGCCAGCTTCATCGGCTGGTGGGTCTTTGCCAGCATCTTCAGCAATATCAACCCCGCCATGCCGCCCGCAAAGGCCAGCAGCGGCAGCACCATCACGGGCAGTGACGGCATAAGAAGTAGAG‑CCCCCACAGAGGCCAGGCTGGCGGCATGGTTAACACCGAGAATATCCGGTGATGCCAGAGGGTTGCGCACAATCCCCTGTATCAGCACGCCCGCCACGGCGAGGGCTGCACCGACAAACAGTGCCAGCAGC  >  NC_000913/4511984‑4512250

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