Predicted mutation
evidence seq id position mutation freq annotation gene description
RA Exported 1,259,838 T→C 18.7% 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
*Exported1,259,8380TC18.7% 191.8 / 36.0 91E362E (GAA→GAGtorZtrimethylamine N‑oxide reductase system III, catalytic subunit
Reads supporting (aligned to +/- strand):  ref base T (47/27);  new base C (10/7);  total (57/34)
Fisher's exact test for biased strand distribution p-value = 7.84e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 8.97e-01

AACCGCCGCCCGGTGTACCAATCTGCCCTAACATCGCCGCCAGCGTCACCAGCATCCAGTGTTTCTGTTCACCATATTGCTGGCGCTGAATTCCCCAGCCTGCCATCAGCATTGTGCGATTCGCCGCCATCAGTTCAG  >  Exported/1259771‑1259908
                                                                   |                                                                      
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)
       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)
             tgtACCAATCTGCCCTAACATCGCCGCCAGCGTCACCAGCATCCAGTGTTTCTGTTCACCATATTGCTGgc                                                        >  1:2667483/1‑71 (MQ=255)
                     tttGCCCCAACATTGCTGCCAGCGTGACGATCATCCAGTGTTTTTGCTCACCAAACTGTTGGCGCTGCAt                                                 >  1:4054543/3‑70 (MQ=255)
                     tttGCCCCAACATTGCTGCCAGCGTGACGATCATCCAGTGTTTTTGCTCACCAAACTGTTGGCGCTGCAt                                                 >  1:1989388/3‑70 (MQ=255)
                     tCTGCCCTAACATCGCCGCCAGCGTCACCAGCATCCAGTGTTTCTGTTCACCATATTGCt                                                           <  1:148507/60‑1 (MQ=255)
                      cTGCCCTAACATCGCCGCCAGCGTCACCAGCATCCAGTGTTTCTGTTCACCATATTGCTGGCGCTGAATTc                                               >  1:761480/1‑71 (MQ=255)
                      cTGCCCTAACATCGCCGCCAGCGTCACCAGCATCCAGTGTTTCTGTTCACCATATTGCTGGCGCTGAATTc                                               >  1:3686733/1‑71 (MQ=255)
                      cTGCCCTAACATCGCCGCCAGCGTCACCAGCATCCAGTGTTTCTGTTCACCATATTGCTGGCGCTGAATTc                                               >  1:4127653/1‑71 (MQ=255)
                       tGCCCTAACATCGCCGCCAGCGTCACCAGCATCCAGTGTTTCTGtt                                                                       >  1:4328969/1‑46 (MQ=255)
                       tGCCCTAACATCGCCGCCAGCGTCACCAGCATCCAGTGTTTCTGtt                                                                       >  1:402258/1‑46 (MQ=255)
                       tGCCCTAACATCGCCGCCAGCGTCACCAGCATCCAGTGTTTCTGtt                                                                       >  1:3432523/1‑46 (MQ=255)
                       tGCCCTAACATCGCCGCCAGCGTCACCAGCATCCAGTGTTTCTGTTCACCATATTGCTGGCGCTGAATTcc                                              >  1:5063077/1‑71 (MQ=255)
                         cccTAACATCGCCGCCAGCGTCACCAGCATCCAGTGTTTCTGTTCACCATATTGCTGGCGCTGAATTcccc                                            <  1:4124712/71‑1 (MQ=255)
                         cccTAACATCGCCGCCAGCGTCACCAGCATCCAGTGTTTCTGTTCACCATATTGCTGGCGCTGAATTcccc                                            <  1:2038243/71‑1 (MQ=255)
                          cccaaCATTGCTGCCAGCGTGACGATCATCCAGTGTTTTTGCTCACCAAACTGTTGGCGCTGCATTccc                                             >  1:1510680/4‑69 (MQ=255)
                          cccaaCATTGCTGCCAGCGTGACGATCATCCAGTGTTTTTGCTCACCAAACTGTTGGCGCTGCATTccc                                             >  1:5694394/4‑69 (MQ=255)
                          cccaaCATTGCTGCCAGCGTGACGATCATCCAGTGTTTTTGCTCACCAAACTGTTGGCGCTGCATTccc                                             >  1:175182/4‑69 (MQ=255)
                          cccaaCATTGCTGCCAGCGTGACGATCATCCAGTGTTTTTGCTCAACAAACTGTTGGCGCTGCATTCCCCa                                           <  1:5584119/68‑1 (MQ=255)
                          ccTAACATCGCCGCCAGCGTCACCAGCATCCAGTGTTTCTGTTCACCATATTGCTGGCGCTGAATTcccc                                            <  1:3988489/70‑1 (MQ=255)
                            tAACATCGCCGCCAGCGTCACCAGCATCCAGTGTTTCTGTTCAc                                                                    <  1:1748286/44‑1 (MQ=255)
                            tAACATCGCCGCCAGCGTCACCAGCATCCAGTGTTTCTGTTCACCatat                                                               <  1:3359860/49‑1 (MQ=255)
                            tAACATCGCCGCCAGCGTCACCAGCATCCAGTGTTTCTGTTCACCATATTGCTGGCGCt                                                     <  1:1665276/59‑1 (MQ=255)
                            tAACATCGCCGCCAGCGTCACCAGCATCCAGTGTTTCTGTTCACCATATTGCTGGCGCTGAATTccc                                             >  1:3949943/1‑67 (MQ=255)
                            tAACATCGCCGCCAGCGTCACCAGCATCCAGTGTTTCTGTTCACCATATTGCTGGCGCTGAATTCCCCAGc                                         >  1:3334757/1‑71 (MQ=255)
                             aaCATCGCCGCCAGCGTCACCAGCATCCAGTGTTTCTGTTCACCATATTGCTGGCGCTGAATTCCCCAGcc                                        >  1:1835596/1‑71 (MQ=255)
                              aCATTGCTGCCAGCGTGACGATCATCCAGTGTTTTTGCTCACCAAACTGTTgg                                                         >  1:4864383/1‑53 (MQ=255)
                              aCATTGCTGCCAGCGTGACGATCATCCAGTGTTTTTGCTCACCAAACTGTTgg                                                         >  1:205362/1‑53 (MQ=255)
                              aCATCGCCGCCAGCGTCACCAGCATCCAGTGTTTCTGTTCACCata                                                                >  1:827036/1‑46 (MQ=255)
                              aCATCGCCGCCAGCGTCACCAGCATCCAGTGTTTCTGTTCACCata                                                                >  1:3982914/1‑46 (MQ=255)
                              aCATCGCCGCCAGCGTCACCAGCATCCAGTGTTTCTGTTCACCata                                                                >  1:5096688/1‑46 (MQ=255)
                              aCATCGCCGCCAGCGTCACCAGCATCCAGTGTTTCTGTTCACCa                                                                  <  1:2800889/44‑1 (MQ=255)
                              aCATCGCCGCCAGCGTCACCAGCATCCAGTGTTTCTGTTCACCa                                                                  <  1:1991530/44‑1 (MQ=255)
                              aCATCGCCGCCAGCGTCACCAGCATCCAGTGTTTCTGTTCACCATATTg                                                             >  1:3103883/1‑49 (MQ=255)
                              aCATCGCCGCCAGCGTCACCAGCATCCAGTGTTTCTGTTCACCATATTg                                                             >  1:1661717/1‑49 (MQ=255)
                               cATCGCCGCCAGCGTCACCAGCATCCAGTGTTTCTGTTCACCATATTGCTGGCGCTGAATTCCCCAGCCTg                                      >  1:2768747/1‑71 (MQ=255)
                                aTCGCCGCCAGCGTCACCAGCATCCAGTGTTTCTGTTCACCa                                                                  <  1:3491325/42‑1 (MQ=255)
                                aTCGCCGCCAGCGTCACCAGCATCCAGTGTTTCTGTTCACCATATTGCTGGCGCTGAATTCCCCAGCCTGc                                     <  1:5090975/71‑1 (MQ=255)
                                  cgccgcCAGCGTCACCAGCATCCAGTGTTTCTGTTCACCata                                                                >  1:1491933/1‑42 (MQ=255)
                                  cgccgcCAGCGTCACCAGCATCCAGTGTTTCTGTTCACCata                                                                >  1:2764681/1‑42 (MQ=255)
                                  cgccgcCAGCGTCACCAGCATCCAGTGTTTCTGTTCACCata                                                                >  1:3638466/1‑42 (MQ=255)
                                  cgccgcCAGCGTCACCAGCATCCAGTGTTTCTGTTCACCata                                                                >  1:4720874/1‑42 (MQ=255)
                                  cgccgcCAGCGTCACCAGCATCCAGTGTTTCTGTTCACCata                                                                >  1:1536825/1‑42 (MQ=255)
                                   gctgccAGCGTGACGATCATCCAGTGTTTTTGCTCACCAAACTGTTGGCGCTGCATTCCCCAGCCTGCCAt                                  >  1:2796622/4‑71 (MQ=255)
                                   gccgccAGCGTCACCAGCATCCAGTGTTTCTGTTCACCa                                                                  >  1:2393787/1‑39 (MQ=255)
                                   gccgccAGCGTCACCAGCATCCAGTGTTTCTGTTCACCATATTGCTGGCGCTGAATTCCCCa                                           >  1:3956571/1‑62 (MQ=255)
                                    ctgccAGCGTGACGATCATCCAGTGTTTTTGCTCACCAAACTGTTGGCGCTGCATTCCCCAGCCTGCCATc                                 <  1:1108008/69‑1 (MQ=255)
                                    ccgccAGCGTCACCAGCATCCAGTGTTTCTGTTCACCATATTGCTGGCGCTGAATTCCCCAGCCTGCCATc                                 <  1:4935187/71‑1 (MQ=255)
                                     tgccAGCGTGACGATCATCCAGTGTTTTTGCTCACCAAACTGTTGGCGCTGCATTc                                               <  1:4673168/55‑1 (MQ=255)
                                      gccAGCGTCACCAGCATCCAGTGTTTCTGTTCACCATATTGCTGGCGCTGAATTCCCCAGCCTGc                                     >  1:2199633/1‑65 (MQ=255)
                                      gccAGCGTCACCAGCATCCAGTGTTTCTGTTCACCATATTGCTGGCGCTGAATTCCCCAGCCTGc                                     >  1:4741654/1‑65 (MQ=255)
                                      gccAGCGTCACCAGCATCCAGTGTTTCTGTTCACCATATTGCTGGCGCTGAATTCCCCAGCCTGCCATCAg                               <  1:2306961/71‑1 (MQ=255)
                                      gccAGCGTCACCAGCATCCAGTGTTTCTGTTAACCATATTGCTGGCGCTGAATTCCCCAGCCTGc                                     >  1:100003/1‑65 (MQ=255)
                                         aGCGTCACCAGCATCCAGTGTTTCTGTTCACCATATTGCTGGCGCTGAATTCCCCAGCCTGCCAt                                  >  1:557917/1‑65 (MQ=255)
                                         aGCGTCACCAGCATCCAGTGTTTCTGTTCACCATATTGCTGGCGCTGAATTCCCCAGCCTGCCAt                                  >  1:3026975/1‑65 (MQ=255)
                                          gCGTCACCAGCATCCAGTGTTTCTGTTCACCATATTGCTGGCGCTGAATTCCCCAGCCTGCCATCAg                               <  1:4246466/67‑1 (MQ=255)
                                           cGTCACCAGCATCCAGTGTTTCTGTTCACCATATTGCTGGCGCTGAATTCCCCAGCCTGCCATCAGCATtg                          <  1:3602227/71‑1 (MQ=255)
                                              cACCAGCATCCAGTGTTTCTGTTCACCATATTGCTGGCGCTGAATTCCCCAGCCTGCCATCAGCATtgtg                        >  1:88160/1‑70 (MQ=255)
                                               aCCAGCATCCAGTGTTTCTGTTCACCATATTGCTGGCGCTGAATTCCCCAGCCTGCCATCAGCATTGTGc                       >  1:5606113/1‑70 (MQ=255)
                                                ccAGCATCCAGTGTTTCTGTTCACCATATTGCTGGCGCTGAATTCCCCAGCCTGCCATCAGCATTGTGCGa                     <  1:1211039/71‑1 (MQ=255)
                                                 cAGCATCCAGTGTTTCTGTTCACCATATTGCTGGCGCTGAATTCCCCAGCCTGCCATCAGCATTGTGCGAt                    >  1:97902/1‑71 (MQ=255)
                                                 cAGCATCCAGTGTTTCTGTTCACCATATTGCTGGCGCTGAATTCCCCAGCCTGCCATCAGCATTGTGCGAt                    >  1:4539136/1‑71 (MQ=255)
                                                 cAGCATCCAGTGTTTCTGTTCACCATATTGCTGGCGCTGAATTCCCCAGCCTGCCATCAGCATTGTGCGAt                    >  1:2999096/1‑71 (MQ=255)
                                                  aGCATCCAGTGTTTCTGTTCACCATATTGCTGGCGCTGAATTCCCCAGCCTGCCATCAGCATTGTGCGAtt                   <  1:5489981/71‑1 (MQ=255)
                                                     aTCCAGTGTTTCTGTTCACCATATTGCTGGCGCTGAATTCCCCAGCCTGCCATCAGCATTGTGCGATTcg                 <  1:5386591/70‑1 (MQ=255)
                                                      tCCAGTGTTTCTGTTCACCATATTGCTGGCGCTGAATTCCCCAGCCTGCCATCAGCATTGTGCGAt                    <  1:5036926/66‑1 (MQ=255)
                                                        cAGTGTTTCTGTTCACCATATTGCTGGCGCTGAATTCCCCAGCCTGCCATCAGCATTGTGCGATTcgccg              <  1:3492910/70‑1 (MQ=255)
                                                         aGTGTTTCTGTTCACCATATTGCTGCCGCTGAATTCCCCAGCCTGCCATCAGCATTGTGCGATTcgccgc             >  1:2695104/1‑70 (MQ=255)
                                                          gtgtTTCTGTTCACCATATTGCTGGCGCTGAATTcc                                              <  1:1023609/36‑1 (MQ=255)
                                                          gtgtTTCTGTTCACCATATTGCTGGCGCTGAATTcc                                              <  1:4637152/36‑1 (MQ=255)
                                                          gtgtTTCTGTTCACCATATTGCTGGCGCTGAATTcc                                              <  1:510371/36‑1 (MQ=255)
                                                            gtTTCTGTTCACCATATTGCTGGCGCTGAATTCCCCAGCCTGCCATCAGCa                             >  1:4169970/1‑51 (MQ=255)
                                                             ttttTGCTCACCAAACTGTTGGGGCTGCATTCCCCAGCCTGCCATCa                                <  1:3389962/47‑1 (MQ=16)
                                                             ttttTGCTCACCAAACTGTTGGCGCTGCATTCCCCAGCCTGCCATCa                                <  1:3626513/47‑1 (MQ=19)
                                                             ttttTGCTCACCAAACTGTTGGCGCTGCATTCCCCAGCCTGCCATCa                                <  1:2860386/47‑1 (MQ=19)
                                                             ttttTGCTCACCAAACTGTTGGCGCTGCATTCCCCAGCCTGCCATCa                                <  1:5131128/47‑1 (MQ=19)
                                                             tttCTGTTCACCATATTGCTGGCGCTGAATTCCCCAGCCTGCCATCAGCATTGTGCGATTCGCCGCCAt          >  1:3645741/1‑69 (MQ=255)
                                                              ttCTGTTCACCATATTGCTGGCGCTGAATTCCCCAGCCTGCCAt                                  >  1:4300713/1‑44 (MQ=255)
                                                               tCTGTTCACCATATTGCTGGCGCTGAATTCCCCAGCCTGCCa                                   >  1:407781/1‑42 (MQ=255)
                                                               tCTGTTCACCATATTGCTGGCGCTGAATTCCCCAGCCTGCCa                                   >  1:5561233/1‑42 (MQ=255)
                                                               tCTGTTCACCATATTGCTGGCGCTGAATTCCCCAGCCTGCCa                                   >  1:3307494/1‑42 (MQ=255)
                                                               tCTGTTCACCATATTGCTGGCGCTGAATTCCCCAGCCTGCCa                                   >  1:3588402/1‑42 (MQ=255)
                                                               tCTGTTCACCATATTGCTGGCGCTGAATTCCCCAGCCTGCCATCAGCATtg                          >  1:937116/1‑51 (MQ=255)
                                                               tCTGTTCACCATATTGCTGGCGCTGAATTCCCCAGCCTGCCATCAGCATTGTGCGATTCGCCGCCATCAGt      >  1:1023915/1‑71 (MQ=255)
                                                                cTGTTCACCATATTGCTGGCGCTGAATTCCCCAGCCTGCCATCAGCATTGTGCGATTCGCCGCCATCAGt      <  1:325743/70‑1 (MQ=255)
                                                                cTGTTCACCATATTGCTGGCGCTGAATTCCCCAGCCTGCCATCAGCATTGTGCGATTCGCCGCCATCAGt      <  1:4283614/70‑1 (MQ=255)
                                                                 tGTTCACCATATTGCTGGCGCTGAATTCCCCAGCCTGCCATCAGCATTGTGCGAtt                   >  1:5607194/1‑56 (MQ=255)
                                                                 tGTTCACCATATTGCTGGCGCTGAATTCCCCAGCCTGCCATCAGCATTGTGCGATTCGCCGCCATCAGtt     >  1:3734576/1‑70 (MQ=255)
                                                                   ttCACCATATTGCTGGCGCTGAATTCCCCAGCCTGCCATCAGCATTGTGc                       <  1:4983520/50‑1 (MQ=255)
                                                                   ttCACCATATTGCTGGCGCTGAATTCCCCAGCCTGCCATCAGCATTGTGCGATTCGCCGCCATCAGTTCAg  >  1:5009449/1‑71 (MQ=255)
                                                                   |                                                                      
AACCGCCGCCCGGTGTACCAATCTGCCCTAACATCGCCGCCAGCGTCACCAGCATCCAGTGTTTCTGTTCACCATATTGCTGGCGCTGAATTCCCCAGCCTGCCATCAGCATTGTGCGATTCGCCGCCATCAGTTCAG  >  Exported/1259771‑1259908

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: