Predicted mutation
evidence seq id position mutation freq annotation gene description
RA minE 190,114 T→C 62.6% F537F (TTT→TTC ldcC → lysine decarboxylase 2, constitutive

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*minE190,1140TC62.6% 48.7 / 70.0 112F537F (TTT→TTCldcClysine decarboxylase 2, constitutive
Reads supporting (aligned to +/- strand):  ref base T (17/25);  new base C (45/25);  total (62/50)
Fisher's exact test for biased strand distribution p-value = 1.86e-02
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.89e-01

Only 100 of 113 total aligned reads displayed.

GCAAAATTCCTCGACGAACGTGGGATCGTAGTAGAGAAAACCGGCCCTTATAACCTGCTGTTTCTCTTTAGTATTGGCATCGATAAAACCAAAGCAATGGGATTATTGCGTGGGTTGACGGAATTCAAACGC  >  minE/190046‑190177
                                                                    |                                                               
gCAAAATTCCTCGACGAACGTGGGATCGTAGTAGAGAAAACCGGCCCTTATAACCTGCTGTTTCTCttt                                                                 <  1:2359602/69‑1 (MQ=255)
    aaTACCTCGACGCACATGGCATCGTTGTTGAGAAAACCGGTCCGTATAACCTGCTGTTCCTGTTCAg                                                               >  1:3303229/1‑67 (MQ=255)
    aaTACCTCGACGAACATTGCATCGTTGTTGAGAAAACCGGTCCGTATAACCTGCTGTTCCTGTTCAgca                                                             <  1:206573/69‑3 (MQ=255)
    aaTACCTCGACGAACATGGCATCGTTGTTGAGAAAACCGGTCCGTATAACCTGCTGTTCCTGTTCAgca                                                             <  1:1468180/69‑3 (MQ=255)
    aaTACCTCGACGAACATGGCATCGTTGTTGAGAAAACCGGTCCGTATAACCTGCTGTTCCTGTTCAgca                                                             <  1:1904420/69‑3 (MQ=255)
    aaTACCTCGACGAACATGGCATCGTTGTTGAGAAAACCGGTCCGTATAACCTGCTGTTCCTGTTCAgca                                                             <  1:1888849/69‑3 (MQ=255)
    aaTACCTCGACGAACATGGCATCGTTGTTGAGAAAACCGGTCCGTATAACCTGCTGTTCCTGTTCAgca                                                             <  1:1846561/69‑3 (MQ=255)
    aaTACCTCGACGAACATGGCATCGTTGTTGAGAAAACCGGTCCGTATAACCTGCTGTTCCTGTTCAgca                                                             <  1:2072512/69‑3 (MQ=255)
    aaTACCTCGACGAACATGGCATCGTTGTTGAGAAAACCGGTCCGTATAACCTGCTGTTCCTGTTCAgca                                                             <  1:3077589/69‑3 (MQ=255)
    aaTACCTCGACGAACATGGCATCGTTGTTGAGAAAACCGGTCCGTATAACCTGCTGTTCCTGTTCAgca                                                             <  1:2926966/69‑3 (MQ=255)
    aaTACCTCGACGAACATGGCATCGTTGTTGAGAAAACCGGTCCGTATAACCTGCTGTTCCTGTTCAgca                                                             <  1:1607158/69‑3 (MQ=255)
    aaTACCTCGACGAACATGGCATCGTTGTTGAGAAAACCGGTCCGTATAACCTGCTGTTCCTGTTCAgca                                                             <  1:3224580/69‑3 (MQ=255)
    aaTACCTCGACGAACATGGCATCGTTGTTGAGAAAACCGGTCCGTATAACCTGCTGTTCCTGTTCAgca                                                             <  1:2151660/69‑3 (MQ=255)
    aaTACCTCGACGAACATGGCATCGTTGTTGAGAAAACCGGTCCGTATAACCTGCTGTTCCTGTTCAgca                                                             <  1:3314474/69‑3 (MQ=255)
    aaTACCTCGACGAACATGGCATCGTTGTTGAGAAAACCGGTCCGTATAACCTGCTGTTCCTGTTCAgca                                                             <  1:953241/69‑3 (MQ=255)
    aaTACCTCGACGAACATGGCATCGTTGTTGAGAAAACCGGTCCGTATAACCTGCTGTTCCTGTTCAgca                                                             <  1:3352934/69‑3 (MQ=255)
    aaTACCTCGACGAACATGGCATCGTTGTTGAGAAAACCGGTCCGTATAACCTGCTGTTCCTGTTCAgca                                                             <  1:1419785/69‑3 (MQ=255)
    aaTACCTCGACGAACATGGCATCGTTGTTGAGAAAACCGGTCCGTATAACCTGCTGTTCCTGTTCAgca                                                             <  1:2854087/69‑3 (MQ=255)
    aaTACCTCGACGAACATGGCATCGTTGTTGAGAAAACCGGTCCGTATAACCTGCTGTTCCTGTTCAgca                                                             <  1:2842915/69‑3 (MQ=255)
    aaTACCTCGACGAACATGGCATCGTTGTTGAGAAAACCGGTCCGTATAACCTGCTGTTCCTGTTCAgca                                                             <  1:2403666/69‑3 (MQ=255)
    aaTACCTCGACGAACATGGCATCGTTGTTGAGAAAACCGGTCCGTATAACCTGCTGTTCCTGTTCAgca                                                             <  1:870975/69‑3 (MQ=255)
    aaTACCTCGACGAACATGGCATCGTTGTTGAGAAAACCGGTCCGTATAACCTGCTGTTCCTGTTCAgca                                                             <  1:1233416/69‑3 (MQ=255)
    aaTACCTCGACGAACATGGCATCGTTGTTGAGAAAACCGGTCCGTATAACCTGCTGTTCCTGTTCAgca                                                             <  1:592151/69‑3 (MQ=255)
    aaTACCTCGACGAACATGGCATCGTTGTTGAGAAAACCGGTCCGTATAACCTGCTGTTCCTGTTCAgca                                                             <  1:616781/69‑3 (MQ=255)
    aaTACCTCGACGAACATGGCATCGTTGTTGAGAAAACCGGTCCGTATAACCTGCTGTTCCTGTTCAgca                                                             <  1:2668939/69‑3 (MQ=255)
    aaTACCTCGACGAACATGGCATCGTTGTTGAGAAAACCGGTCCGTATAACCTGCTGTTCCTGTTCAgca                                                             <  1:2663649/69‑3 (MQ=255)
    aaTACCTCGACGAACATGGCATCGTTGTTGAGAAAACCGGTCCGTATAACCTGCTGTTCCTGTTCAg                                                               >  1:209702/1‑67 (MQ=255)
    aaTACCTCGACGAACATGGCATCGTTGTTGAGAAAACCGGTCCGTATAACCTGCTGTTCCTGTTCAg                                                               >  1:2017844/1‑67 (MQ=255)
    aaTACCTCGACGAACATGGCATCGTTGTTGAGAAAACCGGTCCGTATAACCTGCTGTTCCTGTTCAg                                                               >  1:2225287/1‑67 (MQ=255)
    aaTACCTCGACGAACATGGCATCGTTGTTGAGAAAACCGGTCCGTATAACCTGCTGTTCCTGTTCAg                                                               >  1:211024/1‑67 (MQ=255)
    aaTACCTCGACGAACATGGCATCGTTGTTGAGAAAACCGGTCCGTATAACCTGCTGTTCCTGTTCAg                                                               >  1:2485818/1‑67 (MQ=255)
    aaTACCTCGACGAACATGGCATCGTTGTTGAGAAAACCGGTCCGTATAACCTGCTGTTCCTGTTCAg                                                               >  1:3028215/1‑67 (MQ=255)
    aaTACCTCGACGAACATGGCATCGTTGTTGAGAAAACCGGTCCGTATAACCTGCTGTTCCTGTTCAg                                                               >  1:1791176/1‑67 (MQ=255)
    aaTACCTCGACGAACATGGCATCGTTGTTGAGAAAACCGGTCCGTATAACCTGCTGTTCCTGTTCAg                                                               >  1:3183370/1‑67 (MQ=255)
    aaTACCTCGACGAACATGGCATCGTTGTTGAGAAAACCGGTCCGTATAACCTGCTGTTCCTGTTCAg                                                               >  1:3243803/1‑67 (MQ=255)
    aaTACCTCGACGAACATGGCATCGTTGTTGAGAAAACCGGTCCGTATAACCTGCTGTTCCTGTTCAg                                                               >  1:1528336/1‑67 (MQ=255)
    aaTACCTCGACGAACATGGCATCGTTGTTGAGAAAACCGGTCCGTATAACCTGCTGTTCCTGTTCAg                                                               >  1:3317817/1‑67 (MQ=255)
    aaTACCTCGACGAACATGGCATCGTTGTTGAGAAAACCGGTCCGTATAACCTGCTGTTCCTGTTCAg                                                               >  1:397673/1‑67 (MQ=255)
    aaTACCTCGACGAACATGGCATCGTTGTTGAGAAAACCGGTCCGTATAACCTGCTGTTCCTGTTCAg                                                               >  1:1307497/1‑67 (MQ=255)
    aaTACCTCGACGAACATGGCATCGTTGTTGAGAAAACCGGTCCGTATAACCTGCTGTTCCTGTTCAg                                                               >  1:1278889/1‑67 (MQ=255)
    aaTACCTCGACGAACATGGCATCGTTGTTGAGAAAACCGGTCCGTATAACCTGCTGTTCCTGTTCAg                                                               >  1:838202/1‑67 (MQ=255)
    aaTACCTCGACGAACATGGCATCGTTGTTGAGAAAACCGGTCCGTATAACCTGCTGTTCCTGTTCAg                                                               >  1:114925/1‑67 (MQ=255)
                                agagAAAACCGGCCCTTATAACCTGCTGTTTCTCTTTAGTATTGGCATCGATAAAACCAAAGc                                       <  1:3073395/63‑1 (MQ=255)
                                agagAAAACCGGCCCTTATAACCTGCTGTTTCTCTTTAGTATTGGCATCGATAAAACCAAAGc                                       <  1:2927428/63‑1 (MQ=255)
                                agagAAAACCGGCCCTTATAACCTGCTGTTTCTCTTTAGTATTGGCATCGATAAAACCAAAGc                                       <  1:3201208/63‑1 (MQ=255)
                                agagAAAACCGGCCCTTATAACCTGCTGTTTCTCTTTAGTATTGGCATCGATAAAACCAAAGc                                       <  1:3422104/63‑1 (MQ=255)
                                agagAAAACCGGCCCTTATAACCTGCTGTTTCTCTTTAGTATTGGCATCGATAAAACCAAAGc                                       <  1:1040917/63‑1 (MQ=255)
                                agagAAAACCGGCCCTTATAACCTGCTGTTTCTCTTTAGTATTGGCATCGATAAAACCAAAGc                                       <  1:1097696/63‑1 (MQ=255)
                                agagAAAACCGGCCCTTATAACCTGCTGTTTCTCTTTAGTATTGGCATCGATAAAACCAAAGc                                       <  1:152510/63‑1 (MQ=255)
                                 gagaAAACCGGCCCTTATAACCTGCTGTTTCTCTTTAGTATTGGCATCGATAAAACCAAAGCAATggg                                 >  1:849883/1‑68 (MQ=255)
                                 gagaAAACCGGCCCTTATAACCTGCTGTTTCTCTTTAGTATTGGCATCGATAAAACCAAAGCAATggg                                 >  1:1655277/1‑68 (MQ=255)
                                 gagaAAACCGGCCCTTATAACCTGCTGTTTCTCTTTAGTATTGGCATCGATAAAACCAAAGCAATggg                                 >  1:65620/1‑68 (MQ=255)
                                 gagaAAACCGGCCCTTATAACCTGCTGTTTCTCTTTAGTATTGGCATCGATAAAACCAAAGCAATggg                                 >  1:1583036/1‑68 (MQ=255)
                                 gagaAAACCGGCCCTTATAACCTGCTGTTTCTCTTTAGTATTGGCATCGATAAAACCAAAGCAATggg                                 >  1:580706/1‑68 (MQ=255)
                                   gaAAACCGGCCCTTATAACCTGCTGTTTCTCTTTAGTATTGGCATCGATAAAACCAAAGCAATGGGat                               <  1:3028590/68‑1 (MQ=255)
                                                  taACCTGCTGTTCCTGTTCAGCATCGGTATCGATAAGACCaa                                          >  1:52375/1‑42 (MQ=12)
                                                  taACCTGCTGTTCCTGTTCAGCATCGGTATCGATAAGACCaa                                          >  1:3093203/1‑42 (MQ=12)
                                                  taACCTGCTGTTCCTGTTCAGCATCGGTATCGATAAGACCaa                                          >  1:1729141/1‑42 (MQ=12)
                                                  taACCTGCTGTTCCTGTTCAGCATCGGTATCGATAAGACCaa                                          >  1:1571990/1‑42 (MQ=12)
                                                  taACCTGCTGTTCCTGTTCAGCATCGGTATCGATAAGACCaa                                          >  1:1459043/1‑42 (MQ=12)
                                                  taACCTGCTGTTCCTGTTCAGCATCGGTATCGATAAGACCaa                                          >  1:1443292/1‑42 (MQ=12)
                                                  taACCTGCTGTTCCTGTTCAGCATCGGTATCGATAAGACCaa                                          >  1:1024320/1‑42 (MQ=12)
                                                  taACCTGCTGTTCCTGTTCAGCATCGGTATCGATAAGACCaa                                          >  1:119022/1‑42 (MQ=12)
                                                  taACCTGCTGTTCCTGTTCAGCATCGGTATCGATAAGACCaa                                          >  1:1169000/1‑42 (MQ=12)
                                                  taACCTGCTGTTCCTGTTCAGCATCGGTATCGATAAGACCaa                                          >  1:726003/1‑42 (MQ=12)
                                                  taACCTGCTGTTCCTGTTCAGCATCGGTATCGATAAGACCaa                                          >  1:1152916/1‑42 (MQ=12)
                                                  taACCTGCTGTTCCTGTTCAGCATCGGTATCGATAAGACCaa                                          >  1:2508104/1‑42 (MQ=12)
                                                  taACCTGCTGTTCCTGTTCAGCATCGGTATCGATAAGACCaa                                          >  1:1851622/1‑42 (MQ=12)
                                                  taACCTGCTGTTCCTGTTCAGCATCGGTATCGATAAGACCaa                                          >  1:2631109/1‑42 (MQ=12)
                                                  taACCTGCTGTTCCTGTTCAGCATCGGTATCGATAAGACCaa                                          >  1:2724801/1‑42 (MQ=12)
                                                  taACCTGCTGTTCCTGTTCAGCATCGGTATCGATAAGACCaa                                          >  1:2336588/1‑42 (MQ=12)
                                                  taACCTGCTGTTCCTGTTCAGCATCGGTATCGATAAGACCaa                                          >  1:2889633/1‑42 (MQ=12)
                                                  taACCTGCTGTTCCTGTTCAGCATCGGTATCGATAAGACCaa                                          >  1:1898445/1‑42 (MQ=12)
                                                  taACCTGCTGTTCCTGTTCAGCATCGGTATCGATAAGACCaa                                          >  1:2076834/1‑42 (MQ=12)
                                                           gTTTCTCTTTAGTATTGGCATCGATAAAACCAAAGCAATGGGATTATTGCGTGGGTTGACGGAATTCaa      >  1:457121/1‑69 (MQ=255)
                                                              tctcTTTAGTATTGGCATCGATAAAACCAAAGCAATggg                                 <  1:2389604/39‑1 (MQ=255)
                                                              tctcTTTAGTATTGGCATCGATAAAACCAAAGCAATggg                                 <  1:2519936/39‑1 (MQ=255)
                                                              tctcTTTAGTATTGGCATCGATAAAACCAAAGCAATggg                                 <  1:880346/39‑1 (MQ=255)
                                                              tctcTTTAGTATTGGCATCGATAAAACCAAAGCAATggg                                 <  1:3061172/39‑1 (MQ=255)
                                                              tctcTTTAGTATTGGCATCGATAAAACCAAAGCAATggg                                 <  1:844963/39‑1 (MQ=255)
                                                              tctcTTTAGTATTGGCATCGATAAAACCAAAGCAATggg                                 <  1:720310/39‑1 (MQ=255)
                                                              tctcTTTAGTATTGGCATCGATAAAACCAAAGCAATggg                                 <  1:208881/39‑1 (MQ=255)
                                                              tctcTTTAGTATTGGCATCGATAAAACCAAAGCAATggg                                 <  1:1606807/39‑1 (MQ=255)
                                                              tctcTTTAGTATTGGCATCGATAAAACCAAAGCAATggg                                 <  1:1350954/39‑1 (MQ=255)
                                                              tctcTTTAGTATTGGCATCGATAAAACCAAAGCAATggg                                 <  1:2911369/39‑1 (MQ=255)
                                                              tctcTTTAGTATTGGCATCGATAAAACCAAAGCAATggg                                 <  1:3227026/39‑1 (MQ=255)
                                                              tctcTTTAGTATTGGCATCGATAAAACCAAAGCAATggg                                 <  1:3401470/39‑1 (MQ=255)
                                                              tctcTTTAGTATTGGCATCGATAAAACCAAAGCAATggg                                 <  1:1767828/39‑1 (MQ=255)
                                                              tctcTTTAGTATTGGCATCGATAAAACCAAAGCAATggg                                 <  1:375301/39‑1 (MQ=255)
                                                              tctcTTTAGTATTGGCATCGATAAAACCAAAGCAATggg                                 <  1:1321500/39‑1 (MQ=255)
                                                              tctcTTTAGTATTGGCATCGATAAAACCAAAGCAATGGGATTATTGCGTg                      >  1:2413438/1‑50 (MQ=255)
                                                              tctcTTTAGTATTGGCATCGATAAAACCAAAGCAATGGGATTATTGCGTg                      >  1:680247/1‑50 (MQ=255)
                                                              tctcTTTAGTATTGGCATCGATAAAACCAAAGCAATGGGATTATTGCGTg                      >  1:825791/1‑50 (MQ=255)
                                                              tctcTTTAGTATTGGCATCGATAAAACCAAAGCAATGGGATTATTGCGTg                      >  1:3386938/1‑50 (MQ=255)
                                                               ctctTTAGTATTGGCATCGATAAAACCAAAGCAATGGGATTATTGCGTGGGTTGACGGAATTCAAACGc  >  1:3055574/1‑69 (MQ=255)
                                                               ctctTTAGTATTGGCATCGATAAAACCAAAGCAATGGGATTATTGCGTGGGTTGACGGAATTCAAACGc  >  1:2423757/1‑69 (MQ=255)
                                                               ctctTTAGTATTGGCATCGATAAAACCAAAGCAATGGGATTATTGCGTGGGTTGACGGAATTCAAACGc  >  1:59188/1‑69 (MQ=255)
                                                               ctctTTAGTATTGGCATCGATAAAACCAAAGCAATGGGATTATTGCGTGGGTTGACGGAATTCAAACGc  >  1:2613932/1‑69 (MQ=255)
                                                               ctctTTAGTATTGGCATCGATAAAACCAAAGCAATGGGATTATTGCGTGGGTTGACGGAATTCAAACGc  >  1:2884850/1‑69 (MQ=255)
                                                               ctctTTAGTATTGGCATCGATAAAACCAAAGCAATGGGATTATTGCGTCGGTTGACGGAATTCAAACGc  >  1:2652992/1‑69 (MQ=255)
                                                                    |                                                               
GCAAAATTCCTCGACGAACGTGGGATCGTAGTAGAGAAAACCGGCCCTTATAACCTGCTGTTTCTCTTTAGTATTGGCATCGATAAAACCAAAGCAATGGGATTATTGCGTGGGTTGACGGAATTCAAACGC  >  minE/190046‑190177

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: