Predicted mutation
evidence seq id position mutation freq annotation gene description
RA minE 1,259,610 A→G 5.2% N438N (AAT→AAC 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
*minE1,259,6100AG 5.2% 281.6 / 3.8 96N438N (AAT→AACtorZtrimethylamine N‑oxide reductase system III, catalytic subunit
Reads supporting (aligned to +/- strand):  ref base A (54/37);  new base G (3/2);  total (57/39)
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.45e-01

CGCCCACCAGATCATTTTGATATTCGGATACGTTTGTTCTTTACCGTTATGCTGATATTTGCCGCCAGGATTTTCCAGCGCATCGACGATACGGGCAACAGGAATCGCCGTCATTCCGCCATCATCAGCGGCCTCGCTGG  >  minE/1259541‑1259680
                                                                     |                                                                      
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:287859/1‑71 (MQ=255)
        aGATCATTTTGATATTCGGATACGTTTGTTCTTTACCGTTATGCTGATATTTGCCGCCAGGATTTTCCAgc                                                               >  1:908143/1‑71 (MQ=255)
        aGATCATTTTGATATTCGGATACGTTTGTTCTTTACCGTTATGCTGATATTTGCCGCCAGGATTTTCCAgc                                                               >  1:2913737/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:608410/1‑60 (MQ=255)
            cATTTTGATATTCGGATACGTTTGTTCTTTACCGTTATGCTGATATTTGCCGCCAGGAtt                                                                      >  1:874394/1‑60 (MQ=255)
            cATTTTGATATTCGGATACGTTTGTTCTTTACCGTTATGCTGATATTTGCCGCCAGGAtt                                                                      >  1:2568848/1‑60 (MQ=255)
            cATTTTGATATTCGGATACGTTTGTTCTTTACCGTTATGCTGATATTTGCCGCCAGGATTTTCCAGCGCAt                                                           >  1:3971528/1‑71 (MQ=255)
            cATTTTGATATTCGGATACGTTTGTTCTTTACCGTTATGCTGATATTTGCCGCCAGGATTTTCCAGCGCAt                                                           >  1:2593680/1‑71 (MQ=255)
            cATTTTGATATTCGGATACGTTTGTTCTTTACCGTTATGCTGATATTTGCCGCCAGGATTTTCCAGCGCAt                                                           >  1:2478823/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)
                         ggATACGTTTGTTCTTTACCGTTATGCTGATATTTGCCGCCAGGATTTTCCAGCGCATCGACGATACGGGc                                              <  1:750469/71‑1 (MQ=255)
                            tACGTTTGTTCTTTACCGTTATGCTGATATTTGCCGCCAGGATTTTcc                                                                  <  1:3220034/48‑1 (MQ=255)
                                tttGTTCTTTACCGTTATGCTGATATTTGCCGCCAGGATTTTCCAGCGCATCGACGATa                                                   <  1:1396129/59‑1 (MQ=255)
                                tttGTTCTTTACCGTTATGCTGATATTTGCCGCCAGGATTTTCCAGCGCATCGACGATa                                                   <  1:1846398/59‑1 (MQ=255)
                                 ttGTTCTTTACCGTTATGCTGATATTTGCCGCCAGGATTTTCCAGCGCATcgacga                                                     >  1:3731956/1‑56 (MQ=255)
                                 ttGTTCTTTACCGTTATGCTGATATTTGCCGCCAGGATTTTCCAGCGCATcgacga                                                     >  1:3067841/1‑56 (MQ=255)
                                 ttGTTCTTTACCGTTATGCTGATATTTGCCGCCAGGATTTTCCAGCGCATcgacga                                                     >  1:4249166/1‑56 (MQ=255)
                                 ttGTTCTTTACCGTTATGCTGATATTTGCCGCCAGGATTTTCCAGCGCATcgacga                                                     >  1:3849606/1‑56 (MQ=255)
                                 ttGTTCTTTACCGTTATGCTGATATTTGCCGCCAGGATTTTCCAGCGCATcgacga                                                     >  1:5074405/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:5262237/1‑71 (MQ=255)
                                  tGTTCTTTACCGTTATGCTGATATTTGCCGCCAGGATTTTCCAGCGCATCGACGATACGGGCAACAGGAAt                                     >  1:567388/1‑71 (MQ=255)
                                  tGTTCTTTACCGTTATGCTGATATTTGCCGCCAGGATTTTCCAGCGCATCGACGATACGGGCAACAGGAAt                                     >  1:850539/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:2381079/1‑48 (MQ=255)
                                      cTTTACCGTTATGCTGATATTTGCCGCCAGGATTTTCCAGCGCATcga                                                        <  1:3090719/48‑1 (MQ=255)
                                      cTTTACCGTTATGCTGATATTTGCCGCCAGGATTTTCCAGCGCATCGACGATACGGGCAACAGGAAt                                     >  1:429945/1‑67 (MQ=255)
                                      cTTTACCGTTATGCTGATATTTGCCGCCAGGATTTTCCAGCGCATCGACGATACGGGCAACAGGAAt                                     >  1:5448335/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:4964833/70‑1 (MQ=255)
                                         tACCGTTATGCTGATATTTGCCGCCAGGATTTTCCAGCGCATCGACGATACGGGCAACAGGAATCGCCGt                               >  1:1514098/1‑70 (MQ=255)
                                         tACCGTTATGCTGATATTTGCCGCCAGGATTTTCCAGCGCATCGACGATACGGGCAACAGGAATCGCCGt                               <  1:340776/70‑1 (MQ=255)
                                         tACCGTTATGCTGATATTTGCCGCCAGGATTTTCCAGCGCATCGACGATACGGGCAACAGGAATCGCCGt                               <  1:4261229/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)
                                                  gCTGATATTTGCCGCCAGGATTTTCCAGCGCATcgac                                                       >  1:5494210/1‑37 (MQ=255)
                                                  gCTGATATTTGCCGCCAGGATTTTCCAGCGCATcgac                                                       >  1:2703606/1‑37 (MQ=255)
                                                  gCTGATATTTGCCGCCAGGATTTTCCAGCGCATcgac                                                       >  1:4065946/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:4898602/1‑61 (MQ=255)
                                                       tatTTGCCGCCAGGATTTTCCAGCGCATCGACGATACGGGCAACAGGAATCGCCGTCATTc                          >  1:3113091/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:1653200/71‑1 (MQ=255)
                                                              cgccAGGATTTTCCAGCGCATCGACGATACGGGCAACAGGAATCGCCGTCATTCCGCCATCATCAGCGGcc         >  1:82557/1‑71 (MQ=255)
                                                              cgccAGGATTTTCCAGCGCATCGACGATACGGGCAACAGGAATCGCCGTCATTCCGCCATCATCAGCGGcc         <  1:4989511/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:4889252/43‑1 (MQ=255)
                                                                    gATTTTCCAGCGCATCGACGATACGGGCAACAGGAATCGCCGt                               <  1:5176533/43‑1 (MQ=255)
                                                                    gATTTTCCAGCGCATCGACGATACGGGCAACAGGAATCGCCGt                               <  1:1315152/43‑1 (MQ=255)
                                                                    gATTTTCCAGCGCATCGACGATACGGGCAACAGGAATCGCCGt                               <  1:5182770/43‑1 (MQ=255)
                                                                    gATTTTCCAGCGCATCGACGATACGGGCAACAGGAATCGCCGTCATTCCGCCatcatc                <  1:4849533/58‑1 (MQ=255)
                                                                     aTTTTCCAGCGCATCGACGATACGGGCAACAGGAATCGCCGTc                              <  1:2801507/43‑1 (MQ=255)
                                                                     aTTTTCCAGCGCATCGACGATACGGGCAACAGGAATCGCCGTc                              <  1:782481/43‑1 (MQ=255)
                                                                     aTTTTCCAGCGCATCGACGATACGGGCAACAGGAATCGCCGTCATTCCGCCATCATCAGCGGCCTCGCTgg  <  1:210047/71‑1 (MQ=255)
                                                                     |                                                                      
CGCCCACCAGATCATTTTGATATTCGGATACGTTTGTTCTTTACCGTTATGCTGATATTTGCCGCCAGGATTTTCCAGCGCATCGACGATACGGGCAACAGGAATCGCCGTCATTCCGCCATCATCAGCGGCCTCGCTGG  >  minE/1259541‑1259680

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: