Predicted mutation
evidence seq id position mutation freq annotation gene description
RA minE 1,141,923 T→C 50.0% I122I (ATT→ATC katE → hydroperoxidase HPII(III)

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*minE1,141,9230TC50.0% ‑4.7 / 45.7 52I122I (ATT→ATCkatEhydroperoxidase HPII(III)
Reads supporting (aligned to +/- strand):  ref base T (13/13);  new base C (14/12);  total (27/25)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

CGCTGCTGGAAGATTTTATTCTGCGCGAGAAAATCACCCACTTTGACCATGAGCGCATTCCGGAACGTATTGTTCATGCACGCGGATCAGCCGCTCACGGTTATTTCCAGCCATATAAAAGCTTA  >  minE/1141865‑1141989
                                                          |                                                                  
cGCTGCTGGAAGATTTTATTCTGCGCGAGAAAATCACCCACTTTGACCATGAGCGCATTCCGGAACGTa                                                          <  1:935614/69‑1 (MQ=255)
cGCTGCTGGAAGATTTTATTCTGCGCGAGAAAATCACCCACTTTGACCATGAGCGCATTCCGGAACGTa                                                          <  1:1668748/69‑1 (MQ=255)
   tgctgGAAGATTTTATTCTGCGCGAGAAAATCACTCACTTTGACCATGAGCGCATCCCGGa                                                               >  1:1332185/1‑61 (MQ=255)
   tgctgGAAGATTTTATTCTGCGCGAGAAAATCACTCACTTTGACCATGAGCGCATCCCGGa                                                               >  1:588384/1‑61 (MQ=255)
   tgctgGAAGATTTTATTCTGCGCGAGAAAATCACTCACTTTGACCATGAGCGCATCCCGGa                                                               >  1:440700/1‑61 (MQ=255)
   tgctgGAAGATTTTATTCTGCGCGAGAAAATCACTCACTTTGACCATGAGCGCATCCCGGa                                                               >  1:3307463/1‑61 (MQ=255)
   tgctgGAAGATTTTATTCTGCGCGAGAAAATCACTCACTTTGACCATGAGCGCATCCCGGa                                                               >  1:166430/1‑61 (MQ=255)
   tgctgGAAGATTTTATTCTGCGCGAGAAAATCACTCACTTTGACCATGAGCGCATCCCGGa                                                               >  1:3032586/1‑61 (MQ=255)
   tgctgGAAGATTTTATTCTGCGCGAGAAAATCACTCACTTTGACCATGAGCGCATCCCGGa                                                               >  1:1793139/1‑61 (MQ=255)
   tgctgGAAGATTTTATTCTGCGCGAGAAAATCACTCACTTTGACCATGAGCGCATCCCGGa                                                               >  1:2432466/1‑61 (MQ=255)
   tgctgGAAGATTTTATTCTGCGCGAGAAAATCACTCACTTTGACCATGAGCGCATCCCGGa                                                               >  1:2412212/1‑61 (MQ=255)
   tgctgGAAGATTTTATTCTGCGCGAGAAAATCACTCACTTTGACCATGAGCGCATCCCGGa                                                               >  1:2334/1‑61 (MQ=255)
   tgctgGAAGATTTTATTCTGCGCGAGAAAATCACTCACTTTGACCATGAGCGCATCCCGGa                                                               >  1:2303340/1‑61 (MQ=255)
   tgctgGAAGATTTTATTCTGCGCGAGAAAATCACTCACTTTGACCATGAGCGCATCCCGGa                                                               >  1:2270510/1‑61 (MQ=255)
   tgctgGAAGATTATATTCTGCGCGAGAAAATCACTCACTTTGACCATGAGCGCATCCCGGa                                                               >  1:2513008/1‑61 (MQ=255)
   tgatgGAAGATTTTATTCTGCGCGAGAAAATCACTCACTTTGACCATGAGCGCATCCCGGa                                                               >  1:3250496/4‑61 (MQ=255)
    gctcGAAGATTTTATTCTGCGCGAGAAAATCACCCACTTTGACCATGAGCGCATTCCGGAACGTATTGt                                                      <  1:3301336/69‑1 (MQ=255)
                 aTTCTGCGCGAGAAAATCACCCACTTTGACCATGAGCGCATTCCGGAACGTATTGTTCATGCACGCgg                                          <  1:3391947/68‑1 (MQ=255)
                 aTTCTGCGCGAGAAAATCACCCACTTTGACCATGAGCGCATTCCGGAACGTATTGTTCATGCACGCgg                                          <  1:3344102/68‑1 (MQ=255)
                 aTTCTGCGCGAGAAAATCACCCACTTTGACCATGAGCGCATTCCGGAACGTATTGTTCATGCACGCgg                                          <  1:2365898/68‑1 (MQ=255)
                 aTTCTGCGCGAGAAAATCACCCACTTTGACCATGAGCGCATTCCGGAACGTATTGTTCATGCACGCgg                                          <  1:1290715/68‑1 (MQ=255)
                    cTGCGCGAGAAAATCACTCACTTTGACCATGAGCGCATCCCGGAACGTATTGTTCATGCACGCGGATc                                       <  1:2607961/68‑1 (MQ=255)
                    cTGCGCGAGAAAATCACTCACTTTGACCATGAGCGCATCCCGGAACGTATTGTTCATGCACGCGGATc                                       <  1:151346/68‑1 (MQ=255)
                    cTGCGCGAGAAAATCACTCACTTTGACCATGAGCGCATCCCGGAACGTATTGTTCATGCACGCGGATc                                       <  1:3131786/68‑1 (MQ=255)
                    cTGCGCGAGAAAATCACTCACTTTGACCATGAGCGCATCCCGGAACGTATTGTTCATGCACGCGGATc                                       <  1:3422162/68‑1 (MQ=255)
                    cTGCGCGAGAAAATCACTCACTTTGACCATGAGCGCATCCCGGAACGTATTGTTCATGCACGCGGATc                                       <  1:1773120/68‑1 (MQ=255)
                    cTGCGCGAGAAAATCACTCACTTTGACCATGAGCGCATCCCGGAACGTATTGTTCATGCACGCGGATc                                       <  1:517458/68‑1 (MQ=255)
                    cTGCGCGAGAAAATCACTCACTTTGACCATGAGCGCATCCCGGAACGTATTGTTCATGCACGCGGATc                                       <  1:1831071/68‑1 (MQ=255)
                    cTGCGCGAGAAAATCACTCACTTTGACCATGAGCGCATCCCGGAACGTATTGTTCATGCACGCGGATc                                       <  1:685620/68‑1 (MQ=255)
                    cTGCGCGAGAAAATCACTCACTTTGACCATGAGCGCATCCCGGAACGTATTGTTCATGCACGCGGATc                                       <  1:2400813/68‑1 (MQ=255)
                    cTGCGCGAGAAAATCACTCACTTTGACCATGAGCGCATCCCGGAACGTATTGTTCATGCACGCGGATc                                       <  1:206078/68‑1 (MQ=255)
                    cTGCGCGAGAAAATCACTCACTTTGACCATGAGCGCATCCCGGAACGTATTGTTCATGCACGCGGATc                                       <  1:2079485/68‑1 (MQ=255)
                    cTGCGCGAGAAAATCACTCACTTTGACCATGAGCGCATCCCGGAACGTATTGTTCATGCACGCGGATc                                       <  1:954670/68‑1 (MQ=255)
                                        cTTTGACCATGAGCGCATTCCGGAACGTATTGTTCATGCACGCgg                                          <  1:1407433/45‑1 (MQ=255)
                                          ttGACCATGAGCGCATTCCGGAACGTATTGTTCATGCACGCgg                                          >  1:29925/1‑43 (MQ=255)
                                          ttGACCATGAGCGCATTCCGGAACGTATTGTTCATGCACGCGGATCAGCCGCTCACGGTTATTTCCAg                 >  1:2810952/1‑68 (MQ=255)
                                          ttGACCATGAGCGCATTCCGGAACGTATTGTTCATGCACGCGGATCAGCCGCTCACGGTTATTTCCAg                 >  1:2814911/1‑68 (MQ=255)
                                          ttGACCATGAGCGCATTCCGGAACGTATTGTTCATGCACGCGGATCAGCCGCTCACGGTTATTTCCAGc                >  1:1458027/1‑69 (MQ=255)
                                          ttGACCATGAGCGCATTCCGGAACGTATTGTTCATGCACGCGGATCAGCCGCTCACGGTTATTTCCAGc                >  1:3078272/1‑69 (MQ=255)
                                          ttGACCATGAGCGCATTCCGGAACGTATTGTTCATGCACGCGGATCAGCCGCTCACGGTTATTTCCAGc                >  1:429719/1‑69 (MQ=255)
                                          ttGACCATGAGCGCATTCCGGAACGTATTGTTCATGCACGCGGATCAGCCGCTCACGGTTATTTCCAGc                >  1:700611/1‑69 (MQ=255)
                                          ttGACCATGAGCGCATTCCGGAACGTATTGTTCATGCACGCGGATCAGCCGCTCACGGTTATTTCCAGc                >  1:3349027/1‑69 (MQ=255)
                                          ttGACCATGAGCGCATTCCGGAACGTATTGTTCATGCACGCGGATCAGCCGCTCACGGTTATTTCCAGc                >  1:3082245/1‑69 (MQ=255)
                                          ttGACCATGAGCGCATTCCGGAACGTATTGTTCATGCACGCGGATCAGCCGCTCACGGTTATTTCCAGc                >  1:21351/1‑69 (MQ=255)
                                          ttGACCATGAGCGCATTCCGGAACGTATTGTTCATGCACGCGGATCAGCCGCTCACGGTTATTTCCAGc                >  1:2110591/1‑69 (MQ=255)
                                          ttGACCATGAGCGCATTCCGGAACGTATTGTTCATGCACGCGGATCAGCCGCTCACGGTTATTTCCAGc                >  1:2158500/1‑69 (MQ=255)
                                                   aGCGCATTCCGGAACGTATTGTTCATGCACGCGGATCAGCCGCTCACGGTTATTTCCAGCCATATaaaa       >  1:1399880/1‑69 (MQ=255)
                                                        aTTCCGGAACGTATTGTTCATGCACGCGGATCAGCCGCTCACGGTTATTTCCAGCCATATAAAAGCTTa  <  1:2577203/69‑1 (MQ=255)
                                                        aTTCCGGAACGTATTGTTCATGCACGCGGATCAGCCGCTCACGGTTATTTCCAGCCATATAAAAGCTTa  <  1:2487018/69‑1 (MQ=255)
                                                        aTTCCGGAACGTATTGTTCATGCACGCGGATCAGCCGCTCACGGTTATTTCCAGCCATATAAAAGCTTa  <  1:2478433/69‑1 (MQ=255)
                                                        aTTCCGGAACGTATTGTTCATGCACGCGGATCAGCCGCTCACGGTTATTTCCAGCCATATAAAAGCTTa  <  1:2023492/69‑1 (MQ=255)
                                                        aTTCCGGAACGTATTGTTCATGCACGCGGATCAGCCGCTCACGGTTATTTCCAGCCATATAAAAGCTTa  <  1:723351/69‑1 (MQ=255)
                                                          |                                                                  
CGCTGCTGGAAGATTTTATTCTGCGCGAGAAAATCACCCACTTTGACCATGAGCGCATTCCGGAACGTATTGTTCATGCACGCGGATCAGCCGCTCACGGTTATTTCCAGCCATATAAAAGCTTA  >  minE/1141865‑1141989

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: