Predicted mutation
evidence seq id position mutation freq annotation gene description
RA minE 1,259,838 T→C 22.2% 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,8380TC22.2% 134.8 / 37.4 72E362E (GAA→GAGtorZtrimethylamine N‑oxide reductase system III, catalytic subunit
Reads supporting (aligned to +/- strand):  ref base T (25/31);  new base C (11/5);  total (36/36)
Fisher's exact test for biased strand distribution p-value = 1.55e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 3.83e-01

CAAAACCGCCGCCCGGTGTACCAATCTGCCCTAACATCGCCGCCAGCGTCACCAGCATCCAGTGTTTCTGTTCACCATATTGCTGGCGCTGAATTCCCCAGCCTGCCATCAGCATTGTGCGATTCGCCGCCATCAGTT  >  minE/1259768‑1259905
                                                                      |                                                                   
cAAAACCGCCGCCCGGTGTACCAATCTGCCCTAACATCGCCGCCAGCGTCACCAGCATCCAGTGTTTCTGt                                                                     >  1:237533/1‑71 (MQ=255)
cAAAACCGCCGCCCGGTGTACCAATCTGCCCTAACATCGCCGCCAGCGTCACCAGCATCCAGTGTTTCTGt                                                                     >  1:1587592/1‑71 (MQ=255)
      cgccgcCCGGTGTACCAATCTGCCCTAACATCGCCGCCAGCGTCACCAGCATCCAGTGTTTCTGTTCACCa                                                               >  1:956166/1‑71 (MQ=255)
      cgccgcCCGGTGTACCAATCTGCCCTAACATCGCCGCCAGCGTCACCAGCATCCAGTGTTTCTGTTCACCa                                                               <  1:2071491/71‑1 (MQ=255)
           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:1504740/71‑1 (MQ=255)
              ggTGTACCAATCTGCCCTAACATCGCCGCCAGCGTCACCAGCATCCAGTGTTTCTGTTCACCATATTGCTg                                                       >  1:1156444/1‑71 (MQ=255)
               gtgtACCAATCTGCCCTAACATCGCCGCCAGCGTCACCAGCATCCAGTGTTTCTGtt                                                                    >  1:2807651/1‑57 (MQ=255)
                 gtACCAATCTGCCCTAACATCGCCGCCAGCGTCACCAGCATCCAGTGTTTCTGTTCACCATATTGCTGgcg                                                    >  1:3950268/1‑71 (MQ=255)
                  tACCAATCTGCCCTAACATCGCCTCCAGCGTCACCAGCATCCAGTGTTTCTGTTCACCATATTGCTGgcg                                                    <  1:1955581/70‑1 (MQ=255)
                    ccAATCTGCCCTAACATCGCCGCCAGCGTCACCAGCATCCAGTGTTTCTGtt                                                                    <  1:332172/52‑1 (MQ=255)
                    ccAATCTGCCCTAACATCGCCGCCAGCGTCACCAGCATCCAGTGTTTCTGt                                                                     <  1:1372389/51‑1 (MQ=255)
                       atttGCCCCAACATTGCTGCCAGCGTGACGATCATCCAGTGTTTTTGCTCACCAAACTGTTGGCGCTGCAt                                              >  1:1223111/4‑71 (MQ=255)
                       aTCTGCCCTAACATCGCCGCCAGCGTCACCAGCATCCAGTGTTTCTGtt                                                                    <  1:3436580/49‑1 (MQ=255)
                       aTCTGCCCTAACATCGCCGCCAGCGTCACCAGCATCCAGTGTTTCTGtt                                                                    <  1:2919854/49‑1 (MQ=255)
                        tCTGCCCTAACATCGCCGCCAGCGTCACCAGCATCCAGTGTTTCTGTTCa                                                                  >  1:2617609/1‑50 (MQ=255)
                          tGCCCTAACATCGCCGCCAGCGTCACCAGCATCCAGTGTTTCTGTTCACCATATTGCTGGCGCTg                                                 >  1:3416274/1‑65 (MQ=255)
                          tGCCCTAACATCGCCGCCAGCGTCACCAGCATCCAGTGTTTCTGTTCACCATATTGCTGGCGCTg                                                 >  1:3121549/1‑65 (MQ=255)
                            cccTAACATCGCCGCCAGCGTCACCAGCATCCAGTGTTTCTGTTCACCATAtt                                                           <  1:3882319/53‑1 (MQ=255)
                             ccTAACATCGCCGCCAGCGTCACCAGCATCCAGTGTTTCTGt                                                                     <  1:3317346/42‑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)
                                aaCATTGCTGCCAGCGTGACGATCATCCAGTGTTTTTGCTCACCaaa                                                             >  1:2259757/1‑45 (MQ=12)
                                aaCATCGCCGCCAGCGTCACCAGCATCCAGTGTTTCTGtt                                                                    >  1:3326263/1‑40 (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:2221078/1‑71 (MQ=255)
                                      gccgccAGCGTCACCAGCATCCAGTGTTTCTGTTCACCATATTGCTGGCGCTGAATTCCCCAGCCTGCCAt                               >  1:2216343/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:3621450/1‑71 (MQ=255)
                                           cAGCGTCACCAGCATCCAGTGTTTCTGTTCACCATATTGCTGGCGCTGAATTCCCCAGCCTGCCATCAGCa                          <  1:1507521/71‑1 (MQ=255)
                                           cAGCGTCACCAGCATCCAGTGTTTCTGTTCACCATATTGCTGGCGCTGAATTCCCCAGCCTGCCATCAGCa                          <  1:3841247/71‑1 (MQ=255)
                                           cAGCGTCACCAGCATCCAGTGTTTCTGTTCACCATATTGCTGGCGCTGAATTCCCCAGCCTGCCATCAGCa                          <  1:1910719/71‑1 (MQ=255)
                                           cAGCGTCACCAGCATCCAGTGTTTCTGTTCACCATATTGCTGGCGCTGAATTCCCCAGCCTGCCATCAGCa                          <  1:462589/71‑1 (MQ=255)
                                           cAGCGTCACCAGCATCCAGTGTTTCTGTTCACCATATTGCTGGCGCTGAATTCCCCAGCCTGCCATCAGCa                          <  1:235307/71‑1 (MQ=255)
                                            aGCGTCACCAGCATCCAGTGTTTCTGTTCACCATATTGCTGGCGCTGAATTCCCCAGCCTGCCATCAGCAt                         >  1:792378/1‑71 (MQ=255)
                                             gCGTGACGATCATCCAGTGTTTTTGCTCACCAAACTGTTGGCGCTGCATTCCCCAGCCTGCCATCAGCAtg                        <  1:1881162/71‑2 (MQ=255)
                                             gCGTGACGATCATCCAGTGTTTTTGCTCACCAAACTGTTGGCGCTGCATTCCCCAGCCTGCCATCAGCAtg                        <  1:480872/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:2928067/71‑1 (MQ=255)
                                                      gCATCCAGTGTTTCTGTTCACCATATTGCTGGCGCTGAATTCCCCAGCCTGCCATCAGCATTGTGCGATTc               >  1:3191134/1‑71 (MQ=255)
                                                        aTCCAGTGTTTCTGTTCACCATATTGCTGGCGCTGAATTCCCCAGCCTGCCATCAGCATTGTGCGATTcgc             <  1:802895/71‑1 (MQ=255)
                                                         tCCAGTGTTTCTGTTCACCATATTGCTGGCGCTGAATTCCCCAGCCTGCCATCAGCATTGTGCGATTcgcc            <  1:715872/71‑1 (MQ=255)
                                                           cAGTGTTTTTGCTCACCAAACTGTTGGCGCTGCATTCCCCAGCCTGCCATc                              >  1:1240331/1‑51 (MQ=255)
                                                           cAGTGTTTTTGCTCACCAAACTGTTGGCGCTGCATTCCCCAGCCTGCCATc                              >  1:1226739/1‑51 (MQ=255)
                                                           cAGTGTTTTTGCTCACCAAACTGTTGGCGCTGCATTCCCCAGCCTGCCATc                              >  1:1671019/1‑51 (MQ=255)
                                                           cAGTGTTTTTGCTCACCAAACTGTTGGCGCTGCATTCCCCAGCCTGCCATc                              >  1:3445555/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:1407074/71‑1 (MQ=255)
                                                              tgtTTCTGTTCACCATATTGCTGGCGCTGAATTCCCCAGCCTGCCATCAGCATTGTGCGATTCGCCGCCAt       <  1:1126505/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:3924097/1‑71 (MQ=255)
                                                                   cTGTTCACCATATTGCTGGCGCTGAATTCCCCAGCCTGCCATCAGCATTGTGCGATTCGCCGCCATCAGtt  >  1:310941/1‑71 (MQ=255)
                                                                      |                                                                   
CAAAACCGCCGCCCGGTGTACCAATCTGCCCTAACATCGCCGCCAGCGTCACCAGCATCCAGTGTTTCTGTTCACCATATTGCTGGCGCTGAATTCCCCAGCCTGCCATCAGCATTGTGCGATTCGCCGCCATCAGTT  >  minE/1259768‑1259905

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: