Predicted mutation
evidence seq id position mutation freq annotation gene description
RA minE 1,142,011 T→C 43.0% F152L (TTC→CTC)  katE → hydroperoxidase HPII(III)

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*minE1,142,0110TC43.0% 30.0 / 105.7 79F152L (TTC→CTC) katEhydroperoxidase HPII(III)
Reads supporting (aligned to +/- strand):  ref base T (25/20);  new base C (24/10);  total (49/30)
Fisher's exact test for biased strand distribution p-value = 2.42e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.81e-01

CACGCGGATCAGCCGCTCACGGTTATTTCCAGCCATATAAAAGCTTAAGCGATATTACCAAAGCGGATTTCCTCTCAGATCCGAACAAAATCACCCCAGTATTTGTACGTTTCTCTACCGTTCAGGGTGGTGCTGGCT  >  minE/1141943‑1142080
                                                                    |                                                                     
cACGCGGATCAGCCGCTCACGGTTATTTCCAGCCATATAAAAGCTTAAGCGATATTACCAAAGCGGAttt                                                                      <  1:898549/70‑1 (MQ=255)
   gcgGATCAGCCGCTCACGGTTATTTCCAGCCATATAAAAGCTTAAGCGATATTACCAAAGCGGATTTCctc                                                                  <  1:277302/71‑1 (MQ=255)
     ggATCAGCCGCTCACGGTTATTTCCAGCCATATAAAAGCTTAAGCGATATTACCAAAGCGGATTTCCtctc                                                                <  1:3559812/71‑1 (MQ=255)
     ggATCAGCCGCTCACGGTTATTTCCAGCCATATAAAAGCTTAAGCGATATTACCAAAGCGGATTTCCtctc                                                                <  1:4056383/71‑1 (MQ=255)
      gATCAGCCGCTCACGGTTATTTCCAGCCATATAAAAGCTTAAGCGATATTACCAAAGCGGATTTCCtctc                                                                <  1:4208560/70‑1 (MQ=255)
       aTCAGCCGCTCACGGTTATTTCCAGCCATATAAAAGCTTAAGCGATATTACCAAAGCGGATTTCctct                                                                 >  1:234748/1‑68 (MQ=255)
        tCAGACGCTCACGGTTATTTCCAGCCATATAAAAGCTTAAGCGATATTACCAAAGCGGATCTCCTCTCAGa                                                             <  1:3146878/71‑1 (MQ=255)
           gCCGCTCACGGTTATTTCCAGCCATATAAAAGCTTAAGCGATATTACCAAAGCGGATCTCCTCTCAGATcc                                                          >  1:10923/1‑71 (MQ=255)
                    ggTTATTTCCAGCCATATAAAAGCTTAAGCGATATTACGAAAGCGGATCTCCTCTCAGATCCGAACAAAAt                                                 <  1:3397597/71‑1 (MQ=255)
                    ggTTATTTCCAGCCATATAAAAGCTTAAGCGATATTACCAAAGCGGATCTCCTCTCAGATCCGAACAAAAt                                                 <  1:2424611/71‑1 (MQ=255)
                    ggTTATTTCCAGCCATATAAAAGCTTAAGCGATATTACCAAAGCGGATCTCCTCTCAGATCCGAACAAAAt                                                 <  1:2893070/71‑1 (MQ=255)
                    ggTTATTTCCAGCCATATAAAAGCTTAAGCGATATTACCAAAGCGGATCTCCTCTCAGATCCGAACAAAAt                                                 <  1:1940564/71‑1 (MQ=255)
                    ggTTATTTCCAGCCATATAAAAGCTTAAGCGATATTACCAAAGCGGATCTCCTCTCAGATCCGAACAAAAt                                                 <  1:3506550/71‑1 (MQ=255)
                     gTTATTTCCAGCCATATAAAAGCTTAAGCGATATTACCAAAGCGGATTTCctct                                                                 <  1:1453334/54‑1 (MQ=255)
                     gTTATTTCCAGCCATATAAAAGCTTAAGCGATATTACCAAAGCGGATCTCCTCTCAGATCCGAACAAAAt                                                 >  1:481914/1‑70 (MQ=255)
                      ttATTTCCAGCCATATAAAAGCTTAAGCGATATTAACAAAGCGGATTTCCTCTCAGATCCGAACAAAATCa                                               >  1:559062/1‑71 (MQ=255)
                        aTTTCCAGCCATATAAAAGCTTAAGCGATATTACCAAAGCGGATCTCCTCTCAGATCCGAACAAAATCAc                                              <  1:3404247/70‑1 (MQ=255)
                           tCCAGCCATATAAAAGCTTAAGCGATATTACCAAAGCGGATTTCCTCTCAGATCCGAACAAAATCAccc                                            >  1:3416768/1‑69 (MQ=255)
                           tCCAGCCATATAAAAGCTTAAGCGATATTACCAAAGCGGATTTCCTCTCAGATCCGAACAAAATCAccc                                            >  1:3680554/1‑69 (MQ=255)
                           tCCAGCCATATAAAAGCTTAAGCGATATTACCAAAGCGGATCTCCTCTCAGATCCGAACAAAATCAccc                                            >  1:366770/1‑69 (MQ=255)
                             cAGCCATATAAAAGCTTAAGCGATATTACCAAAGCGGATTTCCTCTCAGATCCGAACAAAATCACCCCAg                                         <  1:828631/70‑1 (MQ=255)
                             cAGCCATATAAAAGCTTAAGCGATATTACCAAAGCGGATTTCCTCTCAGATCCGAACAAAATCACCCCAGt                                        <  1:1033251/71‑1 (MQ=255)
                                ccATATAAAAGCTTAAGCGATATTACCAAAGCGGATCTCCTCTCAGATCCGAACAATATCACCCCAGTAtt                                     <  1:2917237/71‑1 (MQ=255)
                                 cATATAAAAGCTTAAGCGATATTACCAAAGCGGATTTCCTCTCAGATCCGAACAAAATCACCCCAGTAttt                                    >  1:2278063/1‑71 (MQ=255)
                                 cATATAAAAGCTTAAGCGATATTACCAAAGCGGATTTCCTCTCAGATCCGAACAAAATCACCCCAGTAttt                                    >  1:4188012/1‑71 (MQ=255)
                                   tataAAAGCTTAAGCGATATTACCAAAGCGGATTTCCTCTCAGATCCGAACAAAATCACCCCAGTATTTg                                   >  1:2172604/1‑70 (MQ=255)
                                   tataAAAGCTTAAGCGATATTACCAAAGCGGATTTCCTCTCAGATCCGAACAAAATCACCCCAGTATTTg                                   >  1:2934079/1‑70 (MQ=255)
                                   tataAAAGCTTAAGCGATATTACCAAAGCGGATCTCCTCTCAGATCCGAACaa                                                    <  1:658434/53‑1 (MQ=255)
                                   tataAAAGCTTAAGCGATATTACCAAAGCGGATCTCCTCTCAGATCCGAACaa                                                    <  1:923981/53‑1 (MQ=255)
                                    ataAAAGCTTAAGCGATATTACCAAAGCGGATTTCCTCTCAGATCCGAACAAAATCa                                               >  1:3748467/1‑57 (MQ=255)
                                    ataAAAGCTTAAGCGATATTACCAAAGCGGATTTCCTCTCAGATCCGAACAAAATCa                                               >  1:3954321/1‑57 (MQ=255)
                                    ataAAAGCTTAAGCGATATTACCAAAGCGGATTTCCTCTCAGATCCGAACAAAATCa                                               >  1:2580506/1‑57 (MQ=255)
                                    ataAAAGCTTAAGCGATATTACCAAAGCGGATTTCCTCTCAGATCCGAACAAAATCa                                               >  1:4155553/1‑57 (MQ=255)
                                     taAAAGCTTAAGCGATATTACCAAAGCGGATCTCCTCTCAGATCCGAACAAAAt                                                 >  1:4051210/1‑54 (MQ=255)
                                     taAAAGCTTAAGCGATATTACCAAAGCGGATCTCCTCTCAGATCCGAACAAAAt                                                 >  1:1135864/1‑54 (MQ=255)
                                     taAAAGCTTAAGCGATATTACCAAAGCGGATCTCCTCTCAGATCCGAACAAAAt                                                 >  1:2972880/1‑54 (MQ=255)
                                     taAAAGCTTAAGCGATATTACCAAAGCGGATCTCCTCTCAGATCCGAACAAAATCACCCCAGTATTTGTAc                                >  1:1899065/1‑71 (MQ=255)
                                     taAAAGCTTAAGCGATATTACCAAAGCGGATCTCCTCTCAGATCCGAACAAAATCACCCCAGTATTTGTAc                                >  1:1739771/1‑71 (MQ=255)
                                     taAAAGCTTAAGCGATATTACAAAAGCGGATCTCCTCTCAGATCCGAACAAAATCACCCCAGTATTTGTAc                                >  1:952357/1‑71 (MQ=255)
                                       aaaGCTTAAGCGATATTACCAAAGCGGATTTCCTCTCa                                                               >  1:1505057/1‑38 (MQ=255)
                                       aaaGCTTAAGCGATATTACCAAAGCGGATTTCCTCTCAGATCCGAACAAAATCACCCCAGTATTTGTAc                                >  1:1654137/1‑69 (MQ=255)
                                         aGCTTAAGCGATATTACCAAAGCGGATTTCCTCTCAGATCCGAACAAAATCAcccc                                           >  1:834140/1‑56 (MQ=255)
                                         aGCTTAAGCGATATTACCAAAGCGGATTTCCTCTCAGATCCGAACAAAATCACCCCa                                          >  1:2013382/1‑57 (MQ=255)
                                         aGCTTAAGCGATATTACCAAAGCGGATTTCCTCTCAGATCCGAACAAAATCACCCCa                                          >  1:2915751/1‑57 (MQ=255)
                                         aGCTTAAGCGATATTACCAAAGCGGATTTCCTCTCAGATCCGAACAAAATCACCCCa                                          >  1:2111813/1‑57 (MQ=255)
                                         aGCTTAAGCGATATTACCAAAGCGGATTTCCTCTCAGATCCGAACAAAATCACCCCa                                          >  1:4228868/1‑57 (MQ=255)
                                               aGCGATATTACCAAAGCGGATTTCCTCTCAGATCCGAACAAAATCACCCCAGTATTTGTACGTTtc                           <  1:2574749/66‑1 (MQ=255)
                                               aGCGATATTACCAAAGCGGATTTCCTCTCAGATCCGAACAAAATCACCCCAGTATTTGTACGTTtc                           <  1:1921484/66‑1 (MQ=255)
                                               aGCGATATTACCAAAGCGGATTTCCTCTCAGATCCGAACAAAATCACCCCAGTATTTGTACGTTtc                           <  1:521003/66‑1 (MQ=255)
                                                 cGATATTACCAAAGCGGATTTCCTCTCAGATCCGAACAAAATCACCCCAGTATTTGTACGTTtctc                         >  1:2959739/1‑66 (MQ=255)
                                                   atatTACCAAAGCGGATTTCCTCTCAGATCCGAACaa                                                    <  1:3600580/37‑1 (MQ=255)
                                                   atatTACCAAAGCGGATTTCCTCTCAGATCCGAACaa                                                    <  1:4070408/37‑1 (MQ=255)
                                                   atatTACCAAAGCGGATTTCCTCTCAGATCCGAACAAAATCACCCCAGTATTTGTACGTTTCTCTACCGtt                  >  1:2701362/1‑71 (MQ=255)
                                                   atatTACCAAAGCGGATTTCCTCTCAGATCCGAACAAAATCACCCCAGTATTTGTACGTTTCTCTACCGtt                  >  1:4091019/1‑71 (MQ=255)
                                                    tatTACCAAAGCGGATTTCCTCTCAGATCCGAACAAAATCACCCCAGTATTTGTACGTTTCTCTACCGTTc                 <  1:2650730/71‑1 (MQ=255)
                                                    tatTACCAAAGCGGATCTCCTCTCAGATCCGAACaaaa                                                  >  1:4256960/1‑38 (MQ=255)
                                                    tatTACCAAAGCGGATCTCCTCTCAGATCCGAACAAAATCACCCCAGTATTTGTACGTTTCTCTACCGtt                  >  1:3545106/1‑70 (MQ=255)
                                                    tatTACCAAAGCGGATCTCCTCTCAGATCCGAACAAAATCACCCCAGTATTTGTACGTTTCTCTACCGtt                  >  1:1219898/1‑70 (MQ=255)
                                                     atTACCAAAGCGGATTTCCTCTCAGATCCGAACAAAATCACCCCAGTATTTGTACGTTTCTCTACCGtt                  >  1:3765941/1‑69 (MQ=255)
                                                      ttACCAAAGCGGATTTCCTCTCAGATCCGAACAAAATCACCCCAGTATTTGTACGTTTCTCTACCGTTCa                <  1:1876306/70‑1 (MQ=255)
                                                      ttACCAAAGCGGATTTCCTCTCAGATCCGAACAAAATCACCCCAGTATTTGTACGTTTCTCTACCGTTCa                <  1:1503252/70‑1 (MQ=255)
                                                      ttACCAAAGCGGATTTCCTCTCAGATCCGAACAAAATCACCCCAGTATTTGTACGTTTCTCTACCGTTCa                <  1:2727012/70‑1 (MQ=255)
                                                      ttACCAAAGCGGATCTCCTCTCAGATCCGAACAAAATCACCCCAGTATTTGTACGTTTCTCTACCGTTCAg               >  1:3891156/1‑71 (MQ=255)
                                                      ttACCAAAGCGGATCTCCTCTCAGATCCGAACAAAATCACCCCAGTATTTGTACGTTTCTCTACCGTTCAg               >  1:3597227/1‑71 (MQ=255)
                                                       tACCAAAGCGGATCTCCTCTCAGATCCGAACAAAATCACCCCAGTATTTGTa                                 >  1:2231264/1‑52 (MQ=255)
                                                       tACCAAAGCGGATCTCCTCTCAGATCCGAACAAAATCACCCCAGTATTTGTa                                 >  1:3995852/1‑52 (MQ=255)
                                                       tACCAAAGCGGATCTCCTCTCAGATCCGAACAAAATCACCCCAGTATTTGTa                                 >  1:1423111/1‑52 (MQ=255)
                                                        aCCAAAGCGGATTTCCTCTCAGATCCGAACAAAATCACCCCAGTATTTGTACGTTTCTCTACCGTTCAggg             >  1:1363599/1‑71 (MQ=255)
                                                        aCCAAAGCGGATTTCATCTCAGATCCGAACAAAATCACCCCAGTATTTGTACGTTTCTCTACCGTTCAggg             >  1:1078876/1‑71 (MQ=255)
                                                        aCCAAAGCGGATCTCCTCTCAGATCCGAACAAAATCa                                               >  1:2873102/1‑37 (MQ=255)
                                                        aCCAAAGCGGATCTCCTCTCAGATCCGAACAAAATCACCCCAGTATTTGTACGTTTCTCTACCGTTCAggg             >  1:1069772/1‑71 (MQ=255)
                                                        aCCAAAGCGGATCTCCTCTCAGATCCGAACAAAATCACCCCAGTATTTGTACGTTTCTCTACCGTTCAggg             >  1:947060/1‑71 (MQ=255)
                                                        aCCAAAGCGGATCTCCTCTCAGATCCGAACAAAATCACCCCAGTATTTGTACGTTTCTCTACCGTTCAggg             >  1:3640966/1‑71 (MQ=255)
                                                        aCCAAAGCGGATCTCCTCTCAGATCCGAACAAAATCACCCCAGTATTTGTACGTTTCTCTACCGTTCAggg             >  1:2478900/1‑71 (MQ=255)
                                                             aGCGGATCTCCTCTCAGATCCGAACAAAATCACCCCAGTATTTGTACg                               >  1:3198540/1‑48 (MQ=255)
                                                                 gATTTCCTCTCAGATCCGAACAAAATCACCCCAGTATTTGTACGt                              <  1:3251686/45‑1 (MQ=255)
                                                                 gATTTCCTCTCAGATCCGAACAAAATCACCCCAGTATTTGTACGt                              <  1:3039210/45‑1 (MQ=255)
                                                                 gATTTCCTCTCAGATCCGAACAAAATCACCCCAGTATTTGTACGt                              <  1:2646213/45‑1 (MQ=255)
                                                                 gATTTCCTCTCAGATCCGAACAAAATCACCCCAGTATTTGTACGt                              <  1:1709155/45‑1 (MQ=255)
                                                                 gATCTCCTCTCAGATCCGAACAAAATCACCCCAGTATTTGTACGTTTCTCTACCGTTCAGGGTGGTGCTgg    >  1:2145075/1‑71 (MQ=255)
                                                                   tttCCTCTCAGATCCGAACAAAATCACCCCAGTAt                                      <  1:3169136/35‑1 (MQ=255)
                                                                   tttCCTCTCAGATCCGAACAAAATCACCCCAGTAt                                      <  1:340999/35‑1 (MQ=255)
                                                                   tttCCTCTCAGATCCGAACAAAATCACCCCAGTATTTGTACGTTTCTCTACCGTTCAGGGTGGTGCTGGct  >  1:2039035/1‑71 (MQ=255)
                                                                    |                                                                     
CACGCGGATCAGCCGCTCACGGTTATTTCCAGCCATATAAAAGCTTAAGCGATATTACCAAAGCGGATTTCCTCTCAGATCCGAACAAAATCACCCCAGTATTTGTACGTTTCTCTACCGTTCAGGGTGGTGCTGGCT  >  minE/1141943‑1142080

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: