Predicted mutation
evidence seq id position mutation freq annotation gene description
RA minE 372,467 C→A 62.0% S1138S (TCC→TCA entF → enterobactin synthase multienzyme complex component, ATP‑dependent

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*minE372,4670CA62.0% 9.3 / 22.8 29S1138S (TCC→TCAentFenterobactin synthase multienzyme complex component, ATP‑dependent
Reads supporting (aligned to +/- strand):  ref base C (4/7);  new base A (7/11);  total (11/18)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 8.45e-01

TGCGAAGCGCATCTGGCAACGTTACTTGAACAACAACCGCATGGCCCTTATTACCTGCTGGGGTATTCCCTTGGCGGTACGCTGGCGCAGGGTATTGCGGCGCGACTGCGTGCCCGTGGCGAACAGGTGGC  >  minE/372399‑372529
                                                                    |                                                              
tGCGAAGCGCATCTGGCAACGTTACTTGAACAACAACCGCATGGCCCTTATTACCTGCTGGGGTATTcc                                                                <  1:398995/69‑1 (MQ=255)
           tCTGGCAACGTTACTTGAACAACAACCGCATGGCCCTTATTACCTGCTGGGGTATTCCCTTGGCGGTAc                                                     <  1:117501/69‑1 (MQ=255)
           tCTGGCAACGTTACTTGAACAACAACCGCATGGCCCTTATTACCTGCTGGGGTATTCCCTTGGCGGTAc                                                     <  1:3048020/69‑1 (MQ=255)
           tCTGGCAACGTTACTTGAACAACAACCGCATGGCCCTTATTACCTGCTGGGGTATTCCCTTGGCGGTAc                                                     <  1:293269/69‑1 (MQ=255)
             tGGCAACGTTACTTGAACAACAACCGCATGGCCCTTATTACCTGCTGGGGTATTCCCTTGGCGGTACGc                                                   >  1:1460471/1‑69 (MQ=255)
             tGGCAACGTTACTTGAACAACAACCGCATGGCCCTTATTACCTGCTGGGGTATTCCCTTGGCGGTACGc                                                   >  1:111651/1‑69 (MQ=255)
             tGGCAACGTTACTTGAACAACAACCGCATGGCCCTTATTACCTGCTGGGGTATTCCCTTGGCGGTACGc                                                   >  1:751419/1‑69 (MQ=255)
             tGGCAACGTTACTTGAACAACAACCGCATGGCCCTTATTACCTGCTGGGGTATTCCCTTGGCGGTACGc                                                   >  1:2113597/1‑69 (MQ=255)
                       aCTTGAACAACAACCGCATGGCCCTTATTACCTGCTGGGGTATTCCCTTGGCGGTACGCTGGCGCAggg                                         <  1:1400537/69‑1 (MQ=255)
                         ttGAACAACAACCGCATGGCCCTTATTACCTGCTGGGGTATTCCCTTGGCGGTACGCTGGCGCAGGGt                                        <  1:3197962/68‑1 (MQ=255)
                                                   taCCTGCTGGGGTATTCCCTTGGCGGTACGCTGGCGCAGGGTATTGCGGCGCGACTGCGTGCCCGTg               <  1:1652562/67‑1 (MQ=255)
                                                       tgctgGGTTATTCACTTGGCGGTACGCTGGCGCAGGGCATTGCGGCGCGACTGCGTGCCCg                 >  1:1352968/1‑61 (MQ=255)
                                                       tgctgGGTTATTCACTTGGCGGTACGCTGGCGCAGGGCATTGCGGCGCGACTGCGTGCCCGTGGCGAAc         >  1:524889/1‑69 (MQ=39)
                                                       tgctgGGTTATTCACTTGGCGGTACGCTGGCGCAGGGCATTGCGGCGCGACTGCGTGCCCGTGGCGAAc         >  1:1826873/1‑69 (MQ=39)
                                                       tgctgGGTTATTCACTTGGCGGTACGCTGGCGCAGGGCATTGCGGCGCGACTGCGTGCCCGTGGCGAAc         >  1:2272703/1‑69 (MQ=39)
                                                       tgctgGGTTATTCACTTGGCGGTACGCTGGCGCAGGGCATTGCGGCGCGACTGCGTGCCCGTGGCGAAc         >  1:2465860/1‑69 (MQ=39)
                                                       tgctgGGTTATTCACTTGGCGGTACGCTGGCGCAGGGCATTGCGGCGCGACTGCGTGCCCGTGGCGAAc         >  1:2579821/1‑69 (MQ=39)
                                                       tgctgGGTTATTCACTTGGCGGTACGCTGGCGCAGGGCATTGCGGCGCGACTGCGTGCCCGTGGCGAAc         >  1:3429452/1‑69 (MQ=39)
                                                        gctgGGGTATTCCCTTGGCGGTACGCTGGCGCAGGGTATTGCGGCGCGACTGCGTGCCCGTg               <  1:2007268/62‑1 (MQ=255)
                                                              ttATTCACTTGGCGGTACGCTGGCGCAGGGCATTGCGGCGCGACTGCGTGCCCGTGGCGAACAGGTGGc  <  1:583376/68‑1 (MQ=39)
                                                              ttATTCACTTGGCGGTACGCTGGCGCAGGGCATTGCGGCGCGACTGCGTGCCCGTGGCGAACAGGTGGc  <  1:3382497/68‑1 (MQ=39)
                                                              ttATTCACTTGGCGGTACGCTGGCGCAGGGCATTGCGGCGCGACTGCGTGCCCGTGGCGAACAGGTGGc  <  1:92744/68‑1 (MQ=39)
                                                              ttATTCACTTGGCGGTACGCTGGCGCAGGGCATTGCGGCGCGACTGCGTGCCCGTGGCGAACAGGTGGc  <  1:3223445/68‑1 (MQ=39)
                                                              ttATTCACTTGGCGGTACGCTGGCGCAGGGCATTGCGGCGCGACTGCGTGCCCGTGGCGAACAGGTGGc  <  1:3130289/68‑1 (MQ=39)
                                                              ttATTCACTTGGCGGTACGCTGGCGCAGGGCATTGCGGCGCGACTGCGTGCCCGTGGCGAACAGGTGGc  <  1:2970658/68‑1 (MQ=39)
                                                              ttATTCACTTGGCGGTACGCTGGCGCAGGGCATTGCGGCGCGACTGCGTGCCCGTGGCGAACAGGTGGc  <  1:2522928/68‑1 (MQ=39)
                                                              ttATTCACTTGGCGGTACGCTGGCGCAGGGCATTGCGGCGCGACTGCGTGCCCGTGGCGAACAGGTGGc  <  1:1915921/68‑1 (MQ=39)
                                                              ttATTCACTTGGCGGTACGCTGGCGCAGGGCATTGCGGCGCGACTGCGTGCCCGTGGCGAACAGGTGGc  <  1:1812631/68‑1 (MQ=39)
                                                              ttATTCACTTGGCGGTACGCTGGCGCAGGGCATTGCGGCGCGACTGCGTGCCCGTGGCGAACAGGTGGc  <  1:1796389/68‑1 (MQ=39)
                                                              ttATTCACTTGGCGGTACGCTGGCGCAGGGCATTGCGGCGCGACTGCGTGCCCGTGGCGAACAGGTGGc  <  1:1521132/68‑1 (MQ=39)
                                                                    |                                                              
TGCGAAGCGCATCTGGCAACGTTACTTGAACAACAACCGCATGGCCCTTATTACCTGCTGGGGTATTCCCTTGGCGGTACGCTGGCGCAGGGTATTGCGGCGCGACTGCGTGCCCGTGGCGAACAGGTGGC  >  minE/372399‑372529

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: