Predicted mutation
evidence seq id position mutation freq annotation gene description
RA Exported 1,259,787 A→G 27.3% G379G (GGT→GGC 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,7870AG27.3% 79.5 / 33.7 55G379G (GGT→GGCtorZtrimethylamine N‑oxide reductase system III, catalytic subunit
Reads supporting (aligned to +/- strand):  ref base A (16/24);  new base G (10/5);  total (26/29)
Fisher's exact test for biased strand distribution p-value = 1.29e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 4.18e-01

GCCGCCAACACGCGTCGGGTTGCCACCGTTGGAGTAATGATAAGAGAAACCAAAACCGCCGCCCGGTGTACCAATCTGCCCTAACATCGCCGCCAGCGTCACCAGCATCCAGTGTTTCTGTTCACCATATTGCTGGC  >  Exported/1259718‑1259854
                                                                     |                                                                   
gccgccAACACGCGTCGGGTTGCCACCGTTGGAGTAATGATAAGAGAAACCAAAACCGCCGCCCGGTGTAc                                                                    >  1:5665775/1‑71 (MQ=255)
 ccgccAACACGCGTCGGGTTGCCACCGTTGGAGTAATGATAAGAGAAACCAAAACCGCCGCCCGGTGTAcc                                                                   <  1:1749692/71‑1 (MQ=255)
  cgccAACACGCGTCGGGTTGCCACCGTTGGAGTAATGATAAGAGAAACCAAAACCGCCGCCCGGTGTACCa                                                                  <  1:1386230/71‑1 (MQ=255)
  cgccAACACGCGGCGGGTTGCCACCGTTGGAGTAATGATAAGAGAAACCAAAACCGCCGCCCGGTGTACCa                                                                  <  1:5047619/71‑1 (MQ=255)
     cAACACGCGTCGGGTTGCCACCGTTGGAGTAATGATAAGAGAAACCAAAACCGCCGCCCGGTGTACCAATc                                                               >  1:2102408/1‑71 (MQ=255)
      aaCACGCGTCGGGTTGCCACCGTTGGAGTAATGATAAGAGAAACCAAAACCGCCGCCCGGTGTACCAATCt                                                              >  1:109426/1‑71 (MQ=255)
            cgTCGGGTTGCCACCGTTGGAGTAATGATAAGAGAAACCAAAACCGCCGCCCGGTGTACCAATCTGCCCt                                                         <  1:1088197/70‑1 (MQ=255)
            cgTCGGGTTGCCACCGTTGGAGTAATGATAAGAGAAACCAAAACCGCCGCCCGGTGTACCAATCTGCCCTa                                                        <  1:1946550/71‑1 (MQ=255)
              tggggTTACCACCATTGGCAAAATGGTAAGAAAGACCAAAACCGCCGCCGGGTGTGCCGATTTGCCCCaa                                                       >  1:3184431/3‑70 (MQ=255)
              tCGGGTTGCCACCGTTGGAGTAATGATAAGAGAAACCAAAACCGCCGCCCGGTGTACCAATCTGCCCTa                                                        <  1:742430/69‑1 (MQ=255)
                gggTTGCCACCGTTGGAGTAATGATAAGAGAAACCAAAACCGCCGCCCGGTGTa                                                                     >  1:4083805/1‑54 (MQ=255)
                  gTTACCACCATTGGCAAAATGGTAAGAAAGACCAAAACCGCCGCCGGGTGTGCCGATTTGCCCCAACAttg                                                  >  1:1004409/1‑69 (MQ=255)
                   ttaccaccaTTGGCAAAATGGTAAGAAAGACCAAAACCGCCGCCGGGTGTGCCGATTTGCCCCAACATtgc                                                 >  1:2755281/4‑71 (MQ=255)
                   ttGCCACCGTTGGAGTAATGATAAGAGAAACCAAAACCGCCGCCCGGTGTACCAATCTGCCCTAACATcgc                                                 >  1:4433378/1‑71 (MQ=255)
                    taccaccaTTGGCAAAATGGTAAGAAAGACCAAAACCGCCGCCGGGTGTGCCGATTTGCCCCAACATtgc                                                 >  1:239651/3‑70 (MQ=255)
                      ccaccaTTGGCAAAATGGTAAGAAAGACCAAAACCGCCGCCGGGTGTGcc                                                                   <  1:5319517/50‑1 (MQ=255)
                              ggAGTAATTATAAGAGAAACCAAAACCGCCGCCCGGTGTACCa                                                                  <  1:4634668/43‑1 (MQ=255)
                               gAGTAATGATAAGAGAAACCAAAACCGCCGCCCGGTGTACCAATCTGCCCTAACATCGCCGCCAGCGTCAc                                     <  1:2844667/71‑1 (MQ=255)
                               gAGTAATGATAAGAGAAACCAAAACCGCCGCCCGGTGTACCAATCTGCCCTAACATCGCCGCCAGCGTCAc                                     <  1:3032831/71‑1 (MQ=255)
                                aGTAATGATAAGAGAAACCAAAACCGCCGCCCGGTGTACCAATCTGCCCTAACATCGCCGCCAGCGTCAc                                     >  1:4600250/1‑70 (MQ=255)
                                  aaaTGGTAAGAAAGACCAAAACCGCCGCCGGGTGTGCCGATTTGCCCCAACATTGCTGCCAGc                                          <  1:3618007/62‑1 (MQ=255)
                                  aaaTGGTAAGAAAGACCAAAACCGCCGCCGGGTGTGCCGATTTGCCCCAACATTGCTGCCAGCGTGAcga                                   <  1:419678/69‑3 (MQ=255)
                                     tGATAAGAGAAACCAAAACCGCCGCCCGGTGTACCa                                                                  <  1:1196808/36‑1 (MQ=255)
                                       gtAAGAAAGACCAAAACCGCCGCCGGGTGTGCCGATTTGCCCCAACATTGCTGCCAGCGTGACGATCATcc                             <  1:138006/70‑1 (MQ=255)
                                       gtAAGAAAGACCAAAACCGCCGCCGGGTGTGCCGATTTGCCCCAACATTGCTGCCAGCGTGACGATCATcc                             <  1:4525914/70‑1 (MQ=255)
                                        tAAGAGAAACCAAAACCGCCGCCCGGTGTACCAATCTGCCCTAACATCGCCGCCAGCGTCACCAGCATc                              >  1:3898621/1‑69 (MQ=255)
                                           gagaAACCAAAACCGCCGCCCGGTGTACCAATCTGCCCTAACATCGCCGCCAGCGTCACCAGc                                 <  1:1149535/63‑1 (MQ=255)
                                           gagaAACCAAAACCGCCGCACGGTGTACCAATCTGCCCTAACATCGCCGCCAGCGTCACCAGc                                 <  1:248870/63‑1 (MQ=255)
                                             gaAACCAAAACCGCCGCCCGGTGTACCAATCTGCCCTAACAt                                                    >  1:3003638/1‑42 (MQ=255)
                                             gaAACCAAAACCGCCGCCCGGTGTACCAATCTGCCCTAACATCGCCGCCAGCGTCACCAGCATCCAGTGtt                       >  1:2201941/1‑71 (MQ=255)
                                              agaCCAAAACCGCCGCCGGGTGTGCCGATTTGCCCCAACATTGCTGCCAGCGTGACGATCATCCa                            >  1:402649/3‑65 (MQ=255)
                                               gaCCAAAACCGCCGCCGGGTGTGCCGATTTGCCCCAACATTGCTGCCAGCGTGACGATCATCCAGTGtttt                     >  1:4128930/2‑70 (MQ=255)
                                               aaCCAAAACCGCCGCCCGGTGTACCAATCTGCCCTAACATCGCCGCCAGc                                          <  1:3615561/50‑1 (MQ=255)
                                               aaCCAAAACCGCCGCCCGGTGTACCAATCTGCCCTAACATCGCCGCCAGCGTCACCAGCa                                <  1:908069/60‑1 (MQ=255)
                                               aaCCAAAACCGCCGCCCGGTGTACCAATCTGCCCTAACATCGCCGCCAGCGTCACCAGCATCCAGTGTTTc                     <  1:2122035/71‑1 (MQ=255)
                                                aCCAAAACCGCCGCCCGGTGTACCAATCTGCCCTAACATCGCCGCCAGCGTCACCAGCATCCAGTGTTTc                     >  1:2733240/1‑70 (MQ=255)
                                                 ccAAAACCGCCGCCCGGTGTACCAATCTGCCCTAACATCGCCGCCAGCGTCACCa                                   >  1:2902326/1‑55 (MQ=255)
                                                  cAAAACCGCCGCCGGGTGTGCCGATTTGCCCCAACATtgct                                                >  1:3988192/1‑40 (MQ=12)
                                                  cAAAACCGCCGCCGGGTGTGCCGATTTGCCCCAACATTGCTGCCAGCGTGACGATCATCCAGTGtttt                     >  1:3137196/1‑67 (MQ=255)
                                                  cAAAACCGCCGCCCGGTGTACCAATCTGCCCTAACATCGCCGCCAGCGTCACCAGCATcc                             <  1:1660741/60‑1 (MQ=255)
                                                  cAAAACCGCCGCCCGGTGTACCAATCTGCCCTAACATCGCCGCCAGCGTCACCAGCATCCAgt                          <  1:3178490/63‑1 (MQ=255)
                                                     aaCCGCCGCCCGGTGTACCAATCTGCCCTAACATCGCCGCCAGCGTCACCAGCATCCAGTGTTTCTGTTCa               >  1:5384322/1‑71 (MQ=255)
                                                     aaCCGCCGCCCGGTGTACCAATCTGCCCTAACATCGCCGCCAGCGTCACCAGCATCCAGTGTTTCTGTTCa               >  1:2015270/1‑71 (MQ=255)
                                                        cgccgcCCGGTGTACCAATCTGCCCTAACATCGCCGCCAGCGTCACCAGCATCCAGTGTTTCTGt                  <  1:317868/65‑1 (MQ=255)
                                                         gccgccGGGTGTGCCGATTTGCCCCAACATTGCTGCCAGCGTGACGATCATCCAGTGTTTTTGCTCACCa            >  1:2928180/1‑70 (MQ=255)
                                                         gccgccCGGTGTACCAATCTGCCCTAACATCGCCGCCAGCGTCACCAGCATCCAGTGTTTCTGTTCACCa            >  1:3376499/1‑70 (MQ=255)
                                                           cgccCGGTGTACCAATCTGCCCTAACATCGCCGCCAGCGTCACCAGCATCCAGTGTTTCTg                   <  1:5535675/61‑1 (MQ=255)
                                                            gccCGGTGTACCAATCTGCCCTAACATCGCCGCCAGCGTCACCAGCATCCAGTGttt                      <  1:2246028/57‑1 (MQ=255)
                                                            gccCGGTGTACCAATCTGCCCTAACATCGCCGCCAGCGTCACCAGCATCCAGTGTTTCTGTTCACCATAtt        <  1:4789950/71‑1 (MQ=255)
                                                             ccgggTGTGCCGATTTGCCCCAACATTGCTGCCAGCGTGACGATCATCCAGTGTTTTTGCTCa               >  1:3004314/4‑63 (MQ=255)
                                                             cccGGTGTACCAATCTGCCCTAACATCGCCGCCAGCGTCACCAGCATCCAGTGTTTCTGTTCACCATATTg       >  1:2805669/1‑71 (MQ=255)
                                                               cGGTGTACCAATCTGCCCTAACATCGCCGCCAGCGTCACCAGGATCCAGTGTTTCTGTTCACCATATTGCt     <  1:3908737/71‑1 (MQ=255)
                                                               cGGTGTACCAATCTGCCCTAACATCGCCGCCAGCGTCACCAGCATCCAg                           >  1:1119392/1‑49 (MQ=255)
                                                                  tgtACCAATCTGCCCTAACATCGCCGCCAGCGTCACCAGCATCCAGTGTTTCTGTTCACCATATTGCTGgc  >  1:2667483/1‑71 (MQ=255)
                                                                    tACCAATCTGCCCTAACATCGCCGCCAGCGTCACCAGCATCCAGTGTTTCTg                   <  1:1224449/52‑1 (MQ=255)
                                                                    tACCAATCTGCCCTAACATCGCCGCCAGCGTCACCAGCATCCAGTGTTTCTg                   <  1:3432571/52‑1 (MQ=255)
                                                                     |                                                                   
GCCGCCAACACGCGTCGGGTTGCCACCGTTGGAGTAATGATAAGAGAAACCAAAACCGCCGCCCGGTGTACCAATCTGCCCTAACATCGCCGCCAGCGTCACCAGCATCCAGTGTTTCTGTTCACCATATTGCTGGC  >  Exported/1259718‑1259854

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