Predicted mutation
evidence seq id position mutation annotation gene description
RA NC_000913 4,516,115 C→T R188H (CGC→CAC)  fecA ← ferric citrate outer membrane transporter

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009134,516,1150CT92.9% 35.6 / NA 14R188H (CGC→CAC) fecAferric citrate outer membrane transporter
Reads supporting (aligned to +/- strand):  ref base C (0/0);  major base T (7/6);  minor base G (0/1);  total (7/7)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

CGTACCACGTCAATGGCATCCATGTTGCCGAGCGAAACGGGAGCCAGTGAAAGCTGCGGCTGACCGTAAGGGGCGAAGGGGACGGGGATGCCGTCCATCAGGACGGTCGAGCGGCTGGCGAGGCGCGGGTTCAGGCCCCGGATGCCAAAGTTCATCGCCAGGTCGTGGCTGCCGGTGCCGTTGTTTTCCGGCGCGCTGACGCCAGGGATGCGGTTAAGTACCT  >  NC_000913/4515992‑4516214
                                                                                                                           |                                                                                                   
cgTACCACGTCAATGGCATCCATGTTGCCGAGCGAAACGGGAGCCAGTGAAAGCTGCGGCTGACCGTAAGGGGCGAAGGGGACGGGGATGCCGTCCATCAGGACGGTCGAGCGGCTGGCGAGGTGCGGGTTCa                                                                                            >  2:418376/1‑133 (MQ=255)
  tACCACGTCAATGGCATCCATGTTGCCGAGCGAAACGGGAGCCAGTGAAAGCTGCGGCTGACCGTAAGGGGCGAAGGGGACGGGGATGCCGTCCATCAGGACGGTCGAGCGGCTGGCGAGGTGCGGGTTCAgg                                                                                          >  2:494047/1‑133 (MQ=255)
       cGTCAATGGCATCCATGTTGCCGAGCGAAACGGGAGCCAGTGAAAGCTGCGGCTGACCGTAAGGGGCGAAGGGGACGGGGATGCCGTCCATCAGGACGGTCGAGCGGCTGGCGAGGTGCGGGTTCAGGCCCCg                                                                                     >  2:478904/1‑133 (MQ=255)
           aaTGGCATCCATGTTGCCGAGCGAAACGGGAGCCAGTGAAAGCTGCGGCTGACCGTAAGGGGCGAAGGGGACGGGGATGCCGTCCATCAGGACGGTCGAGCGGCTGGCGAGGTGCGGGTTCAGGCCCCGGATg                                                                                 <  2:414403/133‑1 (MQ=255)
                                            cAGTGAAAGCTGCGGCTGACCGTAAGGGGCGAAGGGGACGGGGATGCCGTCCATCAGGACGGTCGAGCGGCTGGCGAGGTGCGGGTTCAGGCCCCGGATGCCAAAGTTCATCGCCAGGTCGTGGCTGCCGGTg                                                <  1:494047/133‑1 (MQ=255)
                                                   aGCTGCGGCTGACCGTAAGGGGCGAAGGGGACGGGGATGCCGTCCATCAGGACGGTCGAGCGGCTGGCGAGGTGCGGGTTCAGGCCCCGGATGCCAAAGTTCATCGCCAGGTCGTGGCTGCCGGTGCCgttgt                                         >  2:363208/1‑133 (MQ=255)
                                                         ggCTGACCGTAAGGGGCGAAGGGGACGGGGATGCCGTCCATCAGGACGGTCGAGCGGCTGGCGAGGTGCGGGTTCAGGTCCCGGATGCCAAAGTTCATCGCCAGGTCGTGGCGGGCGGTGCCGTGCTTTTCCg                                   >  2:347242/1‑133 (MQ=255)
                                                                        gCGAAGGGGACGGGGATGCCGTCCATCAGGACGGTCGAGCGGCTGGCGAGGTGCGGGTTCAGGCCCCGGATGCCAAAGTTCATCGCCAGGTCGTGGCTGCCGGTGCCGTTGTTTTCCGGCGCGCTGACGCCAg                    >  2:24814/1‑133 (MQ=255)
                                                                             ggggACGGGGATGCCGTCCATCAGGACGGTCGAGCGGCTGGCGAGGTGCGGGTTCAGGCCCCGGATGCCAAAGTTCATCGCCAGGTCGTGGCTGCCGGTGCCGTTGTTTTCCGGCGCGCTGACGCCAGGGATg               <  1:151123/133‑1 (MQ=255)
                                                                             ggggACGGGGATGCCGTCCATCAGGACGGTCGAGCGGCTGGCGAGGTGCGGGTTCAGGCCCCGGATGCCAAAGTTCATCGCCAGGTCGTGGCTGCCGGTGCCGTTGTTTTCCGGCGCGCTGACGCCAGGGATg               <  1:418376/133‑1 (MQ=255)
                                                                                  cGGGGATGCCGTCCATCAGGACGGACGATCGGTTGGCGAGGGGCGGGTTCAGGCCCCGGATGCCAAAGTTCATCGCCAGGTCGTGGCTGCCGGTGCCGTTGTTTTCCGGCGCGCTGACGCCAGGGATGCGGtt          <  2:176044/133‑1 (MQ=255)
                                                                                    gggATGCCGTCCATCAGGACGGTCGAGCGGCTGGCGAGGTGCGGGTTCAGGCCCCGGATGCCAAAGTTCATCGCCAGGTCGTGGCTGCCGGTGCCGTTGTTTTCCGGCGCGCTGACGCCAGGGATGCGGTTaa        <  1:478904/133‑1 (MQ=255)
                                                                                        tGCCGTCCATCAGGACGGTCGAGCGGCTGGCGAGGTGCGGGTTCAGGCCCCGGATGCCAAAGTTCATCGCCAGGTCGTGGCTGCCGGTGCCGTTGTTTTCCGGCGCGCTGACGCCAGGGATGCGGTTaa        <  1:308213/129‑1 (MQ=255)
                                                                                        tGCCGTCCATCAGGACGGTCGAGCGGCTGGCGAGGTGCGGGTTCAGGCCCCGGATGCCAAAGTTCATCGCCAGGTCGTGGCTGCCGGTGCCGTTGTTTTCCGGCGCGCTGACGCCAGGGATGCGGTTaa        >  2:308213/1‑129 (MQ=255)
                                                                                          ccGTCCATCAGGACGGTCGAGCGGCTGGCGAGGTGCGGGTTCAGGCCCCGGATGCCAAAGTTCATCGCCAGGTCGTGGCTGCCGGTGCCGTAGTTTTCCGGCGCCCTGCCGCCAGGGATGCGGTTAAGTACCt  >  2:273898/1‑133 (MQ=255)
                                                                                                                           |                                                                                                   
CGTACCACGTCAATGGCATCCATGTTGCCGAGCGAAACGGGAGCCAGTGAAAGCTGCGGCTGACCGTAAGGGGCGAAGGGGACGGGGATGCCGTCCATCAGGACGGTCGAGCGGCTGGCGAGGCGCGGGTTCAGGCCCCGGATGCCAAAGTTCATCGCCAGGTCGTGGCTGCCGGTGCCGTTGTTTTCCGGCGCGCTGACGCCAGGGATGCGGTTAAGTACCT  >  NC_000913/4515992‑4516214

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