Predicted mutation
evidence seq id position mutation freq annotation gene description
RA minE 190,120 T→C 12.3% I539I (ATT→ATC ldcC → lysine decarboxylase 2, constitutive

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*minE190,1200TC12.3% 195.8 / 17.9 73I539I (ATT→ATCldcClysine decarboxylase 2, constitutive
Reads supporting (aligned to +/- strand):  ref base T (35/29);  new base C (5/4);  total (40/33)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.85e-01

CCTCGACGAACGTGGGATCGTAGTAGAGAAAACCGGCCCTTATAACCTGCTGTTTCTCTTTAGTATTGGCATCGATAAAACCAAAGCAATGGGATTATTGCGTGGGTTGACGGAATTCAAACGCTCTTACGATCTC  >  minE/190054‑190189
                                                                  |                                                                     
ccTCGACGAACGTGGGATCGTAGTAGAGAAAACCGGCCCTTATAACCTGCTGTTTCTCTTTAGTATTGGCa                                                                   >  1:2442605/1‑71 (MQ=255)
 cTCGACGAACGTGGGATCGTAGTAGAGAAAACCGGCCCTTATAACCTGCTGTTTCTCTTTAGTATTGGCAt                                                                  >  1:3512026/1‑71 (MQ=255)
 cTCGACGAACGTGGGATCGTAGTAGAGAAAACCGGCCCTTATAACCTGCTGTTTCTCTTTAGTATTGGCAt                                                                  >  1:3564507/1‑71 (MQ=255)
  tCGACGAACGTGGGATCGTAGTAGAGAAAACCGGCCCTTATAACCTGCTGTTTCTCTTTAGTATTGGCATc                                                                 >  1:3450749/1‑71 (MQ=255)
    gacgaACGTGGGATCGTAGTAGAGAAAACCGGCCCTTATAACCTGCTGTTTCTCTTTAGTATTGGCATc                                                                 <  1:1918597/69‑1 (MQ=255)
    gacgaACATGGCATCGTTGTTGAGAAAACCGGTCCGTATAACCTGCTGTTCCTGTTCAGCATCGGTATCGa                                                               >  1:1369840/1‑71 (MQ=255)
          cGTGGGATCGTAGTAGAGAAAACCGGCCCTTATAACCTGCTGTTTCTCTTTAGTATTGGCATCGATaaaaa                                                         <  1:2821963/71‑2 (MQ=255)
             gggATCGTAGTAGAGAAAACCGGCCCTTATAACCTGCTGTTTCTCTTTAGTATTGGCATCGAt                                                              >  1:3267638/1‑63 (MQ=255)
             ggcaTCGTTGTTGAGAAAACCGGTCCGTATAACCTGCTGTTCCTGTTCAGCATCGGTATCGATAAGACCa                                                       >  1:3703444/4‑70 (MQ=255)
             ggcaTCGTTGTTGAGAAAACCGGTCCGTATAACCTGCTGTTCCTGTTCAGCATCGGTATCGATAAGACCa                                                       >  1:1480316/4‑70 (MQ=255)
              ggATCGTAGTAGAGAAAACCGGCCCTTATAACCTGCTGTTTCTCTTTAGTATTGGCATCGATaaaa                                                          >  1:1260951/1‑66 (MQ=255)
                 tCGTAGTAGAGAAAACCGGCCCTTATAACCTGCTGTTTCTCTTTAGTATTGGCATCGATAAAACCAAAGCa                                                  <  1:3279586/71‑1 (MQ=255)
                  cGTAGTAGAGAAAACCGGCCCTTATAACCTGCTGTTTCTCTTTAGTATTGGCATCGATAAAACCAAAGCaa                                                 >  1:4054129/1‑71 (MQ=255)
                     agtagAGAAAACCGGCCCTTATAACCTGCTGTTTCTCTTTAGTATTGGCATCGAt                                                              >  1:4042529/1‑55 (MQ=255)
                         gagaAAACCGGTCCGTATAACCTGCTGTTCCTGTTCAGCATCGGTATCGAt                                                              <  1:1978557/51‑1 (MQ=255)
                            aaaaCCGGTCCGTATAACCTGCTGTTCCTGTTCAGCATCGGTATCGAt                                                              <  1:123114/48‑1 (MQ=255)
                            aaaaCCGGTCCGTATAACCTGCTGTTCCTGTTCAGCATCGGTATCGAt                                                              <  1:1400862/48‑1 (MQ=255)
                            aaaaCCGGCCCTTATAACCTGCTGTTTCTCTTTAGTATTGGCATCGATAAAACCAAAGCAATGGGattatt                                       >  1:2056157/1‑71 (MQ=255)
                            aaaaCCGGCCCTTATAACCTGCTGTTTCTCTTTAGTATTGGCATCGATAAAACCAAAGCAATGGGattatt                                       >  1:3405240/1‑71 (MQ=255)
                            aaaaCCGGCCCTTATAACCTGCTGTTTCTCTTTAGTATTGGCATCGATAAAACCAAAGCAATGGGattat                                        >  1:1062013/1‑70 (MQ=255)
                            aaaaCCGGCCCTTATAACCTGCTGTTTCTCTTTAGTATTGGCATCGATAAAACCAAAGCAATGGGatta                                         >  1:3646824/1‑69 (MQ=255)
                             aaaCCGGTCCGTATAACCTGCTGTTCCTGTTCAGCATCGGTATCg                                                                <  1:793907/45‑1 (MQ=255)
                              aaCCGGCCCTTATAACCTGCTGTTTCTCTTTAGTATTGGCATCGATAAAACCAAAGCAATGGGATTATTg                                      >  1:3524962/1‑70 (MQ=255)
                                 cGGCCCTTATAACCTGCTGTTTCTCTTTAGTATTGGCATCGATAAAACCAAAGCAATGGGATTATTGCGTg                                  <  1:1231794/71‑1 (MQ=255)
                                     ccTTATAACCTGCTGTTTCTCTTTAGTATTGGCATCGATAAAACCAAAGCaa                                                 <  1:3204363/52‑1 (MQ=255)
                                      cTTATAACCTGCTGTTTCTCTTTAGTATTGGCATCGATaaaa                                                          >  1:1858961/1‑42 (MQ=255)
                                      cTTATAACCTGCTGTTTCTCTTTAGTATTGGCATCGATaaaa                                                          >  1:454097/1‑42 (MQ=255)
                                      cTTATAACCTGCTGTTTCTCTTTAGTATTGGCATCGATAAAACCAAAGCAATGGGat                                           <  1:528183/57‑1 (MQ=255)
                                      cTTATAACCTGCTGTTTCTCTTTAGTATTGGCATCGATAAAACCAAAGCAATGGGat                                           <  1:1104122/57‑1 (MQ=255)
                                       ttATAACCTGCTGTTTCTCTTTAGTATTGGCATCGATaa                                                            <  1:601685/39‑1 (MQ=255)
                                       ttATAACCTGCTGTTTCTCTTTAGTATTGGCATCGATAAAACCa                                                       <  1:526092/44‑1 (MQ=255)
                                       ttATAACCTGCTGTTTCTCTTTAGTATTGGCATCGATAAAACCa                                                       <  1:2296245/44‑1 (MQ=255)
                                       ttATAACCTGCTGTTTCTCTTTAGTATTGGCATCGATAAAACCAAAGCAATGGGATTATTg                                      <  1:413607/61‑1 (MQ=255)
                                       ttATAACCTGCTGTTTCTCTTTAGTATTGGCATCGATAAAACCAAAGCAATGGGATTATTg                                      <  1:2585473/61‑1 (MQ=255)
                                        tataACCTGCTGTTTCTCTTTAGTATTGGCATCGATaa                                                            >  1:3303219/1‑38 (MQ=255)
                                        tataACCTGCTGTTTCTCTTTAGTATTGGCATCGATAAAAcc                                                        <  1:1308194/42‑1 (MQ=255)
                                         ataACCTGCTGTTTCTCTTTAGTATTGGCATCGATAAAACCAAAGCAATGGGATTATTGCGTGGGTTGACg                          >  1:494207/1‑71 (MQ=255)
                                         ataACCTGCTGTTTCTCTTTAGTATTGGCATCGATAAAACCAAAGCAATGGGATTATTGCGTGGGTTGACg                          >  1:3735814/1‑71 (MQ=255)
                                         ataACCTGCTGTTTCTCTTTAGTATTGGCATCGATAAAACCAAAGCAATGGGATTATTGCGTGGGTTGACg                          >  1:634149/1‑71 (MQ=255)
                                         ataACCTGCTGTTTCTCTTTAGTATTGGCATCGATAAAACCAAAGCAATGGGATTATTGCGTGGGTTGACg                          >  1:3672948/1‑71 (MQ=255)
                                         ataACCTGCTGTTTCTCTTTAGTATTGGCATCGATAAAACCAAAGCAATGGGATTATTGCGTGGGTTGACg                          >  1:1787046/1‑71 (MQ=255)
                                         ataACCTGCTGTTTCTCTTTAGTATTGGCATCGATAAAACCAAAGCAATGGGATTATTGCGTGGGTTGACg                          >  1:1326514/1‑71 (MQ=255)
                                         ataACCTGCTGTTTCTCTTTAGTATTGGCATCGATAAAACCAAAGCAATGGGATTATTGCGTGGGTTGACg                          >  1:1081613/1‑71 (MQ=255)
                                         ataACCTGCTGTTCCTGTTCAGCATCGGTATCGATAAGACCAAAGCa                                                  >  1:2876014/1‑47 (MQ=255)
                                         ataACCTGCTGTTCCTGTTCAGCATCGGTATCGATAAGACCAAAGCa                                                  >  1:465623/1‑47 (MQ=255)
                                          taACCTGCTGTTTCTCTTTAGTATTGGCATCGATAAAACCa                                                       <  1:1272730/41‑1 (MQ=255)
                                          taACCTGCTGTTTCTCTTTAGTATTGGCATCGATAAAACCAAAGCAATGGGATTATTGCGTggg                                >  1:2022705/1‑64 (MQ=255)
                                            aCCTGCTGTTTCTCTTTAGTATTGGCATCGATAAAACCAAAGCAATGGGattat                                        >  1:2867743/1‑54 (MQ=255)
                                            aCCTGCTGTTTCTCTTTAGTATTGGCATCGATAAAACCAAAGCAATGGGATTATTGCGTggg                                >  1:4085314/1‑62 (MQ=255)
                                             ccTGCTGTTTCTCTTTAGTATTGGCATCGATAAAACCAAAGCAATGGGATTATTGCGTGGGTTGACGGAAt                      <  1:1704434/71‑1 (MQ=255)
                                              ctgctgTTTCTCTTTAGTATTGGCATCGATAAAACCAAAGCGATGGGATTATTGCGTGGGTTGACGGAAtt                     >  1:2887998/1‑71 (MQ=255)
                                                gctgTTTCTCTTTAGTATTGGCATCGATAAAACCAAAGCAATGGGat                                           <  1:2068007/47‑1 (MQ=255)
                                                gctgTTTCTCTTTAGTATTGGCATCGATAAAACCAAAGCAATGGGATTATTGCGTGGGTTGACGGa                        >  1:2560406/1‑66 (MQ=255)
                                                  tgTTTCTCTTTAGTATTGGCATCGATAAAACCAAAGCAATGGGat                                           <  1:2279967/45‑1 (MQ=255)
                                                   gTTTCTCTTTAGTATTGGCATCGATAAAACCATAGc                                                   <  1:1995760/36‑1 (MQ=255)
                                                    tttCTCTTTAGTATTGGCATCGATAAAATCAAAGc                                                   <  1:3519/35‑1 (MQ=255)
                                                    tttCTCTTTAGTATTGGCATCGATAAAACCAAAGc                                                   <  1:631527/35‑1 (MQ=255)
                                                    tttCTCTTTAGTATTGGCATCGATAAAACCAAAGc                                                   <  1:418897/35‑1 (MQ=255)
                                                    tttCTCTTTAGTATTGGCATCGATAAAACCAAAGCaa                                                 <  1:2391942/37‑1 (MQ=255)
                                                    tttCTCTTTAGTATTGGCATCGATAAAACCAAAGCAATgg                                              <  1:1260194/40‑1 (MQ=255)
                                                    tttCTCTTTAGTATTGGCATCGATAAAACCAAAGCAATGGGatt                                          <  1:1382404/44‑1 (MQ=255)
                                                    tttCTCTTTAGTATTGGCATCGATAAAACCAAAGCAATGGGATTATTGCGTGGGTTGACg                          <  1:2674286/60‑1 (MQ=255)
                                                     ttCTCTTTAGTATTGGCATCGATAAAACCAAAGCaa                                                 >  1:2877912/1‑36 (MQ=255)
                                                      tctcTTTAGTATTGGCATCGATAAAACCAAAGCAATg                                               >  1:1369305/1‑37 (MQ=255)
                                                      tctcTTTAGTATTGGCATCGATAAAACCAAAGCAATGGGatta                                         >  1:1219563/1‑43 (MQ=255)
                                                      tctcTTTAGTATTGGCATCGATAAAACCAAAGCAATGGGATTATTGCGTGGGTTGa                            >  1:513882/1‑56 (MQ=255)
                                                      tctcTTTAGTATTGGCATCGATAAAACCAAAGCAATGGGATTATTGCGTGGGTTGACGGAATTCAAACGc              >  1:3164196/1‑70 (MQ=255)
                                                      tctcTTTAGTATTGGCATCGAAAAAGCCAAAGCAATGGGATTATTGCGTGGGTTGACGGAATTCAAACGct             >  1:1624280/1‑71 (MQ=255)
                                                        tctTTAGTATTGGCATCGATAAAACCAAAGCAATGGGATTATTGCGTGGGTTGACGGAATTCAAACGctct           >  1:3296535/1‑71 (MQ=255)
                                                          tttAGTATTGGCATCGATAAAACCAAAGCAATGGGATTATTGCGTGGGTTGACGGAATTCAAACGctct           <  1:1242821/69‑1 (MQ=255)
                                                            tAGTATTGGCATCGATAAAACCAAAGCAATGGGATTATTGCGTGGGTTg                             <  1:3907176/49‑1 (MQ=255)
                                                              gTATTGGCATCGATAAAACCAAAGCAATGGGATTATTGCGTGGGTTGACGGAATTCAAACGCTCTTACg       <  1:1679967/69‑1 (MQ=255)
                                                               tATTGGCATCGATAAAACCAAAGCAATGGGATTATTGCGTGGGTTGACGGAATTCAAACGCTCtt          <  1:2789881/65‑1 (MQ=255)
                                                                 ttGGCATCGATAAAACCAAAGCAATGGGATTATTGCGTGGGTTg                             <  1:333407/44‑1 (MQ=255)
                                                                 ttGGCATCGATAAAACCAAAGCAATGGGATTATTGCGTGGGTTGACGGAATTCAAACGCTCTTACGAtctc  <  1:1350601/71‑1 (MQ=255)
                                                                 ttGGCATCGATAAAACCAAAGCAATGGGATTATTGCGTGGGTTGACGGAATTCAAACGCTCTTACGAtct   >  1:3733548/1‑70 (MQ=255)
                                                                 ttGGCATCGATAAAACCAAAGCAATGGGATTATTGCGTGGGTTGACGGAATTCAAACGCTCTTACGAtct   >  1:1277429/1‑70 (MQ=255)
                                                                 ttGGCATCGATAAAACCAAAGCAATGGGATTATTGCGTGGGTTGACGGAATTCAAACGCTCTTACGAtct   >  1:1986002/1‑70 (MQ=255)
                                                                  |                                                                     
CCTCGACGAACGTGGGATCGTAGTAGAGAAAACCGGCCCTTATAACCTGCTGTTTCTCTTTAGTATTGGCATCGATAAAACCAAAGCAATGGGATTATTGCGTGGGTTGACGGAATTCAAACGCTCTTACGATCTC  >  minE/190054‑190189

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: