Predicted mutation
evidence seq id position mutation freq annotation gene description
RA Exported 1,259,790 A→G 27.4% I378I (ATT→ATC 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,7900AG27.4% 73.5 / 31.1 51I378I (ATT→ATCtorZtrimethylamine N‑oxide reductase system III, catalytic subunit
Reads supporting (aligned to +/- strand):  ref base A (18/19);  new base G (10/4);  total (28/23)
Fisher's exact test for biased strand distribution p-value = 2.10e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.42e-01

CGCCAACACGCGTCGGGTTGCCACCGTTGGAGTAATGATAAGAGAAACCAAAACCGCCGCCCGGTGTACCAATCTGCCCTAACATCGCCGCCAGCGTCACCAGCATCCAGTGTTTCTGTTCACCATATTGCTGGC  >  Exported/1259720‑1259854
                                                                      |                                                                
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)
                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)
                            ggAGTAATTATAAGAGAAACCAAAACCGCCGCCCGGTGTACCa                                                                  <  1:4634668/43‑1 (MQ=255)
                             gAGTAATGATAAGAGAAACCAAAACCGCCGCCCGGTGTACCAATCTGCCCTAACATCGCCGCCAGCGTCAc                                     <  1:3032831/71‑1 (MQ=255)
                             gAGTAATGATAAGAGAAACCAAAACCGCCGCCCGGTGTACCAATCTGCCCTAACATCGCCGCCAGCGTCAc                                     <  1:2844667/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:4525914/70‑1 (MQ=255)
                                     gtAAGAAAGACCAAAACCGCCGCCGGGTGTGCCGATTTGCCCCAACATTGCTGCCAGCGTGACGATCATcc                             <  1:138006/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)
                                                                    ccAATCTGCCCTAACATCGCCGCCAGCGTCACCAGCATCCAgtgt                        >  1:2318236/1‑45 (MQ=255)
                                                                    ccAATCTGCCCTAACATCGCCGCCAGCGTCACCAGCATCCAgtgt                        >  1:1547495/1‑45 (MQ=255)
                                                                    ccAATCTGCCCTAACATCGCCGCCAGCGTCACCAGCATCCAgtgt                        >  1:2839865/1‑45 (MQ=255)
                                                                      |                                                                
CGCCAACACGCGTCGGGTTGCCACCGTTGGAGTAATGATAAGAGAAACCAAAACCGCCGCCCGGTGTACCAATCTGCCCTAACATCGCCGCCAGCGTCACCAGCATCCAGTGTTTCTGTTCACCATATTGCTGGC  >  Exported/1259720‑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: