Predicted mutation
evidence seq id position mutation freq annotation gene description
RA minE 1,259,822 G→T 7.4% L368M (CTG→ATG) ‡ 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,8220GT 7.4% 222.0 / 2.8 67L368M (CTG→ATG) ‡torZtrimethylamine N‑oxide reductase system III, catalytic subunit
Reads supporting (aligned to +/- strand):  ref base G (34/28);  new base T (2/3);  total (36/31)
Fisher's exact test for biased strand distribution p-value = 6.56e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.79e-01

ATGATAAGAGAAACCAAAACCGCCGCCCGGTGTACCAATCTGCCCTAACATCGCCGCCAGCGTCACCAGCATCCAGTGTTTCTGTTCACCATATTGCTGGCGCTGAATTCCCCAGCCTGCCATCAGCATTGTGCGATTC  >  minE/1259754‑1259892
                                                                    |                                                                      
aTGATAAGAGAAACCAAAACCGCCGCCCGGTGTACCAATCTGCCCTAACATCGCCGCCAGCGTCACCAGCa                                                                      >  1:2585453/1‑71 (MQ=255)
    tAAGAGAAACCAAAACCGCCGCCCGGTGTACCAATCTGCCCTAACATCGCCGCCAGCGTCACCAGCATc                                                                    >  1:3670793/1‑69 (MQ=255)
            aCCAAAACCGCCGCCCGGTGTACCAATCTGCCCTAACATCGCCGCCAGCGTCACCAGCATCCAGTGTTTCt                                                          >  1:1320390/1‑71 (MQ=255)
              cAAAACCGCCGCCGGGTGTGCCGATTTGCCCCAACATTGCTGCCAGCGTGACGATCATCCAGTGTTTTTgc                                                        <  1:2528450/71‑2 (MQ=255)
                 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)
                            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:2043001/1‑42 (MQ=255)
                                        tGCCCTAACATCGCCGCCAGCGTCACCAGCATCCAGTGTTTc                                                           >  1:3817044/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:917326/1‑71 (MQ=255)
                                            cTAACATCGCCGCCAGCGTCACCAGCATCCAGTGTTTCTGTTCACCATATTGCTGGCGCTGAATTCCCCAg                          <  1:620954/71‑1 (MQ=255)
                                            cTAACATCGCCGCCAGCGTCACCAGCATCCAGTGTTTCTGTTCACCATATTGCTGGCGCTGAATTCCCCAg                          <  1:825555/71‑1 (MQ=255)
                                              aaCATCGCCTCCAGCGTCACCAGCATCCAGTGTTTCTg                                                         <  1:1663325/38‑1 (MQ=255)
                                              aaCATCGCCGCCAGCGTCACCAGCATCCAGTGTTTCTg                                                         <  1:3236500/38‑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:3206601/41‑1 (MQ=255)
                                                  tCGCCGCCAGCGTCACCAGCATCCAGTGTTTCTGTTCACCa                                                  <  1:1761356/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)
                                                                   aGCATCCAGTGTTTCTGTTCACCATATTGCTGGCGCTGAATTCCCCAGCCTGCCATCAGCATTGTGCGAtt   <  1:1078791/71‑1 (MQ=255)
                                                                    gCATCCAGTGTTTCTGTTCACCATATTGCTGGCGCTGAATTCCCCAGCCTGCCATCAGCATTGTGCGATTc  >  1:1605191/1‑71 (MQ=255)
                                                                    gCATCCAGTGTTTCTGTTCACCATATTGCTGGCGCTGAATTCCCCAGCCTGCCATCAGCATTGTGCGATTc  <  1:3837314/71‑1 (MQ=255)
                                                                    |                                                                      
ATGATAAGAGAAACCAAAACCGCCGCCCGGTGTACCAATCTGCCCTAACATCGCCGCCAGCGTCACCAGCATCCAGTGTTTCTGTTCACCATATTGCTGGCGCTGAATTCCCCAGCCTGCCATCAGCATTGTGCGATTC  >  minE/1259754‑1259892

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