Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 4,515,480 C→T 100% G400S (GGC→AGC)  fecA ← ferric citrate outer membrane transporter

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009134,515,4800CT100.0% 93.1 / NA 26G400S (GGC→AGC) fecAferric citrate outer membrane transporter
Reads supporting (aligned to +/- strand):  ref base C (0/0);  new base T (14/12);  total (14/12)
Rejected as polymorphism: Frequency below/above cutoff threshold.
Rejected as polymorphism: Variant not supported by required number of reads on each strand.

CGCGGTCGTAAGGGCTTGAGCCGGACGGCAACTGCCCGCTGCTGGTGGCGGTGTAGTAACGCATTTCATGCGTTGATTCATTCAAATAGCGATAGCCCACGCCCACTTCGTGCGCGGAAGGGCCGATCATAAAGATCTGGCTGTAGCGTGGCTCAATACCGCGCACCCAGTAGTTACGCGGCGAGAGGGTGATGCGTTTGCCTTGCTCCAGGTAGCCGCTGCGCAGGGTTTGGGTGTAGAACCCCTGAATGTTGAATTTATGCTGGCTGTCTGGCTGGAACTGGTAGCCCAGGCTCGCCAGCTTGCGACGACCCCAGAAGCGGTCATACGGGCGGGTGGATTGCCAGCGATCGGCGTCGTAATCCGCGCGAGACAGGCCACCGGGCATGTCGGCTTCACCGTCGTAATATTGCAG  >  NC_000913/4515264‑4515678
                                                                                                                                                                                                                        |                                                                                                                                                                                                      
cgcgGTCGTAAGTGCTTGAGCCGGACGGCAACTGCCCGCGGCTGGTGGCGGTGTAGTAACGCATTTCATGCGTTGATTCATTCAAATAGCGATAGCCCACGCCCACTTCGTGCGCGGAAGGGCCGATCATAAAGATCTGGCTGTAGCGTGGCTCAATACCGCGCACCCAGTAGTTACGCGGCGAGAGGGTGATGCGTTTGCCTTGCTCCAGGTAGCTGCTGCGCAGGGTTTGGGTGTAGAACCCCTGAATGTTGAATTTATGCTGGCTGTc                                                                                                                                                  <  1:141495/271‑1 (MQ=255)
   ggTCGTAAGGGCTTGAGCCGGACGGCAACTGCCCGCTGCTGGTGGCGGTGTAGTAACGCATTTCATGCGTTGATTCATTCAAATAGCGATAGCCCACGCCCACTTCGTGCGCGGAAGGGCCGATCATAAAGATCTGGCTGTAGCGTGGCTCAATACCGCGCACCCAGTAGTTACGCGGCGAGAGGGTGATGCGTTTGCCTTGCTCCAGGTAGCTGCTGCGCAGGGTTTGGGTGTAGAACCCCTGAATGTTGAATTTATGCTGGCTGTctgg                                                                                                                                               >  1:453748/1‑271 (MQ=255)
           gggCTTGAGCCGGACGGCAACTGCCCGCTGCTGGTGGCGGTGTAGTAACGCATTTCATGCGTTGATTCATTTAAATAGCGATAGCCCACGCCCACTTCGTGCGCGGAAGGGCCGATCCTAAAGATCTGGCTGTAGCGTGGCTCAATACCGCGCACCCAGTAGTTACGCGGCGAGAGGGTGATGCGTTTGCCTTGCTCCAGGTAGCTGCTGCGCAGGGGTTGGGTGTAGAACCCCTGAATGTTGAATTTATGCTGGCTGTCTGGCGGGAACt                                                                                                                                       >  1:430234/1‑271 (MQ=255)
                          ggCAACTGCCCGCTGCTGGTGGCGGTGTAGTAACGCATTTCATGCGTTGATTCATTCAAATAGCGATAGCCCACGCCCACTTCGTGCGCGGAAGGGCCGATCATAAAGATCTGGCTGTAGCGTGGCTCAATACCGCGCACCCAGTAGTTACGCGGCGAGAGGGTGATGCGTTTGCCTTGCTCCAGGTAGCTGCTGCGCAGGGTTTGGGTGTAGAACCCCTGAATGTTGAATTTATGCTGGCTGTCTGGCTGGAACTGGTAGCCCAGGCTCg                                                                                                                        >  1:88673/1‑271 (MQ=255)
                                     gctgctGGTGGCGGTGTAGTAACGCATTTCATGCGTTGATTCATTCAAATAGCGATAGCCCACGCCCACTTCGTGCGCGGAAGGGCCGATCATAAAGATCTGGCTGTAGCGTGGCTCAATACCGCGCACCCAGTAGTTACGCGGCGAGAGGGTGATGCGTTTGCCTTGCTCCAGGTAGCTGCTGCGCa                                                                                                                                                                                                >  2:383848/1‑188 (MQ=255)
                                      ctgctgGTGGCGGTGTAGTAACGCATTTCATGCGTTGATTCATTCAAATAGCGATAGCCCACGCCCACTTCGTGCGCGGAAGGGCCGATCATAAAGATCTGGCTGTAGCGTGGCTCAATACCGCGCACCCAGTAGTTACGCGGCGAGAGGGTGATGCGTTTGCCTTGCTCCAGGTAGCTGCTGCGCAg                                                                                                                                                                                               <  1:383848/188‑1 (MQ=255)
                                                      agtaACGCATTTCATGCGTTGATTCATTCAAATAGCGATAGCCCACGCCCACTTCGTGCGCGGAAGGGCCGATCATAAAGATCTGGCTGTAGCGTGGCTCAATACCGCGCACCCAGTAGTTACGCGGCGAGAGGGTGATGCGTTTGCCTTGCTCCAGGTAGCTGCTGCGCAgg                                                                                                                                                                                              >  1:304963/1‑173 (MQ=255)
                                                       gtaACGCATTTCATGCGTTGATTCATTCAAATAGCGATAGCCCACGCCCACTTCGTGCGCGGAAGGGCCGATCATAAAGATCTGGCTGTAGCGTGGCTCAATACCGCGCACCCAGTAGTTACGCGGCGAGAGGGTGATGCGTTTGCCTTGCTCCAGGTAGCTGCTGCGCAggg                                                                                                                                                                                             <  2:304963/173‑1 (MQ=255)
                                                           cGCATTTCATGCGTTGATTCATTCAAATAGCGATAGCCCACGCCCACTTCGTGCGCGGAAGGGCCGATCATAAAGATCTGGCTGTAGCGTGGCTCAATACCGCGCACCCAGTAGTTACGCGGCGAGAGGGTGATGCGTTTGCCTTGCTCCAGGTAGCTGCTGCGCAGGGTTTGGGTGTAGAACCCCTGAATGTTGAATTTATGCTGGCTGTCTGGCTGGAACTGGTAGCCCAGGCTCGCCAGCTTGCGACGAACCCAGacggcgtcagacag                                                                                      >  2:87554/1‑259 (MQ=255)
                                                                                  caAATAGCGATAGCCCACGCCCACTTCGTGCGCGGAAGGGCCGATCATAAAGATCTGGCTGTAGCGTGGCTCAATACCGCGCACCCAGTAGTTACGCGGCGAGAGGGTGATGCGTTTGCCTTGCTCCAGGTAGCTGc                                                                                                                                                                                                      >  2:130075/1‑137 (MQ=255)
                                                                                   aaaTAGCGATAGCCCACGCCCACTTCGTGCGCGGAAGGGCCGATCATAAAGATCTGGCTGTAGCGTGGCTCAATACCGCGCACCCAGTAGTTACGCGGCGAGAGGGTGATGCGTTTGCCTTGCTCCAGGTAgctgct                                                                                                                                                                                                     <  1:130075/137‑1 (MQ=255)
                                                                                             aGCCCACGCCCACTTCGTGCGCGGAAGGGCCGATAATAAAGATCTGGCTGGAGCGTGGCTCAATACCGCGCACCCAGTAGTTACGCGGCGAGAGGGTGATGCGTTTGCCTTGCTCCAGGTAGCTGCTGCGCAGGGTTTGGGTGTAGAACCCCTGAATGTTGAATTTATGCTGGCTGTCTGGCTGGAACTGGTAGCCCAGGCTCGCCAGCTTGCGACGACCCCAGAAGCGGTCATACGGGCGGGTGGATTGCCAGCGATCGGCGTCGTAATcc                                                    <  1:368202/272‑1 (MQ=255)
                                                                                                            cGTGCGCGGAAGGGCCGATCATAAAGATCTGGCTGTAGCGTGGCTCAATACCGCGCACCCAGTAGTTACGCGGCGAGAGGGTGATGCGTTTGCCTTGCTCCAGGTAGCTGCTGCGCAGGGTTTGGgtgt                                                                                                                                                                                    >  1:87489/1‑129 (MQ=255)
                                                                                                             gTGCGCGGAAGGGCCGATCATAAAGATCTGGCTGTAGCGTGGCTCAATACCGCGCACCCAGTAGTTACGCGGCGAGAGGGTGATGCGTTTGCCTTGCTCCAGGTAGCTGCTGCGCAGGGTTTGGGTGTa                                                                                                                                                                                   <  2:87489/129‑1 (MQ=255)
                                                                                                             gTGCGCGGAAGGGCCGATCATAAAGATCTGGCTGTAGCGTGGCTCAATACCGCGCACCCAGTAGTTACGCGGCGAGAGGGTGATGCGTTTGCCTTGCTCCAGGTAGCTGCTGCGCAGGGTTTGGGTGTAGAACCCCTGAATGTTGAATTTATGCTGGCTGTCTGGCTGGAACTGGTAGCCCAGGCTCGCCAGCTTGCGACGACCCCAGAAGCGGTCATAc                                                                                        <  1:166529/220‑1 (MQ=255)
                                                                                                                     aaGGGCCGATCATAAAGATCTGGCTGTAGCGTGGCTCAATACCGCGCACCCAGTAGTTACGCGGCGAGAGGGTGATGCGTTTGCCTTGCTCCAGGTAGCTGCTGCGCAGGGTTTGGGTGTAGAACCCCTGAATGTTGAATTTATGCTGGCTGTCTGGCTGGAACTGGTAGCCCAGGCTCGCCAGCTTGCGACGACCCCAg                                                                                                    >  1:346462/1‑200 (MQ=255)
                                                                                                                      aGGGCCGATCATAAAGATCTGGCTGTAGCGTGGCTCAATACCGCGCACCCAGTAGTTACGCGGCGAGAGGGTGATGCGTTTGCCTTGCTCCAGGTAGCTGCTGCGCAGGGTTTGGGTGTAGAACCCCTGAATGTTGAATTTATGCTGGCTGTCTGGCTGGAACTGGTAGCCCAGGCTCGCCAGCTTGCGACGACCCCAGa                                                                                                   <  2:346462/200‑1 (MQ=255)
                                                                                                                                           gCTGTAGCGTGGCTCAATACCGCGCACCCAGTAGTTACGCGGCGAGAGGGCGATGCGTTTGCCTTGCTCCAGGTAGCTGCTGCGCAGGGTTTGGGTGTAGAACCCCTGAATGTTGAATTTATGCTGGCTGTCTGGCTGGAACTGGTAGCCCAGGCTCGCCAGCTTGCGACGACCCCAGAAGCGGTCATACGGGCGGGTGGATTGCCAGCg                                                                    >  1:318521/1‑210 (MQ=255)
                                                                                                                                            cTGTAGCGTGGCTCAATACCGCGCACCCAGTAGTTACGCGGCGAGAGGGTGATGCGTTTGCCTTGCTCCAGGTAGCTGCTGCGCAGGGTTTGGGTGTAGAACCCCTGAATGTTGAATTTATGCTGGCTGTCTGGCTGGAACTGGTAGCCCAGGCTCGCCAGCTTGCGACGACCCCAGAAGCGGTCATACGGGCGGGTGGATTGCCAGCGa                                                                   <  2:318521/210‑1 (MQ=255)
                                                                                                                                               tAGCGTGTCTCAATCCCGCGCACCCAGTAGTTACGGGGCGGGGGGGTGATGCTTTGGCCTTGCTCCAGGTAGCTGCTGCGCAGGGTTTGGGTGTAGAACCCCTGATTTTTGAATTTATGCTGGCTGTCTGGCGGGAACTGGTAGCCCAGGCTCGCCAGCTTGCGACGACCCCAGAAGCGGTCATACGGGCGGGTGGATTGCCAGCGATCGGCGTCGTAATCCGCGCGAGACAGGCCACCGGGCATGTCGGCTTCACCGTCGTAATATTgcag  <  1:87554/272‑1 (MQ=255)
                                                                                                                                                                    aCCCAGTAGTTACGCGGCGAGAGGGTGATGCGTTTGCCTTGCTCCAGGTAGCTGTTGCGCAGGGTTTGGGTGTAGAACCCCTGAATGTTGAATTTATGCTGGCTGTCTGGCTGGAACTGGTAGCCCAGGCTCGCCAGCTTGCGACGACCCCAGAAGCGGTCATACGGGCGGGTGGATTGCCAGCGa                                                                   >  2:301422/1‑186 (MQ=255)
                                                                                                                                                                    aCCCAGTAGTTACGCGGCGAGAGGGTGATGCGTTTGCCTTGCTCCAGGTAGCTGCTGCGCAGGGTTTGGGTGTAGAACCCCTGAATGTTGAATTTATGCTGGCTGTCTGGCTGGAACTGGTAGCCCAGGCTCGCCAGCTTGCGACGACCCCAGAAGCGGTCATACGGGCGGGTGGATTGCCAGCGATCGgcggcg                                                          >  1:9657/1‑195 (MQ=255)
                                                                                                                                                                    aCCCAGTAGTTACGCGGCGAGAGGGTGATGCGTTTGCCTTGCTCCAGGTAGCTGCTGCGCAGGGTTTGGGTGTAGAACCCCTGAATGTTGAATTTATGCTGGCTGTCTGGCTGGAACTGGTAGCCCAGGCTCGCCAGCTTGCGACGACCCCAGAAGCGGTCATACGGGCGGGTGGATTGCCAGCGATCGGCGTCGTAATc                                                     >  2:393598/1‑200 (MQ=255)
                                                                                                                                                                     cccAGTAGTTACGCGGCGAGAGGGTGATGCGTTTGCCTTGCTCCAGGTAGCTGCTGCGCAGGGTTTGGGTGTAGAACCCCTGAATGTTGAATTTATGCTGGCTGTCTGGCTGGAACTGGTAGCCCAGGCTCGCCAGCTTGCGACGACCCCAGAAGCGGTCATACGGGCGGGTGGATTGCCAGCGATCGGCGTCGTAATcc                                                    <  1:393598/200‑1 (MQ=255)
                                                                                                                                                                     cccAGTAGTTACGCGGCGAGAGGGTGATGCGGTTGCCTTGCTCCAGGTAGCTGTTGCGCAGGGTTTGGGTGTAGAACCCCTGAATGTTGAATTTATGCTGGCTGTCTGGCTGGAACTGGTAGCCCAGGCTCGCCAGCTTGCGACGACCCCAGAAGCGGTCATACGGGCGGGTGGATTGCCAGCGAt                                                                  <  1:301422/186‑1 (MQ=255)
                                                                                                                                                                                       agagGGTGATGCGTTTGCCTTGCTCCAGGTAGCTGCTGCGCAGGGTTTGGGTGTAGAACCCCTGAATGTTGAATTTATGCTGGCTGTCTGGCTGGAACTGGTAGCCCAGGCTCGCCAGCTTg                                                                                                                >  2:160646/1‑122 (MQ=255)
                                                                                                                                                                                        gagGGTGATGCGTTTGCCTTGCTCCAGGTAGCTGCTGCGCAGGGTTTGGGTGTAGAACCCCTGAATGTTGAATTTATGCTGGCTGTCTGGCTGGAACTGGTAGCCCAGGCTCGCCAGCTTGc                                                                                                               <  1:160646/122‑1 (MQ=255)
                                                                                                                                                                                                                        |                                                                                                                                                                                                      
CGCGGTCGTAAGGGCTTGAGCCGGACGGCAACTGCCCGCTGCTGGTGGCGGTGTAGTAACGCATTTCATGCGTTGATTCATTCAAATAGCGATAGCCCACGCCCACTTCGTGCGCGGAAGGGCCGATCATAAAGATCTGGCTGTAGCGTGGCTCAATACCGCGCACCCAGTAGTTACGCGGCGAGAGGGTGATGCGTTTGCCTTGCTCCAGGTAGCCGCTGCGCAGGGTTTGGGTGTAGAACCCCTGAATGTTGAATTTATGCTGGCTGTCTGGCTGGAACTGGTAGCCCAGGCTCGCCAGCTTGCGACGACCCCAGAAGCGGTCATACGGGCGGGTGGATTGCCAGCGATCGGCGTCGTAATCCGCGCGAGACAGGCCACCGGGCATGTCGGCTTCACCGTCGTAATATTGCAG  >  NC_000913/4515264‑4515678

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