Predicted mutation
evidence seq id position mutation freq annotation gene description
RA minE 2,334,294 C→T 30.5% T581T (ACC→ACT fdoG → formate dehydrogenase‑O, large subunit

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*minE2,334,2940CT30.5% 124.0 / 56.9 121T581T (ACC→ACTfdoGformate dehydrogenase‑O, large subunit
Reads supporting (aligned to +/- strand):  ref base C (34/50);  new base T (26/11);  total (60/61)
Fisher's exact test for biased strand distribution p-value = 3.09e-03
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 8.50e-01

Only 100 of 121 total aligned reads displayed.

GTTGCCTCATTCCCGAACAAAAACAAAGTGATCGGTTGTCTGTCGAAACTGAAGTTCCTCGTCACCATCGACCCACTGAACACTGAAACCTCTAACTTCTGGCAGAACCACGGTGAGCTGAACGA  >  minE/2334229‑2334353
                                                                 |                                                           
gTTGCCTCATTCCCGAACAAAAACAAAGTGATCGGTTGTCTGTCGAAACTGAAGTTCCTCGTCACCATc                                                          >  1:938734/1‑69 (MQ=255)
gTTGCCTCATTCCCGAACAAAAACAAAGTGATCGGTTGTCTGTCGAAACTGAAGTTCCTCGTCACCATc                                                          >  1:214170/1‑69 (MQ=255)
gTTGCCTCATTCCCGAACAAAAACAAAGTGATCGGTTGTCTGTCGAAACTGAAGTTCCTCGTCACCATc                                                          >  1:2203081/1‑69 (MQ=255)
gTTGCCTCATTCCCGAACAAAAACAAAGTGATCGGTTGTCTGTCGAAACTGAAGTTCCTCGTCACCATc                                                          >  1:1923357/1‑69 (MQ=255)
gTTGCCTCATTCCCGAACAAAAACAAAGTGATCGGTTGTCTGTCGAAACTGAAGTTCCTCGTCACCATc                                                          >  1:2788809/1‑69 (MQ=255)
gTTGCCTCATTCCCGAACAAAAACAAAGTGATCGGTTGTCTGTCGAAACTGAAGTTCCTCGTCACCATc                                                          >  1:1314239/1‑69 (MQ=255)
gTTGCCTCATTCCCGAACAAAAACAAAGTGATCGGTTGTCTGTCGAAACTGAAGTTCCTCGTCACCATc                                                          >  1:3102112/1‑69 (MQ=255)
gTTGCCTCATTCCCGAACAAAAACAAAGTGATCGGTTGTCTGTCGAAACTGAAGTTCCTCGTCACCATc                                                          >  1:447790/1‑69 (MQ=255)
   gCCTCATTCCCGAACAAAAACAAAGTGATCGGTTGTCTGTCGAAACTGAAGTTCCTCGTCACCATCGa                                                        >  1:2572183/1‑68 (MQ=255)
   gCCTCATTCCCGAACAAAAACAAAGTGATCGGTTGTCTGTCGAAACTGAAGTTCCTCGTCACCATCGa                                                        >  1:1673502/1‑68 (MQ=255)
   gCCTCATTCCCGAACAAAAACAAAGTGATCGGTTGTCTGTCGAAACTGAAGTTCCTCGTCACCATCGa                                                        >  1:1123520/1‑68 (MQ=255)
   gCCTCATTCCCGAACAAAAACAAAGTGATCGGTTGTCTGTCGAAACTGAAGTTCCTCGTCACCATCGa                                                        >  1:1553607/1‑68 (MQ=255)
   gCCTCATTCCCGAACAAAAACAAAGTGATCGGTTGTCTGTCGAAACTGAAGTTCCTCGTCACCATCGa                                                        >  1:135632/1‑68 (MQ=255)
      tCATTCCCGAACAAAAACAAAGTGATCGGTTGTCTGTCGAAACTGAAGTTCCTCGTCACCATCGACCca                                                    >  1:1225997/1‑69 (MQ=255)
      tCATTCCCGAACAAAAACAAAGTGATCGGTTGTCTGTCGAAACTGAAGTTCCTCGTCACCATCGACCca                                                    >  1:2779816/1‑69 (MQ=255)
      tCATTCCCGAACAAAAACAAAGTGATCGGTTGTCTGTCGAAACTGAAGTTCCTCGTCACCATCGACCca                                                    >  1:2501252/1‑69 (MQ=255)
      tCATTCCCGAACAAAAACAAAGTGATCGGTTGTCTGTCGAAACTGAAGTTCCTCGTCACCATCGACCca                                                    >  1:2791965/1‑69 (MQ=255)
      tCATTCCCGAACAAAAACAAAGTGATCGGTTGTCTGTCGAAACTGAAGTTCCTCGTCACCATCGACCca                                                    >  1:138028/1‑69 (MQ=255)
           cccGAACAAAAACAAAGTGATCGGTTGTCTGTCGAAACTGAAGTTCCTCGTCACCATCGACCcactgaa                                               <  1:2421091/69‑1 (MQ=255)
           cccGAACAAAAACAAAGTGATCGGTTGTCTGTCGAAACTGAAGTTCCTCGTCACCATCGACCcactgaa                                               <  1:94509/69‑1 (MQ=255)
           cccGAACAAAAACAAAGTGATCGGTTGTCTGTCGAAACTGAAGTTCCTCGTCACCATCGACCcactgaa                                               <  1:2411085/69‑1 (MQ=255)
           cccGAACAAAAACAAAGTGATCGGTTGTCTGTCGAAACTGAAGTTCCTCGTCACCATCGACCcactgaa                                               <  1:239232/69‑1 (MQ=255)
           cccGAACAAAAACAAAGTGATCGGTTGTCTGTCGAAACTGAAGTTCCTCGTCACCATCGACCcactgaa                                               <  1:2181087/69‑1 (MQ=255)
           cccGAACAAAAACAAAGTGATCGGTTGTCTGTCGAAACTGAAGTTCCTCGTCACCATCGACCcactgaa                                               <  1:3274371/69‑1 (MQ=255)
           cccGAACAAAAACAAAGTGATCGGTTGTCTGTCGAAACTGAAGTTCCTCGTCACCATCGACCcactgaa                                               <  1:2599246/69‑1 (MQ=255)
           cccGAACAAAAACAAAGTGATCGGTTGTCTGTCGAAACTGAAGTTCCTCGTCACCATCGACCcactgaa                                               <  1:43196/69‑1 (MQ=255)
           cccGAACAAAAACAAAGTGATCGGTTGTCTGTCGAAACTGAAGTTCCTCGTCACCATCGACCcactgaa                                               <  1:2867517/69‑1 (MQ=255)
           cccGAACAAAAACAAAGTGATCGGTTGTCTGTCGAAACTGAAGTTCCTCGTCACCATCGACCcactgaa                                               <  1:3370006/69‑1 (MQ=255)
           cccGAACAAAAACAAAGTGATCGGTTGTCTGTCGAAACTGAAGTTCCTCGTCACCATCGACCcactgaa                                               <  1:2887270/69‑1 (MQ=255)
           cccGAACAAAAACAAAGTGATCGGTTGTCTGTCGAAACTGAAGTTCCTCGTCACCATCGACCcactgaa                                               <  1:2970074/69‑1 (MQ=255)
           cccGAACAAAAACAAAGTGATCGGTTGTCTGTCGAAACTGAAGTTCCTCGTCACCATCGACCcactgaa                                               <  1:3107038/69‑1 (MQ=255)
           cccGAACAAAAACAAAGTGATCGGTTGTCTGTCGAAACTGAAGTTCCTCGTCACCATCGACCcactgaa                                               <  1:3152646/69‑1 (MQ=255)
           cccGAACAAAAACAAAGTGATCGGTTGTCTGTCGAAACTGAAGTTCCTCGTCACCATCGACCcactgaa                                               <  1:3233747/69‑1 (MQ=255)
           cccGAACAAAAACAAAGTGATCGGTTGTCTGTCGAAACTGAAGTTCCTCGTCACCATCGACCcactgaa                                               <  1:14130/69‑1 (MQ=255)
           cccGAACAAAAACAAAGTGATCGGTTGTCTGTCGAAACTGAAGTTCCTCGTCACCATCGACCcactgaa                                               >  1:2041725/1‑69 (MQ=255)
           cccGAACAAAAACAAAGTGATCGGTTGTCTGTCGAAACTGAAGTTCCTCGTCACCATCGACCcactgaa                                               <  1:1282944/69‑1 (MQ=255)
           cccGAACAAAAACAAAGTGATCGGTTGTCTGTCGAAACTGAAGTTCCTCGTCACCATCGACCcactgaa                                               <  1:1105178/69‑1 (MQ=255)
           cccGAACAAAAACAAAGTGATCGGTTGTCTGTCGAAACTGAAGTTCCTCGTCACCATCGACCcactgaa                                               <  1:1227893/69‑1 (MQ=255)
             cGAACAAAAACAAAGTGATCGGTTGTCTGGCGAAACTGAAGTTCCTCGTCACCATCGACCcactgaa                                               <  1:3069008/67‑1 (MQ=255)
                acaaaaacaaaGTGATCGGCTGTCTGTCGAAGCTGAAATTCCTCGTGACTAt                                                           >  1:2320308/1‑52 (MQ=25)
                acaaaaacaaaGTGATCGGCTGTCTGTCGAAGCTGAAATTCCTCGTGACTATc                                                          >  1:1135172/1‑53 (MQ=25)
                acaaaaacaaaGTGATCGGCTGTCTGTCGAAGCTGAAATTCCTCGTGACTATc                                                          >  1:2736589/1‑53 (MQ=25)
                acaaaaacaaaGTGATCGGCTGTCTGTCGAAGCTGAAATTCCTCGTGACTATc                                                          >  1:2781325/1‑53 (MQ=25)
                acaaaaacaaaGTGATCGGCTGTCTGTCGAAGCTGAAATTCCTCGTGACTATc                                                          >  1:1343472/1‑53 (MQ=25)
                acaaaaacaaaGTGATCGGCTGTCTGTCGAAGCTGAAATTCCTCGTGACTATc                                                          >  1:2958475/1‑53 (MQ=25)
                acaaaaacaaaGTGATCGGCTGTCTGTCGAAGCTGAAATTCCTCGTGACTATc                                                          >  1:3342141/1‑53 (MQ=25)
                acaaaaacaaaGTGATCGGCTGTCTGTCGAAGCTGAAATTCCTCGTGACTATc                                                          >  1:3394846/1‑53 (MQ=25)
                acaaaaacaaaGTGATCGGCTGTCTGTCGAAGCTGAAATTCCTCGTGACTATc                                                          >  1:2104780/1‑53 (MQ=25)
                acaaaaacaaaGTGATCGGCTGTCTGTCGAAGCTGAAATTCCTCGTGACTATc                                                          >  1:456776/1‑53 (MQ=25)
                acaaaaacaaaGTGATCGGCTGTCTGTCGAAGCTGAAATTCCTCGTGACTATc                                                          >  1:465582/1‑53 (MQ=25)
                acaaaaacaaaGTGATCGGCTGTCTGTCGAAGCTGAAATTCCTCGTGACTATc                                                          >  1:477556/1‑53 (MQ=25)
                acaaaaacaaaGTGATCGGCTGTCTGTCGAAGCTGAAATTCCTCGTGACTATc                                                          >  1:549253/1‑53 (MQ=25)
                acaaaaacaaaGTGATCGGCTGTCTGTCGAAGCTGAAATTCCTCGTGACTATc                                                          >  1:622798/1‑53 (MQ=25)
                acaaaaacaaaGTGATCGGCTGTCTGTCGAAGCTGAAATTCCTCGTGACTATc                                                          >  1:775449/1‑53 (MQ=25)
                acaaaaacaaaGTGATCGGCTGTCTGTCGAAGCTGAAATTCCTCGTGACTATc                                                          >  1:1063014/1‑53 (MQ=25)
                acaaaaacaaaGTGATCGGCTGTCTGTCGAAGCTGAAATTCCTCGTGACTATc                                                          >  1:2689155/1‑53 (MQ=25)
                acaaaaacaaaGTGATCGGCTGTCTGTCGAAGCTGAAATTCCTCGTGACTATc                                                          >  1:2691943/1‑53 (MQ=25)
                acaaaaacaaaGTGATCGGCTGTCTGTCGAAGCTGAAATTCCTCGTGACTATc                                                          >  1:2124770/1‑53 (MQ=25)
                acaaaaacaaaGTGATCGGCTGTCTGTCGAAGCTGAAATTCCTCGTGACTATc                                                          >  1:2276925/1‑53 (MQ=25)
                acaaaaacaaaGTGATCGGCTGTCTGTCGAAGCTGAAATTCCTCGTGACTATc                                                          >  1:1680172/1‑53 (MQ=25)
                acaaaaacaaaGTGATCGGCTGTCTGTCGAAGCTGAAATTCCTCGTGACTATc                                                          >  1:1727682/1‑53 (MQ=25)
                acaaaaacaaaGTGATAGGCTGTCTGTCGAAGCTGAAATTCCTCGTGACTATc                                                          >  1:1022315/1‑53 (MQ=25)
                 caaaaacaaaGTGATCGGTTGTCTGTCGAAACTGAAGTTCCTCGTCACCATCGACCcactgaacactga                                         <  1:700780/69‑1 (MQ=255)
                 caaaaacaaaGTGATCGGTTGTCTGTCGAAACTGAAGTTCCTCGTCACCATCGACCcactgaacactga                                         <  1:2116289/69‑1 (MQ=255)
                  aaaaacaaaGTGTTCGGTTGTCTGTCGAAACTGAAGTTCCTCGTCACCATCGACCcactgaacactgaa                                        <  1:2449073/69‑1 (MQ=255)
                  aaaaacaaaGTGATCGGTTGTCTGTCGAAACTGAAGTTCCTCGTCACCATCGACCcactgaacactgaa                                        <  1:3272354/69‑1 (MQ=255)
                  aaaaacaaaGTGATCGGTTGTCTGTCGAAACTGAAGTTCCTCGTCACCATCGACCcactgaacactgaa                                        <  1:1471083/69‑1 (MQ=255)
                  aaaaacaaaGTGATCGGTTGTCTGTCGAAACTGAAGTTCCTCGTCACCATCGACCcactgaacactgaa                                        <  1:2061423/69‑1 (MQ=255)
                  aaaaacaaaGTGATCGGTTGTCTGTCGAAACTGAAGTTCCTCGTCACCATCGACCcactgaacactgaa                                        <  1:2677807/69‑1 (MQ=255)
                  aaaaacaaaGTGATCGGTTGTCTGTCGAAACTGAAGTTCCTCGTCACCATCGACCcactgaacactgaa                                        <  1:1353794/69‑1 (MQ=255)
                  aaaaacaaaGTGATCGGTTGTCTGTCGAAACTGAAGTTCCTCGTCACCATCGACCcactgaacactgaa                                        <  1:908086/69‑1 (MQ=255)
                  aaaaacaaaGTGATCGGTTGTCTGTCGAAACTGAAGTTCCTCGTCACCATCGACCcactgaacactgaa                                        <  1:91015/69‑1 (MQ=255)
                  aaaaacaaaGTGATCGGTTGTCTGTCGAAACTGAAGTTCCTCGTCACCATCGACCcactgaacactgaa                                        <  1:1766353/69‑1 (MQ=255)
                  aaaaacaaaGTGATCGGTTGTCTGTCGAAACTGAAGTTCCTCGTCACCATCGACCcactgaacactgaa                                        <  1:2071884/69‑1 (MQ=255)
                  aaaaacaaaGTGATCGGTTGTCTGTCGAAACTGAAGTTCCTCGTCACCATCGACCcactgaacactgaa                                        <  1:2922932/69‑1 (MQ=255)
                  aaaaacaaaGTGATCGGATGTCTGTCGAAACTGAAGTTCCTCGTCACCATCGACCcactgaacactgaa                                        <  1:3217355/69‑1 (MQ=255)
                  aaaaacaaaGCGATCGGTTGTCTGTCGAAACTGAAGTTCCTCGTCACCATCGACCcactgaacactgaa                                        <  1:1284047/69‑1 (MQ=255)
                    aaacaaaGTGATCGGTTGTCTGTCGAAACTGAAGTTCCTCGTCACCATCGACCCACTGAACACTGAAAc                                      <  1:1235507/69‑1 (MQ=255)
                     aacaaaGTGATCGGTTGTCTGTCGAAACTGAAGTTCCTCGTCACCa                                                            >  1:1191924/1‑46 (MQ=255)
                     aacaaaGTGATCGGTTGTCTGTCGAAACTGAAGTTCCTCGTCACCATCGACCCACTGAACACTGAAAcc                                     >  1:1230281/1‑69 (MQ=255)
                     aacaaaGTGATCGGTTGTCTGTCGAAACTGAAGTTCCTCGTCACCATCGACCCACTGAACACTGAAAc                                      <  1:1810709/68‑1 (MQ=255)
                     aacaaaGTGATCGGTTGTCTGTCGAAACTGAAGTTCCTCGTCACCATCGACCCACTGAACACTGAAAc                                      >  1:921091/1‑68 (MQ=255)
                        aaaGTGATCGGTTGTCTGTCGAAACTGAAGTTCCTCGTCACCATCGACCcactgaa                                               <  1:3412973/56‑1 (MQ=255)
                            tGATCGGTTGTCTGTCGAAACTGAAGTTCCTCGTCACCATCGACCCACTGAACACTGAAACCTCTAACt                              >  1:843076/1‑69 (MQ=255)
                            tGATCGGTTGTCTGTCGAAACTGAAGTTCCTCGTCACCATCGACCCACTGAACACTGAAACCTCTAACt                              >  1:1731415/1‑69 (MQ=255)
                                    tgtctgtcGAAACTGAAGTTCCTCGTCACCATCGACCCACTGAACACTGAAACCTCTAACTTCTGGCAg                      <  1:2355332/69‑1 (MQ=255)
                                    tgtctgtcGAAACTGAAGTTCCTCGTCACCATCGACCCACTGAACACTGAAACCTCTAACTTCTGGCAg                      <  1:900187/69‑1 (MQ=255)
                                    tgtctgtcGAAACTGAAGTTCCTCGTCACCATCGACCCACTGAACACTGAAACCTCTAACTTCTGGCAg                      <  1:213446/69‑1 (MQ=255)
                                    tgtctgtcGAAACTGAAGTTCCTCGTCACCATCGACCCACTGAACACTGAAACCTCTAACTTCTGGCAg                      <  1:1800910/69‑1 (MQ=255)
                                    tgtctgtcGAAACTGAAGTTCCTCGTCACCATCGAACCACTGAACACTGAAACCTCTAACTTCTGGCAg                      <  1:828525/69‑1 (MQ=255)
                                              agcTGAAATTCCTCGTGACTATCGACCCGCTGAACACTGAAACCTCTAACTTCTTGCAGAACCACGGt             <  1:2003894/66‑1 (MQ=18)
                                              agcTGAAATTCCTCGTGACTATCGACCCGCTGAACACTGAAACCTCTAACTTCTGGCAGAACCACGGt             <  1:2217535/66‑1 (MQ=18)
                                              agcTGAAATTCCTCGTGACTATCGACCCGCTGAACACTGAAACCTCTAACTTCTGGCAGAACCACGGt             <  1:2373874/66‑1 (MQ=18)
                                              agcTGAAATTCCTCGTGACTATCGACCCGCTGAACACTGAAACCTCTAACTTCTGGCAGAACCACGGt             <  1:847326/66‑1 (MQ=18)
                                              agcTGAAATTCCTCGTGACTATCGACCCGCTGAACACTGAAACCTCTAACTTCTGGCAGAACCACGGt             <  1:1258429/66‑1 (MQ=18)
                                              agcTGAAATTCCTCGTGACTATCGACCCGCTGAACACTGAAACCTCTAACTTCTGGCAGAACCACGGt             <  1:1272905/66‑1 (MQ=18)
                                              agcTGAAATTCCTCGTGACTATCGACCCGCTGAACACTGAAACCTCTAACTTCTGGCAGAACCACGGt             <  1:2751002/66‑1 (MQ=18)
                                                     gTTCCTCGTCACCATCGACCCACTGAACACTGAAACCTCTAACTTCTGGCAGAACCACGGTGAGCTg       <  1:330967/67‑1 (MQ=255)
                                                        ccTCGTCACCATCGACCCACTGAACACTGAAACCTCTAACTTCTGGCAGAACCACGGTGAGCTGAACGa  >  1:681722/1‑69 (MQ=255)
                                                        ccTCGTCACCATCGACCCACTGAACACTGAAACCTCTAACTTCTGGCAGAACCACGGTGAGCTGAACGa  >  1:3042427/1‑69 (MQ=255)
                                                                 |                                                           
GTTGCCTCATTCCCGAACAAAAACAAAGTGATCGGTTGTCTGTCGAAACTGAAGTTCCTCGTCACCATCGACCCACTGAACACTGAAACCTCTAACTTCTGGCAGAACCACGGTGAGCTGAACGA  >  minE/2334229‑2334353

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