Predicted mutation
evidence seq id position mutation freq annotation gene description
RA Exported 1,260,356 A→C 8.3% S190A (TCT→GCT) ‡ 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,260,3560AC 8.3% 204.9 / 11.1 84S190A (TCT→GCT) ‡torZtrimethylamine N‑oxide reductase system III, catalytic subunit
Reads supporting (aligned to +/- strand):  ref base A (40/37);  new base C (3/4);  total (43/41)
Fisher's exact test for biased strand distribution p-value = 7.10e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

GCACCACAACCTGGCTGTTTTCCAGAATCAGCGGCCAGGAGGTCTGCTGTTCATACACCTCTACAGAACCGACCACATGCGGCATGATCACCTGTGCTGCCCCGGTTGAATAATCGCCGCTATGCCCGGAATAA  >  Exported/1260290‑1260423
                                                                  |                                                                   
gCACCACAACCTGGCTGTTTTCCAGAATCAGCGGCCAGGAGGTCTGCTGTTCATACACCTCTACAGAAcc                                                                  >  1:4129664/1‑70 (MQ=255)
    cacaACCTGGCTGTTTTCCAGAATCAGCGGCCAGGAGGTCTGCTGTTCATACACCTCTACAGAACCGACca                                                             <  1:657866/71‑1 (MQ=255)
     acaACCTGGCTGTTTTCCAGAATCAGCGGCCAGGAGGTCTGCTGTTCATACACCTCTACAGAACCGACcac                                                            <  1:1199900/71‑1 (MQ=255)
      caACCTGGCTGTTTTCCAGAATCAGCGGCCAGGAGGTCTGCTGTTCATACACCTCTACAGAACCGACcaca                                                           >  1:2621935/1‑71 (MQ=255)
       aaCCTGGCTGTTTTCCAGAATCAGCGGCCAGGAGGTCTGCTGTTCATACACCTCTACAGAACCGACCACAt                                                          >  1:745218/1‑71 (MQ=255)
        aCCTGGCTGTTTTCCAGAATCAGCGGCCAGGAGGTCTGCTGTTCATACACCTCTACAGAACCGACcaca                                                           <  1:5497130/69‑1 (MQ=255)
        aCCTGGCTGTTTTCCAGAATCAGCGGCCAGGAGGTCTGCTGTTCATACACCTCTACAGAACCGACCACATg                                                         <  1:1657138/71‑1 (MQ=255)
        aCCTGGCTGTTTTCCAGAATCAGCGGCCAGGAGGTCTGCTGTTCATACACCTCTACAGAACCGACCACATg                                                         >  1:1607392/1‑71 (MQ=255)
          cTGGCTGTTTTCCAGAATCAGCGGCCAGGAGGTCTGCTGTTCATACACCTCTACAGAACCGACCACATGCg                                                       <  1:1819726/71‑1 (MQ=255)
            ggCTGTTTTCCAGAATCAGCGGCCAGGAGGTCTGCTGTTCATACACCTCTACAGAACCGACCACATGCgg                                                      >  1:5348665/1‑70 (MQ=255)
                 ttttCCAGAATCAGCGGCCAGGAGGTCTGCTGTTCATACACCTCTACAGAAcc                                                                  <  1:380373/53‑1 (MQ=255)
                 ttttCCAGAATCAGCGGCCAGGAGGTCTGCTGTTCATACACCTCTACAGAAcc                                                                  <  1:3638016/53‑1 (MQ=255)
                 ttttCCAGAATCAGCGGCCAGGAGGTCTGCTGTTCATACACCTCTACAGAACCGACcaca                                                           >  1:2118921/1‑60 (MQ=255)
                 ttttCCAGAATCAGCGGCCAGGAGGTCTGCTGTTCATACACCTCTACAGAACCGACcaca                                                           >  1:2125503/1‑60 (MQ=255)
                 ttttCCAGAATCAGCGGCCAGGAGGTCTGCTGTTCATACACCTCTACAGAACCGACCACATGCGGCa                                                    <  1:4982701/67‑1 (MQ=255)
                 ttttCCAGAATCAGCGGCCAGGAGGTCTGCTGTTCATACACCTCTACAGAACCGACCACATGCGGCa                                                    <  1:3013841/67‑1 (MQ=255)
                 ttttCCAGAATCAGCGGCCAGGAGGTCTGCTGTTCATACACCTCTACAGAACCGACCACATGCGGCATGa                                                 >  1:2053805/1‑70 (MQ=255)
                 tgttCCAGCACCAGCGGCCAACTGGTCTGCTGTTGATAAACTTCACTACCACCCACGACATACGGCATGa                                                 <  1:3837518/68‑1 (MQ=255)
                   ttCCAGCACCAGCGGCCAACTGGTCTGCTGTTGATAAACTTCACTACCACCCACGACATACGGCATGAt                                                >  1:3448164/1‑69 (MQ=255)
                   ttCCAGCACCAGCGGCCAACTGGTCTGCTGTTGATAAACTTCACTACCACCCACGACATACGGCATGATcg                                              <  1:5635217/71‑2 (MQ=255)
                   ttCCAGCACCAGCGGCCAACTGGTCTGCTGTTGATAAACTTCACTACCACCCACGACATACGGCATGATcg                                              <  1:3771724/71‑2 (MQ=255)
                   ttCCAGCACCAGCGGCCAACTGGTCTGCTGTTGATAAACTTCACTACCACCCACGACATACGGCATGATcg                                              >  1:4888777/1‑70 (MQ=255)
                   ttCCAGAATCAGCGGCCAGGAGGTCTGCTGTTCATACACCTCTACAGAACCGACcaca                                                           >  1:5386442/1‑58 (MQ=255)
                   ttCCAGAATCAGCGGCCAGGAGGTCTGCTGTTCATACACCTCTACAGAACCGACcac                                                            >  1:2792368/1‑57 (MQ=255)
                   ttCCAGAATCAGCGGCCAGGAGGTCTGCTGTTCATACACCTCTACAGAACCGACCACATGCGGCATg                                                  >  1:5605029/1‑67 (MQ=255)
                   ttCCAGAATCAGCGGCCAGGAGGTCTGCTGTTCATACACCTCTACAGAACCGACCACATGCGGCATg                                                  >  1:3667326/1‑67 (MQ=255)
                   ttCCAGAATCAGCGGCCAGGAGGTCTGCAGTTCATACACCTCTACAGAACCGACCACATGCGGCATg                                                  >  1:1734160/1‑67 (MQ=255)
                     ccAGAATCAGCGGCCAGGAGGTCTGCTGTTCATACACCTCTACAGAACCGACCACATGCGGCATGATc                                               <  1:4622432/68‑1 (MQ=255)
                     ccAGAATCAGCGGCCAGGAGGTCTGCTGTTCATACACCTCTACAGAACCGACCACATGCGGCATGATc                                               <  1:1905987/68‑1 (MQ=255)
                     ccAGAATCAGCGGCCAGGAGGTCTGCTGTTCATACACCTCTACAGAACCGACCACATGCGGCATGATCAcc                                            >  1:2463294/1‑71 (MQ=255)
                     ccAGAATCAGCGGCCAGGAGGTCTGCTGTTCATACACCTCTACAGAACCGACCACATGCGGCATGATCAc                                             <  1:65811/70‑1 (MQ=255)
                     ccAGAATCAGCGGCCAGGAGGTCTGCTGTTCATACACCTCTACAGAACCGACCACATGCGGCATGATCAc                                             <  1:4802583/70‑1 (MQ=255)
                     ccAGAATCAGCGGCCAGGAGGTCTGCTGTTCATACACCTCTACAGAACCGACCACATGCGGCATGATCAc                                             >  1:1843578/1‑70 (MQ=255)
                     ccAGAATCAGCGGCCAGGAGGTCTGCTGTTCATACACCTCTACAGAACCGACCACATGCGGCATGATCAc                                             <  1:3754421/70‑1 (MQ=255)
                       aGAATCAGCGGCCAGGAGGTCTGCTGTTCATACACCTCTACAGa                                                                     <  1:2130681/44‑1 (MQ=255)
                         aaTCAGCGGCCAGGAGGTCTGCTGTTCATACACCTCTACAGAACCGACCACATGCGGCATGATCACCtgtg                                        >  1:2739886/1‑71 (MQ=255)
                         aaTCAGCGGCCAGGAGGTCTGCTGTTCATACACCTCTACAGAACCGACCACATGCGGCATGATCACCtgtg                                        >  1:636680/1‑71 (MQ=255)
                           ccAGCGGCCAGGAGGTTTGCTGTTCATACACTTCCATCGAACCGACTACGCGCGGCAGGATCAc                                             >  1:4428977/2‑64 (MQ=255)
                           ccAGCGGCCAGGAGGTTTGCTGTTCATACACTTCCATCGAACCGACTACGCGCGGCAGGATCAc                                             >  1:2469112/2‑64 (MQ=255)
                           ccAGCGGCCAGGAGGTTTGCTGTTCATACACTTCCATCGAACCGACTACGCGCGGCAGGATCACCtgcg                                        >  1:2004066/2‑67 (MQ=255)
                               cGGCCAGGAGGTCTGCTGTTCATACACCTCTACAGAACCGACCACATGCGGCATGATCACCtgt                                         <  1:1127941/64‑1 (MQ=255)
                               cGGCCAGGAGGTCTGCTGTTCATACACCTCTACAGAACCGACCACATGCGGCATGATCACCtgt                                         <  1:2812304/64‑1 (MQ=255)
                               cGGCCAGGAGGTCTGCTGTTCATACACCTCTACAGAACCGACCACATGCGGCATGATCACCTGTGCTGccc                                  <  1:3558249/71‑1 (MQ=255)
                                 gCCAGGAGGTCTGCTGTTCATACACCTCTACAGAACCGACCACATGCGGCATGATCACCTGTGCTGCCCCg                                >  1:2095323/1‑71 (MQ=255)
                                 gCCAGGAGGTCTGCTGTTCATACACCTCTACAGAACCGACCACATGCGGCATGATCACCTGTGCTGCCCCg                                >  1:3156104/1‑71 (MQ=255)
                                   cAGGAGGTTTGCTGTTCATACACTTCCATCGAACCGACTACGCGCGGCAGGATCACCAGCGCAGCACCgg                               >  1:1902837/1‑70 (MQ=255)
                                    aggaggTCTGCTGTTCATACACCTCTACAGAACCGACCACATGCGGCATGATCACCTGTGCTGCCCCGGtt                             >  1:2762715/1‑71 (MQ=255)
                                      gaggTCTGCTGTTCATACACCTCTACAGAACCGACCACATGCGGCATGATCAc                                             >  1:2985231/1‑53 (MQ=255)
                                       tggTCTGCTGTTGATAAACTTCACTACCACCCACGACATACGGCATGAt                                                >  1:4584112/2‑49 (MQ=255)
                                        ggTTTGCTGTTCATACACTTCCATCGAACCGACTAAGCGCGGCAGGATCACCTGCGCAGCACAGGTAGag                          >  1:866146/1‑69 (MQ=255)
                                          tCTGCTGTTCATACACCTCTACAGAACCGACCACATGCGGCATGATCACCTGTGCTGCCCCGGTTGaataa                       <  1:3517141/71‑1 (MQ=255)
                                           ctgctgTTCATACACCTCTACAGAACCGACCACATGCGGCATGATCACCTGTGCTGCCCCGGTTGaataat                      >  1:2371063/1‑71 (MQ=255)
                                           ctgctgTTCATACACCTCTACAGAACCGACCACATGCGGCATGATCACCTGTGCTGCCCCGGTTGaataat                      >  1:5492888/1‑71 (MQ=255)
                                            tgctgTTCATACACTTCCATCGAACCGACTACGCGCGGCAGGATCACCTGCGCAGCACCGGTAGagtaat                      <  1:5009371/70‑1 (MQ=255)
                                            tgctgTTCATACACCTCTACAGAACCGACCACATGCGGCATGATc                                               <  1:2666496/45‑1 (MQ=255)
                                            tgctgTTCATACACCTCTACAGAACCGACCACATGCGGCATGATCACCTGTGCTGCCCCGGTTGAATAATc                     >  1:3645477/1‑71 (MQ=255)
                                               tgTTCATACACCTCTACAGAACCGACCACATGCgg                                                      <  1:702544/35‑1 (MQ=255)
                                               tgTTCATACACCTCTACAGAACCGACCACATGCGGCATGATCACCTGTGCTGCCCCGGTTGAATAATcgcc                  >  1:181128/1‑71 (MQ=255)
                                                 ttCATACACCTCTACAGCACCGACCACATGCGGCATGATCACCTGTGCTGCCCCGGTTGaat                         <  1:4559352/62‑1 (MQ=255)
                                                 ttCATACACCTCTACAGAACCGACCACATGCGGCATGATCACCTGTGCTGCCCCGGTTGaat                         <  1:3231834/62‑1 (MQ=255)
                                                 ttCATACACCTCTACAGAACCGACCACATGCGGCATGATCACCTGTGCTGCCCCGGTTGAATAATcgccgc                <  1:1894208/71‑1 (MQ=255)
                                                   cATACACCTCTACAGAACCGACCACATGCGGCATGATCACCTGTGCTGCCCCGGTTGaataa                       >  1:440420/1‑62 (MQ=255)
                                                   cATACACCTCTACAGAACCGACCACATGCGGCATGATCACCTGTGCTGCCCCGGTTGaataa                       >  1:4765835/1‑62 (MQ=255)
                                                     tACACCTCTACAGAACCGACCACATGCGGCATGATCa                                              <  1:1579972/37‑1 (MQ=255)
                                                     tACACCTCTACAGAACCGACCACATGCGGCATGATCACCTGtgctgc                                    <  1:4105586/47‑1 (MQ=255)
                                                       cacCTCTACAGAACCGACCACATGCGGCATGATCACCTGTGCTGccc                                  >  1:1769775/1‑47 (MQ=255)
                                                         ccTCTACAGAACCGACCACATGCGGCATGATCACCTGTGCTGCCCCGGTTGAATAATCGCCGCTATGCCCg        <  1:2810837/71‑1 (MQ=255)
                                                         ccTCTACAGAACCGACCACATGCGGCATGATCACCTGTGCTGCCCCGGTTGAATAATCGCCGCTATGCCCg        <  1:5546190/71‑1 (MQ=255)
                                                          ctctACAGAACCGACCACATGCGGCATGATCACCTGTGCTGCCCCGGTTGAATAATCGCCGCTATGCCCg        >  1:2166735/1‑70 (MQ=255)
                                                           tctACAGAACCGACCACATGCGGCATGATCACCtg                                          <  1:5660438/35‑1 (MQ=255)
                                                            ctACAGAACCGACCACATGCGGCATGATCACCTGTGCTGCCCCGGTTGAATAATCGCCGCTATGCCCg        <  1:949837/68‑1 (MQ=255)
                                                            ctACAGAACCGACCACATGCGGCATGATCACCTGTGCTGCCCCGGTTGAATAATCGCCGCTATGCCCg        <  1:3560212/68‑1 (MQ=255)
                                                            ctACAGAACCGACCACATGCGGCATGATCACCTGTGCTGCCCCGGTTGAATAATCGCCGCTATGCCCGGaa     >  1:4810591/1‑71 (MQ=255)
                                                             tACAGAACCGACCACATGCGGCATGATCACCtgtg                                        >  1:167511/1‑35 (MQ=255)
                                                              actgAACCGACCACATGCGGCATGATCACCTGTGCTGCCCCGGTTGAATAATCGCCGCTATGCCCGGAAt    >  1:4850669/4‑70 (MQ=255)
                                                               cAGAACCGACCACATGCGGCATGATCACCTGTGCTGCCCCGGTTGAATAATCGCCGCTATGCCCGGAAt    <  1:2134174/69‑1 (MQ=255)
                                                               cAGAACCGACCACATGCGGCATGATCACCTGTGCTGCCCCGGTTGAATAATCGCCGCTATGCCCGGAAt    <  1:1682495/69‑1 (MQ=255)
                                                               cAGAACCGACCACATGCGGCATGATCACCTGTGCTGCCCCGGTTGAATAATCGCCGCTATGCCCGGAATaa  >  1:5338316/1‑71 (MQ=255)
                                                               cAGAACCGACCACATGCGGCATGATCACCTGTGCTGCCCCGGTTGAATAATCGCCGCTATGCCCGGAATaa  <  1:4300170/71‑1 (MQ=255)
                                                               cAGAACCGACCACATGCGGCATGATCACCTGTGCTGCCCCGGTTGAATAATCGCCGCTATGCCCGGAATa   <  1:4662064/70‑1 (MQ=255)
                                                               cAGAACCGACCACATGCGGCATGATCACCTGTGCTGCCCCGGTTGAATAATCGCCGCTATGCCCGGAATa   <  1:1919190/70‑1 (MQ=255)
                                                               cAGAACCGACCACATGCGGCATGATCACCTGTGCTGCCCCGGTTGAATAATCGCCGCTATGCCCGGAATa   <  1:1350277/70‑1 (MQ=255)
                                                                aGAACCGACCACATGCGGCATGATCACCTGTGCTGCCCCGGTTGAATAATcgccgc                <  1:3756124/56‑1 (MQ=255)
                                                                 gAACCGACCACATGCGGCATGATCACCTGTGCTGCCCCGGTTGaat                         >  1:5011626/1‑46 (MQ=255)
                                                                 gAACCGACCACATGCGGCATGATCACCTGTGCTGCCCCGGTTGAATAATCGCCGCTATGCCCGGAATaa  >  1:5476154/1‑69 (MQ=255)
                                                                  |                                                                   
GCACCACAACCTGGCTGTTTTCCAGAATCAGCGGCCAGGAGGTCTGCTGTTCATACACCTCTACAGAACCGACCACATGCGGCATGATCACCTGTGCTGCCCCGGTTGAATAATCGCCGCTATGCCCGGAATAA  >  Exported/1260290‑1260423

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