Predicted mutation
evidence seq id position mutation freq annotation gene description
RA minE 1,259,838 T→C 9.4% E362E (GAA→GAG torZ ← trimethylamine N‑oxide reductase system III, catalytic subunit

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*minE1,259,8380TC 9.4% 173.2 / 9.3 64E362E (GAA→GAGtorZtrimethylamine N‑oxide reductase system III, catalytic subunit
Reads supporting (aligned to +/- strand):  ref base T (28/30);  new base C (3/3);  total (31/33)
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.70e-01

AACCGCCGCCCGGTGTACCAATCTGCCCTAACATCGCCGCCAGCGTCACCAGCATCCAGTGTTTCTGTTCACCATATTGCTGGCGCTGAATTCCCCAGCCTGCCATCAGCATTGTGCGATTCGCCGC  >  minE/1259771‑1259897
                                                                   |                                                           
aaCCGCCGCCCGGTGTACCAATCTGCCCTAACATCGCCGCCAGCGTCACCAGCATCCAGTGTTTCTGTTCa                                                          <  1:3617065/71‑1 (MQ=255)
    gccgccCGGTGTACCAATCTGCCCTAACATCGCCGCCAGCGTCACCAGCATCCAGTGTTTCTGTTCACCa                                                       <  1:1517962/70‑1 (MQ=255)
          cGGTGTACCAATCTGCCCTAACATCGCCGCCAGCGTCACCAGCATCCAGTGTTTCTGTTCACCATATTGCt                                                <  1:3016714/71‑1 (MQ=255)
          cGGTGTACCAATCTGCCCTAACATCGCCGCCAGCGTCACCAGCATCCAGTGTTTCTGTTCACCATATTGCt                                                >  1:1710474/1‑71 (MQ=255)
            gtgtACCAATCTGCCCTAACATCGCCGCCAGCGTCACCAGCATCCAGTGTTTCTGTTCACCa                                                       >  1:1072496/1‑62 (MQ=255)
             tgtgCCGATTTGCCCCAACATTGCTGCCAGCGTGACGATCATCCAGTGTTTTTGCTCACCAAACTGTTGgc                                             <  1:2356310/71‑1 (MQ=255)
              gtACCAATCTGCCCTAACATCGCCGCCAGCGTCACCAGCATCCAGTGTTTCTGTTCACCata                                                     >  1:2754179/1‑62 (MQ=255)
                 ccAATCTGCCCTAACATCGCCGCCAGCGTCACCAGCATCCAGTGTTTCTGTTCACCATATTg                                                  >  1:4101501/1‑62 (MQ=255)
                 ccAATCTGCCCTAACATCGCCGCCAGCGTCACCAGCATCCAGTGTTTCTGTTCACCATATTGc                                                 >  1:86129/1‑63 (MQ=255)
                  cAATCTGCCCTAACATCGCCGCCAGCGTCACCAGCATCCAGTGTTTCTGTTCACCATATTGCTGGCGCTg                                         <  1:841231/70‑1 (MQ=255)
                    aTCTGCCCTAACATCGCCGCCAGCGTCACCAGCATCCAGTGTTTCTGTTCACCATATTGCTGGCGCTGAAt                                      <  1:3498550/71‑1 (MQ=255)
                    aTCTGCCCTAACATCGCCCCCAGCGTCACCAGCATCCAGTGTTTCTGTTCACCATATTGCTCGCGCTGa                                        >  1:2035281/1‑69 (MQ=255)
                      ttGCCCCAACATTGCTGCCAGCGTGACGATCATCCAGTGTTTTTGCTCACCAAACTGTTGGCGCTGCAttg                                    >  1:3651372/2‑70 (MQ=255)
                      cTGCCCTAACATCGCCGCCAGCGTCACCAGCATCCAGTGTTTCTGtt                                                            >  1:1286498/1‑47 (MQ=255)
                      cTGCCCTAACATCGCCGCCAGCGTCACCAGCATCCAGTGTTTCTGTTCACCATATTGCTGGCGCTGAATTc                                    >  1:949744/1‑71 (MQ=255)
                       tGCCCTAACATCGCCGCCAGCGTCACCAGCATCCAGTGTTTCTGTTCACCATATTGCTGgc                                             <  1:1924226/61‑1 (MQ=255)
                       tGCCCTAACATCGCCGCCAGCGTCACCAGCATCCAGTGTTTCTGTTCACCATATTGCTGGCGCTGAAtt                                     <  1:3921336/69‑1 (MQ=255)
                       tGCCCTAACATCGCCGCCAGCGTCACCAGCATCCAGTGTTTCTGTTCACCATATTGCTGGCGCTGAATTcc                                   >  1:676311/1‑71 (MQ=255)
                       tGCCCTAACATCGCCGCCAGCGTCACCAGCATCCAGTGTTTCTGTTCACCATATTGCTGGCGCTGAATTcc                                   >  1:489801/1‑71 (MQ=255)
                       tGCCCTAACATCGCCGCCAGCGTCACCAGCATCCAGTGTTTCTGTTCACCATATTGCTGGCGCTGAATTc                                    >  1:2951145/1‑70 (MQ=255)
                         cccTAACATCGCCGCCAGCGTCACCAGCATCCAGTGTTTCTGTTCACCATATTGCTGGCGCTGAATTc                                    <  1:1958911/68‑1 (MQ=255)
                          ccTAACATCGCCGCCAGCGTCACCAGCATCCAGTGTTTCTGTTCACCATATTGCTGGCGCTGAATTCCCCa                                <  1:1306150/71‑1 (MQ=255)
                          ccTAACATCGCCGCCAGCGTCACCAGCATCCAGTGTTTCTGTTCACCATATTGCTGGCGCTGAATTCCCCa                                <  1:2132639/71‑1 (MQ=255)
                           cTAACATCGCCGCCAGCGTCACCAGCATCCAGTGTTTCTGTTCACCATATTGCTGGCGCTGAATTCCCCAg                               <  1:825555/71‑1 (MQ=255)
                           cTAACATCGCCGCCAGCGTCACCAGCATCCAGTGTTTCTGTTCACCATATTGCTGGCGCTGAATTCCCCAg                               >  1:917326/1‑71 (MQ=255)
                           cTAACATCGCCGCCAGCGTCACCAGCATCCAGTGTTTCTGTTCACCATATTGCTGGCGCTGAATTCCCCAg                               <  1:620954/71‑1 (MQ=255)
                             aaCATCGCCGCCAGCGTCACCAGCATCCAGTGTTTCTGtt                                                            >  1:634628/1‑40 (MQ=255)
                             aaCATCGCCGCCAGCGTCACCAGCATCCAGTGTTTCTGtt                                                            >  1:985012/1‑40 (MQ=255)
                             aaCATCGCCGCCAGCGTCACCAGCATCCAGTGTTTCTGTTCAc                                                         >  1:3447520/1‑43 (MQ=255)
                             aaCATCGCCGCCAGCGTCACCAGCATCCAGTGTTTCTGTTCACCATATTGCTGGCGCTGaa                                       >  1:645374/1‑61 (MQ=255)
                              aCATCGCCGCCAGCGTCCCCAGCATCCAGTGTTTCTGTTCACCata                                                     >  1:789654/1‑46 (MQ=255)
                              aCATCGCCGCCAGCGTCACCAGCATCCAGTGTTTCTGTTCACCa                                                       >  1:1615056/1‑44 (MQ=255)
                              aCATCGCCGCCAGCGTCACCAGCATCCAGTGTTTCTGTTCACCATATTGCTGGCGCTGaa                                       >  1:959563/1‑60 (MQ=255)
                              aCATCGCCGCCAGCGTCACCAGCATCCAGTGTTTCTGTTCACCATATTGCTGGCGCTGAATTCCCCAg                               <  1:3682026/68‑1 (MQ=255)
                              aCATCGCCGCCAGCGTCACCAGCATCCAGTGTTTCTGTTCACCATATTGCTGGCGCTGAATTCCCCAGCCt                            >  1:3255222/1‑71 (MQ=255)
                               cATCGCCGCCAGCGTCACCAGCATCCAGTGTTTCTGt                                                             >  1:2809879/1‑37 (MQ=255)
                                aTCGCCGCCAGCGTCACCAGCATCCAGTGTTTCTGTTCACCATATTGCTGGCGCTGAATTCCCCAGCCTGc                          <  1:2724486/71‑1 (MQ=255)
                                aTCGCCGCCAGCGTCACCAGCATCCAGTGTTTCTGTTCACCATATTGCTGGCGCTGAATTCCCCAGCCTGc                          <  1:1289206/71‑1 (MQ=255)
                                 tCGCCGCCAGCGTCACCAGCATCCAGTGTTTCTGTTCACCa                                                       <  1:1761356/41‑1 (MQ=255)
                                 tCGCCGCCAGCGTCACCAGCATCCAGTGTTTCTGTTCACCa                                                       <  1:3206601/41‑1 (MQ=255)
                                 tCGCCGCCAGCGTCACCAGCATCCAGTGTTTCTGTTCACCATATTGCTGGCGCTGAATTCCCCAGCCTGc                          >  1:3702334/1‑70 (MQ=255)
                                   gctgccAGCGTGACGATCATCCAGTGTTTTTGCTCACCAAACTGTTGGCGCTGCATTCCCCAGCCTGCCAt                       >  1:840797/4‑71 (MQ=255)
                                     cgccAGCGTCACCAGCATCCAGTGTTTCTGTTCACCATATTGCTGGCGCTGAATTCCCCAGCCTGCCAt                       >  1:1375322/1‑69 (MQ=255)
                                     cgccAGCGTCACCAGCATCCAGTGTTTCTGTTCACCATATTGCTGGCGCTGAATTCCCCAGCCTGCCAt                       >  1:3751334/1‑69 (MQ=255)
                                         aGCGTGACGATCATCCAGTGTTTTTGCTCACCAAACTGTTGGCGCTg                                         <  1:3275553/47‑1 (MQ=12)
                                          gCGTCACCAGCATCCAGTGTTTCTGTTCACCATATTGCTGGCGCTGAATTCCCCACCCTGCCATCAGCAt                 <  1:245410/70‑1 (MQ=255)
                                           cGTCACCAGCATCCAGTGTTTCTGTTCACCATATTGCTGGCGCTGAATTCCCCAGCCTGCCATCAGCATtg               <  1:2839781/71‑1 (MQ=255)
                                            gTCACCAGCATCCAGTGTTTCTGTTCACCATATTGCTGGCGCTGAATTCCCCAGCCTGCCATCAGc                   >  1:2158958/1‑66 (MQ=255)
                                              cACCAGCATCCAGTGTTTCTGTTCACCATATTGCTGGCGCTGAATTccc                                  <  1:1958147/49‑1 (MQ=255)
                                               aCCAGCATCCAGTGTTTCTGTTCACCATATTGCTGGCGCTGa                                        <  1:1383694/42‑1 (MQ=255)
                                               aCCAGCATCCAGTGTTTCTGTTCACCATATTGCTGGCGCTGAATTCCCCAGCCTGCCATCAGCATTGTGCg           <  1:1953343/71‑1 (MQ=255)
                                                 gatcATCCAGTGTTTTTGCTCACCAAACTGTTGGCGCTGCATTCCCCAGCCTGCCATCAGCa                  >  1:3837548/4‑62 (MQ=255)
                                                  aGCATCCAGTGTTTCTGTTCACCATATTGCTGGCGCTGAATTCCCCAGCCTGCCATCAGCATTGTGCGAtt        <  1:1078791/71‑1 (MQ=255)
                                                   gCATCCAGTGTTTCTGTTCACCATATTGCTGGCGCTGAATTCCCCAGCCTGCCATCAGCATTGTGCGATTc       >  1:1605191/1‑71 (MQ=255)
                                                   gCATCCAGTGTTTCTGTTCACCATATTGCTGGCGCTGAATTCCCCAGCCTGCCATCAGCATTGTGCGATTc       <  1:3837314/71‑1 (MQ=255)
                                                      tCCAGTGTTTCTGTTCACCATATTGCTGGCGCTGAATTCCCCAGc                              >  1:213007/1‑45 (MQ=255)
                                                      tCCAGTGTTTCTGTTCACCATATTGCTGGCGCTGAATTCCCCAGc                              >  1:3160937/1‑45 (MQ=255)
                                                        cAGTGTTTCTGTTCACCATATTGCTGGCGCTGAATTcccc                                 >  1:3446708/1‑40 (MQ=255)
                                                        cAGTGTTTCTGTTCACCATATTGCTGGCGCTGAATTcccc                                 >  1:961537/1‑40 (MQ=255)
                                                        cAGTGTTTCTGTTCACCATATTGCTGGCGCTGAATTCCCCAGCCTGCCATCAGCATTGTGCGATTcgccgc  <  1:1865392/71‑1 (MQ=255)
                                                            gtTTCTGTTCACCATATTGCTGGCGCTGAATTCCCCAGCCTGCCATCAGCAt                 >  1:106811/1‑52 (MQ=255)
                                                             ttttTGCTCACCAAACTGTTGGCGCTGCATTCCCCAGCCTGCCATCa                     <  1:1316465/47‑1 (MQ=19)
                                                             tttCTGTTCACGATATTGCTGGCGCTGAATTCCCCAGCCTg                           <  1:3686287/41‑1 (MQ=255)
                                                             tttCTGTTCACCATATTGCTGGCGCTGAATTCCCCAGCCTg                           <  1:1820502/41‑1 (MQ=255)
                                                             tttCTGTTCACCATATTGCTGGCGCTGAATTCCCCAGCCTGCCATCAGCa                  <  1:808004/50‑1 (MQ=255)
                                                             tttCTGTTCACCATATTGCTGGCGCTGAATTCCCCAGCCTGCCATCAGCATTGTGCGATTc       <  1:1496848/61‑1 (MQ=255)
                                                               tCTGTTCACCATATTGCTGGCGCTGAATTCCCCAGCCTGcc                         >  1:563965/1‑41 (MQ=255)
                                                                 tGTTCACCATATTGCTGGCGCTGAATTCCCCAGCCTGc                          <  1:3151430/38‑1 (MQ=255)
                                                                   |                                                           
AACCGCCGCCCGGTGTACCAATCTGCCCTAACATCGCCGCCAGCGTCACCAGCATCCAGTGTTTCTGTTCACCATATTGCTGGCGCTGAATTCCCCAGCCTGCCATCAGCATTGTGCGATTCGCCGC  >  minE/1259771‑1259897

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