Predicted mutation
evidence seq id position mutation annotation gene description
RA NC_000913 4,515,284 C→T G465D (GGC→GAC)  fecA ← TonB‑dependent outer membrane ferric citrate transporter and signal transducer, ferric citrate extracelluar receptor, FecR‑interacting protein

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009134,515,2840CT96.7% 192.5 / ‑1.0 61coding (1394/2325 nt)fecATonB‑dependent outer membrane ferric citrate transporter and signal transducer, ferric citrate extracelluar receptor, FecR‑interacting protein
Reads supporting (aligned to +/- strand):  ref base C (0/0);  major base T (29/30);  minor base . (1/1);  total (30/31)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 8.04e-01

GACTCGATATGTTCGAAACGCATACCCGGCGTGATGGTCCAGTTGCCGATGTCGATTTTGTCATCCAGATACCAGGCGTGCGCCTCGGTGCCGGAACGCGTATCGCGGTCGTAAGGGCTTGAGCCGGACGGCAACTGCCCGCTGCTGGTGGCGGTGTAGTAACGCATTTCATGCGTTGATTCATTCAAATAGCGATAGCCCACGCCCACTTCGTGCGCGGAAGGGCCGATCATAAAGATCTGGCTGTAGCG  >  NC_000913/4515161‑4515411
                                                                                                                           |                                                                                                                               
gACTCGATATGTTCGAAACGCATACCCGGCGTGATGGTCCAGTTGCCGATGTCGATTTTGTCATCCAGATACCAGGCGTGCGCCTCGGTGCCGGAACGCGTATCGCGGTCGTAAGGGCATGAGTCGGACGGCAACAGcc                                                                                                                  <  2:1111217/139‑1 (MQ=255)
 aCTCGATATGTTCGAAACGCATACCCGGCGTGATGGTCCAGTTGCCGATGTCGATTTTGTCATCCAGATACCAGGCGTGCGCCTCGGTGCCGGAACGCGTATCGCGGTCGTAAGGGCTTGAGTCGGACGGCAACTGccc                                                                                                                 >  1:1101363/1‑139 (MQ=255)
 aCTCGATATGTTCGAAACGCATACCCGGCGTGATGGTCCAGTTGCCGATGTCGATTTTGTCATCCAGATACCAGGCGTGCGCCTCGGTGCCGGAACGCGTATCGCGGTCGTAAGGGCTTGAGTCGGACGGCAACTGccc                                                                                                                 >  1:415789/1‑139 (MQ=255)
 aCTCGATATGTTCGAAACGCATACCCGGCGTGATGGTCCAGTTGCCGATGTCGATTTTGTCATCCAGATACCAGGCGTGCGCCTCGGTGCCGGAACGCGTATCGCGGTCGTAAGGGCTTGAGTCGGACGGCAACTGccc                                                                                                                 >  2:566215/1‑139 (MQ=255)
 aCTCGATATGTTCGAAACGCATACCCGGCGTGATGGTCCAGTTGCCGATGTCGATTTTGTCATCCAGATACCAGGCGTGCGCCTCGGTGCCGGAACGCGTATCGCGGTCGTAAGGGCTTGAGTCGGACGGCAACTGccc                                                                                                                 >  2:764159/1‑139 (MQ=255)
 aCTCGATATGTTCGAAACGCATACCCGGCGTGATGGTCCAGTTACCGATGTCGATTTTGTCATCCAGATACCAGGCGTGCGCCTCGGTGCCGGAACGCGTATCGCGGTCGTAAGGGCTTGAGTGGACGGCAACTGCCCg                                                                                                                >  1:166518/1‑139 (MQ=255)
  cTCGATATGTTCGAAACGCATACCCGGCGTGATGGTCCAGTTGCCGATGTCGATTTTGTCATCCAGATACCAGGCGTGCGCCTCGGTGCCGGAACGCGTATCGCGGTCGTAAGGGCTTGAGTCGGACGGCAACTGCCCg                                                                                                                <  1:1610537/139‑1 (MQ=255)
      atatGTTCGAAACGCATACCCGGCGTGATGGTCCAGTTGCCGATGTCGATTTTGTCATCCAGATACCAGGCGTGCGCCTCGGTGCCGGAACGCGTATCGCGGTCGTAAGGGCTTGAGTCGGACGGCAACTGcc                                                                                                                  <  1:914608/133‑1 (MQ=255)
      atatGTTCGAAACGCATACCCGGCGTGATGGTCCAGTTGCCGATGTCGATTTTGTCATCCAGATACCAGGCGTGCGCCTCGGTGCCGGAACGCGTATCGCGGTCGTAAGGGCTTGAGTCGGACGGCAACTGcc                                                                                                                  >  2:914608/1‑133 (MQ=255)
            tCGAAACGCATACCCGGCGTGATGGTCCAGTTGCCGATGTCGATTTTGTCATCCAGATACCAGGCGTGCGCCTCGGTGCCGGAACGCGTATCGCGGTCGTAAGGGCTTGAGTCGGACGGCAACTGCCCGCTGCtggtgg                                                                                                      <  1:84157/139‑1 (MQ=255)
                             cGTGATGGTCCAGTTGCCGATGTCGATTTTGTCATCCAGATACCAGGCGTGCGCCTCGGTGCCGGAACGCGTATCGCGGTCGTAAGGGCTTGAGTCGGACGGCAACTGCCCGCTGCTGGTGGCGGTGTAGTAACGCAtt                                                                                     >  1:402461/1‑139 (MQ=255)
                               tGATGGTCCAGTTGCCGATGTCGATTTTGTCATCCAGATACCAGGCGTGCGCCTCGGTGCCGGAACGCGTATCGCGGTCGTAAGGGCTTGAGTCGGACGGCAACTGCCCGCTGCTGGTGGCGGTGTAGTAACGCATTTc                                                                                   <  1:764159/139‑1 (MQ=255)
                               tGATGGTCCAGTTGCCGATGTCGATTTTGTCATCCAGATACCAGGCGTGCGCCTCGGTGCCGGAACGCGTATCGCGGTCGTAAGGGCTTGAGTCGGACGGCAACTGCCCGCTGCTGGTGGCGGTGTAGTAACGCATTTc                                                                                   <  2:1226099/139‑1 (MQ=255)
                                 aTGGTCCAGTTGCCGATGTCGATTTTGTCATCCAGATACCAGGCGTGCGCCTCGGTGCCGGAACGCGTATCGCGGTCGTAAGGGCTTGAGTCGGACGGCAACTGCCCGCTGCTGGTGGCGGTgtag                                                                                              <  1:1342188/126‑1 (MQ=255)
                                 aTGGTCCAGTTGCCGATGTCGATTTTGTCATCCAGATACCAGGCGTGCGCCTCGGTGCCGGAACGCGTATCGCGGTCGTAAGGGCTTGAGTCGGACGGCAACTGCCCGCTGCTGGTGGCGGTgtag                                                                                              >  2:1342188/1‑126 (MQ=255)
                                   ggTCCAGTTGCCGATGTCGATTTTGTCATCCAGATACCAGGCGTGCGCCTCGGTGCCGGAACGCGTATCGCGGTCGTAAGGGCTTGAGTCGGACGGCAACTGCCCGCTGCTGGTGGCGGTGTAGTAACGCATTTCATGc                                                                               >  2:1627160/1‑139 (MQ=255)
                                        aGTTGCCGATGTCGATTTTGTCATCCAGATACCAGGCGTACGCCTCGGTGCCGGAACGCGTATCGCCGTCGTAAGGGCTTGAGTCGGACGGCAACTGCTCGCTGCTGGTGGCGGTGTAGTAACGGATTTCATGCGTTGa                                                                          <  1:1029345/139‑1 (MQ=255)
                                            gCCGATGTCGATTTTGTCATCCAGATACCAGGCGTGCGCCTCGGTGCCGGAACGCGTATCGCGGTCGTAAGGGCTTGAGTCGGACGGCAACTGCCCGCTGCTGGTGGCgg                                                                                                   >  1:315952/1‑110 (MQ=255)
                                            gCCGATGTCGATTTTGTCATCCAGATACCAGGCGTGCGCCTCGGTGCCGGAACGCGTATCGCGGTCGTAAGGGCTTGAGTCGGACGGCAACTGCCCGCTGCTGGTGGCgg                                                                                                   <  2:315952/110‑1 (MQ=255)
                                             ccGATGTCGATTTTGTCATCCAGATACCAGGCGTGCGCCTCGGTGCCGGAACGCGTATCGCGGTCGTAAGGGCTTGAGTCGGACGGCAACTGCCCGCTGCTGGTGGCGGTGTAGTAACGCATTTCATGCGTTGattcat                                                                     <  1:1645378/139‑1 (MQ=255)
                                               gATGTCGATTTTGTCATCCAGATACCAGGCGTGCGCCTCGGTGCCGGAACGCGTATCGCGGTCGTAAGGGCTTGAGTCGGACGGCAACTGCCCGCTGCTGGTGGCGGTGTAGTAACGCATTTCATGCGTTGattcaatc                                                                   >  2:631210/1‑136 (MQ=255)
                                                                          ggCGTGCGCCTCGGTGCCGGAACGCGTATCGCGGTCGTAAGGGCTTGAGTCgg                                                                                                                              >  2:1366319/1‑53 (MQ=255)
                                                                          ggCGTGCGCCTCGGTGCCGGAACGCGTATCGCGGTCGTAAGGGCTTGAGTCgg                                                                                                                              <  1:1366319/53‑1 (MQ=255)
                                                                          ggCGTGCGCCTCGGTGCCGGAACGCGTATCGCGGTCGTAAGGGCTTGAGTCGGACGGCAACTGCCCGCTGCTGGTGGCGGTGTAGTAACGCATTTCATGCg                                                                              >  2:917331/1‑101 (MQ=255)
                                                                          ggCGTGCGCCTCGGTGCCGGAACGCGTATCGCGGTCGTAAGGGCTTGAGTCGGACGGCAACTGCCCGCTGCTGGTGGCGGTGTAGTAACGCATTTCATGCg                                                                              <  1:917331/101‑1 (MQ=255)
                                                                          ggCGTGCGCCTCGGTGCCGGAACGCGTATCGCGGTCGTAAGGGCTTGAGTCGGACGGCAACTGCCCGCTGCTGGTGGCGGTGTAGTAACGCATTTCATGCGTTGATTCATTCAAATAGCGATAGCCCACGCCCACTTCg                                        >  2:1174608/1‑139 (MQ=255)
                                                                           gCGTGCGCCTCGGTGCCGGAACGCGTATCGCGGTCGTAAGGGCTTGAGTCGGACGGCAACTGCCCGCTGCTGGTGGCGGTGTAGTAACGCATTTCATGCg                                                                              <  2:427610/100‑1 (MQ=255)
                                                                           gCGTGCGCCTCGGTGCCGGAACGCGTATCGCGGTCGTAAGGGCTTGAGTCGGACGGCAACTGCCCGCTGCTGGTGGCGGTGTAGTAACGCATTTCATGCg                                                                              >  1:427610/1‑100 (MQ=255)
                                                                              tGCGCCTCGGTGCCGGAACGCGTATCGCGGTCGTAAGGGCTTGAGTCGGACGGCAACTGcc                                                                                                                  <  1:1523390/61‑1 (MQ=255)
                                                                              tGCGCCTCGGTGCCGGAACGCGTATCGCGGTCGTAAGGGCTTGAGTCGGACGGCAACTGcc                                                                                                                  >  2:1523390/1‑61 (MQ=255)
                                                                              tGCGCCTCGGTGCCGGAACGCGTATCGCGGTCGTAAGGGCTTGAGTCGGACGGCAACTGCCCGCTGCtggtgg                                                                                                      >  1:1405392/1‑73 (MQ=255)
                                                                              tGCGCCTCGGTGCCGGAACGCGTATCGCGGTCGTAAGGGCTTGAGTCGGACGGCAACTGCCCGCTGCtggtgg                                                                                                      <  2:1405392/73‑1 (MQ=255)
                                                                              tGCGCCTCGGTGCCGGAACGCGTATCGCGGTCGTAAGGGCTTGAGTCGGACGGCAACTGCCCGCTGCTGGTGGCGGTGTAGTAACGCATTTCATGCGTTGATTCATTCAAATAGCGATAGCCCACGCCCACTTCg                                        <  2:1189603/135‑1 (MQ=255)
                                                                              tGCGCCTCGGTGCCGGAACGCGTATCGCGGTCGTAAGGGCTTGAGTCGGACGGCAACTGCCCGCTGCTGGTGGCGGTGTAGTAACGCATTTCATGCGTTGATTCATTCAAATAGCGATAGCCCACGCCCACTTCg                                        >  1:1189603/1‑135 (MQ=255)
                                                                              tGCGCCTCGGTGCCGGAACGCGTATCGCGGTCGTAAGGGCTTGAGTCGGACGGCAACTGCCCGCTGCTGGTGGCGGTGTAGTAACGCATTTCATGCGTTGATTCATTCAAATAGCGATAGCCCACGCCCACTTCGTgcg                                    >  2:1569374/1‑139 (MQ=255)
                                                                              tGCGCCTCGGTGCCGGAACGCGTATCGCGGTCGTAAGGGCTTGAGTCGGACGGCAACTGCCCGCTGCTGGTGGCGGTGTAGTAACGCATTTCATGCGTTGATTCATTCAAATAGCGATAGCCCACGCCCACTTCGTgcg                                    >  1:94082/1‑139 (MQ=255)
                                                                                cgcCTCGGTGCCGGAACGCGTATCGCGGTCGTAAGGGCTTGAGTCGGACGGCAACTGCCCGCTGCtggtgg                                                                                                      >  1:879658/1‑71 (MQ=255)
                                                                                cgcCTCGGTGCCGGAACGCGTATCGCGGTCGTAAGGGCTTGAGTCGGACGGCAACTGCCCGCTGCtggtgg                                                                                                      <  2:879658/71‑1 (MQ=255)
                                                                                cgcCTCGGTGCCGGAACGCGTATCGCGGTCGTAAGGGCTTGAGTCGGACGGCAACTGCCCGCTGCTGGTGGCGGTCTAGTAACGCATTTCATGCGTTGATTCATTCAAATAGCGATAGCCCACGCCCACTTCGTGcgcg                                  <  2:1101363/139‑1 (MQ=255)
                                                                                  ccTCGGTGCCGGAACGCGTATCGCGGTCGTAAGGGCTTGAGTCGGACGGCAACTGCCCGCTGCTGGTGGCGGTGTAGTAACGCATTTCATGCGTTGATTCATTCAAATAGCGATAGCCCACGCCCACTTCGTGCGCGGa                                <  2:94082/139‑1 (MQ=255)
                                                                                         gCCGGAACGCGTATCGCGGTCGTAAGGGCTTGAGTCGGACGGCAACTGCCCGCTGCTGGTGGCGGTGTAGTAACGCATTTCATGCGTTGATTCATTCAAATAGCGATAGCCCACGCCCACTTCGTGCGCGGAAGGGCCg                         <  1:1174608/139‑1 (MQ=255)
                                                                                            ggAACGCGTATCGCGGTCGTAAGGGCTTGAGTCGGACGGCAACTGCCCGCTGCTGGTGGCGGTgta                                                                                               <  1:1583207/66‑1 (MQ=255)
                                                                                            ggAACGCGTATCGCGGTCGTAAGGGCTTGAGTCGGACGGCAACTGCCCGCTGCTGGTGGCGGTgta                                                                                               >  2:1583207/1‑66 (MQ=255)
                                                                                               aCGCGTATCGCGGTCGTAAGGGCTTGAGTCGGACGGCAACTGCCCGCTGCTGGTGGCGGTGTAGTAACGCATTTCATGCGTTGATTCATTCAAATAGCGATAGCCCACGCCCACTTCGTGCGCGGAAGGGCCGATCATa                   <  1:1627160/139‑1 (MQ=255)
                                                                                                  cgTATCGCGGTCGTAAGGGCTTGAGTCGGACGGCAACTGCCCGCTGCtggtgg                                                                                                      <  2:1223130/53‑1 (MQ=255)
                                                                                                  cgTATCGCGGTCGTAAGGGCTTGAGTCGGACGGCAACTGCCCGCTGCtggtgg                                                                                                      >  1:1223130/1‑53 (MQ=255)
                                                                                                  cgTATCGCGGTCGTAAGGGCTTGAGTCGGACGGCAACTGCCCGCTGCTGGTGGCg                                                                                                    <  1:1084197/55‑1 (MQ=255)
                                                                                                  cgTATCGCGGTCGTAAGGGCTTGAGTCGGACGGCAACTGCCCGCTGCTGGTGGCg                                                                                                    >  2:1084197/1‑55 (MQ=255)
                                                                                                  cgTATCGCGGTCGTAAGGGCTTGAGTCGGACGGCAACTGCCCGCTGCTGGTGGCGGTGTAGTAACGCATTTCATGCGTTGATTCATTCaaa                                                                <  1:724948/91‑1 (MQ=255)
                                                                                                  cgTATCGCGGTCGTAAGGGCTTGAGTCGGACGGCAACTGCCCGCTGCTGGTGGCGGTGTAGTAACGCATTTCATGCGTTGATTCATTCaaa                                                                >  2:724948/1‑91 (MQ=255)
                                                                                                  cgTATCGCGGTCGTAAGGGCTTGAGTCGGACGGCAACTGCCCGCTGCTGGTGGCGGTGTAGTAACGCATTTCATGCGTTGATTCATTCAAATAGCGATAGCCCACGCCCACTTCGTGCGCGGAAGGGc                           <  2:557786/128‑1 (MQ=255)
                                                                                                  cgTATCGCGGTCGTAAGGGCTTGAGTCGGACGGCAACTGCCCGCTGCTGGTGGCGGTGTAGTAACGCATTTCATGCGTTGATTCATTCAAATAGCGATAGCCCACGCCCACTTCGTGCGCGGAAGGGc                           >  1:557786/1‑128 (MQ=255)
                                                                                                  cgTATCGCGGTCGTAAGGGCTTGAGTCGGACGGCAACTGCCCGCTGCTGGTGGCGGTGTAGTAACGCATTTCATGCGTTGATTCATTCAAATAGCGATAGCCCACGCCCACTTCGTGCGCGGAAGGGCCGATCATAAAg                >  1:334438/1‑139 (MQ=255)
                                                                                                    tATCGCGGTCGTAAGGGCTTGAGTCGGACGGCAACTGCCCGCTGCTGGTGGCGGTGTAGTAACGCATTTCATGCGTTGATTCATTCAAATAGCGATAGCCCACGCCCACTTCGTGCGCGGAAGGGCCGATCATAAAGAt              >  1:1321596/1‑139 (MQ=255)
                                                                                                     aTCGCGGTCGTAAGGGCTTGAGTCGGACGGCAACTGCCCGCTGCtggtgg                                                                                                      >  1:993261/1‑50 (MQ=255)
                                                                                                     aTCGCGGTCGTAAGGGCTTGAGTCGGACGCCAACTGCCCGCTGCtggtgg                                                                                                      <  2:993261/50‑1 (MQ=255)
                                                                                                               tAAGGGCTTGAGTGGACGGCAACTGCCCGCTGCTGGTGGCGGTGTAGTAACGCATTTCATGCGTTGATTCATTCAAATAGCGATAGCCCACGCCCACTTCGTGCGCGGAAGGGCCGATCATAAAGATCTGGCTGTAGCg  <  2:166518/139‑1 (MQ=255)
                                                                                                                aaGGGCTTGAGTCGGACGGCAACTGCCCGCTGCTGGTGGCGGTGTAGTAACGCATTTCATGCGTTGATTCATTCAAATAGCGATAGCCCACGCCCACTTCGTGCGCGGAAGGGCCGATCATAAAGATCTGGCTGTAGCg  <  2:402461/139‑1 (MQ=255)
                                                                                                                aaGGGCTTGAGTCGGACGGCAACTGCCCGCTGCTGGTGGCGGTGTAGTAACGCATTTCATGCGTTGATTCAATCAAATAGCGATAGCCCACGCCCACTTCGTGCGCGGAAGGGCCGATCATAAAGATCTGGCTGTAGCg  <  1:631210/139‑1 (MQ=255)
                                                                                                                      ttGAGTCGGACGGCAACTGCCCGCTGCTGGTGGCGGTGTAGTAACGCATTTCATGCGTTGATTCATTCaaa                                                                >  2:39277/1‑71 (MQ=255)
                                                                                                                      ttGAGTCGGACGGCAACTGCCCGCTGCTGGTGGCGGTGTAGTAACGCATTTCATGCGTTGATTCATTCaaa                                                                <  1:39277/71‑1 (MQ=255)
                                                                                                                           |                                                                                                                               
GACTCGATATGTTCGAAACGCATACCCGGCGTGATGGTCCAGTTGCCGATGTCGATTTTGTCATCCAGATACCAGGCGTGCGCCTCGGTGCCGGAACGCGTATCGCGGTCGTAAGGGCTTGAGCCGGACGGCAACTGCCCGCTGCTGGTGGCGGTGTAGTAACGCATTTCATGCGTTGATTCATTCAAATAGCGATAGCCCACGCCCACTTCGTGCGCGGAAGGGCCGATCATAAAGATCTGGCTGTAGCG  >  NC_000913/4515161‑4515411

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