Predicted mutation
evidence seq id position mutation freq annotation gene description
RA minE 1,513,411 C→T 58.8% P156P (CCG→CCA ubiX ← 3‑octaprenyl‑4‑hydroxybenzoate carboxy‑lyase

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*minE1,513,4110CT58.8% 19.4 / 44.3 63P156P (CCG→CCAubiX3‑octaprenyl‑4‑hydroxybenzoate carboxy‑lyase
Reads supporting (aligned to +/- strand):  ref base C (9/17);  new base T (7/30);  total (16/47)
Fisher's exact test for biased strand distribution p-value = 2.40e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 8.18e-01

GGTTATCGCAAACTGGTCAAGAACACGATTAACCGTCTGATTTATCACATCATCAAGGGATTGCGGGCGATGATAAAACGCCGGAACGGGAGGCATAATCACCGCACCGATTTCTGCCGCCTGAGTCATTA  >  minE/1513348‑1513478
                                                               |                                                                   
ggTTATCGCAAACTGGTCAAGAACACGATTAACCGTCTGATTTATCACATCATCAAGGGATTGCGGGCg                                                                <  1:352816/69‑1 (MQ=255)
ggTTATCGCAAACTGGTCAAGAACACGATTAACCGTCTGATTTATCACATCATCAAGGGATTGCGGGCg                                                                <  1:603385/69‑1 (MQ=255)
ggTTATCGCAAACTGGTCAAGAACACGATTAACCGTCTGATTTATCACATCATCAAGGGATTGCGGGCg                                                                <  1:2618119/69‑1 (MQ=255)
    aTCGCAAACTGGTCAAGAACACGATTAACCGTCTGATTTATCACATCATCAAGGGATTGCGGGCgatg                                                             >  1:2006216/1‑68 (MQ=255)
    aTCGCAAACTGGTCAAGAACACGATTAACCGTCTGATTTATCACATCATCAAGGGATTGCGGGCgatg                                                             >  1:420520/1‑68 (MQ=255)
    aTCGCAAACTGGTCAAGAACACGATTAACCGTCTGATTTATCACATCATCAAGGGATTGCGGGCgatg                                                             >  1:461509/1‑68 (MQ=255)
    aTCGCAAACTGGTCAAGAACACGATTAACCGTCTGATTTATCACATCATCAAGGGATTGCGGGCgatg                                                             >  1:3393792/1‑68 (MQ=255)
                  aaGAACACGATTAACCGTCTGATTTATCGCATCATCAAGGGATTGTGGGCGATGATAAAACGCCGGAAc                                              >  1:952333/1‑69 (MQ=255)
                     aaCACGATTAACCGTCTGATTTATCACATCATCAAGGTATTGCGGGCGATGATAAAACGCCGGAACggg                                           <  1:3108494/69‑1 (MQ=255)
                      acacGATTAACCGTCTGATTTATCACATCATCAAGGGATTGCGGGCGATGATAAAACGCCGGAACGGGa                                          >  1:654965/1‑69 (MQ=255)
                        acGATTAACCGTCTGATTTATCACATCATCAAGGGATTGCGGGCGATGATAAAACGCCGGAACGGGAgg                                        <  1:2763752/69‑1 (MQ=255)
                        acGATTAACCGTCTGATTTATCACATCATCAAGGGATTGCGGGCGATGATAAAACGCCGGAACGGGAgg                                        <  1:1037285/69‑1 (MQ=255)
                        acGATTAACCGTCTGATTTATCACATCATCAAGGGATTGCGGGCGATGATAAAACGCCGGAACGGGAgg                                        <  1:3244141/69‑1 (MQ=255)
                        acGATTAACCGTCTGATTTATCACATCATCAAGGGATTGCGGGCGATGATAAAACGCCGGAACGGGAgg                                        <  1:3206550/69‑1 (MQ=255)
                        acGATTAACCGTCTGATTTATCACATCATCAAGGGATTGCGGGCGATGATAAAACGCCGGAACGGGAgg                                        <  1:2600354/69‑1 (MQ=255)
                        acGATTAACCGTCTGATTTATCACATCATCAAGGGATTGCGGGCGATGATAAAACGCCGGAACGGGAgg                                        <  1:3194309/69‑1 (MQ=255)
                        acGATTAACCGTCTGATTTATCACATCATCAAGGGATTGCGGGCGATGATAAAACGCCGGAACGGGAgg                                        <  1:974188/69‑1 (MQ=255)
                        acGATTAACCGTCTGATTTATCACATCATCAAGGGATTGCGGGCGATGATAAAACGCCGGAACGGGAgg                                        <  1:1902506/69‑1 (MQ=255)
                        acGATTAACCGTCTGATTTATCACATCATCAAGGGATTGCGGGCGATGATAAAACGCCGGAACGGGAgg                                        <  1:1566944/69‑1 (MQ=255)
                        acGATTAACCGTCTGATTTATCACATCATCAAGGGATTGCGGGCGATGATAAAACGCCGGAACGGGAgg                                        <  1:849542/69‑1 (MQ=255)
                        acGATTAACCGTCTGATTTATCACATCATCAAGGGATTGCGGGCGATGATAAAACGCCGGAACGGGAgg                                        <  1:1414893/69‑1 (MQ=255)
                        acGATTAACCGTCTGATTTATCACATCATCAAGGGATTGCGGGCGATGATAAAACCCCGGAACGGGAgg                                        <  1:181720/69‑1 (MQ=255)
                         cGATTAACCGTCTGATTTATCGCATCATCAAGGGATTGTGGGCGATGATAAAACGCCGGAACGGGAgg                                        >  1:1219997/1‑68 (MQ=255)
                         cGATTAACCGTCTGATTTATCGCATCATCAAGGGATTGTGGGCGATGATAAAACGCCGGAACGGGAgg                                        >  1:1735438/1‑68 (MQ=255)
                         cGATTAACCGTCTGATTTATCGCATCATCAAGGGATTGTGGGCGATGATAAAACGCCGGAACGGGAgg                                        >  1:1713739/1‑68 (MQ=255)
                         cGATTAACCGTCTGATTTATCGCATCATCAAGGGATTGTGGGCGATGATAAAACGCCGGAACGGGAgg                                        >  1:1693055/1‑68 (MQ=255)
                         cGATTAACCGTCTGATTTATCGCATCATCAAGGGATTGTGGGCGATGATAAAACGCCGGAACGGGAgg                                        >  1:3420273/1‑68 (MQ=255)
                         cGATTAACCGTCTGATTTATCGCATCATCAAGGGATTGTGGGCGATGATAAAACGCCGGAACGGGAgg                                        >  1:2703281/1‑68 (MQ=255)
                                ccGTCTGATTTATCGCATCATCAAGTGATTGTGGGCGATGATAAAACGCCGGAACGGGAGGCATAATc                                 <  1:30085/68‑1 (MQ=255)
                                ccGTCTGATTTATCGCATCATCAAGGGATTGTGGGCGATGATAAAACGCCGGAACGGGAGGCATAATc                                 <  1:604963/68‑1 (MQ=255)
                                ccGTCTGATTTATCGCATCATCAAGGGATTGTGGGCGATGATAAAACGCCGGAACGGGAGGCATAATc                                 <  1:2770672/68‑1 (MQ=255)
                                ccGTCTGATTTATCGCATCATCAAGGGATTGTGGGCGATGATAAAACGCCGGAACGGGAGGCATAATc                                 <  1:311815/68‑1 (MQ=255)
                                ccGTCTGATTTATCGCATCATCAAGGGATTGTGGGCGATGATAAAACGCCGGAACGGGAGGCATAATc                                 <  1:884421/68‑1 (MQ=255)
                                ccGTCTGATTTATCGCATCATCAAGGGATTGTGGGCGATGATAAAACGCCGGAACGGGAGGCATAATc                                 <  1:3218808/68‑1 (MQ=255)
                                ccGTCTGATTTATCGCATCATCAAGGGATTGTGGGCGATGATAAAACGCCGGAACGGGAGGCATAATc                                 <  1:3417339/68‑1 (MQ=255)
                                ccGTCTGATTTATCGCATCATCAAGGGATTGTGGGCGATGATAAAACGCCGGAACGGGAGGCATAATc                                 <  1:368957/68‑1 (MQ=255)
                                ccGTCTGATTTATCGCATCATCAAGGGATTGTGGGCGATGATAAAACGCCGGAACGGGAGGCATAATc                                 <  1:373736/68‑1 (MQ=255)
                                ccGTCTGATTTATCGCATCATCAAGGGATTGTGGGCGATGATAAAACGCCGGAACGGGAGGCATAATc                                 <  1:437071/68‑1 (MQ=255)
                                ccGTCTGATTTATCGCATCATCAAGGGATTGTGGGCGATGATAAAACGCCGGAACGGGAGGCATAATc                                 <  1:456770/68‑1 (MQ=255)
                                ccGTCTGATTTATCGCATCATCAAGGGATTGTGGGCGATGATAAAACGCCGGAACGGGAGGCATAATc                                 <  1:250623/68‑1 (MQ=255)
                                ccGTCTGATTTATCGCATCATCAAGGGATTGTGGGCGATGATAAAACGCCGGAACGGGAGGCATAATc                                 <  1:1098860/68‑1 (MQ=255)
                                ccGTCTGATTTATCGCATCATCAAGGGATTGTGGGCGATGATAAAACGCCGGAACGGGAGGCATAATc                                 <  1:1183021/68‑1 (MQ=255)
                                ccGTCTGATTTATCGCATCATCAAGGGATTGTGGGCGATGATAAAACGCCGGAACGGGAGGCATAATc                                 <  1:1483392/68‑1 (MQ=255)
                                ccGTCTGATTTATCGCATCATCAAGGGATTGTGGGCGATGATAAAACGCCGGAACGGGAGGCATAATc                                 <  1:1599718/68‑1 (MQ=255)
                                ccGTCTGATTTATCGCATCATCAAGGGATTGTGGGCGATGATAAAACGCCGGAACGGGAGGCATAATc                                 <  1:1610405/68‑1 (MQ=255)
                                ccGTCTGATTTATCGCATCATCAAGGGATTGTGGGCGATGATAAAACGCCGGAACGGGAGGCATAATc                                 <  1:1713406/68‑1 (MQ=255)
                                ccGTCTGATTTATCGCATCATCAAGGGATTGTGGGCGATGATAAAACGCCGGAACGGGAGGCATAATc                                 <  1:2423944/68‑1 (MQ=255)
                                ccGTCTGATTTATCGCATCATCAAGGGATTGTGGGCGATGATAAAACGCCGGAACGGGAGGCATAATc                                 <  1:2473491/68‑1 (MQ=255)
                                ccGTCTGATTTATCGCATCATCAAGGGATTGTGGGCGATGATAAAACGCCGGAACGGGAGGCATAATc                                 <  1:3094170/68‑1 (MQ=255)
                                ccGTCTGATTTATCGCATCATCAAGGGATTGTGGGCGATGATAAAACGCCGGAACGGGAGGCATAATc                                 <  1:2613813/68‑1 (MQ=255)
                                ccGTCTGATTTATCGCATCATCAAGGGATTGTGGGCGATGATAAAACGCCGGAACGGGAGGCATAATc                                 <  1:2617818/68‑1 (MQ=255)
                                ccGTCTGATTTATCGCATCATCAAGGGATTGTGGGCGATGATAAAACGCCGGAACGGGAGGCATAATc                                 <  1:2654700/68‑1 (MQ=255)
                                ccGTCTGATTTATCGCATCATCAAGGGATTGTGGGCGATGATAAAACGCCGGAACGGGAGGCATAATc                                 <  1:2669330/68‑1 (MQ=255)
                                ccGTCTGATTTATCGCATCATCAAGGGATTGTGGGCGATGATAAAACGCCGGAACGGGAGGCATAATc                                 <  1:293230/68‑1 (MQ=255)
                                ccGTCTGATTTATCGCATCATCAAGGGATTGTGGGCGATGATAAAACGCCGGAACGGGAGGCATAATc                                 <  1:3031257/68‑1 (MQ=255)
                                 cGTCTGATTTATCGCATCATCAAGGGATTGTGGGCGATGATAAAACGCCGGAACGGGAGGCATAATc                                 <  1:3263145/67‑1 (MQ=255)
                                  gTCTGATTTATCGCATCATCAAGGGAGCGTGGGCGATGATAAAACGCCGGAACGGGAGGCATAATc                                 <  1:2917120/66‑1 (MQ=255)
                                     tGATTTATCACATCATCAAGGGATTGCGGGCGATGATAAAACGCCGGAACGGGAGGCATAATcaccgca                           >  1:1353730/1‑69 (MQ=255)
                                     tGATTTATCACATCATCAAGGGATTGCGGGCGATGATAAAACGCCGGAACGGGAGGCATAATcaccgc                            >  1:2673106/1‑68 (MQ=255)
                                     tGATTTATCACATCATCAAGGGATTGCGGGCGATGATAAAACGCCGGAACGGGAGGCATAATcaccgc                            >  1:2965833/1‑68 (MQ=255)
                                                   atcaAGGGATTGCGGGCGATGATAAAACGCCGGAACGGGAGGCATAATCACCGCACCGATTTCTgccgc             >  1:2525897/1‑69 (MQ=255)
                                                           aTTGTGGGCGATGATAAAACGCCGGAACGGGAGGCATAATc                                 <  1:65921/41‑1 (MQ=255)
                                                              gCGGGCGATGATAAAACGCCGGAACGGGAGGCATAATCACCGCACCGATTTCTGCCGCCTGAGTCATTa  <  1:1714219/69‑1 (MQ=255)
                                                               |                                                                   
GGTTATCGCAAACTGGTCAAGAACACGATTAACCGTCTGATTTATCACATCATCAAGGGATTGCGGGCGATGATAAAACGCCGGAACGGGAGGCATAATCACCGCACCGATTTCTGCCGCCTGAGTCATTA  >  minE/1513348‑1513478

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: