Predicted mutation
evidence seq id position mutation freq annotation gene description
RA minE 1,259,781 C→G 28.0% P381P (CCG→CCC 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,7810CG28.0% 72.2 / 32.3 50P381P (CCG→CCCtorZtrimethylamine N‑oxide reductase system III, catalytic subunit
Reads supporting (aligned to +/- strand):  ref base C (15/21);  new base G (9/5);  total (24/26)
Fisher's exact test for biased strand distribution p-value = 2.11e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 6.78e-01

GCACGCCGCCAACACGCGTCGGGTTGCCACCGTTGGAGTAATGATAAGAGAAACCAAAACCGCCGCCCGGTGTACCAATCTGCCCTAACATCGCCGCCAGCGTCACCAGCATCCAGTGTTTCTGTTCACCATATTGCT  >  minE/1259714‑1259851
                                                                   |                                                                      
gcaCGCCGCCAACACGCGTCGGGTTGCCACCGTTGGAGTAATGATAAGAGAAACCAAAACCGCCGCCCGgt                                                                     >  1:3502688/1‑71 (MQ=255)
    gccgccAACACGCGTCGGGTTGCCACCGTTGGAGTAATGATAAGAGAAACCAAAACCGCCGCCCGGTGTAc                                                                 >  1:5665775/1‑71 (MQ=255)
     ccgccAACACGCGTCGGGTTGCCACCGTTGGAGTAATGATAAGAGAAACCAAAACCGCCGCCCGGTGTAcc                                                                <  1:1749692/71‑1 (MQ=255)
      cgccAACACGCGTCGGGTTGCCACCGTTGGAGTAATGATAAGAGAAACCAAAACCGCCGCCCGGTGTACCa                                                               <  1:1386230/71‑1 (MQ=255)
      cgccAACACGCGGCGGGTTGCCACCGTTGGAGTAATGATAAGAGAAACCAAAACCGCCGCCCGGTGTACCa                                                               <  1:5047619/71‑1 (MQ=255)
         cAACACGCGTCGGGTTGCCACCGTTGGAGTAATGATAAGAGAAACCAAAACCGCCGCCCGGTGTACCAATc                                                            >  1:2102408/1‑71 (MQ=255)
          aaCACGCGTCGGGTTGCCACCGTTGGAGTAATGATAAGAGAAACCAAAACCGCCGCCCGGTGTACCAATCt                                                           >  1:109426/1‑71 (MQ=255)
                cgTCGGGTTGCCACCGTTGGAGTAATGATAAGAGAAACCAAAACCGCCGCCCGGTGTACCAATCTGCCCt                                                      <  1:1088197/70‑1 (MQ=255)
                cgTCGGGTTGCCACCGTTGGAGTAATGATAAGAGAAACCAAAACCGCCGCCCGGTGTACCAATCTGCCCTa                                                     <  1:1946550/71‑1 (MQ=255)
                  tggggTTACCACCATTGGCAAAATGGTAAGAAAGACCAAAACCGCCGCCGGGTGTGCCGATTTGCCCCaa                                                    >  1:3184431/3‑70 (MQ=255)
                  tCGGGTTGCCACCGTTGGAGTAATGATAAGAGAAACCAAAACCGCCGCCCGGTGTACCAATCTGCCCTa                                                     <  1:742430/69‑1 (MQ=255)
                    gggTTGCCACCGTTGGAGTAATGATAAGAGAAACCAAAACCGCCGCCCGGTGTa                                                                  >  1:4083805/1‑54 (MQ=255)
                      gTTACCACCATTGGCAAAATGGTAAGAAAGACCAAAACCGCCGCCGGGTGTGCCGATTTGCCCCAACAttg                                               >  1:1004409/1‑69 (MQ=255)
                       ttaccaccaTTGGCAAAATGGTAAGAAAGACCAAAACCGCCGCCGGGTGTGCCGATTTGCCCCAACATtgc                                              >  1:2755281/4‑71 (MQ=255)
                       ttGCCACCGTTGGAGTAATGATAAGAGAAACCAAAACCGCCGCCCGGTGTACCAATCTGCCCTAACATcgc                                              >  1:4433378/1‑71 (MQ=255)
                        taccaccaTTGGCAAAATGGTAAGAAAGACCAAAACCGCCGCCGGGTGTGCCGATTTGCCCCAACATtgc                                              >  1:239651/3‑70 (MQ=255)
                          ccaccaTTGGCAAAATGGTAAGAAAGACCAAAACCGCCGCCGGGTGTGcc                                                                <  1:5319517/50‑1 (MQ=255)
                                  ggAGTAATTATAAGAGAAACCAAAACCGCCGCCCGGTGTACCa                                                               <  1:4634668/43‑1 (MQ=255)
                                   gAGTAATGATAAGAGAAACCAAAACCGCCGCCCGGTGTACCAATCTGCCCTAACATCGCCGCCAGCGTCAc                                  <  1:2844667/71‑1 (MQ=255)
                                   gAGTAATGATAAGAGAAACCAAAACCGCCGCCCGGTGTACCAATCTGCCCTAACATCGCCGCCAGCGTCAc                                  <  1:3032831/71‑1 (MQ=255)
                                    aGTAATGATAAGAGAAACCAAAACCGCCGCCCGGTGTACCAATCTGCCCTAACATCGCCGCCAGCGTCAc                                  >  1:4600250/1‑70 (MQ=255)
                                      aaaTGGTAAGAAAGACCAAAACCGCCGCCGGGTGTGCCGATTTGCCCCAACATTGCTGCCAGc                                       <  1:3618007/62‑1 (MQ=255)
                                      aaaTGGTAAGAAAGACCAAAACCGCCGCCGGGTGTGCCGATTTGCCCCAACATTGCTGCCAGCGTGAcga                                <  1:419678/69‑3 (MQ=255)
                                         tGATAAGAGAAACCAAAACCGCCGCCCGGTGTACCa                                                               <  1:1196808/36‑1 (MQ=255)
                                           gtAAGAAAGACCAAAACCGCCGCCGGGTGTGCCGATTTGCCCCAACATTGCTGCCAGCGTGACGATCATcc                          <  1:138006/70‑1 (MQ=255)
                                           gtAAGAAAGACCAAAACCGCCGCCGGGTGTGCCGATTTGCCCCAACATTGCTGCCAGCGTGACGATCATcc                          <  1:4525914/70‑1 (MQ=255)
                                            tAAGAGAAACCAAAACCGCCGCCCGGTGTACCAATCTGCCCTAACATCGCCGCCAGCGTCACCAGCATc                           >  1:3898621/1‑69 (MQ=255)
                                               gagaAACCAAAACCGCCGCCCGGTGTACCAATCTGCCCTAACATCGCCGCCAGCGTCACCAGc                              <  1:1149535/63‑1 (MQ=255)
                                               gagaAACCAAAACCGCCGCACGGTGTACCAATCTGCCCTAACATCGCCGCCAGCGTCACCAGc                              <  1:248870/63‑1 (MQ=255)
                                                 gaAACCAAAACCGCCGCCCGGTGTACCAATCTGCCCTAACAt                                                 >  1:3003638/1‑42 (MQ=255)
                                                 gaAACCAAAACCGCCGCCCGGTGTACCAATCTGCCCTAACATCGCCGCCAGCGTCACCAGCATCCAGTGtt                    >  1:2201941/1‑71 (MQ=255)
                                                  agaCCAAAACCGCCGCCGGGTGTGCCGATTTGCCCCAACATTGCTGCCAGCGTGACGATCATCCa                         >  1:402649/3‑65 (MQ=255)
                                                   gaCCAAAACCGCCGCCGGGTGTGCCGATTTGCCCCAACATTGCTGCCAGCGTGACGATCATCCAGTGtttt                  >  1:4128930/2‑70 (MQ=255)
                                                   aaCCAAAACCGCCGCCCGGTGTACCAATCTGCCCTAACATCGCCGCCAGc                                       <  1:3615561/50‑1 (MQ=255)
                                                   aaCCAAAACCGCCGCCCGGTGTACCAATCTGCCCTAACATCGCCGCCAGCGTCACCAGCa                             <  1:908069/60‑1 (MQ=255)
                                                   aaCCAAAACCGCCGCCCGGTGTACCAATCTGCCCTAACATCGCCGCCAGCGTCACCAGCATCCAGTGTTTc                  <  1:2122035/71‑1 (MQ=255)
                                                    aCCAAAACCGCCGCCCGGTGTACCAATCTGCCCTAACATCGCCGCCAGCGTCACCAGCATCCAGTGTTTc                  >  1:2733240/1‑70 (MQ=255)
                                                     ccAAAACCGCCGCCCGGTGTACCAATCTGCCCTAACATCGCCGCCAGCGTCACCa                                >  1:2902326/1‑55 (MQ=255)
                                                      cAAAACCGCCGCCGGGTGTGCCGATTTGCCCCAACATtgct                                             >  1:3988192/1‑40 (MQ=12)
                                                      cAAAACCGCCGCCGGGTGTGCCGATTTGCCCCAACATTGCTGCCAGCGTGACGATCATCCAGTGtttt                  >  1:3137196/1‑67 (MQ=255)
                                                      cAAAACCGCCGCCCGGTGTACCAATCTGCCCTAACATCGCCGCCAGCGTCACCAGCATcc                          <  1:1660741/60‑1 (MQ=255)
                                                      cAAAACCGCCGCCCGGTGTACCAATCTGCCCTAACATCGCCGCCAGCGTCACCAGCATCCAgt                       <  1:3178490/63‑1 (MQ=255)
                                                         aaCCGCCGCCCGGTGTACCAATCTGCCCTAACATCGCCGCCAGCGTCACCAGCATCCAGTGTTTCTGTTCa            >  1:5384322/1‑71 (MQ=255)
                                                         aaCCGCCGCCCGGTGTACCAATCTGCCCTAACATCGCCGCCAGCGTCACCAGCATCCAGTGTTTCTGTTCa            >  1:2015270/1‑71 (MQ=255)
                                                            cgccgcCCGGTGTACCAATCTGCCCTAACATCGCCGCCAGCGTCACCAGCATCCAGTGTTTCTGt               <  1:317868/65‑1 (MQ=255)
                                                             gccgccGGGTGTGCCGATTTGCCCCAACATTGCTGCCAGCGTGACGATCATCCAGTGTTTTTGCTCACCa         >  1:2928180/1‑70 (MQ=255)
                                                             gccgccCGGTGTACCAATCTGCCCTAACATCGCCGCCAGCGTCACCAGCATCCAGTGTTTCTGTTCACCa         >  1:3376499/1‑70 (MQ=255)
                                                               cgccCGGTGTACCAATCTGCCCTAACATCGCCGCCAGCGTCACCAGCATCCAGTGTTTCTg                <  1:5535675/61‑1 (MQ=255)
                                                                gccCGGTGTACCAATCTGCCCTAACATCGCCGCCAGCGTCACCAGCATCCAGTGttt                   <  1:2246028/57‑1 (MQ=255)
                                                                gccCGGTGTACCAATCTGCCCTAACATCGCCGCCAGCGTCACCAGCATCCAGTGTTTCTGTTCACCATAtt     <  1:4789950/71‑1 (MQ=255)
                                                                 cccGGTGTACCAATCTGCCCTAACATCGCCGCCAGCGTCACCAGCATCCAGTGTTTCTGTTCACCATATTg    >  1:2805669/1‑71 (MQ=255)
                                                                   cGGTGTACCAATCTGCCCTAACATCGCCGCCAGCGTCACCAGGATCCAGTGTTTCTGTTCACCATATTGCt  <  1:3908737/71‑1 (MQ=255)
                                                                   cGGTGTACCAATCTGCCCTAACATCGCCGCCAGCGTCACCAGCATCCAg                        >  1:1119392/1‑49 (MQ=255)
                                                                   |                                                                      
GCACGCCGCCAACACGCGTCGGGTTGCCACCGTTGGAGTAATGATAAGAGAAACCAAAACCGCCGCCCGGTGTACCAATCTGCCCTAACATCGCCGCCAGCGTCACCAGCATCCAGTGTTTCTGTTCACCATATTGCT  >  minE/1259714‑1259851

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