Predicted mutation
evidence seq id position mutation freq annotation gene description
RA minE 395,051 C→T 20.8% D134D (GAC→GAT ahpF → alkyl hydroperoxide reductase, F52a subunit, FAD/NAD(P)‑binding

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*minE395,0510CT20.8% 117.3 / 19.9 81D134D (GAC→GATahpFalkyl hydroperoxide reductase, F52a subunit, FAD/NAD(P)‑binding
Reads supporting (aligned to +/- strand):  ref base C (24/40);  new base T (3/14);  total (27/54)
Fisher's exact test for biased strand distribution p-value = 1.55e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.68e-01

TTGACGGTGATTTTGAATTCGAAACCTATTACTCGCTCTCTTGCCACAACTGCCCGGACGTGGTGCAGGCGCTGAACCTGATGAGCGTACTGAACCCGCGCATCAAGCACA  >  minE/394993‑395103
                                                          |                                                    
ttGACGGTGATTTTGAATTCGAAACCTATTACTCGCTCTCTTGCCACAACTGCCCGGATgtggtg                                                >  1:2199103/1‑65 (MQ=255)
ttGACGGTGATTTTGAATTCGAAACCTATTACTCGCTCTCTTGCCACAACTGCCCGGATgtggtg                                                >  1:536600/1‑65 (MQ=255)
ttGACGGTGATTTTGAATTCGAAACCTATTACTCGCTCTCTTGCCACAACTGCCCGGATgtggtg                                                >  1:3055539/1‑65 (MQ=255)
  gACGGTGATTTTGAATTCGAAACCTATTACTCGCTCTCTTGCCACAACTGCCCGGACGTGGTGCAGgcg                                          >  1:1361163/1‑69 (MQ=255)
  gACGGTGATTTTGAATTCGAAACCTATTACTCGCTCTCTTGCCACAACTGCCCGGACGTGGTGCAGgcg                                          >  1:2010998/1‑69 (MQ=255)
  gACGGTGATTTTGAATTCGAAACCTATTACTCGCTCTCTTGCCACAACTGCCCGGACGTGGTGCAGgcg                                          >  1:764227/1‑69 (MQ=255)
  gACGGTGATTTTGAATTCGAAACCTATTACTCGCTCTCTTGCCACAACTGCCCGGACGTGGTGCAGgcg                                          >  1:1339214/1‑69 (MQ=255)
  gACGGTGATTTTGAATTCGAAACCTATTACTCGCTCTCTTGCCACAACTGCCCGGACGTGGTGCAGgcg                                          >  1:1337573/1‑69 (MQ=255)
  gACGGTGATTTTGAATTCGAAACCTATTACTCGCTCTCTTGCCACAACTGCCCGGACGTGGTGCAGgcg                                          >  1:111922/1‑69 (MQ=255)
      gTGATTTTGAATTCGAAACCTATTACTCGCTCTCTTGCCACAACTGCCCGGACGTGGTGCAGGCGCTGa                                      <  1:967410/69‑1 (MQ=255)
      gTGATTTTGAATTCGAAACCTATTACTCGCTCTCTTGCCACAACTGCCCGGACGTGGTGCAGGCGCTGa                                      <  1:3337282/69‑1 (MQ=255)
      gTGATTTTGAATTCGAAACCTATTACTCGCTCTCTTGCCACAACTGCCCGGACGTGGTGCAGGCGCTGa                                      <  1:1113413/69‑1 (MQ=255)
      gTGATTTTGAATTCGAAACCTATTACTCGCTCTCTTGCCACAACTGCCCGGACGTGGTGCAGGCGCTGa                                      <  1:2617159/69‑1 (MQ=255)
        gATTTTGAATTCGAAACCTATTACTCGCTCTCTTGCCACAACTGCCCGGACGTGGTGCAGGCGCTGAAc                                    >  1:2899119/1‑69 (MQ=255)
          ttttGAATTCGAAACCTATTACTCGCTCTCTTGCCACAACTGCCCGGATgt                                                    <  1:3221177/51‑1 (MQ=255)
          ttttGAATTCGAAACCTATTACTCGCTCTCTTGCCACAACTGCCCGGATgt                                                    <  1:215493/51‑1 (MQ=255)
          ttttGAATTCGAAACCTATTACTCGCTCTCTTGCCACAACTGCCCGGATgt                                                    <  1:3326473/51‑1 (MQ=255)
          ttttGAATTCGAAACCTATTACTCGCTCTCTTGCCACAACTGCCCGGATgt                                                    <  1:2327190/51‑1 (MQ=255)
          ttttGAATTCGAAACCTATTACTCGCTCTCTTGCCACAACTGCCCGGATgt                                                    <  1:2248295/51‑1 (MQ=255)
          ttttGAATTCGAAACCTATTACTCGCTCTCTTGCCACAACTGCCCGGATgt                                                    <  1:3155220/51‑1 (MQ=255)
          ttttGAATTCGAAACCTATTACTCGCTCTCTTGCCACAACTGCCCGGATgt                                                    <  1:311159/51‑1 (MQ=255)
          ttttGAATTCGAAACCTATTACTCGCTCTCTTGCCACAACTGCCCGGATgt                                                    <  1:461325/51‑1 (MQ=255)
          ttttGAATTCGAAACCTATTACTCGCTCTCTTGCCACAACTGCCCGGATgt                                                    <  1:1081468/51‑1 (MQ=255)
          ttttGAATTCGAAACCTATTACTCGCTCTCTTGCCACAACTGCCCGGATgt                                                    <  1:997816/51‑1 (MQ=255)
          ttttGAATTCGAAACCTATTACTCGCTCTCTTGCCACAACTGCCCGGATgt                                                    <  1:2933330/51‑1 (MQ=255)
          ttttGAATTCGAAACCTATTACTCGCTCTCTTGCCACAACTGCCCGGATgt                                                    <  1:1500825/51‑1 (MQ=255)
          ttttGAATTCGAAACCTATTACTCGCTCTCTTGCCACAACTGCCCGGATgt                                                    <  1:1433445/51‑1 (MQ=255)
           tttGAATTCGAAACCTATTACTCGCTCTCTTGCCACAACTGCCCGGATgt                                                    <  1:558314/50‑1 (MQ=255)
           tttGAATTCGAAACCTATTACTCGCTCTCTTGCCACAACTGCCCGGACGTGGTGCAGGCGCTGAACCt                                  <  1:1004612/68‑1 (MQ=255)
           tttGAATTCGAAACCTATTACTCGCTCTCTTGCCACAACTGCCCGGACGTGGTGCAGGCGCTGAACCt                                  <  1:3127400/68‑1 (MQ=255)
           tttGAATTCGAAACCTATTACTCGCTCTCTTGCCACAACTGCCCGGACGTGGTGCAGGCGCTGAACCt                                  <  1:1001816/68‑1 (MQ=255)
           tttGAATTCGAAACCTATTACTCGCTCTCTTGCCACAACTGCCCGGACGTGGTGCAGGCGCTGAACCt                                  <  1:2756990/68‑1 (MQ=255)
           tttGAATTCGAAACCTATTACTCGCTCTCTTGCCACAACTGCCCGGACGTGGTGCAGGCGCTGAACCt                                  <  1:2977409/68‑1 (MQ=255)
           tttGAATTCGAAACCTATTACTCGCTCTCTTGCCACAACTGCCCGGACGTGGTGCAGGCGCTGAACCt                                  <  1:2276636/68‑1 (MQ=255)
           tttGAATTCGAAACCTATTACTCGCTCTCTTGCCACAACTGCCCGGACGTGGTGCAGGCGCTGAACCt                                  <  1:3018386/68‑1 (MQ=255)
           tttGAATTCGAAACCTATTACTCGCTCTCTTGCCACAACTGCCCGGACGTGGTGCAGGCGCTGAACCt                                  <  1:2245274/68‑1 (MQ=255)
           tttGAATTCGAAACCTATTACTCGCTCTCTTGCCACAACTGCCCGGACGTGGTGCAGGCGCTGAACCt                                  <  1:222592/68‑1 (MQ=255)
                 ttCGAAACCTATTACTCGCTCTCTTGCGACAACTGCCCGGACGTGGTGCAGGCGCTGAACCTGATGAg                            <  1:2535466/68‑1 (MQ=255)
                 ttCGAAACCTATTACTCGCTCTCTTGCCACAACTGCCCGGACGTGGTGCAGGCGCTg                                       <  1:2848960/57‑1 (MQ=255)
                 ttCGAAACCTATTACTCGCTCTCTTGCCACAACTGCCCGGACGTGGTGCAGGCGCTg                                       <  1:2973116/57‑1 (MQ=255)
                 ttCGAAACCTATTACTCGCTCTCTTGCCACAACTGCCCGGACGTGGTGCAGGCGCTg                                       <  1:1178505/57‑1 (MQ=255)
                 ttCGAAACCTATTACTCGCTCTCTTGCCACAACTGCCCGGACGTGGTGCAGGCGCTGAACCTGATGAg                            <  1:2752215/68‑1 (MQ=255)
                  tCGAAACCTATTACTCGCTCTCTTGCCACAACTGCCCGGACGTGGTGCAGGCGCTGAACCTGATGAGCg                          <  1:2573423/69‑1 (MQ=255)
                  tCGAAACCTATTACTCGCTCTCTTGCCACAACTGCCCGGACGTGGTGCAGGCGCTGAACCTGATGAGCg                          <  1:2799371/69‑1 (MQ=255)
                  tCGAAACCTATTACTCGCTCTCTTGCCACAACTGCCCGGACGTGGTGCAGGCGCTGAACCTGATGAGCg                          <  1:2803433/69‑1 (MQ=255)
                         cTATTACTCGCTCTCTTGCCACAACTGCCCGGACgtggt                                                 >  1:403400/1‑39 (MQ=255)
                         cTATTACTCGCTCTCTTGCCACAACTGCCCGGACGTGGTGCAGGCGCTg                                       >  1:2821318/1‑49 (MQ=255)
                         cTATTACTCGCTCTCTTGCCACAACTGCCCGGACGTGGTGCAGGCGCTGAACCTGATGAGCGTACTg                     >  1:2589815/1‑67 (MQ=255)
                         cTATTACTCGCTCTCTTGCCACAACTGCCCGGACGTGGTGCAGGCGCTGAACCTGATGAGCGTACTGaa                   >  1:400489/1‑69 (MQ=255)
                         cTATTACTCGCTCTCTTGCCACAACTGCCCGGACGTGGTGCAGGCGCTGAACCTGATGAGCGTACTGaa                   >  1:67052/1‑69 (MQ=255)
                         cTATTACTCGCTCTCTTGCCACAACTGCCCGGACGTGGTGCAGGCGCTGAACCTGATGAGCGTACTGaa                   >  1:78090/1‑69 (MQ=255)
                         cTATTACTCGCTCTCTTGCCACAACTGCCCGGACGTGGTGCAGGCGCTGAACCTGATGAGCGTACTGaa                   >  1:99838/1‑69 (MQ=255)
                         cTATTACTCGCTCTCTTGCCACAACTGCCCGGACGTGGTGCAGGCGCTGAACCTGATGAGCGTACTGaa                   >  1:2933158/1‑69 (MQ=255)
                         cTATTACTCGCTCTCTTGCCACAACTGCCCGGACGTGGTGCAGGCGCTGAACCTGATGAGCGTACTGa                    >  1:1866649/1‑68 (MQ=255)
                              aCTCGCTCTCTTGCCACAACTGCCCGGACGTGGTGCAGGCGCTGAACCt                                  <  1:2317670/49‑1 (MQ=255)
                              aCTCGCTCTCTTGCCACAACTGCCCGGACGTGGTGCAGGCGCTGAACCt                                  <  1:3116760/49‑1 (MQ=255)
                              aCTCGCTCTCTTGCCACAACTGCCCGGACGTGGTGCAGGCGCTGAACCt                                  <  1:3116809/49‑1 (MQ=255)
                              aCTCGCTCTCTTGCCACAACTGCCCGGACGTGGTGCAGGCGCTGAACCt                                  <  1:2452996/49‑1 (MQ=255)
                                tCGCTCTCTTGCCACAACTGCCCGGACGTGGTGCAGGCGCTGAACCTGATGAGCGTACTGAACCcgcg             <  1:2933676/68‑1 (MQ=255)
                                tCGCTCTCTTGCCACAACTGCCCGGACGTGGTGCAGGCGCTGAACCTGATGAGCGTACTGAACCcgcg             <  1:583159/68‑1 (MQ=255)
                                tCGCTCTCTTGCCACAACTGCCCGGACGTGGTGCAGGCGCTGAACCTGATGAGCGTACTGAACCcgcg             <  1:2186760/68‑1 (MQ=255)
                                tCGCTCTCTTGCCACAACTGCCCGGACGTGGTGCAGGCGCTGAACCTGATGAGCGTACTGAACCcgcg             <  1:114137/68‑1 (MQ=255)
                                tCGCTCTCTTGCCACAACTGCCCGGACGTGGTGCAGGCGCTGAACCTGATGAGCGTACTGAACCcgcg             <  1:2614025/68‑1 (MQ=255)
                                tCGCTCTCTTGCCACAACTGCCCGGACGTGGTGCAGGCGCTGAACCTGATGAGCGTACTGAACCcgcg             <  1:2981697/68‑1 (MQ=255)
                                 cGCTCTCTTGCCACAACTGCCCGGACGTGGTGCAGGCGCTGAACCTGATGAGCg                          <  1:2972117/54‑1 (MQ=255)
                                 cGCTCTCTTGCCACAACTGCCCGGACGTGGTGCAGGCGCTGAACCTGATGAGCg                          <  1:321371/54‑1 (MQ=255)
                                 cGCTCTCTTGCCACAACTGCCCGGACGTGGTGCAGGCGCTGAACCTGATGAGCg                          <  1:3248412/54‑1 (MQ=255)
                                 cGCTCTCTTGCCACAACTGCCCGGACGTGGTGCAGGCGCTGAACCTGATGAGCg                          <  1:2722718/54‑1 (MQ=255)
                                 cGCTCTCTTGCCACAACTGCCCGGACGTGGTGCAGGCGCTGAACCTGATGAGCg                          <  1:1969936/54‑1 (MQ=255)
                                 cGCTCTCTTGCCACAACTGCCCGGACGTGGTGCAGGCGCTGAACCTGATGAGCGTACTGAACCcgcg             <  1:1372556/67‑1 (MQ=255)
                                 cGCTCTCTTGCCACAACTGCCCGGACGTGGTGCAGGCGCTGAACCTGATGAGCGTACTGAACCcgcg             <  1:2629355/67‑1 (MQ=255)
                                          gCCACAACTGCCCGGACGTGGTGCAGGCGCTGAACCTGATGAGCGTACTGAACCCGCGCATc         <  1:234343/62‑1 (MQ=255)
                                          gCCACAACTGCCCGGACGTGGTGCAGGCGCTGAACCTGATGAGCGTACTGAACCCGCGCATc         <  1:1953413/62‑1 (MQ=255)
                                          gCCACAACTGCCCGGACGTGGTGCAGGCGCTGAACCTGATGAGCGTACTGAACCCGCGCATCAAGcaca  >  1:180566/1‑69 (MQ=255)
                                          gCCACAACTGCCCGGACGTGGTGCAGGCGCTGAACCTGATGAGCGTACTGAACCCGCGCATCAAGcaca  >  1:1098600/1‑69 (MQ=255)
                                          gCCACAACTGCCCGGACGTGGTGCAGGCGCTGAACCTGATGAGCGTACTGAACCCGCGCATCAAGcaca  >  1:1932242/1‑69 (MQ=255)
                                          gCCACAACTGCCCGGACGTGGTGCAGGCGCTGAACCTGATGAGCGTACTGAACCCGCGCATCAAGcaca  >  1:2269354/1‑69 (MQ=255)
                                          gCCACAACTGCCCGGACGTGGTGCAGGCGCTGAACCTGATGAGCGTACTGAACCCGCGCATCAAGcaca  >  1:2292119/1‑69 (MQ=255)
                                          gCCACAACTGCCCGGACGTGGTGCAGGCGCTGAACCTGATGAGCGTACTGAACCCGCGCATCAAGcaca  >  1:3323938/1‑69 (MQ=255)
                                          gCCACAACTGCCCGGACGTGGTGCAGGCGCTGAACCTGATGAGCGTACTGAACCCGCGCATCAAGcaca  >  1:3185588/1‑69 (MQ=255)
                                          gCCACAACTGCCCGGACGTGGTGCAGGCGCTGAACCTGATGAGCGTACTGAACCCGCGCATCAAGcac   >  1:2607337/1‑68 (MQ=255)
                                                          |                                                    
TTGACGGTGATTTTGAATTCGAAACCTATTACTCGCTCTCTTGCCACAACTGCCCGGACGTGGTGCAGGCGCTGAACCTGATGAGCGTACTGAACCCGCGCATCAAGCACA  >  minE/394993‑395103

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: