Predicted mutation
evidence seq id position mutation freq annotation gene description
RA minE 190,090 C→T 48.4% G529G (GGC→GGT ldcC → lysine decarboxylase 2, constitutive

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*minE190,0900CT48.4% 13.4 / 79.0 63G529G (GGC→GGTldcClysine decarboxylase 2, constitutive
Reads supporting (aligned to +/- strand):  ref base C (18/14);  new base T (8/22);  total (26/37)
Fisher's exact test for biased strand distribution p-value = 2.26e-02
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

AGGAGGGGATCCCGGCGGCGCTGGTAGCAAAATTCCTCGACGAACGTGGGATCGTAGTAGAGAAAACCGGCCCTTATAACCTGCTGTTTCTCTTTAGTATTGGCATCGATAAAACCAAAGCAATGGGATTATTG  >  minE/190020‑190153
                                                                      |                                                               
aggaggGGATCCCGGCGGCGCTGGTAGCAAAATTCCTCGACGAACGTGGGATCGTAGTAGAGAAAACCGGc                                                                 >  1:3891542/1‑71 (MQ=255)
            cggcggCGCTGGTAGCAAAATTCCTCGACGAACGTGGGATCGTAGTAGAGAAAACCGGCCCTTATAACctg                                                     >  1:4643796/1‑71 (MQ=255)
             ggcggcGCTGGTAGCAAAATTCCTCGACGAACGTGGGATCGTAGTAGAGAAAACCGGCCCTta                                                            <  1:1127544/63‑1 (MQ=255)
             ggcggcGCTGGTAGCAAAATTCCTCGACGAACGTGGGATCGTAGTAGAGAAAACCGGCCCTTATAACctgc                                                    <  1:1885781/71‑1 (MQ=255)
                     tGGTAGCAAAATTCCTCGACGAACGTGGGATCGTAGTAGAGAAAACCGGCCCTTATAACCTGCTGTTtctc                                            >  1:5501753/1‑71 (MQ=255)
                       gTAGCAAAATTCCTCGACGAACGTGGGATCGTAGTAGAGAAAACCGGcc                                                                <  1:3979568/49‑1 (MQ=255)
                       gTAGCAAAATTCCTCGACGAACGTGGGATCGTAGTAGAGAAAACCGGcc                                                                <  1:3856959/49‑1 (MQ=255)
                           cAAAATTCCTCGACGAACGTGGGATCGTAGTAGAGAAAACCGGCCCtt                                                             <  1:4343062/48‑1 (MQ=255)
                           cAAAATTCCTCGACGAACGTGGGATCGTAGTAGAGAAAACCGGCCCtt                                                             <  1:3965627/48‑1 (MQ=255)
                           cAAAATTCCTCGACGAACGTGGGATCGTAGTAGAGAAAACCGGCCCTTATAACCTGCTGTTTCTCTTTa                                        <  1:4080404/69‑1 (MQ=255)
                           cAAAATTCCTCGACGAACGTGGGATCGTAGTAGAGAAAACCGGCCCTTATAACCTGCTGTTTCTCTTTAg                                       >  1:3135692/1‑70 (MQ=255)
                            gaaaTACCTCGACGAACATGGCATCGTTGTTGAGAAAACCGGTCCGTATAACCTGCTGTTCCTGTTCAgca                                     <  1:1921034/70‑3 (MQ=255)
                            gaaaTACCTCGACGAACATGGCATCGTTGTTGAGAAAACCGGTCCGTATAACCTGCTGTTCCTGTTCAgca                                     <  1:2563564/70‑3 (MQ=255)
                            aaaaTTCCTCGACGAACGTGGGATCGTAGTAGAGAAAACCGGCCCTTATAACCTGCTGTTTCTCTTTAGTa                                     >  1:5446262/1‑71 (MQ=255)
                             aaaTACCTCGACGAACATGGCATCGTTGTTGAGAAAACCGGTCCGTATAACCTGCTGTTCCTGttc                                         <  1:212265/66‑2 (MQ=255)
                              aaTACCTCGACGAACATGGCATCGTTGTTGAGAAAACCGGTCCGTATAACCTGCTGTTCCTGTTCAGCAtc                                   >  1:3633037/1‑70 (MQ=255)
                                ttCCTCGACGAACGTGGGATCGTAGTAGAGAAAACCGGCCCTTATAACCTGCTGTTTCTCTTTAGTATTgg                                 >  1:705990/1‑71 (MQ=255)
                                taccTCGACGAACATGGCATCGTTGTTGAGAAAACCGGTCCGTATAACCTGCTGTTCCTGTTCAGCATCgg                                 <  1:164118/69‑1 (MQ=255)
                                   cTCGACGAACATGGCATCGTTGTTGAGAAAACCGGTCCGTATAACCTGCTGTTCCTGTTCAGCATCGGTAt                              >  1:880512/1‑71 (MQ=255)
                                     cgacgaACGTGGGATCGTAGTAGAGAAAACCGGCCCTTATAACCTGCTGTTTCTCTTTAGTAt                                    <  1:973463/63‑1 (MQ=255)
                                     cgacgaACGTGGGATCGTAGTAGAGAAAACCGGCCCTTATAACCTGCTGTTTCTCTTTAGTATTGGCATc                             >  1:682283/1‑70 (MQ=255)
                                      gacgaACGTGGGATCGTAGTAGAGAAAACCGGCCCTTATAAc                                                        >  1:5039417/1‑42 (MQ=255)
                                      gacgaACGTGGGATCGTAGTAGAGAAAACCGGCCCTTATAAc                                                        >  1:2433211/1‑42 (MQ=255)
                                      gacgaACATGGCATCGTTGTTGAGAAAACCGGTCCGTATAACCTGCTGTTCCTGTTCAGCATCGGTATc                             <  1:4136023/69‑1 (MQ=255)
                                      gacgaACATGGCATCGTTGTTGAGAAAACCGGTCCGTATAACCTGCTGTTCCTGTTCAGCATCGGTATc                             <  1:288101/69‑1 (MQ=255)
                                       acgaACATGGCATCGTTGTTGAGAAAACCGGTCCGTATAACCTGCTGTTCCTGTTCAGCATCGGTATCGa                           >  1:3479110/1‑70 (MQ=255)
                                       acgaACATGGCATCGTTGTTGAGAAAACCGGTCCGTATAACCTGCTGTTCCTGTTCAGCATCGGTATCGAt                          <  1:3428638/71‑1 (MQ=255)
                                       acgaACATGGCATCGTTGTTGAGAAAACCGGTCCGTATAACCTGCTGTTCCTGTTCAGCATCGGTATCGAt                          <  1:2344690/71‑1 (MQ=255)
                                        cgaACGTGGGATCGTAGTAGAGAAAACCGGCCCTtata                                                          <  1:3088921/38‑1 (MQ=255)
                                         gaACGTGGGATCGTAGTAGAGAAAACCGGCCCTTATAACCTGCTGTTTCTCTTTAGTATTGGCATCGATa                         >  1:4932308/1‑70 (MQ=255)
                                         gaACATGGCATCGTTGTTGAGAAAACCGGTCCGTATAACCTGCTGTTCCTGTTCAGCAtc                                   >  1:1705181/1‑59 (MQ=12)
                                          aaCATGGCATCGTTGTTGAGAAAACCGGTCCGTATAACctgc                                                    <  1:1434562/42‑1 (MQ=12)
                                           aCGTGGGATCGTAGTAGAGAAAACCGGCCCTTATAACctgc                                                    >  1:176999/1‑41 (MQ=255)
                                           aCGTGGGATCGTAGTAGAGAAAACCGGCCCTTATAACCTGCTGTTTCTCTTTAGTATTGGc                                >  1:4708487/1‑61 (MQ=255)
                                             atGGCATCGTTGTTGAGAAAACCGGTCCGTATAACCTGCTGTTCCTGTTCAGCATCGGTATCGATAAGAcc                    <  1:4827765/70‑1 (MQ=255)
                                             atGGCATCGTTGTTGAGAAAACCGGTCCGTATAACCTGCTGTTCCTGTTCAGCATCGGTATCGATAAGAcc                    <  1:719281/70‑1 (MQ=255)
                                              tGGCATCGTTGTTGAGAAAACCGGTCCGTATAACCTGCTGTTCCTGTTCAGCATCGGTATCg                            <  1:1265055/62‑1 (MQ=12)
                                              tGGCATCGTTGTTGAGAAAACCGGTCCGTATAACCTGCTGTTCCTGTTCAGCATCGGTATCg                            <  1:4706954/62‑1 (MQ=12)
                                              tGGCATCGTTGTTGAGAAAACCGGTCCGTATAACCTGCTGTTCCTGTTCAGCATCGGTATCg                            <  1:1373705/62‑1 (MQ=12)
                                              tGGCATCGTTGTTGAGAAAACCGGTCCGTATAACCTGCTGTTCCTGTTCAGCATCGGTATCGATAAGACCa                   <  1:1566860/71‑1 (MQ=255)
                                                ggATCGTAGTAGAGAAAACCGGCCCTTATAACCTGCTGTTTCTCTTTAGTATTGGCATCGATAAAACCa                   >  1:2260222/1‑69 (MQ=255)
                                                 caTCGTTGTTGAGAAAACCGGTCCGTATAACCTGCTGTTCCTGTTCAGCAtc                                   <  1:608460/51‑2 (MQ=12)
                                                 caTCGTTGTTGAGAAAACCGGTCCGTATAACCTGCTGTTCCTGTTCAGCATCGGTATCGATAAGACCAAAg                <  1:5253021/70‑1 (MQ=255)
                                                 caTCGTTGTTGAGAAAACCGGTCCGTATAACCTGCTGTTCCTGTTCAGCATCGGTATCGATAAGACCAAAg                <  1:2596964/70‑1 (MQ=255)
                                                  aTCGTTGTTGAGAAAACCGGTCCGTATAACCTGCTGTTCCTGTTCAGCATCGGTATCGATAAGACCAAAGc               <  1:2341838/71‑1 (MQ=255)
                                                   tCGTTGTTGAGAAAACCGGTCCGTATAACCTGCTGTTCCTGTTCAg                                       <  1:5287639/46‑1 (MQ=12)
                                                   tCGTAGTAGAGAAAACCGGCCCTTATAACCTGCTGTTTCTCTTTAGTATTGGCATCGATAAAACCAAAGCa              <  1:1695552/71‑1 (MQ=255)
                                                         ttgagaAAACCGGTCCGTATAACCTGCTGTTCCTGTTCAGCATCGGTATCGAt                          <  1:4744508/51‑1 (MQ=255)
                                                          agagAAAACCGGCCCTTATAACCTGCTGTTTCTCTTTAGTATTGGCATCGATAAAACCAAAGCAATGGGat       >  1:2420110/1‑71 (MQ=255)
                                                           gagaAAACCGGCCCTTATAACCTGCTGTTTCTCTTTAGTATTGGCATCGATaaa                       >  1:4133317/1‑54 (MQ=255)
                                                           gagaAAACCGGCCCTTATAACCTGCTGTTTCTCTTTAGTATTGGCATCGATAAAACCaaa                 >  1:2881683/1‑60 (MQ=255)
                                                           gagaAAACCGGCCCTTATAACCTGCTGTTTCTCTTTAGTATTGGCATCGATAAAACCaaa                 >  1:644475/1‑60 (MQ=255)
                                                           gagaAAACCGGCCCTTATAACCTGCTGTTTCTCTTTAGTATTGGCATCGATAAAACCAAAGCAATGGGatt      <  1:2332978/71‑1 (MQ=255)
                                                            agaAAACCGGCCCTTATAACCTGCTGTTTCTCTTTAGTa                                     <  1:4684909/39‑1 (MQ=255)
                                                            agaAAACCGGCCCTTATAACCTGCTGTTTCTCTTTAGTa                                     <  1:4322470/39‑1 (MQ=255)
                                                             gaAAACCGGCCCTTATAACCTGCTGTTTCTCTTTAGTa                                     >  1:4259462/1‑38 (MQ=255)
                                                             gaAAACCGGCCCTTATAACCTGCTGTTTCTCTTTAGTATTGGCATCGATAAAACCAAAGCaa             <  1:314931/62‑1 (MQ=255)
                                                              aaaaCCGGGCCTTATAACCTGCTGTTTCTCTTTAGTa                                     <  1:5341206/37‑1 (MQ=255)
                                                              aaaaCCGGCCCTTATAACCTGCTGTTTCTCTTTAGTAt                                    >  1:5043411/1‑38 (MQ=255)
                                                               aaaCCGGTCCGTATAACCTGCTGTTCCTGTTCAGCATCGGTATCGATAAGAc                     <  1:299019/52‑1 (MQ=255)
                                                               aaaCCGGCCCTTATAACCTGCTGTTTCTCTTTAGTATTGGCATCGATAAAACCAAAGCAATGGGATTATTg  <  1:1636107/71‑1 (MQ=255)
                                                                aaCCGGTCCGTATAACCTGCTGTTCCTGTTCAGCATCGGTATCGa                           >  1:3036080/1‑45 (MQ=255)
                                                                aaCCGGTCCGTATAACCTGCTGTTCCTGTTCAGCATCGGTATCGAt                          >  1:2727915/1‑46 (MQ=255)
                                                                aaCCGGTCCGTATAACCTGCTGTTCCTGTTCAGCATCGGTATCGATAAGACCAAAGCa              >  1:3654314/1‑58 (MQ=255)
                                                                 aCCGGTCCGTATAACCTGCTGTTCCTGTTCAGCATCGGTATCGATAAGACCAAAGCac             >  1:4835337/1‑57 (MQ=255)
                                                                     gCCCTTATAACCTGCTGTTTCTCTTTAGTATTGGCATCGATAAAACCAAAGCAATGGGatt      <  1:4273272/61‑1 (MQ=255)
                                                                      |                                                               
AGGAGGGGATCCCGGCGGCGCTGGTAGCAAAATTCCTCGACGAACGTGGGATCGTAGTAGAGAAAACCGGCCCTTATAACCTGCTGTTTCTCTTTAGTATTGGCATCGATAAAACCAAAGCAATGGGATTATTG  >  minE/190020‑190153

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: