Predicted mutation
evidence seq id position mutation freq annotation gene description
RA minE 369,416 G→A 37.8% Q121Q (CAG→CAA entF → enterobactin synthase multienzyme complex component, ATP‑dependent

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*minE369,4160GA37.8% 42.9 / 46.4 74Q121Q (CAG→CAAentFenterobactin synthase multienzyme complex component, ATP‑dependent
Reads supporting (aligned to +/- strand):  ref base G (26/20);  new base A (13/15);  total (39/35)
Fisher's exact test for biased strand distribution p-value = 4.75e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.27e-01

GCATTAATGCAGGCGGATTTGCAACAAGATCTGCGCGTCGATAGCGGTAAACCACTGGTCTTTCATCAGCTGATACAGGTGGCGGATAACCGCTGGTACTGGTATCAGCGTTATCACCATTTGCTGGTCGAT  >  minE/369348‑369479
                                                                    |                                                               
gCATTAATGCAGGCGGATTTGCAACAAGATCTGCGCGTCGATAGCGGTAAACCGCTGGTCTTTCATCAg                                                                 >  1:382781/1‑69 (MQ=255)
gCATTAATGCAGGCGGATTTGCAACAAGATCTGCGCGTCGATAGCGGTAAACCACTGGTCTTTCATCAg                                                                 >  1:901834/1‑69 (MQ=255)
gCATTAATGCAGGCGGATTTGCAACAAGATCTGCGCGTCGATAGCGGTAAACCACTGGTCTTTCATCAg                                                                 >  1:2012883/1‑69 (MQ=255)
gCATTAATGCAGGCGGATTTGCAACAAGATCTGCGCGTCGATAGCGGTAAACCACTGGTCTTTCATCAg                                                                 >  1:2438089/1‑69 (MQ=255)
gCATTAATGCAGGCGGATTTGCAACAAGATCTGCGCGTCGATAGCGGTAAACCACTGGTCTTTCATCAg                                                                 >  1:1623350/1‑69 (MQ=255)
gCATTAATGCAGGCGGATTTGCAACAAGATCTGCGCGTCGATAGCGGTAAACCACTGGTCTTTCATCAg                                                                 >  1:161141/1‑69 (MQ=255)
gCATTAATGCAGGCGGATTTGCAACAAGATCTGCGCGTCGATAGCGGTAAACCACTGGTCTTTCATCAg                                                                 >  1:1591357/1‑69 (MQ=255)
gCATTAATGCAGGCGGATTTGCAACAAGATCTGCGCGTCGATAGCGGTAAACCACTGGTCTTTCATCAg                                                                 >  1:2497243/1‑69 (MQ=255)
gCATTAATGCAGGCGGATTTGCAACAAGATCTGCGCGTCGATAGCGGTAAACCACTGGTCTTTCATCAg                                                                 >  1:3362790/1‑69 (MQ=255)
gCATTAATGCAGGCGGATTTGCAACAAGATCTGCGCGTCGATAGCGGTAAACCACTGGTCTTTCATCAg                                                                 >  1:3365084/1‑69 (MQ=255)
gCATTAATGCAGGCGGATTTGCAACAAGATCTGCGCGTCGATAGCGGTAAACCACTGGTCTTTCATCAg                                                                 >  1:1210699/1‑69 (MQ=255)
gCATTAATGCAGGCGGATTTGCAACAAGATCTGCGCGTCGATAGCGGTAAACCACTGGTCTTTCATCAg                                                                 >  1:738089/1‑69 (MQ=255)
     aaTGCAGGCGGATTTGCAACAAGATCTGCGCGTCGATAGCGGTAAACCACTGGTCTTTCATCAGCTGa                                                             >  1:2962222/1‑68 (MQ=255)
     aaTGCAGGCGGATTTGCAACAAGATCTGCGCGTCGATAGCGGTAAACCACTGGTCTTTCATCAGCTGa                                                             >  1:2958253/1‑68 (MQ=255)
     aaTGCAGGCGGATTTGCAACAAGATCTGCGCGTCGATAGCGGTAAACCACTGGTCTTTCATCAGCTGa                                                             >  1:2507590/1‑68 (MQ=255)
     aaTGCAGGCGGATTTGCAACAAGATCTGCGCGTCGATAGCGGTAAACCACTGGTCTTTCATCAGCTGa                                                             >  1:849648/1‑68 (MQ=255)
     aaTGCAGGCGGATTTGCAACAAGATCTGCGCGTCGATAGCGGTAAACCACTGGTCTTTCATCAGCTGa                                                             >  1:229619/1‑68 (MQ=255)
     aaTGCAGGCGGATTTGCAACAAGATCTGCGCGTCGATAGCGGTAAACCACTGGTCTTTCATCAGCTGa                                                             >  1:947335/1‑68 (MQ=255)
     aaTGCAGGCGGATTTGCAACAAGATCTGCGCGTCGATAGCGGTAAACCACTGGTCTTTCATCAGCTGa                                                             >  1:1104460/1‑68 (MQ=255)
     aaTGCAGGCGGATTTGCAACAAGATCTGCGCGTCGATAGCGGTAAACCACTGGTCTTTCATCAGCTGa                                                             >  1:1736529/1‑68 (MQ=255)
     aaTGCAGGCGGATATGCAACAAGATCTGCGCGTCGATAGCGGTAAACCACTGGTCTTTCATCAGCTGa                                                             >  1:3364926/1‑68 (MQ=255)
                 tttGCAACAAGATCTGCGCGTCGTTAGCGGTAAACCACTGGTCTTTCATCAGCTGa                                                             <  1:2850255/56‑1 (MQ=255)
                 tttGCAACAAGATCTGCGCGTCGATAGCGGTAAACCACTGGTCTTTCATCAGCTGa                                                             >  1:483857/1‑56 (MQ=255)
                 tttGCAACAAGATCTGCGCGTCGATAGCGGTAAACCACTGGTCTTTCATCAGCTGa                                                             <  1:1655380/56‑1 (MQ=255)
                 tttGCAACAAGATCTGCGCGTCGATAGCGGTAAACCACTGGTCTTTCATCAGCTGa                                                             <  1:3340986/56‑1 (MQ=255)
                 tttGCAACAAGATCTGCGCGTCGATAGCGGTAAACCACTGGTCTTTCATCAGCTGa                                                             >  1:2089763/1‑56 (MQ=255)
                 tttGCAACAAGATCTGCGCGTCGATAGCGGTAAACCACTGGTCTTTCATCAGCTGa                                                             <  1:2087258/56‑1 (MQ=255)
                 tttGCAACAAGATCTGCGCGTCGATAGCGGTAAACCACTGGTCTTTCATCAGCTGa                                                             >  1:2673826/1‑56 (MQ=255)
                 tttGCAACAAGATCTGCGCGTCGATAGCGGTAAACCACTGGTCTTTCATCAGCTGa                                                             >  1:3201522/1‑56 (MQ=255)
                 tttGCAACAAGATCTGCGCGTCGATAGCGGTAAACCACTGGTCTTTCATCAGCTGa                                                             <  1:2734997/56‑1 (MQ=255)
                 tttGCAACAAGATCTGCGCGTCGATAGCGGTAAACCACTGGTCTTTCATCAGCTGa                                                             >  1:2854032/1‑56 (MQ=255)
                 tttGCAACAAGATCTGCGCGTCGATAGCGGTAAACCACTGGTCTTTCATCAGCTGATACAGGTGGCgg                                                 >  1:1191908/1‑68 (MQ=255)
                 tttGCAACAAGATCTGCGCGTCGATAGCGGTAAACCACTGGTCTTTCATCAGCTGATACAGGTGGCgg                                                 >  1:1341545/1‑68 (MQ=255)
                 tttGCAACAAGATCTGCGCGTCGATAGCGGTAAACCACTGGTCTTTCATCAGCTGATACAGGTGGCgg                                                 >  1:2145187/1‑68 (MQ=255)
                   tGCAACAAGATCTGCGCGTCGATAGCGGTAAACCACTGGTCTTTCATCAGCTGATACAGGTGGCGGATa                                              <  1:460495/69‑1 (MQ=255)
                   tGCAACAAGATCTGCGCGTCGATAGCGGTAAACCACTGGTCTTTCATCAGCTGATACAGGTGGCGGATa                                              <  1:3051163/69‑1 (MQ=255)
                                   cgTCGATAGCGGTAAACCACTGGTCTTTCATCAGCTGATACAGGTGGCGGATAACCGctggtactggt                               <  1:334314/68‑1 (MQ=255)
                                   cgTCGATAGCGGTAAACCACTGGTCTTTCATCAGCTGATACAGGTGGCGGATAACCGctggtactggt                               <  1:417331/68‑1 (MQ=255)
                                   cgTCGATAGCGGTAAACCACTGGTCTTTCATCAGCTGATACAGGTGGCGGATAACCGctggtactggt                               <  1:1269772/68‑1 (MQ=255)
                                   cgTCGATAGCGGTAAACCACTGGTCTTTCATCAGCTGATACAGGTGGCGGATAACCGctggtactgg                                <  1:1191115/67‑1 (MQ=255)
                                   cgTCGATAGCGGTAAACCACTGGTCTTTCATCAGCTGATACAGGTGGCGGATAACCGctggtactgg                                <  1:2034790/67‑1 (MQ=255)
                                   cgTCGATAGCGGTAAACCACTGGTCTTTCATCAGCTGATACAGGTGGCGGATAACCGctggtactgg                                <  1:1928213/67‑1 (MQ=255)
                                   cgTCGATAGCGGTAAACCACTGGTCTTTCATCAGCTGATACAGGTGGCGGATAACCGctggtactgg                                <  1:1924318/67‑1 (MQ=255)
                                   cgTCGATAGCGGTAAACCACTGGTCTTTCATCAGCTGATACAGGTGGCGGATAACCGctggtactgg                                <  1:1913596/67‑1 (MQ=255)
                                           gCGGTAAACCACTGGTCTTTCATCAGCTGATACAGGTGGCGGATAACCGCTGGTACTGGTATCAGCGtt                      <  1:1747202/69‑1 (MQ=255)
                                           gCGGTAAACCACTGGTCTTTCATCAGCTGATACAGGTGGCGGATAACCGCTGGTACTGGTATCAGCGtt                      <  1:1029497/69‑1 (MQ=255)
                                           gCGGTAAACCACTGGTCTTTCATCAGCTGATACAGGTGGCGGATAACCGCTGGTACTGGTATCAGCGtt                      <  1:2445941/69‑1 (MQ=255)
                                           gCGGTAAACCACTGGTCTTTCATCAGCTGATACAGGTGGCGGATAACCGCTGGTACTGGTATCAGCGtt                      <  1:2024736/69‑1 (MQ=255)
                                           gCGGTAAACCACTGGTCTTTCATCAGCTGATACAGGTGGCGGATAACCGCTGGTACTGGTATCAGCGtt                      <  1:3001636/69‑1 (MQ=255)
                                                           gtttCATCAACTGATACAGGTGGCGGATAACCGCTGGTACTGGTATCAGCGTTATCACCATTTGCTgg       >  1:2343808/2‑68 (MQ=255)
                                                           cTTTCATCAACTGATACAGGTGGCGGATAACCGCTGGTACTGGTATCAGCGTTATCACCATTTGCTgg       >  1:3373929/1‑68 (MQ=255)
                                                           cTTTCATCAACTGATACAGGTGGCGGATAACCGCTGGTACTGGTATCAGCGTTATCACCATTTGCTgg       >  1:3179941/1‑68 (MQ=255)
                                                           cTTTCATCAACTGATACAGGTGGCGGATAACCGCTGGTACTGGTATCAGCGTTATCACCATTTGCTgg       >  1:1001258/1‑68 (MQ=255)
                                                           cTTTCATCAACTGATACAGGTGGCGGATAACCGCTGGTACTGGTATCAGCGTTATCACCATTTGCTgg       >  1:3234607/1‑68 (MQ=255)
                                                           cTTTCATCAACTGATACAGGTGGCGGATAACCGCTGGTACTGGTATCAGCGTTATCACCATTTGCTgg       >  1:2690454/1‑68 (MQ=255)
                                                           cTTTCATCAACTGATACAGGTGGCGGATAACCGCTGGTACTGGTATCAGCGTTATCACCATTTGCTgg       >  1:2645588/1‑68 (MQ=255)
                                                           cTTTCATCAACTGATACAGGTGGCGGATAACCGCTGGTACTGGTATCAGCGTTATCACCATTTGCTgg       >  1:263148/1‑68 (MQ=255)
                                                           cTTTCATCAACTGATACAGGTGGCGGATAACCGCTGGTACTGGTATCAGCGTTATCACCATTTGCTgg       >  1:581284/1‑68 (MQ=255)
                                                           cTTTCATCAACTGATACAGGTGGCGGATAACCGCTGGTACTGGTATCAGCGTTATCACCATTTGCTgg       >  1:2100053/1‑68 (MQ=255)
                                                           cTTTCATCAACTGATACAGGTGGCGGATAACCGCTGGTACTGGTATCAGCGTTATCACCATTTGCTgg       >  1:2045288/1‑68 (MQ=255)
                                                           cTTTCATCAACTGATACAGGTGGCGGATAACCGCTGGTACTGGTATCAGCGTTATCACCATTTGCTg        >  1:2676331/1‑67 (MQ=255)
                                                           cTTTCATCAACTGATACAGGTGGCGGATAACCGCTGGTACTGGTATCAGCGTTATCACCATTTGATgg       >  1:2161721/1‑68 (MQ=255)
                                                                atcaACTGATACAGGTGGCGGATAACCGCTGGTACTTGTATCAGCGTTATCACCATTTGCTGGTCg    <  1:3337708/66‑1 (MQ=255)
                                                                atcaACTGATACAGGTGGCGGATAACCGCTGGTACTGGTATCAGCGTTATCACCATTTGCTGGTCg    <  1:352587/66‑1 (MQ=255)
                                                                atcaACTGATACAGGTGGCGGATAACCGCTGGTACTGGTATCAGCGTTATCACCATTTGCTGGTCg    <  1:1286595/66‑1 (MQ=255)
                                                                atcaACTGATACAGGTGGCGGATAACCGCTGGTACTGGTATCAGCGTTATCACCATTTGCTGGTCg    <  1:849012/66‑1 (MQ=255)
                                                                atcaACTGATACAGGTGGCGGATAACCGCTGGTACTGGTATCAGCGTTATCACCATTTGCTGGTCg    <  1:285255/66‑1 (MQ=255)
                                                                atcaACTGATACAGGTGGCGGATAACCGCTGGTACTGGTATCAGCGTTATCACCATTTGCTGGTCg    <  1:2495388/66‑1 (MQ=255)
                                                                atcaACTGATACAGGTGGCGGATAACCGCTGGTACTGGTATCAGCGTTATCACCATTTGCTGGTCg    <  1:1761865/66‑1 (MQ=255)
                                                                atcaACTGATACAGGTGGCGGATAACCGCTGGTACTGGTATCAGCGTTATCACCATTTGCTGGTCg    <  1:1680723/66‑1 (MQ=255)
                                                                atcaACTGATACAGGTGGCGGATAACCGCTGGTACTGGTATCAGCGTTATCACCATTTGCTGGTCg    <  1:2444743/66‑1 (MQ=255)
                                                                atcaACTGATACAGGTGGCGGATAACCGCTGGTACTGGTATCAGCGTTATCACCATTTGCTGGTCg    <  1:1443683/66‑1 (MQ=255)
                                                                atcaACTGATACAGGTGGCGGATAACCGCTGGTACTGGTATCAGCGTTATCACCATTTGCTGGTCg    <  1:2880209/66‑1 (MQ=255)
                                                                atcaACTGATACAGGTGGCGGATAACCGCTGGTACTGGTATCAGCGTTATCACCATTTGCTGGTCg    <  1:1314606/66‑1 (MQ=255)
                                                                atcaACTGATACAGGTGGCGGATAACCGCTGGTACTGGTATCAGCGTTATCACCATTTGCTGGTCg    <  1:3019481/66‑1 (MQ=255)
                                                                atcaACTGATACAGGTGGCGGATAACCGCTGGTACTGGTATCAGCGTTATCACCATTTGCTGGTCg    <  1:3153563/66‑1 (MQ=255)
                                                                 tcaACTGATACAGGTGGCGGATAACCGCTGGTACTGGTATCAGCGTTATCACCATTTGCTGGTCg    <  1:3008182/65‑1 (MQ=255)
                                                                   aGCTGATACAGGTGGCGGATAACCGCTGGTACTGGTATCAGCGTTATCACCATTTGCTGGTCGAt  >  1:1382125/1‑65 (MQ=255)
                                                                   aGCTGATACAGGTGGCGGATAACCGCTGGTACTGGTATCAGCGTTATCACCATTTGCTGGTCGAt  >  1:3299055/1‑65 (MQ=255)
                                                                   aGCTGATACAGGTGGCGGATAACCGCTGGTACTGGTATCAGCGTTATCACCATTTGCTGGTCGAt  >  1:2414272/1‑65 (MQ=255)
                                                                   aGCTGATACAGGTGGCGGATAACCGCTGGTACTGGTATCAGCGTTATCACCATTTGCTGGTCGAt  >  1:3093681/1‑65 (MQ=255)
                                                                   aGCTGATACAGGTGGCGGATAACCGCTGGTACTGGTATCAGCGTTATCACCATTTGCTGGTCGAt  >  1:706412/1‑65 (MQ=255)
                                                                   aGCTGATACAGGTGGCGGATAACCGCTGGTACTGGTATCAGCGTTATCACCATTTGCTGGTCGAt  >  1:2273251/1‑65 (MQ=255)
                                                                   aGCTGATACAGGTGGCGGATAACCGCTGGTACTGGTATCAGCGTTATCACCATTTGCTGGTCGAt  >  1:1918825/1‑65 (MQ=255)
                                                                   aGCTGATACAGGTGGCGGATAACCGCTGGTACTGGTATCAGCGTTATCACCATTTGCTGGTCGAt  >  1:1899749/1‑65 (MQ=255)
                                                                   aGCTGATACAGGTGGCGGATAACCGCTGGTACTGGTATCAGCGTTATCACCATTTGCTGGTCGAt  >  1:1802921/1‑65 (MQ=255)
                                                                    |                                                               
GCATTAATGCAGGCGGATTTGCAACAAGATCTGCGCGTCGATAGCGGTAAACCACTGGTCTTTCATCAGCTGATACAGGTGGCGGATAACCGCTGGTACTGGTATCAGCGTTATCACCATTTGCTGGTCGAT  >  minE/369348‑369479

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: