Predicted mutation
evidence seq id position mutation annotation gene description
RA NC_000913 4,516,827 C→T A297A (GCG→GCA fecR ← regulator for fec operon, periplasmic

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009134,516,8270CT96.1% 160.1 / ‑1.4 51A297A (GCG→GCAfecRregulator for fec operon, periplasmic
Reads supporting (aligned to +/- strand):  ref base C (1/1);  new base T (34/15);  total (35/16)
Fisher's exact test for biased strand distribution p-value = 5.33e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 6.30e-01

TTTTTGTTGTGTTCAGCTATGAGTCGAACGAGAAGGGTAAAAAGGACAATCGAAATAAGAATTATTTTCCTTACAGTGGTGAAATGTTTATCCAGTACCGCGTAATAGACTGAATTTTAACGGGAAGCGTTTGCGCGATAACGTTCAGGATCGCATCGGTATTTTTCAGCGGGAACGTCCCGCTCAGGCGCAGCCCGGCAACGGCGGGATCGCAGCGCAGCACGCCGTTGCGGTAACGGGTTAGCGTGGCTATCACCTCACCCAGC  >  NC_000913/4516694‑4516959
                                                                                                                                     |                                                                                                                                    
tttttGTTGTGTTCAGCTATGAGTCGAACGAGAAGGGTAAAAAGGACAATCGAAATAAGAATTATTTTCCTTACAGTGGTGAAATGTTTATCCAGTACCGCGTAATAGACTGAATTTTAACGGGAAGCGTTTGTGCGATa                                                                                                                                <  1:673494/140‑1 (MQ=255)
           ttCAGCTATGAGTCGAACGAGAAGGGTAAAAAGGACAATCGAAATAAGAATTATTTTCCTTACAGTGGTGAAATGTTTATCCAGTACCGCGTAATAGACTGAATTTTAACGGGAAGCGTTTGTGCGATAACGTTCAGGAt                                                                                                                     <  1:422218/140‑1 (MQ=255)
                   tGAGTCGAACGAGAAGGGTAAAAAGGACAATCGAAATAAGAATTATTTTCCTTACAGTGGTGAAATGTTTATCCAGTACCGCGTAATAGACTGAATTTTAACGGGAAGCGTTTGTGCGATAACGTTCAGGATCGCATCgg                                                                                                             <  1:541592/140‑1 (MQ=255)
                          aaCGAGAAGGGTAAAAAGGACAATCGAAATAAGAATTATTTTCCTTACAGTGGTGAAATGTTTATCCAGTACCGCGTAATAGACTGAATTTTAACGGGAAGCGTTTGTGCGATAACGTTCAGGATCGCATCGGTAttttt                                                                                                      >  2:623876/1‑140 (MQ=255)
                          aaCGAGAAGGGTAAAAAGGACAATCGAAATAAGAATTATTTTCCTTACAGTGGTGAAATGTTTATCCAGTACCGCGTAATAGACTGAATTTTAACGGGAAGCGTTTGCGCGATAACGTTCAGGATCGCATCGGTAttttt                                                                                                      >  1:171159/1‑140 (MQ=255)
                                     tAAAAAGGACAATCGAAATAAGAATTATTTTCCTTACAGTGGTGAAATGTTTATCCAGTACCGCGTAATAGACTGAATTTTAACGGGAAGCGTTTGTGCGATAACGTTCAGGATCGCATCGGTATTTTTCAGCGGGAACg                                                                                           <  1:563159/140‑1 (MQ=255)
                                             aCAATCGAAATAAGAATTATTTTCCTTACAGTGGTGAAATGTTTATCCAGTACCGCGTAATAGACTGAATTTTAACGGGAAGCGTTTGTGCGATAACGTTCAGGATCGCATCGGTATTTTTCAGCGGGAACGTCCCGCTc                                                                                   >  1:213878/1‑140 (MQ=255)
                                             aCAATCGAAATAAGAATTATTTTCCTTACAGTGGTGAAATGTTTATCCAGTACCGCGTAATAGACTGAATTTTAACGGGAAGCGTTTGTGCGATAACGTTCAGGATCGCATCGGTATTTTTCAGCGGGAACGTCCCGCTc                                                                                   >  1:447315/1‑140 (MQ=255)
                                             aCAATCGAAATAAGAATTATTTTCCTTACAGTGGTGAAATGTTTATCCAGTACCGCGTAATAGACTGAATTTTAACGGGAAGCGTTTGTGCGATAACGTTCAGGATCGCATCGGTATTTTTCAGCGGGAACGTCCCGCTc                                                                                   <  1:869965/140‑1 (MQ=255)
                                             aCAATCGAAATAAGAATTATTTTCCTTACAGTGGTGAAATGTTTATCCAGTACCGCGTAATAGACTGAATTTTAACGGGAAGCGTTTGTGCGATAACGTTCAGGATCGCATCGGTATTTTTCAGCGGGAACGTCCCGCTc                                                                                   >  2:78082/1‑140 (MQ=255)
                                             aCAATCGAAATAAGAATTATTTTCCGTACAGTGGTGAAATGTTTATCCAGTACCGCGTAATAGACTGAATTTTAACGGGAAGCGTTTGTGCGATAACGTTCAGGATCGCATCGGTATTTTTCAGCGGGAACGTCCCGCTc                                                                                   >  1:862461/1‑140 (MQ=255)
                                              cAATCGAAATAAGAATTATTTTCCTTACAGTGGTGAAATGTTTATCCAGTACCGCGTAATAGACTGAATTTTAACGGGAAGCGTTTGTGCGATAACGTTCAGGATCGCATCGGTATTTTTCAGCGGGAACGTCCCGCTCa                                                                                  >  2:933964/1‑140 (MQ=255)
                                                  cGAAATAAGAATTATTTTCCTTACAGTGGTGAAATGTTTATCCAGTACCGCGTAATAGACTGAATTTTAACGGGAAGCGTTTGTGCGATAACGTTCAGGATCGCATCGGTATTTTTCAGCGGGAACGTCCCGCTCAGgcg                                                                              <  2:95584/140‑1 (MQ=255)
                                                  cGAAATAAGAATTATTTTCCTTACAGTGGTGAAATGTTTATCCAGTACCGCGTAATAGACTGAATTTTAACGGGAAGCGTTTGTGCGATAACGTTCAGGATCGCATCGGTATTTTTCAGCGGGAACGTCCCGCTCAGgcg                                                                              <  2:503939/140‑1 (MQ=255)
                                                                   tCCTTACAGTGGTGAAATGTTTATCCAGTACCGCGTAATAGACTGAATTTTAACGGGAAGCGTTTGTGCGATAACGTTCAGGATCGCATCGGTATTTTTCAGCGGGAACGTCCCGCTCAGGCGCAGCCCGGCAAcggcgg                                                             <  1:539918/140‑1 (MQ=255)
                                                                        aCAGTGGTGAAATGTTTATCCAGTACCGCGTAATAGACTGAATTTTAACGGGAAGCGTTTGTGCGATAACGTTCAGGATCGCATCGGTATTTTTCAGCGGGAACGTCCCGCTCAGGCGCAGCCCGGCAACGGCGGGATcg                                                        >  1:208272/1‑140 (MQ=255)
                                                                              gtgAAATGTTTATCCAGTACCGCGTAATAGACTGAATTTTAACGGGAAGCGTTTGTGCGATAACGTTCAGGATCGCATCGGTATTTTTCAGCGGGAACGTCCCGCTCAGGCGCAGCCCGGCAACGGCGGGATcgcagcgc                                                  >  2:916677/1‑140 (MQ=255)
                                                                              gtgAAATGTTTATCCAGTACCGCGTAATAGACTGAATTTTAACGGGAAGCGTTTGTGCGATAACGTTCAGGATCGCATCGGTATTTTTCAGCGGGAACGTCCCGCTCAGGCGCAGCCCGGCAACGGCGGGATcgcagcgc                                                  >  1:80391/1‑140 (MQ=255)
                                                                               tgAAATGTTTATCCAGTACCGCGTAATAGACTGAATTTTAACGGGAAGCGTTTGTGCGATAACGTTCAGGATCGCATCGGTATTTTTCAGCGGGAACGTCCCGCTCAGGCGCAGCCCGGCAACGGCGGGATcgcagcgca                                                 >  2:309981/1‑140 (MQ=255)
                                                                                  aaTGTTTATCCAGTACCGCGTAATAGACTGAATTTTAACGGGAAGCGTTTGTGCGATAACGTTCAGGATCGCATCGGTATTTTTCAGCGGGAACGTCCCGCTCAGGCGCAGCCCGGCAACGGCGGGATCGCAGCgcagca                                              <  1:594254/140‑1 (MQ=255)
                                                                                     gTTTATCCAGTACCGCGTAATAGACTGAATTTTAACGGGAAGCGTTTGTGCGATAACGTTCAGGATCGCATCGGTATTTTTCAGCGGGAACGTCCCGCTCAGGCGCAGCCCGGCAACGGCGGGATCGCAGCGCAGCACGc                                           >  1:502530/1‑140 (MQ=255)
                                                                                      tttATCCAGTACCGCGTAATAGACTGAATTTTAACGGGAAGCGTTTGTGCGATAGCGTTCAGGATCGCATCGGTATTTTTCAGCGGGAACGTCCCGCTCAGGCGCAGCCCGGCAACGGCGGGATCGCAGCGCAGCACGcc                                          >  2:758863/1‑140 (MQ=255)
                                                                                      tttATCCAGTACCGCGTAATAGACTGAATTTTAACGGGAAGCGTTTGTGCGATAACGTTCAGGATCGCATCGGTATTTTTCAGCGGGAACGTCCCGCTCAGGCGCAGCCCGGCAACGGCGGGATCGCAGCGCAGCACGcc                                          >  2:21403/1‑140 (MQ=255)
                                                                                      tttATCCAGTACCGCGTAATAGACTGAATTTTAACGGGAAGCGTTTGTGCGATAACGTTCAGGATCGCATCGGTATTTTTCAGCGGGAACGTCCCGCTCAGGCGCAGCCCGGCAACGGCGGGATCGCAGCGCAGCACGcc                                          >  1:702632/1‑140 (MQ=255)
                                                                                      tttATCCAGTACCGCGTAATAGACTGAATTTTAACGGGAAGCGTTTGTGCGATAACGTTCAGGATCGCATCGGTATTTTTCAGCGGGAACGTCCCGCTCAGGCGCAGCCCGGCAACGGCGGGATCGCAACGCAGCAcccc                                          >  2:744457/1‑140 (MQ=255)
                                                                                                aCCGCGTAATAGACTGAATTTTAACGGGAAGCGTTTGTGCGATAACGTTCAGGATCGCATCGGTATTTTTCAGCGGGAACGTCCCGCTCAGGCGCAGCCCGGCAACGGCGGGATCGCAGCGCAGCACGCCGTTTCGGTaa                                >  1:279353/1‑140 (MQ=255)
                                                                                                aCCGCGTAATAGACTGAATTTTAACGGGAAGCGTTTGTGCGATAACGTTCAGGATCGCATCGGTATTTTTCAGCGGGAACGTCCCGCTCAGGCGCAGCCCGGCAACGGCGGGATCGCAGCGCAGCACGCCGTTGCGGTaa                                >  1:736305/1‑140 (MQ=255)
                                                                                                aCCGCGTAATAGACTGAATTTTAACGGGAAGCGTTTGTGCGATAACGTTCAGGATCGCATCGGTATTTTTCAGCGGGAACGTCCCGCTCAGGCGCAGCCCGGCAACGGCGGGATCGCAGCCCCGCACGCCGTGGCGGGaa                                >  2:210381/1‑140 (MQ=255)
                                                                                                  cgcgTAATAGACTGAATTTTAACGGGAAGCGTTTGTGCGATAACGTTCAGGATCGCATCGGTATTTTTCAGCGGGAACGTCCCGCTCAGGCGCAGCCCGGCAACGGCGGGACTGCAGCGGAGCACGACGTTGCggggaac                              >  2:726536/1‑135 (MQ=255)
                                                                                                   gcgTAATAGACTGAATTTTAACGGGAAGCGTTTGTGCGATAACGTTCAGGATCGCATCGGTATTTTTCAGCGGGAACGTCCCGCTCAGGCGCAGCCCGGCAACGGCGGGATCGCAGCGCAGCACGCCGTTGCGGTAACgg                             <  1:933964/140‑1 (MQ=255)
                                                                                                    cgTAATAGACTGAATTTTAACGGGAAGCGTTTGTGCGATAACGTTCAGGATCGCATCGGTATTTTTCAGCGGGAACGTCCCGCTCAGGCGCAGCCCGGCAACGGCGGGATCGCAGCGCAGCACGCCGTTGCGGTAACggg                            >  2:998612/1‑140 (MQ=255)
                                                                                                    cgTAATAGACTGAATTTTAACGGGAAGCGTTTGTGCGATAACGTTCAGGATCGCATCGGTATTTTTCAGCGGGAACGTCCCGCTCAGGCGCAGCCCGGCAACGGCGGGATCGCAGCGCAGCACGCCGTTGCGGTAACggg                            >  2:56046/1‑140 (MQ=255)
                                                                                                    cgTAATAGACTGAATTTTAACGGGAAGCGTTTGTGCGATAACGTTCAGGATCGCATCGGTATTTTTCAGCGGGAACGTCCCGCTCAGGCGCAGCCCGGCAACGGCGGGATCGCAGCGCAGCACGCCGTTGCGGTAACggg                            >  1:108822/1‑140 (MQ=255)
                                                                                                    cgTAATAGACTGAATTTTAACGGGAAGCGTTTGTGCGATAACGTTCAGGATCGCATCGGTATTTTTCAGCGGGAACGTCCCGCTCAGGCGCAGCCCGGCAACGGCGGGATCGCAGCGCAGCACGCCGTTGCGGTAACggg                            >  1:909309/1‑140 (MQ=255)
                                                                                                    cgTAATAGACTGAATTTTAACGGGAAGCGTTTGTGCGATAACGTTCAGGATCGCATCGGTATTTTTCAGCGGGAACGTCCCGCTCAGGCGCAGCCCGGCAACGGCGGGATCGCAGCGCAGCACGCCGTTGCGGTAACggg                            >  1:68799/1‑140 (MQ=255)
                                                                                                    cgTAATAGACTGAATTTTAACGGGAAGCGTTTGTGCGATAACGTTCAGGATCGCATCGGTATTTTTCAGCGGGAACGTCCCGCTCAGGCGCAGCCCGGCAACGGCGGGATCGCAGCGCAGCACGCCGTTGCGGTAACggg                            >  1:700366/1‑140 (MQ=255)
                                                                                                    cgTAATAGACTGAATTTTAACGGGAAGCGTTTGTGCGATAACGTTCAGGATCGCATCGGTATTTTTCAGCGGGAACGTCCCGCTCAGGCGCAGCCCGGCAACGGCGGGATCGCAGCGCAGCACGCCGTTGCGGGAACggg                            >  2:799729/1‑140 (MQ=255)
                                                                                                    cgTAATAGACTGAATTTTAACGGGAAGCGTTTGTGCGATAACGTTCAGGATCGCATCGGTATTTTTCAGCGGGAACGTCCCGCTCAGGCGCAGCCCGGCAACGGCGGGATCGCAGCGCAGCACGCCGGTGCGGTAACggg                            >  2:826733/1‑140 (MQ=255)
                                                                                                    cgTAATAGACTGAATTTTAACGGGAAGCGTTTGTGCGATAACGTTCAGGATCGCATCGGTATTTTTCAGCGGGAACGTCCCGCTCAGGCGCAGCCCGGCAACGGCGGGATCGCAGCGCAGCACGCCGGTGCGGTAACggg                            >  2:782509/1‑140 (MQ=255)
                                                                                                    cgTAATAGACTGAATTTTAACGGGAAGCGTTTGTGCGATAACGTTCAGGATCGCATCGGTATTTTTCAGCGGGAACGTCCCGCTCAGGCGCAGCCCGGCAACGGCGGGATCGCAGCGCAGCACGCCGGTGCGGTAACggg                            >  2:748641/1‑140 (MQ=255)
                                                                                                    cgTAATAGACTGAATTTTAACGGGAAGCGTTTGTGCGATAACGTTCAGGATCGCATCGGTATTTTTCAGCGGGAACGTCCCGCTCAGGCGCAGCCCGGCAACGGCGGGATCGCAGCGCAGCACGCCGGTGCGGGAACggg                            >  2:343318/1‑140 (MQ=255)
                                                                                                    cgTAATAGACTGAATTTTAACGGGAAGCGTTTGTGCGATAACGTTCAGGATCGCATCGGTATTTTTCAGCGGGAACGTCCCGCTCAGGCGCAGCCCGGCAACGGCGGGATCGCAGCGCAGCACGCCCTTGCGGTAACggg                            >  1:611036/1‑140 (MQ=255)
                                                                                                          aGACTGAATTTTAACGGGAAGCGTTTGTGCGATAACGTTCAGGATCGCATCGGTATTTTTCAGCGGGAACGTCCCGCTCAGGCGCAGCCCGGCAACGGCGGGATCGCAGCGCAGCACGCCGTTGCGGTAACGGGTTAGCg                      <  1:799729/140‑1 (MQ=255)
                                                                                                          aGACTGAATTTTAACGGGAAGCGTTTGTGCGATAACGTTCAGGATCGCATCGGTATTTTTCAGCGGGAACGTCCCGCTCAGGCGCAGCCCGGCAACGGCGGGATCGCAGCGCAGCACGCCGTTGCGGTAACGGGTTAGCg                      <  1:289623/140‑1 (MQ=255)
                                                                                                           gACTGAATTTTAACGGGAAGCGTTTGTGCGATAACGTTCAGGATCGCATCGGTATTTTTCAGCGGGAACGTCCCGCTCAGGCGCAGCCCGGCAACGGCGGGATCGCAGCGCAGCACGCCGTTGCGGTAACGGGTTAGCGt                     <  1:202804/140‑1 (MQ=255)
                                                                                                            aCTGAATTTTAACGGGAAGCGTTTGTGCGATAACGTTCAGGATCGCATCGGTATTTTTCAGCGGGAACGTCCCGCTCAGGCGCAGCCCGGCAACGGCGGGATCGCAGCGCGGCACGCCGTTGCGGTAAAGGGTTTGCGTg                    >  1:601722/1‑140 (MQ=255)
                                                                                                            aCTGAATTTTAACGGGAAGCGTTTGTGCGATAACGTTCAGGATCGCATCGGTATTTTTCAGCGGGAACGTCCCGCTCAGGCGCAGCCCGGCAACGGCGGGATCGCAGCGCAGCACGCCGTTGCGGTAACGGGTTAGCGTg                    >  1:858189/1‑140 (MQ=255)
                                                                                                                aaTTTTAACGGGAAGCGTTTGTGCGATAACGTTCAGGATCGCATCGGTATTTTTCAGCGGGAACGTCCCGCTCAGGCGCAGCCCGGCAACGGCGGGATCGCAGCGCAGCACGCCGTTGCGGTAACGGGTTAGCGTGGCTa                <  1:826733/140‑1 (MQ=255)
                                                                                                                aaTTTTAACGGGAAGCGTTTGCGCGATAACGTTCAGGATCGCATCGGTATTTTTCAGCGGGAACGTCCCGCTCAGGCGCAGCCCGGCAACGGCGGGATCGCAGCGCAGCACGCCGTTGCGGTAACGGGTTAGCGTGGCTa                <  2:171159/140‑1 (MQ=255)
                                                                                                                            aaGCGTTTGTGCGATAACGTTCAGGATCGCATCGGTATTTTTCAGCGGGAACGTCCCGCTCAGGCGCAGCCCGGCAACGGCGGGGTCGCCGCGCAGCACGCCGTTGCGGTAACGGGGTAACGTGGCTATcaccccacccc    >  2:597456/1‑139 (MQ=255)
                                                                                                                              gCGTTTGTGCGATAACGTTCAGGATCGCATCGGTATTTTTCAGCGGGAACGTCCCGCTCAGGCGCAGCCCGGCAACGGCGGGATCGCAGCGCAGCACGCCGTTGCGGTAACGGGTTAGCGTGGCTATCACCTCACCCAGc  <  1:722180/140‑1 (MQ=255)
                                                                                                                                     |                                                                                                                                    
TTTTTGTTGTGTTCAGCTATGAGTCGAACGAGAAGGGTAAAAAGGACAATCGAAATAAGAATTATTTTCCTTACAGTGGTGAAATGTTTATCCAGTACCGCGTAATAGACTGAATTTTAACGGGAAGCGTTTGCGCGATAACGTTCAGGATCGCATCGGTATTTTTCAGCGGGAACGTCCCGCTCAGGCGCAGCCCGGCAACGGCGGGATCGCAGCGCAGCACGCCGTTGCGGTAACGGGTTAGCGTGGCTATCACCTCACCCAGC  >  NC_000913/4516694‑4516959

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