Predicted mutation
evidence seq id position mutation freq annotation gene description
RA minE 190,117 T→C 14.6% S538S (AGT→AGC ldcC → lysine decarboxylase 2, constitutive

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*minE190,1170TC14.6% 189.0 / 22.4 75S538S (AGT→AGCldcClysine decarboxylase 2, constitutive
Reads supporting (aligned to +/- strand):  ref base T (36/28);  new base C (6/5);  total (42/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 = 1.00e+00

AATTCCTCGACGAACGTGGGATCGTAGTAGAGAAAACCGGCCCTTATAACCTGCTGTTTCTCTTTAGTATTGGCATCGATAAAACCAAAGCAATGGGATTATTGCGTGGGTTGACGGAATTCAAACGCTCTTACG  >  minE/190050‑190184
                                                                   |                                                                   
aaTACCTCGACGAACATGGCATCGTTGTTGAGAAAACCGGTCCGTATAACCTGCTGTTCCTGTTCAGCAtc                                                                  >  1:989448/1‑70 (MQ=255)
 ataccTCGACGAACATGGCATCGTTGTTGAGAAAACCGGTCCGTATAACCTGCTGTTCCTGTTCAGCAt                                                                   <  1:1156892/66‑1 (MQ=255)
    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)
               gTGGGATCGTAGTAGAGAAAACCGGCCCTTATAACCTGCTGTTTCTCTTTAGt                                                                     <  1:2450675/53‑1 (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)
                             gagaAAACCGGTCCGTATAACCTGCTGTTCCTGTTCAGCATCGGTATCGAt                                                         <  1:1978557/51‑1 (MQ=255)
                                aaaaCCGGTCCGTATAACCTGCTGTTCCTGTTCAGCATCGGTATCGAt                                                         <  1:1400862/48‑1 (MQ=255)
                                aaaaCCGGTCCGTATAACCTGCTGTTCCTGTTCAGCATCGGTATCGAt                                                         <  1:123114/48‑1 (MQ=255)
                                aaaaCCGGCCCTTATAACCTGCTGTTTCTCTTTAGTATTGGCATCGATAAAACCAAAGCAATGGGattatt                                  >  1:3405240/1‑71 (MQ=255)
                                aaaaCCGGCCCTTATAACCTGCTGTTTCTCTTTAGTATTGGCATCGATAAAACCAAAGCAATGGGattatt                                  >  1:2056157/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:454097/1‑42 (MQ=255)
                                          cTTATAACCTGCTGTTTCTCTTTAGTATTGGCATCGATaaaa                                                     >  1:1858961/1‑42 (MQ=255)
                                          cTTATAACCTGCTGTTTCTCTTTAGTATTGGCATCGATAAAACCAAAGCAATGGGat                                      <  1:1104122/57‑1 (MQ=255)
                                          cTTATAACCTGCTGTTTCTCTTTAGTATTGGCATCGATAAAACCAAAGCAATGGGat                                      <  1:528183/57‑1 (MQ=255)
                                           ttATAACCTGCTGTTTCTCTTTAGTATTGGCATCGATaa                                                       <  1:601685/39‑1 (MQ=255)
                                           ttATAACCTGCTGTTTCTCTTTAGTATTGGCATCGATAAAACCa                                                  <  1:2296245/44‑1 (MQ=255)
                                           ttATAACCTGCTGTTTCTCTTTAGTATTGGCATCGATAAAACCa                                                  <  1:526092/44‑1 (MQ=255)
                                           ttATAACCTGCTGTTTCTCTTTAGTATTGGCATCGATAAAACCAAAGCAATGGGATTATTg                                 <  1:2585473/61‑1 (MQ=255)
                                           ttATAACCTGCTGTTTCTCTTTAGTATTGGCATCGATAAAACCAAAGCAATGGGATTATTg                                 <  1:413607/61‑1 (MQ=255)
                                            tataACCTGCTGTTTCTCTTTAGTATTGGCATCGATaa                                                       >  1:3303219/1‑38 (MQ=255)
                                            tataACCTGCTGTTTCTCTTTAGTATTGGCATCGATAAAAcc                                                   <  1:1308194/42‑1 (MQ=255)
                                             ataACCTGCTGTTTCTCTTTAGTATTGGCATCGATAAAACCAAAGCAATGGGATTATTGCGTGGGTTGACg                     >  1:634149/1‑71 (MQ=255)
                                             ataACCTGCTGTTTCTCTTTAGTATTGGCATCGATAAAACCAAAGCAATGGGATTATTGCGTGGGTTGACg                     >  1:494207/1‑71 (MQ=255)
                                             ataACCTGCTGTTTCTCTTTAGTATTGGCATCGATAAAACCAAAGCAATGGGATTATTGCGTGGGTTGACg                     >  1:3672948/1‑71 (MQ=255)
                                             ataACCTGCTGTTTCTCTTTAGTATTGGCATCGATAAAACCAAAGCAATGGGATTATTGCGTGGGTTGACg                     >  1:1081613/1‑71 (MQ=255)
                                             ataACCTGCTGTTTCTCTTTAGTATTGGCATCGATAAAACCAAAGCAATGGGATTATTGCGTGGGTTGACg                     >  1:3735814/1‑71 (MQ=255)
                                             ataACCTGCTGTTTCTCTTTAGTATTGGCATCGATAAAACCAAAGCAATGGGATTATTGCGTGGGTTGACg                     >  1:1326514/1‑71 (MQ=255)
                                             ataACCTGCTGTTTCTCTTTAGTATTGGCATCGATAAAACCAAAGCAATGGGATTATTGCGTGGGTTGACg                     >  1:1787046/1‑71 (MQ=255)
                                             ataACCTGCTGTTCCTGTTCAGCATCGGTATCGATAAGACCAAAGCa                                             >  1:465623/1‑47 (MQ=255)
                                             ataACCTGCTGTTCCTGTTCAGCATCGGTATCGATAAGACCAAAGCa                                             >  1:2876014/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)
                                                                   |                                                                   
AATTCCTCGACGAACGTGGGATCGTAGTAGAGAAAACCGGCCCTTATAACCTGCTGTTTCTCTTTAGTATTGGCATCGATAAAACCAAAGCAATGGGATTATTGCGTGGGTTGACGGAATTCAAACGCTCTTACG  >  minE/190050‑190184

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: