Predicted mutation
evidence seq id position mutation freq annotation gene description
RA Exported 1,259,757 A→G 14.0% Y389Y (TAT→TAC 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,7570AG14.0% 117.8 / 10.3 50Y389Y (TAT→TACtorZtrimethylamine N‑oxide reductase system III, catalytic subunit
Reads supporting (aligned to +/- strand):  ref base A (17/26);  new base G (4/3);  total (21/29)
Fisher's exact test for biased strand distribution p-value = 4.34e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

GCAATCGCCGCAGACATTTCCGGCAGCACGCCGCCAACACGCGTCGGGTTGCCACCGTTGGAGTAATGATAAGAGAAACCAAAACCGCCGCCCGGTGTACCAATCTGCCCTAACATCGCCGCCAGCGTCACCA  >  Exported/1259689‑1259821
                                                                    |                                                                
gCAATCGCCGCAGACATTTCCGGCAGCACGCCGCCAACACGCGTCGGGTTGCCACCGTTGGAGTAATGAt                                                                 <  1:923526/70‑1 (MQ=255)
gCAATCGCCGCAGACATTTCCGGCAGCACGCCGCCAACACGCGTCGGGTTGCCACCGTTGGAGTAATGAt                                                                 <  1:916944/70‑1 (MQ=255)
  aaTCGCCGCAGACATTTCCGGCAGCACGCCGCCAACACGCGTCGGGTTGCCACCGTTGGAGTAATGAt                                                                 >  1:5641937/1‑68 (MQ=255)
  aaTCGCCGCAGACATTTCCGGCAGCACGCCGCCAACACGCGTCGGGTTGCCACCGTTGGAGTAATGAt                                                                 >  1:321698/1‑68 (MQ=255)
  aaTCGCCGCAGACATTTCCGGCAGCACGCCGCCAACACGCGTCGGGTTGCCACCGTAGGAGTAATGATAag                                                              <  1:4812935/71‑1 (MQ=255)
  aaTCGCCGCAGACATTTCCGGCAGCAAGCCGCCAACACGCGTCGGGTTGCCACCGTTGGAGTAATGAt                                                                 >  1:381636/1‑68 (MQ=255)
   aTCGCCGCAGACATTTCCGGCAGCACGCCGCCAACACGCGTCGGGTTGCCACCGTTGGAGTAATGATAaga                                                             <  1:3590063/71‑1 (MQ=255)
   aTCGCCGCAGACATTTCCGGCAGCACGCCGCCAACACGCGTCGGGTTGCCACCGTTGGAGTAATGATAaga                                                             <  1:318121/71‑1 (MQ=255)
    tCGCCGCAGACATTTCCGGCAGCACGCCGCCAACACGCGTCGGGTTGCCACCGTTGGAGTAATGATAagag                                                            <  1:1018400/71‑1 (MQ=255)
     cgccgcAGACATTTCCGGCAGCACGCCGCCAACACGCGTCGGGTTGCCACCGTTGGAGTAATGATAag                                                              >  1:1223154/1‑68 (MQ=255)
      gccgcAGACATTTCCGGCAGCACGCCGCCAACACGCGTCGGGTTGCCACCGTTGGAGTAATGATAAgaga                                                           >  1:1541340/1‑70 (MQ=255)
          cAGACATTTCCGGCAGCACGCCGCCAACACGCGTCGGGTTGCCACCGTTGGAGTAATGATAAGAGAAACCa                                                      <  1:1028929/71‑1 (MQ=255)
           aGACATTTCCGGCAGCACGCCGCCAACACGCGTCGGGTTGCCACCGTTGGAGTAATGATAAGAGAAACCaa                                                     <  1:1502292/71‑1 (MQ=255)
           aGACATTTCCGGCAGCACGCCGCCAACACGCGTCGGGTTGCCACCGTTGGAGTAATGATAAGAGAAACCaa                                                     <  1:4042393/71‑1 (MQ=255)
           aGACATTTCCGGCAGCACGCCGCCAACACGCGTCGGGTTGCCACCGTTGGAGTAATGATAAGAGAAACCaa                                                     <  1:5161683/71‑1 (MQ=255)
           aGACATTTCCGGCAGCACGCCGCCAACACGCGTCGGGTTGCCACCGTTGGAGTAATGATAAGAGAAACCaa                                                     <  1:2435816/71‑1 (MQ=255)
           aGACATTTCCGGCAGCACGCCGCCAACACGCGTCGGGTTGCCACCGTTGGAGTAATGATAAGAGAAACCaa                                                     <  1:1930038/71‑1 (MQ=255)
           aGACATTTCCGGCAGCACGCCGCCAACACGCGTCGGGTTGCCACCGTTGGAGTAATGATAAGAGAAACCaa                                                     <  1:1665503/71‑1 (MQ=255)
             aCATTTCCGGCAGCACGCCGCCAACACGCGTCGGTTTGCCACCGTTGGAGTAATGATAAGAGAAACCaaaa                                                   <  1:5297226/71‑1 (MQ=255)
             aCATTTCCGGCAGCACGCCGCCAACACGCGTCGGGTTGCCACCGTTGGAGTAATGATAAGAGAAAcc                                                       <  1:2691006/67‑1 (MQ=255)
              cATTTCCGGCAGCACGCCGCCAACACGCGTCGGGTTGCCACCGTTGGAGTAATGATAAGAGAAACCAAAAc                                                  <  1:5537514/71‑1 (MQ=255)
                 ttCAGGCAGCACGCCGCCAACACGCGTCGGGTTGCCACCGTTGGAGTAATGATAAgaga                                                           >  1:135840/1‑59 (MQ=255)
                         gcaCGCCGCCAACACGCGTCGGGTTGCCACCGTTGGAGTAATGATAAGAGAAACCAAAACCGCCGCCCGgt                                       >  1:3502688/1‑71 (MQ=255)
                           aCGCCGCCAACACGCGTCGGGTTGCCACCGTTGGAGTAATGa                                                                  >  1:3587890/1‑42 (MQ=255)
                           aCGCCGCCAACACGCGTCGGGTTGCCACCGTTGGAGTAATGa                                                                  >  1:3037138/1‑42 (MQ=255)
                           aCGCCGCCAACACGCGTCGGGTTGCCACCGTTGGAGTAATGa                                                                  >  1:3310084/1‑42 (MQ=255)
                             gccgccAACACGCGTCGGGTTGCCACCGTTGGAGTAATGATAAGAGAAACCAAAACCGCCGCCCGGTGTAc                                   >  1:5665775/1‑71 (MQ=255)
                              ccgccAACACGCGTCGGGTTGCCACCGTTGGAGTAATGATAAGAGAAACCAAAACCGCCGCCCGGTGTAcc                                  <  1:1749692/71‑1 (MQ=255)
                               cgccAACACGCGTCGGGTTGCCACCGTTGGAGTAATGATAAgaga                                                           >  1:1222053/1‑45 (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)
                                      acGCGTCGGGTTGCCACCGTTGGAGTAATGATAAgaga                                                           >  1:4261957/1‑38 (MQ=255)
                                      acGCGTCGGGTTGCCACCGTTGGAGTAATGATAAgaga                                                           >  1:4913034/1‑38 (MQ=255)
                                       cgcgTCGGGTTGCCACCGTTGGAGTAATGAAAAGAGAAACCAAAAccgccg                                             >  1:2359074/1‑51 (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)
                                                                    |                                                                
GCAATCGCCGCAGACATTTCCGGCAGCACGCCGCCAACACGCGTCGGGTTGCCACCGTTGGAGTAATGATAAGAGAAACCAAAACCGCCGCCCGGTGTACCAATCTGCCCTAACATCGCCGCCAGCGTCACCA  >  Exported/1259689‑1259821

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: