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,2840CT98.0% 155.2 / ‑4.1 50G465D (GGC→GAC) 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 (1/0);  new base T (16/33);  total (17/33)
Fisher's exact test for biased strand distribution p-value = 3.40e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 8.32e-01

TTCTGGTATGACTCGATATGTTCGAAACGCATACCCGGCGTGATGGTCCAGTTGCCGATGTCGATTTTGTCATCCAGATACCAGGCGTGCGCCTCGGTGCCGGAACGCGTATCGCGGTCGTAAGGGCTTGAGCCGGACGGCAACTGCCCGCTGCTGGTGGCGGTGTAGTAACGCATTTCATGCGTTGATTCATTCAAATAGCGATAGCCCACGCCCACTTCGTGCGCGGAAGGGCCGATCATAAAGATCTGGCTGTAGCG  >  NC_000913/4515152‑4515411
                                                                                                                                    |                                                                                                                               
ttCTGGTATGACTCGATATGTTCGAAACGCATACCCGGCGTGATGGTCCAGTTGCCGATGTCGATTTTGTCATCCAGATACCAGGCGTGCGCCTCGGTGCCGGAACGCGTATCGCGGTCGTAAGGGCTTGAGTCGGACg                                                                                                                           <  1:102207/139‑1 (MQ=255)
          aCTCGATATGTTCGAAACGCATACTCGGCGTGATGGTCCAGCTGCCGATGTCGATTTTGTCATCCAGATACCAGGCGTGCGCCTCGGTGCCGGAACGCGTATCGCGGTCGTAAGGGCTTGAGTCGGACGGCAACTGccc                                                                                                                 >  1:48638/1‑139 (MQ=255)
          aCTCGATATGTTCGAAACGCATACCCGGCGTGATGGTCCAGTTGCCGATGTCGATTTTGTCATCCAGATACCAGGCGTGCGCCTCGGTGCCGGAACGCGTATCGCGGTCGTAAGGGCTTGAGTCGGACGGCAACTGccc                                                                                                                 >  1:159002/1‑139 (MQ=255)
          aCTCGATATGTTCGAAACGCATACCCGGCGTGATGGTCCAGTTGCCGATGTCGATTTTGTCATCCAGATACCAGGCGTGCGCCTCGGTGCCGGAACGCGTATCGCGGTCGTAAGGGCTTGAGTCGGACGGCAACTGccc                                                                                                                 >  2:280419/1‑139 (MQ=255)
          aCTCGATATGTTCGAAACGCATACCCGGCGTGATGGTCCAGTTGCCGATGTCGATTTTGTCATCCAGATACCAGGCGTGCGCCTCGGTGCCGGAACGCGTATCGCGGTCGTAAGGGCTTGAGTCGGACGGCAACTGccc                                                                                                                 >  2:325600/1‑139 (MQ=255)
          aCTCGATATGTTCGAAACGCATACCCGGCGTGATGGTCCAGTTGCCGATGTCGATTTTGTCATCCAGATACCAGGCGTGCGCCTCGGTGCCGGAACGCGTATCGCGGTCGTAAGGGCTTGAGTCGGACGGCAACTGccc                                                                                                                 >  1:766787/1‑139 (MQ=255)
          aCTCGATATGTTCGAAACGCATACCCGGCGTGATGGTCCAGTTGCCGATGTCGATTTTGTCATCCAGATACCAGGCGTGCGCCTCGGTGCCGGAACGCGTATCGCGGTCGTAAGGGCTTGAGTCGGACGGCAACTGccc                                                                                                                 >  1:475337/1‑139 (MQ=255)
             cGATATGTTCGAAACGCATACCCGGCGTGATGGTCCAGTTGCCGATGTCGATTTTGTCATCCAGATACCAGGCGTGCGCCTCGGTGCCGGAACGCGTATCGCGGTCGTAAGGGCTTGAGTCGGACGGCAACTGCCCgct                                                                                                              <  2:637986/139‑1 (MQ=255)
               atatGTTCGAAACGCATACCCGGCGTGATGGTCCAGTTGCCGATGTCGATTTTGTCATCCAGATACCAGGCGTGCGCCTCGGTGCCGGAACGCGTATCGCGGTCGTAAGGGCTTGAGTCGGACGGCAACTGCCCgctgc                                                                                                            >  1:49075/1‑139 (MQ=255)
                  tGTTCGAAACGCATACCCGGCGTGATGGTCCAGTTGCCGATGTCGATTTTGTCATCCAGATACCAGGCGTGCGCCTCGGTGCCGGAACGCGTATCGCGGTCGTAAGGGCTTGAGTCGGACGGCAACTGCCCGCTGCtgg                                                                                                         <  1:557769/139‑1 (MQ=255)
                    ttCGAAACGCATACCCGGCGTGATGGTCCAGTTGCCGATGTCGATTTTGTCATCCAGATACCAGGCGTGCGCCTCGGTGCCGGAACGCGTATCGCGGTCGTAAGGGCTTGAGTCGGACGGCAACTGCCCGCTGCtgg                                                                                                         <  1:48189/137‑1 (MQ=255)
                    ttCGAAACGCATACCCGGCGTGATGGTCCAGTTGCCGATGTCGATTTTGTCATCCAGATACCAGGCGTGCGCCTCGGTGCCGGAACGCGTATCGCGGTCGTAAGGGCTTGAGTCGGACGGCAACTGCCCGCTGCtgg                                                                                                         >  2:48189/1‑137 (MQ=255)
                     tCGAAACGCATACCCGGCGTGATGGTCCAGTTGCCGATGTCGATTTTGTCATCCAGATACCAGGCGTGCGCCTCGGTGCCGGAACGCGTATCGCGGTCGTAAGGGCTTGAGTCGGACGGCAACTGCCCGCTGCtggtgg                                                                                                      <  2:625098/139‑1 (MQ=255)
                     tCGAAACGCATACCCGGCGTGATGGTCCAGTTGCCGATGTCGATTTTGTCATCCAGATACCAGGCGTGCGCCTCGGTGCCGGAACGCGTATCGCGGTCGTAAGGGCTTGAGTCGGACGGCAACTGCCCGCTGCtggtgg                                                                                                      <  2:222652/139‑1 (MQ=255)
                     tCGAAACGCATACCCGGCGTGATGGTCCAGTTGCCGATGTCGATTTTGTCATCCAGATACCAGGCGTGCGCCTCGGTGCCGGAACGCGTATCGCGGTCGTAAGGGCTTGAGTCGGACGGCAACTGCCCGCTGCtggtgg                                                                                                      <  2:284339/139‑1 (MQ=255)
                     tCGAAACGCATACCCGGCGTGATGGTCCAGTTGCCGATGTCGATTTTGTCATCCAGATACCAGGCGTGCGCCTCGGTGCCGGAACGCGTATCGCGGTCGTAAGGGCTTGAGTCGGACGGCAACTGCCCGCTGCtggtgg                                                                                                      >  1:743894/1‑139 (MQ=255)
                     tCGAAACGCATACCCGGCGTGATGGTCCAGTTGCCGATGTCGATTTTGTCATCCAGATACCAGGCGTGCGCCTCGGTGCCGGAACGCGTATCGCGGTCGTAAGGGCTTGAGTCGGACGGCAACTGCCCGCTGCtggtgg                                                                                                      <  2:49075/139‑1 (MQ=255)
                          aCGCATACCCGGCGTGATGGTCCAGTTGCCGATGTCGATTTTGTCATCCAGATACCAGGCGTGCGCCTCGGTGCCGGAACGCGTATCGCGGTCGTAAGGGCTTGAGTCGGACGGCAACTGCCCGCTGCTGGTGGCGgtg                                                                                                 <  1:428949/139‑1 (MQ=255)
                                aCCCGGCGTGATGGTCCAGTTGCCGATGTCGATTTTGTCATCCAGATACCAGGCGTGCGCCTCGGTGCCGGAACGCGTATCGCGGTCGTAAGGGCTTGAGTCGGACGGCAACTGCCCGCTGCTGGTGGCGGTGTAGTaa                                                                                           <  2:273450/139‑1 (MQ=255)
                                aCCCGGCGTGATGGTCCAGTTGCCGATGTCGATTTTGTCATCCAGATACCAGGCGTGCGCCTCGGTGCCGGAACGCGTATCGCGGTCGTAAGGGCTTGAGTCGGACGGCAACTGCCCGCTGCTGGTGGCGGTGTAGTaa                                                                                           <  2:766787/139‑1 (MQ=255)
                                      cGTGATGGTCCAGTTGCCGATGTCGATTTTGTCATCCAGATACCAGGCGTGCGCCTCGGTGCCGGAACGCGTATCGCGGTCGTAAGGGCTTGAGCCGGACGGCAACTGCCCGCTGCTGGTGGCGGTGTAGTAACGCAtt                                                                                     >  2:474714/1‑139 (MQ=255)
                                        tGATGGTCCAGTTGCCGATGTCGATTTTGTCATCCAGATACCAGGCGTGCGCCTCGGTGCCGGAACGCGTATCGCGGTCGTAAGGGCTTGAGTCGGACGGCAACTGCCCGCTGCTGGTGGCGGTGTAGTAACGCATTTc                                                                                   <  1:658997/139‑1 (MQ=255)
                                        tGATGGTCCAGTTGCCGATGTCGATTTTGTCATCCAGATACCAGGCGTGCGCCTCGGTGCCGGAACGCGTATCGCGGTCGTAAGGGCTTGAGTCGGACGGCAACTGCCCGCTGCTGGTGGCGGTGTAGTAACGCATTTc                                                                                   <  2:743894/139‑1 (MQ=255)
                                         gATGGTCCAGTTGCCGATGTCGATTTTGTCATCCAGATACCAGGCGTGCGCCTCGGTGCCGGAACGCGTATCGCGGTCGTAAGGGCTTGAGTCGGACGGCAACTGCCCGCTGCTGGTGGCGGTGTAGTAACGCATTTCa                                                                                  <  2:692708/139‑1 (MQ=255)
                                          aTGGTCCAGTTGCCGATGTCGATTTTGTCATCCAGATACCAGGCGTGCGCCTCGGTGCCGGAACGCGTATCGCGGTCGTAAGGGCTTGAGTCGGACGGCAACTGCCCGCTGCTGGTGGCGGTGTAGTAACGCATTTCAt                                                                                 >  2:49491/1‑139 (MQ=255)
                                             gTCCAGTTGCCGATGTCGATTTTGTCATCCAGATACCAGGCGTGCGCCTCGGTGCCGGAACGCGTATCGCGGTCGTAAGGGCTTGAGTCGGACGGCAACTGCCCGCTGCTGGTGGCGGTGTAGTAACGCATTTCATGCg                                                                              <  2:159751/139‑1 (MQ=255)
                                                   ttGCCGATGTCGATTTTGTCATCCAGATACCAGGCGTGCGCCTCGGTGCCGGAACGCGTATCGCGGTCGTAAGGGCTTGAGTCGGACGGCAACTGCCCGCTGCTGGTGGCGGTGTAGTAACGCATTTCATGCGTTGatt                                                                        <  1:432583/139‑1 (MQ=255)
                                                       cGATGTCGATTTTGTCATCCAGATACCAGGCGTGCGCCTCGGTGCCGGAACGCGTATCGCGGTCGTAAGGGCTTGAGTCGGACGGCAACTGCCCGCTGCTGGTGGCGGTGTAGTAACGCATTTCATGCGTTGattcatt                                                                    <  1:49491/139‑1 (MQ=255)
                                                            tCGATTTTGTCATCCAGATACCAGGCGTGCGCCTCGGTGCCGGAACGCGTATCGCGGTCGTAAGGGCTTGAGTCGGACGGCAACTGCCCGCTGCTGGTGGCGGTGTAGTAACGCATTTCATGCGTTGATTCATTCAAAt                                                               <  1:338087/139‑1 (MQ=255)
                                                                 tttGTCATCCAGATACCAGGCGTGCGCCTCGGTGCCGGAACGCGTATCGCGGTCGTAAGGGCTTGAGTCGGACGGCAACTGCCCGCTGCTGGTGGCGGTGTAGTAACGCATTTCATGCGTTGATTCATTCAAATAGCGa                                                          <  1:280419/139‑1 (MQ=255)
                                                                    gTCATCCAGATACCAGGCGTGCGCCTCGGTGCCGGAACGCGTATCGCGGTCGTAAGGGCTTGAGTCGGACGGCAACTGCCCGCTGCTGGTGGCGGTGTAGTAACGCATTGCATGCGTTGATTCATTCAAATAGCGATAg                                                       <  2:795385/139‑1 (MQ=255)
                                                                      cATCCAGATACCAGGCGTGCGCCTCGGTGCCGGAACGCGTATCGCGGTCGTAAGGGCTTGAGTCGGACGGCAACTGCCCGCTGCTGGTGGCGGTGTAGTAACGCATTTCATGCGTTGATTCATTCAAATAGCGATAgcc                                                     <  1:800542/139‑1 (MQ=255)
                                                                      cATCCAGATACCAGGCGTGCGCCTCGGTGCCGGAACGCGTATCGCGGTCGTAAGGGCTTGAGTCGGACGGCAACTGCCCGCTGCTGGTGGCGGTGTAGTAACGCATTTCATGCGTTGATTCATTCAAATAGCGATAgcc                                                     <  1:325600/139‑1 (MQ=255)
                                                                      cATCCAGATACCAGGCGTGCGCCTCGGTGCCGGAACGCGTATCGCGGTCGTAAGGGCTTGAGTCGGACGGCAACTGCCCGCTGCTGGTGGCGGTGTAGTAACGCATTTCATGCGTTGATTCATTCAAATAGCGATAgcc                                                     <  1:389906/139‑1 (MQ=255)
                                                                          cAGATACCAGGCGTGCGCCTCGGTGCCGGAACGCGTATCGCGGTCGTAAGGGCTTGAGTCGGACGGCAACTGCCCGCTGCTGGTGGCGGTGTAGTAACGCATTTCATGCGTTGATTCATTCAAATAGCGATAgcccacg                                                 <  2:333175/139‑1 (MQ=255)
                                                                          cAGATACCAGGCGTGCGCCTCGGTGCCGGAACGCGTATCGCGGTCGTAAGGGCTTGAGTCGGACGGCAACTGCCCGCTGCTGGTGGCGGTGTAGTAACGCATTTCATGCGTTGATTCATTCAAATAGCGATAgcc                                                     >  1:636726/1‑135 (MQ=255)
                                                                          cAGATACCAGGCGTGCGCCTCGGTGCCGGAACGCGTATCGCGGTCGTAAGGGCTTGAGTCGGACGGCAACTGCCCGCTGCTGGTGGCGGTGTAGTAACGCATTTCATGCGTTGATTCATTCAAATAGCGATAgcc                                                     <  2:636726/135‑1 (MQ=255)
                                                                                     cGTGCGCCTCGGTGCCGGAACGCGTATCGCGGTCGTAAGGGCTTGAGTCGGACGGCAACTGCCCGCTGCTGGTGGCGGTGTAGTAACGCATTTCATGCGTTGATTCATTCAAATAGCGATAGCCCACGCCCACTTCGTg                                      <  1:246328/139‑1 (MQ=255)
                                                                                       tGCGCCTCGGTGCCGGAACGCGTATCGCGGTCGTAAGGGCTTGAGTCGGACGGCAACTGCCCGCTGCTGGTGGCGGTGTAGTAACGCATTTCATGCGTTGattca                                                                      <  2:327277/105‑1 (MQ=255)
                                                                                       tGCGCCTCGGTGCCGGAACGCGTATCGCGGTCGTAAGGGCTTGAGTCGGACGGCAACTGCCCGCTGCTGGTGGCGGTGTAGTAACGCATTTCATGCGTTGattca                                                                      >  1:327277/1‑105 (MQ=255)
                                                                                         cgcCTCGGTGCCGGAACGCGTATCGCGGTCGTAAGGGCTTGAGTCGGACGGCAACTGCCCGCTGCtggtgg                                                                                                      <  2:175040/71‑1 (MQ=255)
                                                                                         cgcCTCGGTGCCGGAACGCGTATCGCGGTCGTAAGGGCTTGAGTCGGACGGCAACTGCCCGCTGCtggtgg                                                                                                      >  1:175040/1‑71 (MQ=255)
                                                                                         cgcCTCGGTGCCGGAACGCGTATCGCGGTCGTAAGGGCTTGAGTCGGACGGCAACTGCCCGCTGCTGGTGGCGGTGTAGTAACGCATTTCATGCGTTGATTCATTCAAATAGCGATAGCCCACGCCCACTTCGTGcgcg                                  >  1:554055/1‑139 (MQ=255)
                                                                                         cgcCTCGGTGCCGGAACGCGTATCGCGGTCGTAAGGGCTTGAGTCGGACGGCAACTGCCCGCTGCTGGTGGCGGTGTAGTAACGCATTTCATGCGTTGATTCATTCAAATAGCGATAGCCCACGCCCACTTCGTGcgcg                                  >  1:641712/1‑139 (MQ=255)
                                                                                         cgcCTCGGTGCCGGAACGCGTATCGCGGTCGTAAGGGCTTGAGTCGGACGGCAACTGCCCGCTGCTGGTGGCGGTGTAGTAACGCATTTCATGCGTTGATTCATTCAAATAGCGATAGCCCACGCCCACTTCGTGcgcg                                  >  1:737798/1‑139 (MQ=255)
                                                                                                   ccGGAACGCGTATCGCGGTCGTAAGGGCTTGAGTCGGACGGCAACTGCCCGCTGCTGGTGGCGGTGTAGTAACGCATTTCATGCGTTGATTCATTCAAATAGCGATAGCCCACGCCCACTTCGTGCGCGGAAGGGCCGa                        <  2:232571/139‑1 (MQ=255)
                                                                                                               tCGCGGTCGTAAGGGCTTGAGTCGGACGGCAACTGCCCGCTGCTGGTGGCGGTGTAGTAACGCATTTCATGCGTTGATTCATTCAAATAGCGATAGCCCACGCCCACTTCGTGCGCGGAAGGGCCGATCATAAAGATCt            <  1:439444/139‑1 (MQ=255)
                                                                                                                       gTAAGGGCTTGAGTCGGACGGCAACTGCCCGCTGCTGGTGGCGGTGTAGTAACGCATTTCATGCGTTGATTCATTCAAATAGCGATAGCCCACGCCCACTTCGTGCGCGGAAGGGCCGATCATAAAGATCTGGCTGTAg    <  2:475337/139‑1 (MQ=255)
                                                                                                                         aaGGGCTTGAGTCGGACGGCAACTGCCCGCTGCTGGTGGCGGTGTAGTAACGCATTTCATGCGTTGATTCATTCAAATAGCGATAGCCCACGCCCACTTCGTGCGCGGAAGGGCCGATCATAAAGATCTGGCTGTAGCg  <  2:641712/139‑1 (MQ=255)
                                                                                                                         aaGGGCTTGAGTCGGACGGCAACTGCCCGCTGCTGGTGGCGGTGTAGTAACGCATTTCATGCGTTGATTCATTCAAATAGCGATAGCCCACGCCCACTTCGTGCGCGGAAGGGCCGATCATAAAGATCTGGCTGTAGCg  <  2:737798/139‑1 (MQ=255)
                                                                                                                                    |                                                                                                                               
TTCTGGTATGACTCGATATGTTCGAAACGCATACCCGGCGTGATGGTCCAGTTGCCGATGTCGATTTTGTCATCCAGATACCAGGCGTGCGCCTCGGTGCCGGAACGCGTATCGCGGTCGTAAGGGCTTGAGCCGGACGGCAACTGCCCGCTGCTGGTGGCGGTGTAGTAACGCATTTCATGCGTTGATTCATTCAAATAGCGATAGCCCACGCCCACTTCGTGCGCGGAAGGGCCGATCATAAAGATCTGGCTGTAGCG  >  NC_000913/4515152‑4515411

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