Predicted mutation
evidence seq id position mutation freq annotation gene description
RA minE 1,259,847 T→C 21.4% Q359Q (CAA→CAG 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,8470TC21.4% 136.6 / 34.4 70Q359Q (CAA→CAGtorZtrimethylamine N‑oxide reductase system III, catalytic subunit
Reads supporting (aligned to +/- strand):  ref base T (24/31);  new base C (10/5);  total (34/36)
Fisher's exact test for biased strand distribution p-value = 1.49e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 6.58e-01

CCCGGTGTACCAATCTGCCCTAACATCGCCGCCAGCGTCACCAGCATCCAGTGTTTCTGTTCACCATATTGCTGGCGCTGAATTCCCCAGCCTGCCATCAGCATTGTGCGATTCGCCGCCATCAGTTCAGCGAGTTT  >  minE/1259779‑1259915
                                                                    |                                                                    
cccGGTGTACCAATCTGCCCTAACATCGCCGCCAGCGTCACCAGCATCCAGTGTTTCTGTTCACCATAtt                                                                     <  1:507695/70‑1 (MQ=255)
  cGGTGTACCAATCTGCCCTAACATCGCCGCCAGCGTCACCAGCATCCAGTGTTTCTGTTCACCATATTGc                                                                   >  1:589207/1‑70 (MQ=255)
  cGGTGTACCAATCTGCCCTAACATCGCCGCCAGCGTCACCAGCATCCAGTGTTTCTGTTCACCATATTGCt                                                                  >  1:1012164/1‑71 (MQ=255)
  cGGTGTACCAATCTGCCCTAACATCGCCGCCAGCGTCACCAGCATCCAGTGTTTCTGTTCACCATATTGCt                                                                  >  1:1824658/1‑71 (MQ=255)
   ggTGTACCAATCTGCCCTAACATCGCCGCCAGCGTCACCAGCATCCAGTGTTTCTGTTCACCATATTGCTg                                                                 >  1:1156444/1‑71 (MQ=255)
   ggTGTACCAATCTGCCCTAACATCGCCGCCAGCGTCACCAGCATCCAGTGTTTCTGTTCACCATATTGCTg                                                                 <  1:1504740/71‑1 (MQ=255)
      gtACCAATCTGCCCTAACATCGCCGCCAGCGTCACCAGCATCCAGTGTTTCTGTTCACCATATTGCTGgcg                                                              >  1:3950268/1‑71 (MQ=255)
       tACCAATCTGCCCTAACATCGCCTCCAGCGTCACCAGCATCCAGTGTTTCTGTTCACCATATTGCTGgcg                                                              <  1:1955581/70‑1 (MQ=255)
            atttGCCCCAACATTGCTGCCAGCGTGACGATCATCCAGTGTTTTTGCTCACCAAACTGTTGGCGCTGCAt                                                        >  1:1223111/4‑71 (MQ=255)
               tGCCCTAACATCGCCGCCAGCGTCACCAGCATCCAGTGTTTCTGTTCACCATATTGCTGGCGCTg                                                           >  1:3416274/1‑65 (MQ=255)
               tGCCCTAACATCGCCGCCAGCGTCACCAGCATCCAGTGTTTCTGTTCACCATATTGCTGGCGCTg                                                           >  1:3121549/1‑65 (MQ=255)
                 cccTAACATCGCCGCCAGCGTCACCAGCATCCAGTGTTTCTGTTCACCATAtt                                                                     <  1:3882319/53‑1 (MQ=255)
                  ccTAACATCGCCGCCAGCGTCACCAGCATCCAGTGTTTCTGTTCACCATATTGCTGacgc                                                             >  1:2492214/1‑60 (MQ=255)
                  ccTAACATCGCCGCCAGCGTCACCAGCATCCAGTGTTTCTGTTCACCATATTGCTGGCGCTGAATTCCCCa                                                  <  1:843594/71‑1 (MQ=255)
                  ccTAACATCGCCGCCAGCGTCACCAGCATCCAGTGTTTCTGTTCACCATATTGCTGGCGCTGAATTCCCCa                                                  <  1:3514365/71‑1 (MQ=255)
                    tAACATCGCCGCCAGCGTCACCAGCATCCAGTGTTTCTGTTCACCATATTg                                                                    <  1:1012511/51‑1 (MQ=255)
                    caaCATTGCTGCCAGCGTGACGATCATCCAGTGTTTTTGCTCACCAAACTGTTGGCGCTg                                                           <  1:2524134/59‑1 (MQ=255)
                    caaCATTGCTGCCAGCGTGACGATCATCCAGTGTTTTTGCTCACCAAACTGTTGGCGCTGCATTCCCCAGc                                                >  1:3033525/2‑71 (MQ=255)
                      aCATTGCTGCCAGCGTGACGATCATCCAGTGTTTTTGCTCACCAAACTGTTGGCGCTgca                                                         >  1:3921167/1‑58 (MQ=255)
                      aCATTGCTGCCAGCGTGACGATCATCCAGTGTTTTTGCTCACCAAACTGTTGGCGCTGCATTCCCCAGCCt                                              >  1:2520160/1‑71 (MQ=255)
                      aCATCGCCGCCAGCGTCACCAGCATCCAGTGTTTCTGTTCACCATATTGCTGGCGCTg                                                           >  1:3673769/1‑58 (MQ=255)
                       cATCGCCGCCAGCGTCACCAGCATCCAGTGTTTCTGTTCACCATATTGCTGGCGCTGAATTCCCCAGCCTg                                             <  1:2655954/71‑1 (MQ=255)
                           gctgccAGCGTGACGATCATCCAGTGTTTTTGCTCACCAAACTGTTGGCGCTGCATTCCCCAGCCTGCCa                                          <  1:3094874/67‑1 (MQ=255)
                           gccgccAGCGTCACCAGCATCCAGTGTTTCTGTTCACCATATTGc                                                                   >  1:1904084/1‑45 (MQ=255)
                           gccgccAGCGTCACCAGCATCCAGTGTTTCTGTTCACCATATTGc                                                                   >  1:3002658/1‑45 (MQ=255)
                           gccgccAGCGTCACCAGCATCCAGTGTTTCTGTTCACCATATTGCTGGCGCTGAATTCCCCAGCCTGCCAt                                         >  1:2216343/1‑71 (MQ=255)
                           gccgccAGCGTCACCAGCATCCAGTGTTTCTGTTCACCATATTGCTGGCGCTGAATTCCCCAGCCTGCCAt                                         >  1:2221078/1‑71 (MQ=255)
                            ctgccAGCGTGACGATCATCCAGTGTTTTTGCTCACCAAACTGTTGGCGCTGCATTCCCCAGCCTGCCAt                                         >  1:3379490/3‑70 (MQ=255)
                            ctgccAGCGTGACGATCATCCAGTGTTTTTGCTCACCAAACTGTTGGCGCTGCATTCCCCAGCCTGCCAt                                         <  1:2241161/68‑1 (MQ=255)
                              gccAGCGTGACGATCATCCAGTGTTTTTGCTCACCAAACTGTTGGCGCTGCATTCCCCAGCCTGCCATCAg                                      >  1:1137738/1‑71 (MQ=255)
                              gccAGCGTCACCAGCATCCAGTGTTTCTGTTCACCATATTGCTGGCGCTGAATTCCCCAGCCTGCCATCAg                                      >  1:2827844/1‑71 (MQ=255)
                               ccAGCCTCACCAGCATCCAGTGTTTCTGTTCACCATATTGCTGGCGCTGAATTCCCCAGCCTg                                             <  1:1938514/63‑1 (MQ=255)
                                cAGCGTCACCAGCATCCAGTGTTTCTGTTCACCATATTGCTGGCGCTGAATTCCCCAGCCTGCCATCAGCa                                    <  1:235307/71‑1 (MQ=255)
                                cAGCGTCACCAGCATCCAGTGTTTCTGTTCACCATATTGCTGGCGCTGAATTCCCCAGCCTGCCATCAGCa                                    <  1:462589/71‑1 (MQ=255)
                                cAGCGTCACCAGCATCCAGTGTTTCTGTTCACCATATTGCTGGCGCTGAATTCCCCAGCCTGCCATCAGCa                                    <  1:3841247/71‑1 (MQ=255)
                                cAGCGTCACCAGCATCCAGTGTTTCTGTTCACCATATTGCTGGCGCTGAATTCCCCAGCCTGCCATCAGCa                                    <  1:1910719/71‑1 (MQ=255)
                                cAGCGTCACCAGCATCCAGTGTTTCTGTTCACCATATTGCTGGCGCTGAATTCCCCAGCCTGCCATCAGCa                                    >  1:3621450/1‑71 (MQ=255)
                                cAGCGTCACCAGCATCCAGTGTTTCTGTTCACCATATTGCTGGCGCTGAATTCCCCAGCCTGCCATCAGCa                                    <  1:1507521/71‑1 (MQ=255)
                                 aGCGTCACCAGCATCCAGTGTTTCTGTTCACCATATTGCTGGCGCTGAATTCCCCAGCCTGCCATCAGCAt                                   >  1:792378/1‑71 (MQ=255)
                                  gCGTGACGATCATCCAGTGTTTTTGCTCACCAAACTGTTGGCGCTGCATTCCCCAGCCTGCCATCAGCAtg                                  <  1:480872/71‑2 (MQ=255)
                                  gCGTGACGATCATCCAGTGTTTTTGCTCACCAAACTGTTGGCGCTGCATTCCCCAGCCTGCCATCAGCAtg                                  <  1:1881162/71‑2 (MQ=255)
                                   cGTCACCAGCATCCAGTGTTTCTGTTCACCATATTGCTGGCGCTGAATTCCCCAGCCTGCCATCAGCATtg                                 >  1:1084140/1‑71 (MQ=255)
                                    gTCACCAGCATCCAGTGTTTCTGTTCACCATATTGCTGGCGCTGAATTCCCCAGCCTGCCATCAGCATtgt                                <  1:415203/71‑1 (MQ=255)
                                    gTCACCAGCATCCAGTGTTTCTGTTCACCATATTGCTGGCGCTGAATTCCCCAGCCTGCCATCAGCATtg                                 <  1:371696/70‑1 (MQ=255)
                                      cACCAGCATCCAGTGTTTCTGTTCACCATATTGCTGGCGCTGAATTCCCCAGCCTGCCATCAGCATTGTGc                              >  1:2726561/1‑71 (MQ=255)
                                          aGCATCCAGTGTTTCTGTTCACCATATTGCTGGCGCTGAATTCCCCAGCCTGCCATCAGCATTGTGCGAtt                          >  1:3160833/1‑71 (MQ=255)
                                           gCATCCAGTGTTTCTGTTCACCATATTGCTGGCGCTGAATTCCCCAGCCTGCCATCAGCATTGTGCGATTc                         >  1:3191134/1‑71 (MQ=255)
                                           gCATCCAGTGTTTCTGTTCACCATATTGCTGGCGCTGAATTCCCCAGCCTGCCATCAGCATTGTGCGATTc                         <  1:2928067/71‑1 (MQ=255)
                                             aTCCAGTGTTTCTGTTCACCATATTGCTGGCGCTGAATTCCCCAGCCTGCCATCAGCATTGTGCGATTcgc                       <  1:802895/71‑1 (MQ=255)
                                              tCCAGTGTTTCTGTTCACCATATTGCTGGCGCTGAATTCCCCAGCCTGCCATCAGCATTGTGCGATTcgcc                      <  1:715872/71‑1 (MQ=255)
                                                cAGTGTTTTTGCTCACCAAACTGTTGGCGCTGCATTCCCCAGCCTGCCATc                                        >  1:1671019/1‑51 (MQ=255)
                                                cAGTGTTTTTGCTCACCAAACTGTTGGCGCTGCATTCCCCAGCCTGCCATc                                        >  1:3445555/1‑51 (MQ=255)
                                                cAGTGTTTTTGCTCACCAAACTGTTGGCGCTGCATTCCCCAGCCTGCCATc                                        >  1:1226739/1‑51 (MQ=255)
                                                cAGTGTTTTTGCTCACCAAACTGTTGGCGCTGCATTCCCCAGCCTGCCATc                                        >  1:1240331/1‑51 (MQ=255)
                                                  gtgtTTCTGTTCACCATATTGCTGGCGCTGAATTcccc                                                   <  1:352911/38‑1 (MQ=255)
                                                  gtgtTTCTGTTCACCATATTGCTGGCGCTGAATTCCCCAGCCTGCCATCAGCATTGTGCGAt                           <  1:1949623/62‑1 (MQ=255)
                                                   tgtTTCTGTTCACCATATTGCTGGCGCTGAATTCCCCAGCCTGCCATCAGCATTGTGCGATTc                         <  1:3684044/63‑1 (MQ=255)
                                                   tgtTTCTGTTCACCATATTGCTGGCGCTGAATTCCCCAGCCTGCCATCAGCATTGTGCGATTCGCCGCCAt                 <  1:1126505/71‑1 (MQ=255)
                                                   tgtTTCTGTTCACCATATTGCTGGCGCTGAATTCCCCAGCCTGCCATCAGCATTGTGCGATTCGCCGCCAt                 <  1:1407074/71‑1 (MQ=255)
                                                    gtTTCTGTTCACCATATTGCTGGCGCTGAATTCCCCAGCCTGGCATCAGCATTGTGCGATTCGCCGCCAt                 >  1:34580/1‑70 (MQ=255)
                                                    gtTTCTGTTCACCATATTGCTGGCGCTGAATTCCCCAGCCTGCCATCAGCATTGTGc                              <  1:494234/57‑1 (MQ=255)
                                                     tttCTGTTCACCATATTGCTGGCGCTGAATTCCCCAGCCTGCCATCAGCATTGTGCGATTcgccgc                    <  1:295695/66‑1 (MQ=255)
                                                     tttCTGTTCACCATATTGCTGGCGCTGAATTCCCCAGCCTGCCATCAGCATTGTGCGATTCgccgcc                   <  1:3369877/67‑1 (MQ=255)
                                                     tttCTGTTCACCATATTGCTGGCGCTGAATTCCCCAGCCTGCCATCAGCATTGTGCGATTCGCCGCCATc                <  1:1356276/70‑1 (MQ=255)
                                                     tttCTGTTCACCATATTGCTGGCGCTGAATTCCCCAGCCTGCCATCAGCATTGTGCGATTCGCCGCCATc                <  1:3370159/70‑1 (MQ=255)
                                                       tCTGTTCACCATATTGCTGGCGCTGAATTCCCCAGCCTGCCATCAGCATTGTGCGATTCGCCGCCATCAg              <  1:57242/70‑1 (MQ=255)
                                                        cTGTTCACCATATTGCTGGCGCTGAATTCCCCAGCCTGCCATCAGCATTGTGCGATTCGCCGCCATCAGtt            >  1:310941/1‑71 (MQ=255)
                                                        cTGTTCACCATATTGCTGGCGCTGAATTCCCCAGCCTGCCATCAGCATTGTGCGATTCGCCGCCATCAGtt            >  1:3924097/1‑71 (MQ=255)
                                                             cACCATATTGCTGGCGCTGAATTCCCCAGCCTGCCATCAGCATTGTGCGATTCGCCGCCATCAGTTCAGc        <  1:3130496/70‑1 (MQ=255)
                                                              aCCATATTGCTGGCGCTGAATTCCCCAGCCTGCCATCAGCATTGTGCGATTCGCCGCCATCAGTTCAGCGa      <  1:2943778/71‑1 (MQ=255)
                                                               ccATATTGCTGGCGCTGAATTCCCCAGCCTGCCATCAGCATTGTGCGATTCGCCGCCATCAGTTCAGCGAg     <  1:2178134/71‑1 (MQ=255)
                                                                cATATTGCTGGCGCTGAATTCCCCAGCCTGCCATCAGCATTGTGCGATTCGCCGCCATCAGTTCAGCGAGt    >  1:3100847/1‑71 (MQ=255)
                                                                  tatTGCTGGCGCTGAATTCCCCAGCCTGCCATCAGCATTGTGCGATTCGCCGCCATCAGTTCAGCGAGttt  <  1:3767300/71‑1 (MQ=255)
                                                                    |                                                                    
CCCGGTGTACCAATCTGCCCTAACATCGCCGCCAGCGTCACCAGCATCCAGTGTTTCTGTTCACCATATTGCTGGCGCTGAATTCCCCAGCCTGCCATCAGCATTGTGCGATTCGCCGCCATCAGTTCAGCGAGTTT  >  minE/1259779‑1259915

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: