Predicted mutation
evidence seq id position mutation freq annotation gene description
RA minE 1,259,859 A→C 20.3% I355M (ATT→ATG 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,8590AC20.3% 104.0 / 32.4 64I355M (ATT→ATGtorZtrimethylamine N‑oxide reductase system III, catalytic subunit
Reads supporting (aligned to +/- strand):  ref base A (17/34);  new base C (9/4);  total (26/38)
Fisher's exact test for biased strand distribution p-value = 2.72e-02
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

ATCTGCCCTAACATCGCCGCCAGCGTCACCAGCATCCAGTGTTTCTGTTCACCATATTGCTGGCGCTGAATTCCCCAGCCTGCCATCAGCATTGTGCGATTCGCCGCCATCAGTTCAGCGAGTTTGACTATCTGCGCT  >  minE/1259791‑1259928
                                                                    |                                                                     
atttGCCCCAACATTGCTGCCAGCGTGACGATCATCCAGTGTTTTTGCTCACCAAACTGTTGGCGCTGCAt                                                                     >  1:1223111/4‑71 (MQ=255)
      ccTAACATCGCCGCCAGCGTCACCAGCATCCAGTGTTTCTGTTCACCATATTGCTGGCGCTGAATTCCCCa                                                               <  1:843594/71‑1 (MQ=255)
      ccTAACATCGCCGCCAGCGTCACCAGCATCCAGTGTTTCTGTTCACCATATTGCTGGCGCTGAATTCCCCa                                                               <  1:3514365/71‑1 (MQ=255)
        caaCATTGCTGCCAGCGTGACGATCATCCAGTGTTTTTGCTCACCAAACTGTTGGCGCTGCATTCCCCAGc                                                             >  1:3033525/2‑71 (MQ=255)
          aCATTGCTGCCAGCGTGACGATCATCCAGTGTTTTTGCTCACCAAACTGTTGGCGCTGCATTCCCCAGCCt                                                           >  1:2520160/1‑71 (MQ=255)
           cATCGCCGCCAGCGTCACCAGCATCCAGTGTTTCTGTTCACCATATTGCTGGCGCTGAATTCCCCAGCCTg                                                          <  1:2655954/71‑1 (MQ=255)
               gctgccAGCGTGACGATCATCCAGTGTTTTTGCTCACCAAACTGTTGGCGCTGCATTCCCCAGCCTGCCa                                                       <  1:3094874/67‑1 (MQ=255)
               gccgccAGCGTCACCAGCATCCAGTGTTTCTGTTCACCATATTGCTGGCGCTGAATTCCCCAGCCTGCCAt                                                      >  1:2221078/1‑71 (MQ=255)
               gccgccAGCGTCACCAGCATCCAGTGTTTCTGTTCACCATATTGCTGGCGCTGAATTCCCCAGCCTGCCAt                                                      >  1:2216343/1‑71 (MQ=255)
                ctgccAGCGTGACGATCATCCAGTGTTTTTGCTCACCAAACTGTTGGCGCTGCATTCCCCAGCCTGCCAt                                                      >  1:3379490/3‑70 (MQ=255)
                ctgccAGCGTGACGATCATCCAGTGTTTTTGCTCACCAAACTGTTGGCGCTGCATTCCCCAGCCTGCCAt                                                      <  1:2241161/68‑1 (MQ=255)
                  gccAGCGTGACGATCATCCAGTGTTTTTGCTCACCAAACTGTTGGCGCTGCATTCCCCAGCCTGCCATCAg                                                   >  1:1137738/1‑71 (MQ=255)
                  gccAGCGTCACCAGCATCCAGTGTTTCTGTTCACCATATTGCTGGCGCTGAATTCCCCAGCCTGCCATCAg                                                   >  1:2827844/1‑71 (MQ=255)
                   ccAGCCTCACCAGCATCCAGTGTTTCTGTTCACCATATTGCTGGCGCTGAATTCCCCAGCCTg                                                          <  1:1938514/63‑1 (MQ=255)
                    cAGCGTCACCAGCATCCAGTGTTTCTGTTCACCATATTGCTGGCGCTGAATTCCCCAGCCTGCCATCAGCa                                                 <  1:462589/71‑1 (MQ=255)
                    cAGCGTCACCAGCATCCAGTGTTTCTGTTCACCATATTGCTGGCGCTGAATTCCCCAGCCTGCCATCAGCa                                                 >  1:3621450/1‑71 (MQ=255)
                    cAGCGTCACCAGCATCCAGTGTTTCTGTTCACCATATTGCTGGCGCTGAATTCCCCAGCCTGCCATCAGCa                                                 <  1:3841247/71‑1 (MQ=255)
                    cAGCGTCACCAGCATCCAGTGTTTCTGTTCACCATATTGCTGGCGCTGAATTCCCCAGCCTGCCATCAGCa                                                 <  1:235307/71‑1 (MQ=255)
                    cAGCGTCACCAGCATCCAGTGTTTCTGTTCACCATATTGCTGGCGCTGAATTCCCCAGCCTGCCATCAGCa                                                 <  1:1507521/71‑1 (MQ=255)
                    cAGCGTCACCAGCATCCAGTGTTTCTGTTCACCATATTGCTGGCGCTGAATTCCCCAGCCTGCCATCAGCa                                                 <  1:1910719/71‑1 (MQ=255)
                     aGCGTCACCAGCATCCAGTGTTTCTGTTCACCATATTGCTGGCGCTGAATTCCCCAGCCTGCCATCAGCAt                                                >  1:792378/1‑71 (MQ=255)
                      gCGTGACGATCATCCAGTGTTTTTGCTCACCAAACTGTTGGCGCTGCATTCCCCAGCCTGCCATCAGCAtg                                               <  1:480872/71‑2 (MQ=255)
                      gCGTGACGATCATCCAGTGTTTTTGCTCACCAAACTGTTGGCGCTGCATTCCCCAGCCTGCCATCAGCAtg                                               <  1:1881162/71‑2 (MQ=255)
                       cGTCACCAGCATCCAGTGTTTCTGTTCACCATATTGCTGGCGCTGAATTCCCCAGCCTGCCATCAGCATtg                                              >  1:1084140/1‑71 (MQ=255)
                        gTCACCAGCATCCAGTGTTTCTGTTCACCATATTGCTGGCGCTGAATTCCCCAGCCTGCCATCAGCATtgt                                             <  1:415203/71‑1 (MQ=255)
                        gTCACCAGCATCCAGTGTTTCTGTTCACCATATTGCTGGCGCTGAATTCCCCAGCCTGCCATCAGCATtg                                              <  1:371696/70‑1 (MQ=255)
                          cACCAGCATCCAGTGTTTCTGTTCACCATATTGCTGGCGCTGAATTCCCCAGCCTGCCATCAGCATTGTGc                                           >  1:2726561/1‑71 (MQ=255)
                              aGCATCCAGTGTTTCTGTTCACCATATTGCTGGCGCTGAATTCCCCAGCCTGCCATCAGCATTGTGCGAtt                                       >  1:3160833/1‑71 (MQ=255)
                               gCATCCAGTGTTTCTGTTCACCATATTGCTGGCGCTGAATTCCCCAGCCTGCCATCAGCATTGTGCGATTc                                      <  1:2928067/71‑1 (MQ=255)
                               gCATCCAGTGTTTCTGTTCACCATATTGCTGGCGCTGAATTCCCCAGCCTGCCATCAGCATTGTGCGATTc                                      >  1:3191134/1‑71 (MQ=255)
                                 aTCCAGTGTTTCTGTTCACCATATTGCTGGCGCTGAATTCCCCAGCCTGCCATCAGCATTGTGCGATTcgc                                    <  1:802895/71‑1 (MQ=255)
                                  tCCAGTGTTTCTGTTCACCATATTGCTGGCGCTGAATTCCCCAGCCTGCCATCAGCATTGTGCGATTcgcc                                   <  1:715872/71‑1 (MQ=255)
                                    cAGTGTTTTTGCTCACCAAACTGTTGGCGCTGCATTCCCCAGCCTGCCATc                                                     >  1:1671019/1‑51 (MQ=255)
                                    cAGTGTTTTTGCTCACCAAACTGTTGGCGCTGCATTCCCCAGCCTGCCATc                                                     >  1:3445555/1‑51 (MQ=255)
                                    cAGTGTTTTTGCTCACCAAACTGTTGGCGCTGCATTCCCCAGCCTGCCATc                                                     >  1:1226739/1‑51 (MQ=255)
                                    cAGTGTTTTTGCTCACCAAACTGTTGGCGCTGCATTCCCCAGCCTGCCATc                                                     >  1:1240331/1‑51 (MQ=255)
                                      gtgtTTCTGTTCACCATATTGCTGGCGCTGAATTcccc                                                                <  1:352911/38‑1 (MQ=255)
                                      gtgtTTCTGTTCACCATATTGCTGGCGCTGAATTCCCCAGCCTGCCATCAGCATTGTGCGAt                                        <  1:1949623/62‑1 (MQ=255)
                                       tgtTTCTGTTCACCATATTGCTGGCGCTGAATTCCCCAGCCTGCCATCAGCATTGTGCGATTc                                      <  1:3684044/63‑1 (MQ=255)
                                       tgtTTCTGTTCACCATATTGCTGGCGCTGAATTCCCCAGCCTGCCATCAGCATTGTGCGATTCGCCGCCAt                              <  1:1407074/71‑1 (MQ=255)
                                       tgtTTCTGTTCACCATATTGCTGGCGCTGAATTCCCCAGCCTGCCATCAGCATTGTGCGATTCGCCGCCAt                              <  1:1126505/71‑1 (MQ=255)
                                        gtTTCTGTTCACCATATTGCTGGCGCTGAATTCCCCAGCCTGGCATCAGCATTGTGCGATTCGCCGCCAt                              >  1:34580/1‑70 (MQ=255)
                                        gtTTCTGTTCACCATATTGCTGGCGCTGAATTCCCCAGCCTGCCATCAGCATTGTGc                                           <  1:494234/57‑1 (MQ=255)
                                         tttCTGTTCACCATATTGCTGGCGCTGAATTCCCCAGCCTGCCATCAGCATTGTGCGATTcgccgc                                 <  1:295695/66‑1 (MQ=255)
                                         tttCTGTTCACCATATTGCTGGCGCTGAATTCCCCAGCCTGCCATCAGCATTGTGCGATTCgccgcc                                <  1:3369877/67‑1 (MQ=255)
                                         tttCTGTTCACCATATTGCTGGCGCTGAATTCCCCAGCCTGCCATCAGCATTGTGCGATTCGCCGCCATc                             <  1:1356276/70‑1 (MQ=255)
                                         tttCTGTTCACCATATTGCTGGCGCTGAATTCCCCAGCCTGCCATCAGCATTGTGCGATTCGCCGCCATc                             <  1:3370159/70‑1 (MQ=255)
                                           tCTGTTCACCATATTGCTGGCGCTGAATTCCCCAGCCTGCCATCAGCATTGTGCGATTCGCCGCCATCAg                           <  1:57242/70‑1 (MQ=255)
                                            cTGTTCACCATATTGCTGGCGCTGAATTCCCCAGCCTGCCATCAGCATTGTGCGATTCGCCGCCATCAGtt                         >  1:310941/1‑71 (MQ=255)
                                            cTGTTCACCATATTGCTGGCGCTGAATTCCCCAGCCTGCCATCAGCATTGTGCGATTCGCCGCCATCAGtt                         >  1:3924097/1‑71 (MQ=255)
                                                 cACCATATTGCTGGCGCTGAATTCCCCAGCCTGCCATCAGCATTGTGCGATTCGCCGCCATCAGTTCAGc                     <  1:3130496/70‑1 (MQ=255)
                                                  aCCATATTGCTGGCGCTGAATTCCCCAGCCTGCCATCAGCATTGTGCGATTCGCCGCCATCAGTTCAGCGa                   <  1:2943778/71‑1 (MQ=255)
                                                   ccATATTGCTGGCGCTGAATTCCCCAGCCTGCCATCAGCATTGTGCGATTCGCCGCCATCAGTTCAGCGAg                  <  1:2178134/71‑1 (MQ=255)
                                                    cATATTGCTGGCGCTGAATTCCCCAGCCTGCCATCAGCATTGTGCGATTCGCCGCCATCAGTTCAGCGAGt                 >  1:3100847/1‑71 (MQ=255)
                                                      tatTGCTGGCGCTGAATTCCCCAGCCTGCCATCAGCATTGTGCGATTCGCCGCCATCAGTTCAGCGAGttt               <  1:3767300/71‑1 (MQ=255)
                                                         tGCTGGCGCTGAATTCCCCAGCCTGCCATCAGCATTGTGCGATTCGCCGCCAt                              >  1:1236801/1‑53 (MQ=255)
                                                         tGCTGGCGCTGAATTCCCCAGCCTGCCATCAGCATTGTGCGATTCGCCGCCAACAGTTCAGCGAGTTTGAc            >  1:1866554/1‑71 (MQ=255)
                                                          gCTGGCGCTGAATTCCCCAGCCTGCCATCAGCATTGTGCGATTCGCCGCCATCAGTTCAGCGAGTTTGAc            <  1:780591/70‑1 (MQ=255)
                                                              gcgcTGAATTCCCCAGCCTGCCATCAGCATTGTGCGATTCTCCGCCATCAGTTCAGCGAGTTTGACTATCt       <  1:816315/71‑1 (MQ=255)
                                                              gcgcTGAATTCCCCAGCCTGCCATCAGCATTGTGCGATTCGCCGCCATCAGTTCAGCGAGTTTGACTAt         <  1:3090433/69‑1 (MQ=255)
                                                                gcTGAATTCCCCAGCCTGCCATCAGCATTGTGCGATTCGCCGCCATCAGTTCAGCGAGTTTGACTATCTgc     <  1:1413161/71‑1 (MQ=255)
                                                                 cTGAATTCCCCAGCCTGCCATCAGCATTGTGCGATTCGCCGCCATCAGTTCAGCGAGTTTGACTATCTgcg    <  1:1911283/71‑1 (MQ=255)
                                                                   gAATTCCCCAGCCTGCCATCAGCATTGTGCGATTCGCCGCCATCAGTTCAGCGAGTTTGa             >  1:2623334/1‑60 (MQ=255)
                                                                   gAATTCCCCAGCCTGCCATCAGCATTGTGCGATTCGCCGCCATCAGTTCAGCGAGTTTGACTATCTGCGCt  >  1:720879/1‑71 (MQ=255)
                                                                    aaTTCCCCAGCCTGCCATCAGCATTGTGCGATTCGCCGCCATCAGTTCAGc                     <  1:842486/51‑1 (MQ=255)
                                                                    aaTTCCCCAGCCTGCCATCAGCATTGTGCGATTCGCCGCCATCAGTTCAGCGAg                  <  1:3915133/54‑1 (MQ=255)
                                                                    |                                                                     
ATCTGCCCTAACATCGCCGCCAGCGTCACCAGCATCCAGTGTTTCTGTTCACCATATTGCTGGCGCTGAATTCCCCAGCCTGCCATCAGCATTGTGCGATTCGCCGCCATCAGTTCAGCGAGTTTGACTATCTGCGCT  >  minE/1259791‑1259928

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: