Predicted mutation
evidence seq id position mutation freq annotation gene description
RA Exported 1,259,622 A→T 5.2% D434E (GAT→GAA 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
*Exported1,259,6220AT 5.2% 274.8 / 5.4 96D434E (GAT→GAAtorZtrimethylamine N‑oxide reductase system III, catalytic subunit
Reads supporting (aligned to +/- strand):  ref base A (51/40);  new base T (3/2);  total (54/42)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.11e-01

CATTTTGATATTCGGATACGTTTGTTCTTTACCGTTATGCTGATATTTGCCGCCAGGATTTTCCAGCGCATCGACGATACGGGCAACAGGAATCGCCGTCATTCCGCCATCATCAGCGGCCTCGCTGGCATGTCC  >  Exported/1259553‑1259687
                                                                     |                                                                 
cATTTTGATATTCGGATACGTTTGTTCTTTACCGTTATGCTGATATTTGCCGCCAGGATTTTCCAGCGCAt                                                                  >  1:3971528/1‑71 (MQ=255)
cATTTTGATATTCGGATACGTTTGTTCTTTACCGTTATGCTGATATTTGCCGCCAGGATTTTCCAGCGCAt                                                                  >  1:2593680/1‑71 (MQ=255)
cATTTTGATATTCGGATACGTTTGTTCTTTACCGTTATGCTGATATTTGCCGCCAGGATTTTCCAGCGCAt                                                                  >  1:2478823/1‑71 (MQ=255)
 aTTTTGATATTCGGATACGTTTGTTCTTTACCGTTATGCTGATATTTGCCGCCAGGATTTTCCAGCGCATc                                                                 <  1:5390111/71‑1 (MQ=255)
 aTTTTGATATTCGGATACGTTTGTTCTTTACCGTTATGCTGATATTTGCCGCCAGGATTTTCCAGCGCATc                                                                 <  1:1507028/71‑1 (MQ=255)
  ttttGATATTCGGATACGTTTGTTCTTTACCGTTATGCTGATATTTGCCGCCAGGATTTTCCAGCGCATc                                                                 >  1:3886777/1‑70 (MQ=255)
   tttGATATTCGGATACGTTTGTTCTTTACCGTTATGCTGATATTTGCCGCCAGGATTTTCCAGCGCATcga                                                               >  1:3963912/1‑71 (MQ=255)
    ttGATATTCGGATACGTTTGTTCTTTACCGTTATGCTGATATTTGCCGCCAGGATTTTCCAGCGCATcgac                                                              >  1:2727893/1‑71 (MQ=255)
    ttGATATTCGGATACGTTTGTTCTTTACCGTTATGCTGATATTTGCCGCCAGGATTTTCCAGCGCATcgac                                                              >  1:1933358/1‑71 (MQ=255)
     tGATATTCGGATACGTTTGTTCTTTACCGTTATGCTGATATTTGCCGCCAGGATTTTCCAGCGCATcgac                                                              <  1:4129135/70‑1 (MQ=255)
       atatTCGGATACGTTTGTTCTTTACCGTTATGCTGATATTTGCCGCCAGGATTTTCCAGCGCATcgacga                                                            <  1:3541201/70‑1 (MQ=255)
          ttCGGATACGTTTGTTCTTTACCGTTATGCTGATATTTGCCGCCAGGATTTTCCAGCGCATCGACGATAc                                                         <  1:3431187/70‑1 (MQ=255)
           tCGGATACGTTTGTTCTTTACCGTTATGCTGATATTTGCCGCCAGGATTTTCCAGCGCATcga                                                               <  1:1494828/63‑1 (MQ=255)
             ggATACGTTTGTTCTTTACCGTTATGCTGATATTTGCCGCCAGGATTTTCCAGCGCATCGACGATACGGGc                                                     <  1:750469/71‑1 (MQ=255)
                    tttGTTCTTTACCGTTATGCTGATATTTGCCGCCAGGATTTTCCAGCGCATCGACGATa                                                          <  1:1396129/59‑1 (MQ=255)
                    tttGTTCTTTACCGTTATGCTGATATTTGCCGCCAGGATTTTCCAGCGCATCGACGATa                                                          <  1:1846398/59‑1 (MQ=255)
                     ttGTTCTTTACCGTTATGCTGATATTTGCCGCCAGGATTTTCCAGCGCATcgacga                                                            >  1:3067841/1‑56 (MQ=255)
                     ttGTTCTTTACCGTTATGCTGATATTTGCCGCCAGGATTTTCCAGCGCATcgacga                                                            >  1:5074405/1‑56 (MQ=255)
                     ttGTTCTTTACCGTTATGCTGATATTTGCCGCCAGGATTTTCCAGCGCATcgacga                                                            >  1:3849606/1‑56 (MQ=255)
                     ttGTTCTTTACCGTTATGCTGATATTTGCCGCCAGGATTTTCCAGCGCATcgacga                                                            >  1:3731956/1‑56 (MQ=255)
                     ttGTTCTTTACCGTTATGCTGATATTTGCCGCCAGGATTTTCCAGCGCATcgacga                                                            >  1:4249166/1‑56 (MQ=255)
                     ttGTTCTTTACCGTTATGCTGATATTTGCCGCCAGGATTTTCCAGCGCATCGACGATACGGGCAACAGGaa                                             >  1:2301660/1‑71 (MQ=255)
                      tGTTCTTTACCGTTATGCTGATATTTGCCGCCAGGATTTTCCAGCGCATcga                                                               >  1:1171895/1‑52 (MQ=255)
                      tGTTCTTTACCGTTATGCTGATATTTGCCGCCAGGATTTTCCAGCGCATCGACGATACGGGCaa                                                   >  1:4416331/1‑64 (MQ=255)
                      tGTTCTTTACCGTTATGCTGATATTTGCCGCCAGGATTTTCCAGCGCATCGACGATACGGGCAACAGGAAt                                            >  1:567388/1‑71 (MQ=255)
                      tGTTCTTTACCGTTATGCTGATATTTGCCGCCAGGATTTTCCAGCGCATCGACGATACGGGCAACAGGAAt                                            >  1:5262237/1‑71 (MQ=255)
                      tGTTCTTTACCGTTATGCTGATATTTGCCGCCAGGATTTTCCAGCGCATCGACGATACGGGCAACAGGAAt                                            >  1:850539/1‑71 (MQ=255)
                       gTTCTTTACCGTTATGCTGATATTTGCCGCCAGGATTTTCCAGCGCATCGACGATACGGGCATCAGGAATc                                           >  1:168942/1‑71 (MQ=255)
                       gTTCTTTACCGTTATGCTGATATTTGCCGCCAGGATTTTCCAGCGCATCGACGATACGGGCAACAGGAAt                                            <  1:4522051/70‑1 (MQ=255)
                       gTTCTTTACCGTTATGCTGATATTTGCCGCCAGGATTTTCCAGCGCATCGACGATACGGGCAACAGGAATc                                           >  1:3700420/1‑71 (MQ=255)
                        ttCTTTACCGTTATGCTGATATTTGCCGCCAGGATTTTCCAGCGCATcgacga                                                            >  1:2298950/1‑53 (MQ=255)
                        ttCTTTACCGTTATGCTGATATTTGCCGCCAGGATTTTCCAGCGCATc                                                                 <  1:4758252/48‑1 (MQ=255)
                        ttCTTTACCGTTATGCTGATATTTGCCGCCAGGATTTTCCAGCGCATc                                                                 <  1:1679952/48‑1 (MQ=255)
                        ttCTTTACCGTTATGCTGATATTTGCCGCCAGGATTTTCCAGCGCATCGACGATACGGGCAACAGGa                                              <  1:5199668/67‑1 (MQ=255)
                          cTTTACCGTTATGCTGATATTTGCCGCCAGGATTTTCCAGCGCATcga                                                               <  1:3090719/48‑1 (MQ=255)
                          cTTTACCGTTATGCTGATATTTGCCGCCAGGATTTTCCAGCGCATcga                                                               >  1:2381079/1‑48 (MQ=255)
                          cTTTACCGTTATGCTGATATTTGCCGCCAGGATTTTCCAGCGCATCGACGATACGGGCAACAGGAAt                                            >  1:5448335/1‑67 (MQ=255)
                          cTTTACCGTTATGCTGATATTTGCCGCCAGGATTTTCCAGCGCATCGACGATACGGGCAACAGGAAt                                            >  1:429945/1‑67 (MQ=255)
                          cTTTACCGTTATGCTGATATTTGCCGCCAGGATTTTCCAGCGCATCGACGATACGGGCAACAGGAATCGcc                                        >  1:2998410/1‑71 (MQ=255)
                           tttACCGTTATGCTGATATTTGCCGCCAGGATTTTCCAGCGCATc                                                                 <  1:3247156/45‑1 (MQ=255)
                            ttACCGTTATGCTGATATTTGCCGCCAGGATTTTCCAGCGCATCGACGATACGGGCAACAGGAATCGCCGt                                      >  1:128213/1‑71 (MQ=255)
                             tACCGTTATGCTGATATTTGCCGCCAGGATTTTCCAGCGCATCGACGATACGGGCAACAGGAATCGCCGt                                      <  1:4964833/70‑1 (MQ=255)
                             tACCGTTATGCTGATATTTGCCGCCAGGATTTTCCAGCGCATCGACGATACGGGCAACAGGAATCGCCGt                                      <  1:4261229/70‑1 (MQ=255)
                             tACCGTTATGCTGATATTTGCCGCCAGGATTTTCCAGCGCATCGACGATACGGGCAACAGGAATCGCCGt                                      <  1:340776/70‑1 (MQ=255)
                             tACCGTTATGCTGATATTTGCCGCCAGGATTTTCCAGCGCATCGACGATACGGGCAACAGGAATCGCCGt                                      >  1:1514098/1‑70 (MQ=255)
                               ccGTTGTGTTGATATGCGCCACCAGGGTTTTCCAGTGCTTCAACAATGCGGGCAACAggg                                              >  1:3145899/1‑59 (MQ=255)
                               ccGTTGTGTTGATATGCGCCACCAGGGTTTTCCAGTGCTTCAACAATGCGGGCAACAggg                                              >  1:5087717/1‑59 (MQ=255)
                               ccGTTGTGTTGATATGCGCCACCAGGGTTTTCCAGTGCTTCAACAATGCGGGCAACAGGGAttt                                          >  1:1695375/1‑62 (MQ=255)
                                    aTGCTGATATTTGCCGCCAGGATTTTCCAGCGCATCGACGATACGGGCAACAGGAATCGCCGTCATTCCg                               <  1:5195118/70‑1 (MQ=255)
                                     tGCTGATATTTGCCGCCAGGATTTTCCAGCGCATCGACGATACGGGCAACAGGAATCGCCGTCATTCCGc                              >  1:1056701/1‑70 (MQ=255)
                                     tGCTGATATTTGCCGCCAGGATTTTCCAGCGCATCGACGATACGGGCAACAGGAATCGCCGTCATTCCGc                              >  1:4152849/1‑70 (MQ=255)
                                      gCTGATATTTGCCGCCAGGATTTTCCAGCGCATcgac                                                              >  1:2703606/1‑37 (MQ=255)
                                      gCTGATATTTGCCGCCAGGATTTTCCAGCGCATcgac                                                              >  1:4065946/1‑37 (MQ=255)
                                      gCTGATATTTGCCGCCAGGATTTTCCAGCGCATcgac                                                              >  1:5494210/1‑37 (MQ=255)
                                       ttGATATGCGCCACCAGGGTTTTCCAGTGCTTCAACAATGCGGGCAACAGGGAt                                            <  1:5642868/53‑1 (MQ=255)
                                        tGATATTTGCCGCCAGGATTTTCCAGCGCATCGACGATa                                                          <  1:3528062/39‑1 (MQ=255)
                                        tGATATGCGCCACCAGGGTTTTCCAGTGCTTCAACAATGCGGGCAACAGGGAt                                            <  1:5274536/53‑1 (MQ=255)
                                          atatTTGCCGCCAGGATTTTCCAGCGCATCGACGATACGGGCAACAGGAATCGCCGTCATTCCGCcatca                         <  1:5301244/70‑1 (MQ=255)
                                           tatTTGCCGCCAGGATTTTCCAGCGCATCGACGATACGGGCAACAGGAATCGCCGTCATTc                                 >  1:4898602/1‑61 (MQ=255)
                                           tatTTGCCGCCAGGATTTTCCAGCGCATCGACGATACGGGCAACAGGAATCGCCGTCATTc                                 >  1:3113091/1‑61 (MQ=255)
                                           tatTTGCCGCCAGGATTTTCCAGCGCATCGACGATACGGGCAACAGGAATCGCCGTCATTc                                 >  1:4371361/1‑61 (MQ=255)
                                           tatTTGCCGCCAGGATTTTCCAGCGCATCGACGATACGGGCAACAGGAATCGCCGTCATTCCGCcatcat                        >  1:489677/1‑70 (MQ=255)
                                             tttGCCGCCAGGATTTTCCAGCGCATCGACGATACGGGCAACAGGaa                                             >  1:4554316/1‑47 (MQ=255)
                                             tttGCCGCCAGGATTTTCCAGCGCATCGACGATACGGGCAACAGGa                                              <  1:1910830/46‑1 (MQ=255)
                                               tGCCGCCAGGATTTTCCAGCGCATCGACGATACGGGCAACAGGAATCGCCGTCATTCCGCCATCATCAGCg                   >  1:332939/1‑71 (MQ=255)
                                                 ccgccAGGATTTTCCAGCGCATCGACGATACGGGCAACAGGAATCGCCGTCATTCCGc                              <  1:476519/58‑1 (MQ=255)
                                                  cgccAGGATTTTCCAGCGCATCGACGATACGGGCAACAGGAATCGCCGTCATTCCGCcatcat                        >  1:1364440/1‑63 (MQ=255)
                                                  cgccAGGATTTTCCAGCGCATCGACGATACGGGCAACAGGAATCGCCGTCATTCCGCCATCATCAGCGGcc                >  1:82557/1‑71 (MQ=255)
                                                  cgccAGGATTTTCCAGCGCATCGACGATACGGGCAACAGGAATCGCCGTCATTCCGCCATCATCAGCGGcc                <  1:1653200/71‑1 (MQ=255)
                                                  cgccAGGATTTTCCAGCGCATCGACGATACGGGCAACAGGAATCGCCGTCATTCCGCCATCATCAGCGGcc                <  1:4989511/71‑1 (MQ=255)
                                                      aGGATTTTCCAGCGCATCGACGATACGGGCAACAGGAATCGcc                                        <  1:4784928/43‑1 (MQ=255)
                                                      aGGATTTTCCAGCGCATCGACGATACGGGCAACAGGAATCGCCGTCATTCCGCCATCATCAGCGGCCTCGc            >  1:2855123/1‑71 (MQ=255)
                                                        gATTTTCCAGCGCATCGACGATACGGGCAACAGGa                                              <  1:1339194/35‑1 (MQ=255)
                                                        gATTTTCCAGCGCATCGACGATACGGGCAACAGGa                                              <  1:3863570/35‑1 (MQ=255)
                                                        gATTTTCCAGCGCATCGACGATACGGGCAACAGGAATCGCCGt                                      <  1:1315152/43‑1 (MQ=255)
                                                        gATTTTCCAGCGCATCGACGATACGGGCAACAGGAATCGCCGt                                      <  1:5182770/43‑1 (MQ=255)
                                                        gATTTTCCAGCGCATCGACGATACGGGCAACAGGAATCGCCGt                                      <  1:5176533/43‑1 (MQ=255)
                                                        gATTTTCCAGCGCATCGACGATACGGGCAACAGGAATCGCCGt                                      <  1:4889252/43‑1 (MQ=255)
                                                        gATTTTCCAGCGCATCGACGATACGGGCAACAGGAATCGCCGTCATTCCGCCatcatc                       <  1:4849533/58‑1 (MQ=255)
                                                         aTTTTCCAGCGCATCGACGATACGGGCAACAGGAATCGCCGTc                                     <  1:782481/43‑1 (MQ=255)
                                                         aTTTTCCAGCGCATCGACGATACGGGCAACAGGAATCGCCGTc                                     <  1:2801507/43‑1 (MQ=255)
                                                         aTTTTCCAGCGCATCGACGATACGGGCAACAGGAATCGCCGTCATTCCGCCATCATCAGCGGCCTCGCTgg         <  1:210047/71‑1 (MQ=255)
                                                          ttttCCAGCGCATCGACGATACGGGCAACAGGAAt                                            >  1:1115151/1‑35 (MQ=255)
                                                          ttttCCAGCGCATCGACGATACGGGCAACAGGAAt                                            >  1:4680944/1‑35 (MQ=255)
                                                          ttttCCAGCGCATCGACGATACGGGCAACAGGAAt                                            >  1:1220249/1‑35 (MQ=255)
                                                          ttttCCAGCGCATCGACGATACGGGCAACAGGAATCGCCGTCATTCCGCCATCATCAg                     >  1:4939809/1‑58 (MQ=255)
                                                          ttttCCAGCGCATCGACGATACGGGCAACAGGAATCGCCGTCATTCCGCCATCATCAg                     >  1:1573795/1‑58 (MQ=255)
                                                           tttCCAGCGCATCGACGATACGGGCAACAGGAATCGCCGTCATTCCGCcatcat                        <  1:3976185/54‑1 (MQ=255)
                                                           tttCCAGCGCATCGACGATACGGGCAACAGGAATCGCCGTCATTCCGCcatcat                        <  1:1938859/54‑1 (MQ=255)
                                                           tttCCAGCGCATCGACGATACGGGCAACAGGAATCGCCGTCATTCCGCcatcat                        <  1:634346/54‑1 (MQ=255)
                                                           tttCCAGCGCATCGACGATACGGGCAACAGGAATCGCCGTCATTCCGCcatcat                        <  1:895565/54‑1 (MQ=255)
                                                             tCCAGCGCATCGACGATACGGGCAACAGGAATCGCCGTCa                                    >  1:3915666/1‑40 (MQ=255)
                                                             tCCAGCGCATCGACGATACGGGCAACAGGAATCGCCGTCa                                    >  1:1238459/1‑40 (MQ=255)
                                                                aGCGCATCGACGATACGGGCAACAGGAATCGCCGTCATTCCGCCATCATCAGCGGCCTCGCTGGCATGTcc  >  1:1461784/1‑71 (MQ=255)
                                                                aGCGCATCGACGATACGGGCAACAGGAATCGCCGTCATTCCGCCATCATCAGCGGCCTCGCTGGCATGTcc  >  1:1176889/1‑71 (MQ=255)
                                                                aGCGCATCGACGATACGGGCAACAGGAATCGCCGTCATTCCGCCATCATCAGCGGCCTCGCTGGCATGTc   >  1:4970305/1‑70 (MQ=255)
                                                                     |                                                                 
CATTTTGATATTCGGATACGTTTGTTCTTTACCGTTATGCTGATATTTGCCGCCAGGATTTTCCAGCGCATCGACGATACGGGCAACAGGAATCGCCGTCATTCCGCCATCATCAGCGGCCTCGCTGGCATGTCC  >  Exported/1259553‑1259687

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