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

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*minE190,1080TC80.6% 111.7 / 19.7 87F535F (TTT→TTCldcClysine decarboxylase 2, constitutive
Reads supporting (aligned to +/- strand):  ref base T (7/10);  new base C (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 = 9.15e-01

Only 100 of 113 total aligned reads displayed.

AGCAAAATTCCTCGACGAACGTGGGATCGTAGTAGAGAAAACCGGCCCTTATAACCTGCTGTTTCTCTTTAGTATTGGCATCGATAAAACCAAAGCAATGGGATTATTGCGTGGGTTGACGGAATTCAA  >  minE/190045‑190173
                                                               |                                                                 
aGCAAAATTCCTCGACGAACGTGGGATCGTAGTAGAGAAAACCGGCCCTTATAACCTGCTGTTtctc                                                                >  1:872027/1‑67 (MQ=255)
aGCAAAATTCCTCGACGAACGTGGGATCGTAGTAGAGAAAACCGGCCCTTATAACCTGCTGTTtctc                                                                >  1:1747043/1‑67 (MQ=255)
aGCAAAATTCCTCGACGAACGTGGGATCGTAGTAGAGAAAACCGGCCCTTATAACCTGCTGTTtctc                                                                >  1:431183/1‑67 (MQ=255)
aGCAAAATTCCTCGACGAACGTGGGATCGTAGTAGAGAAAACCGGCCCTTATAACCTGCTGTTtct                                                                 >  1:3346748/1‑66 (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:2072512/69‑3 (MQ=255)
     aaTACCTCGACGAACATGGCATCGTTGTTGAGAAAACCGGTCCGTATAACCTGCTGTTCCTGTTCAgca                                                         <  1:1607158/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:3077589/69‑3 (MQ=255)
     aaTACCTCGACGAACATGGCATCGTTGTTGAGAAAACCGGTCCGTATAACCTGCTGTTCCTGTTCAgca                                                         <  1:953241/69‑3 (MQ=255)
     aaTACCTCGACGAACATGGCATCGTTGTTGAGAAAACCGGTCCGTATAACCTGCTGTTCCTGTTCAgca                                                         <  1:2151660/69‑3 (MQ=255)
     aaTACCTCGACGAACATGGCATCGTTGTTGAGAAAACCGGTCCGTATAACCTGCTGTTCCTGTTCAgca                                                         <  1:2926966/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:2668939/69‑3 (MQ=255)
     aaTACCTCGACGAACATGGCATCGTTGTTGAGAAAACCGGTCCGTATAACCTGCTGTTCCTGTTCAgca                                                         <  1:2663649/69‑3 (MQ=255)
     aaTACCTCGACGAACATGGCATCGTTGTTGAGAAAACCGGTCCGTATAACCTGCTGTTCCTGTTCAgca                                                         <  1:3314474/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:1419785/69‑3 (MQ=255)
     aaTACCTCGACGAACATGGCATCGTTGTTGAGAAAACCGGTCCGTATAACCTGCTGTTCCTGTTCAgca                                                         <  1:592151/69‑3 (MQ=255)
     aaTACCTCGACGAACATGGCATCGTTGTTGAGAAAACCGGTCCGTATAACCTGCTGTTCCTGTTCAg                                                           >  1:3243803/1‑67 (MQ=255)
     aaTACCTCGACGAACATGGCATCGTTGTTGAGAAAACCGGTCCGTATAACCTGCTGTTCCTGTTCAg                                                           >  1:2485818/1‑67 (MQ=255)
     aaTACCTCGACGAACATGGCATCGTTGTTGAGAAAACCGGTCCGTATAACCTGCTGTTCCTGTTCAg                                                           >  1:114925/1‑67 (MQ=255)
     aaTACCTCGACGAACATGGCATCGTTGTTGAGAAAACCGGTCCGTATAACCTGCTGTTCCTGTTCAg                                                           >  1:838202/1‑67 (MQ=255)
     aaTACCTCGACGAACATGGCATCGTTGTTGAGAAAACCGGTCCGTATAACCTGCTGTTCCTGTTCAg                                                           >  1:1278889/1‑67 (MQ=255)
     aaTACCTCGACGAACATGGCATCGTTGTTGAGAAAACCGGTCCGTATAACCTGCTGTTCCTGTTCAg                                                           >  1:2225287/1‑67 (MQ=255)
     aaTACCTCGACGAACATGGCATCGTTGTTGAGAAAACCGGTCCGTATAACCTGCTGTTCCTGTTCAg                                                           >  1:1307497/1‑67 (MQ=255)
     aaTACCTCGACGAACATGGCATCGTTGTTGAGAAAACCGGTCCGTATAACCTGCTGTTCCTGTTCAg                                                           >  1:209702/1‑67 (MQ=255)
     aaTACCTCGACGAACATGGCATCGTTGTTGAGAAAACCGGTCCGTATAACCTGCTGTTCCTGTTCAg                                                           >  1:2017844/1‑67 (MQ=255)
     aaTACCTCGACGAACATGGCATCGTTGTTGAGAAAACCGGTCCGTATAACCTGCTGTTCCTGTTCAg                                                           >  1:397673/1‑67 (MQ=255)
     aaTACCTCGACGAACATGGCATCGTTGTTGAGAAAACCGGTCCGTATAACCTGCTGTTCCTGTTCAg                                                           >  1:3183370/1‑67 (MQ=255)
     aaTACCTCGACGAACATGGCATCGTTGTTGAGAAAACCGGTCCGTATAACCTGCTGTTCCTGTTCAg                                                           >  1:1528336/1‑67 (MQ=255)
     aaTACCTCGACGAACATGGCATCGTTGTTGAGAAAACCGGTCCGTATAACCTGCTGTTCCTGTTCAg                                                           >  1:1791176/1‑67 (MQ=255)
     aaTACCTCGACGAACATGGCATCGTTGTTGAGAAAACCGGTCCGTATAACCTGCTGTTCCTGTTCAg                                                           >  1:3196280/1‑67 (MQ=255)
     aaTACCTCGACGAACATGGCATCGTTGTTGAGAAAACCGGTCCGTATAACCTGCTGTTCCTGTTCAg                                                           >  1:3317817/1‑67 (MQ=255)
                                 agagAAAACCGGCCCTTATAACCTGCTGTTTCTCTTTAGTATTGGCATCGATAAAACCAAAGc                                   <  1:2927428/63‑1 (MQ=255)
                                 agagAAAACCGGCCCTTATAACCTGCTGTTTCTCTTTAGTATTGGCATCGATAAAACCAAAGc                                   <  1:3073395/63‑1 (MQ=255)
                                 agagAAAACCGGCCCTTATAACCTGCTGTTTCTCTTTAGTATTGGCATCGATAAAACCAAAGc                                   <  1:3422104/63‑1 (MQ=255)
                                 agagAAAACCGGCCCTTATAACCTGCTGTTTCTCTTTAGTATTGGCATCGATAAAACCAAAGc                                   <  1:3201208/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:580706/1‑68 (MQ=255)
                                  gagaAAACCGGCCCTTATAACCTGCTGTTTCTCTTTAGTATTGGCATCGATAAAACCAAAGCAATggg                             >  1:1655277/1‑68 (MQ=255)
                                  gagaAAACCGGCCCTTATAACCTGCTGTTTCTCTTTAGTATTGGCATCGATAAAACCAAAGCAATggg                             >  1:1583036/1‑68 (MQ=255)
                                  gagaAAACCGGCCCTTATAACCTGCTGTTTCTCTTTAGTATTGGCATCGATAAAACCAAAGCAATggg                             >  1:65620/1‑68 (MQ=255)
                                  gagaAAACCGGCCCTTATAACCTGCTGTTTCTCTTTAGTATTGGCATCGATAAAACCAAAGCAATggg                             >  1:849883/1‑68 (MQ=255)
                                  gagaAAACCGGCCCTTATAACCTGCTGTTTCTCTTTAGTATTGGCATCGATAAAACCAAAGCAATggg                             >  1:2499447/1‑68 (MQ=255)
                                    gaAAACCGGCCCTTATAACCTGCTGTTTCTCTTTAGTATTGGCATCGATAAAACCAAAGCAATGGGat                           <  1:3028590/68‑1 (MQ=255)
                                                   taACCTGCTGTTCCTGTTCAGCATCGGTATCGATAAGACCaa                                      >  1:1169000/1‑42 (MQ=12)
                                                   taACCTGCTGTTCCTGTTCAGCATCGGTATCGATAAGACCaa                                      >  1:119022/1‑42 (MQ=12)
                                                   taACCTGCTGTTCCTGTTCAGCATCGGTATCGATAAGACCaa                                      >  1:726003/1‑42 (MQ=12)
                                                   taACCTGCTGTTCCTGTTCAGCATCGGTATCGATAAGACCaa                                      >  1:1571990/1‑42 (MQ=12)
                                                   taACCTGCTGTTCCTGTTCAGCATCGGTATCGATAAGACCaa                                      >  1:52375/1‑42 (MQ=12)
                                                   taACCTGCTGTTCCTGTTCAGCATCGGTATCGATAAGACCaa                                      >  1:1443292/1‑42 (MQ=12)
                                                   taACCTGCTGTTCCTGTTCAGCATCGGTATCGATAAGACCaa                                      >  1:1459043/1‑42 (MQ=12)
                                                   taACCTGCTGTTCCTGTTCAGCATCGGTATCGATAAGACCaa                                      >  1:1524466/1‑42 (MQ=12)
                                                   taACCTGCTGTTCCTGTTCAGCATCGGTATCGATAAGACCaa                                      >  1:1152916/1‑42 (MQ=12)
                                                   taACCTGCTGTTCCTGTTCAGCATCGGTATCGATAAGACCaa                                      >  1:3330938/1‑42 (MQ=12)
                                                   taACCTGCTGTTCCTGTTCAGCATCGGTATCGATAAGACCaa                                      >  1:2508104/1‑42 (MQ=12)
                                                   taACCTGCTGTTCCTGTTCAGCATCGGTATCGATAAGACCaa                                      >  1:1024320/1‑42 (MQ=12)
                                                   taACCTGCTGTTCCTGTTCAGCATCGGTATCGATAAGACCaa                                      >  1:1743818/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:2207819/1‑42 (MQ=12)
                                                   taACCTGCTGTTCCTGTTCAGCATCGGTATCGATAAGACCaa                                      >  1:2947134/1‑42 (MQ=12)
                                                   taACCTGCTGTTCCTGTTCAGCATCGGTATCGATAAGACCaa                                      >  1:2076834/1‑42 (MQ=12)
                                                   taACCTGCTGTTCCTGTTCAGCATCGGTATCGATAAGACCaa                                      >  1:3093203/1‑42 (MQ=12)
                                                   taACCTGCTGTTCCTGTTCAGCATCGGTATCGATAAGACCaa                                      >  1:1851622/1‑42 (MQ=12)
                                                            gTTTCTCTTTAGTATTGGCATCGATAAAACCAAAGCAATGGGATTATTGCGTGGGTTGACGGAATTCaa  >  1:457121/1‑69 (MQ=255)
                                                               tctcTTTAGTATTGGCATCGATAAAACCAAAGCAATggg                             <  1:2911369/39‑1 (MQ=255)
                                                               tctcTTTAGTATTGGCATCGATAAAACCAAAGCAATggg                             <  1:3227026/39‑1 (MQ=255)
                                                               tctcTTTAGTATTGGCATCGATAAAACCAAAGCAATggg                             <  1:880346/39‑1 (MQ=255)
                                                               tctcTTTAGTATTGGCATCGATAAAACCAAAGCAATggg                             <  1:1606807/39‑1 (MQ=255)
                                                               tctcTTTAGTATTGGCATCGATAAAACCAAAGCAATggg                             <  1:844963/39‑1 (MQ=255)
                                                               tctcTTTAGTATTGGCATCGATAAAACCAAAGCAATggg                             <  1:2389604/39‑1 (MQ=255)
                                                               tctcTTTAGTATTGGCATCGATAAAACCAAAGCAATggg                             <  1:1767828/39‑1 (MQ=255)
                                                               tctcTTTAGTATTGGCATCGATAAAACCAAAGCAATggg                             <  1:720310/39‑1 (MQ=255)
                                                               tctcTTTAGTATTGGCATCGATAAAACCAAAGCAATggg                             <  1:3401470/39‑1 (MQ=255)
                                                               tctcTTTAGTATTGGCATCGATAAAACCAAAGCAATggg                             <  1:1321500/39‑1 (MQ=255)
                                                               tctcTTTAGTATTGGCATCGATAAAACCAAAGCAATggg                             <  1:1350954/39‑1 (MQ=255)
                                                               tctcTTTAGTATTGGCATCGATAAAACCAAAGCAATggg                             <  1:208881/39‑1 (MQ=255)
                                                               tctcTTTAGTATTGGCATCGATAAAACCAAAGCAATggg                             <  1:1984854/39‑1 (MQ=255)
                                                               tctcTTTAGTATTGGCATCGATAAAACCAAAGCAATggg                             <  1:3061172/39‑1 (MQ=255)
                                                               tctcTTTAGTATTGGCATCGATAAAACCAAAGCAATggg                             <  1:375301/39‑1 (MQ=255)
                                                               tctcTTTAGTATTGGCATCGATAAAACCAAAGCAATggg                             <  1:2519936/39‑1 (MQ=255)
                                                               tctcTTTAGTATTGGCATCGATAAAACCAAAGCAATGGGATTATTGCGTg                  >  1:680247/1‑50 (MQ=255)
                                                               tctcTTTAGTATTGGCATCGATAAAACCAAAGCAATGGGATTATTGCGTg                  >  1:2413438/1‑50 (MQ=255)
                                                               tctcTTTAGTATTGGCATCGATAAAACCAAAGCAATGGGATTATTGCGTg                  >  1:3386938/1‑50 (MQ=255)
                                                               |                                                                 
AGCAAAATTCCTCGACGAACGTGGGATCGTAGTAGAGAAAACCGGCCCTTATAACCTGCTGTTTCTCTTTAGTATTGGCATCGATAAAACCAAAGCAATGGGATTATTGCGTGGGTTGACGGAATTCAA  >  minE/190045‑190173

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: