Predicted mutation
evidence seq id position mutation freq annotation gene description
RA minE 1,087,363 G→A 54.1% I284I (ATC→ATT sufS ← selenocysteine lyase, PLP‑dependent

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*minE1,087,3630GA54.1% 3.6 / 50.3 61I284I (ATC→ATTsufSselenocysteine lyase, PLP‑dependent
Reads supporting (aligned to +/- strand):  ref base G (13/15);  new base A (21/12);  total (34/27)
Fisher's exact test for biased strand distribution p-value = 2.05e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 8.45e-01

CATACTCGGCTATGTTATTAAGCCCCAGCGCCGAAACATACTCCAGCGCCGCGCCAAGACCAATGATGCCCCCGGTATTGGGTGTACCGGCTTCAAACCGCCATGGTGCTTTGGTCCAGGTAGTGCCTTCA  >  minE/1087299‑1087429
                                                                |                                                                  
cATACTCGGCTATGTTATTAAGCCCCAGCGCCGAAACATACTCCAGCGCCGCGCCAATACCAATGATGc                                                                >  1:2288325/1‑69 (MQ=255)
    cTCGGCTATGTTATTAAGCCCCAGCGCCGAAACATACTCCAGCGCCGCGCCAAGACCAATGATGccccc                                                            >  1:1260866/1‑69 (MQ=255)
    cTCGGCTATGTTATTAAGCCCCAGCGCCGAAACATACTCCAGCGCCGCGCCAAGACCAATGATGccccc                                                            >  1:1702736/1‑69 (MQ=255)
      cGGCAATGTTGTTAAGCCCCAGCGCCGAAACATACTCCAGCGCCGCGCCAAGGCCAATAATACCTCCgg                                                          <  1:967305/69‑1 (MQ=255)
      cGGCAATGTTGTTAAGCCCCAGCGCCGAAACATACTCCAGCGCCGCGCCAAGGCCAATAATACCTCCgg                                                          <  1:2885617/69‑1 (MQ=255)
      cGGCAATGTTGTTAAGCCCCAGCGCCGAAACATACTCCAGCGCCGCGCCAAGGCCAATAATACCTCCgg                                                          <  1:3171479/69‑1 (MQ=255)
      cGGCAATGTTGTTAAGCCCCAGCGCCGAAACATACTCCAGCGCCGCGCCAAGGCCAATAATACCTCCgg                                                          <  1:2581662/69‑1 (MQ=255)
      cGGCAATGTTGTTAAGCCCCAGCGCCGAAACATACTCCAGCGCCGCGCCAAGGCCAATAATACCTCCgg                                                          <  1:2560684/69‑1 (MQ=255)
      cGGCAATGTTGTTAAGCCCCAGCGCCGAAACATACTCCAGCGCCGCGCCAAGGCCAATAATACCTCCgg                                                          <  1:2207943/69‑1 (MQ=255)
      cGGCAATGTTGTTAAGCCCCAGCGCCGAAACATACTCCAGCGCCGCGCCAAGGCCAATAATACCTCCgg                                                          <  1:2167880/69‑1 (MQ=255)
      cGGCAATGTTGTTAAGCCCCAGCGCCGAAACATACTCCAGCGCCGCGCCAAGGCCAATAATACCTCCgg                                                          <  1:456891/69‑1 (MQ=255)
      cGGCAATGTTGTTAAGCCCCAGCGCCGAAACATACTCCAGCGCCGCGCCAAGGCCAATAATACCTCCgg                                                          <  1:187646/69‑1 (MQ=255)
      cGGCAATGTTGTTAAGCCCCAGCGCCGAAACATACTCCAGCGCCGCGCCAAGGCCAATAATACCTCCgg                                                          <  1:1845643/69‑1 (MQ=255)
      cGGCAATGTTGTTAAGCCCCAGCGCCGAAACATACTCCAGCGCCGCGCCAAGGCCAATAATACCTCCgg                                                          <  1:1300905/69‑1 (MQ=255)
      cGGCAATGTTGTTAAGCCCCAGCGCCGAAACATACTCCAGCGCCGCGCCAAGGCCAATAATACCTCCgg                                                          <  1:1736633/69‑1 (MQ=255)
         cTATGTTATTAAGCCCCAGCGCCGAAACATACTCCAGCGCCGCGCCAAGACCAATGATGCCCCCGGt                                                         <  1:3201161/67‑1 (MQ=255)
           aTGTTATTAAGCCCCAGCGCCGAAACATACTCCAGCGCCGCGCCAAGACCAATGATGCCCCCGGTATTg                                                     <  1:2471674/69‑1 (MQ=255)
           aTGTTATTAAGCCCCAGCGCCGAAACATACTCCAGCGCCGCGCCAAGACCAATGATGCCCCCGGTATTg                                                     <  1:2641708/69‑1 (MQ=255)
            tgttgtTAAGCCCCAGCGCCGAAACATACTCCAGCGCCGCGCCAAGGCCAATAATACCTCCGGTGTTgg                                                    >  1:3261112/1‑69 (MQ=255)
            tgttgtTAAGCCCCAGCGCCGAAACATACTCCAGCGCCGCGCCAAGGCCAATAATACCTCCGGTATTgg                                                    >  1:1905795/1‑69 (MQ=255)
            tgttgtTAAGCCCCAGCGCCGAAACATACTCCAGCGCCGCGCCAAGGCCAATAATACCTCCGGTATTgg                                                    >  1:2365854/1‑69 (MQ=255)
            tgttgtTAAGCCCCAGCGCCGAAACATACTCCAGCGCCGCGCCAAGGCCAATAATACCTCCGGTATTgg                                                    >  1:1849093/1‑69 (MQ=255)
            tgttgtTAAGCCCCAGCGCCGAAACATACTCCAGCGCCGCGCCAAGGCCAATAATACCTCCGGTATTgg                                                    >  1:80228/1‑69 (MQ=255)
            tgttgtTAAGCCCCAGCGCCGAAACATACTCCAGCGCCGCGCCAAGGCCAATAATACCTCCGGTATTgg                                                    >  1:546436/1‑69 (MQ=255)
            tgttgtTAAGCCCCAGCGCCGAAACATACTCCAGCGCCGCGCCAAGGCCAATAATACCTCCGGTATTgg                                                    >  1:2127022/1‑69 (MQ=255)
            tgttgtTAAGCCCCAGCGCCGAAACATACTCCAGCGCCGCGCCAAGGCCAATAATACCTCCGGTATTgg                                                    >  1:3327917/1‑69 (MQ=255)
            tgttgtTAAGCCCCAGCGCCGAAACATACTCCAGCGCCGCGCCAAGGCCAATAATACCTCCGGTATTgg                                                    >  1:2022915/1‑69 (MQ=255)
            tgttgtTAAGCCCCAGCGCCGAAACATACTCCAGCGCCGCGCCAAGGCCAATAATACCTCCGGTATTgg                                                    >  1:1981867/1‑69 (MQ=255)
            tgttgtTAAGCCCCAGCGCCGAAACATACTCCAGCGCCGCGCCAAGGCCAATAATACCTCCGGTATTgg                                                    >  1:2841767/1‑69 (MQ=255)
            tgttgtTAAGCCCCAGCGCCGAAACATACTCCAGCGCCGCGCCAAGGCCAATAATACCTCCGGTATTgg                                                    >  1:708081/1‑69 (MQ=255)
            tgttgtTAAGCCCCAGCGCCGAAACATACTCCAGCGCCGCGCCAAGGCCAATAATACCTCCGGTATTg                                                     >  1:899167/1‑68 (MQ=255)
            tgttgtTAAGCCCCAGCGCCGAAACATACTCCAGCGCCGCGCCAAGGCCAATAATACCTCCGGTATTg                                                     >  1:1989057/1‑68 (MQ=255)
              ttattaAGCCCCAGCGCCGAAACATACTCCAGCGCCGCGCCAAGACCAATGATGCCCCCGGTATTggg                                                   >  1:293727/1‑68 (MQ=255)
              ttattaAGCCCCAGCGCCGAAACATACTCCAGCGCCGCGCCAAGACCAATGATGCCCCCGGTATTggg                                                   >  1:1599803/1‑68 (MQ=255)
              ttattaAGCCCCAGCGCCGAAACATACTCCAGCGCCGCGCCAAGACCAATGATGCCCCCGGTATTgg                                                    <  1:198686/67‑1 (MQ=255)
              ttattaAGCCCCAGCGCCGAAACATACTCCAGCGCCGCGCCAAGACCAATGATGCCCCCGGTATTgg                                                    <  1:2766339/67‑1 (MQ=255)
                attaAGCCCCAGCGCCGAAACATACTCCAGCGCCGCGCCAAGACCAATGATGCCCCCGGTATTGGgtgt                                                <  1:1255715/69‑1 (MQ=255)
                attaAGCCCCAGCGCCGAAACATACTCCAGCGCCGCGCCAAGACCAATGATGCCCCCGGTATTGGgtgt                                                <  1:1911358/69‑1 (MQ=255)
                attaAGCCCCAGCGCCGAAACATACTCCAGCGCCGCGCCAAGACCAATGATGCCCCCGGTATTGGgtg                                                 <  1:292250/68‑1 (MQ=255)
                 ttaAGCCCCAGCGCCGAAACATACTCCAGCGCCGCGCCAAGACCAATGATGCCCCCGGTATTGGGTGTa                                               >  1:3320990/1‑69 (MQ=255)
                 ttaAGCCCCAGCGCCGAAACATACTCCAGCGCCGCGCCAAGACCAATGATGCCCCCGGTATTGGGTGTa                                               >  1:363924/1‑69 (MQ=255)
                 ttaAGCCCCAGCGCCGAAACATACTCCAGCGCCGCGCCAAGACCAATGATGCCCCCGGTATTGGGTGTa                                               >  1:980712/1‑69 (MQ=255)
                 ttaAGCCCCAGCGCCGAAACATACTCCAGCGCCGCGCCAAGACCAATGATGCCCCCGGTATTGGGTGTa                                               >  1:1780484/1‑69 (MQ=255)
                 ttaAGCACCAGCGCCGAAACATACTCCAGCGCCGCGCCAAGACCAATGATGCCCCCGGTATTGGGTGTa                                               >  1:442433/1‑69 (MQ=255)
                          aGCGCCGAAACATACTCCAGCGCCGCGCCAAGGCCAATAATACCTCCGGTATTGGGCGTACCGGCTTCa                                      >  1:412356/1‑69 (MQ=255)
                          aGCGCCGAAACATACTCCAGCGCCGCGCCAAGGCCAATAATACCTCCGGTATTGGGCGTACCGGCTTCa                                      >  1:2101926/1‑69 (MQ=255)
                                       aCTCCAGCGCCGCGCCAAGGCCAATAATACCTCCGGTATTGGGCGTACCGGCTTCAAACCGCCACggcg                         >  1:102034/1‑67 (MQ=255)
                                       aCTCCAGCGCCGCGCCAAGGCCAATAATACCTCCGGTATTGGGCGTACCGGCTTCAAACCGCCACggcg                         >  1:759155/1‑67 (MQ=255)
                                       aCTCCAGCGCCGCGCCAAGGCCAATAATACCTCCGGTATTGGGCGTACCGGCTTCAAACCGCCACggcg                         >  1:3345539/1‑67 (MQ=255)
                                       aCTCCAGCGCCGCGCCAAGGCCAATAATACCTCCGGTATTGGGCGTACCGGCTTCAAACCGCCACggcg                         >  1:2612804/1‑67 (MQ=255)
                                       aCTCCAGCGCCGCGCCAAGACCAATGATGCCCCCGGTATTGGGTGTACCGGCTTCAAACCGCCATGGTg                         <  1:3237951/69‑1 (MQ=255)
                                       aCTCCAGCGCCGCGCCAAGACCAATGATGCCCCCGGTATTGGGTGTACCGGCTTCAAACCGCCATGGTg                         <  1:1621249/69‑1 (MQ=255)
                                       aCTCAAGCGCCGCGCCAAGGCCAATAATACCTCCGGTATTGGGCGTACCGGCTTCAAACCGCCACggcg                         >  1:1223514/1‑67 (MQ=255)
                                        cTCCAGCGCCGCGCCAAGACCAATGATGCCCCCGTTATTGGGTGTACCGGCTTCAAACCGCCATGGTGc                        >  1:2070141/1‑69 (MQ=255)
                                        cTCCAGCGCCGCGCCAAGACCAATGATGCCCCCGGTATTGGGTGTACCGGCTTCAAACCGCCATGGTg                         <  1:2623709/68‑1 (MQ=255)
                                        cTCCAGCGCCGCGCCAAGACCAATGATGCCCCCGGTATTGGGTGTACCGGCTTCAAACCGCCATGGTGc                        >  1:826535/1‑69 (MQ=255)
                                        cTCCAGCGCCGCGCCAAGACCAATGATGCCCCCGGTATTGGTGTACCGGCTTCAAACCGCCATGGTGCt                       >  1:188798/1‑69 (MQ=255)
                                                        aGACCAATGATGCCCCCGGTATTGGGTGTACCGGCTTCAAACCGCCATGGTGCTTTGGTCCAGGTAGt         <  1:2344288/68‑1 (MQ=255)
                                                        aGACCAATGATGCCCCCGGTATTGGGTGTACCGGCTTCAAACCGCCATGGTGCTTTGGTCCAGGTAGt         <  1:522720/68‑1 (MQ=255)
                                                        aGACCAATGATGCCCCCGGTATTGGGTGTACCGGCTTCAAACCGCCATGGTGCTTTGGTCCAGGTAGt         <  1:854177/68‑1 (MQ=255)
                                                            cAATGATGCCCCCGGTATTGGGTGTACCGGCTTCAAACCGCCATGGTGCTTTGGTCCAGGTAGTGCCtt    <  1:1227848/69‑1 (MQ=255)
                                                              atgatgCCCCCGGTATTGGGTGTACCGGCTTCAAACCGCCATGGTGCTTTGGTCCAGGTAGTGCCTTc   >  1:2316863/1‑68 (MQ=255)
                                                              atgatgCCCCCGGTATTGGGTGTACCGGCTTCAAACCGCCATGGTGCTTTGGTCCAGGTAGTGCCTTc   >  1:2527661/1‑68 (MQ=255)
                                                              atgatgCCCCCGGTATTGGGTGTACCGGCTTCAAACCGCCATGGTGCTTTGGTCCAGGTAGTGCCTTCa  >  1:2729932/1‑69 (MQ=255)
                                                              atgatgCCCCCGGTATTGGGTGTACCGGCTTCAAACCGCCATGGTGCTTTGGTCCAGGTAGTGCCTTCa  >  1:1899092/1‑69 (MQ=255)
                                                              atgatgCCCCCGGTATTGGGTGTACCGGCTTCAAACCGCCATGGTGCTTTGGTCCAGGTAGTGCCTTCa  >  1:3278674/1‑69 (MQ=255)
                                                                gatgCCCCCGGTATTGGGTGTACCGGCTTCAAACCGCCATGGTGCTTTGGTCCAGGTAGTGCCTTCa  >  1:743892/1‑67 (MQ=255)
                                                                |                                                                  
CATACTCGGCTATGTTATTAAGCCCCAGCGCCGAAACATACTCCAGCGCCGCGCCAAGACCAATGATGCCCCCGGTATTGGGTGTACCGGCTTCAAACCGCCATGGTGCTTTGGTCCAGGTAGTGCCTTCA  >  minE/1087299‑1087429

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: