Predicted mutation
evidence seq id position mutation freq annotation gene description
RA minE 2,872,350 C→G 22.2% Q201E (CAA→GAA)  ulaF → L‑ribulose 5‑phosphate 4‑epimerase

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*minE2,872,3500CG22.2% 207.2 / 53.5 117Q201E (CAA→GAA) ulaFL‑ribulose 5‑phosphate 4‑epimerase
Reads supporting (aligned to +/- strand):  ref base C (54/37);  new base G (19/7);  total (73/44)
Fisher's exact test for biased strand distribution p-value = 2.54e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.95e-01

Only 100 of 117 total aligned reads displayed.

TAACGCGGTGGTGATGGAAGAAGTGGCGAAAATGGCGTGGATTGCCCGCGGCATTAACCCACAACTCAATCACATCGACAGCTTCCTGATGAATAAACACTTCATGCGTAAACACGGTCCTAAC  >  minE/2872284‑2872407
                                                                  |                                                         
tAACGCGGTGGTGATGGAAGAAGTGGCGAAAATGGCGTGGATTGCCCGCGGCATTAACCCACAACTCa                                                          >  1:2082913/1‑68 (MQ=255)
tAACGCGGTGGTGATGGAAGAAGTGGCGAAAATGGCGTGGATTGCCCGCGGCATTAACCCACAACTCa                                                          >  1:2534476/1‑68 (MQ=255)
tAACGCGGTGGTGATGGAAGAAGTGGCGAAAATGGCGTGGATTGCCCGCGGCATTAACCCACAACTCa                                                          >  1:2620190/1‑68 (MQ=255)
tAACGCGGTGGTGATGGAAGAAGTGGCGAAAATGGCGTGGATTGCCCGCGGCATTAACCCACAACTCa                                                          >  1:2622199/1‑68 (MQ=255)
tAACGCGGTGGTGATGGAAGAAGTGGCGAAAATGGCGTGGATTGCCCGCGGCATTAACCCACAACTCa                                                          >  1:2669394/1‑68 (MQ=255)
tAACGCGGTGGTGATGGAAGAAGTGGCGAAAATGGCGTGGATTGCCCGCGGCATTAACCCACAACTCa                                                          >  1:3082803/1‑68 (MQ=255)
tAACGCGGTGGTGATGGAAGAAGTGGCGAAAATGGCGTGGATTGCCCGCGGCATTAACCCACAACTCa                                                          >  1:325488/1‑68 (MQ=255)
tAACGCGGTGGTGATGGAAGAAGTGGCGAAAATGGCGTGGATTGCCCGCGGCATTAACCCACAACTCa                                                          >  1:373853/1‑68 (MQ=255)
tAACGCGGTGGTGATGGAAGAAGTGGCGAAAATGGCGTGGATTGCCCGCGGCATTAACCCACAACTCa                                                          >  1:507081/1‑68 (MQ=255)
tAACGCGGTGGTGATGGAAGAAGTGGCGAAAATGGCGTGGATTGCCCGCGGCATTAACCCACAACTCa                                                          >  1:1253460/1‑68 (MQ=255)
tAACGCGGTGGTGATGGAAGAAGTGGCGAAAATGGCGTGGATTGCCCGCGGCATTAACCCACAACTCa                                                          >  1:532612/1‑68 (MQ=255)
tAACGCGGTGGTGATGGAAGAAGTGGCGAAAATGGCGTGGATTGCCCGCGGCATTAACCCACAACTCa                                                          >  1:684853/1‑68 (MQ=255)
  aCGCGGTGGTGATGGAAGAAGTGGCTAAAATGGCGTGGATTGCCCGCGGCATTAACCCACAACTCAAt                                                        <  1:1541384/68‑1 (MQ=255)
  aCGCGGTGGTGATGGAAGAAGTGGCGAAAATGGCGTGGATTGCCCGCGGCATTAACCCACAACTCAAt                                                        <  1:2335985/68‑1 (MQ=255)
  aCGCGGTGGTGATGGAAGAAGTGGCGAAAATGGCGTGGATTGCCCGCGGCATTAACCCACAACTCAAt                                                        <  1:1159805/68‑1 (MQ=255)
  aCGCGGTGGTGATGGAAGAAGTGGCGAAAATGGCGTGGATTGCCCGCGGCATTAACCCACAACTCAAt                                                        <  1:586113/68‑1 (MQ=255)
  aCGCGGTGGTGATGGAAGAAGTGGCGAAAATGGCGTGGATTGCCCGCGGCATTAACCCACAACTCAAt                                                        <  1:632898/68‑1 (MQ=255)
       gtggtgATGGAAGAAGTGGCGAAAATGGCGTGGATTGCCCGCGGCATTAACCCACAACTCAATCACATc                                                  <  1:898143/69‑1 (MQ=255)
       gtggtgATGGAAGAAGTGGCGAAAATGGCGTGGATTGCCCGCGGCATTAACCCACAACTCAATCACATc                                                  <  1:1896392/69‑1 (MQ=255)
       gtggtgATGGAAGAAGTGGCGAAAATGGCGTGGATTGCCCGCGGCATTAACCCACAACTCAATCACATc                                                  <  1:3264715/69‑1 (MQ=255)
       gtggtgATGGAAGAAGTGGCGAAAATGGCGTGGATTGCCCGCGGCATTAACCCACAACTCAATCACATc                                                  <  1:385272/69‑1 (MQ=255)
       gtggtgATGGAAGAAGTGGCGAAAATGGCGTGGATTGCCCGCGGCATTAACCCACAACTCAATCACATc                                                  <  1:464791/69‑1 (MQ=255)
       gtggtgATGGAAGAAGTGGCGAAAATGGCGTGGATTGCCCGCGGCATTAACCCACAACTCAATCACATc                                                  <  1:1435984/69‑1 (MQ=255)
          gtgATGGAAGAAGTGGCGAAAATGGCGTGGATTGCCCGCGGCATTAACCCACAACTCAATCACATCGAc                                               >  1:874158/1‑69 (MQ=255)
          gtgATGGAAGAAGTGGCGAAAATGGCGTGGATTGCCCGCGGCATTAACCCACAACTCAATCACATCGAc                                               >  1:1961525/1‑69 (MQ=255)
               ggAAGAAGTGGCGAAAATGGCGTGGATTGCGCGTAGTCTTAACCCACAACTGAATCACATCGACAGCt                                           <  1:760691/68‑1 (MQ=255)
               ggAAGAAGTGGCGAAAATGGCGTGGATTGCGCGTAGTATTAACCCACAACTGAATCACATCGACAGCt                                           <  1:1667743/68‑1 (MQ=255)
               ggAAGAAGTGGCGAAAATGGCGTGGATTGCGCGTAGTATTAACCCACAACTGAATCACATCGACAGCt                                           <  1:1527108/68‑1 (MQ=255)
               ggAAGAAGTGGCGAAAATGGCGTGGATTGCGCGTAGTATTAACCCACAACTGAATCACATCGACAGCt                                           <  1:426996/68‑1 (MQ=255)
               ggAAGAAGTGGCGAAAATGGCGTGGATTGCGCGTAGTATTAACCCACAACTGAATCACATCGACAGCt                                           <  1:1432560/68‑1 (MQ=255)
               ggAAGAAGTGGCGAAAATGGCGTGGATTGCCCGCGGCATTAACCCACAACTCAATCACATCGACAGCtt                                          >  1:3194032/1‑69 (MQ=255)
               ggAAGAAGTGGCGAAAATGGCGTGGATTGCCCGCGGCATTAACCCACAACTCAATCACATCGACAGCtt                                          >  1:319836/1‑69 (MQ=255)
               ggAAGAAGTGGCGAAAATGGCGTGGATTGCCCGCGGCATTAACCCACAACTCAATCACATCGACAGCtt                                          >  1:2890721/1‑69 (MQ=255)
                  agaagTGGCGAAAATGGCGTTGATTGCCCGCGGCATTAACCCACAACTCAATCACATCGACAGCTTCCt                                       <  1:511440/69‑1 (MQ=255)
                  agaagTGGCGAAAATGGCGTGGATTGCCCGCGGCATTAACCCACAACTCAATCACATCGACAGCTTCCt                                       <  1:2115037/69‑1 (MQ=255)
                      gTGGCGAAAATGGCGTGGATTGCCCGCGGCATTAACCCACAACTCAATCACATCGACAGCTTCCtgatg                                   <  1:2280354/69‑1 (MQ=255)
                      gTGGCGAAAATGGCGTGGATTGCCCGCGGCATTAACCCACAACTCAATCACATCGACAGCTTCCtgatg                                   <  1:2393675/69‑1 (MQ=255)
                       tGGCGAAAATGGCGTGGATTGCCCGCGGCATTAACCCACAACTCAATCACATCGACAGCTTCCtgatg                                   <  1:1786254/68‑1 (MQ=255)
                            aaaaTGGCGTGGATTGCGCGTAGTATTAACCCACAACTGAATCACATCGACAGCt                                           <  1:865353/55‑1 (MQ=255)
                            aaaaTGGCGTGGATTGCCCGCGGCATTAACCCACAACTCAATCACATCGACAGCTTCCTGATGAATaaa                             <  1:1939814/69‑1 (MQ=255)
                            aaaaTGGCGTGGATTGCCCGCGGCATTAACCCACAACTCAATCACATCGACAGCTTCCTGATGAATaaa                             <  1:335799/69‑1 (MQ=255)
                            aaaaTGGCGTGGATTGCCCGCGGCATTAACCCACAACTCAATCACATCGACAGCTTCCTGATGAATaaa                             <  1:383840/69‑1 (MQ=255)
                            aaaaTGGCGTGGATTGCCCGCGGCATTAACCCACAACTCAATCACATCGACAGCTTCCTGATGAATaaa                             <  1:1214866/69‑1 (MQ=255)
                             aaaTGGCGTGGATTGCGCGTAGTATTAACCCACAACTGAATCACATCGACAGCt                                           >  1:1961135/1‑54 (MQ=255)
                             aaaTGGCGTGGATTGCGCGTAGTATTAACCCACAACTGAATCACATCGACAGCTTCCTGATGAATaaa                             >  1:353862/1‑68 (MQ=255)
                             aaaTGGCGTGGATTGCGCGTAGTATTAACCCACAACTGAATCACATCGACAGCTTCCTGATGAATAAac                            >  1:2834302/1‑69 (MQ=255)
                             aaaTGGCGTGGATTGCGCGTAGTATTAACCCACAACTGAATCACATCGACAGCTTCCTGATGAATAAac                            >  1:3114208/1‑69 (MQ=255)
                             aaaTGGCGTGGATTGCGCGTAGTATTAACCCACAACTGAATCACATCGACAGCTTCCTGATGAATAAac                            >  1:457761/1‑69 (MQ=255)
                             aaaTGGCGTGGATTGCGCGTAGTATTAACCCACAACTGAATCACATCGACAGCTTCCTGATGAATAAac                            >  1:548105/1‑69 (MQ=255)
                             aaaTGGCGTGGATTGCGCGTAGTATTAACCCACAACTGAATCACATCGACAGCTTCCTGATGAATAAac                            >  1:610579/1‑69 (MQ=255)
                             aaaTGGCGTGGATTGCGCGTAGTATTAACCCACAACTGAATCACATCGACAGCTTCCTGATGAATAAac                            >  1:800711/1‑69 (MQ=255)
                             aaaTGGCGTGGATTGCGCGTAGTATTAACCCACAACTGAATCACATCGACAGCTTCCTGATGAATAAac                            >  1:1003361/1‑69 (MQ=255)
                             aaaTGGCGTGGATTGCGCGTAGTATTAACCCACAACTGAATCACATCGACAGCTTCCTGATGAATAAac                            >  1:1871950/1‑69 (MQ=255)
                             aaaTGGCGTGGATTGCGCGTAGTATTAACCCACAACTGAATCACATCGACAGCTTCCTGATGAATAAac                            >  1:1771166/1‑69 (MQ=255)
                             aaaTGGCGTGGATTGCGCGTAGTATTAACCCACAACTGAATCACATCGACAGCTTCCTGATGAATAAac                            >  1:2355464/1‑69 (MQ=255)
                             aaaTGGCGTGGATTGCGCGTAGTATTAACCCACAACTGAATCACATCGACAGCTTCCTGATGAATAAac                            >  1:2168973/1‑69 (MQ=255)
                             aaaTGGCGTGGATTGCGCGTAGTATTAACCCACAACTGAATCACATCGACAGCTTCCTGATGAATAAac                            >  1:2532198/1‑69 (MQ=255)
                             aaaTGGCGTGGATTGCGCGTAGTATTAACCCACAACTGAATCACATCGACAGCTTCCTGATGAATAAac                            >  1:2083952/1‑69 (MQ=255)
                             aaaTGGCGTGGATTGCGCGTAGTATTAACCCACAACTGAATCACATCGACAGCTTCCTGATGAATAAac                            >  1:1906182/1‑69 (MQ=255)
                             aaaTGGCGTGGATTGCGCGTAGTATTAACCCACAACTGAATCACATCGACAGCTTCCTGATGAATAAac                            >  1:2614956/1‑69 (MQ=255)
                             aaaTGGCGTGGATTGCGCGTAGTATTAACCCACAACTGAATCACATCGACAGCTTCCTGATGAATAAac                            >  1:2029631/1‑69 (MQ=255)
                             aaaTGGCGTGGATTGCCCGCGGCATTAACCCACAACTCAATCACATCGACAGCTTCCTGATGAATaaa                             <  1:1559420/68‑1 (MQ=255)
                                 ggCGTGGATTGCCCGCGGCATTAACCCACAACTCAATCACATCGACAGCTTCCTGATGAATAAACACtt                        <  1:1765879/69‑1 (MQ=255)
                                  gCGTGGATTGCCCGCGGCATTAACCCACAACTCAATCACATCGACAGCTTCCTGATGAATAAACACTTc                       >  1:2859824/1‑69 (MQ=255)
                                       gATTGCCCGCGGCATTAACCCACAACTCAATCACATCGACAGCTTCCTGATGAATAAACACTTCATGc                   >  1:764542/1‑68 (MQ=255)
                                       gATTGCCCGCGGCATTAACCCACAACTCAATCACATCGACAGCTTCCTGATGAATAAACACTTCATGc                   >  1:1456525/1‑68 (MQ=255)
                                       gATTGCCCGCGGCATTAACCCACAACTCAATCACATCGACAGCTTCCTGATGAATAAACACTTCATGc                   >  1:1092255/1‑68 (MQ=255)
                                       gATTGCCCGCGGCATTAACCCACAACTCAATCACATCGACAGCTTCCTGATGAATAAACACTTCATGc                   >  1:827475/1‑68 (MQ=255)
                                       gATTGCCCGCGGCATTAACCCACAACTCAATCACATCGACAGCTTCCTGATGAATAAACACTTCATGc                   >  1:843800/1‑68 (MQ=255)
                                       gATTGCCCGCGGCATTAACCCACAACTCAATCACATCGACAGCTTCCTGATGAATAAACACTTCATGc                   >  1:3074752/1‑68 (MQ=255)
                                       gATTGCCCGCGGCATTAACCCACAACTCAATCACATCGACAGCTTCCTGATGAATAAACACTTCATGc                   >  1:1582685/1‑68 (MQ=255)
                                       gATTGCCCGCGGCATTAACCCACAACTCAATCACATCGACAGCTTCCTGATGAATAAACACTTCATGc                   >  1:2852120/1‑68 (MQ=255)
                                        aTTGCCCGCGGCATTAACCCACAACTCAATCACATc                                                  <  1:1808642/36‑1 (MQ=255)
                                              cgcgGCATTAACCCACAACTCAATCACATCGACAGCTTCCTGATGAATAAACACTTCATGCGTAAaca            <  1:2776126/68‑1 (MQ=255)
                                                 ggCATTAACCCACAACTCAATCACATCGACAGCTTCCtgatg                                   <  1:3398463/42‑1 (MQ=255)
                                                     ttAACCCACAACTCAATCACATCGACAGCTTCCTGATGAATAAACACTTCATGCGTAAACACGGTCCt     >  1:2577050/1‑68 (MQ=255)
                                                     ttAACCCACAACTCAATCACATCGACAGCTTCCTGATGAATAAACACTTCATGCGTAAACACGGTCCTa    >  1:2435986/1‑69 (MQ=255)
                                                     ttAACCCACAACTCAATCACATCGACAGCTTCCTGATGAATAAACACTTCATGCGTAAACACGGTCCTa    >  1:1125122/1‑69 (MQ=255)
                                                     ttAACCCACAACTCAATCACATCGACAGCTTCCTGATGAATAAACACTTCATGCGTAAACACGGTCCTa    >  1:869752/1‑69 (MQ=255)
                                                     ttAACCCACAACTCAATCACATCGACAGCTTCCTGATGAATAAACACTTCATGCGTAAACACGGTCCTa    >  1:2171385/1‑69 (MQ=255)
                                                     ttAACCCACAACTCAATCACATCGACAGCTTCCTGATGAATAAACACTTCATGCGTAAACACGGTCCTa    >  1:516814/1‑69 (MQ=255)
                                                     ttAACCCACAACTCAATCACATCGACAGCTTCCTGATGAATAAACACTTCATGCGTAAACACGGTCCTa    >  1:1463683/1‑69 (MQ=255)
                                                     ttAACCCACAACTCAATCACATCGACAGCTTCCTGATGAATAAACACTTCATGCGTAAACACGGTCCTa    >  1:2766257/1‑69 (MQ=255)
                                                     ttAACCCACAACTCAATCACATCGACAGCTTCCTGATGAATAAACACTTCATGCGTAAACACGGTCCTa    >  1:2998610/1‑69 (MQ=255)
                                                     ttAACCCACAACTCAATCACATCGACAGCTTCCTGATGAATAAACACTTCATGCGTAAACACGGTCCTa    >  1:3117882/1‑69 (MQ=255)
                                                     ttAACCCACAACTCAATCACATCGACAGCTTCCTGATGAATAAACACTTCATGCGTAAACACGGTCCTa    >  1:3263422/1‑69 (MQ=255)
                                                     ttAACCCACAACTCAATCACATCGACAGCTTCCTGATGAATAAACACTTCATGCGTAAACACGGTCCTa    >  1:2559160/1‑69 (MQ=255)
                                                     ttAACCCACAACTCAATCACATCGACAGCTTCCTGATGAATAAACACTTCATGCGTAAACACGGTCCTa    >  1:3423577/1‑69 (MQ=255)
                                                     ttAACCCACAACTCAATCACATCGACAGCTTCCTGATGAATAAACACTTCATGCGTAAACACGGTCCTa    >  1:2463852/1‑69 (MQ=255)
                                                     ttAACCCACAACTCAATCACATCGACAGCTTCCTGATGAATAAACACTTCATGCGTAAACACGGTCCTa    >  1:445489/1‑69 (MQ=255)
                                                     ttAACCCACAACTCAATCACATCGACAGCTTCCTGATGAATAAACACTTCATGCGTAAACACGGTCCTa    >  1:1350783/1‑69 (MQ=255)
                                                       aaCCCACAACTCAATCACATCGACAGCTTCCTGATGAAt                                >  1:871131/1‑39 (MQ=255)
                                                       aaCCCACAACTCAATCACATCGACAGCTTCCTGATGAATAAACACTTCATGCGTAAACACGGTCCTAAc  >  1:639315/1‑69 (MQ=255)
                                                       aaCCCACAACTCAATCACATCGACAGCTTCCTGATGAATAAACACTTCATGCGTAAACACGGTCCTAAc  >  1:373733/1‑69 (MQ=255)
                                                       aaCCCACAACTCAATCACATCGACAGCTTCCTGATGAATAAACACTTCATGCGTAAACACGGTCCTAAc  >  1:2187940/1‑69 (MQ=255)
                                                       aaCCCACAACTCAATCACATCGACAGCTTCCTGATGAATAAACACTTCATGCGTAAACACGGTCCTAAc  >  1:1137340/1‑69 (MQ=255)
                                                       aaCCCACAACTCAATCACATCGACAGCTTCCTGATGAATAAACACTTCATGCGTAAACACGGTCCTAAc  >  1:1842881/1‑69 (MQ=255)
                                                       aaCCCACAACTCAATCACATCGACAGCTTCCTGATGAATAAACACTTCATGCGTAAACACGGTCCTAAc  >  1:1898744/1‑69 (MQ=255)
                                                       aaCCCACAACTCAATCACATCGACAGCTTCCTGATGAATAAACACTTCATGCGTAAACACGGTCCTAAc  >  1:2948805/1‑69 (MQ=255)
                                                       aaCCCACAACTCAATCACATCGACAGCATCCTGATGAATAAACACTTCATGCGTAAACACGGTCCTAAc  >  1:1616755/1‑69 (MQ=255)
                                                                  |                                                         
TAACGCGGTGGTGATGGAAGAAGTGGCGAAAATGGCGTGGATTGCCCGCGGCATTAACCCACAACTCAATCACATCGACAGCTTCCTGATGAATAAACACTTCATGCGTAAACACGGTCCTAAC  >  minE/2872284‑2872407

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: