Predicted mutation
evidence seq id position mutation freq annotation gene description
RA minE 1,525,931 C→T 84.8% E309E (GAG→GAA fabB ← 3‑oxoacyl‑[acyl‑carrier‑protein] synthase I

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*minE1,525,9310CT84.8% 61.0 / 3.2 39E309E (GAG→GAAfabB3‑oxoacyl‑[acyl‑carrier‑protein] synthase I
Reads supporting (aligned to +/- strand):  ref base C (4/2);  new base T (13/20);  total (17/22)
Fisher's exact test for biased strand distribution p-value = 3.74e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 7.92e-01

CGGTCATGGCTTTGGTTGCAGAAATCGCCGGGCTCTTATCGCCGAACACTTCACGGATAGCTGCCAGCTCTTTCACGTCGCCAACCGGAGTCGAAGTACCGTGGGAGTTCAGGTAATCGATTGGGGTATCAAC  >  minE/1525864‑1525996
                                                                   |                                                                 
cGGTCATGGCTTTGGTTGCAGAAATCGCCGGGCTCTTATCGCCGAACACTTCACGGATAGCTGCCAGct                                                                  >  1:1436401/1‑69 (MQ=255)
cGGTCATGGCTTTGGTTGCAGAAATCGCCGGGCTCTTATCGCCGAACACTTCACGGATAGCTGCCAGct                                                                  >  1:3350954/1‑69 (MQ=255)
       ggCTTTGGTTGCAGAAATCGCCGGGCTCTTATCGCCGAACACTTCACGGATAGCTGCCAGCTCTTTCAc                                                           >  1:934928/1‑69 (MQ=255)
                             gggCTCTTATCGCCGAACACTTCACGGATAGCTGCCAGCTCTTTCACGTCGCCAACCGGAGTCGAAGTa                                     <  1:1212570/69‑1 (MQ=255)
                                ctctTATCGCCGAACACTTCACGGATAGCTGCCAGTTCTTTCa                                                            >  1:2904024/1‑43 (MQ=255)
                                ctctTATCGCCGAACACTTCACGGATAGCTGCCAGTTCTTTCACGTCGCCAACCGGAGTCGAAGTAcc                                   >  1:1034582/1‑68 (MQ=255)
                                ctctTATCGCCGAACACTTCACGGATAGCTGCCAGTTCTTTCACGTCGCCAACCGGAGTCGAAGTAcc                                   >  1:2768859/1‑68 (MQ=255)
                                ctctTATCGCCGAACACTTCACGGATAGCTGCCAGTTCTTTCACGTCGCCAACCGGAGTCGAAGTAcc                                   >  1:231724/1‑68 (MQ=255)
                                ctctTATCGCCGAACACTTCACGGATAGCTGCCAGTTCTTTCACGTCGCCAACCGGAGTCGAAGTAcc                                   >  1:2174945/1‑68 (MQ=255)
                                ctctTATCGCCGAACACTTCACGGATAGCTGCCAGTTCTTTCACGTCGCCAACCGGAGTCGAAGTAcc                                   >  1:1978799/1‑68 (MQ=255)
                                ctctTATCGCCGAACACTTCACGGATAGCTGCCAGTTCTTTCACGTCGCCAACCGGAGTCGAAGTAcc                                   >  1:188908/1‑68 (MQ=255)
                                ctctTATCGCCGAACACTTCACGGATAGCTGCCAGTTCTTTCACGTCGCCAACCGGAGTCGAAGTAcc                                   >  1:298593/1‑68 (MQ=255)
                                ctctTATCGCCGAACACTTCACGGATAGCTGCCAGTTCTTTCACGTCGCCAACCGGAGTCGAAGTAcc                                   >  1:3290757/1‑68 (MQ=255)
                                ctctTATCGCCGAACACTTCACGGATAGCTGCCAGTTCTTTCACGTCGCCAACCGGAGTCGAAGTAcc                                   >  1:11365/1‑68 (MQ=255)
                                ctctTATCGCCGAACACTTCACGGATAGCTGCCAGTTCTTTCACGTCGCCAACCGGAGTCGAAGTAcc                                   >  1:1087341/1‑68 (MQ=255)
                                ctctTATCGCCGAACACTTCACGGATAGCTGCAAGTTAttt                                                              >  1:3249330/1‑41 (MQ=37)
                                ctctTATCGCCGAACACTTCACGGATAGCAGCCAGTTCTTTCACGTCGCCAACCGGAGTCGAAGTAcc                                   >  1:1974557/1‑68 (MQ=255)
                                        gCCGAACACTTCACGGATAGCTGCCAGCTCTTTCACGTCGCCAACCGGAGTCGAAGTACCGTggg                              >  1:98699/1‑65 (MQ=255)
                                              cacTTCACGGATAGCTGCCAGCTCTTTCACGTCGCCAACCGGAGTCGAAGTACCGTGGGAGTTCAGGTa                    <  1:661929/69‑1 (MQ=255)
                                               acTTCACGGATAGCTGCCAGCTCTTTCACGTCGCCAACCGGAGTCGAAGTACCGTGGGAGTTCAGGTaa                   >  1:3093953/1‑69 (MQ=255)
                                                     cGGATAGCTGCCAGTTCTTTCACGTCGCCAAGCGGAGTCGAAGTACGATGAGAGTTCa                        <  1:103747/58‑1 (MQ=255)
                                                     cGGATAGCTGCCAGTTCTTTCACGTCGCCAACCGGAGTCGAAGTACCATGAGAGTTCa                        <  1:338688/58‑1 (MQ=255)
                                                     cGGATAGCTGCCAGTTCTTTCACGTCGCCAACCGGAGTCGAAGTACCATGAGAGTTCa                        <  1:1121523/58‑1 (MQ=255)
                                                     cGGATAGCTGCCAGTTCTTTCACGTCGCCAACCGGAGTCGAAGTACCATGAGAGTTCa                        <  1:899351/58‑1 (MQ=255)
                                                     cGGATAGCTGCCAGTTCTTTCACGTCGCCAACCGGAGTCGAAGTACCATGAGAGTTCa                        <  1:879397/58‑1 (MQ=255)
                                                     cGGATAGCTGCCAGTTCTTTCACGTCGCCAACCGGAGTCGAAGTACCATGAGAGTTCa                        <  1:687456/58‑1 (MQ=255)
                                                     cGGATAGCTGCCAGTTCTTTCACGTCGCCAACCGGAGTCGAAGTACCATGAGAGTTCa                        <  1:621921/58‑1 (MQ=255)
                                                     cGGATAGCTGCCAGTTCTTTCACGTCGCCAACCGGAGTCGAAGTACCATGAGAGTTCa                        <  1:369628/58‑1 (MQ=255)
                                                     cGGATAGCTGCCAGTTCTTTCACGTCGCCAACCGGAGTCGAAGTACCATGAGAGTTCa                        <  1:356029/58‑1 (MQ=255)
                                                     cGGATAGCTGCCAGTTCTTTCACGTCGCCAACCGGAGTCGAAGTACCATGAGAGTTCa                        <  1:2413806/58‑1 (MQ=255)
                                                     cGGATAGCTGCCAGTTCTTTCACGTCGCCAACCGGAGTCGAAGTACCATGAGAGTTCa                        <  1:1449210/58‑1 (MQ=255)
                                                     cGGATAGCTGCCAGTTCTTTCACGTCGCCAACCGGAGTCGAAGTACCATGAGAGTTCa                        <  1:1512845/58‑1 (MQ=255)
                                                     cGGATAGCTGCCAGTTCTTTCACGTCGCCAACCGGAGTCGAAGTACCATGAGAGTTCa                        <  1:1561344/58‑1 (MQ=255)
                                                     cGGATAGCTGCCAGTTCTTTCACGTCGCCAACCGGAGTCGAAGTACCATGAGAGTTCa                        <  1:1580828/58‑1 (MQ=255)
                                                     cGGATAGCTGCCAGTTCTTTCACGTCGCCAACCGGAGTCGAAGTACCATGAGAGTTCa                        <  1:2905933/58‑1 (MQ=255)
                                                     cGGATAGCTGCCAGTTCTTTCACGTCGCCAACCGGAGTCGAAGTACCATGAGAGTTCa                        <  1:1913049/58‑1 (MQ=255)
                                                     cGGATAGCTGCCAGTTCTTTCACGTCGCCAACCGGAGTCGAAGTACCATGAGAGTTCa                        <  1:288190/58‑1 (MQ=255)
                                                     cGGATAGCTGCCAGTTCTTTCACGTCGCCAACCGGAGTCGAAGTACCATGAGAGTTCa                        <  1:2640326/58‑1 (MQ=255)
                                                     cGGATAGCTGCCAGTTCTTTCACGTCGCCAACCGGAGTCGAAGTACAATGAGAGTTCa                        <  1:536045/58‑1 (MQ=255)
                                                     cGGATAGCTCCCAGTTCTTTCACGTCGCCAACCGGAGTCGAAGTACCATGAGAGTTCa                        <  1:629220/58‑1 (MQ=255)
                                                                 aGCTCTTTCACGTCGCCAACCGGAGTCGAAGTACCGTGGGAGTTCAGGTAATCGATTGGGGTATCAAc  >  1:2390754/1‑68 (MQ=255)
                                                                   |                                                                 
CGGTCATGGCTTTGGTTGCAGAAATCGCCGGGCTCTTATCGCCGAACACTTCACGGATAGCTGCCAGCTCTTTCACGTCGCCAACCGGAGTCGAAGTACCGTGGGAGTTCAGGTAATCGATTGGGGTATCAAC  >  minE/1525864‑1525996

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