Predicted mutation
evidence seq id position mutation freq annotation gene description
RA minE 2,334,282 G→A 28.1% K577K (AAG→AAA fdoG → formate dehydrogenase‑O, large subunit

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*minE2,334,2820GA28.1% 151.8 / 55.8 131K577K (AAG→AAAfdoGformate dehydrogenase‑O, large subunit
Reads supporting (aligned to +/- strand):  ref base G (38/56);  new base A (26/11);  total (64/67)
Fisher's exact test for biased strand distribution p-value = 3.30e-03
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

Only 100 of 132 total aligned reads displayed.

AGGGCTTTAACCCTGTTGCCTCATTCCCGAACAAAAACAAAGTGATCGGTTGTCTGTCGAAACTGAAGTTCCTCGTCACCATCGACCCACTGAACACTGAAACCTCTAACTTCTGGCAGAACCACGGTGAGCTG  >  minE/2334215‑2334348
                                                                   |                                                                  
aGGGCTTTAACCCTGTTGCCTCATTCCCGAACAAAAACAAAGTGATCGGTTGTCTGTCGAAACTGAAGt                                                                   >  1:2366620/1‑69 (MQ=255)
aGGGCTTTAACCCTGTTGCCTCATTCCCGAACAAAAACAAAGTGATCGGTTGTCTGTCGAAACTGAAGt                                                                   >  1:2065589/1‑69 (MQ=255)
aGGGCTTTAACCCTGTTGCCTCATTCCCGAACAAAAACAAAGTGATCGGTTGTCTGTCGAAACTGAAGt                                                                   >  1:2418010/1‑69 (MQ=255)
aGGGCTTTAACCCTGTTGCCTCATTCCCGAACAAAAACAAAGTGATCGGTTGTCTGTCGAAACTGAAGt                                                                   >  1:1910522/1‑69 (MQ=255)
  ggCTTTAACCCTGTTGCCTCATTCCCGAACAAAAACAAAGTGATCGGTTGTCTGTCGAAACTGAAGTTc                                                                 <  1:2438209/69‑1 (MQ=255)
  ggCTTTAACCCTGTTGCCTCATTCCCGAACAAAAACAAAGTGATCGGTTGTCTGTCGAAACTGAAGTTc                                                                 <  1:1403861/69‑1 (MQ=255)
              gTTGCCTCATTCCCGAACAAAAACAAAGTGATCGGTTGTCTGTCGAAACTGAAGTTCCTCGTCACCATc                                                     >  1:2788809/1‑69 (MQ=255)
              gTTGCCTCATTCCCGAACAAAAACAAAGTGATCGGTTGTCTGTCGAAACTGAAGTTCCTCGTCACCATc                                                     >  1:214170/1‑69 (MQ=255)
              gTTGCCTCATTCCCGAACAAAAACAAAGTGATCGGTTGTCTGTCGAAACTGAAGTTCCTCGTCACCATc                                                     >  1:2203081/1‑69 (MQ=255)
              gTTGCCTCATTCCCGAACAAAAACAAAGTGATCGGTTGTCTGTCGAAACTGAAGTTCCTCGTCACCATc                                                     >  1:241970/1‑69 (MQ=255)
              gTTGCCTCATTCCCGAACAAAAACAAAGTGATCGGTTGTCTGTCGAAACTGAAGTTCCTCGTCACCATc                                                     >  1:1923357/1‑69 (MQ=255)
              gTTGCCTCATTCCCGAACAAAAACAAAGTGATCGGTTGTCTGTCGAAACTGAAGTTCCTCGTCACCATc                                                     >  1:3102112/1‑69 (MQ=255)
              gTTGCCTCATTCCCGAACAAAAACAAAGTGATCGGTTGTCTGTCGAAACTGAAGTTCCTCGTCACCATc                                                     >  1:3178869/1‑69 (MQ=255)
              gTTGCCTCATTCCCGAACAAAAACAAAGTGATCGGTTGTCTGTCGAAACTGAAGTTCCTCGTCACCATc                                                     >  1:938734/1‑69 (MQ=255)
              gTTGCCTCATTCCCGAACAAAAACAAAGTGATCGGTTGTCTGTCGAAACTGAAGTTCCTCGTCACCATc                                                     >  1:816154/1‑69 (MQ=255)
              gTTGCCTCATTCCCGAACAAAAACAAAGTGATCGGTTGTCTGTCGAAACTGAAGTTCCTCGTCACCATc                                                     >  1:447790/1‑69 (MQ=255)
               ttGCCTCATTCCCGAACAAAAACAAAGTGATCGGTTGTCTGTCGAAACTGAAGTTc                                                                 <  1:1510461/56‑1 (MQ=255)
                 gCCTCATTCCCGAACAAAAACAAAGTGATCGGTTGTCTGTCGAAACTGAAGTTCCTCGTCACCATCGa                                                   >  1:1553607/1‑68 (MQ=255)
                 gCCTCATTCCCGAACAAAAACAAAGTGATCGGTTGTCTGTCGAAACTGAAGTTCCTCGTCACCATCGa                                                   >  1:1673502/1‑68 (MQ=255)
                 gCCTCATTCCCGAACAAAAACAAAGTGATCGGTTGTCTGTCGAAACTGAAGTTCCTCGTCACCATCGa                                                   >  1:908160/1‑68 (MQ=255)
                 gCCTCATTCCCGAACAAAAACAAAGTGATCGGTTGTCTGTCGAAACTGAAGTTCCTCGTCACCATCGa                                                   >  1:2572183/1‑68 (MQ=255)
                 gCCTCATTCCCGAACAAAAACAAAGTGATCGGTTGTCTGTCGAAACTGAAGTTCCTCGTCACCATCGa                                                   >  1:1123520/1‑68 (MQ=255)
                    tCATTCCCGAACAAAAACAAAGTGATCGGTTGTCTGTCGAAACTGAAGTTCCTCGTCACCATCGACCca                                               >  1:2791965/1‑69 (MQ=255)
                    tCATTCCCGAACAAAAACAAAGTGATCGGTTGTCTGTCGAAACTGAAGTTCCTCGTCACCATCGACCca                                               >  1:2779816/1‑69 (MQ=255)
                    tCATTCCCGAACAAAAACAAAGTGATCGGTTGTCTGTCGAAACTGAAGTTCCTCGTCACCATCGACCca                                               >  1:2501252/1‑69 (MQ=255)
                    tCATTCCCGAACAAAAACAAAGTGATCGGTTGTCTGTCGAAACTGAAGTTCCTCGTCACCATCGACCca                                               >  1:138028/1‑69 (MQ=255)
                    tCATTCCCGAACAAAAACAAAGTGATCGGTTGTCTGTCGAAACTGAAGTTCCTCGTCACCATCGACCca                                               >  1:1225997/1‑69 (MQ=255)
                         cccGAACAAAAACAAAGTGATCGGTTGTCTGTCGAAACTGAAGTTCCTCGTCACCATCGACCcactgaa                                          >  1:3331391/1‑69 (MQ=255)
                         cccGAACAAAAACAAAGTGATCGGTTGTCTGTCGAAACTGAAGTTCCTCGTCACCATCGACCcactgaa                                          <  1:3152646/69‑1 (MQ=255)
                         cccGAACAAAAACAAAGTGATCGGTTGTCTGTCGAAACTGAAGTTCCTCGTCACCATCGACCcactgaa                                          <  1:2181087/69‑1 (MQ=255)
                         cccGAACAAAAACAAAGTGATCGGTTGTCTGTCGAAACTGAAGTTCCTCGTCACCATCGACCcactgaa                                          <  1:239232/69‑1 (MQ=255)
                         cccGAACAAAAACAAAGTGATCGGTTGTCTGTCGAAACTGAAGTTCCTCGTCACCATCGACCcactgaa                                          <  1:2980693/69‑1 (MQ=255)
                         cccGAACAAAAACAAAGTGATCGGTTGTCTGTCGAAACTGAAGTTCCTCGTCACCATCGACCcactgaa                                          <  1:2887270/69‑1 (MQ=255)
                         cccGAACAAAAACAAAGTGATCGGTTGTCTGTCGAAACTGAAGTTCCTCGTCACCATCGACCcactgaa                                          <  1:43196/69‑1 (MQ=255)
                         cccGAACAAAAACAAAGTGATCGGTTGTCTGTCGAAACTGAAGTTCCTCGTCACCATCGACCcactgaa                                          <  1:3370006/69‑1 (MQ=255)
                         cccGAACAAAAACAAAGTGATCGGTTGTCTGTCGAAACTGAAGTTCCTCGTCACCATCGACCcactgaa                                          <  1:2421091/69‑1 (MQ=255)
                         cccGAACAAAAACAAAGTGATCGGTTGTCTGTCGAAACTGAAGTTCCTCGTCACCATCGACCcactgaa                                          <  1:2599246/69‑1 (MQ=255)
                         cccGAACAAAAACAAAGTGATCGGTTGTCTGTCGAAACTGAAGTTCCTCGTCACCATCGACCcactgaa                                          <  1:94509/69‑1 (MQ=255)
                         cccGAACAAAAACAAAGTGATCGGTTGTCTGTCGAAACTGAAGTTCCTCGTCACCATCGACCcactgaa                                          <  1:3274371/69‑1 (MQ=255)
                         cccGAACAAAAACAAAGTGATCGGTTGTCTGTCGAAACTGAAGTTCCTCGTCACCATCGACCcactgaa                                          <  1:14130/69‑1 (MQ=255)
                         cccGAACAAAAACAAAGTGATCGGTTGTCTGTCGAAACTGAAGTTCCTCGTCACCATCGACCcactgaa                                          <  1:1105178/69‑1 (MQ=255)
                         cccGAACAAAAACAAAGTGATCGGTTGTCTGTCGAAACTGAAGTTCCTCGTCACCATCGACCcactgaa                                          <  1:118051/69‑1 (MQ=255)
                         cccGAACAAAAACAAAGTGATCGGTTGTCTGTCGAAACTGAAGTTCCTCGTCACCATCGACCcactgaa                                          <  1:1282944/69‑1 (MQ=255)
                           cGAACAAAAACAAAGTGATCGGTTGTCTGGCGAAACTGAAGTTCCTCGTCACCATCGACCcactgaa                                          <  1:3069008/67‑1 (MQ=255)
                              acaaaaacaaaGTGATCGGCTGTCTGTCGAAGCTGAAATTCCTCGTGACTAt                                                      >  1:2320308/1‑52 (MQ=25)
                              acaaaaacaaaGTGATCGGCTGTCTGTCGAAGCTGAAATTCCTCGTGACTATc                                                     >  1:456776/1‑53 (MQ=25)
                              acaaaaacaaaGTGATCGGCTGTCTGTCGAAGCTGAAATTCCTCGTGACTATc                                                     >  1:1680172/1‑53 (MQ=25)
                              acaaaaacaaaGTGATCGGCTGTCTGTCGAAGCTGAAATTCCTCGTGACTATc                                                     >  1:2856490/1‑53 (MQ=25)
                              acaaaaacaaaGTGATCGGCTGTCTGTCGAAGCTGAAATTCCTCGTGACTATc                                                     >  1:2691943/1‑53 (MQ=25)
                              acaaaaacaaaGTGATCGGCTGTCTGTCGAAGCTGAAATTCCTCGTGACTATc                                                     >  1:2689155/1‑53 (MQ=25)
                              acaaaaacaaaGTGATCGGCTGTCTGTCGAAGCTGAAATTCCTCGTGACTATc                                                     >  1:2958475/1‑53 (MQ=25)
                              acaaaaacaaaGTGATCGGCTGTCTGTCGAAGCTGAAATTCCTCGTGACTATc                                                     >  1:1343472/1‑53 (MQ=25)
                              acaaaaacaaaGTGATCGGCTGTCTGTCGAAGCTGAAATTCCTCGTGACTATc                                                     >  1:1316340/1‑53 (MQ=25)
                              acaaaaacaaaGTGATCGGCTGTCTGTCGAAGCTGAAATTCCTCGTGACTATc                                                     >  1:3394846/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:2104780/1‑53 (MQ=25)
                              acaaaaacaaaGTGATCGGCTGTCTGTCGAAGCTGAAATTCCTCGTGACTATc                                                     >  1:775449/1‑53 (MQ=25)
                              acaaaaacaaaGTGATCGGCTGTCTGTCGAAGCTGAAATTCCTCGTGACTATc                                                     >  1:2276925/1‑53 (MQ=25)
                              acaaaaacaaaGTGATCGGCTGTCTGTCGAAGCTGAAATTCCTCGTGACTATc                                                     >  1:2124770/1‑53 (MQ=25)
                              acaaaaacaaaGTGATCGGCTGTCTGTCGAAGCTGAAATTCCTCGTGACTATc                                                     >  1:1063014/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)
                                aaaaacaaaGTGATCGGTTGTCTGTCGAAACTGAAGTTCCTCGTCACCATCGACCcactgaacactgaa                                   <  1:908086/69‑1 (MQ=255)
                                aaaaacaaaGTGATCGGTTGTCTGTCGAAACTGAAGTTCCTCGTCACCATCGACCcactgaacactgaa                                   <  1:2071884/69‑1 (MQ=255)
                                aaaaacaaaGTGATCGGTTGTCTGTCGAAACTGAAGTTCCTCGTCACCATCGACCcactgaacactgaa                                   <  1:2061423/69‑1 (MQ=255)
                                aaaaacaaaGTGATCGGTTGTCTGTCGAAACTGAAGTTCCTCGTCACCATCGACCcactgaacactgaa                                   <  1:1241137/69‑1 (MQ=255)
                                aaaaacaaaGTGATCGGTTGTCTGTCGAAACTGAAGTTCCTCGTCACCATCGACCcactgaacactgaa                                   <  1:3272354/69‑1 (MQ=255)
                                aaaaacaaaGTGATCGGTTGTCTGTCGAAACTGAAGTTCCTCGTCACCATCGACCcactgaacactgaa                                   <  1:1766353/69‑1 (MQ=255)
                                aaaaacaaaGTGATCGGTTGTCTGTCGAAACTGAAGTTCCTCGTCACCATCGACCcactgaacactgaa                                   <  1:1471083/69‑1 (MQ=255)
                                aaaaacaaaGTGATCGGTTGTCTGTCGAAACTGAAGTTCCTCGTCACCATCGACCcactgaacactgaa                                   <  1:1694688/69‑1 (MQ=255)
                                aaaaacaaaGTGATCGGTTGTCTGTCGAAACTGAAGTTCCTCGTCACCATCGACCcactgaacactgaa                                   <  1:2677807/69‑1 (MQ=255)
                                aaaaacaaaGTGATCGGTTGTCTGTCGAAACTGAAGTTCCTCGTCACCATCGACCcactgaacactgaa                                   <  1:1353794/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)
                                   aacaaaGTGATCGGTTGTCTGTCGAAACTGAAGTTCCTCGTCACCATCGACCCACTGAACACTGAAAcc                                >  1:2304058/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:900187/69‑1 (MQ=255)
                                                  tgtctgtcGAAACTGAAGTTCCTCGTCACCATCGACCCACTGAACACTGAAACCTCTAACTTCTGGCAg                 <  1:213446/69‑1 (MQ=255)
                                                  tgtctgtcGAAACTGAAGTTCCTCGTCACCATCGACCCACTGAACACTGAAACCTCTAACTTCTGGCAg                 <  1:482673/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:2373874/66‑1 (MQ=18)
                                                            agcTGAAATTCCTCGTGACTATCGACCCGCTGAACACTGAAACCTCTAACTTCTGGCAGAACCACGGt        <  1:2217535/66‑1 (MQ=18)
                                                            agcTGAAATTCCTCGTGACTATCGACCCGCTGAACACTGAAACCTCTAACTTCTGGCAGAACCACGGt        <  1:1272905/66‑1 (MQ=18)
                                                            agcTGAAATTCCTCGTGACTATCGACCCGCTGAACACTGAAACCTCTAACTTCTGGCAGAACCACGGt        <  1:2751002/66‑1 (MQ=18)
                                                            agcTGAAATTCCTCGTGACTATCGACCCGCTGAACACTGAAACCTCTAACTTCTGGCAGAACCACGGt        <  1:2152581/66‑1 (MQ=18)
                                                            agcTGAAATTCCTCGTGACTATCGACCCGCTGAACACTGAAACCTCTAACTTCTGGCAGAACCACGGt        <  1:2755713/66‑1 (MQ=18)
                                                                   gTTCCTCGTCACCATCGACCCACTGAACACTGAAACCTCTAACTTCTGGCAGAACCACGGTGAGCTg  <  1:330967/67‑1 (MQ=255)
                                                                   |                                                                  
AGGGCTTTAACCCTGTTGCCTCATTCCCGAACAAAAACAAAGTGATCGGTTGTCTGTCGAAACTGAAGTTCCTCGTCACCATCGACCCACTGAACACTGAAACCTCTAACTTCTGGCAGAACCACGGTGAGCTG  >  minE/2334215‑2334348

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: