Predicted mutation
evidence seq id position mutation freq annotation gene description
RA minE 190,075 A→T 16.3% V524V (GTA→GTT ldcC → lysine decarboxylase 2, constitutive

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*minE190,0750AT16.3% 121.5 / 15.7 49V524V (GTA→GTTldcClysine decarboxylase 2, constitutive
Reads supporting (aligned to +/- strand):  ref base A (27/14);  new base T (6/2);  total (33/16)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

GGGCAATATGAGCGAGGAGGGGATCCCGGCGGCGCTGGTAGCAAAATTCCTCGACGAACGTGGGATCGTAGTAGAGAAAACCGGCCCTTATAACCTGCTGTTTCTCTTTAGTATTGGCATCGATAAAACCAAAGCAA  >  minE/190006‑190142
                                                                     |                                                                   
gggCAATATGAGCGAGGAGGGGATCCCGGCGGCGCTGGTAGCAAAATTCCTCGACGAACGTGGGATCgtag                                                                    >  1:1280975/1‑71 (MQ=255)
  gCAATATGAGCGAGGAGGGGATCCCGGCGGCGCTGGTAGCAAAATTCCTCGACGAACGTGGGATCgtagta                                                                  <  1:2844800/71‑1 (MQ=255)
      tatGAGCGAGGAGGGGATCCCGGCGGCGCTGGTAGCAAAATTCCTCGACGAACGTGGGATCGTAGTagag                                                               >  1:3342698/1‑70 (MQ=255)
        tGAGCGAGGAGGGGATCCCGGCGGCGCTGGTAGCAAAATTCCTCGACGAACGTGGGATCgta                                                                     >  1:995718/1‑62 (MQ=255)
        tGAGCGAGGAGGGGATCCCGGCGGCGCTGGTAGCAAAATTCCTCGACGAACGTGGGATCgta                                                                     >  1:2987016/1‑62 (MQ=255)
        tGAGCGAGGAGGGGATCCCGGCGGCGCTGGTAGCAAAATTCCTCGACGAACGTGGGATCGTAGTAGAGaa                                                             >  1:927767/1‑70 (MQ=255)
                 aggGGATCCCGGCGGCGCTGGTAGCAAAATTCCTCGACGAACGTGGGATCGTAGTAGAGaaaa                                                           >  1:2771031/1‑63 (MQ=255)
                   gggATCCCGGCGGCGCTGGTAGCAAAATTCCTCGACGAACGTGGGATCGTAGTAGAGaaa                                                            <  1:2156612/60‑1 (MQ=255)
                      aTCCCGGCGGCGCTGGTAGCAAAATTCCTCGACGAACGTGGGATCGTAGTAGAGAAAACCGGCCCTTATaa                                              >  1:3132968/1‑71 (MQ=255)
                        cccGGCGGCGCTGGTAGCAAAATTCCTCGACGAACGTGGGATCgta                                                                     >  1:3784601/1‑46 (MQ=255)
                        cccGGCGGCGCTGGTAGCAAAATTCCTCGACGAACGTGGGATCGTAGTAGAGAAAAcc                                                         >  1:3721888/1‑58 (MQ=255)
                        cccGGCGGCGCTGGTAGCAAAATTCCTCGACGAACGTGGGATCGTAGTAGAGAAAACCGGCCCTTATAAcc                                            >  1:1751251/1‑71 (MQ=255)
                        cccGGCGGCGCTGGTAGCAAAATTCCTCGACGAACGTGGGATCGTAGTAGAGAAAACCGGCCCTTATAAcc                                            >  1:397263/1‑71 (MQ=255)
                        cccGGCGGCGCTGGTAGCAAAATTCCTCGACGAACGTGGGATCGTAGTAGAGAAAACCGGCCCTTATAAc                                             >  1:3751403/1‑70 (MQ=255)
                         ccGGCGGCGCTGGTAGCAAAATTCCTCGACGAACGTGGGATCGTAGTAGAGAAAAc                                                          >  1:1664714/1‑56 (MQ=255)
                           ggcggcGCTGGTAGCAAAATTCCTCGACGAACGTGGGATCGTAGTAGAGAAAACCGGCCCTTATAACctgc                                         <  1:3767392/71‑1 (MQ=255)
                            gcggagCTGGTAGCAAAATTCCTCGACGAACGTGGGATCGTAGTAGAGAAAACCGGCCCTTATAACctgct                                        >  1:2416242/1‑71 (MQ=255)
                                 gcTGGTAGCAAAATTCCTCGACGAACGTGGGATCGTAGTAGAGAAAACCGGCCCTTATAACCTGCTGTTtc                                   <  1:2875556/71‑1 (MQ=255)
                                    ggTAGCAAAATTCCTCGACGAACGTGGGATCGTAGTAGAGAAAACCGGCCCTTATAAc                                             >  1:3534410/1‑58 (MQ=255)
                                     gTAGCAAAATTCCTCGACGAACGTGGGATCGTAGTaga                                                                <  1:628541/38‑1 (MQ=255)
                                     gTAGCAAAATTCCTCGACGAACGTGGGATCGTAGTAGAGAAAACCGGCCCTTATAACCTGCTGTTtctc                                 <  1:4049635/69‑1 (MQ=255)
                                     gTAGCAAAATTCCTCGACGAACGTGGGATCGTAGTAGAGAAAACCGGCCCTTATAACCTGCTGTTTCTCtt                               <  1:2393807/71‑1 (MQ=255)
                                     gTAGCAAAATTCCTCGACGAACGTGGGATCGTAGTAGAGAAAACCGGCCCTTATAACCTGCTGTTTCTCtt                               <  1:1680610/71‑1 (MQ=255)
                                     gTAGCAAAATTCCTCGACGAACGTGGGATCGTAGTAGAGAAAACCGGCCCTTATAACCTGCTGTTTCTCtt                               >  1:2480545/1‑71 (MQ=255)
                                     gTAGCAAAATTCCTCGACGAACGTGGGATCGTAGTAGAGAAAACCGGCCCTTATAACCTGCTGTTTCTCtt                               <  1:1466335/71‑1 (MQ=255)
                                       aGCAAAATTCCTCGACGAACGTGGGATCGTAGTAGAGaa                                                             >  1:1306369/1‑39 (MQ=255)
                                       aGCAAAATTCCTCGACGAACGTGGGATCGTAGTAGAGAAAACCGGCCCTTATAACCTGCTGTTTCTCtt                               >  1:1335423/1‑69 (MQ=255)
                                        gCAAAATTCCTCGACGAACGTGGGATCGTAGTAGAGAAAACCGGCCCTTATAACCTGCTGTTTCTCTTTAg                            <  1:738377/71‑1 (MQ=255)
                                            aaTACCTCGACGAACATGGCATCGTTGTTGAGAAAACCGGTCCGTATAACCTGCTGTTCCTGTTCAGCAtc                        >  1:989448/1‑70 (MQ=255)
                                             ataccTCGACGAACATGGCATCGTTGTTGAGAAAACCGGTCCGTATAACCTGCTGTTCCTGTTCAGCAt                         <  1:1156892/66‑1 (MQ=255)
                                              ttCCTCGACGAACGTGGGATCGTAGTAGAGAAAAcc                                                         >  1:3688156/1‑36 (MQ=255)
                                              ttCCTCGACGAACGTGGGATCGTAGTAGAGAAAAcc                                                         >  1:1837572/1‑36 (MQ=255)
                                               accTCGACGAACATGGCATCGTTGTTGAGAAAACCGGTCCGTATAACCTGCTGTTcctg                                 <  1:2034573/58‑2 (MQ=19)
                                                ccTCGACGAACGTGGGATCGTAGTAGAGAAAACCGGCCCTTATAACCTGCTGTTTCTCTTTAGTATTGGCa                    >  1:2442605/1‑71 (MQ=255)
                                                 cTCGACGAACGTGGGATCGTAGTAGAGAAAACCGGCCCTTATAACCTGCTGTTTCTCTTTAGTATTGGCAt                   >  1:3564507/1‑71 (MQ=255)
                                                 cTCGACGAACGTGGGATCGTAGTAGAGAAAACCGGCCCTTATAACCTGCTGTTTCTCTTTAGTATTGGCAt                   >  1:3512026/1‑71 (MQ=255)
                                                  tCGACGAACGTGGGATCGTAGTAGAGAAAACCGGCCCTTATAACCTGCTGTTTCTCTTTAGTATTGGCATc                  >  1:3450749/1‑71 (MQ=255)
                                                    gacgaACGTGGGATCGTAGTAGAGAAAACCGGCCCTta                                                 >  1:2721051/1‑38 (MQ=255)
                                                    gacgaACGTGGGATCGTAGTAGAGAAAACCGGCCCTTATAACCTGCTGTTTCTCTTTAGTATTGGCATc                  <  1:1918597/69‑1 (MQ=255)
                                                    gacgaACATGGCATCGTTGTTGAGAAAACCGGTCCGTATAACCTGCTGTTCCTGTTCAGCATCGGTATCGa                >  1:1369840/1‑71 (MQ=255)
                                                       gaACGTGGGATCGTAGTAGAGAAAACCGGCCCTTATAACCTGCTGTTTCTCTTTAg                            <  1:178726/56‑1 (MQ=255)
                                                       gaACATGGCATCGTTGTTGAGAAAACCGGTCCGTATAACCTGCTGTTcct                                  >  1:1916936/1‑50 (MQ=16)
                                                       gaACATGGCATCGTTGTTGAGAAAACCGGTCCGTATAACCTGCTGTTcct                                  >  1:2411829/1‑50 (MQ=16)
                                                          cGTGGGATCGTAGTAGAGAAAACCGGCCCTTATAAc                                             >  1:1728748/1‑36 (MQ=255)
                                                          cGTGGGATCGTAGTAGAGAAAACCGGCCCTTATAACCTGCTGTTTCTCTTTAGTATTGGCATCGATaaaaa          <  1:2821963/71‑2 (MQ=255)
                                                           gTGGGATCGTAGTAGAGAAAACCGGCCCTTATAACCTGCTGTTTCTCTTTAGt                           <  1:2450675/53‑1 (MQ=255)
                                                             gggATCGTAGTAGAGAAAACCGGCCCTTATAACCTGCTGTTTCTCTTTa                             >  1:244781/1‑49 (MQ=255)
                                                             gggATCGTAGTAGAGAAAACCGGCCCTTATAACCTGCTGTTTCTCTTTAGTATTGGCATCGAt               >  1:3267638/1‑63 (MQ=255)
                                                             ggcaTCGTTGTTGAGAAAACCGGTCCGTATAACCTGCTGTTCCTGTTCAGCATCGGTATCGATAAGACCa        >  1:3703444/4‑70 (MQ=255)
                                                             ggcaTCGTTGTTGAGAAAACCGGTCCGTATAACCTGCTGTTCCTGTTCAGCATCGGTATCGATAAGACCa        >  1:1480316/4‑70 (MQ=255)
                                                              ggATCGTAGTAGAGAAAACCGGCCCTTATAACCTGCTGTTTCTCTTTAGTAt                         >  1:3199325/1‑52 (MQ=255)
                                                              ggATCGTAGTAGAGAAAACCGGCCCTTATAACCTGCTGTTTCTCTTTAGTATTGGCATCGATaaaa           >  1:1260951/1‑66 (MQ=255)
                                                                 tCGTAGTAGAGAAAACCGGCCCTTATAACCTGCTGTTTCTCTTTAGTATTGGCATCGATAAAACCAAAGCa   <  1:3279586/71‑1 (MQ=255)
                                                                  cGTAGTAGAGAAAACCGGCCCTTATAACCTGCTGTTTCTCTTTAGTATTGGCATCGATAAAACCAAAGCaa  >  1:4054129/1‑71 (MQ=255)
                                                                    tagtagAGAAAACCGGCCCTTATAACCTGCTGTTTCTCTTTAGt                           <  1:3837901/44‑1 (MQ=255)
                                                                     agtagAGAAAACCGGCCCTTATAACCTGCTGTTTCTCTTTAGTATTGGCATCGAt               >  1:4042529/1‑55 (MQ=255)
                                                                     |                                                                   
GGGCAATATGAGCGAGGAGGGGATCCCGGCGGCGCTGGTAGCAAAATTCCTCGACGAACGTGGGATCGTAGTAGAGAAAACCGGCCCTTATAACCTGCTGTTTCTCTTTAGTATTGGCATCGATAAAACCAAAGCAA  >  minE/190006‑190142

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: