Predicted mutation
evidence seq id position mutation freq annotation gene description
RA minE 1,259,850 C→T 21.4% Q358Q (CAG→CAA 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,8500CT21.4% 141.0 / 32.7 70Q358Q (CAG→CAAtorZtrimethylamine N‑oxide reductase system III, catalytic subunit
Reads supporting (aligned to +/- strand):  ref base C (23/32);  new base T (10/5);  total (33/37)
Fisher's exact test for biased strand distribution p-value = 1.43e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

CGGTGTACCAATCTGCCCTAACATCGCCGCCAGCGTCACCAGCATCCAGTGTTTCTGTTCACCATATTGCTGGCGCTGAATTCCCCAGCCTGCCATCAGCATTGTGCGATTCGCCGCCATCAGTTCAGCGAGTTTGAC  >  minE/1259781‑1259918
                                                                     |                                                                    
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:3121549/1‑65 (MQ=255)
             tGCCCTAACATCGCCGCCAGCGTCACCAGCATCCAGTGTTTCTGTTCACCATATTGCTGGCGCTg                                                              >  1:3416274/1‑65 (MQ=255)
                ccTAACATCGCCGCCAGCGTCACCAGCATCCAGTGTTTCTGTTCACCATATTGCTGacgc                                                                >  1:2492214/1‑60 (MQ=255)
                ccTAACATCGCCGCCAGCGTCACCAGCATCCAGTGTTTCTGTTCACCATATTGCTGGCGCTGAATTCCCCa                                                     <  1:3514365/71‑1 (MQ=255)
                ccTAACATCGCCGCCAGCGTCACCAGCATCCAGTGTTTCTGTTCACCATATTGCTGGCGCTGAATTCCCCa                                                     <  1:843594/71‑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:3002658/1‑45 (MQ=255)
                         gccgccAGCGTCACCAGCATCCAGTGTTTCTGTTCACCATATTGc                                                                      >  1:1904084/1‑45 (MQ=255)
                         gccgccAGCGTCACCAGCATCCAGTGTTTCTGTTCACCATATTGCTGGCGCTGAATTCCCCAGCCTGCCAt                                            >  1:2216343/1‑71 (MQ=255)
                         gccgccAGCGTCACCAGCATCCAGTGTTTCTGTTCACCATATTGCTGGCGCTGAATTCCCCAGCCTGCCAt                                            >  1:2221078/1‑71 (MQ=255)
                          ctgccAGCGTGACGATCATCCAGTGTTTTTGCTCACCAAACTGTTGGCGCTGCATTCCCCAGCCTGCCAt                                            <  1:2241161/68‑1 (MQ=255)
                          ctgccAGCGTGACGATCATCCAGTGTTTTTGCTCACCAAACTGTTGGCGCTGCATTCCCCAGCCTGCCAt                                            >  1:3379490/3‑70 (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:1507521/71‑1 (MQ=255)
                              cAGCGTCACCAGCATCCAGTGTTTCTGTTCACCATATTGCTGGCGCTGAATTCCCCAGCCTGCCATCAGCa                                       <  1:235307/71‑1 (MQ=255)
                              cAGCGTCACCAGCATCCAGTGTTTCTGTTCACCATATTGCTGGCGCTGAATTCCCCAGCCTGCCATCAGCa                                       <  1:1910719/71‑1 (MQ=255)
                              cAGCGTCACCAGCATCCAGTGTTTCTGTTCACCATATTGCTGGCGCTGAATTCCCCAGCCTGCCATCAGCa                                       >  1:3621450/1‑71 (MQ=255)
                              cAGCGTCACCAGCATCCAGTGTTTCTGTTCACCATATTGCTGGCGCTGAATTCCCCAGCCTGCCATCAGCa                                       <  1:462589/71‑1 (MQ=255)
                              cAGCGTCACCAGCATCCAGTGTTTCTGTTCACCATATTGCTGGCGCTGAATTCCCCAGCCTGCCATCAGCa                                       <  1:3841247/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:1226739/1‑51 (MQ=255)
                                              cAGTGTTTTTGCTCACCAAACTGTTGGCGCTGCATTCCCCAGCCTGCCATc                                           >  1:3445555/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)
                                                                   tGCTGGCGCTGAATTCCCCAGCCTGCCATCAGCATTGTGCGATTCGCCGCCAt                    >  1:1236801/1‑53 (MQ=255)
                                                                   tGCTGGCGCTGAATTCCCCAGCCTGCCATCAGCATTGTGCGATTCGCCGCCAACAGTTCAGCGAGTTTGAc  >  1:1866554/1‑71 (MQ=255)
                                                                    gCTGGCGCTGAATTCCCCAGCCTGCCATCAGCATTGTGCGATTCGCCGCCATCAGTTCAGCGAGTTTGAc  <  1:780591/70‑1 (MQ=255)
                                                                     |                                                                    
CGGTGTACCAATCTGCCCTAACATCGCCGCCAGCGTCACCAGCATCCAGTGTTTCTGTTCACCATATTGCTGGCGCTGAATTCCCCAGCCTGCCATCAGCATTGTGCGATTCGCCGCCATCAGTTCAGCGAGTTTGAC  >  minE/1259781‑1259918

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: