Predicted mutation
evidence seq id position mutation freq annotation gene description
RA minE 1,259,805 C→T 23.7% A373A (GCG→GCA 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,8050CT23.7% 135.0 / 39.7 76A373A (GCG→GCAtorZtrimethylamine N‑oxide reductase system III, catalytic subunit
Reads supporting (aligned to +/- strand):  ref base C (30/28);  new base T (13/5);  total (43/33)
Fisher's exact test for biased strand distribution p-value = 1.75e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 8.78e-01

TTGCCACCGTTGGAGTAATGATAAGAGAAACCAAAACCGCCGCCCGGTGTACCAATCTGCCCTAACATCGCCGCCAGCGTCACCAGCATCCAGTGTTTCTGTTCACCATATTGCTGGCGCTGAATTCCCCAGCCTGC  >  minE/1259737‑1259873
                                                                    |                                                                    
ttaccaccaTTGGCAAAATGGTAAGAAAGACCAAAACCGCCGCCGGGTGTGCCGATTTGCCCCAACATtgc                                                                    >  1:2755281/4‑71 (MQ=255)
ttGCCACCGTTGGAGTAATGATAAGAGAAACCAAAACCGCCGCCCGGTGTACCAATCTGCCCTAACATcgc                                                                    >  1:4433378/1‑71 (MQ=255)
 taccaccaTTGGCAAAATGGTAAGAAAGACCAAAACCGCCGCCGGGTGTGCCGATTTGCCCCAACATtgc                                                                    >  1:239651/3‑70 (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)
                    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)
                          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)
                                          ccgggTGTGCCGATTTGCCCCAACATTGCTGCCAGCGTGACGATCATCCAGTGTTTTTGCTCa                                  >  1:3004314/4‑63 (MQ=255)
                                          cccGGTGTACCAATCTGCCCTAACATCGCCGCCAGCGTCACCAGCATCCAGTGTTTCTGTTCACCATATTg                          >  1:2805669/1‑71 (MQ=255)
                                            cGGTGTACCAATCTGCCCTAACATCGCCGCCAGCGTCACCAGGATCCAGTGTTTCTGTTCACCATATTGCt                        <  1:3908737/71‑1 (MQ=255)
                                            cGGTGTACCAATCTGCCCTAACATCGCCGCCAGCGTCACCAGCATCCAg                                              >  1:1119392/1‑49 (MQ=255)
                                               tgtACCAATCTGCCCTAACATCGCCGCCAGCGTCACCAGCATCCAGTGTTTCTGTTCACCATATTGCTGgc                     >  1:2667483/1‑71 (MQ=255)
                                                 tACCAATCTGCCCTAACATCGCCGCCAGCGTCACCAGCATCCAGTGTTTCTg                                      <  1:1224449/52‑1 (MQ=255)
                                                 tACCAATCTGCCCTAACATCGCCGCCAGCGTCACCAGCATCCAGTGTTTCTg                                      <  1:3432571/52‑1 (MQ=255)
                                                   ccAATCTGCCCTAACATCGCCGCCAGCGTCACCAGCATCCAgtgt                                           >  1:1547495/1‑45 (MQ=255)
                                                   ccAATCTGCCCTAACATCGCCGCCAGCGTCACCAGCATCCAgtgt                                           >  1:2318236/1‑45 (MQ=255)
                                                   ccAATCTGCCCTAACATCGCCGCCAGCGTCACCAGCATCCAgtgt                                           >  1:2839865/1‑45 (MQ=255)
                                                       tttGCCCCAACATTGCTGCCAGCGTGACGATCATCCAGTGTTTTTGCTCACCAAACTGTTGGCGCTGCAt              >  1:4054543/3‑70 (MQ=255)
                                                       tttGCCCCAACATTGCTGCCAGCGTGACGATCATCCAGTGTTTTTGCTCACCAAACTGTTGGCGCTGCAt              >  1:1989388/3‑70 (MQ=255)
                                                       tCTGCCCTAACATCGCCGCCAGCGTCACCAGCATCCAGTGTTTCTGTTCACCATATTGCt                        <  1:148507/60‑1 (MQ=255)
                                                        cTGCCCTAACATCGCCGCCAGCGTCACCAGCATCCAGTGTTTCTGTTCACCATATTGCTGGCGCTGAATTc            >  1:4127653/1‑71 (MQ=255)
                                                        cTGCCCTAACATCGCCGCCAGCGTCACCAGCATCCAGTGTTTCTGTTCACCATATTGCTGGCGCTGAATTc            >  1:761480/1‑71 (MQ=255)
                                                        cTGCCCTAACATCGCCGCCAGCGTCACCAGCATCCAGTGTTTCTGTTCACCATATTGCTGGCGCTGAATTc            >  1:3686733/1‑71 (MQ=255)
                                                         tGCCCTAACATCGCCGCCAGCGTCACCAGCATCCAGTGTTTCTGtt                                    >  1:4328969/1‑46 (MQ=255)
                                                         tGCCCTAACATCGCCGCCAGCGTCACCAGCATCCAGTGTTTCTGtt                                    >  1:402258/1‑46 (MQ=255)
                                                         tGCCCTAACATCGCCGCCAGCGTCACCAGCATCCAGTGTTTCTGtt                                    >  1:3432523/1‑46 (MQ=255)
                                                         tGCCCTAACATCGCCGCCAGCGTCACCAGCATCCAGTGTTTCTGTTCACCATATTGCTGGCGCTGAATTcc           >  1:5063077/1‑71 (MQ=255)
                                                           cccTAACATCGCCGCCAGCGTCACCAGCATCCAGTGTTTCTGTTCACCATATTGCTGGCGCTGAATTcccc         <  1:4124712/71‑1 (MQ=255)
                                                           cccTAACATCGCCGCCAGCGTCACCAGCATCCAGTGTTTCTGTTCACCATATTGCTGGCGCTGAATTcccc         <  1:2038243/71‑1 (MQ=255)
                                                            cccaaCATTGCTGCCAGCGTGACGATCATCCAGTGTTTTTGCTCACCAAACTGTTGGCGCTGCATTccc          >  1:5694394/4‑69 (MQ=255)
                                                            cccaaCATTGCTGCCAGCGTGACGATCATCCAGTGTTTTTGCTCACCAAACTGTTGGCGCTGCATTccc          >  1:1510680/4‑69 (MQ=255)
                                                            cccaaCATTGCTGCCAGCGTGACGATCATCCAGTGTTTTTGCTCACCAAACTGTTGGCGCTGCATTccc          >  1:175182/4‑69 (MQ=255)
                                                            cccaaCATTGCTGCCAGCGTGACGATCATCCAGTGTTTTTGCTCAACAAACTGTTGGCGCTGCATTCCCCa        <  1:5584119/68‑1 (MQ=255)
                                                            ccTAACATCGCCGCCAGCGTCACCAGCATCCAGTGTTTCTGTTCACCATATTGCTGGCGCTGAATTcccc         <  1:3988489/70‑1 (MQ=255)
                                                              tAACATCGCCGCCAGCGTCACCAGCATCCAGTGtt                                          <  1:5567807/35‑1 (MQ=255)
                                                              tAACATCGCCGCCAGCGTCACCAGCATCCAGTGTTTCTGTTCAc                                 <  1:1748286/44‑1 (MQ=255)
                                                              tAACATCGCCGCCAGCGTCACCAGCATCCAGTGTTTCTGTTCACCatat                            <  1:3359860/49‑1 (MQ=255)
                                                              tAACATCGCCGCCAGCGTCACCAGCATCCAGTGTTTCTGTTCACCATATTGCTGGCGCt                  <  1:1665276/59‑1 (MQ=255)
                                                              tAACATCGCCGCCAGCGTCACCAGCATCCAGTGTTTCTGTTCACCATATTGCTGGCGCTGAATTccc          >  1:3949943/1‑67 (MQ=255)
                                                              tAACATCGCCGCCAGCGTCACCAGCATCCAGTGTTTCTGTTCACCATATTGCTGGCGCTGAATTCCCCAGc      >  1:3334757/1‑71 (MQ=255)
                                                               aaCATCGCCGCCAGCGTCACCAGCATCCAGTGTTTCTGTTCACCATATTGCTGGCGCTGAATTCCCCAGcc     >  1:1835596/1‑71 (MQ=255)
                                                                aCATTGCTGCCAGCGTGACGATCATCCAGTGtttt                                        >  1:1916089/1‑34 (MQ=12)
                                                                aCATTGCTGCCAGCGTGACGATCATCCAGTGTTTTTGCTCACCAAACTGTTgg                      >  1:205362/1‑53 (MQ=255)
                                                                aCATTGCTGCCAGCGTGACGATCATCCAGTGTTTTTGCTCACCAAACTGTTgg                      >  1:4864383/1‑53 (MQ=255)
                                                                aCATCGCCGCCAGCGTCACCAGCATCCAGTGTTTCTGTTCACCata                             >  1:5096688/1‑46 (MQ=255)
                                                                aCATCGCCGCCAGCGTCACCAGCATCCAGTGTTTCTGTTCACCata                             >  1:3982914/1‑46 (MQ=255)
                                                                aCATCGCCGCCAGCGTCACCAGCATCCAGTGTTTCTGTTCACCata                             >  1:827036/1‑46 (MQ=255)
                                                                aCATCGCCGCCAGCGTCACCAGCATCCAGTGTTTCTGTTCACCa                               <  1:2800889/44‑1 (MQ=255)
                                                                aCATCGCCGCCAGCGTCACCAGCATCCAGTGTTTCTGTTCACCa                               <  1:1991530/44‑1 (MQ=255)
                                                                aCATCGCCGCCAGCGTCACCAGCATCCAGTGTTTCTGTTCACCATATTg                          >  1:3103883/1‑49 (MQ=255)
                                                                aCATCGCCGCCAGCGTCACCAGCATCCAGTGTTTCTGTTCACCATATTg                          >  1:1661717/1‑49 (MQ=255)
                                                                 cATCGCCGCCAGCGTCACCAGCATCCAGTGTTTCTGTTCACCATATTGCTGGCGCTGAATTCCCCAGCCTg   >  1:2768747/1‑71 (MQ=255)
                                                                  aTCGCCGCCAGCGTCACCAGCATCCAGTGTTTCTGTTCACCa                               <  1:3491325/42‑1 (MQ=255)
                                                                  aTCGCCGCCAGCGTCACCAGCATCCAGTGTTTCTGTTCACCATATTGCTGGCGCTGAATTCCCCAGCCTGc  <  1:5090975/71‑1 (MQ=255)
                                                                    cgccgcCAGCGTCACCAGCATCCAGTGTTTCTGTTCACCata                             >  1:4720874/1‑42 (MQ=255)
                                                                    cgccgcCAGCGTCACCAGCATCCAGTGTTTCTGTTCACCata                             >  1:1491933/1‑42 (MQ=255)
                                                                    cgccgcCAGCGTCACCAGCATCCAGTGTTTCTGTTCACCata                             >  1:2764681/1‑42 (MQ=255)
                                                                    cgccgcCAGCGTCACCAGCATCCAGTGTTTCTGTTCACCata                             >  1:1536825/1‑42 (MQ=255)
                                                                    cgccgcCAGCGTCACCAGCATCCAGTGTTTCTGTTCACCata                             >  1:3638466/1‑42 (MQ=255)
                                                                    |                                                                    
TTGCCACCGTTGGAGTAATGATAAGAGAAACCAAAACCGCCGCCCGGTGTACCAATCTGCCCTAACATCGCCGCCAGCGTCACCAGCATCCAGTGTTTCTGTTCACCATATTGCTGGCGCTGAATTCCCCAGCCTGC  >  minE/1259737‑1259873

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: