Predicted mutation
evidence seq id position mutation freq annotation gene description
RA minE 1,260,748 G→A 9.7% S59L (TCA→TTA)  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,260,7480GA 9.7% 115.7 / 4.1 41S59L (TCA→TTA) torZtrimethylamine N‑oxide reductase system III, catalytic subunit
Reads supporting (aligned to +/- strand):  ref base G (15/22);  new base A (2/2);  total (17/24)
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.95e-01

TGTGTACCTGATCCGCCGCCGTAGACTGTAAGGAATTCGGTATGGTTTTCGCCAGCGCGCCTGTTGAAGAAACAATCTTGCCGTCCTTCACTTCTACGTTCATCGCTCCCCAACGTCCCGCGGTGAGGATTTTACC  >  minE/1260683‑1260818
                                                                 |                                                                      
tgtgtACCTGATCCGCCGCCGTAGACTGTAAGGAATTCGGTATTGTTTTCGCCAGCGCGCCTGTTgaagaa                                                                   <  1:1659084/71‑1 (MQ=255)
tgtgtACCTGATCCGCCGCCGTAGACTGTAAGGAATTCGGTATGGTTTTCGCCAGCGCGCCTGTTgaagaa                                                                   <  1:2159942/71‑1 (MQ=255)
    tACCTGATCCGCCGCCGTAGACTGTAAGGAATTCGGTATGGTTTTCGCCAGCGCGCCTGTTGAAGAAACa                                                                <  1:1653838/70‑1 (MQ=255)
      ccTGATCCGCCGCCGTAGACTGTAAGGAATTCGGTATGGTTTTCGCCAGCGCGCCTGTTGAAGAAACAAt                                                              <  1:803287/70‑1 (MQ=255)
       cTGATCCGCCGCCGTAGACTGTAAGGAATTCGGTATGGTTTTCGCCAGCGCGCCTGTTGAAGAAACAATCt                                                            <  1:1296703/71‑1 (MQ=255)
         gATCCGCCGCCGTAGACTGTAAGGAATTCGGTATGGTTTTCGCCAGCGCGCCTGTTGAAGAAACAATCTTg                                                          >  1:1343445/1‑71 (MQ=255)
         gATCCGCCGCCGTAGACTGTAAGGAATTCGGTATGGTTTTCGCCAGCGCGCCTGTTGAAGAAACAATCTTg                                                          >  1:673566/1‑71 (MQ=255)
          aTCCGCCGCCGTAGACTGTAAGGAATTCGGTATGGTTTTCGCCAGCGCGCCTGTTGAAGAAACAATCtt                                                           <  1:1080835/69‑1 (MQ=255)
          aTCCGCCGCCGTAGACTGTAAGGAATTCGGTATGGTTTTCGCCAGCGCGCCTGTTGAAGAAACAATCtt                                                           <  1:1437612/69‑1 (MQ=255)
          aTCCGCCGCCGTAGACTGTAAGGAATTCGGTATGGTTTTCGCCAGCGCGCCTGTTGAAGAAACAATCTTGc                                                         <  1:863163/71‑1 (MQ=255)
              gccgccgTAGACTGTAAGGAATTCGGTATGGTTTTCGCCAGCGCGCCTGTTGAAGAAACAATCTTGCCGTc                                                     <  1:2037837/71‑1 (MQ=255)
                cgccgTAGACTGTAAGGAATTCGGTATGGTTTTCGCCAGCGCGCCTGTTGAAGAAACAATCTTGCCGTCCt                                                   <  1:1329605/71‑1 (MQ=255)
                    gTAGACTGTAAGGAATTCGGTATGGTTTTCGCCAGCGCGCCTGTTGAAGAAAGAATCTTGCCGTCCtt                                                  <  1:1288142/68‑1 (MQ=255)
                     tAGACTGTAAGGAATTCGGTATGGTTTTCGCCAGCGCGCCTGTTGAAGAAACAATCTTGCCGTCCTTCACt                                              <  1:129701/71‑1 (MQ=255)
                          tGTAAGGAATTCGGTATGGTTTTCGCCAGCGCGCCTGTTGAAGAAACAATCTTGCCGTCCTTCACTTCTa                                          >  1:1426765/1‑70 (MQ=255)
                          tGTAAGGAATTCGGTATGGTTTTCGCCAGCGCGCCTGTTGAAGAAACAATCTTGCCGTCCTTCACTTCTAc                                         >  1:2013647/1‑71 (MQ=255)
                            tAAGGAATTCGGTATGGTTTTCGCCAGCGCGCCTGTTGAAGAAACaa                                                               <  1:1590217/47‑1 (MQ=255)
                            tAAGGAATTCGGTATGGTTTTCGCCAGCGCGCCTGTTGAAGAAACaa                                                               <  1:1348259/47‑1 (MQ=255)
                            tAAGGAATTCGGTATGGTTTTCGCCAGCGCGCCTGTTGAAGAAACaa                                                               <  1:1128939/47‑1 (MQ=255)
                            tAAGGAATTCGGTATGGTTTTCGCCAGCGCGCCTGTTGAAGAAACAATCTTGCCGTCCTTCACTTCTACGt                                       >  1:172541/1‑71 (MQ=255)
                              aGGAATTCGGTATGGTTTTCGCCAGCGCGCCTGTTGAAGAAACAATCTTGCCGTCCTTCACTTCTACGtt                                      >  1:549128/1‑70 (MQ=255)
                               ggAATTCGGTATGGTTTTCGCCAGCGCGCCTGTTGAAGAAACAATCTTGCCGTCCTTCACTTCTACGTTCa                                    <  1:1547844/71‑1 (MQ=255)
                                  aTTCGGTATGGTTTTCGCCAGCGCGCCTGTTGAAGAAACAATCTTGCCGTc                                                     <  1:26751/51‑1 (MQ=255)
                                   ttCGGTATGGTTTTCGCCAGCGCGCCTGTTgaaga                                                                    <  1:901420/35‑1 (MQ=255)
                                   ttCGGTATGGTTTTCGCCAGCGCGCCTGTTGAAGAAACa                                                                <  1:773835/39‑1 (MQ=255)
                                   ttCGGTATGGTTTTCGCCAGCGCGCCTGTTGAAGAAACAATCTTGCCg                                                       <  1:1315889/48‑1 (MQ=255)
                                     cGGTATGGTTTTCGCCAGCGCGCCTGTTGAAGaaa                                                                  >  1:1794627/1‑35 (MQ=255)
                                         aTGGTTTTCGCCAGCGCGCCTGTTGAAGAAACAATCTTGCCGTCCTTCACTTCTACGTTCa                                    >  1:2193782/1‑61 (MQ=255)
                                         aTGGTTTTCGCCAGCGCGCCTGTTGAAGAAACAATCTTGCCGTCCTTCACTTCTACGTTCa                                    >  1:994084/1‑61 (MQ=255)
                                         aTGGTTTTCGCCAGCGCGCCTGTTGAAGAAACAATCTTGCCGTCCTTCACTTCTACGTTCATCGCTcccc                           >  1:1990284/1‑70 (MQ=255)
                                         aTGGTTTTCGCCAGCGCGCCTGTTGAAGAAACAATCTTGCCGTCCTTCACTTCTACGTTCATCGCTcccc                           >  1:1901704/1‑70 (MQ=255)
                                              tttCGCCAGCGCGCCTGTTGAAGAAACAATCTTGc                                                         <  1:1645690/35‑1 (MQ=255)
                                               ttCGCCAGCGCGCCTGTTGAAGAAACAATCTTGCCGTCCTTCACTTCTACGTTCATc                                  >  1:1829246/1‑57 (MQ=255)
                                                  gCCAGCGCGCCTGTTGAAGAAACAATCTTGCCGTCCTTCACTTCTACGTTCATCGCTCCCCAACGTCCcgc                 <  1:1961687/71‑1 (MQ=255)
                                                  gCCAGCGCGCCTGTTGAAGAAACAATCTTGCCGTCCTTCACTTCTACGTTCATCGCTCCCCAACGTCCcgc                 <  1:289399/71‑1 (MQ=255)
                                                  gCCAGCGCGCCTGTTGAAGAAACAATCTTGCCGTCCTTCACTTCTACGTTCATCGCTCCCCAACGTCCcg                  <  1:943327/70‑1 (MQ=255)
                                                   ccAGCGCGCCTGTTGAAGAAACAATCTTGCCGTCCTTCACTTCTACGTTCATCGCTCCCCAACGTCCcgcg                >  1:537066/1‑71 (MQ=255)
                                                             tGTTGAAGAAACAATCTTGCCGTCCTTCACTTCTACGTTCATCGCTCCCCAACGTCCCGCGGTGAGGAttt      >  1:285055/1‑71 (MQ=255)
                                                             tGTTAAAGAAACAATCTTGCCGTCCTTCACTTCTa                                          >  1:1604748/1‑35 (MQ=255)
                                                             tGTTAAAGAAACAATCTTGCCGTCCTTCACTTCTAc                                         >  1:2092778/1‑36 (MQ=255)
                                                             tGTTAAAGAAACAATCTTGCCGTCCTTCACTTCTACGTTCATCGCt                               <  1:1233119/46‑1 (MQ=255)
                                                             tGTTAAAGAAACAATCTTGCCGTCCTTCACTTCTACGTTCATCGCt                               <  1:578193/46‑1 (MQ=255)
                                                               ttGAAGAAACAATCTTGCCGTCCTTCACTTCTACGTTCATCGCTCCCCAACGTCCCGCGGTGAGGAtttt     <  1:131927/70‑1 (MQ=255)
                                                               ttGAAGAAACAATCTTGCCGTCCTTCACTTCTACGTTCATCGCTCCCCAACGTCCCGCGGTGAGGATTTTa    >  1:984904/1‑71 (MQ=255)
                                                                 gaagaaACAATCTTGCCGTCCTTCACTTCTACGTTCATCGCTCCCCAACGTCCCGCGGTGAGGAt        <  1:1393888/65‑1 (MQ=255)
                                                                 gaagaaACAATCTTGCCGTCCTTCACTTCTACGTTCATCGCTCCCCAACGTCCCGCGGTGAGGATTTTa    >  1:90284/1‑69 (MQ=255)
                                                                 gaagaaACAATCTTGCCGTCCTTCACTTCTACGTTCATCGCTCCCCAACGTCCCGCGGTGAGGATTTTAcc  >  1:718415/1‑71 (MQ=255)
                                                                 gaagaaACAATCTTGCCGTCCTTCACTTCTACGTTCATCGCTCCCCAACGTCCCGCGGTGAGGATTTTAcc  >  1:1860410/1‑71 (MQ=255)
                                                                 gaagaaACAATCTTGCCGTCCTTCACTTCTACGTTCATCGCTCCCCAACGTCCCGCGGTGAGGATTTTAc   <  1:1271956/70‑1 (MQ=255)
                                                                 gaagaaACAATCTTGCCGTCCTTCACTTCTACGTTCATCGCTCCCCAACGTCCCGCGGTGAGGATTTTAc   <  1:210535/70‑1 (MQ=255)
                                                                 |                                                                      
TGTGTACCTGATCCGCCGCCGTAGACTGTAAGGAATTCGGTATGGTTTTCGCCAGCGCGCCTGTTGAAGAAACAATCTTGCCGTCCTTCACTTCTACGTTCATCGCTCCCCAACGTCCCGCGGTGAGGATTTTACC  >  minE/1260683‑1260818

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: