Predicted mutation
evidence seq id position mutation freq annotation gene description
RA minE 2,679,619 A→G 44.8% I357I (ATT→ATC feoB ← fused ferrous iron transporter, protein B

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*minE2,679,6190AG44.8% 13.7 / 70.8 87I357I (ATT→ATCfeoBfused ferrous iron transporter, protein B
Reads supporting (aligned to +/- strand):  ref base A (37/11);  new base G (13/26);  total (50/37)
Fisher's exact test for biased strand distribution p-value = 6.75e-05
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

ACAAACGCCGCACGCGCCATATACCCGGAGTCCTCAAGGAAGGAGAGGAACAGGTACATCATGCCAATCTGCGGCACCAGTGGCAGCACGGTGTTAATGCCGCCACCCAGGCCCTGGGCGAGGAAGAT  >  minE/2679554‑2679681
                                                                 |                                                              
acaAACGCCGCACGCGCCATATACCCGGAGTCCTCAAGGAAGGAGAGGAACAGGTACATCATGCCAATc                                                             >  1:261505/1‑69 (MQ=255)
acaAACGCCGCACGCGCCATATACCCGGAGTCCTCAAGGAAGGAGAGGAACAGGTACATCATGCCAATc                                                             >  1:838782/1‑69 (MQ=255)
acaAACGCCGCACGCGCCATATACCCGGAGTCCTCAAGGAAGGAGAGGAACAGGTACATCATGCCAATc                                                             >  1:711699/1‑69 (MQ=255)
acaAACGCCGCACGCGCCATATACCCGGAGTCCTCAAGGAAGGAGAGGAACAGGTACATCATGCCAATc                                                             >  1:2937850/1‑69 (MQ=255)
     cgccgcACGCGCCATATACCCGGAGTCCTCAAGGAAGGAGAGGAACAGGTACATCATGCCAATCTGCgg                                                        <  1:42341/69‑1 (MQ=255)
          cACGCGCCATATAGCCGGAGTCCTCAAGGAAGGAGAGGAACAGGTACATCATGCCGATCTGCGGCACCa                                                   <  1:2186889/69‑1 (MQ=255)
          cACGCGCCATATAGCCGGAGTCCTCAAGGAAGGAGAGGAACAGGTACATCATGCCGATCTGCGGCACCa                                                   <  1:2611129/69‑1 (MQ=255)
          cACGCGCCATATAGCCGGAGTCCTCAAGGAAGGAGAGGAACAGGTACATCATGCCGATCTGCGGCACCa                                                   <  1:295408/69‑1 (MQ=255)
          cACGCGCCATATAGCCGGAGTCCTCAAGGAAGGAGAGGAACAGGTACATCATGCCGATCTGCGGCACCa                                                   <  1:65359/69‑1 (MQ=255)
          cACGCGCCATATAGCCGGAGTCCTCAAGGAAGGAGAGGAACAGGTACATCATGCCGATCTGCGGCACCa                                                   <  1:1867300/69‑1 (MQ=255)
          cACGCGCCATATAGCCGGAGTCCTCAAGGAAGGAGAGGAACAGGTACATCATGCCGATCTGCGGCACCa                                                   <  1:551639/69‑1 (MQ=255)
          cACGCGCCATATAGCCGGAGTCCTCAAGGAAGGAGAGGAACAGGTACATCATGCCGATCTGCGGCACCa                                                   <  1:415054/69‑1 (MQ=255)
          cACGCGCCATATAGCCGGAGTCCTCAAGGAAGGAGAGGAACAGGTACATCATGCCGATCTGCGGCACCa                                                   <  1:1702344/69‑1 (MQ=255)
                 cATATAGCCGGAGTCCTCAAGGAAGGAGAGGAACAGGTACATCATGCCGATCTg                                                           >  1:3274847/1‑54 (MQ=255)
                 cATATAGCCGGAGTCCTCAAGGAAGGAGAGGAACAGGTACATCATGCCGATCTGc                                                          >  1:1588826/1‑55 (MQ=255)
                 cATATAGCCGGAGTCCTCAAGGAAGGAGAGGAACAGGTACATCATGCCGATCTGc                                                          >  1:1656895/1‑55 (MQ=255)
                 cATATAGCCGGAGTCCTCAAGGAAGGAGAGGAACAGGTACATCATGCCGATCTGc                                                          >  1:3127087/1‑55 (MQ=255)
                 cATATAGCCGGAGTCCTCAAGGAAGGAGAGGAACAGGTACATCATGCCGATCTGc                                                          >  1:3076382/1‑55 (MQ=255)
                 cATATAGCCGGAGTCCTCAAGGAAGGAGAGGAACAGGTACATCATGCCGATCTGc                                                          >  1:3072381/1‑55 (MQ=255)
                 cATATAGCCGGAGTCCTCAAGGAAGGAGAGGAACAGGTACATCATGCCGATCTGc                                                          >  1:1404514/1‑55 (MQ=255)
                 cATATAGCCGGAGTCCTCAAGGAAGGAGAGGAACAGGTACATCATGCCGATCTGc                                                          >  1:1023979/1‑55 (MQ=255)
                 cATATAGCCGGAGTCCTCAAGGAAGGAGAGGAACAGGTACATCATGCCGATCTGc                                                          >  1:709508/1‑55 (MQ=255)
                 cATATAGCCGGAGTCCTCAAGGAAGGAGAGGAACAGGTACATCATGCCGATCTGc                                                          >  1:2109449/1‑55 (MQ=255)
                 cATATAGCCGGAGTCCTCAAGGAAGGAGAGGAACAGGTACATCATGCCGATCTGc                                                          >  1:1963994/1‑55 (MQ=255)
                 cATATAGCCGGAGTCCTCAAGGAAGGAGAGGAACAGGTACATCATGCCGATCTGc                                                          >  1:2389807/1‑55 (MQ=255)
                 cATATAGCCGGAGTCCTCAAGGAAGGAGAGCAACAGGTACATCATGCCGATCTGc                                                          >  1:3361197/1‑55 (MQ=255)
                     taCCCGGAGTCCTCAAGGAAGGAGAGGAGCAGGTACATCATGCCAATCTGCGGCACCAGTGGCAGCACg                                        >  1:3134076/1‑69 (MQ=255)
                     taCCCGGAGTCCTCAAGGAAGGAGAGGAACAGGTACATCATGCCa                                                                >  1:2274723/1‑45 (MQ=255)
                     taCCCGGAGTCCTCAAGGAAGGAGAGGAACAGGTACATCATGCCa                                                                >  1:672854/1‑45 (MQ=255)
                     taCCCGGAGTCCTCAAGGAAGGAGAGGAACAGGTACATCATGCCAATCTGCGGCAcc                                                    >  1:2167868/1‑57 (MQ=255)
                     taCCCGGAGTCCTCAAGGAAGGAGAGGAACAGGTACATCATGCCAATCTGCGGCAcc                                                    >  1:991609/1‑57 (MQ=255)
                     taCCCGGAGTCCTCAAGGAAGGAGAGGAACAGGTACATCATGCCAATCTGCGGCACCAGTGGCAGCACg                                        >  1:271626/1‑69 (MQ=255)
                     taCCCGGAGTCCTCAAGGAAGGAGAGGAACAGGTACATCATGCCAATCTGCGGCACCAGTGGCAGCACg                                        >  1:2836649/1‑69 (MQ=255)
                     taCCCGGAGTCCTCAAGGAAGGAGAGGAACAGGTACATCATGCCAATCTGCGGCACCAGTGGCAGCACg                                        >  1:883486/1‑69 (MQ=255)
                     taCCCGGAGTCCTCAAGGAAGGAGAGGAACAGGTACATCATGCCAATCTGCGGCACCAGTGGCAGCACg                                        >  1:922928/1‑69 (MQ=255)
                     taCCCGGAGTCCTCAAGGAAGGAGAGGAACAGGTACATCATGCCAATCTGCGGCACCAGTGGCAGCACg                                        >  1:2452904/1‑69 (MQ=255)
                     taCCCGGAGTCCTCAAGGAAGGAGAGGAACAGGTACATCATGCCAATCTGCGGCACCAGTGGCAGCACg                                        >  1:2260473/1‑69 (MQ=255)
                     taCCCGGAGTCCTCAAGGAAGGAGAGGAACAGGTACATCATGCCAATCTGCGGCACCAGTGGCAGCACg                                        >  1:1039580/1‑69 (MQ=255)
                     taCCCGGAGTCCTCAAGGAAGGAGAGGAACAGGTACATCATGCCAATCTGCGGCACCAGTGGCAGCACg                                        >  1:1100477/1‑69 (MQ=255)
                     taCCCGGAGTCCTCAAGGAAGGAGAGGAACAGGTACATCATGCCAATCTGCGGCACCAGTGGCAGCACg                                        >  1:1935720/1‑69 (MQ=255)
                     taCCCGGAGTCCTCAAGGAAGGAGAGGAACAGGTACATCATGCCAATCTGCGGCACCAGTGGCAGCACg                                        >  1:2265745/1‑69 (MQ=255)
                     taCCCGGAGTCCTCAAGGAAGGAGAGGAACAGGTACATCATGCCAATCTGCGGCACCA‑‑GGCAGCACGgt                                      >  1:2880284/1‑69 (MQ=255)
                       gccGGAGTCCTCAAGGAAGGTGAGGAACAGGTACATCATGCCGATCTGCGGCACCAGTGGCAGCACgg                                       <  1:1155594/67‑1 (MQ=255)
                       gccGGAGTCCTCAAGGAAGGATAGGAACAGGTACATCATGCCGATCTGCGGCACCAGTGGCAGCACgg                                       <  1:1715950/67‑1 (MQ=255)
                       gccGGAGTCCTCAAGGAAGGAGAGGAACATGTACATCATGCCGATCTGCGGCACCAGTGGCAGCACgg                                       <  1:689153/67‑1 (MQ=255)
                       gccGGAGTCCTCAAGGAAGGAGAGGAACAGGTACATCATGCCGATCTGCGGCACCAGTGGCAGCACgg                                       <  1:1409541/67‑1 (MQ=255)
                       gccGGAGTCCTCAAGGAAGGAGAGGAACAGGTACATCATGCCGATCTGCGGCACCAGTGGCAGCACgg                                       <  1:3333606/67‑1 (MQ=255)
                       gccGGAGTCCTCAAGGAAGGAGAGGAACAGGTACATCATGCCGATCTGCGGCACCAGTGGCAGCACgg                                       <  1:1555959/67‑1 (MQ=255)
                       gccGGAGTCCTCAAGGAAGGAGAGGAACAGGTACATCATGCCGATCTGCGGCACCAGTGGCAGCACgg                                       <  1:739573/67‑1 (MQ=255)
                       gccGGAGTCCTCAAGGAAGGAGAGGAACAGGTACATCATGCCGATCTGCGGCACCAGTGGCAGCACgg                                       <  1:184939/67‑1 (MQ=255)
                       gccGGAGTCCTCAAGGAAGGAGAGGAACAGGTACATCATGCCGATCTGCGGCACCAGTGGCAGCACgg                                       <  1:2875851/67‑1 (MQ=255)
                       gccGGAGTCCTCAAGGAAGGAGAGGAACAGGTACATCATGCCGATCTGCGGCACCAGTGGCAGCACgg                                       <  1:2744203/67‑1 (MQ=255)
                       gccGGAGTCCTCAAGGAAGGAGAGGAACAGGTACATCATGCCGATCTGCGGCACCAGTGGCAGCACgg                                       <  1:2836033/67‑1 (MQ=255)
                       gccGGAGTCCTCAAGGAAGGAGAGGAACAGGTACATCATGCCGATCTGCGGCACCAGTGGCAGCACgg                                       <  1:27260/67‑1 (MQ=255)
                       gccGGAGTCCTCAAGGAAGGAGAGGAACAGGTACATCATGCCGATCTGCGGCACCAGTGGCAGCACgg                                       <  1:2206973/67‑1 (MQ=255)
                       gccGGAGTCCTCAAGGAAGGAGAGGAACAGGTACATCATGCCGATCTGCGGCACCAGTGGCAGCACgg                                       <  1:2049378/67‑1 (MQ=255)
                       gccGGAGTCCTCAAGGAAGGAGAGGAACAGGTACATCATGCCGATCTGCGGCACCAGTGGCAGCACgg                                       <  1:1150759/67‑1 (MQ=255)
                       gccGGAGTCCTCAAGGAAGGAGAGGAACAGGTACATCATGCCGATCTGCGGCACCAGTGGCAGCACgg                                       <  1:1979968/67‑1 (MQ=255)
                       gccGGAGTCCTCAAGGAAGGAGAGGAACAGGTACATCATGCCGATCTGCGGCACCAGTGGCAGCACgg                                       <  1:2031972/67‑1 (MQ=255)
                       cccGGAGTCCTCAAGGAAGGAGAGGAACAGGTACATCATGCCAATCTGCGGCACCAGTGGCAGCACGgt                                      <  1:862521/69‑1 (MQ=255)
                       cccGGAGTCCTCAAGGAAGGAGAGGAACAGGTACATCATGCCAATCTGCGGCACCAGTGGCAGCACGgt                                      <  1:2933759/69‑1 (MQ=255)
                       cccGGAGTCCTCAAGGAAGGAGAGGAACAGGTACATCATGCCAATCTGCGGCACCAGTGGCAGCACGgt                                      <  1:2294592/69‑1 (MQ=255)
                       cccGGAGTCCTCAAGGAAGGAGAGGAACAGGTACATCATGCCAATCTGCGGCACCAGTGGCAGCACGgt                                      <  1:31126/69‑1 (MQ=255)
                          ggAGTCCTCAAGGAAGGAGAGTAACAGGTACATCATGCCAATCTGCGGCACCAGTGGCAGCACGGTGtt                                   <  1:1168983/69‑1 (MQ=255)
                          ggAGTCCTCAAGGAAGGAGAGGAACAGGTACATCATGCCGATCTGCGGCACCAGTGGCAGCACgg                                       <  1:497655/65‑1 (MQ=255)
                          ggAGTCCTCAAGGAAGGAGAGGAACAGGTACATCATGCCAATCTGCGGCACCAGTGGCAGCACGGTGtt                                   <  1:379903/69‑1 (MQ=255)
                          ggAGTCCTCAAGGAAGGAGAGGAACAGGTACATCATGCCAATCTGCGGCACCAGTGGCAGCACGGTGtt                                   <  1:122007/69‑1 (MQ=255)
                          ggAGTCCTCAAGGAAGGAGAGGAACAGGTACATCATGCCAATCTGCGGCACCAGTGGCAGCACGGTGtt                                   <  1:1935462/69‑1 (MQ=255)
                          ggAGTCCTCAAGGAAGGAGAGGAACAGGTACATCATGCCAATCTGCGGCACCAGTGGCAGCACGGTGtt                                   <  1:1857811/69‑1 (MQ=255)
                          ggAGTCCTCAAGGAAGGAGAGGAACAGGTACATCATGCCAATCTGCGGCACCAGTGGCAGCACGGTGtt                                   <  1:2883650/69‑1 (MQ=255)
                                      gaaggaGAGGAACAGGTACATCATGCCAATCTGCGGCACCAGTGGCAGCAc                                         >  1:13485/1‑51 (MQ=255)
                                      gaaggaGAGGAACAGGTACATCATGCCAATCTGCGGCACCAGTGGCAGCACGGTGTTAATgccgcc                          >  1:1456250/1‑66 (MQ=255)
                                      gaaggaGAGGAACAGGTACATCATGCCAATCTGCGGCACCAGTGGCAGCACGGTGTTAATGCCGCCa                         >  1:3229529/1‑67 (MQ=255)
                                      gaaggaGAGGAACAGGTACATCATGCCAATCTGCGGCACCAGTGGCAGCACGGTGTTAATGCCGCCa                         >  1:882160/1‑67 (MQ=255)
                                      gaaggaGAGGAACAGGTACATCATGCCAATCTGCGGCACCAGTGGCAGCACGGTGTTAATGCCGCCa                         >  1:2284138/1‑67 (MQ=255)
                                      gaaggaGAGGAACAGGTACATCATGCCAATCTGCGGCACCAGTGGCAGCACGGTGTTAATGCCGCCa                         >  1:1289965/1‑67 (MQ=255)
                                      gaaggaGAGGAACAGGTACATCATGCCAATCTGCGGCACCAGTGGCAGCACGGTGTTAATGCCGCCa                         >  1:587844/1‑67 (MQ=255)
                                      gaaggaGAGGAACAGGTACATCATGCCAATCTGCGGCACCAGTGGCAGCACGGTGTTAATGCCGCCa                         >  1:430117/1‑67 (MQ=255)
                                      gaaggaGAGGAACAGGTACATCATGCCAATCTGCGGCACCAGTGGCAGCACGGTGTTAATGCCGCCa                         >  1:1452639/1‑67 (MQ=255)
                                      gaaggaGAGGAACAGGTACATCATGCCAATCTGCGGCACCAGTGGCAGCACGGTGTTAATGCCGCCa                         >  1:1527547/1‑67 (MQ=255)
                                      gaaggaGAGGAACAGGTACATCATGCCAATCTGCGGCACCAGTGGCAGCACGGTGTTAATGCCGCCa                         >  1:2315749/1‑67 (MQ=255)
                                      gaaggaGAGGAACAGGTACATCATGCCAATCTGCGGCACCAGTGGCAGCACGGTGTTAATGCCGCCa                         >  1:195868/1‑67 (MQ=255)
                                      gaaggaGAGGAACAGGTACATCATGCCAATCTGCGGCACCAGTGGCAGCACGGTGTTAATGCCGCCa                         >  1:2489535/1‑67 (MQ=255)
                                      gaaggaGAGGAACAGGTACATCATGCCAATCTGCGGCACCAGTGGCAGCACGGTGTTAATGCCGCCa                         >  1:2397669/1‑67 (MQ=255)
                                      gaaggaGAGGAACAGGTACATCATGCCAATCTGCGGCACCAGTGGCAGCACGGTGTTAATGCCGCCa                         >  1:2345224/1‑67 (MQ=255)
                                            gagGAACAGGTACATCATGCCAATCTGCGGCACCAGTGGCAGCACGGTGTTAATGCCGCCACCCAGGcc                 >  1:163913/1‑69 (MQ=255)
                                               gAACAGGTACATCATGCCAATCTGCGGCACCAGTGGCAGCACGGTGTTAATGCCGCCACCCAGGCCCTg              >  1:1975433/1‑69 (MQ=255)
                                                           catGCCAATCTGCGGCACCAGTGGCAGCACGGTGTTAATGCCGCCACCCAGGCCCTGGGCGAGGAAGAt  >  1:680014/1‑69 (MQ=255)
                                                           catGCCAATCTGCGGCACCAGTGGCAGCACGGTGTTAATGCCGCCACCCAGGCCCTGGGCGAGGAAGAt  >  1:937005/1‑69 (MQ=255)
                                                                 |                                                              
ACAAACGCCGCACGCGCCATATACCCGGAGTCCTCAAGGAAGGAGAGGAACAGGTACATCATGCCAATCTGCGGCACCAGTGGCAGCACGGTGTTAATGCCGCCACCCAGGCCCTGGGCGAGGAAGAT  >  minE/2679554‑2679681

Alignment Legend
Aligned base mismatch/match (shaded by quality score): ATCG/ATCG < 3 ≤ ATCG/ATCG < 21 ≤ ATCG/ATCG < 32 ≤ ATCG/ATCG < 36 ≤ ATCG/ATCG
Unaligned base: atcg    Masked matching base: atcg    Alignment gap:     Deleted base: