Predicted mutation
evidence seq id position mutation freq annotation gene description
RA minE 1,259,599 T→G 5.7% K442T (AAA→ACA) ‡ 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,5990TG 5.7% 233.1 / 6.5 87K442T (AAA→ACA) ‡torZtrimethylamine N‑oxide reductase system III, catalytic subunit
Reads supporting (aligned to +/- strand):  ref base T (54/28);  new base G (3/2);  total (57/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 = 9.66e-01

TTACCGCCACCCGCCCACCAGATCATTTTGATATTCGGATACGTTTGTTCTTTACCGTTATGCTGATATTTGCCGCCAGGATTTTCCAGCGCATCGACGATACGGGCAACAGGAATCGCCGTCATTCCGCCATCAT  >  minE/1259530‑1259665
                                                                     |                                                                  
ttACCGCCACCCGCCCACCAGATCATTTTGATATTCGGATACGTTTGTTCTTTACCGTTATGCTGATAttt                                                                   <  1:2548149/71‑1 (MQ=255)
 tACCGCCACCCGCCCACCAGATCATTTTGATATTCGGATACGTTTGTTCTTTACCGTTATGCTGATATTTg                                                                  >  1:4060732/1‑71 (MQ=255)
  aCCGCCACCCGCCCACCAGATCATTTTGATATTCGGATACGTTTGTTCTTTACCGTTATGCTGATATTTgc                                                                 >  1:4235548/1‑71 (MQ=255)
    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:287859/1‑71 (MQ=255)
                   aGATCATTTTGATATTCGGATACGTTTGTTCTTTACCGTTATGCTGATATTTGCCGCCAGGATTTTCCAgc                                                >  1:2913737/1‑71 (MQ=255)
                   aGATCATTTTGATATTCGGATACGTTTGTTCTTTACCGTTATGCTGATATTTGCCGCCAGGATTTTCCAgc                                                >  1:908143/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)
                     aTCATTTTGATATTCGGATACGTTTGTTCTTTACCGTTATGCTGATATTTgcc                                                                <  1:5461144/53‑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:2478823/1‑71 (MQ=255)
                       cATTTTGATATTCGGATACGTTTGTTCTTTACCGTTATGCTGATATTTGCCGCCAGGATTTTCCAGCGCAt                                            >  1:2593680/1‑71 (MQ=255)
                        aTTTTGATATTCGGATACGTTTGTTCTTTACCGTTATGCTGATATTTGCCGCCAGGATTTTCCAGCGCATc                                           <  1:5390111/71‑1 (MQ=255)
                        aTTTTGATATTCGGATACGTTTGTTCTTTACCGTTATGCTGATATTTGCCGCCAGGATTTTCCAGCGCATc                                           <  1:1507028/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:2727893/1‑71 (MQ=255)
                           ttGATATTCGGATACGTTTGTTCTTTACCGTTATGCTGATATTTGCCGCCAGGATTTTCCAGCGCATcgac                                        >  1:1933358/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)
                                  tCGGATACGTTTGTTCTTTACCGTTATGCTGATATTTg                                                                  >  1:1226866/1‑38 (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:5074405/1‑56 (MQ=255)
                                            ttGTTCTTTACCGTTATGCTGATATTTGCCGCCAGGATTTTCCAGCGCATcgacga                                      >  1:4249166/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: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:1679952/48‑1 (MQ=255)
                                               ttCTTTACCGTTATGCTGATATTTGCCGCCAGGATTTTCCAGCGCATc                                           <  1:4758252/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:1514098/1‑70 (MQ=255)
                                                    tACCGTTATGCTGATATTTGCCGCCAGGATTTTCCAGCGCATCGACGATACGGGCAACAGGAATCGCCGt                <  1:4964833/70‑1 (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:1056701/1‑70 (MQ=255)
                                                            tGCTGATATTTGCCGCCAGGATTTTCCAGCGCATCGACGATACGGGCAACAGGAATCGCCGTCATTCCGc        >  1:4152849/1‑70 (MQ=255)
                                                             gCTGATATTTGCCGCCAGGATTTTCCAGCGCATcgac                                        >  1:4065946/1‑37 (MQ=255)
                                                             gCTGATATTTGCCGCCAGGATTTTCCAGCGCATcgac                                        >  1:2703606/1‑37 (MQ=255)
                                                             gCTGATATTTGCCGCCAGGATTTTCCAGCGCATcgac                                        >  1:5494210/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)
                                                                     |                                                                  
TTACCGCCACCCGCCCACCAGATCATTTTGATATTCGGATACGTTTGTTCTTTACCGTTATGCTGATATTTGCCGCCAGGATTTTCCAGCGCATCGACGATACGGGCAACAGGAATCGCCGTCATTCCGCCATCAT  >  minE/1259530‑1259665

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: