Predicted mutation
evidence seq id position mutation freq annotation gene description
RA minE 856,734 G→C 32.6% T212T (ACC→ACG trpB ← tryptophan synthase, beta subunit

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*minE856,7340GC32.6% 104.9 / 72.5 98T212T (ACC→ACGtrpBtryptophan synthase, beta subunit
Reads supporting (aligned to +/- strand):  ref base G (22/44);  new base C (5/27);  total (27/71)
Fisher's exact test for biased strand distribution p-value = 9.14e-02
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 4.54e-01

GCCGCCAACACAGGCGATAACGGCATCCGGCAGGCGACCTTCTCTTTCCAGAATCTGCGCTTTGGTTTCTTCGCCAATCATCCGCTGAAACTCACGCACAATGGTCGGATAAGGATGCGGGCCAGC  >  minE/856671‑856796
                                                               |                                                              
gccgccAACACAGGCGATAACGGCATCCGGCAGGCGACCTTCTCTTTCCAGAATCTGCGCTTTGGTttc                                                           <  1:2305201/69‑1 (MQ=255)
gccgccAACACAGGCGATAACGGCATCCGGCAGGCGACCTTCTCTTTCCAGAATCTGCGCTTTGGTttc                                                           <  1:1269287/69‑1 (MQ=255)
  cgccAACACAGGCGATAACGGCATCCGGCAGGCGACCTTCTCTTTCCAGAATCTGCGCTTTGGTttc                                                           <  1:3121241/67‑1 (MQ=255)
              cGATAACGGCATCCGGCAGGCGACCTTCTCTTTCCAGAATCTGCGCTTTGGTTTCTTCGCCAatcatc                                              >  1:1657291/1‑68 (MQ=255)
              cGATAACGGCATCCGGCAGGCGACCTTCTCTTTCCAGAATCTGCGCTTTGGTTTCTTCGCCAATCATcc                                             >  1:1110794/1‑69 (MQ=255)
              cGATAACGGCATCCGGCAGGCGACCTTCTCTTTCCAGAATCTGCGCTTTGGTTTCTTCGCCAATCATcc                                             >  1:2903473/1‑69 (MQ=255)
              cGATAACGGCATCCGGCAGGCGACCTTCTCTTTCCAGAATCTGCGCTTTGGTTTCTTCGCCAATCATcc                                             >  1:2417268/1‑69 (MQ=255)
              cGATAACGGCATCCGGCAGGCGACCTTCTCTTTCCAGAATCTGCGCTTTGGTTTCTTCGCCAATCATcc                                             >  1:573694/1‑69 (MQ=255)
              cGATAACGGCATCCGGCAGGCGACCTTCTCTTTCCAGAATCTGCGCTTTGGTTTCTTCGCCAATCATcc                                             >  1:1418747/1‑69 (MQ=255)
              cGATAACGGCATCCGGCAGGCGACCTTCTCTTTCCAGAATCTGCGCTTTGGTTTCTTCGCCAATCATcc                                             >  1:1299299/1‑69 (MQ=255)
                   aCGGCATCCGGCAGGCGACCTTCTCTTTCCAGAATCTGCGCTTTGGTTTCTTCGCCAATCATCCGCTGa                                        <  1:1180255/69‑1 (MQ=255)
                   aCGGCATCCGGCAGGCGACCTTCTCTTTCCAGAATCTGCGCTTTGGTTTCTTCGCCAATCATCCGCTGa                                        <  1:2351667/69‑1 (MQ=255)
                   aCGGCATCCGGCAGGCGACCTTCTCTTTCCAGAATCTGCGCTTTGGTTTCTTCGCCAATCATCCGCTGa                                        <  1:2582494/69‑1 (MQ=255)
                   aCGGCATCCGGCAGGCGACCTTCTCTTTCCAGAATCTGCGCTTTGGTTTCTTCGCCAATCATCCGCTGa                                        <  1:2742906/69‑1 (MQ=255)
                   aCGGCATCCGGCAGGCGACCTTCTCTTTCCAGAATCTGCGCTTTGGTTTCTTCGCCAATCATCCGCTGa                                        <  1:3188511/69‑1 (MQ=255)
                   aCGGCATCCGGCAGGCGACCTTCTCTTTCCAGAATCTGCGCTTTGGTTTCTTCGCCAATCATCCGCTGa                                        <  1:774855/69‑1 (MQ=255)
                   aCGGCATCCGGCAGGCGACCTTCTCTTTCCAGAATCTGCGCTTTGGTTTCTTCGCCAATCATCCGCTGa                                        <  1:1293873/69‑1 (MQ=255)
                      gCATCCGGCAGGCGACCTTCTCTTTCCAGAATCTGCGCTTTGGTTTCTTCGCCAATCATCCGCTGAAAc                                     <  1:3047319/69‑1 (MQ=255)
                          ccGGCAGGCGACCTTCTCTTTCCAGAATCTGCGCTTTCGt                                                              >  1:3158515/1‑40 (MQ=255)
                          ccGGCAGGCGACCTTCTCTTTCCAGAATCTGCGCTTTCGTTTCTTCGCCAATCATCCGCTGAAACTCa                                  >  1:3053034/1‑68 (MQ=255)
                          ccGGCAGGCGACCTTCTCTTTCCAGAATCTGCGCTTTCGTTTCTTCGCCAATCATCCGCTGAAACTCa                                  >  1:3105703/1‑68 (MQ=255)
                          ccGGCAGGCGACCTTCTCTTTCCAGAATCTGCGCTTTCGTTTCTTCGCCAATCATCCGCTGAAACTCa                                  >  1:2419289/1‑68 (MQ=255)
                          ccGGCAGGCGACCTTCTATTTCCAGAATCTGCGCTTTCGTttct                                                          >  1:1161683/1‑44 (MQ=38)
                           cGGCAGGCGACCTTCTCTTTCCAGAATCTGCGCTTTGGTTTCTTCGCCAATCATCCGCTGAAACTCACg                                <  1:837926/69‑1 (MQ=255)
                               aGGCGACCTTCTCTTTCCAGAATCTGCGCTTTGGTTTCTTCGCCAATCATCCGCTGAAACTCACGcaca                            <  1:1609575/69‑1 (MQ=255)
                               aGGCGACCTTCTCTTTCCAGAATCTGCGCTTTGGTTTCTTCGCCAATCATCCGCTGAAACTCACGcaca                            <  1:2376479/69‑1 (MQ=255)
                               aGGCGACCTTCTCTTTCCAGAATCTGCGCTTTGGTTTCTTCGCCAATCATCCGCTGAAACTCACGcaca                            <  1:396190/69‑1 (MQ=255)
                               aGGCGACCTTCTCTTTCCAGAATCTGCGCTTTGGTTTCTTCGCCAATCATCCGCTGAAACTCACGcaca                            <  1:1917522/69‑1 (MQ=255)
                               aGGCGACCTTCTCTTTCCAGAATCTGCGCTTTGGTTTCTTCGCCAATCATCCGCTGAAACTCACGcaca                            <  1:3219491/69‑1 (MQ=255)
                               aGGCGACCTTCTCTTTCCAGAATCTGCGCTTTGGTTTCTTCGCCAATCATCCGCTGAAACTCACGcaca                            <  1:1468422/69‑1 (MQ=255)
                                  cGACCTTCTCTTTCCAGAATCTGCGCTTTGGTTTCTTCGCCAATCATCCGCTGAAACt                                    <  1:3262423/58‑1 (MQ=255)
                                  cGACCTTCTCTTTCCAGAATCTGCGCTTTGGTTTCTTCGCCAATCATCCGCTGAAACt                                    <  1:872425/58‑1 (MQ=255)
                                      cTTCTCTTTCCAGAATCTGCGCTTTGGTTTCTTCGCCAATCATCCGCTGAAACTCACGCACAATGGTCg                     <  1:314400/69‑1 (MQ=255)
                                      cTTCTCTTTCCAGAATCTGCGCTTTGGTTTCTTCGCCAATCATCCGCTGAAACTCACGCACAATGGTCg                     >  1:2869807/1‑69 (MQ=255)
                                         ctctTTCCAGAATCTGCGCTTTGGTTTCTTCGCCAATCATCCGCTGAAACTCACGCACAATGGTCgg                    >  1:785355/1‑67 (MQ=255)
                                         ctctTTCCAGAATCTGCGCTTTGGTTTCTTCGCCAATCATCCGCTGAAACTCACGCACAATGGTCgg                    >  1:1795412/1‑67 (MQ=255)
                                         ctctTTCCAGAATCTGCGCTTTGGTTTCTTCGCCAATCATCCGCTGAAACTCACGCACAATGGTCgg                    >  1:2549761/1‑67 (MQ=255)
                                         ctctTTCCAGAATCTGCGCTTTGGTTTCTTCGCCAATCATCCGCTGAAACTCACGCACAATGGTCgg                    >  1:3223489/1‑67 (MQ=255)
                                         ctctTTCCAGAATCTGCGCTTTGGTTTCTTCGCCAATCATCCGCTGAAACTCACGCACAATGGTCgg                    >  1:3123649/1‑67 (MQ=255)
                                         ctctTTCCAGAATCTGCGCTTTGGTTTCTTCGCCAATCATCCGCTGAAACTCACGCACAATGGTCgg                    >  1:669346/1‑67 (MQ=255)
                                         ctctTTCCAGAATCTGCGCTTTGGTTTCTTCGCCAATCATCCGCTGAAACTCACGCACAATGGTCg                     >  1:611930/1‑66 (MQ=255)
                                         ctctTTCCAGAATCTGCGCTTTGGTTTCTTCGCCAATCATCCGCTGAAACTCACGCACAATGGTCg                     >  1:3173014/1‑66 (MQ=255)
                                           ctTTCCAGAATCTGCGCTTTGGTTTCTTCGCCAATCATCCGCTGAAACTCACGCACAATGGTCGGATa                 >  1:1026702/1‑68 (MQ=255)
                                           ctTTCCAGAATCTGCGCTTTGGTTTCTTCGCCAATCATCCGCTGAAACTCACGCACAATGGTCGGATa                 >  1:354164/1‑68 (MQ=255)
                                           ctTTCCAGAATCTGCGCTTTGGTTTCTTCGCCAATCATCCGCTGAAACTCACGCACAATGGTCGGATa                 >  1:32825/1‑68 (MQ=255)
                                           ctTTCCAGAATCTGCGCTTTGGTTTCTTCGCCAATCATCCGCTGAAACTCACGCACAATGGTCGGATa                 >  1:2547592/1‑68 (MQ=255)
                                           ctTTCCAGAATCTGCGCTTTGGTTTCTTCGCCAATCATCCGCTGAAACTCACGCACAATGGTCGGATa                 >  1:121515/1‑68 (MQ=255)
                                           ctTTCCAGAATCTGCGCTTTCGTTTCTTCGCCAATCATCCGCTGAAACTCACGCACaa                           <  1:687834/58‑1 (MQ=255)
                                           ctTTCCAGAATCTGCGCTTTCGTTTCTTCGCCAATCATCCGCTGAAACTCACGCACaa                           <  1:668135/58‑1 (MQ=255)
                                           ctTTCCAGAATCTGCGCTTTCGTTTCTTCGCCAATCATCCGCTGAAACTCACGCACaa                           <  1:717710/58‑1 (MQ=255)
                                           ctTTCCAGAATCTGCGCTTTCGTTTCTTCGCCAATCATCCGCTGAAACTCACGCACaa                           <  1:957599/58‑1 (MQ=255)
                                           ctTTCCAGAATCTGCGCTTTCGTTTCTTCGCCAATCATCCGCTGAAACTCACGCACaa                           <  1:314055/58‑1 (MQ=255)
                                           ctTTCCAGAATCTGCGCTTTCGTTTCTTCGCCAATCATCCGCTGAAACTCACGCACaa                           <  1:3338803/58‑1 (MQ=255)
                                           ctTTCCAGAATCTGCGCTTTCGTTTCTTCGCCAATCATCCGCTGAAACTCACGCACaa                           <  1:898347/58‑1 (MQ=255)
                                           ctTTCCAGAATCTGCGCTTTCGTTTCTTCGCCAATCATCCGCTGAAACTCACGCACaa                           <  1:498831/58‑1 (MQ=255)
                                           ctTTCCAGAATCTGCGCTTTCGTTTCTTCGCCAATCATCCGCTGAAACTCACGCACaa                           <  1:457583/58‑1 (MQ=255)
                                           ctTTCCAGAATCTGCGCTTTCGTTTCTTCGCCAATCATCCGCTGAAACTCACGCACaa                           <  1:3361122/58‑1 (MQ=255)
                                           ctTTCCAGAATCTGCGCTTTCGTTTCTTCGCCAATCATCCGCTGAAACTCACGCACaa                           <  1:14691/58‑1 (MQ=255)
                                           ctTTCCAGAATCTGCGCTTTCGTTTCTTCGCCAATCATCCGCTGAAACTCACGCACaa                           <  1:1801552/58‑1 (MQ=255)
                                           ctTTCCAGAATCTGCGCTTTCGTTTCTTCGCCAATCATCCGCTGAAACTCACGCACaa                           <  1:3006082/58‑1 (MQ=255)
                                           ctTTCCAGAATCTGCGCTTTCGTTTCTTCGCCAATCATCCGCTGAAACTCACGCACaa                           <  1:1618605/58‑1 (MQ=255)
                                           ctTTCCAGAATCTGCGCTTTCGTTTCTTCGCCAATCATCCGCTGAAACTCACGCACaa                           <  1:1617496/58‑1 (MQ=255)
                                           ctTTCCAGAATCTGCGCTTTCGTTTCTTCGCCAATCATCCGCTGAAACTCACGCACaa                           <  1:2511593/58‑1 (MQ=255)
                                           ctTTCCAGAATCTGCGCTTTCGTTTCTTCGCCAATCATCCGCTGAAACTCACGCACaa                           <  1:2524222/58‑1 (MQ=255)
                                           ctTTCCAGAATCTGCGCTTTCGTTTCTTCGCCAATCATCCGCTGAAACTCACGCACaa                           <  1:1468057/58‑1 (MQ=255)
                                           ctTTCCAGAATCTGCGCTTTCGTTTCTTCGCCAATCATCCGCTGAAACTCACGCACaa                           <  1:2606030/58‑1 (MQ=255)
                                           ctTTCCAGAATCTGCGCTTTCGTTTCTTCGCCAATCATCCGCTGAAACTCACGCACaa                           <  1:1389541/58‑1 (MQ=255)
                                           ctTTCCAGAATCTGCGCTTTCGTTTCTTCGCCAATCATCCGCTGAAACTCACGCACaa                           <  1:2768263/58‑1 (MQ=255)
                                           ctTTCCAGAATCTGCGCTTTCGTTTCTTCGCCAATCATCCGCTGAAACTCACGCACaa                           <  1:2785519/58‑1 (MQ=255)
                                           ctTTCCAGAATCTGCGCTTTCGTTTCTTCGCCAATCATCCGCTGAAACTCACGCACaa                           <  1:1182423/58‑1 (MQ=255)
                                           ctTTCCAGAATCTGCGCTTTCGTTTCTTCGCCAATCATCCGCTGAAACTCACGCACaa                           <  1:2880203/58‑1 (MQ=255)
                                           ctTTCCAGAATCTGCGCTTTCGTTTCTTCGCCAATCATCCGCTGAAACTCACGCACaa                           <  1:2953942/58‑1 (MQ=255)
                                           ctTTCCAGAATCTGCCCTTTCGTTTCTTCGCCAATCATCCGCTGAAACTCACGCACaa                           <  1:410019/58‑1 (MQ=255)
                                            tttCCAGAATCTGCGCTTTCGTTTCTTCGCCAATCATCCGCTGAAACTCACGCACaa                           <  1:3141090/57‑1 (MQ=255)
                                              tCCAGAATCTGCGCTTTGGTTTCTTCGCCAATCATCCGCTGAAACTCACGCACAATGGTCGGATAAg               <  1:428405/67‑1 (MQ=255)
                                              tCCAGAATCTGCGCTTTGGTTTCTTCGCCAATCATCCGCTGAAACTCACGCACAATGGTCGGATAAg               <  1:734811/67‑1 (MQ=255)
                                              tCCAGAATCTGCGCTTTGGTTTCTTCGCCAATCATCCGCTGAAACTCACGCACAATGGTCGGATAAg               <  1:1622674/67‑1 (MQ=255)
                                              tCCAGAATCTGCGCTTTGGTTTCTTCGCCAATCATCCGCTGAAACTCACGCACAATGGTCGGATAAg               <  1:1670819/67‑1 (MQ=255)
                                              tCCAGAATCTGCGCTTTGGTTTCTTCGCCAATCATCCGCTGAAACTCACGCACAATGGTCGGATAAg               <  1:666592/67‑1 (MQ=255)
                                              tCCAGAATCTGCGCTTTGGTTTCTTCGCCAATCATCCGCTGAAACTCACGCACAATGGTCGGATAAg               <  1:1944067/67‑1 (MQ=255)
                                              tCCAGAATCTGCGCTTTGGTTTCTTCGCCAATCATCCGCTGAAACTCACGCACAATGGTCGGATAAg               <  1:2788297/67‑1 (MQ=255)
                                              tCCAGAATCTGCGCTTTGGTTTCTTCGCCAATCATCCGCTGAAACTCACGCACAATGGTCGGATAAg               <  1:2726247/67‑1 (MQ=255)
                                              tCCAGAATCTGCGCTTTGGTTTCTTCGCCAATCATCCGCTGAAACTCACGCACAATGGTCGGATAAg               <  1:3259529/67‑1 (MQ=255)
                                              tCCAGAATCTGCGCTTTGGTTTCTTCGCCAATCATCCGCTGAAACTCACGCACAATGGTCGGATAAg               <  1:3272845/67‑1 (MQ=255)
                                              tCCAGAATCTGCGCTTTGGTTTCTTCGCCAATCATCCGCTGAAACTCACGCACAATGGTCGGATAAg               <  1:2199364/67‑1 (MQ=255)
                                              tCCAGAATCTGCGCTTTGGTTTCTTCGCCAATCATCCGCTGAAACTCACGCACAATGGTCGGATAAg               <  1:3393929/67‑1 (MQ=255)
                                              tCCAGAATCTGCGCTTTGGTTTCTTCGCCAATCATCCGCTGAAACTCACGCACAATGGTCGGATAAg               <  1:2101824/67‑1 (MQ=255)
                                              tCCAGAATCTGCGCTTTGGTTTCTTCGCCAATCATCCGCTGAAACTCACGCACAATGGTCGGATAAg               <  1:459865/67‑1 (MQ=255)
                                                    aTCTGCGCTTTGGTTTCTTCGCCAATCATCCGCTGAAACt                                    <  1:3275166/40‑1 (MQ=255)
                                                    aTCTGCGCTTTGGTTTCTTCGCCAATCATCCGCTGAAACTCACGCACAATGGTCg                     >  1:1420537/1‑55 (MQ=255)
                                                      cTGCGCTTTGGTTTCTTCGCCAATCATCCGCTGAAACTCACGCACAATGGTCGGTTAAGGATGCGGGcc     <  1:2167632/69‑1 (MQ=255)
                                                      cTGCGCTTTGGTTTCTTCGCCAATCATCCGCTGAAACTCACGCACAATGGTCGGATAAGGATGCGGGcc     <  1:642825/69‑1 (MQ=255)
                                                      cTGCGCTTTGGTTTCTTCGCCAATCATCCGCTGAAACTCACGCACAATGGTCGGATAAGGATGCGGGcc     <  1:1765625/69‑1 (MQ=255)
                                                      cTGCGCTTTGGTTTCTTCGCCAATCATCCGCTGAAACTCACGCACAATGGTCGGATAAGGATGCGGGcc     <  1:597885/69‑1 (MQ=255)
                                                      cTGCGCTTTGGTTTCTTCGCCAATCATCCGCTGAAACTCACGCACAATGGTCGGATAAGGATGCGGGcc     <  1:1769587/69‑1 (MQ=255)
                                                      cTGCGCTTTGGTTTCTTCGCCAATCATCCGCTGAAACTCACGCACAATGGTCGGATAAGGATGCGGGcc     <  1:511705/69‑1 (MQ=255)
                                                      cTGCGCTTTGGTTTCTTCGCCAATCATCCGCTGAAACTCACGCACAATGGTCGGATAAGGATGCGGGcc     <  1:1128569/69‑1 (MQ=255)
                                                           cTTTGGTTTCTTCGCCAATCATCCGCTGAAACTCACGCACAATGGTCGGATAAGGATGCGGGCCAGc  <  1:2262576/67‑1 (MQ=255)
                                                               |                                                              
GCCGCCAACACAGGCGATAACGGCATCCGGCAGGCGACCTTCTCTTTCCAGAATCTGCGCTTTGGTTTCTTCGCCAATCATCCGCTGAAACTCACGCACAATGGTCGGATAAGGATGCGGGCCAGC  >  minE/856671‑856796

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: