Predicted mutation
evidence seq id position mutation freq annotation gene description
RA minE 190,111 C→G 80.5% L536L (CTC→CTG ldcC → lysine decarboxylase 2, constitutive

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*minE190,1110CG80.5% 131.0 / 32.0 87L536L (CTC→CTGldcClysine decarboxylase 2, constitutive
Reads supporting (aligned to +/- strand):  ref base C (7/10);  new base G (45/25);  total (52/35)
Fisher's exact test for biased strand distribution p-value = 1.02e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 8.70e-01

Only 100 of 118 total aligned reads displayed.

AGCAAAATTCCTCGACGAACGTGGGATCGTAGTAGAGAAAACCGGCCCTTATAACCTGCTGTTTCTCTTTAGTATTGGCATCGATAAAACCAAAGCAATGGGATTATTGCGTGGGTTGACGGAATTCAAACGC  >  minE/190045‑190177
                                                                  |                                                                  
aGCAAAATTCCTCGACGAACGTGGGATCGTAGTAGAGAAAACCGGCCCTTATAACCTGCTGTTtctc                                                                    >  1:872027/1‑67 (MQ=255)
aGCAAAATTCCTCGACGAACGTGGGATCGTAGTAGAGAAAACCGGCCCTTATAACCTGCTGTTtctc                                                                    >  1:1747043/1‑67 (MQ=255)
aGCAAAATTCCTCGACGAACGTGGGATCGTAGTAGAGAAAACCGGCCCTTATAACCTGCTGTTtctc                                                                    >  1:2442649/1‑67 (MQ=255)
aGCAAAATTCCTCGACGAACGTGGGATCGTAGTAGAGAAAACCGGCCCTTATAACCTGCTGTTtctc                                                                    >  1:1491608/1‑67 (MQ=255)
aGCAAAATTCCTCGACGAACGTGGGATCGTAGTAGAGAAAACCGGCCCTTATAACCTGCTGTTtctc                                                                    >  1:431183/1‑67 (MQ=255)
 gCAAAATTCCTCGACGAACGTGGGATCGTAGTAGAGAAAACCGGCCCTTATAACCTGCTGTTTCTCttt                                                                 <  1:2359602/69‑1 (MQ=255)
    aaaTTCCTCGACGAACGTGGGATCGTAGTAGAGAAAACCGGCCCTTATAACCTGCTGTTTCTCTTTAGt                                                              <  1:786128/69‑1 (MQ=255)
     aaTACCTCGACGCACATGGCATCGTTGTTGAGAAAACCGGTCCGTATAACCTGCTGTTCCTGTTCAg                                                               >  1:3303229/1‑67 (MQ=255)
     aaTACCTCGACGAACATTGCATCGTTGTTGAGAAAACCGGTCCGTATAACCTGCTGTTCCTGTTCAgca                                                             <  1:206573/69‑3 (MQ=255)
     aaTACCTCGACGAACATGGCATCGTTGTTGAGAAAACCGGTCCGTATAACCTGCTGTTCCTGTTCAgca                                                             <  1:3314474/69‑3 (MQ=255)
     aaTACCTCGACGAACATGGCATCGTTGTTGAGAAAACCGGTCCGTATAACCTGCTGTTCCTGTTCAgca                                                             <  1:3224580/69‑3 (MQ=255)
     aaTACCTCGACGAACATGGCATCGTTGTTGAGAAAACCGGTCCGTATAACCTGCTGTTCCTGTTCAgca                                                             <  1:1846561/69‑3 (MQ=255)
     aaTACCTCGACGAACATGGCATCGTTGTTGAGAAAACCGGTCCGTATAACCTGCTGTTCCTGTTCAgca                                                             <  1:1888849/69‑3 (MQ=255)
     aaTACCTCGACGAACATGGCATCGTTGTTGAGAAAACCGGTCCGTATAACCTGCTGTTCCTGTTCAgca                                                             <  1:1904420/69‑3 (MQ=255)
     aaTACCTCGACGAACATGGCATCGTTGTTGAGAAAACCGGTCCGTATAACCTGCTGTTCCTGTTCAgca                                                             <  1:2151660/69‑3 (MQ=255)
     aaTACCTCGACGAACATGGCATCGTTGTTGAGAAAACCGGTCCGTATAACCTGCTGTTCCTGTTCAgca                                                             <  1:2403666/69‑3 (MQ=255)
     aaTACCTCGACGAACATGGCATCGTTGTTGAGAAAACCGGTCCGTATAACCTGCTGTTCCTGTTCAgca                                                             <  1:2926966/69‑3 (MQ=255)
     aaTACCTCGACGAACATGGCATCGTTGTTGAGAAAACCGGTCCGTATAACCTGCTGTTCCTGTTCAgca                                                             <  1:2854087/69‑3 (MQ=255)
     aaTACCTCGACGAACATGGCATCGTTGTTGAGAAAACCGGTCCGTATAACCTGCTGTTCCTGTTCAgca                                                             <  1:2680412/69‑3 (MQ=255)
     aaTACCTCGACGAACATGGCATCGTTGTTGAGAAAACCGGTCCGTATAACCTGCTGTTCCTGTTCAgca                                                             <  1:2668939/69‑3 (MQ=255)
     aaTACCTCGACGAACATGGCATCGTTGTTGAGAAAACCGGTCCGTATAACCTGCTGTTCCTGTTCAgca                                                             <  1:2663649/69‑3 (MQ=255)
     aaTACCTCGACGAACATGGCATCGTTGTTGAGAAAACCGGTCCGTATAACCTGCTGTTCCTGTTCAgca                                                             <  1:953241/69‑3 (MQ=255)
     aaTACCTCGACGAACATGGCATCGTTGTTGAGAAAACCGGTCCGTATAACCTGCTGTTCCTGTTCAgca                                                             <  1:3352934/69‑3 (MQ=255)
     aaTACCTCGACGAACATGGCATCGTTGTTGAGAAAACCGGTCCGTATAACCTGCTGTTCCTGTTCAgca                                                             <  1:870975/69‑3 (MQ=255)
     aaTACCTCGACGAACATGGCATCGTTGTTGAGAAAACCGGTCCGTATAACCTGCTGTTCCTGTTCAgca                                                             <  1:1233416/69‑3 (MQ=255)
     aaTACCTCGACGAACATGGCATCGTTGTTGAGAAAACCGGTCCGTATAACCTGCTGTTCCTGTTCAgca                                                             <  1:1468180/69‑3 (MQ=255)
     aaTACCTCGACGAACATGGCATCGTTGTTGAGAAAACCGGTCCGTATAACCTGCTGTTCCTGTTCAgca                                                             <  1:616781/69‑3 (MQ=255)
     aaTACCTCGACGAACATGGCATCGTTGTTGAGAAAACCGGTCCGTATAACCTGCTGTTCCTGTTCAgca                                                             <  1:592151/69‑3 (MQ=255)
     aaTACCTCGACGAACATGGCATCGTTGTTGAGAAAACCGGTCCGTATAACCTGCTGTTCCTGTTCAgca                                                             <  1:1607158/69‑3 (MQ=255)
     aaTACCTCGACGAACATGGCATCGTTGTTGAGAAAACCGGTCCGTATAACCTGCTGTTCCTGTTCAg                                                               >  1:209702/1‑67 (MQ=255)
     aaTACCTCGACGAACATGGCATCGTTGTTGAGAAAACCGGTCCGTATAACCTGCTGTTCCTGTTCAg                                                               >  1:114925/1‑67 (MQ=255)
     aaTACCTCGACGAACATGGCATCGTTGTTGAGAAAACCGGTCCGTATAACCTGCTGTTCCTGTTCAg                                                               >  1:838202/1‑67 (MQ=255)
     aaTACCTCGACGAACATGGCATCGTTGTTGAGAAAACCGGTCCGTATAACCTGCTGTTCCTGTTCAg                                                               >  1:2485818/1‑67 (MQ=255)
     aaTACCTCGACGAACATGGCATCGTTGTTGAGAAAACCGGTCCGTATAACCTGCTGTTCCTGTTCAg                                                               >  1:1278889/1‑67 (MQ=255)
     aaTACCTCGACGAACATGGCATCGTTGTTGAGAAAACCGGTCCGTATAACCTGCTGTTCCTGTTCAg                                                               >  1:1307497/1‑67 (MQ=255)
     aaTACCTCGACGAACATGGCATCGTTGTTGAGAAAACCGGTCCGTATAACCTGCTGTTCCTGTTCAg                                                               >  1:3028215/1‑67 (MQ=255)
     aaTACCTCGACGAACATGGCATCGTTGTTGAGAAAACCGGTCCGTATAACCTGCTGTTCCTGTTCAg                                                               >  1:2225287/1‑67 (MQ=255)
     aaTACCTCGACGAACATGGCATCGTTGTTGAGAAAACCGGTCCGTATAACCTGCTGTTCCTGTTCAg                                                               >  1:211024/1‑67 (MQ=255)
     aaTACCTCGACGAACATGGCATCGTTGTTGAGAAAACCGGTCCGTATAACCTGCTGTTCCTGTTCAg                                                               >  1:3183370/1‑67 (MQ=255)
     aaTACCTCGACGAACATGGCATCGTTGTTGAGAAAACCGGTCCGTATAACCTGCTGTTCCTGTTCAg                                                               >  1:566402/1‑67 (MQ=255)
     aaTACCTCGACGAACATGGCATCGTTGTTGAGAAAACCGGTCCGTATAACCTGCTGTTCCTGTTCAg                                                               >  1:2017844/1‑67 (MQ=255)
     aaTACCTCGACGAACATGGCATCGTTGTTGAGAAAACCGGTCCGTATAACCTGCTGTTCCTGTTCAg                                                               >  1:3196280/1‑67 (MQ=255)
     aaTACCTCGACGAACATGGCATCGTTGTTGAGAAAACCGGTCCGTATAACCTGCTGTTCCTGTTCAg                                                               >  1:397673/1‑67 (MQ=255)
     aaTACCTCGACGAACATGGCATCGTTGTTGAGAAAACCGGTCCGTATAACCTGCTGTTCCTGTTCAg                                                               >  1:1528336/1‑67 (MQ=255)
     aaTACCTCGACGAACATGGCATCGTTGTTGAGAAAACCGGTCCGTATAACCTGCTGTTCCTGTTCAg                                                               >  1:1791176/1‑67 (MQ=255)
     aaTACCTCGACGAACATGGCATCGTTGTTGAGAAAACCGGTCCGTATAACCTGCTGTTCCTGTTCAg                                                               >  1:3317817/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:1040917/63‑1 (MQ=255)
                                 agagAAAACCGGCCCTTATAACCTGCTGTTTCTCTTTAGTATTGGCATCGATAAAACCAAAGc                                       <  1:3422104/63‑1 (MQ=255)
                                 agagAAAACCGGCCCTTATAACCTGCTGTTTCTCTTTAGTATTGGCATCGATAAAACCAAAGc                                       <  1:1097696/63‑1 (MQ=255)
                                  gagaAAACCGGCCCTTATAACCTGCTGTTTCTCTTTAGTATTGGCATCGATAAAACCAAAGCAATggg                                 >  1:1655277/1‑68 (MQ=255)
                                  gagaAAACCGGCCCTTATAACCTGCTGTTTCTCTTTAGTATTGGCATCGATAAAACCAAAGCAATggg                                 >  1:849883/1‑68 (MQ=255)
                                  gagaAAACCGGCCCTTATAACCTGCTGTTTCTCTTTAGTATTGGCATCGATAAAACCAAAGCAATggg                                 >  1:2499447/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:1743818/1‑42 (MQ=12)
                                                   taACCTGCTGTTCCTGTTCAGCATCGGTATCGATAAGACCaa                                          >  1:2631109/1‑42 (MQ=12)
                                                   taACCTGCTGTTCCTGTTCAGCATCGGTATCGATAAGACCaa                                          >  1:3330938/1‑42 (MQ=12)
                                                   taACCTGCTGTTCCTGTTCAGCATCGGTATCGATAAGACCaa                                          >  1:1571990/1‑42 (MQ=12)
                                                   taACCTGCTGTTCCTGTTCAGCATCGGTATCGATAAGACCaa                                          >  1:1524466/1‑42 (MQ=12)
                                                   taACCTGCTGTTCCTGTTCAGCATCGGTATCGATAAGACCaa                                          >  1:1459043/1‑42 (MQ=12)
                                                   taACCTGCTGTTCCTGTTCAGCATCGGTATCGATAAGACCaa                                          >  1:1443292/1‑42 (MQ=12)
                                                   taACCTGCTGTTCCTGTTCAGCATCGGTATCGATAAGACCaa                                          >  1:726003/1‑42 (MQ=12)
                                                   taACCTGCTGTTCCTGTTCAGCATCGGTATCGATAAGACCaa                                          >  1:1241732/1‑42 (MQ=12)
                                                   taACCTGCTGTTCCTGTTCAGCATCGGTATCGATAAGACCaa                                          >  1:1169000/1‑42 (MQ=12)
                                                   taACCTGCTGTTCCTGTTCAGCATCGGTATCGATAAGACCaa                                          >  1:1152916/1‑42 (MQ=12)
                                                   taACCTGCTGTTCCTGTTCAGCATCGGTATCGATAAGACCaa                                          >  1:1024320/1‑42 (MQ=12)
                                                   taACCTGCTGTTCCTGTTCAGCATCGGTATCGATAAGACCaa                                          >  1:2336588/1‑42 (MQ=12)
                                                   taACCTGCTGTTCCTGTTCAGCATCGGTATCGATAAGACCaa                                          >  1:2724801/1‑42 (MQ=12)
                                                   taACCTGCTGTTCCTGTTCAGCATCGGTATCGATAAGACCaa                                          >  1:2889633/1‑42 (MQ=12)
                                                   taACCTGCTGTTCCTGTTCAGCATCGGTATCGATAAGACCaa                                          >  1:2508104/1‑42 (MQ=12)
                                                   taACCTGCTGTTCCTGTTCAGCATCGGTATCGATAAGACCaa                                          >  1:1898445/1‑42 (MQ=12)
                                                   taACCTGCTGTTCCTGTTCAGCATCGGTATCGATAAGACCaa                                          >  1:3093203/1‑42 (MQ=12)
                                                   taACCTGCTGTTCCTGTTCAGCATCGGTATCGATAAGACCaa                                          >  1:2076834/1‑42 (MQ=12)
                                                            gTTTCTCTTTAGTATTGGCATCGATAAAACCAAAGCAATGGGATTATTGCGTGGGTTGACGGAATTCaa      >  1:457121/1‑69 (MQ=255)
                                                               tctcTTTAGTATTGGCATCGATAAAACCAAAGCAATggg                                 <  1:1767828/39‑1 (MQ=255)
                                                               tctcTTTAGTATTGGCATCGATAAAACCAAAGCAATggg                                 <  1:1606807/39‑1 (MQ=255)
                                                               tctcTTTAGTATTGGCATCGATAAAACCAAAGCAATggg                                 <  1:3227026/39‑1 (MQ=255)
                                                               tctcTTTAGTATTGGCATCGATAAAACCAAAGCAATggg                                 <  1:2911369/39‑1 (MQ=255)
                                                               tctcTTTAGTATTGGCATCGATAAAACCAAAGCAATggg                                 <  1:720310/39‑1 (MQ=255)
                                                               tctcTTTAGTATTGGCATCGATAAAACCAAAGCAATggg                                 <  1:1984854/39‑1 (MQ=255)
                                                               tctcTTTAGTATTGGCATCGATAAAACCAAAGCAATggg                                 <  1:1321500/39‑1 (MQ=255)
                                                               tctcTTTAGTATTGGCATCGATAAAACCAAAGCAATggg                                 <  1:1350954/39‑1 (MQ=255)
                                                               tctcTTTAGTATTGGCATCGATAAAACCAAAGCAATggg                                 <  1:2519936/39‑1 (MQ=255)
                                                               tctcTTTAGTATTGGCATCGATAAAACCAAAGCAATggg                                 <  1:2389604/39‑1 (MQ=255)
                                                               tctcTTTAGTATTGGCATCGATAAAACCAAAGCAATggg                                 <  1:880346/39‑1 (MQ=255)
                                                               tctcTTTAGTATTGGCATCGATAAAACCAAAGCAATggg                                 <  1:375301/39‑1 (MQ=255)
                                                               tctcTTTAGTATTGGCATCGATAAAACCAAAGCAATggg                                 <  1:3061172/39‑1 (MQ=255)
                                                               tctcTTTAGTATTGGCATCGATAAAACCAAAGCAATggg                                 <  1:3401470/39‑1 (MQ=255)
                                                               tctcTTTAGTATTGGCATCGATAAAACCAAAGCAATggg                                 <  1:208881/39‑1 (MQ=255)
                                                               tctcTTTAGTATTGGCATCGATAAAACCAAAGCAATGGGATTATTGCGTg                      >  1:2413438/1‑50 (MQ=255)
                                                               tctcTTTAGTATTGGCATCGATAAAACCAAAGCAATGGGATTATTGCGTg                      >  1:825791/1‑50 (MQ=255)
                                                               tctcTTTAGTATTGGCATCGATAAAACCAAAGCAATGGGATTATTGCGTg                      >  1:680247/1‑50 (MQ=255)
                                                                ctctTTAGTATTGGCATCGATAAAACCAAAGCAATGGGATTATTGCGTGGGTTGACGGAATTCAAACGc  >  1:2884850/1‑69 (MQ=255)
                                                                ctctTTAGTATTGGCATCGATAAAACCAAAGCAATGGGATTATTGCGTGGGTTGACGGAATTCAAACGc  >  1:59188/1‑69 (MQ=255)
                                                                ctctTTAGTATTGGCATCGATAAAACCAAAGCAATGGGATTATTGCGTGGGTTGACGGAATTCAAACGc  >  1:3055574/1‑69 (MQ=255)
                                                                ctctTTAGTATTGGCATCGATAAAACCAAAGCAATGGGATTATTGCGTCGGTTGACGGAATTCAAACGc  >  1:2652992/1‑69 (MQ=255)
                                                                  |                                                                  
AGCAAAATTCCTCGACGAACGTGGGATCGTAGTAGAGAAAACCGGCCCTTATAACCTGCTGTTTCTCTTTAGTATTGGCATCGATAAAACCAAAGCAATGGGATTATTGCGTGGGTTGACGGAATTCAAACGC  >  minE/190045‑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: