Predicted mutation
evidence seq id position mutation freq annotation gene description
RA minE 1,142,011 T→C 26.2% 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,0110TC26.2% 65.0 / 27.7 42F152L (TTC→CTC) katEhydroperoxidase HPII(III)
Reads supporting (aligned to +/- strand):  ref base T (13/18);  new base C (5/6);  total (18/24)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 7.63e-01

GCCGCTCACGGTTATTTCCAGCCATATAAAAGCTTAAGCGATATTACCAAAGCGGATTTCCTCTCAGATCCGAACAAAATCACCCCAGTATTTGTACGTTTCTCTACCGTTCAGGGTGGTGCTGG  >  minE/1141954‑1142078
                                                         |                                                                   
gCCGCTCACGGTTATTTCCAGCCATATAAAAGCTTAAGCGATATTACCAAAGCGGATTTCCTCTCAGATcc                                                        >  1:932440/1‑71 (MQ=255)
gCCGCTCACGGTTATTTCCAGCCATATAAAAGCTTAAGCGATATTACCAAAGCGGATCTCCTCTCAGATcc                                                        <  1:50362/71‑1 (MQ=255)
 ccGCTCACGGTTATTTCCAGCCATATAAAAGCTTAAGCGATATTACCAAAGCGGATTTCCTCTCAg                                                            <  1:2204566/66‑1 (MQ=255)
    cTCACGGTTATTTCCAGCCATATAAAAGCTTAAGCGATATTACCAAAGCGGATTTCCtctc                                                              <  1:659687/61‑1 (MQ=255)
       aCGGTTATTTCCAGCCATATAAAAGCTTAAGCGATATTACCAAAGCGGATCTCCTCTCAGATCCGAACaaa                                                 <  1:1104871/71‑1 (MQ=255)
       aCGGTTATTTCCAGCCATATAAAAGCTTAAGCGATATTACCAAAGCGGATCTCCTCTCAGATCCGAACaa                                                  <  1:1188207/70‑1 (MQ=255)
        cGGTTATTTCCAGCCATATAAAAGCTTAAGCGATATTACCAAAGCGGATTTCCTCTCAGATCCGAACaaaa                                                >  1:1439121/1‑71 (MQ=255)
         ggTTATTTCCAGCCATATAAAAGCTTAAGCGATATTACCAAAGCGGATTTCCTCTCAGATCCGAACAAAAt                                               <  1:2027904/71‑1 (MQ=255)
             aTTTCCAGCCATATAAAAGCTTAAGCGATATTACCAAAGCGGATTTc                                                                   >  1:662444/1‑47 (MQ=255)
                tCCAGCCATATAAAAGCTTAAGCGATATTACCAAAGCGGATTTCCTCTCAGATCCGAACAAAATCAccc                                          >  1:2143475/1‑69 (MQ=255)
                  cAGCCATATAAAAGCTTAAGCGATATTACCAAAGCGGATTTCCTCTCAGATCCGAACAAAATCACCCCAg                                       <  1:2127698/70‑1 (MQ=255)
                    gCCATATAAAAGCTTAAGCGATATTACCAAAGCGGATTTCCTCTCAGATCCGAACAAAAt                                               <  1:1633224/60‑1 (MQ=255)
                        tataAAAGCTTAAGCGATATTACCAAAGCTGATTTCCTCTCAg                                                            <  1:398269/43‑1 (MQ=255)
                        tataAAAGCTTAAGCGATATTACCAAAGCGGATTTCctct                                                               <  1:514017/40‑1 (MQ=255)
                        tataAAAGCTTAAGCGATATTACCAAAGCGGATTTCctct                                                               <  1:942588/40‑1 (MQ=255)
                        tataAAAGCTTAAGCGATATTACCAAAGCGGATTTCCTCTCAg                                                            <  1:776631/43‑1 (MQ=255)
                        tataAAAGCTTAAGCGATATTACCAAAGCGGATTTCCTCTCAg                                                            <  1:148479/43‑1 (MQ=255)
                        tataAAAGCTTAAGCGATATTACCAAAGCGGATTTCCTCTCAg                                                            <  1:1778255/43‑1 (MQ=255)
                          taAAAGCTTAAGCGATATTACCAAAGCGGATTTCctct                                                               >  1:99136/1‑38 (MQ=255)
                          taAAAGCTTAAGCGATATTACCAAAGCGGATTTCCTCTCAGATCCGAACAAAATCAc                                            >  1:1924583/1‑57 (MQ=255)
                           aaaaGCTTAAGCGATATTACCAAAGCGGATTTCCTCTCAGATCCGAAc                                                    >  1:1349565/1‑48 (MQ=255)
                             aaGCTTAAGCGATATTACCAAAGCGGATCTCCTCTCAGATCCGAACAAAATCACCCCAGTATTTGTACGtt                           >  1:1468295/1‑71 (MQ=255)
                              aGCTTAAGCGATATTACCAAAGCGGATTTCCTCTCAGATCCGAACAAAATCACCCCAg                                       <  1:234172/58‑1 (MQ=255)
                               gCTTAAGCGATATTACCAAAGCGGATCTCCTCTCAGATCCGAACa                                                   >  1:555137/1‑45 (MQ=255)
                                 ttAAGCGATATTACCAAAGCGGATTTCCTCTCAGATCCGAACAAAAt                                               <  1:2086282/47‑1 (MQ=255)
                                     gCGATATTACCAAAGCGGATTTCCTCTCAGATCCGAACAAAATCACCCCAGTATTTGTACGTTTCTCTAcc                   >  1:2087631/1‑71 (MQ=255)
                                      cGATATTACCAAAGCGGATTTCCTCTCAGATCCGAACAAAATCACCCCAGTATTTGTACGTTTCTCTACCg                  <  1:358086/71‑1 (MQ=255)
                                          atTACCAAAGCGGATTTCCTCTCAGATCCGAACAAAATCACCCCAGTATTTGTACGTTTCTCTACCGtt                >  1:21753/1‑69 (MQ=255)
                                           ttACCAAAGCGGATTTCCTCTCAGATCCGAACAAAAt                                               <  1:1111339/37‑1 (MQ=255)
                                           ttACCAAAGCGGATTTCCTCTCAGATCCGAACAAAAt                                               <  1:411177/37‑1 (MQ=255)
                                           ttACCAAAGCGGATTTCCTCTCAGATCCGAACAAAAt                                               <  1:327390/37‑1 (MQ=255)
                                           ttACCAAAGCGGATTTCCTCTCAGATCCGAACAAAATCACCCCAGTATTTGTACGTTTCTCTACCGTTCAg             >  1:144684/1‑71 (MQ=255)
                                           ttACCAAAGCGGATTTCCTCTCAGATCCGAACAAAATCACCCCAGTATTTGTACGTTTCTCTACCGTTCAg             >  1:2030461/1‑71 (MQ=255)
                                           ttACCAAAGCGGATCTCCTCTCAGATCCGAACAAAATCA‑‑CCAGTATTTGTACGTTTCTCTACCGTTCAggg           >  1:422037/1‑71 (MQ=255)
                                              ccAAAGCGGATTTCCTCTCAGATCCGAACAAAATCa                                             >  1:64402/1‑36 (MQ=255)
                                              ccAAAGCGGATTTCCTCTCAGATCCGAACAAAATCa                                             >  1:36124/1‑36 (MQ=255)
                                              ccAAAGCGGATTTCCTCTCAGATCCGAACAAAATCACCCCAGTATTTg                                 <  1:115072/48‑1 (MQ=255)
                                               cAAAGCGGATCTCCTCTCAGATCCGAACAAAATCACCCCAGTATTTGTACGTTTCTCTACCGTTCAGggtg         <  1:1917510/71‑1 (MQ=255)
                                                 aaGCGGATCTCCTCTCAGATCCGAACAAAATCACCCCAGTATTTGTACGTTTCTCTACCGTTCAGggtggt       >  1:1815519/1‑71 (MQ=255)
                                                  aGCGGATCTCCTCTCAGATCCGAACAAAATCACCCCAGTATTTGTACGTTTCTCTACCGTTCAGGgtggtg      >  1:378540/1‑71 (MQ=255)
                                                      gATCTCCTCTCAGATCCGAACAAAATCACCCCAGTATTTGTACGTTTCTCTACCGTTCAGGGTGGTGCTgg  <  1:1590141/71‑1 (MQ=255)
                                                      gATCTCCTCTCAGATCCGAACAAAATCACCCCAGTATTTGTACGTTTCTCTACCGTTCAGGGTGGTGCTgg  <  1:1572166/71‑1 (MQ=255)
                                                        tttCCTCTCAGATCCGAACAAAATCACCCCAGTATTTGTACGTTtctc                       >  1:621141/1‑48 (MQ=255)
                                                         |                                                                   
GCCGCTCACGGTTATTTCCAGCCATATAAAAGCTTAAGCGATATTACCAAAGCGGATTTCCTCTCAGATCCGAACAAAATCACCCCAGTATTTGTACGTTTCTCTACCGTTCAGGGTGGTGCTGG  >  minE/1141954‑1142078

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: