Predicted mutation
evidence seq id position mutation freq annotation gene description
RA Exported 1,259,628 G→A 7.2% I432I (ATC→ATT 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,6280GA 7.2% 200.4 / 4.1 69I432I (ATC→ATTtorZtrimethylamine N‑oxide reductase system III, catalytic subunit
Reads supporting (aligned to +/- strand):  ref base G (33/31);  new base A (3/2);  total (36/33)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.80e-01

ATATTCGGATACGTTTGTTCTTTACCGTTATGCTGATATTTGCCGCCAGGATTTTCCAGCGCATCGACGATACGGGCAACAGGAATCGCCGTCATTCCGCCATCATCAGCGGCCTCGCTGGCATGTCCGGCAATCG  >  Exported/1259560‑1259695
                                                                    |                                                                   
atatTCGGATACGTTTGTTCTTTACCGTTATGCTGATATTTGCCGCCAGGATTTTCCAGCGCATcgacga                                                                    <  1:3541201/70‑1 (MQ=255)
   ttCGGATACGTTTGTTCTTTACCGTTATGCTGATATTTGCCGCCAGGATTTTCCAGCGCATCGACGATAc                                                                 <  1:3431187/70‑1 (MQ=255)
      ggATACGTTTGTTCTTTACCGTTATGCTGATATTTGCCGCCAGGATTTTCCAGCGCATCGACGATACGGGc                                                             <  1:750469/71‑1 (MQ=255)
             tttGTTCTTTACCGTTATGCTGATATTTGCCGCCAGGATTTTCCAGCGCATCGACGATa                                                                  <  1:1396129/59‑1 (MQ=255)
             tttGTTCTTTACCGTTATGCTGATATTTGCCGCCAGGATTTTCCAGCGCATCGACGATa                                                                  <  1:1846398/59‑1 (MQ=255)
              ttGTTCTTTACCGTTATGCTGATATTTGCCGCCAGGATTTTCCAGCGCATcgacga                                                                    >  1:3849606/1‑56 (MQ=255)
              ttGTTCTTTACCGTTATGCTGATATTTGCCGCCAGGATTTTCCAGCGCATcgacga                                                                    >  1:3731956/1‑56 (MQ=255)
              ttGTTCTTTACCGTTATGCTGATATTTGCCGCCAGGATTTTCCAGCGCATcgacga                                                                    >  1:3067841/1‑56 (MQ=255)
              ttGTTCTTTACCGTTATGCTGATATTTGCCGCCAGGATTTTCCAGCGCATcgacga                                                                    >  1:5074405/1‑56 (MQ=255)
              ttGTTCTTTACCGTTATGCTGATATTTGCCGCCAGGATTTTCCAGCGCATcgacga                                                                    >  1:4249166/1‑56 (MQ=255)
              ttGTTCTTTACCGTTATGCTGATATTTGCCGCCAGGATTTTCCAGCGCATCGACGATACGGGCAACAGGaa                                                     >  1:2301660/1‑71 (MQ=255)
               tGTTCTTTACCGTTATGCTGATATTTGCCGCCAGGATTTTCCAGCGCATCGACGATACGGGCaa                                                           >  1:4416331/1‑64 (MQ=255)
               tGTTCTTTACCGTTATGCTGATATTTGCCGCCAGGATTTTCCAGCGCATCGACGATACGGGCAACAGGAAt                                                    >  1:5262237/1‑71 (MQ=255)
               tGTTCTTTACCGTTATGCTGATATTTGCCGCCAGGATTTTCCAGCGCATCGACGATACGGGCAACAGGAAt                                                    >  1:567388/1‑71 (MQ=255)
               tGTTCTTTACCGTTATGCTGATATTTGCCGCCAGGATTTTCCAGCGCATCGACGATACGGGCAACAGGAAt                                                    >  1:850539/1‑71 (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)
                 ttCTTTACCGTTATGCTGATATTTGCCGCCAGGATTTTCCAGCGCATCGACGATACGGGCAACAGGa                                                      <  1:5199668/67‑1 (MQ=255)
                   cTTTACCGTTATGCTGATATTTGCCGCCAGGATTTTCCAGCGCATCGACGATACGGGCAACAGGAAt                                                    >  1:5448335/1‑67 (MQ=255)
                   cTTTACCGTTATGCTGATATTTGCCGCCAGGATTTTCCAGCGCATCGACGATACGGGCAACAGGAAt                                                    >  1:429945/1‑67 (MQ=255)
                   cTTTACCGTTATGCTGATATTTGCCGCCAGGATTTTCCAGCGCATCGACGATACGGGCAACAGGAATCGcc                                                >  1:2998410/1‑71 (MQ=255)
                     ttACCGTTATGCTGATATTTGCCGCCAGGATTTTCCAGCGCATCGACGATACGGGCAACAGGAATCGCCGt                                              >  1:128213/1‑71 (MQ=255)
                      tACCGTTATGCTGATATTTGCCGCCAGGATTTTCCAGCGCATCGACGATACGGGCAACAGGAATCGCCGt                                              >  1:1514098/1‑70 (MQ=255)
                      tACCGTTATGCTGATATTTGCCGCCAGGATTTTCCAGCGCATCGACGATACGGGCAACAGGAATCGCCGt                                              <  1:4261229/70‑1 (MQ=255)
                      tACCGTTATGCTGATATTTGCCGCCAGGATTTTCCAGCGCATCGACGATACGGGCAACAGGAATCGCCGt                                              <  1:4964833/70‑1 (MQ=255)
                      tACCGTTATGCTGATATTTGCCGCCAGGATTTTCCAGCGCATCGACGATACGGGCAACAGGAATCGCCGt                                              <  1:340776/70‑1 (MQ=255)
                        ccGTTGTGTTGATATGCGCCACCAGGGTTTTCCAGTGCTTCAACAATGCGGGCAACAggg                                                      >  1:5087717/1‑59 (MQ=255)
                        ccGTTGTGTTGATATGCGCCACCAGGGTTTTCCAGTGCTTCAACAATGCGGGCAACAggg                                                      >  1:3145899/1‑59 (MQ=255)
                        ccGTTGTGTTGATATGCGCCACCAGGGTTTTCCAGTGCTTCAACAATGCGGGCAACAGGGAttt                                                  >  1:1695375/1‑62 (MQ=255)
                             aTGCTGATATTTGCCGCCAGGATTTTCCAGCGCATCGACGATACGGGCAACAGGAATCGCCGTCATTCCg                                       <  1:5195118/70‑1 (MQ=255)
                              tGCTGATATTTGCCGCCAGGATTTTCCAGCGCATCGACGATACGGGCAACAGGAATCGCCGTCATTCCGc                                      >  1:4152849/1‑70 (MQ=255)
                              tGCTGATATTTGCCGCCAGGATTTTCCAGCGCATCGACGATACGGGCAACAGGAATCGCCGTCATTCCGc                                      >  1:1056701/1‑70 (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:3113091/1‑61 (MQ=255)
                                    tatTTGCCGCCAGGATTTTCCAGCGCATCGACGATACGGGCAACAGGAATCGCCGTCATTc                                         >  1:4898602/1‑61 (MQ=255)
                                    tatTTGCCGCCAGGATTTTCCAGCGCATCGACGATACGGGCAACAGGAATCGCCGTCATTc                                         >  1:4371361/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:82557/1‑71 (MQ=255)
                                           cgccAGGATTTTCCAGCGCATCGACGATACGGGCAACAGGAATCGCCGTCATTCCGCCATCATCAGCGGcc                        <  1:4989511/71‑1 (MQ=255)
                                           cgccAGGATTTTCCAGCGCATCGACGATACGGGCAACAGGAATCGCCGTCATTCCGCCATCATCAGCGGcc                        <  1:1653200/71‑1 (MQ=255)
                                               aGGATTTTCCAGCGCATCGACGATACGGGCAACAGGAATCGcc                                                <  1:4784928/43‑1 (MQ=255)
                                               aGGATTTTCCAGCGCATCGACGATACGGGCAACAGGAATCGCCGTCATTCCGCCATCATCAGCGGCCTCGc                    >  1:2855123/1‑71 (MQ=255)
                                                 gATTTTCCAGCGCATCGACGATACGGGCAACAGGa                                                      <  1:1339194/35‑1 (MQ=255)
                                                 gATTTTCCAGCGCATCGACGATACGGGCAACAGGa                                                      <  1:3863570/35‑1 (MQ=255)
                                                 gATTTTCCAGCGCATCGACGATACGGGCAACAGGAATCGCCGt                                              <  1:5176533/43‑1 (MQ=255)
                                                 gATTTTCCAGCGCATCGACGATACGGGCAACAGGAATCGCCGt                                              <  1:4889252/43‑1 (MQ=255)
                                                 gATTTTCCAGCGCATCGACGATACGGGCAACAGGAATCGCCGt                                              <  1:5182770/43‑1 (MQ=255)
                                                 gATTTTCCAGCGCATCGACGATACGGGCAACAGGAATCGCCGt                                              <  1:1315152/43‑1 (MQ=255)
                                                 gATTTTCCAGCGCATCGACGATACGGGCAACAGGAATCGCCGTCATTCCGCCatcatc                               <  1:4849533/58‑1 (MQ=255)
                                                  aTTTTCCAGCGCATCGACGATACGGGCAACAGGAATCGCCGTc                                             <  1:782481/43‑1 (MQ=255)
                                                  aTTTTCCAGCGCATCGACGATACGGGCAACAGGAATCGCCGTc                                             <  1:2801507/43‑1 (MQ=255)
                                                  aTTTTCCAGCGCATCGACGATACGGGCAACAGGAATCGCCGTCATTCCGCCATCATCAGCGGCCTCGCTgg                 <  1:210047/71‑1 (MQ=255)
                                                   ttttCCAGCGCATCGACGATACGGGCAACAGGAAt                                                    >  1:4680944/1‑35 (MQ=255)
                                                   ttttCCAGCGCATCGACGATACGGGCAACAGGAAt                                                    >  1:1220249/1‑35 (MQ=255)
                                                   ttttCCAGCGCATCGACGATACGGGCAACAGGAAt                                                    >  1:1115151/1‑35 (MQ=255)
                                                   ttttCCAGCGCATCGACGATACGGGCAACAGGAATCGCCGTCATTCCGCCATCATCAg                             >  1:1573795/1‑58 (MQ=255)
                                                   ttttCCAGCGCATCGACGATACGGGCAACAGGAATCGCCGTCATTCCGCCATCATCAg                             >  1:4939809/1‑58 (MQ=255)
                                                    tttCCAGCGCATCGACGATACGGGCAACAGGAATCGCCGTCATTCCGCcatcat                                <  1:895565/54‑1 (MQ=255)
                                                    tttCCAGCGCATCGACGATACGGGCAACAGGAATCGCCGTCATTCCGCcatcat                                <  1:634346/54‑1 (MQ=255)
                                                    tttCCAGCGCATCGACGATACGGGCAACAGGAATCGCCGTCATTCCGCcatcat                                <  1:3976185/54‑1 (MQ=255)
                                                    tttCCAGCGCATCGACGATACGGGCAACAGGAATCGCCGTCATTCCGCcatcat                                <  1:1938859/54‑1 (MQ=255)
                                                      tCCAGCGCATCGACGATACGGGCAACAGGAATCGCCGTCa                                            >  1:3915666/1‑40 (MQ=255)
                                                      tCCAGCGCATCGACGATACGGGCAACAGGAATCGCCGTCa                                            >  1:1238459/1‑40 (MQ=255)
                                                         aGCGCATCGACGATACGGGCAACAGGAATCGCCGTCATTCCGCCATCATCAGCGGCCTCGCTGGCATGTcc          >  1:1461784/1‑71 (MQ=255)
                                                         aGCGCATCGACGATACGGGCAACAGGAATCGCCGTCATTCCGCCATCATCAGCGGCCTCGCTGGCATGTcc          >  1:1176889/1‑71 (MQ=255)
                                                         aGCGCATCGACGATACGGGCAACAGGAATCGCCGTCATTCCGCCATCATCAGCGGCCTCGCTGGCATGTc           >  1:4970305/1‑70 (MQ=255)
                                                                 gacgaTACGGGCAACAGGAATCGCCGTCATTCCGCCATCATCAGCGGCCTCGCTGGCATGTCCGGCAATcg  >  1:4138658/1‑71 (MQ=255)
                                                                  acgaTACGGGCAACAGGAATCGCCGTCATTCCGCCAtcatca                              >  1:5117048/1‑42 (MQ=255)
                                                                  acgaTACGGGCAACAGGAATCGCCGTCATTCCGCCATCATCAGCGGCCTCGCTGGCATGTCCGGCAAt    >  1:4219019/1‑68 (MQ=255)
                                                                    |                                                                   
ATATTCGGATACGTTTGTTCTTTACCGTTATGCTGATATTTGCCGCCAGGATTTTCCAGCGCATCGACGATACGGGCAACAGGAATCGCCGTCATTCCGCCATCATCAGCGGCCTCGCTGGCATGTCCGGCAATCG  >  Exported/1259560‑1259695

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