Predicted mutation
evidence seq id position mutation freq annotation gene description
RA minE 1,259,799 T→C 10.2% L375L (TTA→TTG 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,7990TC10.2% 161.2 / 9.5 59L375L (TTA→TTGtorZtrimethylamine N‑oxide reductase system III, catalytic subunit
Reads supporting (aligned to +/- strand):  ref base T (27/26);  new base C (2/4);  total (29/30)
Fisher's exact test for biased strand distribution p-value = 6.71e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 6.29e-01

GCGTCGGGTTGCCACCGTTGGAGTAATGATAAGAGAAACCAAAACCGCCGCCCGGTGTACCAATCTGCCCTAACATCGCCGCCAGCGTCACCAGCATCCAGTGTTTCTGTTCACCATATTGCTGGCGCTGAATTCCCCAG  >  minE/1259729‑1259868
                                                                      |                                                                     
gcgTCGGGTTGCCACCGTTGGAGTAATGATAAGAGAAACCAAAACCGCCGCCCGGTGTACCAATCTGCCCt                                                                       <  1:2445320/71‑1 (MQ=255)
   tggggTTACCACCATTGGCAAAATGGTAAGAAAGACCAAAACCGCCGCCGGGTGTGCCGATTTGCCCCaa                                                                     <  1:3521322/68‑1 (MQ=255)
    cGGGTTGCCACCGTTGGAGTAATGATAAGAGAAACCAAAACCGCCGCCCGGTGTACCAATCTGCCCTAACa                                                                   >  1:2400682/1‑71 (MQ=255)
          gCCACCGTTGGAGTAATGATAAGAGAAACCAAAACCGCCGCCCGGTGTACCAATCTGCCCTAACATcgcc                                                              <  1:3819074/70‑1 (MQ=255)
          gCCACCGTTGGAGTAATGATAAGAGAAACCAAAACCGCCGCCCGGTGTACCAATCTGCCCTAACATcgcc                                                              >  1:2857677/1‑70 (MQ=255)
             aCCGTTGGAGTAATGATAAGAGAAACCAAAACCGCCGCCCTGTGTACCAATCTGCCCTAACATcgccgc                                                            >  1:4054426/1‑69 (MQ=255)
             aCCGTTGGAGTAATGATAAGAGAAACCAAAACCGCCGCCCGGTGTACCAATCTGCCCTAACATCGCCGCCa                                                          <  1:2978550/71‑1 (MQ=255)
             aCCGTTGGAGTAATGATAAGAGAAACCAAAACCGCCGCCCGGTGTACCAATCTGCCCTAACATCGCCGCCa                                                          <  1:177382/71‑1 (MQ=255)
               cGTTGGAGTAATGATAAGAGAAACCAAAACCGCCGCCCGGTGTACCAATCTGCCCTAACATCGCCGCCAGc                                                        >  1:2220293/1‑71 (MQ=255)
                gTTGGAGTAATGATAAGAGAAACCAAAACCGCCGCCCGGTGTACCAATCTGCCCTAAc                                                                    >  1:494220/1‑58 (MQ=255)
                 ttGGAGTAATGATAAGAGAAACCAAAACCGCCGCCCGGTGTACCAATCTGCCCTAACATCGCCGCCAGCGt                                                      <  1:1597373/71‑1 (MQ=255)
                    gAGTAATTATAAGAGAAACCAAAACCGCCGCCCGGTGTACCAATCTGCCCTAACATCGCCGCCAGCGTCa                                                    <  1:475474/70‑1 (MQ=255)
                        aaTGGTAAGAAAGACCAAAACCGCCGCCGGGTGTGCCGATTTGCCCCAACATTGCTGCCAGCGTGAcgat                                                >  1:3313564/1‑67 (MQ=255)
                         aTGATAAGAGAAACCAAAACCGCCGCCCGGTGTACCAATCTGCCCTAACATCGCCGCCAGCGTCACCAGCa                                              >  1:2585453/1‑71 (MQ=255)
                             tAAGAGAAACCAAAACCGCCGCCCGGTGTACCAATCTGCCCTa                                                                      <  1:1611531/43‑1 (MQ=255)
                             tAAGAGAAACCAAAACCGCCGCCCGGTGTACCAATCTGCCCTAACATCGCCGCCAGCGTCACCAGCATc                                            >  1:3670793/1‑69 (MQ=255)
                             tAAGAAAGACCAAAACCGCCGCCGGGTGTGCCGATTTGCCCCAACATTGCTGCCAGCGTGAcga                                                 <  1:1886895/64‑3 (MQ=255)
                                    aaCCAAAACCGCCGCCCGGTGTACCAATCTGCCCTAACATcgccgc                                                            <  1:2095070/46‑1 (MQ=255)
                                     aCCAAAACCGCCGCCCGGTGTACCAATCTGCCCTAACATCGCCGCCAGCGTCACCAGCATCCAGTGTTTCt                                  >  1:1320390/1‑71 (MQ=255)
                                      ccAAAACCGCCGCCCGGTGTACCAATCTGCCCTAACATCgccgcc                                                           <  1:1370525/45‑1 (MQ=255)
                                      ccAAAACCGCCGCCCGGTGTACCAATCTGCCCTAACATCgccgcc                                                           <  1:638572/45‑1 (MQ=255)
                                      ccAAAACCGCCGCCCGGTGTACCAATCTGCCCTAACATCGCCGCCAGCGt                                                      <  1:2427535/50‑1 (MQ=255)
                                       cAAAACCGCCGCCGGGTGTGCCGATTTGCCCCAACATTGCTGCCAGCGTGACGATCATCCAGTGTTTTTgc                                <  1:2528450/71‑2 (MQ=255)
                                       aaaaaCCGCCGCCCGGTGTACCAATCTGGCCTAACATCGCCGCCAGc                                                        <  1:3523507/46‑1 (MQ=255)
                                        aaaaCCGCCGCCCGGTGTACCAATCTGCCCTAACATCgccgcc                                                           >  1:937192/1‑43 (MQ=255)
                                          aaCCGCCGCCCGGTGTACCAATCTGCCCTAACATCGCCGCCAGCGTCACCAGCATCCAGTGTTTCTGTTCa                             <  1:3617065/71‑1 (MQ=255)
                                              gccgccCGGTGTACCAATCTGCCCTAACATCGCCGCCAGCGTCACCAGCATCCAGTGTTTCTGTTCACCa                          <  1:1517962/70‑1 (MQ=255)
                                                    cGGTGTACCAATCTGCCCTAACATCGCCGCCAGCGTCACCAGCATCCAGTGTTTCTGTTCACCATATTGCt                   >  1:1710474/1‑71 (MQ=255)
                                                    cGGTGTACCAATCTGCCCTAACATCGCCGCCAGCGTCACCAGCATCCAGTGTTTCTGTTCACCATATTGCt                   <  1:3016714/71‑1 (MQ=255)
                                                     ggTGTACCAATCTGCCCTAACATCGCCGCCAGCGTCACCAGCATCCAGTGTTTCTg                                 <  1:3907139/56‑1 (MQ=255)
                                                     ggTGTACCAATCTGCCCTAACATCGCCGCCAGCGTCACCAGCATCCAGTGTTTCTg                                 <  1:1351023/56‑1 (MQ=255)
                                                      gtgtACCAATCTGCCCTAACATCGCCGCCAGCGTCACCAGCATCCAGTGTTTCTGTTCACCa                          >  1:1072496/1‑62 (MQ=255)
                                                       tgtgCCGATTTGCCCCAACATTGCTGCCAGCGTGACGATCATCCAGTGTTTTTGCTCACCAAACTGTTGgc                <  1:2356310/71‑1 (MQ=255)
                                                        gtACCAATCTGCCCTAACATCGCCGCCAGCGTCACCAGCATCCAGTGTTTCTg                                 >  1:2074946/1‑53 (MQ=255)
                                                        gtACCAATCTGCCCTAACATCGCCGCCAGCGTCACCAGCATCCAGTGTTTCTGTTCACCata                        >  1:2754179/1‑62 (MQ=255)
                                                           ccAATCTGCCCTAACATCGCCGCCAGCGTCACCAGCATCCAGTGttt                                    >  1:1906000/1‑47 (MQ=255)
                                                           ccAATCTGCCCTAACATCGCCGCCAGCGTCACCAGCATCCAGTGTTTCTGTTCACCATATTg                     >  1:4101501/1‑62 (MQ=255)
                                                           ccAATCTGCCCTAACATCGCCGCCAGCGTCACCAGCATCCAGTGTTTCTGTTCACCATATTGc                    >  1:86129/1‑63 (MQ=255)
                                                            cAATCTGCCCTAACATCGCCGCCAGCGTCACCAGCATCCAGTGTTTCTGTTCACCATATTGCTGGCGCTg            <  1:841231/70‑1 (MQ=255)
                                                             aaTCTGCCCTAACATCGCCGCCAGCGTCACCAGCATCCAGTGtt                                     <  1:3706085/44‑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)
                                                                 tGCCCTAACATCGCCGCCAGCGTCACCAGCATCCAGTGttt                                    >  1:611931/1‑41 (MQ=255)
                                                                 tGCCCTAACATCGCCGCCAGCGTCACCAGCATCCAGTGttt                                    >  1:1722870/1‑41 (MQ=255)
                                                                 tGCCCTAACATCGCCGCCAGCGTCACCAGCATCCAGTGTTTc                                   >  1:3817044/1‑42 (MQ=255)
                                                                 tGCCCTAACATCGCCGCCAGCGTCACCAGCATCCAGTGTTTc                                   >  1:2043001/1‑42 (MQ=255)
                                                                 tGCCCTAACATCGCCGCCAGCGTCACCAGCATCCAGTGTTTCTGTTCACCATATTGCTGgc                <  1:1924226/61‑1 (MQ=255)
                                                                 tGCCCTAACATCGCCGCCAGCGTCACCAGCATCCAGTGTTTCTGTTCACCATATTGCTGGCGCTGAAtt        <  1:3921336/69‑1 (MQ=255)
                                                                 tGCCCTAACATCGCCGCCAGCGTCACCAGCATCCAGTGTTTCTGTTCACCATATTGCTGGCGCTGAATTcc      >  1:489801/1‑71 (MQ=255)
                                                                 tGCCCTAACATCGCCGCCAGCGTCACCAGCATCCAGTGTTTCTGTTCACCATATTGCTGGCGCTGAATTcc      >  1:676311/1‑71 (MQ=255)
                                                                 tGCCCTAACATCGCCGCCAGCGTCACCAGCATCCAGTGTTTCTGTTCACCATATTGCTGGCGCTGAATTc       >  1:2951145/1‑70 (MQ=255)
                                                                   cccTAACATCGCCGCCAGCGTCACCAGCATCCAGTGTTTCTGTTCACCATATTGCTGGCGCTGAATTc       <  1:1958911/68‑1 (MQ=255)
                                                                    ccTAACATCGCCGCCAGCGTCACCAGCATCCAGTGTTTCTGTTCACCATATTGCTGGCGCTGAATTCCCCa   <  1:2132639/71‑1 (MQ=255)
                                                                    ccTAACATCGCCGCCAGCGTCACCAGCATCCAGTGTTTCTGTTCACCATATTGCTGGCGCTGAATTCCCCa   <  1:1306150/71‑1 (MQ=255)
                                                                     cTAACATCGCCGCCAGCGTCACCAGCATCCAGTGTTTCTGTTCACCATATTGCTGGCGCTGAATTCCCCAg  <  1:620954/71‑1 (MQ=255)
                                                                     cTAACATCGCCGCCAGCGTCACCAGCATCCAGTGTTTCTGTTCACCATATTGCTGGCGCTGAATTCCCCAg  <  1:825555/71‑1 (MQ=255)
                                                                     cTAACATCGCCGCCAGCGTCACCAGCATCCAGTGTTTCTGTTCACCATATTGCTGGCGCTGAATTCCCCAg  >  1:917326/1‑71 (MQ=255)
                                                                      |                                                                     
GCGTCGGGTTGCCACCGTTGGAGTAATGATAAGAGAAACCAAAACCGCCGCCCGGTGTACCAATCTGCCCTAACATCGCCGCCAGCGTCACCAGCATCCAGTGTTTCTGTTCACCATATTGCTGGCGCTGAATTCCCCAG  >  minE/1259729‑1259868

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: