Predicted mutation
evidence seq id position mutation freq annotation gene description
RA minE 190,273 A→G 34.9% A590A (GCA→GCG ldcC → lysine decarboxylase 2, constitutive

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*minE190,2730AG34.9% 85.9 / 75.9 126A590A (GCA→GCGldcClysine decarboxylase 2, constitutive
Reads supporting (aligned to +/- strand):  ref base A (52/30);  new base G (17/27);  total (69/57)
Fisher's exact test for biased strand distribution p-value = 9.05e-03
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.67e-01

Only 100 of 126 total aligned reads displayed.

ATGCTACCCGATCTCTATGCAGAAGATCCCGATTTCTACCGCAATATGCGTATTCAGGATCTGGCACAAGGGATCCATAAGCTGATTCGTAAACACGATCTTCCCGGTTTGATGTTGCGGGCATTC  >  minE/190208‑190333
                                                                 |                                                            
atGCTACCCGATCTCTATGCAGAAGATCCCGATTTCTACCGCAATATGCGTATTCAGGATCTGGCACaa                                                           >  1:400939/1‑69 (MQ=255)
atGCTACCCGATCTCTATGCAGAAGATCCCGATTTCTACCGCAATATGCGTATTCAGGATCTGGCACaa                                                           >  1:464870/1‑69 (MQ=255)
atGCTACCCGATCTCTATGCAGAAGATCCCGATTTCTACCGCAATATGCGTATTCAGGATCTGGCACaa                                                           <  1:1518737/69‑1 (MQ=255)
 tGCTACCCGATCTCTATGCAGAAGATCCCGATTTCTACCGCAATATGCGTATTCAGGATCTGGCACaa                                                           <  1:1833731/68‑1 (MQ=255)
   cTACCCGATCTCTATGCAGAAGATCCCGATTTCTACCGCAATATGCGTATTCAGGATCTGGCACAAggg                                                        <  1:344972/69‑1 (MQ=255)
   cTACCCGATCTCTATGCAGAAGATCCCGATTTCTACCGCAATATGCGTATTCAGGATCTGGCACAAggg                                                        <  1:2677615/69‑1 (MQ=255)
   cTACCCGATCTCTATGCAGAAGATCCCGATTTCTACCGCAATATGCGTATTCAGGATCTGGCACAAggg                                                        <  1:1118247/69‑1 (MQ=255)
   cTACCCGATCTCTATGCAGAAGATCCCGATTTCTACCGCAATATGCGTATTCAGGATCTGGCACAAggg                                                        <  1:2570377/69‑1 (MQ=255)
   cTACCCGATCTCTATGCAGAAGATCCCGATTTCTACCGCAATATGCGTATTCAGGATCTGGCACAAggg                                                        <  1:2495121/69‑1 (MQ=255)
    tACCCGATCTCTATGCAGAAGATCCCGATTTCTACCGCAATATGCGTATTCAGGATCTGGCACAAggg                                                        <  1:859711/68‑1 (MQ=255)
            ctctATGCAGAAGATCCCGATTTCTACCGCAATATGCGTATTCAGGATCTGGCACAAGGGATCCATAAg                                               >  1:1584889/1‑69 (MQ=255)
            ctctATGCAGAAGATCCCGATTTCTACCGCAATATGCGTATTCAGGATCTGGCACAAGGGATCCATAAg                                               >  1:1423621/1‑69 (MQ=255)
              ctATGCAGAAGATCCCGATTTCTACCGCAATATGCGTATTCAGGATCTGGCACAAGGGATCCATAAGCt                                             >  1:2760614/1‑69 (MQ=255)
              ctATGCAGAAGATCCCGATTTCTACCGCAATATGCGTATTCAGGATCTGGCACAAGGGATCCATAAGCt                                             >  1:1847208/1‑69 (MQ=255)
              ctATGCAGAAGATCCCGATTTCTACCGCAATATGCGTATTCAGGATCTGGCACAAGGGATCCATAAGCt                                             >  1:3237454/1‑69 (MQ=255)
              ctATGCAGAAGATCCCGATTTCTACCGCAATATGCGTATTCAGGATCTGGCACAAGGGATCCATAAGCt                                             >  1:3036684/1‑69 (MQ=255)
              ctATGCAGAAGATCCCGATTTCTACCGCAATATGCGTATTCAGGATCTGGCACAAGGGATCCATAAGCt                                             >  1:1446452/1‑69 (MQ=255)
              ctATGCAGAAGATCCCGATTTCTACCGCAATATGCGTATTCAGGATCTGGCACAAGGGATCCATAAGCt                                             >  1:2992233/1‑69 (MQ=255)
              ctATGCAGAAGATCCCGATTTCTACCGCAATATGCGTATTCAGGATCTGGCACAAGGGATCCATAAGCt                                             >  1:1068758/1‑69 (MQ=255)
                 tGCAGAAGATCCCGATTTCTACCGCAATATGCGTATTCAGGATCTGGCACAAGGGATCCATAAGCt                                             <  1:3389400/66‑1 (MQ=255)
                 tGCAGAAGATCCCGATTTCTACCGCAATATGCGTATTCAGGATCTGGCACAAGGGATCCATAAGCt                                             <  1:2676414/66‑1 (MQ=255)
                 tGCAGAAGATCCCGATTTCTACCGCAATATGCGTATTCAGGATCTGGCACAAGGGATCCATAAGCt                                             <  1:2592164/66‑1 (MQ=255)
                 tGCAGAAGATCCCGATTTCTACCGCAATATGCGTATTCAGGATCTGGCACAAGGGATCCATAAGCt                                             <  1:1152762/66‑1 (MQ=255)
                            ccGATTTCTACCGCAATATGCGTATTCAGGATCTGGCACAAGGGATCCATAAGCTGATTCGTAAacac                                >  1:3430005/1‑68 (MQ=255)
                            ccGATTTCTACCGCAATATGCGTATTCAGGATCTGGCACAAGGGATCCATAAGCTGATTCGTAAacac                                >  1:773085/1‑68 (MQ=255)
                            ccGATTTCTACCGCAATATGCGTATTCAGGATCTGGCACAAGGGATCCATAAGCTGATTCGTAAACACg                               >  1:2340818/1‑69 (MQ=255)
                                 ttCTACCGCAATATGCGTATTCAGGATCTGGCACAAGGGATCCATAAGCTGATTCGTAAACACGATCtt                          >  1:2269702/1‑69 (MQ=255)
                                 ttCTACCGCAATATGCGTATTCAGGATCTGGCACAAGGGATCCATAAGCTGATTCGTAAACACGATCtt                          >  1:2129342/1‑69 (MQ=255)
                                         cAATATGCGTATTCAGGATCTGGCGCAAGGGATCCATAAGCTGATTCGTAAACACGATCTTCCCGGttt                  >  1:2671164/1‑69 (MQ=255)
                                         cAATATGCGTATTCAGGATCTGGCGCAAGGGATCCATAAGCTGATTCGTAAACACGATCTTCCCGGttt                  >  1:2580735/1‑69 (MQ=255)
                                         cAATATGCGTATTCAGGATCTGGCGCAAGGGATCCATAAGCTGATTCGTAAACACGATCTTCCCGGttt                  >  1:2372515/1‑69 (MQ=255)
                                         cAATATGCGTATTCAGGATCTGGCGCAAGGGATCCATAAGCTGATTCGTAAACACGATCTTCCCGGttt                  >  1:2374597/1‑69 (MQ=255)
                                         cAATATGCGTATTCAGGATCTGGCGCAAGGGATCCATAAGCTGATTCGTAAACACGATCTTCCCGGttt                  >  1:2793038/1‑69 (MQ=255)
                                         cAATATGCGTATTCAGGATCTGGCGCAAGGGATCCATAAGCTGATTCGTAAACACGATCTTCCCGGttt                  >  1:207220/1‑69 (MQ=255)
                                         cAATATGCGTATTCAGGATCTGGCGCAAGGGATCCATAAGCTGATTCGTAAACACGATCTTCCCGGttt                  >  1:1840559/1‑69 (MQ=255)
                                         cAATATGCGTATTCAGGATCTGGCGCAAGGGATCCATAAGCTGATTCGTAAACACGATCTTCCCGGttt                  >  1:1840267/1‑69 (MQ=255)
                                         cAATATGCGTATTCAGGATCTGGCGCAAGGGATCCATAAGCTGATTCGTAAACACGATCTTCCCGGttt                  >  1:11311/1‑69 (MQ=255)
                                         cAATATGCGTATTCAGGATCTGGCGCAAGGGATCCATAAGCTGATTCGTAAACACGATCTTCCCGGttt                  >  1:1203807/1‑69 (MQ=255)
                                         cAATATGCGTATTCAGGATCTGGCGCAAGGGATCCATAAGCTGATTCGTAAACACGATCTTCCCGGttt                  >  1:1786230/1‑69 (MQ=255)
                                         cAATATGCGTATTCAGGATCTGGCGCAAGGGATCCATAAGCTGATTCGTAAACACGATCTTCCCGGttt                  >  1:1642012/1‑69 (MQ=255)
                                         cAATATGCGTATTCAGGATCTGGCGCAAGGGATCCATAAGCTGATTCGTAAACACGATCTTCCCGGttt                  >  1:939919/1‑69 (MQ=255)
                                         cAATATGCGTATTCAGGATCTGGCGCAAGGGATCCATAAGCTGATTCGTAAACACGATCTTCCCGGttt                  >  1:1580771/1‑69 (MQ=255)
                                         cAATATGCGTATTCAGGATCTGGCACAAGGGATCCATAAGCTGATTCGTAAACACGATCTTCCCGGttt                  <  1:940172/69‑1 (MQ=255)
                                           atatGCGTATTCAGGATCTGGCACAAGGGATCCATAAGCTGATTCGTAAACACGATCTTCCCGGTTTGa                >  1:3322236/1‑69 (MQ=255)
                                            tatGCGTATTCAGGATCTGGCACAAGGGATCCATAAGCTGATTCGTAAACACGATCTTCCCGGTTTGAt               <  1:3391591/69‑1 (MQ=255)
                                            tatGCGTATTCAGGATCTGGCACAAGGGATCCATAAGCTGATTCGTAAACACGATCTTCCCGGTTTGAt               <  1:2919192/69‑1 (MQ=255)
                                            tatGCGTATTCAGGATCTGGCACAAGGGATCCATAAGCTGATTCGTAAACACGATCTTCCCGGTTTGAt               <  1:74208/69‑1 (MQ=255)
                                            tatGCGTATTCAGGATCTGGCACAAGGGATCCATAAGCTGATTCGTAAACACGATCTTCCCGGTTTGAt               <  1:596605/69‑1 (MQ=255)
                                            tatGCGTATTCAGGATCTGGCACAAGGGATCCATAAGCTGATTCGTAAACACGATCTTCCCGGTTTGAt               <  1:2182004/69‑1 (MQ=255)
                                            tatGCGTATTCAGGATCTGGCACAAGGGATCCATAAGCTGATTCGTAAACACGATCTTCCCGCTTTGAt               <  1:438313/69‑1 (MQ=255)
                                            tatGCGTATTCAGGATCTATCACAAGGGATCCATAAGCTGATTCGTAAACACGATCTTCCCGGTTTGAt               <  1:817400/69‑1 (MQ=255)
                                              tGCGTATTCAGGATCTGGCACAAGGGATCCATAAGCTg                                            >  1:2687959/1‑38 (MQ=255)
                                              tGCGTATTCAGGATCTGGCACAAGGGATCCATAAGCTGATTCGTaaa                                   >  1:704387/1‑47 (MQ=255)
                                              tGCGTATTCAGGATCTGGCACAAGGGATCCATAAGCTGATTCGTAAACACGATCTTc                         >  1:3364905/1‑57 (MQ=255)
                                              tGCGTATTCAGGATCTGGCACAAGGGATCCATAAGCTGATTCGTAAACACGATCTTCCCGGTTTGATg              >  1:1093129/1‑68 (MQ=255)
                                              tGCGTATTCAGGATCTGGCACAAGGGATCCATAAGCTGATTCGTAAACACGATCTTCCCGGTTTGATg              >  1:2101893/1‑68 (MQ=255)
                                              tGCGTATTCAGGATCTGGCACAAGGGATCCATAAGCTGATTCGTAAACACGATCTTCCCGGTTTGATGt             >  1:2999728/1‑69 (MQ=255)
                                              tGCGTATTCAGGATCTGGCACAAGGGATCCATAAGCTGATTCGTAAACACGATCTTCCCGGTTTGATGt             >  1:682406/1‑69 (MQ=255)
                                              tGCGTATTCAGGATCTGGCACAAGGGATCCATAAGCTGATTCGTAAACACGATCTTCCCGGTTTGATGt             >  1:548204/1‑69 (MQ=255)
                                              tGCGTATTCAGGATCTGGCACAAGGGATCCATAAGCTGATTCGTAAACACGATCTTCCCGGTTTGATGt             >  1:3245885/1‑69 (MQ=255)
                                              tGCGTATTCAGGATCTGGCACAAGGGATCCATAAGCTGATTCGTAAACACGATCTTCCCGGTTTGATGt             >  1:2895109/1‑69 (MQ=255)
                                              tGCGTATTCAGGATCTGGCACAAGGGATCCATAAGCTGATTCGTAAACACGATCTTCCCGGTTTGATGt             >  1:2540216/1‑69 (MQ=255)
                                              tGCGTATTCAGGATCTGGCACAAGGGATCCATAAGCTGATTCGTAAACACGATCTTCCCGGTTTGATGt             >  1:1999496/1‑69 (MQ=255)
                                              tGCGTATTCAGGATCTGGCACAAGGGATCCATAAGCTGATTCGTAAACACGATCTTCCCGGTTTGATGt             >  1:1031380/1‑69 (MQ=255)
                                              tGCGTATTCAGGATCTGGCACAAGGGATCCATAAGCTGATTCGTAAACACGATCTTCCCGGTTTGATGt             >  1:1303644/1‑69 (MQ=255)
                                              tGCGTATTCAGGATCTGGCACAAGGGATCCATAAGCTGATTCGTAAACACGATCTTCCCGGTTTGATGt             >  1:1410277/1‑69 (MQ=255)
                                              tGCGTATTCAGGATCTGGCACAAGGGATCCATAAGCTGATTCGTAAACACGATCTTCCCGGTTTGATGt             >  1:1482493/1‑69 (MQ=255)
                                              tGCGTATTCAGGATCTGGCACAAGGGATCCATAAGCTGATTCGTAAACACGATCTTCCCGGTTTGATGt             >  1:2115210/1‑69 (MQ=255)
                                              tGCGTATTCAGGATCTGGCACAAGGGATCCATAAGCTGATTCGTAAACACGATCTTCCCGGTTTGATGt             >  1:2130087/1‑69 (MQ=255)
                                              tGCGTATTCAGGATCTGGCACAAGGGATCCATAAGCTGATTCGTAAACACGATCTTCCCGGTTTGATGt             >  1:2721333/1‑69 (MQ=255)
                                              tGCGTATTCAGGATCTGGCACAAGGGATCCATAAGCTGATTCGTAAACACGATCTTCCCGGTTTGATGt             >  1:2410641/1‑69 (MQ=255)
                                                cGTATTCAGGATCTGGCGCAAGGGATCCATAAGCTGATTCGTAAACACGATCTTCCCGGTTTGATGTTg           <  1:1594606/69‑1 (MQ=255)
                                                cGTATTCAGGATCTGGCGCAAGGGATCCATAAGCTGATTCGTAAACACGATCTTCCCGGTTTGATGTTg           <  1:2252667/69‑1 (MQ=255)
                                                cGTATTCAGGATCTGGCGCAAGGGATCCATAAGCTGATTCGTAAACACGATCTTCCCGGTTTGATGTTg           <  1:144571/69‑1 (MQ=255)
                                                cGTATTCAGGATCTGGCGCAAGGGATCCATAAGCTGATTCGTAAACACGATCTTCCCGGTTTGATGTTg           <  1:1297054/69‑1 (MQ=255)
                                                cGTATTCAGGATCTGGCGCAAGGGATCCATAAGCTGATTCGTAAACACGATCTTCCCGGTTTGATGTTg           <  1:2463030/69‑1 (MQ=255)
                                                 gTATTCAGGATCTGGCGCAAGGGATCCATAAGCTGATTCGTAAACACGATCTTCCCTGTTTGATGTTg           <  1:2734357/68‑1 (MQ=255)
                                                 gTATTCAGGATCTGGCGCAAGGGATCCATAAGCTGATTCGTAAACACGATCTTCCCGGTTTGATGTTg           <  1:1388090/68‑1 (MQ=255)
                                                  tATTCAGGATCTGGCACAAGGGATCCATAAGCTGATTCGTAAACACGATCTTCCCGGTTTGATGTTGCg         >  1:1493102/1‑69 (MQ=255)
                                                         gATCTGGCACAAGGGATCCATAAGCTGATTCGTAAACACGATCTTCCCGGTTTGATGTTGCGGGCATTc  >  1:2072335/1‑69 (MQ=255)
                                                         gATCTGGCACAAGGGATCCATAAGCTGATTCGTAAACACGATCTTCCCGGTTTGATGTTGCGGGCATTc  >  1:814995/1‑69 (MQ=255)
                                                         gATCTGGCACAAGGGATCCATAAGCTGATTCGTAAACACGATCTTCCCGGTTTGATGTTGCGGGCATTc  >  1:221428/1‑69 (MQ=255)
                                                         gATCTGGCACAAGGGATCCATAAGCTGATTCGTAAACACGATCTTCCCGGTTTGATGTTGCGGGCATTc  >  1:2230185/1‑69 (MQ=255)
                                                         gATCTGGCACAAGGGATCCATAAGCTGATTCGTAAACACGATCTTCCCGGTTTGATGTTGCGGGCATTc  >  1:1463675/1‑69 (MQ=255)
                                                           tCTGGCGCAAGGGATCCATAAGCTGATTCGTAAACACGATCTTCCCgg                     <  1:2900252/48‑1 (MQ=255)
                                                           tCTGGCGCAAGGGATCCATAAGCTGATTCGTAAACACGATCTTCCCgg                     <  1:918428/48‑1 (MQ=255)
                                                           tCTGGCGCAAGGGATCCATAAGCTGATTCGTAAACACGATCTTCCCgg                     <  1:3203002/48‑1 (MQ=255)
                                                           tCTGGCGCAAGGGATCCATAAGCTGATTCGTAAACACGATCTTCCCgg                     <  1:1018527/48‑1 (MQ=255)
                                                           tCTGGCGCAAGGGATCCATAAGCTGATTCGTAAACACGATCTTCCCgg                     <  1:1554158/48‑1 (MQ=255)
                                                           tCTGGCGCAAGGGATCCATAAGCTGATTCGTAAACACGATCTTCCCgg                     <  1:507541/48‑1 (MQ=255)
                                                           tCTGGCGCAAGGGATCCATAAGCTGATTCGTAAACACGATCTTCCCgg                     <  1:1712420/48‑1 (MQ=255)
                                                           tCTGGCGCAAGGGATCCATAAGCTGATTCGTAAACACGATCTTCCCgg                     <  1:1938270/48‑1 (MQ=255)
                                                           tCTGGCGCAAGGGATCCATAAGCTGATTCGTAAACACGATCTTCCCgg                     <  1:2893850/48‑1 (MQ=255)
                                                           tCTGGCGCAAGGGATCCATAAGCTGATTCGTAAACACGATCTTCCCgg                     <  1:3402950/48‑1 (MQ=255)
                                                           tCTGGCGCAAGGGATCCATAAGCTGATTCGTAAACACGATCTTCCCgg                     <  1:3369499/48‑1 (MQ=255)
                                                           tCTGGCGCAAGGGATCCATAAGCTGATTCGTAAACACGATCTTCCCgg                     <  1:3339878/48‑1 (MQ=255)
                                                           tCTGGCGCAAGGGATCCATAAGCTGATTCGTAAACACGATCTTCCCgg                     <  1:3294002/48‑1 (MQ=255)
                                                           tCTGGCGCAAGGGATCCATAAGCTGATTCGTAAACACGATCTTCCCgg                     <  1:305759/48‑1 (MQ=255)
                                                             tGGCGCAAGGGATCCATAAGCTGATTCGTAAACACGATCTTCCCgg                     <  1:476076/46‑1 (MQ=255)
                                                              ggCACAAGGGATCCATAAGCTGATTCGTAAACAAGATCTTCCCGGTTTGAt               <  1:1628340/51‑1 (MQ=255)
                                                                 |                                                            
ATGCTACCCGATCTCTATGCAGAAGATCCCGATTTCTACCGCAATATGCGTATTCAGGATCTGGCACAAGGGATCCATAAGCTGATTCGTAAACACGATCTTCCCGGTTTGATGTTGCGGGCATTC  >  minE/190208‑190333

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: