Predicted mutation
evidence seq id position mutation freq annotation gene description
RA minE 1,259,805 C→T 7.3% 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,8050CT 7.3% 216.8 / 5.1 68A373A (GCG→GCAtorZtrimethylamine N‑oxide reductase system III, catalytic subunit
Reads supporting (aligned to +/- strand):  ref base C (33/30);  new base T (2/3);  total (35/33)
Fisher's exact test for biased strand distribution p-value = 6.68e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

GCCACCGTTGGAGTAATGATAAGAGAAACCAAAACCGCCGCCCGGTGTACCAATCTGCCCTAACATCGCCGCCAGCGTCACCAGCATCCAGTGTTTCTGTTCACCATATTGCTGGCGCTGAATTCCCCAGCCTGC  >  minE/1259739‑1259873
                                                                  |                                                                    
gCCACCGTTGGAGTAATGATAAGAGAAACCAAAACCGCCGCCCGGTGTACCAATCTGCCCTAACATcgcc                                                                   >  1:2857677/1‑70 (MQ=255)
gCCACCGTTGGAGTAATGATAAGAGAAACCAAAACCGCCGCCCGGTGTACCAATCTGCCCTAACATcgcc                                                                   <  1:3819074/70‑1 (MQ=255)
   aCCGTTGGAGTAATGATAAGAGAAACCAAAACCGCCGCCCTGTGTACCAATCTGCCCTAACATcgccgc                                                                 >  1:4054426/1‑69 (MQ=255)
   aCCGTTGGAGTAATGATAAGAGAAACCAAAACCGCCGCCCGGTGTACCAATCTGCCCTAACATCGCCGCCa                                                               <  1:2978550/71‑1 (MQ=255)
   aCCGTTGGAGTAATGATAAGAGAAACCAAAACCGCCGCCCGGTGTACCAATCTGCCCTAACATCGCCGCCa                                                               <  1:177382/71‑1 (MQ=255)
     cGTTGGAGTAATGATAAGAGAAACCAAAACCGCCGCCCGGTGTACCAATCTGCCCTAACATCGCCGCCAGc                                                             >  1:2220293/1‑71 (MQ=255)
       ttGGAGTAATGATAAGAGAAACCAAAACCGCCGCCCGGTGTACCAATCTGCCCTAACATCGCCGCCAGCGt                                                           <  1:1597373/71‑1 (MQ=255)
          gAGTAATTATAAGAGAAACCAAAACCGCCGCCCGGTGTACCAATCTGCCCTAACATCGCCGCCAGCGTCa                                                         <  1:475474/70‑1 (MQ=255)
              aaTGGTAAGAAAGACCAAAACCGCCGCCGGGTGTGCCGATTTGCCCCAACATTGCTGCCAGCGTGAcgat                                                     >  1:3313564/1‑67 (MQ=255)
               aTGATAAGAGAAACCAAAACCGCCGCCCGGTGTACCAATCTGCCCTAACATCGCCGCCAGCGTCACCAGCa                                                   >  1:2585453/1‑71 (MQ=255)
                   tAAGAGAAACCAAAACCGCCGCCCGGTGTACCAATCTGCCCTAACATCGCCGCCAGCGTCACCAGCATc                                                 >  1:3670793/1‑69 (MQ=255)
                   tAAGAAAGACCAAAACCGCCGCCGGGTGTGCCGATTTGCCCCAACATTGCTGCCAGCGTGAcga                                                      <  1:1886895/64‑3 (MQ=255)
                          aaCCAAAACCGCCGCCCGGTGTACCAATCTGCCCTAACATcgccgc                                                                 <  1:2095070/46‑1 (MQ=255)
                           aCCAAAACCGCCGCCCGGTGTACCAATCTGCCCTAACATCGCCGCCAGCGTCACCAGCATCCAGTGTTTCt                                       >  1:1320390/1‑71 (MQ=255)
                            ccAAAACCGCCGCCCGGTGTACCAATCTGCCCTAACATCgccgcc                                                                <  1:1370525/45‑1 (MQ=255)
                            ccAAAACCGCCGCCCGGTGTACCAATCTGCCCTAACATCgccgcc                                                                <  1:638572/45‑1 (MQ=255)
                            ccAAAACCGCCGCCCGGTGTACCAATCTGCCCTAACATCGCCGCCAGCGt                                                           <  1:2427535/50‑1 (MQ=255)
                             cAAAACCGCCGCCGGGTGTGCCGATTTGCCCCAACATTGCTGCCAGCGTGACGATCATCCAGTGTTTTTgc                                     <  1:2528450/71‑2 (MQ=255)
                             aaaaaCCGCCGCCCGGTGTACCAATCTGGCCTAACATCGCCGCCAGc                                                             <  1:3523507/46‑1 (MQ=255)
                              aaaaCCGCCGCCCGGTGTACCAATCTGCCCTAACATCgccgcc                                                                >  1:937192/1‑43 (MQ=255)
                                aaCCGCCGCCCGGTGTACCAATCTGCCCTAACATCGCCGCCAGCGTCACCAGCATCCAGTGTTTCTGTTCa                                  <  1:3617065/71‑1 (MQ=255)
                                    gccgccCGGTGTACCAATCTGCCCTAACATCGCCGCCAGCGTCACCAGCATCCAGTGTTTCTGTTCACCa                               <  1:1517962/70‑1 (MQ=255)
                                          cGGTGTACCAATCTGCCCTAACATCGCCGCCAGCGTCACCAGCATCCAGTGTTTCTGTTCACCATATTGCt                        <  1:3016714/71‑1 (MQ=255)
                                          cGGTGTACCAATCTGCCCTAACATCGCCGCCAGCGTCACCAGCATCCAGTGTTTCTGTTCACCATATTGCt                        >  1:1710474/1‑71 (MQ=255)
                                           ggTGTACCAATCTGCCCTAACATCGCCGCCAGCGTCACCAGCATCCAGTGTTTCTg                                      <  1:3907139/56‑1 (MQ=255)
                                           ggTGTACCAATCTGCCCTAACATCGCCGCCAGCGTCACCAGCATCCAGTGTTTCTg                                      <  1:1351023/56‑1 (MQ=255)
                                            gtgtACCAATCTGCCCTAACATCGCCGCCAGCGTCACCAGCATCCAGTGTTTCTGTTCACCa                               >  1:1072496/1‑62 (MQ=255)
                                             tgtgCCGATTTGCCCCAACATTGCTGCCAGCGTGACGATCATCCAGTGTTTTTGCTCACCAAACTGTTGgc                     <  1:2356310/71‑1 (MQ=255)
                                              gtACCAATCTGCCCTAACATCGCCGCCAGCGTCACCAGCATCCAGTGTTTCTg                                      >  1:2074946/1‑53 (MQ=255)
                                              gtACCAATCTGCCCTAACATCGCCGCCAGCGTCACCAGCATCCAGTGTTTCTGTTCACCata                             >  1:2754179/1‑62 (MQ=255)
                                                 ccAATCTGCCCTAACATCGCCGCCAGCGTCACCAGCATCCAGTGttt                                         >  1:1906000/1‑47 (MQ=255)
                                                 ccAATCTGCCCTAACATCGCCGCCAGCGTCACCAGCATCCAGTGTTTCTGTTCACCATATTg                          >  1:4101501/1‑62 (MQ=255)
                                                 ccAATCTGCCCTAACATCGCCGCCAGCGTCACCAGCATCCAGTGTTTCTGTTCACCATATTGc                         >  1:86129/1‑63 (MQ=255)
                                                  cAATCTGCCCTAACATCGCCGCCAGCGTCACCAGCATCCAGTGTTTCTGTTCACCATATTGCTGGCGCTg                 <  1:841231/70‑1 (MQ=255)
                                                   aaTCTGCCCTAACATCGCCGCCAGCGTCACCAGCATCCAGTGtt                                          <  1:3706085/44‑1 (MQ=255)
                                                    aTCTGCCCTAACATCGCCGCCAGCGTCACCAGCATCCAGTGTTTCTGTTCACCATATTGCTGGCGCTGAAt              <  1:3498550/71‑1 (MQ=255)
                                                    aTCTGCCCTAACATCGCCCCCAGCGTCACCAGCATCCAGTGTTTCTGTTCACCATATTGCTCGCGCTGa                >  1:2035281/1‑69 (MQ=255)
                                                      ttGCCCCAACATTGCTGCCAGCGTGACGATCATCCAGTGTTTTTGCTCACCAAACTGTTGGCGCTGCAttg            >  1:3651372/2‑70 (MQ=255)
                                                      cTGCCCTAACATCGCCGCCAGCGTCACCAGCATCCAGTGTTTCTGtt                                    >  1:1286498/1‑47 (MQ=255)
                                                      cTGCCCTAACATCGCCGCCAGCGTCACCAGCATCCAGTGTTTCTGTTCACCATATTGCTGGCGCTGAATTc            >  1:949744/1‑71 (MQ=255)
                                                       tGCCCTAACATCGCCGCCAGCGTCACCAGCATCCAGTGttt                                         >  1:1722870/1‑41 (MQ=255)
                                                       tGCCCTAACATCGCCGCCAGCGTCACCAGCATCCAGTGttt                                         >  1:611931/1‑41 (MQ=255)
                                                       tGCCCTAACATCGCCGCCAGCGTCACCAGCATCCAGTGTTTc                                        >  1:2043001/1‑42 (MQ=255)
                                                       tGCCCTAACATCGCCGCCAGCGTCACCAGCATCCAGTGTTTc                                        >  1:3817044/1‑42 (MQ=255)
                                                       tGCCCTAACATCGCCGCCAGCGTCACCAGCATCCAGTGTTTCTGTTCACCATATTGCTGgc                     <  1:1924226/61‑1 (MQ=255)
                                                       tGCCCTAACATCGCCGCCAGCGTCACCAGCATCCAGTGTTTCTGTTCACCATATTGCTGGCGCTGAAtt             <  1:3921336/69‑1 (MQ=255)
                                                       tGCCCTAACATCGCCGCCAGCGTCACCAGCATCCAGTGTTTCTGTTCACCATATTGCTGGCGCTGAATTcc           >  1:676311/1‑71 (MQ=255)
                                                       tGCCCTAACATCGCCGCCAGCGTCACCAGCATCCAGTGTTTCTGTTCACCATATTGCTGGCGCTGAATTcc           >  1:489801/1‑71 (MQ=255)
                                                       tGCCCTAACATCGCCGCCAGCGTCACCAGCATCCAGTGTTTCTGTTCACCATATTGCTGGCGCTGAATTc            >  1:2951145/1‑70 (MQ=255)
                                                         cccTAACATCGCCGCCAGCGTCACCAGCATCCAGTGTTTCTGTTCACCATATTGCTGGCGCTGAATTc            <  1:1958911/68‑1 (MQ=255)
                                                          ccTAACATCGCCGCCAGCGTCACCAGCATCCAGTGTTTCTGTTCACCATATTGCTGGCGCTGAATTCCCCa        <  1:2132639/71‑1 (MQ=255)
                                                          ccTAACATCGCCGCCAGCGTCACCAGCATCCAGTGTTTCTGTTCACCATATTGCTGGCGCTGAATTCCCCa        <  1:1306150/71‑1 (MQ=255)
                                                           cTAACATCGCCGCCAGCGTCACCAGCATCCAGTGTTTCTGTTCACCATATTGCTGGCGCTGAATTCCCCAg       <  1:620954/71‑1 (MQ=255)
                                                           cTAACATCGCCGCCAGCGTCACCAGCATCCAGTGTTTCTGTTCACCATATTGCTGGCGCTGAATTCCCCAg       <  1:825555/71‑1 (MQ=255)
                                                           cTAACATCGCCGCCAGCGTCACCAGCATCCAGTGTTTCTGTTCACCATATTGCTGGCGCTGAATTCCCCAg       >  1:917326/1‑71 (MQ=255)
                                                             aaCATCGCCTCCAGCGTCACCAGCATCCAGTGTTTCTg                                      <  1:1663325/38‑1 (MQ=255)
                                                             aaCATCGCCGCCAGCGTCACCAGCATCCAGTGTTTCTg                                      <  1:3236500/38‑1 (MQ=255)
                                                             aaCATCGCCGCCAGCGTCACCAGCATCCAGTGTTTCTGtt                                    >  1:634628/1‑40 (MQ=255)
                                                             aaCATCGCCGCCAGCGTCACCAGCATCCAGTGTTTCTGtt                                    >  1:985012/1‑40 (MQ=255)
                                                             aaCATCGCCGCCAGCGTCACCAGCATCCAGTGTTTCTGTTCAc                                 >  1:3447520/1‑43 (MQ=255)
                                                             aaCATCGCCGCCAGCGTCACCAGCATCCAGTGTTTCTGTTCACCATATTGCTGGCGCTGaa               >  1:645374/1‑61 (MQ=255)
                                                              aCATCGCCGCCAGCGTCCCCAGCATCCAGTGTTTCTGTTCACCata                             >  1:789654/1‑46 (MQ=255)
                                                              aCATCGCCGCCAGCGTCACCAGCATCCAGTGTTTCTGTTCACCa                               >  1:1615056/1‑44 (MQ=255)
                                                              aCATCGCCGCCAGCGTCACCAGCATCCAGTGTTTCTGTTCACCATATTGCTGGCGCTGaa               >  1:959563/1‑60 (MQ=255)
                                                              aCATCGCCGCCAGCGTCACCAGCATCCAGTGTTTCTGTTCACCATATTGCTGGCGCTGAATTCCCCAg       <  1:3682026/68‑1 (MQ=255)
                                                              aCATCGCCGCCAGCGTCACCAGCATCCAGTGTTTCTGTTCACCATATTGCTGGCGCTGAATTCCCCAGCCt    >  1:3255222/1‑71 (MQ=255)
                                                               cATCGCCGCCAGCGTCACCAGCATCCAGTGTTTCTGt                                     >  1:2809879/1‑37 (MQ=255)
                                                                aTCGCCGCCAGCGTCACCAGCATCCAGTGTTTCTGTTCACCATATTGCTGGCGCTGAATTCCCCAGCCTGc  <  1:2724486/71‑1 (MQ=255)
                                                                aTCGCCGCCAGCGTCACCAGCATCCAGTGTTTCTGTTCACCATATTGCTGGCGCTGAATTCCCCAGCCTGc  <  1:1289206/71‑1 (MQ=255)
                                                                 tCGCCGCCAGCGTCACCAGCATCCAGTGTTTCTGTTCACCa                               <  1:3206601/41‑1 (MQ=255)
                                                                 tCGCCGCCAGCGTCACCAGCATCCAGTGTTTCTGTTCACCa                               <  1:1761356/41‑1 (MQ=255)
                                                                 tCGCCGCCAGCGTCACCAGCATCCAGTGTTTCTGTTCACCATATTGCTGGCGCTGAATTCCCCAGCCTGc  >  1:3702334/1‑70 (MQ=255)
                                                                  |                                                                    
GCCACCGTTGGAGTAATGATAAGAGAAACCAAAACCGCCGCCCGGTGTACCAATCTGCCCTAACATCGCCGCCAGCGTCACCAGCATCCAGTGTTTCTGTTCACCATATTGCTGGCGCTGAATTCCCCAGCCTGC  >  minE/1259739‑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: