Predicted mutation
evidence seq id position mutation freq annotation gene description
RA minE 138,172 A→G 24.7% N665N (AAT→AAC htrE ← predicted outer membrane usher protein

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*minE138,1720AG24.7% 130.0 / 48.1 73N665N (AAT→AAChtrEpredicted outer membrane usher protein
Reads supporting (aligned to +/- strand):  ref base A (26/29);  new base G (2/16);  total (28/45)
Fisher's exact test for biased strand distribution p-value = 5.99e-03
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 7.24e-01

AGTCCACCGCTATGCAATACAAAACCACCGTCGGTGCTGAGAGAAACTTGACGGCTGTTATCGCTATTTGCAGAAATTGAACCTGCCAGCGTTCCCCATGGTGACTCATAGCTGGCATAACCCCCAACA  >  minE/138107‑138235
                                                                 |                                                               
agtcCACCGCTATGCAATACAAAACCACCGTCGGTGCTGAGAGAAACTTGACGGCTGTTATCGCTATTTg                                                             >  1:1322617/1‑70 (MQ=255)
agtcAACCGCTATGCAATACAAAACCACCGTCGGTGCTGAGAGAAACTTGACGGCTGTTATCGCTATTTg                                                             >  1:1024353/1‑70 (MQ=255)
    cACCGCTATGCAATACAAAACCACCGTCGGTGCTGAGAGAAACTTGACGGCTGTTATCGCTGTTTGCAGaa                                                        <  1:360727/71‑1 (MQ=255)
    cACCGCTATGCAATACAAAACCACCGTCGGTGCTGAGAGAAACTTGACGGCTGTTATCGCTGTTTGCAGaa                                                        <  1:1568412/71‑1 (MQ=255)
    cACCGCTATGCAATACAAAACCACCGTCGGTGCTGAGAGAAACTTGACGGCTGTTATCGCTATTTGCAGa                                                         <  1:1004943/70‑1 (MQ=255)
    cACCGCTATGCAATACAAAACCACCGTCGGTGCTGAGAGAAACTTGACGGCTGTTATCGCTATTTGCAGa                                                         <  1:4011949/70‑1 (MQ=255)
    cACCGCTATGCAATACAAAACCACCGTCGGTGCTGAGAGAAACTTGACGGCTGTTATCGCTATTTGCAGa                                                         <  1:3447639/70‑1 (MQ=255)
    cACCGCTATGCAATACAAAACCACCGTCGGTGCTGAGAGAAACTTGACGGCTGTTATCGCTATTTGCAGa                                                         <  1:2252189/70‑1 (MQ=255)
         cTATGCAATACAAAACCACCGTCGGTGCTGAGAGAAACTTGACGGCTGTTATCGCTGTTTGCAGAAATTGa                                                   <  1:909559/71‑1 (MQ=255)
         cTATGCAATACAAAACCACCGTCGGTGCTGAGAGAAACTTGACGGCTGTTATCGCTGTTTGCAGAAATTGa                                                   <  1:1949889/71‑1 (MQ=255)
         cTATGCAATACAAAACCACCGTCGGTGCTGAGAGAAACTTGACGGCTGTTATCGCTGTTTGCAGAAATTGa                                                   <  1:1021381/71‑1 (MQ=255)
         cTATGCAATACAAAACCACCGTCGGTGCTGAGAGAAACTTGACGGCTGTTATCGCTATTTGCAGAAATTGa                                                   <  1:3988662/71‑1 (MQ=255)
         cTATGCAATACAAAACCACCGTCGGTGCTGAGAGAAACTTGACGGCTGTTATCGCTATTTGCAGAAATTGa                                                   <  1:1666579/71‑1 (MQ=255)
         cTATGCAATACAAAACCACCGTCGGTGCTGAGAGAAACTTGACGGCTGTTATCGCTATTTGCAGAAATTGa                                                   <  1:2888979/71‑1 (MQ=255)
           aTGCAATACAAAACCACCGTCGGTGCTGAGAGAAACTTGGCGGCTGTTATCGCTGTTTGCAGAAATTg                                                    <  1:3359115/68‑1 (MQ=255)
           aTGCAATACAAAACCACCGTCGGTGCTGAGAGAAACTTGACGGCTGTTATCGCTGTTTGCAGAAATTg                                                    <  1:3731516/68‑1 (MQ=255)
           aTGCAATACAAAACCACCGTCGGTGCTGAGAGAAACTTGACGGCTGTTATCGCTGTTTGCAGAAATTg                                                    <  1:2230944/68‑1 (MQ=255)
           aTGCAATACAAAACCACCGTCGGTGCTGAGAGAAACTTGACGGCTGTTATCGCTGTTTGCAGAAATTg                                                    <  1:2512131/68‑1 (MQ=255)
            tGCAATACAAAACCACCGTCGGTGCTGAGAGAAACTTGACGGCTGTTATCGCTATTTGCAGAAATTGa                                                   <  1:1462069/68‑1 (MQ=255)
            tGCAATACAAAACCACCGTCGGTGCTGAGAGAAACTTGACGGCTGTTATCGCTATTTGCAGAAATTGa                                                   <  1:1269570/68‑1 (MQ=255)
            tGCAATACAAAACCACCGTCGGTGCTGAGAGAAACTTGACGGCTGTTATCGCTATTTGCAGAAATTGAAcc                                                <  1:36881/71‑1 (MQ=255)
             gCAATACAAAACCACCGTCGGTGCTGAGAGAAACTTGACGGCTGTTATCGCTATTTGCAGAAATTGAAc                                                 >  1:480915/1‑69 (MQ=255)
              cAATACAAAACCACCGTCGGTGCTGAGAGAAACTTGACGGCTGTTATCGCTATTTGCAGAAATTGAAc                                                 >  1:2224246/1‑68 (MQ=255)
              cAATACAAAACCACCGTCGGTGCTGAGAGAAACTTGACGGCTGTTATCGCTATTTGCAGAAATTGAAc                                                 >  1:298784/1‑68 (MQ=255)
                aTACAAAACCACCGTCGGTGCTGAGAGAAACTTGACGGCTGTTATCGCTAt                                                                >  1:2995693/1‑51 (MQ=255)
                aTACAAAACCACCGTCGGTGCTGAGAGAAACTTGACGGCTGTTATCGCTAt                                                                >  1:3056014/1‑51 (MQ=255)
                aTACAAAACCACCGTCGGTGCTGAGAGAAACTTGACGGCTGTTATCGCTAt                                                                >  1:593930/1‑51 (MQ=255)
                  aCAAAACCACCGTCGGTGCTGAGAGAAACTTGACGGCTGTTATCGCTATTTGCAGaaa                                                       >  1:3105584/1‑58 (MQ=255)
                   cAAAACCACCGTCGGTGCTGAGAGAAACTTGACGGCTGTTATCGCTATTTGCAGAAATTGAACCTGCCAGc                                         >  1:1383551/1‑71 (MQ=255)
                     aaaCCACCGTCGGTGCTGAGAGAAACTTGACGGCTGTTATCGCTATTTGCAGAAATTGAACCTGCCAGCGt                                       >  1:2946375/1‑71 (MQ=255)
                     aaaCCACCGTCGGTGCTGAGAGAAACTTGACGGCTGTTATCGCTATTTGCAGAAATTGAACCTGCCAGCGt                                       >  1:2444408/1‑71 (MQ=255)
                     aaaCCACCGTCGGTGCTGAGAGAAACTTGACGGCTGTTATCGCTATTTGCAGAAATTGAACCTGCCAGCGt                                       >  1:2481397/1‑71 (MQ=255)
                     aaaCCACCGTCGGTGCTGAGAGAAACTTGACGGCTGTTATCGCTATTTGCAGAAATTGAACCTGCCAGCGt                                       >  1:1089079/1‑71 (MQ=255)
                     aaaCCACCGTCGGTGCTGAGAGAAACTTGACGGCTGTTATCGCTATTTGCAGAAATTGAACCTGCCAGCGt                                       >  1:2449067/1‑71 (MQ=255)
                     aaaCCACCGTCGGTGCTGAGAGAAACTTGACGGCTGTTATCGCTATTTGCAGAAATTGAACCTGCCAGCGt                                       >  1:440664/1‑71 (MQ=255)
                       accaccGTCGGTGCTGAGAGAAACTTGACGGCTGTTATCGCTATTTGCAGAAATTGAACCTGCCAGCGTTc                                     <  1:401499/71‑1 (MQ=255)
                       accaccGTCGGTGCTGAGAGAAACTTGACGGCTGTTATCGCTATTTGCAGAAATTGAACCTGCCAGCGTTc                                     <  1:4110424/71‑1 (MQ=255)
                       accaccGGCGGTGCTGAGAGAAACTTGACGGCTGTTATCGCTATTTGCAGAAATTGAACCTGCCAGCGTTc                                     <  1:3221799/71‑1 (MQ=255)
                         caccGTCGGTGCTGAGAGAAACTTGACGGCTGTTATCGCTATTTg                                                             <  1:3751406/45‑1 (MQ=255)
                          accGTCGGTGCTGAGAGAAACTTGACGGCTGTTATCGCTATTTGCAGAAATTGa                                                   <  1:2422110/54‑1 (MQ=255)
                           ccGTCGGTGCTGAGAGAAACTTGACGGCTGTTATCGCTATTTGCAGAAATTGAACCTGCCAGCGTACCCCa                                 >  1:259752/1‑71 (MQ=255)
                            cGTCGGTGCTGAGAGAAACTTGACGGCTGTTATCGCTGTTTGCAGAAATTGAACCTGCCAGCGTTCCCCa                                 <  1:1517451/70‑1 (MQ=255)
                            cGTCGGTGCTGAGAGAAACTTGACGGCTGTTATCGCTGTTTGCAGAAATTGAACCTGCCAGCGTTCCCCAt                                <  1:2984975/71‑1 (MQ=255)
                                ggTGCTGAGAGAAACTTGACGGCTGTTATCGCTATTTGCAGAAATTGAACCTGCCAGCGTTCCCCATGGt                             <  1:497823/70‑1 (MQ=255)
                                    cTGAGAGAAACTTGACGGCTGTTATCGCTGTTTGCAGAAATTGAACCTGCCAGCGTTCCCCATGGTGACt                         <  1:3005864/70‑1 (MQ=255)
                                    cTGAGAGAAACTTGACGGCTGTTATCGCTGTTTGCAGAAATTGAACCTGCCAGCGTTCCCCATGGTGACt                         <  1:826059/70‑1 (MQ=255)
                                    cTGAGAGAAACTTGACGGCTGTTATCGCTGTTTGCAGAAATTGAACCTGCCAGCGTTCCCCATGGTGACt                         <  1:126325/70‑1 (MQ=255)
                                      gagagaAACTTGACGGCTGTTATCGCTATTTGCAGAAATTGAAc                                                 >  1:1860855/1‑44 (MQ=255)
                                      gagagaAACTTGACGGCTGTTATCGCTATTTGCAGAAATTGAAc                                                 >  1:2024778/1‑44 (MQ=255)
                                      gagagaAACTTGACGGCTGTTATCGCTATTTGCAGAAATTGAAc                                                 >  1:2839821/1‑44 (MQ=255)
                                      gagagaAACTTGACGGCTGTTATCGCTATTTGCAGAAATTGAAc                                                 >  1:3359931/1‑44 (MQ=255)
                                      gagagaAACTTGACGGCTGTTATCGCTATTTGCAGAAATTGAAc                                                 >  1:3568417/1‑44 (MQ=255)
                                       agagaAACTTGACGGCTGTTATCGCTGTTTGCAGAAATTGAACCTGCCAGCGTTCCCCATGGTGa                           >  1:1424877/1‑65 (MQ=255)
                                        gagaAACTTGACGGCTGTTATCGCTGTTTGCAGAAATTGAACCTGCCAGCGTTCCCCATGGTGACTCATAg                    <  1:1471033/71‑1 (MQ=255)
                                        gagaAACTTGACGGCTGTTATCGCTATTTGCAGAAATTGAACCTGCCCGCGTTCCCCATGGTGACTCAt                      <  1:1792978/69‑1 (MQ=255)
                                         agaAACTTGACGGCTGTTATCGCTATTTGCAGAAATTGAACCTGCCAGCGTTCCCCATGGTGACTCATAGc                   <  1:3382725/71‑1 (MQ=255)
                                          gaAACTTGACGGCTGTTATCGCTATTTGCAGAAATTGAACCTGCCAGCGTTCCCCATGGTGACTCATAGc                   <  1:3735753/70‑1 (MQ=255)
                                           aaaCTTGACGGCTGTTATCGCTATTTGCAGAAATTGAACCTGCCAGCGTTCCCCATGGTGACTCATAGCTg                 >  1:1004109/1‑71 (MQ=255)
                                           aaaCTTGACGGCTGTTATCGCTATTTGCAGAAATTGAACCTGCCAGCGTTCCCCATGGTGACTCATAGCTg                 >  1:2390268/1‑71 (MQ=255)
                                           aaaCTTGACGGCTGTTATCGCTATTTGCAGAAATTGAACCTGCCAGCGTTCCCCATGGTGACTCATAGCTg                 >  1:1489118/1‑71 (MQ=255)
                                            aaCTTGACGGCTGTTATCGCTATTTGCAGAAATTg                                                    <  1:1700463/35‑1 (MQ=255)
                                            aaCTTGACGGCTGTTATCGCTATTTGCAGAAATTGAACCTGCCAGCGTTCCCCATGGTGACTCATAGCTgg                <  1:1809256/71‑1 (MQ=255)
                                            aaCTTGACGGCTGTTATCGCTATTTGCAGAAATTGAACCTGCCAGCGTTCCCCATGGTGACTCATAGCTgg                <  1:1461584/71‑1 (MQ=255)
                                                 gACGGCTGTTATCGCTATTTGCAGAAATTGAACCTGCCAGCGTTCCCCATGGTGACTCATAGCTGGCAt             <  1:1062092/69‑1 (MQ=255)
                                                 gACGGCTGTTATCGCTATTTGCAGAAATTGAACCTGCCAGCGTTCCCCATGGTGACTCATAGCTGGCAt             <  1:3371099/69‑1 (MQ=255)
                                                 gACGGCTGTTATCGCTATTTGCAGAAATTGAACCTGCCAGCGTTCCCCATGGTGACTCATAGCTGGCAt             <  1:2517512/69‑1 (MQ=255)
                                                  aCGGCTGTTATCGCTATTTGCAGAAATTGAACCTGCCAGCGTTCCCCATGGTGACTCATAGCTGGCATaa           <  1:3340995/70‑1 (MQ=255)
                                                  aCGGCTGTTATCGCTATTTGCAGAAATTGAACCTGCCAGCGTTCCCCATGGTGACTCATAGCTGGCATaa           <  1:3807190/70‑1 (MQ=255)
                                                  aCGGCTGTTATCGCTATTTGCAGAAATTGAACCTGCCAGCGTTCCCCATGGTGACTCATAGCTGGCATaa           <  1:2105360/70‑1 (MQ=255)
                                                      cTGTTATCGCTGTTTGCAGAAATTGAACCTGCCAGCGTTCCCCATGGTGACTCATAGCTGGCATAAccccc      <  1:1461360/71‑1 (MQ=255)
                                                          tATCGCTGTTTGCAGAAATTGAACCTGCCAGCGTTCCCCATGGTGACTCATAGCTGGCATAACCCCCAACa  >  1:1541597/1‑71 (MQ=255)
                                                          tATCGCTATTTGCAGAAATTGAACCTTCCAGCGTTCCCCATGGTGACTCATa                     <  1:2384224/52‑1 (MQ=255)
                                                           aTCGCTATTTGCAGAAATTGAACCTGCCAGCGTTCCCCATGGTGACTCATAGCTGGCATAACCCCCAACa  >  1:195522/1‑70 (MQ=255)
                                                                 |                                                               
AGTCCACCGCTATGCAATACAAAACCACCGTCGGTGCTGAGAGAAACTTGACGGCTGTTATCGCTATTTGCAGAAATTGAACCTGCCAGCGTTCCCCATGGTGACTCATAGCTGGCATAACCCCCAACA  >  minE/138107‑138235

Alignment Legend
Aligned base mismatch/match (shaded by quality score): ATCG/ATCG < 3 ≤ ATCG/ATCG < 27 ≤ ATCG/ATCG < 32 ≤ ATCG/ATCG < 36 ≤ ATCG/ATCG
Unaligned base: atcg    Masked matching base: atcg    Alignment gap:     Deleted base: