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_0009134,511,1460CT100.0% 211.6 / NA 65W104* (TGG→TGA) fecEferric citrate ABC transporter ATPase
Reads supporting (aligned to +/- strand):  ref base C (0/0);  new base T (34/31);  total (34/31)

AGAAATGCGCGCTGGCGCTGACCGCCGGAAAGCTCGGTTAACCGACGAACGGCAAGATGATTGATCCGGGTCTGGTTCATGGCGACATTAACTCGTGCATTGTCTTCAGCGGAGAGACGCCCCCAGAGTGACAGCCAGGGATTACGACCATACGAAACCAGCTCCTGGACTGTGATCCCCTCTGGCGTTAAATGGTGCTGAGGCAGCAGCGAAAGCCTGCGGGCCAACTGGCGCGATGAGAGCATATTTATGGGATTATCG  >  NC_000913/4511012‑4511272
                                                                                                                                      |                                                                                                                              
aGAAATGCGCGCTGGCGCTGACCGCCGGAAAGCTCGGTTAACCGACGAACGGCAAGATGATTGATCCGGGTCTGGTTCATGGCGACATTAACTCGTGCATTGTCTTCAGCGGAGAGACGCCCCCAGAGTGACAGTCAgg                                                                                                                            >  2:31359/1‑139 (MQ=255)
  aaaTGCGCGCTGGCGCTGACCGCCGGAAAGCTCGGTTAACCGACGAACGGCAAGATGATTGATCCGGGTCTGGTTCATGGCGACATTAACTCGTGCATTGTCTTCAGCGGAGAGACGCCCCCAGAGTGACAGTCAGGGa                                                                                                                          >  1:1410356/1‑139 (MQ=255)
       cgcgcTGGCGCTGACCGCCGGAAAGCTCGGTTAACCGACGAACGGCAAGATGATTGATCCGGGTCTGGTTCATGGCGACATTAACTCGTGCATTGTCTTCAGCGGAGAGACGCCCCCAGAGTGACAGTCAGGGATTACg                                                                                                                     >  2:1584038/1‑139 (MQ=255)
         cgttggcgctgACCGCCGGAAAGCTCGGTTAACCGACGAACGGCAAGATGATTGATCCGGGTCTGGTTCATGGCGACATTAACTCGTGCATTGTCTTCAGCGGAGAGACGCCCCCAGAGTGACAGTCAGGGATTACGAc                                                                                                                   <  2:256361/136‑1 (MQ=255)
            tggcgctgACCGCCGGAAAGCTCGGTTAACCGACGAACGGCAAGATGATTGATCCGGGTCTGGTTCATGGCGACATTAACTCGTGCATTGTCTTCAGCGGAGAGACGCCCCCAGAGTGACAGTCa                                                                                                                              >  1:316743/1‑125 (MQ=255)
                    aCCGCCGGAAAGCTCGGTTAACCGACGAACGGCAAGATGATTGATCCGGGTCTGGTTCATGGCGACATTAACTCGTGCATTGTCTTCAGCGGAGAGACGCCCCCAGAGTGACAGTCAGGGATTACGACCATACGAAAcc                                                                                                        <  2:1723629/139‑1 (MQ=255)
                    aCCGCCGGAAAGCTCGGTTAACCGACGAACGGCAAGATGATTGATCCGGGTCTGGTTCATGGCGACATTAACTCGTGCATTGTCTTCAACGGAGAGACGCCCCCAGAGTGACAGTCAGGGATTACGACCATACGAAAcc                                                                                                        <  2:1039360/139‑1 (MQ=255)
                     ccgccgGAAAGCTCGGTTAACCGACGAACGGCAAGATGATTGATCCGGGTCTGGTTCATGGCGACATTAACTCGTGCATTGTCTTCAGCGGAGAGACGCCCCCAGAGTGACAGTCAGGGATTACGACCATACGAAACCa                                                                                                       <  1:1482362/139‑1 (MQ=255)
                      cgccgGAAAGCTCGGTTAACCGACGAACGGCAAGATGATTGATCCGGGTCTGGTTCATGGCGACATTAACTCGTGCATTGTCTTCAGCGGAGAGACGCCCCCAGAGTGACAGTCa                                                                                                                              >  1:1092356/1‑115 (MQ=255)
                      cgccgGAAAGCTCGGTTAACCGACGAACGGCAAGATGATTGATCCGGGTCTGGTTCATGGCGACATTAACTCGTGCATTGTCTTCAGCGGAGAGACGCCCCCAGAGTGACAGTCAGGGATTACGACCATACGAAACCAg                                                                                                      <  2:1410356/139‑1 (MQ=255)
                        ccgGAAAGCTCGGTTAACCGACGAACGGCAAGATGATTGATCCGGGTCTGGTTCATGGCGACATTAACTCGTGCATTGTCTTCAGCGGAGAGACGCCCCCAGAGTGACAGTCa                                                                                                                              >  2:80604/1‑113 (MQ=255)
                        ccgGAAAGCTCGGTTAACCGACGAACGGCAAGATGATTGATCCGGGTCTGGTTCATGGCGACATTAACTCGTGCATTGTCTTCAGCGGAGAGACGCCCCCAGAGTGACAGTCa                                                                                                                              <  1:80604/113‑1 (MQ=255)
                           gAAAGCTCGGTTAACCGACGAACGGCAAGATGATTGATCCGGGTCTGGTTCATGGCGACATTAACTCGTGCATTGTCTTCAGCGGAGAGACGCCCCCAGAGTGACAGTCAGGGATTACGACCATACGAAACCAGCTCCt                                                                                                 <  1:31359/139‑1 (MQ=255)
                              aGCTCGGTTAACCGACGAACGGCAAGATGATTGATCCGGGTCTGGTTCATGGCGACATTAACTCGTGCATTGTCTTCAGCGGAGAGACGCCCCCAGAGTGACAGTCAGGGATTACGACCATACGAAACCAGCTCCTGGa                                                                                              <  1:449047/139‑1 (MQ=255)
                                   ggTTAACCGACGAACGGCAAGATGATTGATCCGGGTCTGGTTCATGGCGACATTAACTCGTGCATTGTCTTCAGCGGAGAGACGCCCCCAGAGTGACAGTCAGGGATTACGACCATACGAAACCAGCTCCTGGACtgtg                                                                                         >  2:202982/1‑139 (MQ=255)
                                        aCCGACGAACGGCAAGATGATTGATCCGGGTCTGGTTCATGGCGACATTAACTCGTGCATTGTCTTCAGCGGAGAGACGCCCCCAGAGTGACAGTCAGGGATTACGACCATACGAAAcc                                                                                                        <  1:268315/119‑1 (MQ=255)
                                        aCCGACGAACGGCAAGATGATTGATCCGGGTCTGGTTCATGGCGACATTAACTCGTGCATTGTCTTCAGCGGAGAGACGCCCCCAGAGTGACAGTCAGGGATTACGACCATACGAAAcc                                                                                                        >  2:268315/1‑119 (MQ=255)
                                        aCCGACGAACGGCAAGATGATTGATCCGGGTCTGGTTCATGGCGACATTAACTCGTGCATTGTCTTCAGCGGAGAGACGCCCCCAGAGTGACAGTCAGGGATTACGACCATACGAAACCAGCTCCTGGACTGTGATccc                                                                                    >  2:1726183/1‑139 (MQ=255)
                                               aaCGGCAAGATGATTGATCCGGGTCTGGTTCATGGCGACATTAACTCGTGCATTGTCTTCAGCGGAGAGACGCCCCCAGAGTGACAGTCAGGGATTACGACCATACGAAACCAGCTCCTGGACTGTGATCCCCTCTGGc                                                                             >  2:1270046/1‑139 (MQ=255)
                                                 cGGCAAGATGATTGATCCGGGTCTGGTTCATGGCGACATTAACTCGTGCATTGTCTTCAGCGGAGAGACGCCCCCAGAGTGACAGTCa                                                                                                                              >  2:550652/1‑88 (MQ=255)
                                                     aaGATGATTGATCCGGGTCTGGTTCATGGCGACATTAACTCGTGCATTGTCTTCAGCGGAGAGACGCCCCCAGAGTGACAGTCAGGGATTACGACCATACGAAACCAGCTCCTGGACTGTGATCCCCTCTGGCGTTaaa                                                                       <  1:1726183/139‑1 (MQ=255)
                                                                 ccGGGTCTGGTTCATGGCGACATTAACTCGTGCATTGTCTTCAGCGGAGAGACGCCCCCAGAGTGACAGTCAGGGATTACGACCATACGAAACCAGCTCCTGGACTGTGATCCCCTCTGGCGTTAAATGGTGCTGAGgc                                                           <  2:659204/139‑1 (MQ=255)
                                                                    ggTCTGGTTCATGGCGACATTAACTCGTGCATTGTCTTCAGCGGAGAGACGCCCCCAGAGTGACAGTCAGGGATTACGACCATACGAAACCAGCTCCTGGACTGTGATCCCCTCTGGCGTTAAATGGTGCTGAGgcagc                                                        >  1:470272/1‑139 (MQ=255)
                                                                     gTCTGGTTCATGGCGACATTAACTCGTGCATTGTCTTCAGCGGAGAGACGCCCCCAGAGTGACAGTCAGGGATTACGACCATACGAAACCAGCTCCTGGACTGTGATCCCCTCTGGCGTTAAATGGTGCTGAGgcagca                                                       >  1:993443/1‑139 (MQ=255)
                                                                          gTTCATGGCGACATTAACTCGTGCATTGTCTTCAGCGGAGAGACGCCCCCAGAGTGACAGTCAGGGATTACGACCATACGAAACCAGCTCCTGGa                                                                                              <  2:1510624/95‑1 (MQ=255)
                                                                          gTTCATGGCGACATTAACTCGTGCATTGTCTTCAGCGGAGAGACGCCCCCAGAGTGACAGTCAGGGATTACGACCATACGAAACCAGCTCCTGGa                                                                                              >  1:1510624/1‑95 (MQ=255)
                                                                          gTTCATGGCGACATTAACTCGTGCATTGTCTTCAGCGGAGAGACGCCCCCAGAGTGACAGTCAGGGATTACGACCATACGAAACCAGCTCCTGGACTGTGATCCCCTCTGGCGTTAAATGGTGCTGAGGCAGCAGCGaa                                                  >  2:630649/1‑139 (MQ=255)
                                                                          gTTCATGGCGACATTAACTCGTGCATTGTCTTCAGCGGAGAGACGCCCCCAGAGTGACAGTCAAGGATTACGACCATACGAAACCAGCTCCTGGAATGTGATCCCCTCTGGCGTTAAATGGTGCTGAGGCAGCAGCGaa                                                  >  2:1126152/1‑139 (MQ=255)
                                                                               tGGCGACATTAACTCGTGCATTGTCTTCAGCGGAGAGACGCCCCCAGAGTGACAGTCAGGGATTACGACCATACGAAACCAGCTCCTGGACTGTGATCCCCTCTGGCGTTAAATGGTGCTGAGGCAGCAGCGAAAGCCt                                             <  1:1270046/139‑1 (MQ=255)
                                                                                 gCGACATTAACTCGTGCATTGTCTTCAGCGGAGAGACGCCCCCAGAGTGACAGTCAGGGATTACGACCATACGAAACCAGCTCCTGGACTGTGATCCCCTCTGGCGTTAAATGGTGCTGAGGCAGCAGCGAAAGCCTGc                                           <  2:993443/139‑1 (MQ=255)
                                                                                               tGCATTGTCTTCAGCGGAGAGACGCCCCCAGAGTGACAGTCAGGGATTACGACCATACGAAAcc                                                                                                        <  2:586150/64‑1 (MQ=255)
                                                                                               tGCATTGTCTTCAGCGGAGAGACGCCCCCAGAGTGACAGTCAGGGATTACGACCATACGAAAcc                                                                                                        >  1:586150/1‑64 (MQ=255)
                                                                                                gCATTGTCTTCAGCGGAGAGACGCCCCCAGAGTGACAGTCAGGGATTACGACCATACGAAACCAGCTCCTGGACTGTGATCCCCTCTg                                                                               <  1:622391/88‑1 (MQ=255)
                                                                                                gCATTGTCTTCAGCGGAGAGACGCCCCCAGAGTGACAGTCAGGGATTACGACCATACGAAACCAGCTCCTGGACTGTGATCCCCTCTg                                                                               >  2:622391/1‑88 (MQ=255)
                                                                                                  aTTGTCTTCAGCGGAGAGACGCCCCCAGAGTGACAGTCAGGGATTACGACCATACGAAACCAGCTCCTGGACTGTGATCCCCTCTGGCGTTAAATGGTGCTGAGGCAGCAGCGAAAGCCTGCGGGCCAACTGGCGCGAt                          <  1:630649/139‑1 (MQ=255)
                                                                                                    tGTCTTCAGCGGAGAGACGCCCCCAGAGTGACAGTCAGGGATTACGACCATACGAAAcc                                                                                                        <  2:446421/59‑1 (MQ=255)
                                                                                                    tGTCTTCAGCGGAGAGACGCCCCCAGAGTGACAGTCAGGGATTACGACCATACGAAAcc                                                                                                        >  1:446421/1‑59 (MQ=255)
                                                                                                         tCAGCGGAGAGACGCCCCCAGAGTGACAGTCAGGGATTACGACCATACGAAACCAGCTCCTGGACTGTGATCCCCTCTGGCGTTAAATGGTGCTGAGGCAGCAGCGAAAGCCTGCGGGCCAACTGGCGCGATGagag                     <  2:976362/137‑1 (MQ=255)
                                                                                                         tCAGCGGAGAGACGCCCCCAGAGTGACAGTCAGGGATTACGACCATACGAAACCAGCTCCTGGACTGTGATCCCCTCTGGCGTTAAATGGTGCTGAGGCAGCAGCGAAAGCCTGCGGGCCAACTGGCGCGATGagag                     >  1:976362/1‑137 (MQ=255)
                                                                                                          cAGCGGAGAGACGCCCCCAGAGTGACAGTCAGGGATTACGACCATACGAAACCAGCTCCTGGACtg                                                                                           <  2:885378/66‑1 (MQ=255)
                                                                                                          cAGCGGAGAGACGCCCCCAGAGTGACAGTCAGGGATTACGACCATACGAAACCAGCTCCTGGACtg                                                                                           >  1:885378/1‑66 (MQ=255)
                                                                                                             cGGAGAGACGCCCCCAGAGTGACAGTCAGGGATTACGAc                                                                                                                   >  2:1209569/1‑39 (MQ=255)
                                                                                                             cGGAGAGACGCCCCCAGAGTGACAGTCAGGGATTACGAc                                                                                                                   <  1:1209569/39‑1 (MQ=255)
                                                                                                                agagaCGCCCCCAGAGTGACAGTCAGGGATTACGACCATACGAAAcc                                                                                                        <  2:461552/47‑1 (MQ=255)
                                                                                                                agagaCGCCCCCAGAGTGACAGTCAGGGATTACGACCATACGAAAcc                                                                                                        >  1:461552/1‑47 (MQ=255)
                                                                                                                agagaCGCCCCCAGAGTGACAGTCAGGGATTACGACCATACGAAACCAGCTCCTGGACTGTGATCCCCTCTGGCGTTAAATGGTGCTGAGGCAGCAGCGAAAGCCTGCGGGCCAACTGGCGCGATGAGAGCATATTTAt            >  1:145331/1‑139 (MQ=255)
                                                                                                                agagaCGCCCCCAGAGTGACAGTCAGGGATTACGACCATACGAAACCAGCTCCTGGACTGTGATCCCCTCTGGCGTTAAATGGTGCTGAGGCAGCAGCGAAAGCCTGCGGGCCAACTGGCGCGATGAGAGCATATTTAt            >  1:760597/1‑139 (MQ=255)
                                                                                                                  agaCGCCCCCAGAGTGACAGTCAGGGATTACGACCATACGAAACCAGCTc                                                                                                   >  2:29386/1‑50 (MQ=255)
                                                                                                                  agaCGCCCCCAGAGTGACAGTCAGGGATTACGACCATACGAAACCAGCTc                                                                                                   <  1:29386/50‑1 (MQ=255)
                                                                                                                  agaCGCCCCCAGAGTGACAGTCAGGGATTACGACCATACGAAACCAGCTCCTGGACTGTGATCCCCTCTGGCGTTAAATGGTGCTGAGGCAGCAGCGAAAGCCTGCGGGCCaa                                    >  2:1061241/1‑113 (MQ=255)
                                                                                                                  agaCGCCCCCAGAGTGACAGTCAGGGATTACGACCATACGAAACCAGCTCCTGGACTGTGATCCCCTCTGGCGTTAAATGGTGCTGAGGCAGCAGCGAAAGCCTGCGGGCCaa                                    <  1:1061241/113‑1 (MQ=255)
                                                                                                                  agaCGCCCCCAGAGTGACAGTCAGGGATTACGACCATACGAAACCAGCTCCTGGACTGTGATCCCCTCTGGCGTTAAATGGTGCTGAGGCAGCAGCGAAAGCCTGCGGGCCAACTGGCGCGATGAGAGCATATTTATgg          >  1:830199/1‑139 (MQ=255)
                                                                                                                       cccccAGAGTGACAGTCAGGGATTACGACCATACGAAACCAGCTCCTGGACtgt                                                                                          <  2:389558/54‑1 (MQ=255)
                                                                                                                       cccccAGAGTGACAGTCAGGGATTACGACCATACGAAACCAGCTCCTGGACtgt                                                                                          >  1:389558/1‑54 (MQ=255)
                                                                                                                          ccAGAGTGACAGTCAGGGATTACGACCATACGAAACCAGCTCCTGGACTGTGATCCCCTCTGGCGTTAAATGGTGCTGAGGCAGCAGCGAAAGCCTGCGGGCCAACTGGCGCGATGAGAGCATATTTATGGGATTATCg  <  2:470272/139‑1 (MQ=255)
                                                                                                                              agTGACAGTCAGGGATTACGACCATACGAAACCAGCTCCt                                                                                                 >  2:1404032/1‑40 (MQ=255)
                                                                                                                              agTGACAGTCAGGGATTACGACCATACGAAACCAGCTCCt                                                                                                 <  1:1404032/40‑1 (MQ=255)
                                                                                                                                  aCAGTCAGGGATTACGACCATACGAAACCAGCTCCTGGa                                                                                              >  1:1507888/1‑39 (MQ=255)
                                                                                                                                  aCAGTCAGGGATTACGACCATACGAAACCAGCTCCTGGa                                                                                              <  2:1507888/39‑1 (MQ=255)
                                                                                                                                  aCAGTCAGGGATTACGACCATACGAAACCAGCTCCTGGACt                                                                                            >  2:608145/1‑41 (MQ=255)
                                                                                                                                  aCAGTCAGGGATTACGACCATACGAAACCAGCTCCTGGACt                                                                                            <  1:608145/41‑1 (MQ=255)
                                                                                                                                  aCAGTCAGGGATTACGACCATACGAAACCAGCTCCTGGACTGTGATCCCCTCTGGCGTTAAATGGTGCTGAGGCAGCAGCGAAAGCCTGCGGGCCAACTGGCGCGAt                          <  2:1160266/107‑1 (MQ=255)
                                                                                                                                  aCAGTCAGGGATTACGACCATACGAAACCAGCTCCTGGACTGTGATCCCCTCTGGCGTTAAATGGTGCTGAGGCAGCAGCGAAAGCCTGCGGGCCAACTGGCGCGATGAGAGCa                   <  1:1216995/114‑1 (MQ=255)
                                                                                                                                  aCAGTCAGGGATTACGACCATACGAAACCAGCTCCTGGACTGTGATCCCCTCTGGCGTTAAATGGTGCTGAGGCAGCAGCGAAAGCCTGCGGGCCAACTGGCGCGATGAGAGCa                   >  2:1216995/1‑114 (MQ=255)
                                                                                                                                  aCAGTCAGGGATTACGACCATACGAAACCAGCTCCTGCACTGTAATCCCCTCTGGCGTTAAATGGTGCTGAGGCAGCAGCGAAAGCCTGCGGGCCAACAGGCGCGAt                          >  1:1160266/1‑107 (MQ=255)
                                                                                                                                      |                                                                                                                              
AGAAATGCGCGCTGGCGCTGACCGCCGGAAAGCTCGGTTAACCGACGAACGGCAAGATGATTGATCCGGGTCTGGTTCATGGCGACATTAACTCGTGCATTGTCTTCAGCGGAGAGACGCCCCCAGAGTGACAGCCAGGGATTACGACCATACGAAACCAGCTCCTGGACTGTGATCCCCTCTGGCGTTAAATGGTGCTGAGGCAGCAGCGAAAGCCTGCGGGCCAACTGGCGCGATGAGAGCATATTTATGGGATTATCG  >  NC_000913/4511012‑4511272

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