Predicted mutation
evidence seq id position mutation freq annotation gene description
RA minE 2,182,550 T→G 56.3% A942A (GCT→GCG gltB → glutamate synthase, large subunit

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*minE2,182,5500TG56.3% 13.1 / 26.2 32A942A (GCT→GCGgltBglutamate synthase, large subunit
Reads supporting (aligned to +/- strand):  ref base T (4/10);  new base G (16/2);  total (20/12)
Fisher's exact test for biased strand distribution p-value = 7.90e-04
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

GTGAAGGCGGCGAAGACCCGGCGCGCTACGGCACCAACAAAGTGTCGCGCATCAAGCAGGTGGCTTCCGGTCGCTTTGGGGTTACTCCGGCGTATCTGGTCAATGCCGACGTCATTCAGATTAA  >  minE/2182486‑2182609
                                                                |                                                           
gcgAAGGCGGTGAAGACCCGGCGCGCTATGGCACTAACAAAGTGTCGCGCATCAAGCAGGTGGCGTc                                                           <  1:1109174/65‑1 (MQ=255)
   aaGGCGGTGAAGACCCGGCGCGCTATGGCACTAACACAGTGTCGCGCATCAAGCAGGTGGCGTc                                                           <  1:1653985/64‑1 (MQ=255)
           gAAGACCCGGCGCGCTATGGCACTAACAAAGTGTCGCGCATGAAGCAGGTGGCGTCCGGTCGCTTCggc                                              >  1:1555900/1‑68 (MQ=255)
           gAAGACCCGGCGCGCTATGGCACTAACAAAGTGTCGCGCATCAAGCAGGTGGCGTCCGGTCGCTTCggc                                              >  1:962314/1‑68 (MQ=255)
           gAAGACCCGGCGCGCTATGGCACTAACAAAGTGTCGCGCATCAAGCAGGTGGCGTCCGGTCGCTTCggc                                              >  1:79530/1‑68 (MQ=255)
           gAAGACCCGGCGCGCTATGGCACTAACAAAGTGTCGCGCATCAAGCAGGTGGCGTCCGGTCGCTTCggc                                              >  1:672859/1‑68 (MQ=255)
           gAAGACCCGGCGCGCTATGGCACTAACAAAGTGTCGCGCATCAAGCAGGTGGCGTCCGGTCGCTTCggc                                              >  1:544200/1‑68 (MQ=255)
           gAAGACCCGGCGCGCTATGGCACTAACAAAGTGTCGCGCATCAAGCAGGTGGCGTCCGGTCGCTTCggc                                              >  1:423835/1‑68 (MQ=255)
           gAAGACCCGGCGCGCTATGGCACTAACAAAGTGTCGCGCATCAAGCAGGTGGCGTCCGGTCGCTTCggc                                              >  1:400753/1‑68 (MQ=255)
           gAAGACCCGGCGCGCTATGGCACTAACAAAGTGTCGCGCATCAAGCAGGTGGCGTCCGGTCGCTTCggc                                              >  1:3409152/1‑68 (MQ=255)
           gAAGACCCGGCGCGCTATGGCACTAACAAAGTGTCGCGCATCAAGCAGGTGGCGTCCGGTCGCTTCggc                                              >  1:3127595/1‑68 (MQ=255)
           gAAGACCCGGCGCGCTATGGCACTAACAAAGTGTCGCGCATCAAGCAGGTGGCGTCCGGTCGCTTCggc                                              >  1:1128306/1‑68 (MQ=255)
           gAAGACCCGGCGCGCTATGGCACTAACAAAGTGTCGCGCATCAAGCAGGTGGCGTCCGGTCGCTTCggc                                              >  1:1260217/1‑68 (MQ=255)
           gAAGACCCGGCGCGCTATGGCACTAACAAAGTGTCGCGCATCAAGCAGGTGGCGTCCGGTCGCTTCggc                                              >  1:1568388/1‑68 (MQ=255)
           gAAGACCCGGCGCGCTATGGCACTAACAAAGTGTCGCGCATCAAGCAGGTGGCGTCCGGTCGCTTCggc                                              >  1:1618116/1‑68 (MQ=255)
           gAAGACCCAGCGCGCTATGGCACTAACAAAGTGTCGCGCATCAAGCAGGTGGCGTCCGGTCGCTTCggc                                              >  1:3171711/1‑68 (MQ=255)
           gAAGACCAGGCGCGCTATGGCACTAACAAAGTGTCGCGCATCAAGCAGGTGGCGTCCGGTCGCTTCggc                                              >  1:2184351/1‑68 (MQ=255)
           gAAGACACGGCGCGCTATGGCACTAACAAAGTGTCGCGCATCAAGCAGGTGGCGTCCGGTCGCTTCggc                                              >  1:6142/1‑68 (MQ=255)
                cccGGCGCGCTACGGCACCAACAAAGTGTCGCGCATCAAGCAGGTGGCTTCCGGTCGCTTTGGGGTTAc                                         <  1:1300974/69‑1 (MQ=255)
                 ccGGCGCGCTACGGCACCAACAAAGTGTCGCGCATCAAGCAGGTGGCTTCCGGTCGCTTTGGGGTTAc                                         <  1:1671492/68‑1 (MQ=255)
                        gcTACGGCACCAACAAAGTGTCGCGCATCAAGCAGGTGGCTTCCGGTCGCTTTGGGGTTACTCCGGCGt                                 <  1:1055417/69‑1 (MQ=255)
                                aCCAACAAAGTGTCGCGCATCAAGCAGGTGGCTTCCGGTCGCTTTGGGGTTACTCCGGCGTATCTGGTc                         <  1:3207788/69‑1 (MQ=255)
                                              gcgcATCAAGCAGGTGGCTTCCGGTCGCTTTGGGGTTACTCCGGCGTa                                >  1:3420742/1‑48 (MQ=255)
                                              gcgcATCAAGCAGGTGGCTTCCGGTCGCTTTGGGGTTACTCCGGCGTa                                >  1:2604035/1‑48 (MQ=255)
                                                        cAGGTGGCTTCCGGTCGCTTTGGGGTTACTCCGGCGTATCTGGTCAATGCCGACGTCATTCAGATTaa  <  1:803340/68‑1 (MQ=255)
                                                             ggCTTCCGGTCGCTTTGGGGTTACTCCGGCGTATCTGGTCAATGCCGACGTCATTCa        >  1:3221418/1‑57 (MQ=255)
                                                             ggCTTCCGGTCGCTTTGGGGTTACTCCGGCGTATCTGGTCAATGCCGACGTCATTCa        <  1:1370530/57‑1 (MQ=255)
                                                             ggCTTCCGGTCGCTTTGGGGTTACTCCGGCGTATCTGGTCAATGCCGACGTCATTCa        >  1:24952/1‑57 (MQ=255)
                                                             ggCTTCCGGTCGCTTTGGGGTTACTCCGGCGTATCTGGTCAATGCCGACGTCATTCa        <  1:630221/57‑1 (MQ=255)
                                                             ggCTTCCGGTCGCTTTGGGGTTACTCCGGCGTATCTGGTCAATGCCGACGTCATTCa        <  1:2062014/57‑1 (MQ=255)
                                                             ggCTTCCGGTCGCTTTGGGGTTACTCCGGCGTATCTGGTCAATGCCGACGTCATTCa        <  1:1898538/57‑1 (MQ=255)
                                                             ggCTTCCGGTCGCTTTGGGGTTACTCCGGCGTATCTGGTCAATGCAGACGTCATTCa        <  1:2464472/57‑1 (MQ=255)
                                                                |                                                           
GTGAAGGCGGCGAAGACCCGGCGCGCTACGGCACCAACAAAGTGTCGCGCATCAAGCAGGTGGCTTCCGGTCGCTTTGGGGTTACTCCGGCGTATCTGGTCAATGCCGACGTCATTCAGATTAA  >  minE/2182486‑2182609

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: