Predicted mutation
evidence seq id position mutation freq annotation gene description
RA minE 419,065:1 +G 100% coding (664/1176 nt) ubiF → 2‑octaprenyl‑3‑methyl‑6‑methoxy‑1,4‑benzoquinol oxygenase

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*minE419,0651.G96.1% 189.6 / 0.2 51coding (664/1176 nt)ubiF2‑octaprenyl‑3‑methyl‑6‑methoxy‑1,4‑benzoquinol oxygenase
Reads supporting (aligned to +/- strand):  ref base . (0/2);  new base G (0/49);  total (0/51)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.86e-01
Rejected as polymorphism: E-value score below prediction cutoff.
Rejected as polymorphism: Frequency below/above cutoff threshold.
Rejected as polymorphism: Variant not supported by required number of reads on each strand.
Rejected as polymorphism: Polymorphic indel expands or contracts a homopolymer stretch.

GCAATTTACTCCGGACGGACCGCGTGCGTTTCTGCCGTTGTTTGATAACTGGGCATCGCTGGTGTGGTATGACTCTCCGGCGCGTATTCGCCAGTTGCAGAATATGAAT  >  minE/419016‑419124
                                                  |                                                           
gcaATTTACTCGGACGGACCGCGTGCGTTTCTGCCGTTGTTTGATAACTGGGCATCGCTGGTGTGGt                                           <  1:1821961/67‑1 (MQ=255)
 caATTTACTCCGGACGGACCGCGTGCGTTTCTGCCGTTGTTTGATAACTGGGCATCGCTGGTGTGGt                                           <  1:442435/67‑1 (MQ=255)
                                            aTAACTGGGGCATCGCTGGTGTGGTATGACTCTCCGGCGCGTATTCGCCAGTTGCAgaatatgaat  <  1:571033/66‑1 (MQ=255)
                                            aTAACTGGGGCATCGCTGGTGTGGTATGACTCTCCGGCGCGTATTCGCCAGTTGCAgaatatgaat  <  1:2076136/66‑1 (MQ=255)
                                            aTAACTGGGGCATCGCTGGTGTGGTATGACTCTCCGGCGCGTATTCGCCAGTTGCAgaatatgaat  <  1:2114142/66‑1 (MQ=255)
                                            aTAACTGGGGCATCGCTGGTGTGGTATGACTCTCCGGCGCGTATTCGCCAGTTGCAgaatatgaat  <  1:2119613/66‑1 (MQ=255)
                                            aTAACTGGGGCATCGCTGGTGTGGTATGACTCTCCGGCGCGTATTCGCCAGTTGCAgaatatgaat  <  1:213076/66‑1 (MQ=255)
                                            aTAACTGGGGCATCGCTGGTGTGGTATGACTCTCCGGCGCGTATTCGCCAGTTGCAgaatatgaat  <  1:236609/66‑1 (MQ=255)
                                            aTAACTGGGGCATCGCTGGTGTGGTATGACTCTCCGGCGCGTATTCGCCAGTTGCAgaatatgaat  <  1:327095/66‑1 (MQ=255)
                                            aTAACTGGGGCATCGCTGGTGTGGTATGACTCTCCGGCGCGTATTCGCCAGTTGCAgaatatgaat  <  1:330537/66‑1 (MQ=255)
                                            aTAACTGGGGCATCGCTGGTGTGGTATGACTCTCCGGCGCGTATTCGCCAGTTGCAgaatatgaat  <  1:371969/66‑1 (MQ=255)
                                            aTAACTGGGGCATCGCTGGTGTGGTATGACTCTCCGGCGCGTATTCGCCAGTTGCAgaatatgaat  <  1:397015/66‑1 (MQ=255)
                                            aTAACTGGGGCATCGCTGGTGTGGTATGACTCTCCGGCGCGTATTCGCCAGTTGCAgaatatgaat  <  1:54496/66‑1 (MQ=255)
                                            aTAACTGGGGCATCGCTGGTGTGGTATGACTCTCCGGCGCGTATTCGCCAGTTGCAgaatatgaat  <  1:2024687/66‑1 (MQ=255)
                                            aTAACTGGGGCATCGCTGGTGTGGTATGACTCTCCGGCGCGTATTCGCCAGTTGCAgaatatgaat  <  1:57191/66‑1 (MQ=255)
                                            aTAACTGGGGCATCGCTGGTGTGGTATGACTCTCCGGCGCGTATTCGCCAGTTGCAgaatatgaat  <  1:680169/66‑1 (MQ=255)
                                            aTAACTGGGGCATCGCTGGTGTGGTATGACTCTCCGGCGCGTATTCGCCAGTTGCAgaatatgaat  <  1:683846/66‑1 (MQ=255)
                                            aTAACTGGGGCATCGCTGGTGTGGTATGACTCTCCGGCGCGTATTCGCCAGTTGCAgaatatgaat  <  1:686793/66‑1 (MQ=255)
                                            aTAACTGGGGCATCGCTGGTGTGGTATGACTCTCCGGCGCGTATTCGCCAGTTGCAgaatatgaat  <  1:715727/66‑1 (MQ=255)
                                            aTAACTGGGGCATCGCTGGTGTGGTATGACTCTCCGGCGCGTATTCGCCAGTTGCAgaatatgaat  <  1:729651/66‑1 (MQ=255)
                                            aTAACTGGGGCATCGCTGGTGTGGTATGACTCTCCGGCGCGTATTCGCCAGTTGCAgaatatgaat  <  1:78831/66‑1 (MQ=255)
                                            aTAACTGGGGCATCGCTGGTGTGGTATGACTCTCCGGCGCGTATTCGCCAGTTGCAgaatatgaat  <  1:835879/66‑1 (MQ=255)
                                            aTAACTGGGGCATCGCTGGTGTGGTATGACTCTCCGGCGCGTATTCGCCAGTTGCAgaatatgaat  <  1:866201/66‑1 (MQ=255)
                                            aTAACTGGGGCATCGCTGGTGTGGTATGACTCTCCGGCGCGTATTCGCCAGTTGCAgaatatgaat  <  1:900620/66‑1 (MQ=255)
                                            aTAACTGGGGCATCGCTGGTGTGGTATGACTCTCCGGCGCGTATTCGCCAGTTGCAgaatatgaat  <  1:938841/66‑1 (MQ=255)
                                            aTAACTGGGGCATCGCTGGTGTGGTATGACTCTCCGGCGCGTATTCGCCAGTTGCAgaatatgaat  <  1:973730/66‑1 (MQ=255)
                                            aTAACTGGGGCATCGCTGGTGTGGTATGACTCTCCGGCGCGTATTCGCCAGTTGCAgaatatgaat  <  1:1682199/66‑1 (MQ=255)
                                            aTAACTGGGGCATCGCTGGTGTGGTATGACTCTCCGGCGCGTATTCGCCAGTTGCAgaatatgaat  <  1:1038040/66‑1 (MQ=255)
                                            aTAACTGGGGCATCGCTGGTGTGGTATGACTCTCCGGCGCGTATTCGCCAGTTGCAgaatatgaat  <  1:1179146/66‑1 (MQ=255)
                                            aTAACTGGGGCATCGCTGGTGTGGTATGACTCTCCGGCGCGTATTCGCCAGTTGCAgaatatgaat  <  1:1195011/66‑1 (MQ=255)
                                            aTAACTGGGGCATCGCTGGTGTGGTATGACTCTCCGGCGCGTATTCGCCAGTTGCAgaatatgaat  <  1:1250447/66‑1 (MQ=255)
                                            aTAACTGGGGCATCGCTGGTGTGGTATGACTCTCCGGCGCGTATTCGCCAGTTGCAgaatatgaat  <  1:1275368/66‑1 (MQ=255)
                                            aTAACTGGGGCATCGCTGGTGTGGTATGACTCTCCGGCGCGTATTCGCCAGTTGCAgaatatgaat  <  1:1279684/66‑1 (MQ=255)
                                            aTAACTGGGGCATCGCTGGTGTGGTATGACTCTCCGGCGCGTATTCGCCAGTTGCAgaatatgaat  <  1:1351919/66‑1 (MQ=255)
                                            aTAACTGGGGCATCGCTGGTGTGGTATGACTCTCCGGCGCGTATTCGCCAGTTGCAgaatatgaat  <  1:1408099/66‑1 (MQ=255)
                                            aTAACTGGGGCATCGCTGGTGTGGTATGACTCTCCGGCGCGTATTCGCCAGTTGCAgaatatgaat  <  1:1530256/66‑1 (MQ=255)
                                            aTAACTGGGGCATCGCTGGTGTGGTATGACTCTCCGGCGCGTATTCGCCAGTTGCAgaatatgaat  <  1:1593361/66‑1 (MQ=255)
                                            aTAACTGGGGCATCGCTGGTGTGGTATGACTCTCCGGCGCGTATTCGCCAGTTGCAgaatatgaat  <  1:1673574/66‑1 (MQ=255)
                                            aTAACTGGGGCATCGCTGGTGTGGTATGACTCTCCGGCGCGTATTCGCCAGTTGCAgaatatgaat  <  1:1022551/66‑1 (MQ=255)
                                            aTAACTGGGGCATCGCTGGTGTGGTATGACTCTCCGGCGCGTATTCGCCAGTTGCAgaatatgaat  <  1:1753360/66‑1 (MQ=255)
                                            aTAACTGGGGCATCGCTGGTGTGGTATGACTCTCCGGCGCGTATTCGCCAGTTGCAgaatatgaat  <  1:1830838/66‑1 (MQ=255)
                                            aTAACTGGGGCATCGCTGGTGTGGTATGACTCTCCGGCGCGTATTCGCCAGTTGCAgaatatgaat  <  1:1837419/66‑1 (MQ=255)
                                            aTAACTGGGGCATCGCTGGTGTGGTATGACTCTCCGGCGCGTATTCGCCAGTTGCAgaatatgaat  <  1:1896339/66‑1 (MQ=255)
                                            aTAACTGGGGCATCGCTGGTGTGGTATGACTCTCCGGCGCGTATTCGCCAGTTGCAgaatatgaat  <  1:1919084/66‑1 (MQ=255)
                                            aTAACTGGGGCATCGCTGGTGTGGTATGACTCTCCGGCGCGTATTCGCCAGTTGCAgaatatgaat  <  1:1922124/66‑1 (MQ=255)
                                            aTAACTGGGGCATCGCTGGTGTGGTATGACTCTCCGGCGCGTATTCGCCAGTTGCAgaatatgaat  <  1:1963916/66‑1 (MQ=255)
                                            aTAACTGGGGCATCGCTGGTGTGGTATGACTCTCCGGCGCGTATTCGCCAGTTGCAgaatatgaat  <  1:1984607/66‑1 (MQ=255)
                                            aTAACTGGGGCATCGCTGGTGTGGTATGACTCTCCGGCGCGTATTCGCCAGTTGCAgaatatgaat  <  1:1997662/66‑1 (MQ=255)
                                             tAACTGGGGCATCGCTGGTGTGGTATGACTCTCCGGCGCGTATTCGCCAGTTGCAgaatatgaat  <  1:412445/65‑1 (MQ=255)
                                             tAACTGGGGCATCGCTGGTGTGGTATGACTCTCCGGCGCGTATTCGCCAGTTGCAgaatatgaat  <  1:1471037/65‑1 (MQ=255)
                                             tAACTGGGGCATCGCTGGTGTGGTATGACTCTCCGGCGCGTATTCGCCAGTTGCAgaatatgaat  <  1:2002336/65‑1 (MQ=255)
                                                  |                                                           
GCAATTTACTCCGGACGGACCGCGTGCGTTTCTGCCGTTGTTTGATAACTGGGCATCGCTGGTGTGGTATGACTCTCCGGCGCGTATTCGCCAGTTGCAGAATATGAAT  >  minE/419016‑419124

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