Predicted mutation
evidence seq id position mutation freq annotation gene description
RA minE 72,273 G→A 20.8% Q59Q (CAG→CAA 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,2730GA20.8% 117.8 / 20.2 81Q59Q (CAG→CAAilvHacetolactate synthase III, thiamin‑dependent, small subunit
Reads supporting (aligned to +/- strand):  ref base G (17/47);  new base A (11/6);  total (28/53)
Fisher's exact test for biased strand distribution p-value = 8.10e-03
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.43e-01

CCAACCGACGATCCGACATTATCGCGTATGACCATCCAGACCGTGGGCGATGAAAAAGTACTTGAGCAGATCGAAAAGCAATTACACAAACTGGTCGATGTCTTGCGCGTGAGTG  >  minE/72205‑72319
                                                                    |                                              
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:1033627/69‑1 (MQ=255)
                   tATCGCGTATGACCATCCAGACCGTGGGCGATGAAAAAGTACTTGAGCAGATCGAAAAGCAATTAcaca                             <  1:792328/69‑1 (MQ=255)
                   tATCGCGTATGACCATCCAGACCGTGGGCGATGAAAAAGTACTTGAGCAGATCGAAAAGCAATTAcaca                             <  1:868278/69‑1 (MQ=255)
                   tATCGCGTATGACCATCCAGACCGTGGGCGATGAAAAAGTACTTGAGCAGATCGAAAAGCAAGTAcaca                             <  1:436352/69‑1 (MQ=255)
                    aTCGCGTATGACCATCCAGACCGTGGGCGATGAAAAAGTACTTGAGCAGATCGAAAAGCa                                     <  1:1976093/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:2727572/60‑1 (MQ=255)
                    aTCGCGTATGACCATCCAGACCGTGGGCGATGAAAAAGTACTTGAGCAGATCGAAAAGCa                                     <  1:710663/60‑1 (MQ=255)
                    aTCGCGTATGACCATCCAGACCGTGGGCGATGAAAAAGTACTTGAGCAGATCGAAAAGCa                                     <  1:3375350/60‑1 (MQ=255)
                    aTCGCGTATGACCATCCAGACCGTGGGCGATGAAAAAGTACTTGAGCAGATCGAAAAGCa                                     <  1:696290/60‑1 (MQ=255)
                    aTCGCGTATGACCATCCAGACCGTGGGCGATGAAAAAGTACTTGAGCAGATCGAAAAGCa                                     <  1:1478889/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:3035506/1‑67 (MQ=255)
                      cgcgTATGACCATTCAGACCGTGGGCGATGAAAAAGTGCTTGAGCAAATCGAAAAGCAATTACACaagc                          >  1:2701975/1‑67 (MQ=255)
                      cgcgTATGACCATTCAGACCGTGGGCGATGAAAAAGTGCTTGAGCAAATCGAAAAGCAATTACACaagc                          >  1:917746/1‑67 (MQ=255)
                      cgcgTATGACCATTCAGACCGTGGGCGATGAAAAAGTGCTTGAGCAAATCGAAAAGCAATTACACaagc                          >  1:2521235/1‑67 (MQ=255)
                      cgcgTATGACCATTCAGACCGTGGGCGATGAAAAAGTGCTTGAGCAAATCGAAAAGCAATTACACaagc                          >  1:2458665/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)
                                                                    |                                              
CCAACCGACGATCCGACATTATCGCGTATGACCATCCAGACCGTGGGCGATGAAAAAGTACTTGAGCAGATCGAAAAGCAATTACACAAACTGGTCGATGTCTTGCGCGTGAGTG  >  minE/72205‑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: