Predicted mutation
evidence seq id position mutation freq annotation gene description
RA minE 616,470 G→C 16.6% T320T (ACG→ACC aroA → 5‑enolpyruvylshikimate‑3‑phosphate synthetase

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*minE616,4700GC16.6% 136.3 / 18.2 66T320T (ACG→ACCaroA5‑enolpyruvylshikimate‑3‑phosphate synthetase
Reads supporting (aligned to +/- strand):  ref base G (31/24);  new base C (6/5);  total (37/29)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.08e-01

GGTGAACTGAACGCTATTGATATGGATATGAACCATATTCCTGATGCGGCGATGACCATTGCCACGGCGGCGTTATTTGCAAAAGGCACCACCACGCTGCGCAATATCTATAACTGGCGTGTTA  >  minE/616405‑616528
                                                                 |                                                          
ggtgAACTGAACGCTATTGATATGGATATGAACCATATTCCTGATGCGGCGATGACCATTGCCAcggcg                                                         <  1:121094/69‑1 (MQ=255)
ggtgAACTGAACGCTATTGATATGGATATGAACCATATTCCTGATGCGGCGATGACCATTGCCAcggcg                                                         <  1:2955133/69‑1 (MQ=255)
ggtgAACTGAACGCTATTGATATGGATATGAACCATATTCCTGATGCGGCGATGACCATTGCCAcggcg                                                         <  1:2879556/69‑1 (MQ=255)
                  gatatggatatgAACCATATTCCTGATGCGGCGATGACCATTGCCAcggc                                                          >  1:1705314/1‑50 (MQ=255)
                  gatatggatatgAACCATATTCCTGATGCGGCGATGACCATTGCCACGGCGGCGTTATTTGCAAAAGGc                                       >  1:2464912/1‑69 (MQ=255)
                  gatatggatatgAACCATATTCCTGATGCGGCGATGACCATTGCCACGGCGGCGTTATTTGCAAAAGGc                                       >  1:1248966/1‑69 (MQ=255)
                  gatatggatatgAACCATATTCCTGATGCGGCGATGACCATTGCCACGGCGGCGTTATTTGCAAAAGGc                                       >  1:3170025/1‑69 (MQ=255)
                  gatatggatatgAACCATATTCCTGATGCGGCGATGACCATTGCCACGGCGGCGTTATTTGCAAAAGGc                                       >  1:2551734/1‑69 (MQ=255)
                  gatatggatatgAACCATATTCCTGATGCGGCGATGACCATTGCCACGGCGGCGTTATTTGCAAAAGGc                                       >  1:2044445/1‑69 (MQ=255)
                  gatatggatatgAACCATATTCCTGATGCGGCGATGACCATTGCCACGGCGGCGTTATTTGCAAAAGGc                                       >  1:2820561/1‑69 (MQ=255)
                     atggatatgAACCATATTCCTGATGCGGCGATGACCATTGCCACGGCGGCGTTATTTGCAAAAGGcac                                     <  1:3211729/68‑1 (MQ=255)
                     atggatatgAACCATATTCCTGATGCGGCGATGACCATTGCCACGGCGGCGTTATTTGCAAAAGGcac                                     <  1:2015613/68‑1 (MQ=255)
                      tggatatgAACCATATTCCTGATGCGGCGATGACCATTGCCACGGCGGCGTTATTTGCAAAAGGcacc                                    >  1:2176677/1‑68 (MQ=255)
                           atgAACCATATTCCTGATGCGGCGATGACCATTGCCACGGCGGCGTTATTTGCAAAAGGCACCACCACg                              >  1:2421223/1‑69 (MQ=255)
                           atgAACCATATTCCTGATGCGGCGATGACCATTGCCACGGCGGCGTTATTTGCAAAAGGCACCACCACg                              >  1:2252865/1‑69 (MQ=255)
                           atgAACCATATTCCTGATGCGGCGATGACCATTGCCACGGCGGCGTTATTTGCAAAAGGCACCACCACg                              >  1:2083545/1‑69 (MQ=255)
                           atgAACCATATTCCTGATGCGGCGATGACCATTGCCACGGCGGCGTTATTTGCAAAAGGCACCACCACg                              >  1:3056668/1‑69 (MQ=255)
                           atgAACCATATTCCTGATGCGGCGATGACCATTGCCACGGCGGCGTTATTTGCAAAAGGCACCACCACg                              >  1:3146337/1‑69 (MQ=255)
                           atgAACCATATTCCTGATGCGGCGATGACCATTGCCACGGCGGCGTTATTTGCAAAAGGCACCACCACg                              >  1:3250854/1‑69 (MQ=255)
                                 cATATTCCTGATGCGGCGATGACCATTGCCACGGCGGCGTTATTTGCAAAAGGTACCACCACGCTgcgc                        <  1:2821153/69‑1 (MQ=255)
                                 cATATTCCTGATGCGGCGATGACCATTGCCACGGCGGCGTTATTTGCAAAAGGCACCACCACGCTgcgc                        <  1:2732468/69‑1 (MQ=255)
                                 cATATTCCTGATGCGGCGATGACCATTGCCACGGCGGCGTTATTTGCAAAAGGCACCACCACGCTgcgc                        <  1:2545286/69‑1 (MQ=255)
                                 cATATTCCTGATGCGGCGATGACCATTGCCACGGCGGCGTTATTTGCAAAAGGCACCACCACGCTgcgc                        <  1:463665/69‑1 (MQ=255)
                                 cATATTCCTGATGCGGCGATGACCATTGCCACGGCGGCGTTATTTGCAAAAGGCACCACCACGCTgcgc                        <  1:1107130/69‑1 (MQ=255)
                                 cATATTCCTGATGCGGCGATGACCATTGCCACGGCGGCGTTATTTGCAAAAGGCACCACCACGCTgcgc                        <  1:944106/69‑1 (MQ=255)
                                 cATATTCCTGATGCGGCGATGACCATTGCCACGGCGGCGTTATTTGCAAAAGGCACCACCACGCTgcgc                        <  1:1378985/69‑1 (MQ=255)
                                 cATATTCCTGATGCGGCGATGACCATTGCCACGGCGGCGTTATTTGCAAAAGGCACCACCACGCTgcgc                        <  1:1963344/69‑1 (MQ=255)
                                     ttCCTGATGCGGCGATGACCATTGCCACGGCGGCGTTATTTGCAAAAGGCACCACCACGCTGCGCAata                    >  1:2342102/1‑69 (MQ=255)
                                     ttCCTGATGCGGCGATGACCATTGCCACGGCGGCGTTATTTGCAAAAGGCACCACCACGCTGCGCAata                    >  1:367342/1‑69 (MQ=255)
                                     ttCCTGATGCGGCGATGACCATTGCCACGGCGGCGTTATTTGCAAAAGGCACCACCACGCTGCGCAata                    >  1:3063197/1‑69 (MQ=255)
                                     ttCCTGATGCGGCGATGACCATTGCCACGGCGGCGTTATTTGCAAAAGGCACCACCACGCTGCGCAata                    >  1:2117229/1‑69 (MQ=255)
                                       ccTGATGCGGCGATGACCATTGCCACGGCGGCGTTATTTGCAAAAGGcaccac                                  <  1:12313/53‑1 (MQ=255)
                                       ccTGATGAGGCGATGACCATTGCCACGGCGGCGTTATTTGCAAAAGGcaccac                                  <  1:2815959/53‑1 (MQ=255)
                                         tGATGCGGCGATGACCATTGCCACGGCGGCGTTATTTGCAAAAGGCACCACCACGCTGCGCAATATCt                 >  1:1950562/1‑68 (MQ=255)
                                          gATGCGGCGATGACCATTGCCACCGCGGCGTTATTTGCAAAAGGCACCACCACGCTGCGCAATATCtat               <  1:3062553/69‑1 (MQ=255)
                                          gATGCGGCGATGACCATTGCCACCGCGGCGTTATTTGCAAAAGGCACCACCACGCTGCGCAATATCtat               <  1:916225/69‑1 (MQ=255)
                                          gATGCGGCGATGACCATTGCCACCGCGGCGTTATTTGCAAAAGGCACCACCACGCTGCGCAATATCtat               <  1:528475/69‑1 (MQ=255)
                                          gATGCGGCGATGACCATTGCCACCGCGGCGTTATTTGCAAAAGGCACCACCACGCTGCGCAATATCtat               <  1:1410616/69‑1 (MQ=255)
                                          gATGCGGCGATGACCATTGCCACCGCGGCGTTATTTGCAAAAGGCACCACCACGCTGCGCAATATCtat               <  1:1748455/69‑1 (MQ=255)
                                             gcggcgATGACCATTGCCACGGCGGCGTTATTTGCAAAAGGCACCACCACGCTGCGCAATATCtat               <  1:2818598/66‑1 (MQ=255)
                                             gcggcgATGACCATTGCCACGGCGGCGTTATTTGCAAAAGGCACCACCACGCTGCGCAATATCtat               <  1:2827039/66‑1 (MQ=255)
                                             gcggcgATGACCATTGCCACGGCGGCGTTATTTGCAAAAGGCACCACCACGCTGCGCAATATCtat               <  1:2272168/66‑1 (MQ=255)
                                             gcggcgATGACCATTGCCACGGCGGCGTTATTTGCAAAAGGCACCACCACGCTGCGCAATATCtat               <  1:723310/66‑1 (MQ=255)
                                             gcggcgATGACCATTGCCACGGCGGCGTTATTTGCAAAAGGCACCACCACGCTGCGCAATATCtat               <  1:1747727/66‑1 (MQ=255)
                                             gcggcgATGACCATTGCCACGGCGGCGTTATTTGCAAAAGGCACCACCACGCTGCGCAATATCtat               <  1:2422458/66‑1 (MQ=255)
                                             gcggcgATGACCATTGCCACGGCGGCGTTATTTGCAAAAGGCACCACCACGCTGCGCAATATCtat               <  1:3231897/66‑1 (MQ=255)
                                             gcggcgATGACCATTGCCACCGCGGCGTTATTTGCAAAAgg                                        >  1:3400283/1‑41 (MQ=255)
                                             gcggcgATGACCATTGCCACCGCGGCGTTATTTGCAAAAGGCACCACCACGCTGCGCAATATCTATAAc            >  1:642147/1‑69 (MQ=255)
                                             gcggcgATGACCATTGCCACCGCGGCGTTATTTGCAAAAGGCACCACCACGCTGCGCAATATCTATAAc            >  1:471756/1‑69 (MQ=255)
                                             gcggcgATGACCATTGCCACCGCGGCGTTATTTGCAAAAGGCACCACCACGCTGCGCAATATCTATAAc            >  1:2503921/1‑69 (MQ=255)
                                             gcggcgATGACCATTGCCACCGCGGCGTTATTTGCAAAAGGCACCACCACGCTGCGCAATATCTATAAc            >  1:314432/1‑69 (MQ=255)
                                             gcggcgATGACCATTGCCACCGCGGCGTTATTTGCAAAAGGCACCACCACGCTGCGCAATATCTATAAc            >  1:3149538/1‑69 (MQ=255)
                                              cggcgATGACCATTGCCACGGCGGCGTTATTTGCAAAAGGCACCACCACGCTGCGCAATATCTATAACt           >  1:2009193/1‑69 (MQ=255)
                                              cggcgATGACCATTGCCACGGCGGCGTTATTTGCAAAAGGCACCACCACGCTGCGCAATATCTATAACt           >  1:869037/1‑69 (MQ=255)
                                              cggcgATGACCATTGCCACGGCGGCGTTATTTGCAAAAGGCACCACCACGCTGCGCAATATCTATAACt           >  1:2474235/1‑69 (MQ=255)
                                              cggcgATGACCATTGCCACGGCGGCGTTATTTGCAAAAGGCACCACCACGCTGCGCAATATCTATAACt           >  1:2474952/1‑69 (MQ=255)
                                              cggcgATGACCATTGCCACGGCGGCGTTATTTGCAAAAGGCACCACCACGCTGCGCAATATCTATAACt           >  1:289628/1‑69 (MQ=255)
                                              cggcgATGACCATTGCCACGGCGGCGTTATTTGCAAAAGGCACCACCACGCTGCGCAATATCTATAACt           >  1:1845961/1‑69 (MQ=255)
                                                      aCCATTGCCACGGGGGCGTTATTTGCAAAAGGCACCACCACGCTGCGCAat                     >  1:3399239/1‑51 (MQ=255)
                                                      aCCATTGCCACGGCGGCGTTATTTGCAAAAGGcaccaccac                               >  1:888024/1‑41 (MQ=255)
                                                      aCCATTGCCACGGCGGCGTTATTTGCAAAAGGCACCACCACGCTGCGCAat                     >  1:1350443/1‑51 (MQ=255)
                                                      aCCATTGCCACGGCGGCGTTATTTGCAAAAGGCACCACCACGCTGCGCAat                     >  1:1640915/1‑51 (MQ=255)
                                                      aCCATTGCCACGGCGGCGTTATTTGCAAAAGGCACCACCACGCTGCGCAat                     >  1:226325/1‑51 (MQ=255)
                                                      aCCATTGCCACGGCGGCGTTATTTGCAAAAGGCACCACCACGCTGCGCAat                     >  1:3062912/1‑51 (MQ=255)
                                                      aCCATTGCCACGGCGGCGTTATTTGCAAAAGGCACCACAACGCTGCGCAat                     >  1:2972260/1‑51 (MQ=255)
                                                        cATTGCCACGGCGGCGTTATTTGCAAAAGGCACCACCACGCTGCGCAATATCTATAACTGGCGTGTTa  <  1:244656/68‑1 (MQ=255)
                                                        cATTGCCACGGCGGCGTTATTTGCAAAAGGCACCACCACGCTGCGCAATATCTATAACTGGCGTGTTa  <  1:2323208/68‑1 (MQ=255)
                                                        cATTGCCACGGCGGCGTTATTTGCAAAAGGCACCACCACGCTGCGCAATATCTATAACTGGCGTGTTa  <  1:2145564/68‑1 (MQ=255)
                                                        cATTGCCACGGCGGCGTTATTTGCAAAAGGCACCACCACGCTGCGCAATATCTATAACTGGCGTGTTa  <  1:780973/68‑1 (MQ=255)
                                                        cATTGCCACGGCGGCGTTATTTGCAAAAGGCACCACCACGCTGCGCAATATCTATAACTGGCGTGTTa  <  1:1195610/68‑1 (MQ=255)
                                                                 |                                                          
GGTGAACTGAACGCTATTGATATGGATATGAACCATATTCCTGATGCGGCGATGACCATTGCCACGGCGGCGTTATTTGCAAAAGGCACCACCACGCTGCGCAATATCTATAACTGGCGTGTTA  >  minE/616405‑616528

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: