Predicted mutation
evidence seq id position mutation freq annotation gene description
RA minE 1,142,011 T→C 16.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,0110TC16.0% 183.2 / 27.8 75F152L (TTC→CTC) katEhydroperoxidase HPII(III)
Reads supporting (aligned to +/- strand):  ref base T (33/30);  new base C (9/3);  total (42/33)
Fisher's exact test for biased strand distribution p-value = 2.09e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.55e-01

TGCACGCGGATCAGCCGCTCACGGTTATTTCCAGCCATATAAAAGCTTAAGCGATATTACCAAAGCGGATTTCCTCTCAGATCCGAACAAAATCACCCCAGTATTTGTACGTTTCTCTACCGTTCAGGGTGGTGCTGGCTC  >  minE/1141941‑1142081
                                                                      |                                                                      
tGCACGCGGATCAGCCGCTCACGGTTATTTCCAGCCATATAAAAGCTTAAGCGATATTACCAAAGCGGAtt                                                                        <  1:2226744/71‑1 (MQ=255)
                   cACGGTTATTTCCAGCCATATAAAAGCTTAAGCGATATTACCAAAGCGGATCTCCTCTCAGATCCGAACa                                                      >  1:3344302/1‑70 (MQ=255)
                    aCGGTTATTTCCAGCCATATAAAAGCTTAAGCGATATTACCAAAGCGGATTTCCTCTCAGATCCGAACaaa                                                    <  1:3113133/71‑1 (MQ=255)
                     cGGTTATTTCCAGCCATATAAAAGCTTAAGCGATATTACCAAAGCGGAtt                                                                        <  1:3663981/50‑1 (MQ=255)
                     cGGTTATTTCCAGCCATATAAAAGCTTAAGCGATATTACCAAAGCGGAtt                                                                        <  1:1364146/50‑1 (MQ=255)
                     cGGTTATTTCCAGCCATATAAAAGCTTAAGCGATATTACCAAAGCGGAtt                                                                        <  1:3880051/50‑1 (MQ=255)
                     cGGTTATTTCCAGCCATATAAAAGCTTAAGCGATATTACCAAAGCGGAtt                                                                        <  1:615430/50‑1 (MQ=255)
                     cGGTTATTTCCAGCCATATAAAAGCTTAAGCGATATTACCAAAGCGGAtt                                                                        <  1:2710121/50‑1 (MQ=255)
                       gTTATTTCCAGCCATATAAAAGCTTAAGCGATATTACCAAAGCGGATTTCCTCTCAGATCCGAACAAAAt                                                  >  1:2248304/1‑70 (MQ=255)
                       gTTATTTCCAGCCATATAAAAGCTTAAGCGATATTACCAAAGCGGATTTCCTCTCAGATCCGAACAAAAt                                                  >  1:1106455/1‑70 (MQ=255)
                        ttATTTCCAGCCATATAAAAGCTTAAGCGATATTACCAAAGCGGAtt                                                                        <  1:4104017/47‑1 (MQ=255)
                        ttATTTCCAGCCATATAAAAGCTTAAGCGATATTACCAAAGCGGAtt                                                                        <  1:3910457/47‑1 (MQ=255)
                        ttATTTCCAGCCATATAAAAGCTTAAGCGATATTACCAAAGCGGATTTc                                                                      <  1:3037436/49‑1 (MQ=255)
                        ttATTTCCAGCCATATAAAAGCTTAAGCGATATTACCAAAGCGGATTTc                                                                      <  1:689182/49‑1 (MQ=255)
                        ttATTTCCAGCCATATAAAAGCTTAAGCGATATTACCAAAGCGGATTTc                                                                      <  1:2169696/49‑1 (MQ=255)
                          aTTTCCAGCCATATAAAAGCTTAAGCGATATTACCAAAGCGGATTTCCTCTCAGATCCGAACAAAATCAcc                                              >  1:3858612/1‑71 (MQ=255)
                          aTTTCCAGCCATATAAAAGCTTAAGCGATATTACCAAAGCGGATTTCCTCTCAGATCCGAACAAAATCAcc                                              >  1:1281412/1‑71 (MQ=255)
                          aTTTCCAGCCATATAAAAGCTTAAGCGATATTACCAAAGCGGATTTCCTCTCAGATCCGAACAAAATCAcc                                              >  1:2015859/1‑71 (MQ=255)
                            ttCCAGCCATATAAAAGCTTAAGCGATATTACCAAAGCGGATTTCCTCTCAGATCCGAACAAAATCAccc                                             >  1:1696123/1‑70 (MQ=255)
                            ttCCAGCCATATAAAAGCTTAAGCGATATTACCAAAGCGGATTTCCTCTCAGATCCGAACAAAATCAccc                                             >  1:230475/1‑70 (MQ=255)
                            ttCCAGCCATATAAAAGCTTAAGCGATATTACCAAAGCGGATTTCCTCTCAGATCCGAACAAAATCAcc                                              >  1:1611650/1‑69 (MQ=255)
                               cAGCCATATAAAAGCTTAAGCGATATTACCAAAGCGGATTTCCTCTCAGATCCGAACAAAATCACCCCAg                                          <  1:3445696/70‑1 (MQ=255)
                               cAGCCATATAAAAGCTTAAGCGATATTACCAAAGCGGATTTCCTCTCAGATCCGAACAAAATCACCCCAGt                                         <  1:2777581/71‑1 (MQ=255)
                               cAGCCATATAAAAGCTTAAGCGATATTACCAAAGCGGATTTCCTCTCAGATCCGAACAAAATCACCCCAGt                                         <  1:1982637/71‑1 (MQ=255)
                                    atatAAAAGCTTAAGCGATATTACCAAAGCGGATTTCCTCTCAGATCCGAACAAAATCACCCCAGTATTTg                                    <  1:467282/71‑1 (MQ=255)
                                    atatAAAAGCTTAAGCGATATTACCAAAGCGGATTTCCTCTCAGATCCGAACAAAATCACCCCAGTATTTg                                    <  1:2313301/71‑1 (MQ=255)
                                    atatAAAAGCTTAAGCGATATTACCAAAGCGGATCTCCTCTCAGATCCGAACAAAATCACCCCAGTATTTg                                    >  1:2096218/1‑71 (MQ=255)
                                    atatAAAAGCTTAAGCGATATTACCAAAGCGGATCTCCTCTCAGATCCGAACAAAATCACCCCAGTATTTg                                    >  1:2109987/1‑71 (MQ=255)
                                     tataAAAGCTTAAGCGATATTACCAAAGCGTATTTCCTCTCAGATCCGAACaa                                                     <  1:817155/53‑1 (MQ=255)
                                       taAAAGCTTAAGCGATATTACCAAAGCGGATTTCCTCTCAGATCCGAACaaa                                                    <  1:1189049/52‑1 (MQ=255)
                                       taAAAGCTTAAGCGATATTACCAAAGCGGATTTCCTCTCAGATCCGAACaaa                                                    <  1:749353/52‑1 (MQ=255)
                                       taAAAGCTTAAGCGATATTACCAAAGCGGATTTCCTCTCAGATCCGAACAAAATCa                                                >  1:2811299/1‑56 (MQ=255)
                                       taAAAGCTTAAGCGATATTACCAAAGCGGATTTCCTCTCAGATCCGAACAAAATCa                                                >  1:2791158/1‑56 (MQ=255)
                                       taAAAGCTTAAGCGATATTACCAAAGCGGATCTcc                                                                     >  1:2248934/1‑35 (MQ=255)
                                       taAAAGCTTAAGCGATATTACCAAAGCGGATCTcc                                                                     >  1:598112/1‑35 (MQ=255)
                                       taAAAGCTTAAGCGATATTACCAAAGCGGATCTcc                                                                     >  1:2602205/1‑35 (MQ=255)
                                        aaaaGCTTAAGCGATATTACCAAAGCGGATTTCCTCTCAGATCCGAACaaaa                                                   >  1:2226754/1‑52 (MQ=255)
                                        aaaaGCTTAAGCGATATTACCAAAGCGGATTTCCTCTCAGATCCGAACaaaa                                                   >  1:1951658/1‑52 (MQ=255)
                                         aaaGCTTAAGCGATATTACCAAAGCGGATTTCCTCTCa                                                                >  1:2785450/1‑38 (MQ=255)
                                          aaGCTTAAGCGATATTACCAAAGCGGATTTCCTCTCAGATCCGAACaaa                                                    <  1:3249084/49‑1 (MQ=255)
                                          aaGCTTAAGCGATATTACCAAAGCGGATTTCCTCTCAGATCCGAACaaa                                                    <  1:717213/49‑1 (MQ=255)
                                          aaGCTTAAGCGATATTACCAAAGCGGATTTCCTCTCAGATCCGAACAAAATCACCCCAGTATTTGTACGtt                              >  1:1520677/1‑71 (MQ=255)
                                            gCTTAAGCGATATTACCAAAGCGGATTTCCTCTCAGATCCGAACAAAAt                                                  >  1:2772710/1‑49 (MQ=255)
                                            gCTTAAGCGATATTACCAAAGCGGATTTCCTCTCAGATCCGAACAAAAt                                                  >  1:2036370/1‑49 (MQ=255)
                                            gCTTAAGCGATATTACCAAAGCGGATTTCCTCTCAGATCCGAACAAAAt                                                  >  1:646954/1‑49 (MQ=255)
                                            gCTTAAGCGATATTACCAAAGCGGATTTCCTCTCAGATCCGAACAAAAt                                                  >  1:3030536/1‑49 (MQ=255)
                                            gCTTAAGCGATATTACCAAAGCGGATTTCCTCTCAGATCCGAACAAAAt                                                  >  1:2173918/1‑49 (MQ=255)
                                            gCTTAAGCGATATTACCAAAGCGGATTTCCTCTCAGATCCGAACAAAAt                                                  >  1:1621134/1‑49 (MQ=255)
                                            gCTTAAGCGATATTACCAAAGCGGATTTCCTCTCAGATCCGAACAAAAt                                                  >  1:146008/1‑49 (MQ=255)
                                            gCTTAAGCGATATTACCAAAGCGGATTTCCTCTCAGATCCGAACAAAATCACCCCAGTATTTGTACGTTtc                            <  1:442209/71‑1 (MQ=255)
                                            gCTTAAGCGATATTACCAAAGCGGATTTCCTCTCAGATCCGAACAAAATCACCCCAGTATTTGTACGTTtc                            <  1:1489613/71‑1 (MQ=255)
                                                aaGCGATATTACCAAAGCGGATCTCCTCTCAGATCCGAACAAAATCACCCCAGTATTTGTACGTTTCTCTa                        >  1:260372/1‑71 (MQ=255)
                                                aaGCGATATTACCAAAGCGGATCTCCTCTCAGATCCGAACAAAATCACCCCAGTATTTGTACGTTTCTCTa                        >  1:2545275/1‑71 (MQ=255)
                                                  tcGATATTACCAAAGCGGATTTCCTCTCAGATCCGAACAAAATCACCCCAGTATTTGTACGTTTCTCTAcc                      <  1:3038394/70‑1 (MQ=255)
                                                  tcGATATTACCAAAGCGGATTTCCTCTCAGATCCGAACAAAATCACCCCAGTATTTGTACGTTTCTCTAcc                      <  1:820546/70‑1 (MQ=255)
                                                   cGATATTACCAAAGCGGATTTCCTCTCAGATCCGAACAAAATCACCCCAGTATTTGTACGTTTctct                         <  1:3260334/67‑1 (MQ=255)
                                                   cGATATTACCAAAGCGGATTTCCTCTCAGATCCGAACAAAATCACCCCAGTATTTGTACGTTTctct                         <  1:3711430/67‑1 (MQ=255)
                                                   cGATATTACCAAAGCGGATTTCCTCTCAGATCCGAACAAAATCACCCCAGTATTTGTACGTTTctct                         <  1:642482/67‑1 (MQ=255)
                                                   cGATATTACCAAAGCGGATTTCCTCTCAGATCCGAACAAAATCACCCCAGTATTTGTACGTTTctct                         <  1:1954169/67‑1 (MQ=255)
                                                   cGATATTACCAAAGCGGATTTCCTCTCAGATCCGAACAAAATCACCCCAGTATTTGTACGTTTctct                         <  1:3442808/67‑1 (MQ=255)
                                                      tatTACCAAAGCGGATTTCCTCTTAGATCCGAACAAAATCACCCCAGTAt                                       >  1:1695673/1‑50 (MQ=255)
                                                      tatTACCAAAGCGGATTTCCTCTCAGATCCGAACAAAATCACCCCAGTAt                                       >  1:278948/1‑50 (MQ=255)
                                                      tatTACCAAAGCGGATTTCCTCTCAGATCCGAACAAAATCACCCCAGTATTTGTACGTTTCTCTACCGTTc                  <  1:813205/71‑1 (MQ=255)
                                                      tatTACCAAAGCGGATTTCCTCTCAGATCCGAACAAAATCACCCCAGTATTTGTACGTTTCTCTACCGTTc                  <  1:1892853/71‑1 (MQ=255)
                                                       atTACCAAAGCGGATTTCCTCTCAGATCCGAACAAAATCACCCCAGTATTTGTACGt                               <  1:331141/57‑1 (MQ=255)
                                                       atTACCAAAGCGGATTTCCTCTCAGATCCGAACAAAATCACCCCAGTATTTGTACGTTTCTCTACCGTTc                  <  1:3904648/70‑1 (MQ=255)
                                                          aCCAAAGCGGATTTCCTCTCAGATCCGAACAAAATCCCCCCGGttt                                       >  1:3519018/1‑44 (MQ=255)
                                                          aCCAAAGCGGATTTCCTCTCAGATCCGAACAAAATCACCCCAGTATTTGTACGTTTCTCTACCGTTCAggg              >  1:3487546/1‑71 (MQ=255)
                                                          aCCAAAGCGGATTTCCTCTCAGATCCGAACAAAATCACCCCAGTATTTGTACGTTTCTCTACCGTTCAggg              >  1:391804/1‑71 (MQ=255)
                                                          aCCAAAGCGGATTTCCTCTCAGATCCGAACAAAATCACCCCAGTATTTGTACGTTTCTCTACCGTTCAgg               >  1:1142566/1‑70 (MQ=255)
                                                          aCCAAAGCGGATCTCCTCTCAGATCCGAACAAAATCACCCCAGTATTTGTACGTTTCTCTACCGTTCAggg              >  1:2385110/1‑71 (MQ=255)
                                                           ccAAAGCGGATTTCCTCTCAGATCCGAACAAAATCACCCCAGTATTTGTACGTTTCTCTACCGTTCAggg              <  1:3948875/70‑1 (MQ=255)
                                                           ccAAAGCGGATTTCCTCTCAGATCCGAACAAAATCACCCCAGTATTTGTACGTTTCTCTACCGTTCAggg              <  1:2004118/70‑1 (MQ=255)
                                                           ccAAAGCGGATTTCCTCTCAGATCCGAACAAAATCACCCCAGTATTTGTACGTTTCTCTACCGTTCAggg              <  1:2883748/70‑1 (MQ=255)
                                                              aaGCGGATTTCCTCTCAGATCCGAACAAAATCACCCCAGTATTTGTACGTTTCTCTACCGTTCAGggtggt          >  1:3191594/1‑71 (MQ=255)
                                                              aaGCGGATTTCCTCTCAGATCCGAACAAAATCACCCCAGTATTTGTACGTTTCTCTACCGTTCAGggtggt          >  1:3236460/1‑71 (MQ=255)
                                                              aaGCGGATTTCCTCTCAGATCCGAACAAAATCACACCAGTATTTGTACGTTTCTCTACCGTTCAGggtggt          >  1:4044851/1‑71 (MQ=255)
                                                               aGCGGATTTCCTCTCAGATCCGAACAAAATCACCCCAGTATTTGTACGTTTCTCTACCGTTCAGggtg            >  1:2646017/1‑68 (MQ=255)
                                                               aGCGGATTTCCTCTCAGATCCGAACAAAATCACCCCAGTATTTGTACGTTTCTCTACCGTTCAGggtg            >  1:3130134/1‑68 (MQ=255)
                                                               aGCGGATTTCCTCTCAGATCCGAACAAAATAACCCCAGTATTTGTACGTTTCTCTACCGTTCAGggtg            >  1:3534319/1‑68 (MQ=255)
                                                                 cGGATCTCCTCTCAGATCCGAACAAAATCACCCCAGTATTTGTACGt                               <  1:2528939/47‑1 (MQ=255)
                                                                 cGGATCTCCTCTCAGATCCGAACAAAATCACCCCAGTATTTGTACGt                               <  1:1920325/47‑1 (MQ=255)
                                                                 cGGATCTCCTCTCAGATCCGAACAAAATCACCCCAGTATTTGTACGt                               <  1:1643349/47‑1 (MQ=255)
                                                                      ttCCTCTCAGATCCGAACAAAATCACCCCAGTATTTGTACGTTTCTCTACCGTTCAGGGTGGTGCTGGctc  >  1:2937223/1‑71 (MQ=255)
                                                                      ttCCTCTCAGATCCGAACAAAATCACCCCAGTATTTGTACGTTTCTCTACCGTTCAGGGTGGTGCTGGctc  >  1:1881744/1‑71 (MQ=255)
                                                                      ttCCTCTCAGATCCGAACAAAATCACCCCAGTATTTGTACGTTTCTCTACCGTTCAGGGTGGTGCTGGctc  >  1:1782531/1‑71 (MQ=255)
                                                                      ttCCTCTCAGATCCGAACAAAATCACCCCAGTATTTGTACGTTTCTCTACCGTTCAGGGTGGTGCTGGctc  >  1:3326630/1‑71 (MQ=255)
                                                                      |                                                                      
TGCACGCGGATCAGCCGCTCACGGTTATTTCCAGCCATATAAAAGCTTAAGCGATATTACCAAAGCGGATTTCCTCTCAGATCCGAACAAAATCACCCCAGTATTTGTACGTTTCTCTACCGTTCAGGGTGGTGCTGGCTC  >  minE/1141941‑1142081

Alignment Legend
Aligned base mismatch/match (shaded by quality score): ATCG/ATCG < 3 ≤ ATCG/ATCG < 32 ≤ ATCG/ATCG < 36 ≤ ATCG/ATCG
Unaligned base: atcg    Masked matching base: atcg    Alignment gap:     Deleted base: