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% 116.4 / NA 37W104* (TGG→TGAfecEferric citrate ABC transporter ATPase
Reads supporting (aligned to +/- strand):  ref base C (0/0);  new base T (18/19);  total (18/19)

CCAGAAATGCGCGCTGGCGCTGACCGCCGGAAAGCTCGGTTAACCGACGAACGGCAAGATGATTGATCCGGGTCTGGTTCATGGCGACATTAACTCGTGCATTGTCTTCAGCGGAGAGACGCCCCCAGAGTGACAGCCAGGGATTACGACCATACGAAACCAGCTCCTGGACTGTGATCCCCTCTGGCGTTAAATGGTGCTGAGGCAGCAGCGAAAGCCTGCGGGCCAACTGGCGCGATGAGAGCATATTTATGGGA  >  NC_000913/4511010‑4511266
                                                                                                                                        |                                                                                                                        
ccAGAAATGCGCGCTGGCGCTGACCGCCGGAAAGCTCGGTTAACCGACGAACGGCAAGATGATTGATCCGGGTCTGGTTCATGGCGACATTAACTCGTGCATTGTCTTCAGCGGAGAGACGCCCCCAGAGTGACAGTCa                                                                                                                        <  1:771183/139‑1 (MQ=255)
                              aaaGCTCGGTTAACCGACGAACGGCAAGATGATTGATCCGGGTCTGGTTCATGGCGACATTAACTCGTGCATTGTCTTCAGCGGAGAGACGCCCCCAGAGTGACAGTCAGGGATTACGACCATACGAAAACAGCTCCTg                                                                                          <  1:211828/139‑1 (MQ=255)
                                aGCTCGGTTAACCGACGAACGGCAAGATGATTGATCCGGGTCTGGTTCATGGCGACATTAACTCGTGCATTGTCTTCAGCGGAGAGACGCCCCCAGAGTGACAGTCAGGGATTACGACCATACGAAACCAGCTCCTGGa                                                                                        <  2:337384/139‑1 (MQ=255)
                                               cgaACGGCAAGATGATTGATCCGGGTCTGGTTCATGGCGACATTAACTCGTGCATTGTCTTCAGCGGAGAGACGCCCCCAGAGTGACAGTCAGGGATTACGAGCATACGAAAcc                                                                                                  >  2:858532/1‑114 (MQ=255)
                                               cgaACGGCAAGATGATTGATCCGGGTCTGGTTCATGGCGACATTAACTCGTGCATTGTCTTCAGCGGAGAGACGCCCCCAGAGTGACAGTCAGGGATTACGACCATACGAAAcc                                                                                                  <  1:858532/114‑1 (MQ=255)
                                               cgaACGGCAAGATGATTGATCCGGGTCTGGTTCATGGCGACATTAACTCGTGCATTGTCTTCAGCGGAGAGACGCCCCCAGAGTGACAGTCAGGGATTACGACCATACGAAACCAGCTCCTGGACTGTGATCCCCTCTg                                                                         >  2:54777/1‑139 (MQ=255)
                                                                    cGGGTCTGGTTCATGGCGACATTAACTCGTGCATTGTCTTCAGCGGAGAGACGCCCCCAGAGTGACAGTCAGGGATTACGACCATACGAAACCAGCTCCTGGACTGTGATCCCCTCTGGCGTTAAATGGTGCTGAGgca                                                    >  2:1308118/1‑139 (MQ=255)
                                                                            gTTCATGGCGACATTAACTCGTGCATTGTCTTCAGCGGAGAGACGCCCCCAGAGTGACAGTCAGGGATTACGACCATACg                                                                                                       <  1:1047940/80‑1 (MQ=255)
                                                                            gTTCATGGCGACATTAACTCGTGCATTGTCTTCAGCGGAGAGACGCCCCCAGAGTGACAGTCAGGGATTACGACCATACg                                                                                                       >  2:1047940/1‑80 (MQ=255)
                                                                                 tGGCGACATTAACTCGTGCATTGTCTTCAGCGGAGAGACGCCCCCAGAGTGACAGTCAGGGATTACGACCATACGAAACCAGCTCCTGGACTGTGATCCCCTCTGGCGTTAAATGGTGCTGAGGCAGCAGCGAAAGCCt                                       >  1:471464/1‑139 (MQ=255)
                                                                                      aCATTAACTCGTGCATTGTCTTCAGCGGAGAGACGCCCCCAGAGTGACAGTCAGGGATTACGACCATACGAAAcc                                                                                                  <  1:219476/75‑1 (MQ=255)
                                                                                      aCATTAACTCGTGCATTGTCTTCAGCGGAGAGACGCCCCCAGAGTGACAGTCAGGGATTACGACCATACGAAAcc                                                                                                  >  2:219476/1‑75 (MQ=255)
                                                                                                 tGCATTGTCTTCAGCGGAGAGACGCCCCCAGAGTGACAGTCa                                                                                                                        >  1:403614/1‑42 (MQ=255)
                                                                                                 tGCATTGTCTTCAGCGGAGAGACGCCCCCAGAGTGACAGTCa                                                                                                                        <  2:403614/42‑1 (MQ=255)
                                                                                                 tGCATTGTCTTCAGCGGAGAGACGCCCCCAGAGTGACAGTCAGGGATTACGACCATACGAAACCAGCTCCTGGACTGTGATCCCCTCTGGCGTTaa                                                                  <  1:892374/96‑1 (MQ=255)
                                                                                                 tGCATTGTCTTCAGCGGAGAGACGCCCCCAGAGTGACAGTCAGGGATTACGACCATACGAAACCAGCTCCTGGACTGTGATCCCCTCTGGCGTTaa                                                                  >  2:892374/1‑96 (MQ=255)
                                                                                                 tGCATTGTCTTCAGCGGAGAGACGCCCCCAGAGTGACAGTCAGGGATTACGACCATACGAAACCAGCTCCTGGACTGTGATCCCCTCTGGCGTTAAATGGTGCTGAGGCAGCAGCGAAAGCCTGCGGGCCAACTGgcgc                       <  1:54777/139‑1 (MQ=255)
                                                                                                  gCATTGTCTTCAGCGGAGAGACGCCCCCAGAGTGACAGTCAGGGATTACGAc                                                                                                             >  1:544034/1‑52 (MQ=255)
                                                                                                  gCATTGTCTTCAGCGGAGAGACGCCCCCAGAGTGACAGTCAGGGATTACGAc                                                                                                             <  2:544034/52‑1 (MQ=255)
                                                                                                  gCATTGTCTTCAGCGGAGAGACGCCCCCAGAGTGACAGTCAGGGATTACGACCATACGAAACCAGCTCCTGGACTGTGATCCCCTCTGGCGTTAAATGGTGCTGAGGCAGCAGCGAAAGCCTGCGGGCCAACTGGcgcg                      <  2:471464/139‑1 (MQ=255)
                                                                                                        tCTTCAGCGGAGAGACGCCCCCAGAGTGACAGTCAGGGATTACGACCATACGAAACCAGCTCCTGGACTGTGATc                                                                                <  2:945793/75‑1 (MQ=255)
                                                                                                        tCTTCAGCGGAGAGACGCCCCCAGAGTGACAGTCAGGGATTACGACCATACGAAACCAGCTCCTGGACTGTGATc                                                                                >  1:945793/1‑75 (MQ=255)
                                                                                                            cAGCGGAGAGACGCCCCCAGAGTGACAGTCAGGGATTACGACCATACGAAACCAGCTCCTGGACTGTGATCCCCTCTGGCGTTAAATGGTGCTGAGgc                                                     <  1:86967/98‑1 (MQ=255)
                                                                                                            cAGCGGAGAGACGCCCCCAGAGTGACAGTCAGGGATTACGACCATACGAAACCAGCTCCTGGACTGTGATCCCCTCTGGCGTTAAATGGTGCTGAGgc                                                     >  2:86967/1‑98 (MQ=255)
                                                                                                             aGCGGAGAGACGCCCCCAGAGTGACAGTCAGGGATTACGACCATACGAAACCAGCTCCTGGa                                                                                        >  1:897301/1‑62 (MQ=255)
                                                                                                             aGCGGAGAGACGCCCCCAGAGTGACAGTCAGGGATTACGACCATACGAAACCAGCTCCTGGa                                                                                        <  2:897301/62‑1 (MQ=255)
                                                                                                                      aCGCCCCCAGAGTGACAGTCAGGGATTACGACCATACGAAACCAGCTCCTGGACTGTGATCCCCTCTGGCGTTAAATGGTGCTGAGGCAGCAGCGAAAGCCTGCGGGCCAACTGGCGCGATGAGAGCATATTTATGGGa  >  2:338946/1‑139 (MQ=255)
                                                                                                                         cccccAGAGTGACAGTCAGGGATTACGACCATACGAAACCAGCTCCTGGACTGTGATCCCCTCTGGCGtt                                                                    <  2:526090/70‑1 (MQ=255)
                                                                                                                         cccccAGAGTGACAGTCAGGGATTACGACCATACGAAACCAGCTCCTGGACTGTGATCCCCTCTGGCGtt                                                                    >  1:526090/1‑70 (MQ=255)
                                                                                                                              agagTGACAGTCAGGGATTACGACCATACGAAAcc                                                                                                  >  2:380810/1‑35 (MQ=255)
                                                                                                                              agagTGACAGTCAGGGATTACGACCATACGAAAcc                                                                                                  <  1:380810/35‑1 (MQ=255)
                                                                                                                              agagTGACAGTCAGGGATTACGACCATACGAAACCAGCTCCTGGACTGTGATCCCCTCTGGCGTTaaa                                                                 <  2:1347202/68‑1 (MQ=255)
                                                                                                                              agagTGACAGTCAGGGATTACGACCATACGAAACCAGCTCCTGGACTGTGATCCCCTCTGGCGTTaaa                                                                 >  1:1347202/1‑68 (MQ=255)
                                                                                                                                agTGACAGTCAGGGATTACGACCATACGCAACCAGCTCCTGGACTGTGATCCCCTCTGGCGTTAAATGGTGCTGAGGCAGCAGCGAAAGCCTGCGGGCCAACTGGCGCGAt                    >  1:884623/1‑111 (MQ=255)
                                                                                                                                agTGACAGTCAGGGATTACGACCATACGAAACCAGCTCCTGGACTGTGATCCCCTCTGGCGTTAAATGGTGCTGAGGCAGCAGCGAAAGCCTGCGGGCCAACTGGCGCGAt                    <  2:884623/111‑1 (MQ=255)
                                                                                                                                   gACAGTCAGGGATTACGACCATACGAAACCAGCTCCTGGACTGTGATCCCct                                                                            >  2:1098540/1‑52 (MQ=255)
                                                                                                                                   gACAGTCAGGGATTACGACCATACGAAACCAGCTCCTGGACTGTGATCCCct                                                                            <  1:1098540/52‑1 (MQ=255)
                                                                                                                                        |                                                                                                                        
CCAGAAATGCGCGCTGGCGCTGACCGCCGGAAAGCTCGGTTAACCGACGAACGGCAAGATGATTGATCCGGGTCTGGTTCATGGCGACATTAACTCGTGCATTGTCTTCAGCGGAGAGACGCCCCCAGAGTGACAGCCAGGGATTACGACCATACGAAACCAGCTCCTGGACTGTGATCCCCTCTGGCGTTAAATGGTGCTGAGGCAGCAGCGAAAGCCTGCGGGCCAACTGGCGCGATGAGAGCATATTTATGGGA  >  NC_000913/4511010‑4511266

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