Predicted mutation
evidence seq id position mutation freq annotation gene description
RA minE 190,093 T→G 53.7% P530P (CCT→CCG ldcC → lysine decarboxylase 2, constitutive

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*minE190,0930TG53.7% 28.4 / 82.5 82P530P (CCT→CCGldcClysine decarboxylase 2, constitutive
Reads supporting (aligned to +/- strand):  ref base T (13/25);  new base G (19/25);  total (32/50)
Fisher's exact test for biased strand distribution p-value = 4.98e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

ATCCCGGCGGCGCTGGTAGCAAAATTCCTCGACGAACGTGGGATCGTAGTAGAGAAAACCGGCCCTTATAACCTGCTGTTTCTCTTTAGTATTGGCATCGATAAAACCAAAGCAATGGGAT  >  minE/190028‑190148
                                                                 |                                                       
aTCCCGGCGGCGCTGGTAGCAAAATTCCTCGACGAACGTGGGATCGTAGTAGAGAAAACCGGCCCTtat                                                      <  1:51934/69‑1 (MQ=255)
 tCCCGGCGGCGCTGGTAGCAAAATTCCTCGACGAACGTGGGATCGTAGTAGAGAAAACCGGCCCTtata                                                     <  1:196227/69‑1 (MQ=255)
 tCCCGGCGGCGCTGGTAGCAAAATTCCTCGACGAACGTGGGATCGTAGTAGAGAAAACCGGCCCTtata                                                     <  1:618687/69‑1 (MQ=255)
  cccGGCGGCGCTGGTAGCAAAATTCCTCGACGAACGTGGGATCGTAGTAGAGAAAACCGGCCCTTATaa                                                    >  1:1257167/1‑69 (MQ=255)
     ggcggcGCTGGTAGCAAAATTCCTCGACGAACGTGTGATCGTAGTAGAGAAAACCGGCCCTtata                                                     <  1:894249/65‑1 (MQ=255)
     ggcggcGCTGGTAGCAAAATTCCTCGACGAACGTGGGATCGTAGTAGAGAAAACCGGCCCTtata                                                     <  1:106897/65‑1 (MQ=255)
     ggcggcGCTGGTAGCAAAATTCCTCGACGAACGTGGGATCGTAGTAGAGAAAACCGGCCCTtata                                                     <  1:2644101/65‑1 (MQ=255)
     ggcggcGCTGGTAGCAAAATTCCTCGACGAACGTGGGATCGTAGTAGAGAAAACCGGCCCTtata                                                     <  1:463346/65‑1 (MQ=255)
         gcgcTGGTAGCAAAATTCCTCGACGAACGTGGGATCGTAGTAGAGAAAACCGGCCCTTATAACctgctt                                             <  1:1834328/69‑2 (MQ=255)
         gcgcTGGTAGCAAAATTCCTCGACGAACGTGGGATCGTAGTAGAGAAAACCGGCCCTTATAACctgctg                                             <  1:149214/69‑1 (MQ=255)
         gcgcTGGTAGCAAAATTCCTCGACGAACGTGGGATCGTAGTAGAGAAAACCGGCCCTTATAACctgctg                                             <  1:2675369/69‑1 (MQ=255)
         gcgcTGGTAGCAAAATTCCTCGACGAACGTGGGATCGTAGTAGAGAAAACCGGCCCTTATAACctgctg                                             <  1:1821718/69‑1 (MQ=255)
         gcgcTGGTAGCAAAATTCCTCGACGAACGTGGGATCGTAGTAGAGAAAACCGGCCCTTATAACctgctg                                             <  1:621761/69‑1 (MQ=255)
          cgcTGGTAGCAAAATTCCTCGACGAACGTGGGATCGTAGTAGAGAAAACCGGCCCTTATAACctgctg                                             <  1:228627/68‑1 (MQ=255)
                 aGCAAAATTCCTCGACGAACGTGGGATCGTAGTAGAGAAAACCGGCCCTTATAACctgctg                                             <  1:2144436/61‑1 (MQ=255)
                 aGCAAAATTCCTCGACGAACGTGGGATCGTAGTAGAGAAAACCGGCCCTTATAACCTGCTGTTtctc                                       >  1:1491608/1‑67 (MQ=255)
                 aGCAAAATTCCTCGACGAACGTGGGATCGTAGTAGAGAAAACCGGCCCTTATAACCTGCTGTTtctc                                       >  1:2442649/1‑67 (MQ=255)
                 aGCAAAATTCCTCGACGAACGTGGGATCGTAGTAGAGAAAACCGGCCCTTATAACCTGCTGTTtctc                                       >  1:872027/1‑67 (MQ=255)
                 aGCAAAATTCCTCGACGAACGTGGGATCGTAGTAGAGAAAACCGGCCCTTATAACCTGCTGTTtctc                                       >  1:431183/1‑67 (MQ=255)
                 aGCAAAATTCCTCGACGAACGTGGGATCGTAGTAGAGAAAACCGGCCCTTATAACCTGCTGTTtctc                                       >  1:1747043/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:953241/69‑3 (MQ=255)
                      aaTACCTCGACGAACATGGCATCGTTGTTGAGAAAACCGGTCCGTATAACCTGCTGTTCCTGTTCAgca                                <  1:3224580/69‑3 (MQ=255)
                      aaTACCTCGACGAACATGGCATCGTTGTTGAGAAAACCGGTCCGTATAACCTGCTGTTCCTGTTCAgca                                <  1:3077589/69‑3 (MQ=255)
                      aaTACCTCGACGAACATGGCATCGTTGTTGAGAAAACCGGTCCGTATAACCTGCTGTTCCTGTTCAgca                                <  1:2926966/69‑3 (MQ=255)
                      aaTACCTCGACGAACATGGCATCGTTGTTGAGAAAACCGGTCCGTATAACCTGCTGTTCCTGTTCAgca                                <  1:2854087/69‑3 (MQ=255)
                      aaTACCTCGACGAACATGGCATCGTTGTTGAGAAAACCGGTCCGTATAACCTGCTGTTCCTGTTCAgca                                <  1:3314474/69‑3 (MQ=255)
                      aaTACCTCGACGAACATGGCATCGTTGTTGAGAAAACCGGTCCGTATAACCTGCTGTTCCTGTTCAgca                                <  1:2842915/69‑3 (MQ=255)
                      aaTACCTCGACGAACATGGCATCGTTGTTGAGAAAACCGGTCCGTATAACCTGCTGTTCCTGTTCAgca                                <  1:2680412/69‑3 (MQ=255)
                      aaTACCTCGACGAACATGGCATCGTTGTTGAGAAAACCGGTCCGTATAACCTGCTGTTCCTGTTCAgca                                <  1:3352934/69‑3 (MQ=255)
                      aaTACCTCGACGAACATGGCATCGTTGTTGAGAAAACCGGTCCGTATAACCTGCTGTTCCTGTTCAgca                                <  1:870975/69‑3 (MQ=255)
                      aaTACCTCGACGAACATGGCATCGTTGTTGAGAAAACCGGTCCGTATAACCTGCTGTTCCTGTTCAgca                                <  1:592151/69‑3 (MQ=255)
                      aaTACCTCGACGAACATGGCATCGTTGTTGAGAAAACCGGTCCGTATAACCTGCTGTTCCTGTTCAgca                                <  1:616781/69‑3 (MQ=255)
                      aaTACCTCGACGAACATGGCATCGTTGTTGAGAAAACCGGTCCGTATAACCTGCTGTTCCTGTTCAgca                                <  1:2663649/69‑3 (MQ=255)
                      aaTACCTCGACGAACATGGCATCGTTGTTGAGAAAACCGGTCCGTATAACCTGCTGTTCCTGTTCAgca                                <  1:2668939/69‑3 (MQ=255)
                      aaTACCTCGACGAACATGGCATCGTTGTTGAGAAAACCGGTCCGTATAACCTGCTGTTCCTGTTCAgca                                <  1:1888849/69‑3 (MQ=255)
                      aaTACCTCGACGAACATGGCATCGTTGTTGAGAAAACCGGTCCGTATAACCTGCTGTTCCTGTTCAgca                                <  1:1233416/69‑3 (MQ=255)
                      aaTACCTCGACGAACATGGCATCGTTGTTGAGAAAACCGGTCCGTATAACCTGCTGTTCCTGTTCAgca                                <  1:2403666/69‑3 (MQ=255)
                      aaTACCTCGACGAACATGGCATCGTTGTTGAGAAAACCGGTCCGTATAACCTGCTGTTCCTGTTCAgca                                <  1:1904420/69‑3 (MQ=255)
                      aaTACCTCGACGAACATGGCATCGTTGTTGAGAAAACCGGTCCGTATAACCTGCTGTTCCTGTTCAgca                                <  1:2151660/69‑3 (MQ=255)
                      aaTACCTCGACGAACATGGCATCGTTGTTGAGAAAACCGGTCCGTATAACCTGCTGTTCCTGTTCAgca                                <  1:1846561/69‑3 (MQ=255)
                      aaTACCTCGACGAACATGGCATCGTTGTTGAGAAAACCGGTCCGTATAACCTGCTGTTCCTGTTCAgca                                <  1:1419785/69‑3 (MQ=255)
                      aaTACCTCGACGAACATGGCATCGTTGTTGAGAAAACCGGTCCGTATAACCTGCTGTTCCTGTTCAgca                                <  1:1468180/69‑3 (MQ=255)
                      aaTACCTCGACGAACATGGCATCGTTGTTGAGAAAACCGGTCCGTATAACCTGCTGTTCCTGTTCAgca                                <  1:2072512/69‑3 (MQ=255)
                      aaTACCTCGACGAACATGGCATCGTTGTTGAGAAAACCGGTCCGTATAACCTGCTGTTCCTGTTCAgca                                <  1:1607158/69‑3 (MQ=255)
                      aaTACCTCGACGAACATGGCATCGTTGTTGAGAAAACCGGTCCGTATAACCTGCTGTTCCTGTTCAg                                  >  1:1528336/1‑67 (MQ=255)
                      aaTACCTCGACGAACATGGCATCGTTGTTGAGAAAACCGGTCCGTATAACCTGCTGTTCCTGTTCAg                                  >  1:3317817/1‑67 (MQ=255)
                      aaTACCTCGACGAACATGGCATCGTTGTTGAGAAAACCGGTCCGTATAACCTGCTGTTCCTGTTCAg                                  >  1:1791176/1‑67 (MQ=255)
                      aaTACCTCGACGAACATGGCATCGTTGTTGAGAAAACCGGTCCGTATAACCTGCTGTTCCTGTTCAg                                  >  1:397673/1‑67 (MQ=255)
                      aaTACCTCGACGAACATGGCATCGTTGTTGAGAAAACCGGTCCGTATAACCTGCTGTTCCTGTTCAg                                  >  1:566402/1‑67 (MQ=255)
                      aaTACCTCGACGAACATGGCATCGTTGTTGAGAAAACCGGTCCGTATAACCTGCTGTTCCTGTTCAg                                  >  1:2485818/1‑67 (MQ=255)
                      aaTACCTCGACGAACATGGCATCGTTGTTGAGAAAACCGGTCCGTATAACCTGCTGTTCCTGTTCAg                                  >  1:1307497/1‑67 (MQ=255)
                      aaTACCTCGACGAACATGGCATCGTTGTTGAGAAAACCGGTCCGTATAACCTGCTGTTCCTGTTCAg                                  >  1:838202/1‑67 (MQ=255)
                      aaTACCTCGACGAACATGGCATCGTTGTTGAGAAAACCGGTCCGTATAACCTGCTGTTCCTGTTCAg                                  >  1:1278889/1‑67 (MQ=255)
                      aaTACCTCGACGAACATGGCATCGTTGTTGAGAAAACCGGTCCGTATAACCTGCTGTTCCTGTTCAg                                  >  1:114925/1‑67 (MQ=255)
                      aaTACCTCGACGAACATGGCATCGTTGTTGAGAAAACCGGTCCGTATAACCTGCTGTTCCTGTTCAg                                  >  1:2017844/1‑67 (MQ=255)
                      aaTACCTCGACGAACATGGCATCGTTGTTGAGAAAACCGGTCCGTATAACCTGCTGTTCCTGTTCAg                                  >  1:3243803/1‑67 (MQ=255)
                      aaTACCTCGACGAACATGGCATCGTTGTTGAGAAAACCGGTCCGTATAACCTGCTGTTCCTGTTCAg                                  >  1:2225287/1‑67 (MQ=255)
                      aaTACCTCGACGAACATGGCATCGTTGTTGAGAAAACCGGTCCGTATAACCTGCTGTTCCTGTTCAg                                  >  1:3196280/1‑67 (MQ=255)
                      aaTACCTCGACGAACATGGCATCGTTGTTGAGAAAACCGGTCCGTATAACCTGCTGTTCCTGTTCAg                                  >  1:211024/1‑67 (MQ=255)
                      aaTACCTCGACGAACATGGCATCGTTGTTGAGAAAACCGGTCCGTATAACCTGCTGTTCCTGTTCAg                                  >  1:3183370/1‑67 (MQ=255)
                      aaTACCTCGACGAACATGGCATCGTTGTTGAGAAAACCGGTCCGTATAACCTGCTGTTCCTGTTCAg                                  >  1:209702/1‑67 (MQ=255)
                      aaTACCTCGACGAACATGGCATCGTTGTTGAGAAAACCGGTCCGTATAACCTGCTGTTCCTGTTCAg                                  >  1:3028215/1‑67 (MQ=255)
                                           tCGTAGTAGAGAAAACCGGCCCTTATAACctgctg                                             <  1:972656/35‑1 (MQ=255)
                                                  agagAAAACCGGCCCTTATAACCTGCTGTTTCTCTTTAGTATTGGCATCGATAAAACCAAAGc          <  1:152510/63‑1 (MQ=255)
                                                  agagAAAACCGGCCCTTATAACCTGCTGTTTCTCTTTAGTATTGGCATCGATAAAACCAAAGc          <  1:3073395/63‑1 (MQ=255)
                                                  agagAAAACCGGCCCTTATAACCTGCTGTTTCTCTTTAGTATTGGCATCGATAAAACCAAAGc          <  1:2927428/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:3201208/63‑1 (MQ=255)
                                                   gagaAAACCGGCCCTTATAACCTGCTGTTTCTCTTTAGTATTGGCATCGATAAAACCAAAGCAATggg    >  1:580706/1‑68 (MQ=255)
                                                   gagaAAACCGGCCCTTATAACCTGCTGTTTCTCTTTAGTATTGGCATCGATAAAACCAAAGCAATggg    >  1:65620/1‑68 (MQ=255)
                                                   gagaAAACCGGCCCTTATAACCTGCTGTTTCTCTTTAGTATTGGCATCGATAAAACCAAAGCAATggg    >  1:1583036/1‑68 (MQ=255)
                                                   gagaAAACCGGCCCTTATAACCTGCTGTTTCTCTTTAGTATTGGCATCGATAAAACCAAAGCAATggg    >  1:1655277/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)
                                                                 |                                                       
ATCCCGGCGGCGCTGGTAGCAAAATTCCTCGACGAACGTGGGATCGTAGTAGAGAAAACCGGCCCTTATAACCTGCTGTTTCTCTTTAGTATTGGCATCGATAAAACCAAAGCAATGGGAT  >  minE/190028‑190148

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: