Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	NC_000913	4,511,146	C→T	W104* (TGG→TGA)	fecE ←	ferric citrate ABC transporter ATPase

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	4,511,146	0	C	T	100.0%	177.2 / NA	56	W104* (TGG→TGA)	fecE	ferric citrate ABC transporter ATPase
Reads supporting (aligned to +/- strand):  ref base C (0/0);  new base T (28/28);  total (28/28)

CCAGAAATGCGCGCTGGCGCTGACCGCCGGAAAGCTCGGTTAACCGACGAACGGCAAGATGATTGATCCGGGTCTGGTTCATGGCGACATTAACTCGTGCATTGTCTTCAGCGGAGAGACGCCCCCAGAGTGACAGCCAGGGATTACGACCATACGAAACCAGCTCCTGGACTGTGATCCCCTCTGGCGTTAAATGGTGCTGAGGCAGCAGCGAAAGCCTGCGGGCCAACTGGCGCGATGAGAGCATATTTATGGGATTATCGCCGAGAAA  >  NC_000913/4511010‑4511280                                                                                                                                         |                                                                                                                                       ccAGAAATGCGCGCTGGCGCTGACCGCCGGAAAGCTCGGTTAACCGACGAACGGCAAGATGATTGCTCCGGGTCTGGTTCATGGCGACATTAACTCGTGCATTGTCTTCAGCGGAGAGACGCCCCCAGAGTGACAGTCa                                                                                                                                      <  2:289776/139‑1 (MQ=255)   aGAAATGCGCGCTGGCGCTGACCGCCGGAAAGCTCGGTTAACCGACGAACGGCAAGATGATTGATCCGGGTCTGGTTCATGGCGACATTAACTCGTGCATTGTCTTCAGCGGAGAGACGCCCCCAGAGTGACAGTCAgg                                                                                                                                    >  1:1657621/1‑139 (MQ=255)                       aCCGCCGGAAAGCTCGGTTAACCGACGAACGGCAAGATGATTGATCCGGGTCTGGTTCATGGCGACATTAACTCGTGCATTGTCTTCAGCGGAGAGACGCCCCCAGAGTGACAGTCAGGGATTACGACCATACGAAAcc                                                                                                                <  1:577148/139‑1 (MQ=255)                       aCCGCCGGAAAGCTCGGTTAACCGACGAACGGCAAGATGATTGATCCGGGTCTGGTTCATGGCGACATTAACTCGTGCATTGTCTTCAGCGGAGAGACGCCCCCAGAGTGACAGTCAGGGATTACGACCATACGAAAcc                                                                                                                <  2:459035/139‑1 (MQ=255)                         cgccgGAAAGCTCGGTTACCCGACGAACGGCAAGATGATTGGTCCGGGTCTGGTTCATGGCGACATGAACTCGTGCATTGTCTTCAGCGGAGAGTCGCCCCCGGAGTGACAGTCAGGGATTACGACCATACGAAGCCAg                                                                                                              <  1:520432/139‑1 (MQ=255)                         cgccgGAAAGCTCGGTTAACCGACGAACGGCAAGATGATTGATCCGGGTCTGGTTCATGGCGACATTAACTCGTGCATTGTCTTCAGCGGAGAGACGCCCCCAGAGTGACAGTCAGGGATTACGAc                                                                                                                           <  1:535277/126‑1 (MQ=255)                         cgccgGAAAGCTCGGTTAACCGACGAACGGCAAGATGATTGATCCGGGTCTGGTTCATGGCGACATTAACTCGTGCATTGTCTTCAGCGGAGAGACGCCCCCAGAGTGACAGTCAGGGATTACGAc                                                                                                                           >  2:535277/1‑126 (MQ=255)                            cgGAAAGCTCGGTTAACCGACGAACGGCAAGATGATTGATCCGGGTCTGGTTCATGCCGACATTAACTCGTGCATTGTCTTCAGCGGAGAGACGCCCCCAGAGTGACAGTCAGGGATTACGACCATACGAAACCAGCTc                                                                                                           >  1:476023/1‑139 (MQ=255)                                   cTCGGTTAACCGACGAACGGCAA‑‑‑GATTGATCCGGGTCTGGTTCATGGCGACATTAACTCGTGCATTGTCTTCAGCGGAGAGACGCCCCCAGAGTGACAGTCAGGGATTACGACCATACGAAACCAGCTCCTGGACtgtg                                                                                                 <  2:270726/139‑1 (MQ=255)                                                cgaACGGCAAGATGATTGATCCGGGTCTGGTTCATGGCGACATTAACTCGTGCATTGTCTTCAGCGGAGAGACGCCCCCAGAGTGACAGTCAGGGATTACGACCATACGAAACCAGCTCCTGGACTGTGATCCCCTCTg                                                                                       <  2:1657621/139‑1 (MQ=255)                                                  aaCGGCAAGATGATTGATCCGGGTCTGGTTTATGGCGACATTAACTCGTGCATTGTCTTCAGCGGAGAGACGCCCCCAGAGTGACAGTCAgg                                                                                                                                    <  1:425109/92‑1 (MQ=255)                                                  aaCGGCAAGATGATTGATCCGGGTCTGGTTCATGGCGACATTAACTCGTGCATTGTCTTCAGCGGAGAGACGCCCCCAGAGTGACAGTCAgg                                                                                                                                    >  2:425109/1‑92 (MQ=255)                                                                    ccGGGTCTGGTTCATGGCGACATTAACTCGTGCATTGTCTTCAGCGGAGAGACGCCCCCAGAGTGACAGTCAGGGATTACGACCATACGAAACCAGCTCCTGGACTGTGATCCCCTCTGGCGTTAAATGGTGCTGAGgc                                                                   >  2:353716/1‑139 (MQ=255)                                                                    ccGGGTCTGGTTCATGGCGACATTAACTCGTGCATTGTCTTCAGCGGAGAGACGCCCCCAGAGTGACAGTCAGGGATCACGACCATACGAAACCAGCTCCTGGACTGTGATCCCCTCTGGCGTTAAATGGTGCTGAGgc                                                                   >  1:1055738/1‑139 (MQ=255)                                                                         tCTGGTTCATGGCGACATTAACTCGTGCATTGTCTTCAGCGGAGAGACGCCCCCAGAGTGACAGTCAGGGATTACGACCATACGAAACCAGCTc                                                                                                           >  1:1347033/1‑94 (MQ=255)                                                                         tCTGGTTCATGGCGACATTAACTCGTGCATTGTCTTCAGCGGAGAGACGCCCCCAGAGTGACAGTCAGGGATTACGACCATACGAAACCAGCTc                                                                                                           <  2:1347033/94‑1 (MQ=255)                                                                            ggTTCATGGCGACATTAACTCGTGCATTGTCTTCAGCGGAGAGACGCCCCCAGAGTGACAGTCAGGGATTACGACCATACGAAACCAGCTCCTGGACTGTGATCCCCTCTGGCGTTAAATGGTGCTGAGGCAGCAGCGa                                                           >  2:201435/1‑139 (MQ=255)                                                                             gTTCATGGCGACATTAACTCGTGCATTGTCTTCAGCGGAGAGACGCCCCCAGAGTGACAGTCAGGGATTACGACCATACGaaa                                                                                                                  <  2:850521/83‑1 (MQ=255)                                                                             gTTCATGGCGACATTAACTCGTGCATTGTCTTCAGCGGAGAGACGCCCCCAGAGTGACAGTCAGGGATTACGACCATACGaaa                                                                                                                  >  1:850521/1‑83 (MQ=255)                                                                             gTTCATGGCGACATTAACTCGTGCATTGTCTTCAGCGGAGAGACGCCCCCAGAGTGACAGTCAGGGATTACGACCATACGAAAcc                                                                                                                >  2:1065987/1‑85 (MQ=255)                                                                             gTTCATGGCGACATTAACTCGTGCATTGTCTTCAGCGGAGAGACGCCCCCAGAGTGACAGTCAGGGATTACGACCATACGAAAcc                                                                                                                <  1:1065987/85‑1 (MQ=255)                                                                             gTTCATGGCGACATTAACTCGTGCATTGTCTTCAGCGGAGAGACGCCCCCAGAGTGACAGTCAGGGATTACGACCATACGAAACCAGCTCCTGGAc                                                                                                     >  1:245095/1‑96 (MQ=255)                                                                             gTTCATGGCGACATTAACTCGTGCATTGTCTTCAGCGGAGAGACGCCCCCAGAGTGACAGTCAGGGATTACGACCATACGAAACCAGCTCCTGGAc                                                                                                     <  2:245095/96‑1 (MQ=255)                                                                             gTTCATGGCGACATTAACTCGTGCATTGTCTTCAGCGGAGAGACGCCCCCAGAGTGACAGTCAGGGATTACGACCATACGAAACCAGCTCCTGGACTGTGATCCCCTCTGGCGTTAAATGGTGCTGAGGCAGCAGCGaa                                                          >  2:1581623/1‑139 (MQ=255)                                                                                  tGGCGACATTAACTCGTGCATTGTCTTCAGCGGAGAGACGCCCCCAGAGTGACAGTCAGGGATTACGACCATACGAAACCAGCTCCTGGACTGTGATCCCCTCTGGCGTTAAATGGTGCTGAGGCAGCAGCGAAAGCCt                                                     <  2:1055738/139‑1 (MQ=255)                                                                                       aCATTAACTCGTGCATTGTCTTCAGCGGAGAGACGCCCCCAGAGTGACAGTCAGGGATTACGACCATACGAAACCAGCTCCTGGACTGTGATCCCCTCTGGCGTTAAATGGTGCTGAGGCAGCAGCGAAAGCCTGCggg                                                >  2:1374334/1‑139 (MQ=255)                                                                                       aCATTAACTCGTGCATTGTCTTCAGCGGAGAGACGCCCCCAGAGTGACAGTCAGGGATTACGACCATACGAAACCAGCTCCTGGACTGTGATCCCCTCTGGCGTTAAATGGTGCTGAGGCAGCAGCGAAAGCCTGCggg                                                >  1:652052/1‑139 (MQ=255)                                                                                                 gTGCATTGTCTTCAGCGGAGAGACGCCCCCAGAGTGACAGTCAGGGATTACGACCATACGAAACCAGCTCCTGGACTGTGATCCCCTCTGGCGTTAAATGGTGCTGAGGCAGCAGCGAAAGCCTGCGGGCCAACTGgcg                                      <  1:353716/139‑1 (MQ=255)                                                                                                  tGCATTGTCTTCAGCGGAGAGACGCCCCCAGAGTGACAGTCAGGGATTACGACCATACGAAACCAGCTCCt                                                                                                         <  2:613684/71‑1 (MQ=255)                                                                                                  tGCATTGTCTTCAGCGGAGAGACGCCCCCAGAGTGACAGTCAGGGATTACGACCATACGAAACCAGCTCCt                                                                                                         >  1:613684/1‑71 (MQ=255)                                                                                                  tGCATTGCCTTCAGCGGAGAGACGCCCCCAGAGTGACAGTCAGGGATTACGACCATACGAAACCAGCTCCTGGACTGTGATCCCCTCTGGCGTTAAATGGTGCTGAGGCAGCAGCGAAAGCCTGCGGGCCAACTGgcgc                                     >  2:1534049/1‑139 (MQ=255)                                                                                                   gCATTGTCTTCAGCGGAGAGACGCCCCCAGAGTGACAGTCAGGGCTTACGATCATACGAAACACGCTCCTGGACTGTGATCTCCTCTGCTGTCAAATGGTGGTGAGGCAGAAGCGAAAGCCTGCGGGCCAACTGtcgca                                    >  1:1164260/1‑138 (MQ=255)                                                                                                   gCATTGTCTTCAGCGGAGAGACGCCCCCAGAGTGACAGTCAGGGATTACGACCATACGAAACCAGCTCCTGGACTGTGAt                                                                                               <  2:488650/80‑1 (MQ=255)                                                                                                   gCATTGTCTTCAGCGGAGAGACGCCCCCAGAGTGACAGTCAGGGATTACGACCATACGAAACCAGCTCCTGGACTGTGAt                                                                                               >  1:488650/1‑80 (MQ=255)                                                                                                    cATTGTCTTCAGCGGAGAGACGCCCCCAGAGTGACAGTCAGGGATTACGACCATACGAAACCAGCTCCTg                                                                                                        <  2:1379056/70‑1 (MQ=255)                                                                                                    cATTGTCTTCAGCGGAGAGACGCCCCCAGAGTGACAGTCAGGGATTACGACCATACGAAACCAGCTCCTg                                                                                                        >  2:1336587/1‑70 (MQ=255)                                                                                                    cATTGTCTTCAGCGGAGAGACGCCCCCAGAGTGACAGTCAGGGATTACGACCATACGAAACCAGCTCCTg                                                                                                        >  1:1379056/1‑70 (MQ=255)                                                                                                    cATTGTCTTCAGCGGAGAGACGCCCCCAGAGTGACAGTCAGGGATTACGACCATACGAAACCAGCTCCTg                                                                                                        <  1:1336587/70‑1 (MQ=255)                                                                                                     aTTGTCTTCAGCGGAGAGACGCCCCCAGAGTGACAGTCAGGGATTACGACCATACGAAACCAGCTCCTGGACTGTGATCCCCTCTGGCGTTAAATGGTGCTGAGGCAGCAGCGAAAGCCTGCGGGCCAACTGGCGCGAt                                  <  2:652052/139‑1 (MQ=255)                                                                                                         tCTTCAGCGGAGAGACGCCCCCAGAGTGACAGTCAGGGATTACGACCATACGAAACCAGCTCCTGGACTGTGATCCCct                                                                                          >  1:415243/1‑79 (MQ=255)                                                                                                         tCTTCAGCGGAGAGACGCCCCCAGAGTGACAGTCAGGGATTACGACCATACGAAACCAGCTCCTGGACTGTGATCCCct                                                                                          <  2:415243/79‑1 (MQ=255)                                                                                                         tCTTCAGCGGAGAGACGCCCCCAGAGTGACAGTCAGGGATTACGACCATACGAAACCAGCTCCTGGACTGTGATCCCCTCTGGCGTTAAATGGTGCTGAGGCAGCAGCGAAAGCCTGCGGGCCAACTGGCGCGATgaga                              >  1:45914/1‑139 (MQ=255)                                                                                                           ttCAGCGGAGAGACGCCCCCAGAGTGACAGTCAGGGATTACGACCATACGAAAcc                                                                                                                >  2:889423/1‑55 (MQ=255)                                                                                                           ttCAGCGGAGAGACGCCCCCAGAGTGACAGTCAGGGATTACGACCATACGAAAcc                                                                                                                <  1:889423/55‑1 (MQ=255)                                                                                                           ttCAGCGGAGAGACGCCCCCAGAGTGACAGTCAGGGATTACGACCATACGAAACCAGCTCCTGGACTGTGATCCCCTCTGGCGTTAAATGGTGCTGAGGCAGCAGCGAAAGCCTGCGGGCCAACTGGCGCGATGAGAGc                            <  1:201435/139‑1 (MQ=255)                                                                                                           ttCAGCGGAGAGACGCCCCCAGAGTGACAGTCAGGGATTACGACCATACGAAACCAGCTCCTGGACTGTGATCCCCTCTGGCGTTAAATGGTGCTGAGGCAGCAGCGAAAGCCTGCGGGCCAACTGGCGCGATGAGAGc                            <  1:1534049/139‑1 (MQ=255)                                                                                                           ttCAGCGGAGAGACGCCCCCAGAGTGACAGTCAGGGATTACGACCATACGAAACCAGCTCCTGGACTGTGATCCCCTCTGGCGTTAAATGGTGCTGAGGCAGCAGCGAAAGCCTGCGGGCCAACTGGCGCGATGAGAGc                            <  2:476023/139‑1 (MQ=255)                                                                                                            tCAGCGGAGAGACGCCCCCAGAGTGACAGTCAGGGATTACGACCATACGAAACCAGCTCCTGGACTGTGATCCCCTCTGGCGTTAAATGGTGCTGAGGCAGCAGCGAAAGCCTGCGGGCCAACTGGCGCGATGAGAGCa                           <  1:1581623/139‑1 (MQ=255)                                                                                                                   agagaCGCCCCCAGAGTGACAGTCAGGGATTACGAc                                                                                                                           >  2:1010741/1‑36 (MQ=255)                                                                                                                   agagaCGCCCCCAGAGTGACAGTCAGGGATTACGAc                                                                                                                           <  1:1010741/36‑1 (MQ=255)                                                                                                                   agagaCGCCCCCAGAGTGACAGTCAGGGATTACGACCATACGAAACCATCTCCTGGACTGTGATCCCCTCTGGCGTTAAATGGTGCTGAGGCAGCAGCGAAAGCCTGCGGGCCAACTGGCGCGATGAGAGCATATTTAt                    >  2:1388659/1‑139 (MQ=255)                                                                                                                   agagaCGCCCCCAGAGTGACAGTCAGGGATTACGACCATACGAAACCAGCTCCTGGACTGTGATCCCCTCTGGCGTTAAATGGTGCTGAGGCAGCAGCGAAAGCCTg                                                    <  2:966266/107‑1 (MQ=255)                                                                                                                   agagaCGCCCCCAGAGTGACAGTCAGGGATTACGACCATACGAAACCAGCTCCTGGACTGTGATCCCCTCTGGCGTTAAATGGTGCTGAGGCAGCAGCGAAAGCCTGCGGGCCAACTGGCGCGATGAGAGCATATTTAt                    >  1:157/1‑139 (MQ=255)                                                                                                                   agagaCGCCCCCAGAGTGACAGTCAGGGAGTACGACCATACGAAACCAGCTCGTGGACTGTGATCCCCTCTGGCGTTAAATGGTGCTGAGGCAGCAGCGAAAGCCTg                                                    >  1:966266/1‑107 (MQ=255)                                                                                                                             ccAGAGTGACAGTCAGGGATTACGACCATACGAAACCAGCTCCTGGACTGTGATCCCCTCTGGCGTTAAATGGTGCTGAGGCAGCAGCGAAAGCCTGCGGGCCAACTGGCGCGATGAGAGCATATTTATGGGATTATCg          <  1:1374334/139‑1 (MQ=255)                                                                                                                                     aCAGTCAGGGATTACGACCATACGAAACCAGCTCCTGGa                                                                                                      >  2:160529/1‑39 (MQ=255)                                                                                                                                     aCAGTCAGGGATTACGACCATACGAAACCAGCTCCTGGACTGTGATCCCCTCTGGCGTTAAATGGTGCTGAGGCAGCAGCGAAAGCCTGCGGGCCAACTGGCGCGATGAGAGCATATTTATGGGATTATCGCCGAGaaa  >  2:912339/1‑139 (MQ=255)                                                                                                                                         |                                                                                                                                       CCAGAAATGCGCGCTGGCGCTGACCGCCGGAAAGCTCGGTTAACCGACGAACGGCAAGATGATTGATCCGGGTCTGGTTCATGGCGACATTAACTCGTGCATTGTCTTCAGCGGAGAGACGCCCCCAGAGTGACAGCCAGGGATTACGACCATACGAAACCAGCTCCTGGACTGTGATCCCCTCTGGCGTTAAATGGTGCTGAGGCAGCAGCGAAAGCCTGCGGGCCAACTGGCGCGATGAGAGCATATTTATGGGATTATCGCCGAGAAA  >  NC_000913/4511010‑4511280

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