Predicted mutation
evidence seq id position mutation freq annotation gene description
RA minE 1,259,765 A→G 21.4% F387L (TTC→CTC)  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,7650AG21.4% 75.0 / 16.4 42F387L (TTC→CTC) torZtrimethylamine N‑oxide reductase system III, catalytic subunit
Reads supporting (aligned to +/- strand):  ref base A (11/22);  new base G (4/5);  total (15/27)
Fisher's exact test for biased strand distribution p-value = 6.98e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

CAGACATTTCCGGCAGCACGCCGCCAACACGCGTCGGGTTGCCACCGTTGGAGTAATGATAAGAGAAACCAAAACCGCCGCCCGGTGTACCAATCTGCCCTAACATCGCCGCCAGCGTCACCAGCATCCAGTGTTTC  >  minE/1259699‑1259835
                                                                  |                                                                      
cAGACATTTCCGGCAGCACGCCGCCAACACGCGTCGGGTTGCCACCGTTGGAGTAATGATAAGAGAAACCa                                                                    <  1:1028929/71‑1 (MQ=255)
 aGACATTTCCGGCAGCACGCCGCCAACACGCGTCGGGTTGCCACCGTTGGAGTAATGATAAGAGAAACCaa                                                                   <  1:4042393/71‑1 (MQ=255)
 aGACATTTCCGGCAGCACGCCGCCAACACGCGTCGGGTTGCCACCGTTGGAGTAATGATAAGAGAAACCaa                                                                   <  1:5161683/71‑1 (MQ=255)
 aGACATTTCCGGCAGCACGCCGCCAACACGCGTCGGGTTGCCACCGTTGGAGTAATGATAAGAGAAACCaa                                                                   <  1:2435816/71‑1 (MQ=255)
 aGACATTTCCGGCAGCACGCCGCCAACACGCGTCGGGTTGCCACCGTTGGAGTAATGATAAGAGAAACCaa                                                                   <  1:1930038/71‑1 (MQ=255)
 aGACATTTCCGGCAGCACGCCGCCAACACGCGTCGGGTTGCCACCGTTGGAGTAATGATAAGAGAAACCaa                                                                   <  1:1665503/71‑1 (MQ=255)
 aGACATTTCCGGCAGCACGCCGCCAACACGCGTCGGGTTGCCACCGTTGGAGTAATGATAAGAGAAACCaa                                                                   <  1:1502292/71‑1 (MQ=255)
   aCATTTCCGGCAGCACGCCGCCAACACGCGTCGGTTTGCCACCGTTGGAGTAATGATAAGAGAAACCaaaa                                                                 <  1:5297226/71‑1 (MQ=255)
   aCATTTCCGGCAGCACGCCGCCAACACGCGTCGGGTTGCCACCGTTGGAGTAATGATAAGAGAAAcc                                                                     <  1:2691006/67‑1 (MQ=255)
    cATTTCCGGCAGCACGCCGCCAACACGCGTCGGGTTGCCACCGTTGGAGTAATGATAAGAGAAACCAAAAc                                                                <  1:5537514/71‑1 (MQ=255)
               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)
                             cgcgTCGGGTTGCCACCGTTGGAGTAATGAAAAGAGAAACCAAAAccgccg                                                           >  1:2359074/1‑51 (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:4525914/70‑1 (MQ=255)
                                                          gtAAGAAAGACCAAAACCGCCGCCGGGTGTGCCGATTTGCCCCAACATTGCTGCCAGCGTGACGATCATcc          <  1:138006/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)
                                                                  aaCCAAAACCGCCGCCCGGTGTACCAATCTGCCCTAACATCGCCGCCAGc                       <  1:3615561/50‑1 (MQ=255)
                                                                  aaCCAAAACCGCCGCCCGGTGTACCAATCTGCCCTAACATCGCCGCCAGCGTCACCAGCa             <  1:908069/60‑1 (MQ=255)
                                                                  aaCCAAAACCGCCGCCCGGTGTACCAATCTGCCCTAACATCGCCGCCAGCGTCACCAGCATCCAGTGTTTc  <  1:2122035/71‑1 (MQ=255)
                                                                  |                                                                      
CAGACATTTCCGGCAGCACGCCGCCAACACGCGTCGGGTTGCCACCGTTGGAGTAATGATAAGAGAAACCAAAACCGCCGCCCGGTGTACCAATCTGCCCTAACATCGCCGCCAGCGTCACCAGCATCCAGTGTTTC  >  minE/1259699‑1259835

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: