Predicted mutation
evidence seq id position mutation freq annotation gene description
RA minE 72,264 A→G 20.1% V56V (GTA→GTG ilvH → acetolactate synthase III, thiamin‑dependent, small subunit

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*minE72,2640AG20.1% 120.1 / 21.0 84V56V (GTA→GTGilvHacetolactate synthase III, thiamin‑dependent, small subunit
Reads supporting (aligned to +/- strand):  ref base A (17/50);  new base G (11/6);  total (28/56)
Fisher's exact test for biased strand distribution p-value = 3.63e-03
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 8.05e-01

CGTTGCGCCAACCGACGATCCGACATTATCGCGTATGACCATCCAGACCGTGGGCGATGAAAAAGTACTTGAGCAGATCGAAAAGCAATTACACAAACTGGTCGATGTCTTGCGCGTGAGTG  >  minE/72198‑72319
                                                                  |                                                       
cGTTGCGCCAACCGACGATCCGACATTATCGCGTATGACCATCCAGACCGTGGGCGATGAAAAAGTa                                                         <  1:1626109/67‑1 (MQ=255)
cGTTGCGCCAACCGACGATCCGACATTATCGCGTATGACCATCCAGACCGTGGGCGATGAAAAAGTa                                                         <  1:1534840/67‑1 (MQ=255)
  ttGCGCCAACCGACGATCCGACATTATCGGGTATGACCATCCAGACCGTGGGCGATGAAAAAGTa                                                         <  1:2361896/65‑1 (MQ=255)
       ccAACCGACGATCCGACATTATCGCGTATGACCATCCAGACCGTGGGCGATGAAAAAGTACTTGAGCAg                                                <  1:2371562/69‑1 (MQ=255)
       ccAACCGACGATCCGACATTATCGCGTATGACCATCCAGACCGTGGGCGATGAAAAAGTACTTGAGCAg                                                <  1:2955671/69‑1 (MQ=255)
       ccAACCGACGATCCGACATTATCGCGTATGACCATCCAGACCGTGGGCGATGAAAAAGTACTTGAGCAg                                                <  1:1383938/69‑1 (MQ=255)
                   ccGACATTATCGCGTATGACCATCCAGACCGTGGGCGATGAAAAAGTACTTGAGCAGATCGAAAAGCa                                     >  1:528761/1‑68 (MQ=255)
                          tATCGCGTATGACCATCCAGACCGTGGGCGATGAAAAAGTACTTGAGCAGATCGAAAAGCAATTAcaca                             <  1:3306126/69‑1 (MQ=255)
                          tATCGCGTATGACCATCCAGACCGTGGGCGATGAAAAAGTACTTGAGCAGATCGAAAAGCAATTAcaca                             <  1:232615/69‑1 (MQ=255)
                          tATCGCGTATGACCATCCAGACCGTGGGCGATGAAAAAGTACTTGAGCAGATCGAAAAGCAATTAcaca                             <  1:3046537/69‑1 (MQ=255)
                          tATCGCGTATGACCATCCAGACCGTGGGCGATGAAAAAGTACTTGAGCAGATCGAAAAGCAATTAcaca                             <  1:3058856/69‑1 (MQ=255)
                          tATCGCGTATGACCATCCAGACCGTGGGCGATGAAAAAGTACTTGAGCAGATCGAAAAGCAATTAcaca                             <  1:1652068/69‑1 (MQ=255)
                          tATCGCGTATGACCATCCAGACCGTGGGCGATGAAAAAGTACTTGAGCAGATCGAAAAGCAATTAcaca                             <  1:868278/69‑1 (MQ=255)
                          tATCGCGTATGACCATCCAGACCGTGGGCGATGAAAAAGTACTTGAGCAGATCGAAAAGCAATTAcaca                             <  1:1033627/69‑1 (MQ=255)
                          tATCGCGTATGACCATCCAGACCGTGGGCGATGAAAAAGTACTTGAGCAGATCGAAAAGCAATTAcaca                             <  1:792328/69‑1 (MQ=255)
                          tATCGCGTATGACCATCCAGACCGTGGGCGATGAAAAAGTACTTGAGCAGATCGAAAAGCAAGTAcaca                             <  1:436352/69‑1 (MQ=255)
                           aTCGCGTATGACCATCCAGACCGTGGGCGATGAAAAAGTACTTGAGCAGATCGAAAAGCa                                     <  1:2727572/60‑1 (MQ=255)
                           aTCGCGTATGACCATCCAGACCGTGGGCGATGAAAAAGTACTTGAGCAGATCGAAAAGCa                                     <  1:1132875/60‑1 (MQ=255)
                           aTCGCGTATGACCATCCAGACCGTGGGCGATGAAAAAGTACTTGAGCAGATCGAAAAGCa                                     <  1:2558096/60‑1 (MQ=255)
                           aTCGCGTATGACCATCCAGACCGTGGGCGATGAAAAAGTACTTGAGCAGATCGAAAAGCa                                     <  1:2073134/60‑1 (MQ=255)
                           aTCGCGTATGACCATCCAGACCGTGGGCGATGAAAAAGTACTTGAGCAGATCGAAAAGCa                                     <  1:2049829/60‑1 (MQ=255)
                           aTCGCGTATGACCATCCAGACCGTGGGCGATGAAAAAGTACTTGAGCAGATCGAAAAGCa                                     <  1:1478889/60‑1 (MQ=255)
                           aTCGCGTATGACCATCCAGACCGTGGGCGATGAAAAAGTACTTGAGCAGATCGAAAAGCa                                     <  1:1976093/60‑1 (MQ=255)
                           aTCGCGTATGACCATCCAGACCGTGGGCGATGAAAAAGTACTTGAGCAGATCGAAAAGCa                                     <  1:3375350/60‑1 (MQ=255)
                           aTCGCGTATGACCATCCAGACCGTGGGCGATGAAAAAGTACTTGAGCAGATCGAAAAGCa                                     <  1:710663/60‑1 (MQ=255)
                           aTCGCGTATGACCATCCAGACCGTGGGCGATGAAAAAGTACTTGAGCAGATCGAAAAGCa                                     <  1:696290/60‑1 (MQ=255)
                           aTCGCGTATGACCATCCAGACCGTGGGCGATGAAAAAGTACTTGAGCAGATCGAAAAGCAATTAcaca                             <  1:1147301/68‑1 (MQ=255)
                           aTCGCGTATGACCATCCAGACCGTGGGCGATGAAAAAGTACTTGAGCAGATCGAAAAGCAATTACACaa                            >  1:3000425/1‑69 (MQ=255)
                             cgcgTATGACCATTCAGACCGTGGGCGATGAAAAAGTGCTTGAGCAAATCGAAAAGCAATTACACaagc                          >  1:617625/1‑67 (MQ=255)
                             cgcgTATGACCATTCAGACCGTGGGCGATGAAAAAGTGCTTGAGCAAATCGAAAAGCAATTACACaagc                          >  1:2701975/1‑67 (MQ=255)
                             cgcgTATGACCATTCAGACCGTGGGCGATGAAAAAGTGCTTGAGCAAATCGAAAAGCAATTACACaagc                          >  1:3035506/1‑67 (MQ=255)
                             cgcgTATGACCATTCAGACCGTGGGCGATGAAAAAGTGCTTGAGCAAATCGAAAAGCAATTACACaagc                          >  1:2521235/1‑67 (MQ=255)
                             cgcgTATGACCATTCAGACCGTGGGCGATGAAAAAGTGCTTGAGCAAATCGAAAAGCAATTACACaagc                          >  1:2458665/1‑67 (MQ=255)
                             cgcgTATGACCATTCAGACCGTGGGCGATGAAAAAGTGCTTGAGCAAATCGAAAAGCAATTACACaagc                          >  1:917746/1‑67 (MQ=255)
                             cgcgTATGACCATTCAGACCGTGGGCGATGAAAAAGTGCTTGAGCAAATCGAAAAGCAATTACACaagc                          >  1:2390363/1‑67 (MQ=255)
                             cgcgTATGACCATTCAGACCGTGGGCGATGAAAAAGTGCTTGAGCAAATCGAAAAGCAATTACACaagc                          >  1:1322115/1‑67 (MQ=255)
                             cgcgTATGACCATTCAGACCGTGGGCGATGAAAAAGTGCTTGAGCAAATCGAAAAGCAATTACACaagc                          >  1:1690102/1‑67 (MQ=255)
                             cgcgTATGACCATTCAGACCGTGGGCGATGAAAAAGTGCTTGAGCAAATCGAAAAGCAATTACACaagc                          >  1:1673635/1‑67 (MQ=255)
                             cgcgTATGACCATTCAGACCGTGGGCGATGAAAAAGTGCTTGAGCAAATCGAAAAGCAATTACACaagc                          >  1:1672489/1‑67 (MQ=255)
                                           cAGACCGTGGGCGATGTAAAAGTACTTGAGCAGATCGAAAAGc                                      <  1:1256888/43‑1 (MQ=255)
                                           cAGACCGTGGGCGATGAAAAAGTACTTGAGCAGATCGAAAAg                                       >  1:2721919/1‑42 (MQ=255)
                                           cAGACCGTGGGCGATGAAAAAGTACTTGAGCAGATCGAAAAGc                                      <  1:2932216/43‑1 (MQ=255)
                                           cAGACCGTGGGCGATGAAAAAGTACTTGAGCAGATCGAAAAGc                                      >  1:101736/1‑43 (MQ=255)
                                           cAGACCGTGGGCGATGAAAAAGTACTTGAGCAGATCGAAAAGc                                      >  1:1864476/1‑43 (MQ=255)
                                           cAGACCGTGGGCGATGAAAAAGTACTTGAGCAGATCGAAAAGc                                      <  1:3226127/43‑1 (MQ=255)
                                           cAGACCGTGGGCGATGAAAAAGTACTTGAGCAGATCGAAAAGc                                      >  1:324354/1‑43 (MQ=255)
                                           cAGACCGTGGGCGATGAAAAAGTACTTGAGCAGATCGAAAAGc                                      >  1:3286275/1‑43 (MQ=255)
                                           cAGACCGTGGGCGATGAAAAAGTACTTGAGCAGATCGAAAAGc                                      >  1:3386681/1‑43 (MQ=255)
                                           cAGACCGTGGGCGATGAAAAAGTACTTGAGCAGATCGAAAAGc                                      >  1:3411051/1‑43 (MQ=255)
                                           cAGACCGTGGGCGATGAAAAAGTACTTGAGCAGATCGAAAAGc                                      <  1:365078/43‑1 (MQ=255)
                                           cAGACCGTGGGCGATGAAAAAGTACTTGAGCAGATCGAAAAGc                                      <  1:1547050/43‑1 (MQ=255)
                                           cAGACCGTGGGCGATGAAAAAGTACTTGAGCAGATCGAAAAGc                                      >  1:531409/1‑43 (MQ=255)
                                           cAGACCGTGGGCGATGAAAAAGTACTTGAGCAGATCGAAAAGc                                      <  1:577490/43‑1 (MQ=255)
                                           cAGACCGTGGGCGATGAAAAAGTACTTGAGCAGATCGAAAAGc                                      <  1:2857629/43‑1 (MQ=255)
                                           cAGACCGTGGGCGATGAAAAAGTACTTGAGCAGATCGAAAAGc                                      >  1:2450598/1‑43 (MQ=255)
                                           cAGACCGTGGGCGATGAAAAAGTACTTGAGCAGATCGAAAAGc                                      >  1:2795841/1‑43 (MQ=255)
                                           cAGACCGTGGGCGATGAAAAAGTACTTGAGCAGATCGAAAAGc                                      >  1:2192624/1‑43 (MQ=255)
                                           cAGACCGTGGGCGATGAAAAAGTACTTGAGCAGATCGAAAAGc                                      <  1:2525912/43‑1 (MQ=255)
                                           cAGACCGTGGGCGATGAAAAAGTACTTGAGCAGATCGAAAAGc                                      >  1:2016129/1‑43 (MQ=255)
                                           cAGACCGTGGGCGATGAAAAAGTACTTGAGCAGATCGAAAAGCAAt                                   <  1:290609/46‑1 (MQ=255)
                                           cAGACCGTGGGCGATGAAAAAGTACTTGAGCAGATCGAAAAGCAAt                                   <  1:1217100/46‑1 (MQ=255)
                                           cAGACCGTGGGCGATGAAAAAGTACTTGAGCAGATCGAAAAGCAAt                                   <  1:823788/46‑1 (MQ=255)
                                           cAGACCGTGGGCGATGAAAAAGTACTTGAGCAGATCGAAAAGCAAt                                   <  1:780732/46‑1 (MQ=255)
                                           cAGACCGTGGGCGATGAAAAAGTACTTGAGCAGATCGAAAAGCAAt                                   <  1:522526/46‑1 (MQ=255)
                                           cAGACCGTGGGCGATGAAAAAGTACTTGAGCAGATCGAAAAGCAAt                                   <  1:2185648/46‑1 (MQ=255)
                                           cAGACCGTGGGCGATGAAAAAGTACTTGAGCAGATCGAAAAGCAAt                                   <  1:2627929/46‑1 (MQ=255)
                                           cAGACCGTGGGCGATGAAAAAGTACTTGAGCAGATCGAAAAGCAAt                                   <  1:1590033/46‑1 (MQ=255)
                                           cAGACCGTGGGCGATGAAAAAGTACTTGAGCAGATCGAAAAGCAAt                                   <  1:2685743/46‑1 (MQ=255)
                                           cAGACCGTGGGCGATGAAAAAGTACTTGAGCAGATCGAAAAGCAAt                                   <  1:2709957/46‑1 (MQ=255)
                                           cAGACCGTGGGCGATGAAAAAGTACTTGAGCAGATCGAAAAGCAAt                                   <  1:1996868/46‑1 (MQ=255)
                                           cAGACCGTGGGCGATGAAAAAGTACTTGAGCAGATCGAAAAGCAAt                                   <  1:1918805/46‑1 (MQ=255)
                                           cAGACCGTGGGCGATGAAAAAGTACTTGAGCAGATCGAAAAGCAAt                                   <  1:305974/46‑1 (MQ=255)
                                           cAGACCGTGGGCGATGAAAAAGTACTTGAGCAGATCGAAAAGCAAt                                   <  1:2750901/46‑1 (MQ=255)
                                           cAGACCGTGGGCGATGAAAAAGTACTTGAGCAGATCGAAAAGCAATTAcaca                             >  1:1221136/1‑52 (MQ=255)
                                           cAGACCGTGGGCGATGAAAAAGTACTTCAGCAGATCGAAAAGc                                      <  1:68660/43‑1 (MQ=255)
                                           cAGAACGTGGGCGATGAAAAAGTACTTGAGCAGATCGAAAAGCAAt                                   <  1:1344154/46‑1 (MQ=255)
                                               ccGTGGGCGATGAAAAAGTGCTTGAGCAAATCGAAAAGCAATTa                                 <  1:2821057/44‑1 (MQ=38)
                                               ccGTGGGCGATGAAAAAGTGCTTGAGCAAATCGAAAAGCAATTa                                 <  1:2461551/44‑1 (MQ=38)
                                               ccGTGGGCGATGAAAAAGTGCTTGAGCAAATCGAAAAGCAATTa                                 <  1:1489222/44‑1 (MQ=38)
                                               ccGTGGGCGATGAAAAAGTGCTTGAGCAAATCGAAAAGCAATTa                                 <  1:2737230/44‑1 (MQ=38)
                                               ccGTGGGCGATGAAAAAGTGCTTGAGCAAATCGAAAAGCAATTa                                 <  1:380561/44‑1 (MQ=38)
                                               ccGTGGGCGATGAAAAAGTGCTTGAGCAAATCGAAAAGCAATTa                                 <  1:2032478/44‑1 (MQ=38)
                                                    ggCGATGAAAAAGTACTTGAGCAGCTCGAAAAGCAATTACACAAACTGGTCGATGTCTTGCGCgtga     >  1:2567327/1‑67 (MQ=255)
                                                    ggCGATGAAAAAGTACTTGAGCAGATCGAAAAGCa                                     <  1:3139828/35‑1 (MQ=255)
                                                    ggCGATGAAAAAGTACTTGAGCAGATCGAAAAGCAATTACACAAACTGGTCGATGTCTTGCGCgtga     >  1:297367/1‑67 (MQ=255)
                                                      cGATGAAAAAGTACTTGAGCAGATCGAAAAGCAATTACACAAACTGGTCGATGTCTTGCGCgtgagtg  <  1:2525366/68‑1 (MQ=255)
                                                      cGATGAAAAAGTACTTGAGCAGATCGAAAAGCAATTACACAAACTGGTCGATGTCTTGCGCgtgagtg  <  1:2311116/68‑1 (MQ=255)
                                                                  |                                                       
CGTTGCGCCAACCGACGATCCGACATTATCGCGTATGACCATCCAGACCGTGGGCGATGAAAAAGTACTTGAGCAGATCGAAAAGCAATTACACAAACTGGTCGATGTCTTGCGCGTGAGTG  >  minE/72198‑72319

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: