Predicted mutation
evidence seq id position mutation freq annotation gene description
RA minE 1,260,368 G→A 8.1% H186Y (CAT→TAT)  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,260,3680GA 8.1% 210.3 / 6.7 74H186Y (CAT→TAT) torZtrimethylamine N‑oxide reductase system III, catalytic subunit
Reads supporting (aligned to +/- strand):  ref base G (34/34);  new base A (3/3);  total (37/37)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.18e-01

ACCTGGCTGTTTTCCAGAATCAGCGGCCAGGAGGTCTGCTGTTCATACACCTCTACAGAACCGACCACATGCGGCATGATCACCTGTGCTGCCCCGGTTGAATAATCGCCGCTATGCCCGGAATAACCGCCCGCCAGGTT  >  minE/1260298‑1260437
                                                                      |                                                                     
aCCTGGCTGTTTTCCAGAATCAGCGGCCAGGAGGTCTGCTGTTCATACACCTCTACAGAACCGACCACATg                                                                       >  1:1607392/1‑71 (MQ=255)
aCCTGGCTGTTTTCCAGAATCAGCGGCCAGGAGGTCTGCTGTTCATACACCTCTACAGAACCGACCACATg                                                                       <  1:1657138/71‑1 (MQ=255)
  cTGGCTGTTTTCCAGAATCAGCGGCCAGGAGGTCTGCTGTTCATACACCTCTACAGAACCGACCACATGCg                                                                     <  1:1819726/71‑1 (MQ=255)
    ggCTGTTTTCCAGAATCAGCGGCCAGGAGGTCTGCTGTTCATACACCTCTACAGAACCGACCACATGCgg                                                                    >  1:5348665/1‑70 (MQ=255)
         ttttCCAGAATCAGCGGCCAGGAGGTCTGCTGTTCATACACCTCTACAGAACCGACCACATGCGGCa                                                                  <  1:4982701/67‑1 (MQ=255)
         ttttCCAGAATCAGCGGCCAGGAGGTCTGCTGTTCATACACCTCTACAGAACCGACCACATGCGGCa                                                                  <  1:3013841/67‑1 (MQ=255)
         ttttCCAGAATCAGCGGCCAGGAGGTCTGCTGTTCATACACCTCTACAGAACCGACCACATGCGGCATGa                                                               >  1:2053805/1‑70 (MQ=255)
         tgttCCAGCACCAGCGGCCAACTGGTCTGCTGTTGATAAACTTCACTACCACCCACGACATACGGCATGa                                                               <  1:3837518/68‑1 (MQ=255)
           ttCCAGCACCAGCGGCCAACTGGTCTGCTGTTGATAAACTTCACTACCACCCACGACATACGGCATGAt                                                              >  1:3448164/1‑69 (MQ=255)
           ttCCAGCACCAGCGGCCAACTGGTCTGCTGTTGATAAACTTCACTACCACCCACGACATACGGCATGATcg                                                            <  1:3771724/71‑2 (MQ=255)
           ttCCAGCACCAGCGGCCAACTGGTCTGCTGTTGATAAACTTCACTACCACCCACGACATACGGCATGATcg                                                            <  1:5635217/71‑2 (MQ=255)
           ttCCAGCACCAGCGGCCAACTGGTCTGCTGTTGATAAACTTCACTACCACCCACGACATACGGCATGATcg                                                            >  1:4888777/1‑70 (MQ=255)
           ttCCAGAATCAGCGGCCAGGAGGTCTGCTGTTCATACACCTCTACAGAACCGACCACATGCGGCATg                                                                >  1:3667326/1‑67 (MQ=255)
           ttCCAGAATCAGCGGCCAGGAGGTCTGCTGTTCATACACCTCTACAGAACCGACCACATGCGGCATg                                                                >  1:5605029/1‑67 (MQ=255)
           ttCCAGAATCAGCGGCCAGGAGGTCTGCAGTTCATACACCTCTACAGAACCGACCACATGCGGCATg                                                                >  1:1734160/1‑67 (MQ=255)
             ccAGAATCAGCGGCCAGGAGGTCTGCTGTTCATACACCTCTACAGAACCGACCACATGCGGCATGATc                                                             <  1:1905987/68‑1 (MQ=255)
             ccAGAATCAGCGGCCAGGAGGTCTGCTGTTCATACACCTCTACAGAACCGACCACATGCGGCATGATc                                                             <  1:4622432/68‑1 (MQ=255)
             ccAGAATCAGCGGCCAGGAGGTCTGCTGTTCATACACCTCTACAGAACCGACCACATGCGGCATGATCAcc                                                          >  1:2463294/1‑71 (MQ=255)
             ccAGAATCAGCGGCCAGGAGGTCTGCTGTTCATACACCTCTACAGAACCGACCACATGCGGCATGATCAc                                                           <  1:3754421/70‑1 (MQ=255)
             ccAGAATCAGCGGCCAGGAGGTCTGCTGTTCATACACCTCTACAGAACCGACCACATGCGGCATGATCAc                                                           <  1:4802583/70‑1 (MQ=255)
             ccAGAATCAGCGGCCAGGAGGTCTGCTGTTCATACACCTCTACAGAACCGACCACATGCGGCATGATCAc                                                           >  1:1843578/1‑70 (MQ=255)
             ccAGAATCAGCGGCCAGGAGGTCTGCTGTTCATACACCTCTACAGAACCGACCACATGCGGCATGATCAc                                                           <  1:65811/70‑1 (MQ=255)
                 aaTCAGCGGCCAGGAGGTCTGCTGTTCATACACCTCTACAGAACCGACCACATGCGGCATGATCACCtgtg                                                      >  1:636680/1‑71 (MQ=255)
                 aaTCAGCGGCCAGGAGGTCTGCTGTTCATACACCTCTACAGAACCGACCACATGCGGCATGATCACCtgtg                                                      >  1:2739886/1‑71 (MQ=255)
                   ccAGCGGCCAGGAGGTTTGCTGTTCATACACTTCCATCGAACCGACTACGCGCGGCAGGATCAc                                                           >  1:2469112/2‑64 (MQ=255)
                   ccAGCGGCCAGGAGGTTTGCTGTTCATACACTTCCATCGAACCGACTACGCGCGGCAGGATCAc                                                           >  1:4428977/2‑64 (MQ=255)
                   ccAGCGGCCAGGAGGTTTGCTGTTCATACACTTCCATCGAACCGACTACGCGCGGCAGGATCACCtgcg                                                      >  1:2004066/2‑67 (MQ=255)
                       cGGCCAGGAGGTCTGCTGTTCATACACCTCTACAGAACCGACCACATGCGGCATGATCACCtgt                                                       <  1:1127941/64‑1 (MQ=255)
                       cGGCCAGGAGGTCTGCTGTTCATACACCTCTACAGAACCGACCACATGCGGCATGATCACCtgt                                                       <  1:2812304/64‑1 (MQ=255)
                       cGGCCAGGAGGTCTGCTGTTCATACACCTCTACAGAACCGACCACATGCGGCATGATCACCTGTGCTGccc                                                <  1:3558249/71‑1 (MQ=255)
                         gCCAGGAGGTCTGCTGTTCATACACCTCTACAGAACCGACCACATGCGGCATGATCACCTGTGCTGCCCCg                                              >  1:3156104/1‑71 (MQ=255)
                         gCCAGGAGGTCTGCTGTTCATACACCTCTACAGAACCGACCACATGCGGCATGATCACCTGTGCTGCCCCg                                              >  1:2095323/1‑71 (MQ=255)
                           cAGGAGGTTTGCTGTTCATACACTTCCATCGAACCGACTACGCGCGGCAGGATCACCAGCGCAGCACCgg                                             >  1:1902837/1‑70 (MQ=255)
                            aggaggTCTGCTGTTCATACACCTCTACAGAACCGACCACATGCGGCATGATCACCTGTGCTGCCCCGGtt                                           >  1:2762715/1‑71 (MQ=255)
                              gaggTCTGCTGTTCATACACCTCTACAGAACCGACCACATGCGGCATGATCAc                                                           >  1:2985231/1‑53 (MQ=255)
                               tggTCTGCTGTTGATAAACTTCACTACCACCCACGACATACGGCATGAt                                                              >  1:4584112/2‑49 (MQ=255)
                                ggTTTGCTGTTCATACACTTCCATCGAACCGACTAAGCGCGGCAGGATCACCTGCGCAGCACAGGTAGag                                        >  1:866146/1‑69 (MQ=255)
                                  tCTGCTGTTCATACACCTCTACAGAACCGACCACATGCGGCATGATCACCTGTGCTGCCCCGGTTGaataa                                     <  1:3517141/71‑1 (MQ=255)
                                   ctgctgTTCATACACCTCTACAGAACCGACCACATGCGGCATGATCACCTGTGCTGCCCCGGTTGaataat                                    >  1:5492888/1‑71 (MQ=255)
                                   ctgctgTTCATACACCTCTACAGAACCGACCACATGCGGCATGATCACCTGTGCTGCCCCGGTTGaataat                                    >  1:2371063/1‑71 (MQ=255)
                                    tgctgTTCATACACTTCCATCGAACCGACTACGCGCGGCAGGATCACCTGCGCAGCACCGGTAGagtaat                                    <  1:5009371/70‑1 (MQ=255)
                                    tgctgTTCATACACCTCTACAGAACCGACCACATGCGGCATGATc                                                             <  1:2666496/45‑1 (MQ=255)
                                    tgctgTTCATACACCTCTACAGAACCGACCACATGCGGCATGATCACCTGTGCTGCCCCGGTTGAATAATc                                   >  1:3645477/1‑71 (MQ=255)
                                       tgTTCATACACCTCTACAGAACCGACCACATGCgg                                                                    <  1:702544/35‑1 (MQ=255)
                                       tgTTCATACACCTCTACAGAACCGACCACATGCGGCATGATCACCTGTGCTGCCCCGGTTGAATAATcgcc                                >  1:181128/1‑71 (MQ=255)
                                         ttCATACACCTCTACAGCACCGACCACATGCGGCATGATCACCTGTGCTGCCCCGGTTGaat                                       <  1:4559352/62‑1 (MQ=255)
                                         ttCATACACCTCTACAGAACCGACCACATGCGGCATGATCACCTGTGCTGCCCCGGTTGaat                                       <  1:3231834/62‑1 (MQ=255)
                                         ttCATACACCTCTACAGAACCGACCACATGCGGCATGATCACCTGTGCTGCCCCGGTTGAATAATcgccgc                              <  1:1894208/71‑1 (MQ=255)
                                           cATACACCTCTACAGAACCGACCACATGCGGCATGATCACCTGTGCTGCCCCGGTTGaataa                                     >  1:4765835/1‑62 (MQ=255)
                                           cATACACCTCTACAGAACCGACCACATGCGGCATGATCACCTGTGCTGCCCCGGTTGaataa                                     >  1:440420/1‑62 (MQ=255)
                                             tACACCTCTACAGAACCGACCACATGCGGCATGATCa                                                            <  1:1579972/37‑1 (MQ=255)
                                             tACACCTCTACAGAACCGACCACATGCGGCATGATCACCTGtgctgc                                                  <  1:4105586/47‑1 (MQ=255)
                                               cacCTCTACAGAACCGACCACATGCGGCATGATCACCTGTGCTGccc                                                >  1:1769775/1‑47 (MQ=255)
                                                 ccTCTACAGAACCGACCACATGCGGCATGATCACCTGTGCTGCCCCGGTTGAATAATCGCCGCTATGCCCg                      <  1:2810837/71‑1 (MQ=255)
                                                 ccTCTACAGAACCGACCACATGCGGCATGATCACCTGTGCTGCCCCGGTTGAATAATCGCCGCTATGCCCg                      <  1:5546190/71‑1 (MQ=255)
                                                  ctctACAGAACCGACCACATGCGGCATGATCACCTGTGCTGCCCCGGTTGAATAATCGCCGCTATGCCCg                      >  1:2166735/1‑70 (MQ=255)
                                                   tctACAGAACCGACCACATGCGGCATGATCACCtg                                                        <  1:5660438/35‑1 (MQ=255)
                                                    ctACAGAACCGACCACATGCGGCATGATCACCTGTGCTGCCCCGGTTGAATAATCGCCGCTATGCCCg                      <  1:3560212/68‑1 (MQ=255)
                                                    ctACAGAACCGACCACATGCGGCATGATCACCTGTGCTGCCCCGGTTGAATAATCGCCGCTATGCCCg                      <  1:949837/68‑1 (MQ=255)
                                                    ctACAGAACCGACCACATGCGGCATGATCACCTGTGCTGCCCCGGTTGAATAATCGCCGCTATGCCCGGaa                   >  1:4810591/1‑71 (MQ=255)
                                                     tACAGAACCGACCACATGCGGCATGATCACCtgtg                                                      >  1:167511/1‑35 (MQ=255)
                                                      actgAACCGACCACATGCGGCATGATCACCTGTGCTGCCCCGGTTGAATAATCGCCGCTATGCCCGGAAt                  >  1:4850669/4‑70 (MQ=255)
                                                       cAGAACCGACCACATGCGGCATGATCACCTGTGCTGCCCCGGTTGAATAATCGCCGCTATGCCCGGAAt                  <  1:1682495/69‑1 (MQ=255)
                                                       cAGAACCGACCACATGCGGCATGATCACCTGTGCTGCCCCGGTTGAATAATCGCCGCTATGCCCGGAAt                  <  1:2134174/69‑1 (MQ=255)
                                                       cAGAACCGACCACATGCGGCATGATCACCTGTGCTGCCCCGGTTGAATAATCGCCGCTATGCCCGGAATaa                <  1:4300170/71‑1 (MQ=255)
                                                       cAGAACCGACCACATGCGGCATGATCACCTGTGCTGCCCCGGTTGAATAATCGCCGCTATGCCCGGAATaa                >  1:5338316/1‑71 (MQ=255)
                                                       cAGAACCGACCACATGCGGCATGATCACCTGTGCTGCCCCGGTTGAATAATCGCCGCTATGCCCGGAATa                 <  1:4662064/70‑1 (MQ=255)
                                                       cAGAACCGACCACATGCGGCATGATCACCTGTGCTGCCCCGGTTGAATAATCGCCGCTATGCCCGGAATa                 <  1:1919190/70‑1 (MQ=255)
                                                       cAGAACCGACCACATGCGGCATGATCACCTGTGCTGCCCCGGTTGAATAATCGCCGCTATGCCCGGAATa                 <  1:1350277/70‑1 (MQ=255)
                                                        aGAACCGACCACATGCGGCATGATCACCTGTGCTGCCCCGGTTGAATAATcgccgc                              <  1:3756124/56‑1 (MQ=255)
                                                         gAACCGACCACATGCGGCATGATCACCTGTGCTGCCCCGGTTGaat                                       >  1:5011626/1‑46 (MQ=255)
                                                         gAACCGACCACATGCGGCATGATCACCTGTGCTGCCCCGGTTGAATAATCGCCGCTATGCCCGGAATaa                >  1:5476154/1‑69 (MQ=255)
                                                                ccACATGCGGCATGATCACCTGTGCTGCCCCGGTTGAATAATCGCCGCTATGCCCGGAATAAccgcccgc        <  1:1874402/70‑1 (MQ=255)
                                                                ccACATGCGGCATGATCACCTGTGCTGCCCCGGTTGAATAATCGCCGCTATGCCCGGAATAACcgcccgcc       <  1:2946147/71‑1 (MQ=255)
                                                                   caTGCGGCATGATCACCTGTGCTGCCCCGGTTGAATAATCGCCGCTATGCCCGGAATAACCGCCCGCCAgg    >  1:4516885/1‑71 (MQ=255)
                                                                     tGCGGCATGATCACCTGTGCTGCCCCGGTTGAATAATCGCCGCTATGCCCGGAATAACCGCCCGCCAGGtt  >  1:5319022/1‑71 (MQ=255)
                                                                      |                                                                     
ACCTGGCTGTTTTCCAGAATCAGCGGCCAGGAGGTCTGCTGTTCATACACCTCTACAGAACCGACCACATGCGGCATGATCACCTGTGCTGCCCCGGTTGAATAATCGCCGCTATGCCCGGAATAACCGCCCGCCAGGTT  >  minE/1260298‑1260437

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: