Predicted mutation
evidence seq id position mutation freq annotation gene description
RA Exported 1,259,604 G→A 5.8% G440G (GGC→GGT 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,6040GA 5.8% 257.1 / 3.7 85G440G (GGC→GGTtorZtrimethylamine N‑oxide reductase system III, catalytic subunit
Reads supporting (aligned to +/- strand):  ref base G (52/28);  new base A (3/2);  total (55/30)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 8.37e-01

CGCCACCCGCCCACCAGATCATTTTGATATTCGGATACGTTTGTTCTTTACCGTTATGCTGATATTTGCCGCCAGGATTTTCCAGCGCATCGACGATACGGGCAACAGGAATCGCCGTCATTCCGCCATCATCAGCGGCC  >  Exported/1259534‑1259673
                                                                      |                                                                     
cGCCACCCGCCCACCAGATCATTTTGATATTCGGATACGTTTGTTCTTTACCGTTATGCTGATATTTgccg                                                                       >  1:1185557/1‑71 (MQ=255)
     cccGCCCACCAGATCATTTTGATATTCGGATACGTTTGTTCTTTACCGTTATGCTGATATTTGCCGCCAgg                                                                  <  1:3596924/71‑1 (MQ=255)
      ccGCCAACCAGATCATTTTGATATTCGGATACGTTTGTTCTTTACCGTTATGCTGATATTTGCCGCCAgg                                                                  <  1:2294095/70‑1 (MQ=255)
       cGCCCACCGGATCATTTTGATATTCGGATACGTTTGTTCTTTACCGTTATGCTGATATTTGCCGCCAGGAt                                                                <  1:525012/71‑1 (MQ=255)
         cccACCAGATCATTTTGATATTCGGATACGTTTGTTCTTTACCGTTATGCTGATATTTGCCGCCAGGAttt                                                              >  1:5520131/1‑71 (MQ=255)
             ccaGATCATTTTGATATTCGGATACGTTTGTTCTTTACCGTTATGCTGATATTTGCCGCCAGGATTTTc                                                            >  1:3812701/1‑69 (MQ=255)
              caGATCATTTTGATATTCGGATACGTTTGTTCTTTACCGTTATGCTGATATTTGCCGCCAGGATTTTCCAg                                                         >  1:4086117/1‑71 (MQ=255)
               aGATCATTTTGATATTCGGATACGTTTGTTCTTTACCGTTATGCTGATATTTGCCGCCAGGATTTTCCAgc                                                        >  1:2913737/1‑71 (MQ=255)
               aGATCATTTTGATATTCGGATACGTTTGTTCTTTACCGTTATGCTGATATTTGCCGCCAGGATTTTCCAgc                                                        >  1:908143/1‑71 (MQ=255)
               aGATCATTTTGATATTCGGATACGTTTGTTCTTTACCGTTATGCTGATATTTGCCGCCAGGATTTTCCAgc                                                        >  1:287859/1‑71 (MQ=255)
               aGATCATTTTGATATTCGGATACGTTTGTTCTTTACCGTTATGCTGATATTTGCCGCCAGGATTTTCCAg                                                         >  1:1742802/1‑70 (MQ=255)
               aGATCATTTTGATATTCGGATACGTTTGTTCTTTACCGTTATGCTGATATTTGCCGCCAGGATTTTCCAg                                                         >  1:1473343/1‑70 (MQ=255)
                gATCATTTTGATATTCGGATACGTTTGTTCTTTACCGTTATGCTGATATTTGCCGCCAGGATTTTCCAgcg                                                       <  1:4043687/71‑1 (MQ=255)
                  tCATTTTGATATTCGGATACGTTTGTTCTTTACCGTTATGCTGATATTTGCCGCCAGGAtttt                                                             >  1:972954/1‑63 (MQ=255)
                   cATTTTGATATTCGGATACGTTTGTTCTTTACCGTTATGCTGATATTTGCCGCCAGGAtt                                                               >  1:2568848/1‑60 (MQ=255)
                   cATTTTGATATTCGGATACGTTTGTTCTTTACCGTTATGCTGATATTTGCCGCCAGGAtt                                                               >  1:874394/1‑60 (MQ=255)
                   cATTTTGATATTCGGATACGTTTGTTCTTTACCGTTATGCTGATATTTGCCGCCAGGAtt                                                               >  1:608410/1‑60 (MQ=255)
                   cATTTTGATATTCGGATACGTTTGTTCTTTACCGTTATGCTGATATTTGCCGCCAGGATTTTCCAGCGCAt                                                    >  1:3971528/1‑71 (MQ=255)
                   cATTTTGATATTCGGATACGTTTGTTCTTTACCGTTATGCTGATATTTGCCGCCAGGATTTTCCAGCGCAt                                                    >  1:2478823/1‑71 (MQ=255)
                   cATTTTGATATTCGGATACGTTTGTTCTTTACCGTTATGCTGATATTTGCCGCCAGGATTTTCCAGCGCAt                                                    >  1:2593680/1‑71 (MQ=255)
                    aTTTTGATATTCGGATACGTTTGTTCTTTACCGTTATGCTGATATTTGCCGCCAGGATTTTCCAGCGCATc                                                   <  1:1507028/71‑1 (MQ=255)
                    aTTTTGATATTCGGATACGTTTGTTCTTTACCGTTATGCTGATATTTGCCGCCAGGATTTTCCAGCGCATc                                                   <  1:5390111/71‑1 (MQ=255)
                     ttttGATATTCGGATACGTTTGTTCTTTACCGTTATGCTGATATTTGCCGCCAGGATTTTc                                                            >  1:5292358/1‑61 (MQ=255)
                     ttttGATATTCGGATACGTTTGTTCTTTACCGTTATGCTGATATTTGCCGCCAGGATTTTCCAGCGCATc                                                   >  1:3886777/1‑70 (MQ=255)
                      tttGATATTCGGATACGTTTGTTCTTTACCGTTATGCTGATATTTGCCGCCAGGATTTTCCAGCGCATcga                                                 >  1:3963912/1‑71 (MQ=255)
                       ttGATATTCGGATACGTTTGTTCTTTACCGTTATGCTGATATTTGCCGCCAGGATTTTCCAGCGCATcgac                                                >  1:1933358/1‑71 (MQ=255)
                       ttGATATTCGGATACGTTTGTTCTTTACCGTTATGCTGATATTTGCCGCCAGGATTTTCCAGCGCATcgac                                                >  1:2727893/1‑71 (MQ=255)
                        tGATATTCGGATACGTTTGTTCTTTACCGTTATGCTGATATTTGCCGCCAGGATTTTCCAGCGCATcgac                                                <  1:4129135/70‑1 (MQ=255)
                          atatTCGGATACGTTTGTTCTTTACCGTTATGCTGATATTTGCCGCCAGGATTTTCCAGCGCATcgacga                                              <  1:3541201/70‑1 (MQ=255)
                             ttCGGATACGTTTGTTCTTTACCGTTATGCTGATATTTGCCGCCAGGATTTTCCAGCGCATCGACGATAc                                           <  1:3431187/70‑1 (MQ=255)
                              tCGGATACGTTTGTTCTTTACCGTTATGCTGATATTTGCCGCCAGGATTTTCCAGCGCATcga                                                 <  1:1494828/63‑1 (MQ=255)
                               cGGATACGTTTGTTCTTTACCGTTATGCTGATATTTGCCGCCAg                                                                   >  1:4405822/1‑44 (MQ=255)
                                ggATACGTTTGTTCTTTACCGTTATGCTGATATTTGCCGCCAGGATTTTCCAGCGCATCGACGATACGGGc                                       <  1:750469/71‑1 (MQ=255)
                                   tACGTTTGTTCTTTACCGTTATGCTGATATTTGCCGCCAGGATTTTcc                                                           <  1:3220034/48‑1 (MQ=255)
                                       tttGTTCTTTACCGTTATGCTGATATTTGCCGCCAGGATTTTCCAGCGCATCGACGATa                                            <  1:1846398/59‑1 (MQ=255)
                                       tttGTTCTTTACCGTTATGCTGATATTTGCCGCCAGGATTTTCCAGCGCATCGACGATa                                            <  1:1396129/59‑1 (MQ=255)
                                        ttGTTCTTTACCGTTATGCTGATATTTGCCGCCAGGATTTTCCAGCGCATcgacga                                              >  1:4249166/1‑56 (MQ=255)
                                        ttGTTCTTTACCGTTATGCTGATATTTGCCGCCAGGATTTTCCAGCGCATcgacga                                              >  1:5074405/1‑56 (MQ=255)
                                        ttGTTCTTTACCGTTATGCTGATATTTGCCGCCAGGATTTTCCAGCGCATcgacga                                              >  1:3731956/1‑56 (MQ=255)
                                        ttGTTCTTTACCGTTATGCTGATATTTGCCGCCAGGATTTTCCAGCGCATcgacga                                              >  1:3849606/1‑56 (MQ=255)
                                        ttGTTCTTTACCGTTATGCTGATATTTGCCGCCAGGATTTTCCAGCGCATcgacga                                              >  1:3067841/1‑56 (MQ=255)
                                        ttGTTCTTTACCGTTATGCTGATATTTGCCGCCAGGATTTTCCAGCGCATCGACGATACGGGCAACAGGaa                               >  1:2301660/1‑71 (MQ=255)
                                         tGTTCTTTACCGTTATGCTGATATTTGCCGCCAGGATTTTCCAGCGCATcga                                                 >  1:1171895/1‑52 (MQ=255)
                                         tGTTCTTTACCGTTATGCTGATATTTGCCGCCAGGATTTTCCAGCGCATCGACGATACGGGCaa                                     >  1:4416331/1‑64 (MQ=255)
                                         tGTTCTTTACCGTTATGCTGATATTTGCCGCCAGGATTTTCCAGCGCATCGACGATACGGGCAACAGGAAt                              >  1:850539/1‑71 (MQ=255)
                                         tGTTCTTTACCGTTATGCTGATATTTGCCGCCAGGATTTTCCAGCGCATCGACGATACGGGCAACAGGAAt                              >  1:567388/1‑71 (MQ=255)
                                         tGTTCTTTACCGTTATGCTGATATTTGCCGCCAGGATTTTCCAGCGCATCGACGATACGGGCAACAGGAAt                              >  1:5262237/1‑71 (MQ=255)
                                          gTTCTTTACCGTTATGCTGATATTTGCCGCCAGGAt                                                                <  1:4406132/36‑1 (MQ=255)
                                          gTTCTTTACCGTTATGCTGATATTTGCCGCCAGGATTTTCCAGCGCATCGACGATACGGGCATCAGGAATc                             >  1:168942/1‑71 (MQ=255)
                                          gTTCTTTACCGTTATGCTGATATTTGCCGCCAGGATTTTCCAGCGCATCGACGATACGGGCAACAGGAAt                              <  1:4522051/70‑1 (MQ=255)
                                          gTTCTTTACCGTTATGCTGATATTTGCCGCCAGGATTTTCCAGCGCATCGACGATACGGGCAACAGGAATc                             >  1:3700420/1‑71 (MQ=255)
                                           ttCTTTACCGTTATGCTGATATTTGCCGCCAGGATTTTCCAGCGCATcgacga                                              >  1:2298950/1‑53 (MQ=255)
                                           ttCTTTACCGTTATGCTGATATTTGCCGCCAGGATTTTCCAGCGCATc                                                   <  1:4758252/48‑1 (MQ=255)
                                           ttCTTTACCGTTATGCTGATATTTGCCGCCAGGATTTTCCAGCGCATc                                                   <  1:1679952/48‑1 (MQ=255)
                                           ttCTTTACCGTTATGCTGATATTTGCCGCCAGGATTTTCCAGCGCATCGACGATACGGGCAACAGGa                                <  1:5199668/67‑1 (MQ=255)
                                             cTTTACCGTTATGCTGATATTTGCCGCCAGGATTTTCCAGCGCATcga                                                 <  1:3090719/48‑1 (MQ=255)
                                             cTTTACCGTTATGCTGATATTTGCCGCCAGGATTTTCCAGCGCATcga                                                 >  1:2381079/1‑48 (MQ=255)
                                             cTTTACCGTTATGCTGATATTTGCCGCCAGGATTTTCCAGCGCATCGACGATACGGGCAACAGGAAt                              >  1:5448335/1‑67 (MQ=255)
                                             cTTTACCGTTATGCTGATATTTGCCGCCAGGATTTTCCAGCGCATCGACGATACGGGCAACAGGAAt                              >  1:429945/1‑67 (MQ=255)
                                             cTTTACCGTTATGCTGATATTTGCCGCCAGGATTTTCCAGCGCATCGACGATACGGGCAACAGGAATCGcc                          >  1:2998410/1‑71 (MQ=255)
                                              tttACCGTTATGCTGATATTTGCCGCCAGGATTTTCCAGCGCATc                                                   <  1:3247156/45‑1 (MQ=255)
                                               ttACCGTTATGCTGATATTTGCCGCCAGGATTTTCCAGCGCATCGACGATACGGGCAACAGGAATCGCCGt                        >  1:128213/1‑71 (MQ=255)
                                                tACCGTTATGCTGATATTTGCCGCCAGGATTTTCCAGCGCATCGACGATACGGGCAACAGGAATCGCCGt                        <  1:340776/70‑1 (MQ=255)
                                                tACCGTTATGCTGATATTTGCCGCCAGGATTTTCCAGCGCATCGACGATACGGGCAACAGGAATCGCCGt                        >  1:1514098/1‑70 (MQ=255)
                                                tACCGTTATGCTGATATTTGCCGCCAGGATTTTCCAGCGCATCGACGATACGGGCAACAGGAATCGCCGt                        <  1:4964833/70‑1 (MQ=255)
                                                tACCGTTATGCTGATATTTGCCGCCAGGATTTTCCAGCGCATCGACGATACGGGCAACAGGAATCGCCGt                        <  1:4261229/70‑1 (MQ=255)
                                                  ccGTTGTGTTGATATGCGCCACCAGGGTTTTCCAGTGCTTCAACAATGCGGGCAACAggg                                >  1:3145899/1‑59 (MQ=255)
                                                  ccGTTGTGTTGATATGCGCCACCAGGGTTTTCCAGTGCTTCAACAATGCGGGCAACAggg                                >  1:5087717/1‑59 (MQ=255)
                                                  ccGTTGTGTTGATATGCGCCACCAGGGTTTTCCAGTGCTTCAACAATGCGGGCAACAGGGAttt                            >  1:1695375/1‑62 (MQ=255)
                                                       aTGCTGATATTTGCCGCCAGGATTTTCCAGCGCATCGACGATACGGGCAACAGGAATCGCCGTCATTCCg                 <  1:5195118/70‑1 (MQ=255)
                                                        tGCTGATATTTGCCGCCAGGATTTTCCAGCGCATCGACGATACGGGCAACAGGAATCGCCGTCATTCCGc                >  1:1056701/1‑70 (MQ=255)
                                                        tGCTGATATTTGCCGCCAGGATTTTCCAGCGCATCGACGATACGGGCAACAGGAATCGCCGTCATTCCGc                >  1:4152849/1‑70 (MQ=255)
                                                         gCTGATATTTGCCGCCAGGATTTTCCAGCGCATcgac                                                >  1:4065946/1‑37 (MQ=255)
                                                         gCTGATATTTGCCGCCAGGATTTTCCAGCGCATcgac                                                >  1:5494210/1‑37 (MQ=255)
                                                         gCTGATATTTGCCGCCAGGATTTTCCAGCGCATcgac                                                >  1:2703606/1‑37 (MQ=255)
                                                          ttGATATGCGCCACCAGGGTTTTCCAGTGCTTCAACAATGCGGGCAACAGGGAt                              <  1:5642868/53‑1 (MQ=255)
                                                           tGATATTTGCCGCCAGGATTTTCCAGCGCATCGACGATa                                            <  1:3528062/39‑1 (MQ=255)
                                                           tGATATGCGCCACCAGGGTTTTCCAGTGCTTCAACAATGCGGGCAACAGGGAt                              <  1:5274536/53‑1 (MQ=255)
                                                             atatTTGCCGCCAGGATTTTCCAGCGCATCGACGATACGGGCAACAGGAATCGCCGTCATTCCGCcatca           <  1:5301244/70‑1 (MQ=255)
                                                              tatTTGCCGCCAGGATTTTCCAGCGCATCGACGATACGGGCAACAGGAATCGCCGTCATTc                   >  1:4371361/1‑61 (MQ=255)
                                                              tatTTGCCGCCAGGATTTTCCAGCGCATCGACGATACGGGCAACAGGAATCGCCGTCATTc                   >  1:3113091/1‑61 (MQ=255)
                                                              tatTTGCCGCCAGGATTTTCCAGCGCATCGACGATACGGGCAACAGGAATCGCCGTCATTc                   >  1:4898602/1‑61 (MQ=255)
                                                              tatTTGCCGCCAGGATTTTCCAGCGCATCGACGATACGGGCAACAGGAATCGCCGTCATTCCGCcatcat          >  1:489677/1‑70 (MQ=255)
                                                                tttGCCGCCAGGATTTTCCAGCGCATCGACGATACGGGCAACAGGaa                               >  1:4554316/1‑47 (MQ=255)
                                                                tttGCCGCCAGGATTTTCCAGCGCATCGACGATACGGGCAACAGGa                                <  1:1910830/46‑1 (MQ=255)
                                                                  tGCCGCCAGGATTTTCCAGCGCATCGACGATACGGGCAACAGGAATCGCCGTCATTCCGCCATCATCAGCg     >  1:332939/1‑71 (MQ=255)
                                                                    ccgccAGGATTTTCCAGCGCATCGACGATACGGGCAACAGGAATCGCCGTCATTCCGc                <  1:476519/58‑1 (MQ=255)
                                                                     cgccAGGATTTTCCAGCGCATCGACGATACGGGCAACAGGAATCGCCGTCATTCCGCcatcat          >  1:1364440/1‑63 (MQ=255)
                                                                     cgccAGGATTTTCCAGCGCATCGACGATACGGGCAACAGGAATCGCCGTCATTCCGCCATCATCAGCGGcc  <  1:4989511/71‑1 (MQ=255)
                                                                     cgccAGGATTTTCCAGCGCATCGACGATACGGGCAACAGGAATCGCCGTCATTCCGCCATCATCAGCGGcc  <  1:1653200/71‑1 (MQ=255)
                                                                     cgccAGGATTTTCCAGCGCATCGACGATACGGGCAACAGGAATCGCCGTCATTCCGCCATCATCAGCGGcc  >  1:82557/1‑71 (MQ=255)
                                                                      |                                                                     
CGCCACCCGCCCACCAGATCATTTTGATATTCGGATACGTTTGTTCTTTACCGTTATGCTGATATTTGCCGCCAGGATTTTCCAGCGCATCGACGATACGGGCAACAGGAATCGCCGTCATTCCGCCATCATCAGCGGCC  >  Exported/1259534‑1259673

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: