Predicted mutation
evidence seq id position mutation freq annotation gene description
RA minE 499,428 A→T 30.7% D272E (GAT→GAA modF ← fused subunits of molybdate transporter and ATP‑binding components of ABC superfamily

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*minE499,4280AT30.7% 59.9 / 32.0 52D272E (GAT→GAAmodFfused subunits of molybdate transporter and ATP‑binding components of ABC superfamily
Reads supporting (aligned to +/- strand):  ref base A (20/16);  new base T (5/11);  total (25/27)
Fisher's exact test for biased strand distribution p-value = 1.38e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 4.86e-01

TTGGCCCGACAATTTGCCAGTGTTCGCCTGGATTCACCTGCCAGCTAAGGTTATTAAGAATGGGGCGATCGTTATAAGAAACCACGCCATTGTTCAGCACAATGCGCGGTTCGTTGGCGGGTAAGGCGTGACG  >  minE/499361‑499493
                                                                   |                                                                 
ttGGCCCGACAATTTGCCAGTGTTCGCCTGGATTCACCTGCCAGCTAAGGTTATTAAGAATGGGGCGAt                                                                  <  1:3009081/69‑1 (MQ=255)
ttGGCCCGACAATTTGCCAGTGTTCGCCTGGATTCACCTGCCAGCTAAGGTTATTAAGAATGGGGCGAt                                                                  <  1:1224635/69‑1 (MQ=255)
ttGGCCCGACAATTTGCCAGTGTTCGCCTGGATTCACCTGCCAGCTAAGGTTATTAAGAATGGGGCGAt                                                                  <  1:2810520/69‑1 (MQ=255)
ttGGCCCGACAATTTGCCAGTGTTCGCCTGGATTCACCTGCCAGCTAAGGTTATTAAGAATGGGGCGAt                                                                  <  1:2335415/69‑1 (MQ=255)
ttGGCCCGACAATTTGCCAGTGTTCGCCTGGATTCACCTGCCAGCTAAGGTTATTAAGAATGGGGCGAt                                                                  <  1:2326481/69‑1 (MQ=255)
 tGGCCCGACAATTTGCCAGTGTTCGCCTGGATTCACCTGCCAGCTAAGGTTATTAAGAATGGGGCGAt                                                                  <  1:2693659/68‑1 (MQ=255)
    cccGACAATTTGCCAGTGTTCGCCTGGATTCACCTGCCAGCTAAGGTTATTAAGAATGGGGCGATCGt                                                               <  1:2931611/68‑1 (MQ=255)
    cccGACAATTTGCCAGTGTTCGCCTGGATTCACCTGCCAGCTAAGGTTATTAAGAATGGGGCGATCGt                                                               <  1:102963/68‑1 (MQ=255)
    cccGACAATTTGCCAGTGTTCGCCTGGATTCACCTGCCAGCTAAGGTTATTAAGAATGGGGCGATCGt                                                               <  1:2216279/68‑1 (MQ=255)
    cccGACAATTTGCCAGTGTTCGCCTGGATTCACCTGCCAGCTAAGGTTATTAAGAATGGGGCGATCGt                                                               <  1:3342208/68‑1 (MQ=255)
                ccAGTGTTCGCCTGGATTCACCTGCCAGCTAAGGTTATTAAGAATGGGGCGTTCGTTATAAGAAACCa                                                   <  1:3309338/68‑1 (MQ=255)
                ccAGTGTTCGCCTGGATTCACCTGCCAGCTAAGGTTATTAAGAATGGGGCGTTCGTTATAAGAAACCa                                                   <  1:2941201/68‑1 (MQ=255)
                ccAGTGTTCGCCTGGATTCACCTGCCAGCTAAGGTTATTAAGAATGGGGCGTTCGTTATAAGAAACCa                                                   <  1:2865919/68‑1 (MQ=255)
                ccAGTGTTCGCCTGGATTCACCTGCCAGCTAAGGTTATTAAGAATGGGGCGTTCGTTATAAGAAACCa                                                   <  1:2851888/68‑1 (MQ=255)
                ccAGTGTTCGCCTGGATTCACCTGCCAGCTAAGGTTATTAAGAATGGGGCGTTCGTTATAAGAAACCa                                                   <  1:3207691/68‑1 (MQ=255)
                ccAGTGTTCGCCTGGATTCACCTGCCAGCTAAGGTTATTAAGAATGGGGCGTTCGTTATAAGAAACCa                                                   <  1:2580133/68‑1 (MQ=255)
                ccAGTGTTCGCCTGGATTCACCTGCCAGCTAAGGTTATTAAGAATGGGGCGTTCGTTATAAGAAACCa                                                   <  1:577934/68‑1 (MQ=255)
                ccAGTGTTCGCCTGGATTCACCTGCCAGCTAAGGTTATTAAGAATGGGGCGTTCGTTATAAGAAACCa                                                   <  1:1868379/68‑1 (MQ=255)
                ccAGTGTTCGCCTGGATTCACCTGCCAGCTAAGGTTATTAAGAATGGGGCGTTCGTTATAAGAAACCa                                                   <  1:1519809/68‑1 (MQ=255)
                ccAGTGTTCGCCTGGATTCACCTGCCAGCTAAGGTTATTAAGAATGGGGCGTTCGTTATAAGAAACCa                                                   <  1:584477/68‑1 (MQ=255)
                      ttCGCCTGGATTCACCTGCCAGCTAAGGTTATTAAGAATGGGGCGATCGTTATAAGAAACCACGCCAtt                                            >  1:1964662/1‑69 (MQ=255)
                            tGGATTCACCTGCCAGCTAAGGTTATTAAGAATGGGGCGTTCGTTATAAGAAACCa                                                   <  1:1515219/56‑1 (MQ=255)
                                  cACCTGCCAGCTAAGGTTATTAAGAATGGGGCGTTCGTTATAAGAAACCACGCCATTGTTCAGCACaa                                 >  1:943160/1‑68 (MQ=255)
                                  cACCTGCCAGCTAAGGTTATTAAGAATGGGGCGTTCGTTATAAGAAACCACGCCATTGTTCAGCACaa                                 >  1:1139517/1‑68 (MQ=255)
                                  cACCTGCCAGCTAAGGTTATTAAGAATGGGGCGTTCGTTATAAGAAACCACGCCATTGTTCAGCACAAt                                >  1:2536582/1‑69 (MQ=255)
                                  cACCTGCCAGCTAAGGTTATTAAGAATGGGGCGTTCGTTATAAGAAACCACGCCATTGTTCAGCACAAt                                >  1:127817/1‑69 (MQ=255)
                                  cACCTGCCAGCTAAGGTTATTAAGAATGGGGCGTTCGTTATAAGAAACCACGCCATTGTTCAGCACAAt                                >  1:3129111/1‑69 (MQ=255)
                                                    attaAGAATGGGGCGATCGTTATAAGAAACCACGCCATTGTTCAGCACAATGCGCGGTTCGTTGGCgg               >  1:1709456/1‑68 (MQ=255)
                                                    attaAGAATGGGGCGATCGTTATAAGAAACCACGCCATTGTTCAGCACAATGCGCGGTTCGTTGGCgg               >  1:719680/1‑68 (MQ=255)
                                                    attaAGAATGGGGCGATCGTTATAAGAAACCACGCCATTGTTCAGCACAATGCGCGGTTCGTTGGCgg               >  1:18164/1‑68 (MQ=255)
                                                    attaAGAATGGGGCGATCGTTATAAGAAACCACGCCATTGTTCAGCACAATGCGCGGTTCGTTGGCgg               >  1:2267626/1‑68 (MQ=255)
                                                    attaAGAATGGGGCGATCGTTATAAGAAACCACGCCATTGTTCAGCACAATGCGCGGTTCGTTGGCgg               >  1:2842060/1‑68 (MQ=255)
                                                    attaAGAATGGGGCGATCGTTATAAGAAACCACGCCATTGTTCAGCACAATGCGCGGTTCGTTGGCgg               >  1:941139/1‑68 (MQ=255)
                                                    attaAGAATGGGGCGATCGTTATAAGAAACCACGCCATTGTTCAGCACAATGCGCGGTTCGTTGGCgg               >  1:1340481/1‑68 (MQ=255)
                                                    attaAGAATGGGGCGATCGTTATAAGAAACCACGCCATTGTTCAGCACAATGCGCGGTTCGTTGGCg                >  1:2886221/1‑67 (MQ=255)
                                                            tGGGGCGATCGTTATAAGAAACCACGCCATTGTTCAGCACAATGCGCGGTTCGTTGGCGGGTAAGGCGt      >  1:643890/1‑69 (MQ=255)
                                                            tGGGGCGATCGTTATAAGAAACCACGCCATTGTTCAGCACAATGCGCGGTTCGTTGGCGGGTAAGGCGt      >  1:2654820/1‑69 (MQ=255)
                                                            tGGGGCGATCGTTATAAGAAACCACGCCATTGTTCAGCACAATGCGCGGTTCGTTGGCGGGTAAGGCGt      >  1:400710/1‑69 (MQ=255)
                                                            tGGGGCGATCGTTATAAGAAACCACGCCATTGTTCAGCACAATGCGCGGTTCGTTGGCGGGTAAGGCGt      >  1:2191606/1‑69 (MQ=255)
                                                            tGGGGCGATCGTTATAAGAAACCACGCCATTGTTCAGCACAATGCGCGGTTCGTTGGCGGGTAAGGCGt      >  1:2150731/1‑69 (MQ=255)
                                                            tGGGGCGATCGTTATAAGAAACCACGCCATTGTTCAGCACAATGCGCGGTTCGTTGGCGGGTAAGGCGt      >  1:1933295/1‑69 (MQ=255)
                                                            tGGGGCGATCGTTATAAGAAACCACGCCATTGTTCAGCACAATGCGCGGTTCGTTGGCGGGTAAGGCGt      >  1:1929204/1‑69 (MQ=255)
                                                            tGGGGCGATCGTTATAAGAAACCACGCCATTGTTCAGCACAATGCGCGGTTCGTTGGCGGGTAAGGCGt      >  1:1722224/1‑69 (MQ=255)
                                                            tGGGGCGATCGTTATAAGAAACCACGCCATTGTTCAGCACAATGCGCGGTTCGTTGGCGGGTAAGGCGt      >  1:1660309/1‑69 (MQ=255)
                                                            tGGGGCGATCGTTATAAGAAACCACGCCATTGTTCAGCACAATGCGCGGTTCGTTGGCGGGTAAGGCGt      >  1:1388378/1‑69 (MQ=255)
                                                            tGGGGCGATCGTTATAAGAAACCACGACATTGTTCAGCACAATGCGCGGTTCGTTGGCGGGTAAGGCGt      >  1:165248/1‑69 (MQ=255)
                                                                gCGATCGTTATAAGATACCACGCCATTGTTCAGCACAATGCGCGGTTCGTTGGCGGGTAAGGCGTGACg  <  1:1073557/69‑1 (MQ=255)
                                                                gCGATCGTTATAAGAAACCACGCCATTGTTCAGCACAATGCGCGGTTCGTTGGCGGGTAAGGCGTGACg  <  1:2992578/69‑1 (MQ=255)
                                                                gCGATCGTTATAAGAAACCACGCCATTGTTCAGCACAATGCGCGGTTCGTTGGCGGGTAAGGCGTGACg  <  1:3180808/69‑1 (MQ=255)
                                                                gCGATCGTTATAAGAAACCACGCCATTGTTCAGCACAATGCGCGGTTCGTTGGCGGGTAAGGCGTGACg  <  1:2841230/69‑1 (MQ=255)
                                                                gCGATCGTTATAAGAAACCACGCCATTGTTCAGCACAATGCGCGGTTCGTTGGCGGGTAAGGCGTGACg  <  1:3432793/69‑1 (MQ=255)
                                                                gCGATCGTTATAAGAAACCACGCCATTGTTCAGCACAATGCGCGGTTCGTTGGCGGGTAAGGCGTGACg  <  1:675538/69‑1 (MQ=255)
                                                                   |                                                                 
TTGGCCCGACAATTTGCCAGTGTTCGCCTGGATTCACCTGCCAGCTAAGGTTATTAAGAATGGGGCGATCGTTATAAGAAACCACGCCATTGTTCAGCACAATGCGCGGTTCGTTGGCGGGTAAGGCGTGACG  >  minE/499361‑499493

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: