Predicted mutation
evidence seq id position mutation freq annotation gene description
RA minE 940,401 G→A 44.6% G124G (GGC→GGT sfcA ← malate dehydrogenase, NAD‑requiring

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*minE940,4010GA44.6% 13.6 / 49.9 65G124G (GGC→GGTsfcAmalate dehydrogenase, NAD‑requiring
Reads supporting (aligned to +/- strand):  ref base G (10/26);  new base A (19/10);  total (29/36)
Fisher's exact test for biased strand distribution p-value = 2.97e-03
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.60e-01

GTTCGGCACGTTTTGCAGAATATCGTCCATATTGTGCCGGTTCTGGTAAGAGATAAACACGCCGCGTGAACGGCGGTAGATCTCAGAAAAACGCTCACAGGCTGCGCCGACGGTTGGGGTATAAATAA  >  minE/940341‑940468
                                                            |                                                                   
gTTCGGCACGTTTTGCAGAATATCGTCCATATTGTGCCGGTTCTGGTAAGAGATAAACACGCCGCGTg                                                              <  1:2684202/68‑1 (MQ=255)
          ttttGCAGAATATCGTCCATATTGTGCCGGTTCTGGTAAGAGATAAACACGCCGCGTGAACGGCGGTAg                                                   <  1:1779543/69‑1 (MQ=255)
          ttttGCAGAATATCGTCCATATTGTGCCGGTTCTGGTAAGAGATAAACACGCCGCGTGAACGGCGGTAg                                                   <  1:297035/69‑1 (MQ=255)
          ttttGCAGAATATCGTCCATATTGTGCCGGTTCTGGTAAGAGATAAACACGCCGCGTGAACGGCGGTAg                                                   <  1:306775/69‑1 (MQ=255)
           tttGCAGAATATCGTCCATATTGTGCCGGTTCTGGTAAGAGATAAACACGCCGCGTGAACGGCGGTAg                                                   <  1:1122384/68‑1 (MQ=255)
           tttGCAGAATATCGTCCATATTGTGCCGGTTCTGGTAAGAGATAAACACGCCGCGTGAACGGCGGTAg                                                   <  1:3325951/68‑1 (MQ=255)
                  aaTATCGTCCATATTGTGCCGGTGCTGGTAAGAGATAAACACGCCGCGTGAACGGCGGTAg                                                   <  1:1371675/61‑1 (MQ=255)
                     atCGTCCATATTGTGCCGGTTCTGGTAAGAGATAAACACGCCGCGTGAACGGCGGTAGATCTCAGaaa                                         <  1:771632/68‑1 (MQ=255)
                      tCGTCCATATTGTGCCGGTTGTGGTAAGAGATAAACACGCCGCGTg                                                              <  1:1464977/46‑1 (MQ=255)
                      tCGTCCATATTGTGCCGGTTCTGGTAAGAGATAAACACGTCTCGTg                                                              <  1:5290/46‑1 (MQ=255)
                      tCGTCCATATTGTGCCGGTTCTGGTAAGAGATAAACACGCCGCGTg                                                              <  1:3136363/46‑1 (MQ=255)
                      tCGTCCATATTGTGCCGGTTCTGGTAAGAGATAAACACGCCGCGTg                                                              <  1:2663942/46‑1 (MQ=255)
                      tCGTCCATATTGTGCCGGTTCTGGTAAGAGATAAACACGCCGCGTg                                                              <  1:1950365/46‑1 (MQ=255)
                      tCGTCCATATTGTGCCGGTTCTGGTAAGAGATAAACACGCCGCGTg                                                              <  1:1885964/46‑1 (MQ=255)
                      tCGTCCATATTGTGCCGGTTCTGGTAAGAGATAAACACGCCGCGTg                                                              <  1:22132/46‑1 (MQ=255)
                      tCGTCCATATTGTGCCGGTTCTGGTAAGAGATAAACACGCCGCGTGAACGGCGGTAGATCTCAGaaaaa                                       >  1:2031424/1‑69 (MQ=255)
                      tCGTCCATATTGTGCCGGTTCTGGTAAGAGATAAACACGCCGCCTg                                                              <  1:764529/46‑1 (MQ=255)
                        gTCCATATTGTGCCGGTTCTGGTATGAGATAAACACGCCGCGTg                                                              <  1:2140125/44‑1 (MQ=255)
                             tatTGTGCCGGTTCTGGTAAGAGATAACCACGCCGCGTGAACGGCGGT‑‑ATCTCAGAAAAACGCTCACAg                              <  1:3348734/69‑1 (MQ=255)
                             tatTGTGCCGGTTCTGGTAAGAGATAAACACGCCGCGTGAACGGCGGTAGATCTCAGAAAAACGCTca                                 <  1:2122315/68‑1 (MQ=255)
                               ttGTGCCGGTTCTGGTAAGAGATAAACACGCCGCGTg                                                              <  1:2902707/37‑1 (MQ=255)
                               ttGTGCCGGTTCTGGTAAGAGATAAACACGCCGCGTGAACGGCGGTAg                                                   >  1:1033790/1‑48 (MQ=255)
                               ttGTGCCGGTTCTGGTAAGAGATAAACACGCCGCGTGAACGGCGGTAg                                                   >  1:2322455/1‑48 (MQ=255)
                               ttGTGCCGGTTCTGGTAAGAGATAAACACGCCGCGTGAACGGCGGTAg                                                   >  1:2366640/1‑48 (MQ=255)
                               ttGTGCCGGTTCTGGTAAGAGATAAACACGCCGCGTGAACGGCGGTAg                                                   >  1:504910/1‑48 (MQ=255)
                               ttGTGCCGGTTCTGGTAAGAGATAAACACGCCGCGTGAACGGCGGTAg                                                   >  1:1675100/1‑48 (MQ=255)
                               ttGTGCCGGTTCTGGTAAGAGATAAACACGCCGCGTGAACGGCGGTAg                                                   >  1:1610270/1‑48 (MQ=255)
                                        ttCTGGTAAGAGATAAACACACCGCGTGAACGGCGGTAGATCTCGGAAAAACGCTCACAGGCTGCGCCg                     <  1:1782487/69‑1 (MQ=255)
                                        ttCTGGTAAGAGATAAACACACCGCGTGAACGGCGGTAGATCTCAGAAAAACGCTCACAGGCTGCGCCg                     <  1:378241/69‑1 (MQ=255)
                                        ttCTGGTAAGAGATAAACACACCGCGTGAACGGCGGTAGATCTCAGAAAAACGCTCACAGGCTGCGCCg                     <  1:3377313/69‑1 (MQ=255)
                                        ttCTGGTAAGAGATAAACACACCGCGTGAACGGCGGTAGATCTCAGAAAAACGCTCACAGGCTGCGCCg                     <  1:3297795/69‑1 (MQ=255)
                                        ttCTGGTAAGAGATAAACACACCGCGTGAACGGCGGTAGATCTCAGAAAAACGCTCACAGGCTGCGCCg                     <  1:3272769/69‑1 (MQ=255)
                                        ttCTGGTAAGAGATAAACACACCGCGTGAACGGCGGTAGATCTCAGAAAAACGCTCACAGGCTGCGCCg                     <  1:2928610/69‑1 (MQ=255)
                                        ttCTGGTAAGAGATAAACACACCGCGTGAACGGCGGTAGATCTCAGAAAAACGCTCACAGGCTGCGCCg                     <  1:1578210/69‑1 (MQ=255)
                                         tCTGGTAAGAGATAAACACGCCGCGTGAACGGCGGTAGATCTCAGAAAAACGCt                                   <  1:2138121/54‑1 (MQ=255)
                                         tCTGGTAAGAGATAAACACGCCGCGTGAACGGCGGTAGATCTCAGAAAAACGCt                                   <  1:683760/54‑1 (MQ=255)
                                         tCTGGTAAGAGATAAACACGCCGCGTGAACGGCGGTAGATCTCAGAAAAACGCt                                   <  1:1136567/54‑1 (MQ=255)
                                         tCTGGTAAGAGATAAACACGCCGCGTGAACGGCGGTAGATCTCAGAAAAACGCt                                   <  1:3083983/54‑1 (MQ=255)
                                         tCTGGTAAGAGATAAACACACCGCGTGAACGGCGGTAGATCTCAGAAAAACGCTCACAGGCTGCGCCg                     <  1:1323048/68‑1 (MQ=255)
                                         tCTGGTAAGAGATAAACACACCGCGTGAACGGCGGTAGATCTCAGAAAAACGCTCACAGGCTGCGCCg                     <  1:2349789/68‑1 (MQ=255)
                                           tGGTAAGAGATAAACACGCCGCGTGAACGGCGGTAGATCTCAGaaa                                         <  1:2567149/46‑1 (MQ=255)
                                                agagATAAACACGCCGCGTGAACGGCGGTAGATCTCAGAAAAACGCTCACAGGCTGCGCCGACGGTTgg             >  1:1622742/1‑69 (MQ=255)
                                                agagATAAACACACCGCGTGAACGGCGGTAGATCTCAGAAAAACGCTCACAGGCTGCGCCg                     <  1:1700107/61‑1 (MQ=255)
                                                 gagaTAAACACACCGCGTGAACGGCGGTAGATCTCAGAAAAACGCTCACAGGCTGCGCCGACGGTTggg            >  1:1021384/1‑69 (MQ=255)
                                                 gagaTAAACACACCGCGTGAACGGCGGTAGATCTCAGAAAAACGCTCACAGGCTGCGCCGACGGTTggg            >  1:2301314/1‑69 (MQ=255)
                                                 gagaTAAACACACCGCGTGAACGGCGGTAGATCTCAGAAAAACGCTCACAGGCTGCGCCGACGGTTggg            >  1:2300423/1‑69 (MQ=255)
                                                 gagaTAAACACACCGCGTGAACGGCGGTAGATCTCAGAAAAACGCTCACAGGCTGCGCCGACGGTTggg            >  1:440232/1‑69 (MQ=255)
                                                 gagaTAAACACACCGCGTGAACGGCGGTAGATCTCAGAAAAACGCTCACAGGCTGCGCCGACGGTTggg            >  1:1160123/1‑69 (MQ=255)
                                                 gagaTAAACACACCGCGTGAACGGCGGTAGATCTCAGAAAAACGCTCACAGGCTGCGCCGACGGTTggg            >  1:3225969/1‑69 (MQ=255)
                                                 gagaTAAACACACCGCGTGAACGGCGGTAGATCTCAGAAAAACGCTCACAGGCTGCGCCGACGGTTggg            >  1:3182210/1‑69 (MQ=255)
                                                 gagaTAAACACACCGCGTGAACGGCGGTAGATCTCAGAAAAACGCTCACAGGCTGCGCCGACGGTTggg            >  1:2489315/1‑69 (MQ=255)
                                                 gagaTAAACACACCGCGTGAACGGCGGTAGATCTCAGAAAAACGCTCACAGGCTGCGCCGACGGTTggg            >  1:3040192/1‑69 (MQ=255)
                                                  agaTAAACACGCCGCGTGAACGGCGGTAGATCTCAGaaa                                         <  1:1085631/39‑1 (MQ=255)
                                                   gaTAAACACACCGCGTGAACGGCGGTAGATCTCAGAAAAACGCTCACAGGCTGCGCCGACGGTTggg            >  1:1867391/1‑67 (MQ=255)
                                                   gaTAAACACACCGCGTGAACGGCGGTAGATCTCAGAAAAACGCTCACAGGCTGCGCCGACGGTTggg            >  1:540292/1‑67 (MQ=255)
                                                   gaTAAACACACCGCGTGAACGGCGGTAGATCTCAGAAAAACGCTCACAGGCTGCGCCGACGGTTggg            >  1:607113/1‑67 (MQ=255)
                                                   gaTAAACACACCGCGTGAACGGCGGTAGATCTCAGAAAAACGCTCACAGGCTGCGCCGACGGTTggg            >  1:2866682/1‑67 (MQ=255)
                                                   gaTAAACACACCGCGTGAACGGCGGTAGATCTCAGAAAAACGCTCACAGGCTGCGCCGACGGTTggg            >  1:336988/1‑67 (MQ=255)
                                                   gaTAAACACACCGCGTGAACGGCGGTAGATCTCAGAAAAACGCTCACAGGCTGCGCCGACGGTTggg            >  1:3349484/1‑67 (MQ=255)
                                                   gaTAAACACACCGCGTGAACGGCGGTAGATCTCAGAAAAACGCTCACAGGCTGCGCCGACGGTTggg            >  1:136055/1‑67 (MQ=255)
                                                   gaTAAACACACCGCGTGAACGGCGGTAGATCTCAGAAAAACGCTCACAGGCTGCGCCGACGGTTggg            >  1:1083420/1‑67 (MQ=255)
                                                   gaTAAACACACCGCGTGAACGGCGGTAGATCTCAGAAAAACGATCACAGGCTGAGCCGACGGTTggg            >  1:1534223/1‑67 (MQ=255)
                                                   gaTAAAACACCGCGTGAACGGCGGTAGATCTCAGAAAAACGCTCACAGGCTGCGCCGACGGTTggg            >  1:1552546/1‑66 (MQ=255)
                                                       aaCACGCCGCGTGAACGGCGGTAGATCTCAGAAAAACGCTCACAGGCTGCGCCGACGGTTGGGGtata       >  1:3091033/1‑68 (MQ=255)
                                                       aaCACGCCGCGTGAACGGCGGTAGATCTCAGAAAAACGCTCACAGGCTGCGCCGACGGTTGGGGtata       >  1:280366/1‑68 (MQ=255)
                                                           cgccgcGTGAACGGCGGTAGATCTCAGAAAAACGCTCACAGGCTGCGCCGACGGTTGGGGTataaataa  <  1:203366/69‑1 (MQ=255)
                                                           cgccgcGTGAACGGCGGTAGATCTCAGAAAAACGCTCACAGGCTGCGCCGACGGTTGGGGTataaataa  <  1:892407/69‑1 (MQ=255)
                                                            |                                                                   
GTTCGGCACGTTTTGCAGAATATCGTCCATATTGTGCCGGTTCTGGTAAGAGATAAACACGCCGCGTGAACGGCGGTAGATCTCAGAAAAACGCTCACAGGCTGCGCCGACGGTTGGGGTATAAATAA  >  minE/940341‑940468

Alignment Legend
Aligned base mismatch/match (shaded by quality score): ATCG/ATCG < 3 ≤ ATCG/ATCG < 14 ≤ ATCG/ATCG < 32 ≤ ATCG/ATCG < 36 ≤ ATCG/ATCG
Unaligned base: atcg    Masked matching base: atcg    Alignment gap:     Deleted base: