Predicted mutation
evidence seq id position mutation freq annotation gene description
RA minE 190,108 T→C 15.9% F535F (TTT→TTC ldcC → lysine decarboxylase 2, constitutive

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*minE190,1080TC15.9% 167.0 / 24.3 69F535F (TTT→TTCldcClysine decarboxylase 2, constitutive
Reads supporting (aligned to +/- strand):  ref base T (32/26);  new base C (6/5);  total (38/31)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 7.80e-01

GCTGGTAGCAAAATTCCTCGACGAACGTGGGATCGTAGTAGAGAAAACCGGCCCTTATAACCTGCTGTTTCTCTTTAGTATTGGCATCGATAAAACCAAAGCAATGGGATTATTGCGTGGGTTGACGGAATTCAAACGCT  >  minE/190039‑190178
                                                                     |                                                                      
gcTGGTAGCAAAATTCCTCGACGAACGTGGGATCGTAGTAGAGAAAACCGGCCCTTATAACCTGCTGTTtc                                                                       <  1:2875556/71‑1 (MQ=255)
    gTAGCAAAATTCCTCGACGAACGTGGGATCGTAGTAGAGAAAACCGGCCCTTATAACCTGCTGTTtctc                                                                     <  1:4049635/69‑1 (MQ=255)
    gTAGCAAAATTCCTCGACGAACGTGGGATCGTAGTAGAGAAAACCGGCCCTTATAACCTGCTGTTTCTCtt                                                                   >  1:2480545/1‑71 (MQ=255)
    gTAGCAAAATTCCTCGACGAACGTGGGATCGTAGTAGAGAAAACCGGCCCTTATAACCTGCTGTTTCTCtt                                                                   <  1:2393807/71‑1 (MQ=255)
    gTAGCAAAATTCCTCGACGAACGTGGGATCGTAGTAGAGAAAACCGGCCCTTATAACCTGCTGTTTCTCtt                                                                   <  1:1680610/71‑1 (MQ=255)
    gTAGCAAAATTCCTCGACGAACGTGGGATCGTAGTAGAGAAAACCGGCCCTTATAACCTGCTGTTTCTCtt                                                                   <  1:1466335/71‑1 (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)
              accTCGACGAACATGGCATCGTTGTTGAGAAAACCGGTCCGTATAACCTGCTGTTcctg                                                                     <  1:2034573/58‑2 (MQ=19)
               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)
                      gaACGTGGGATCGTAGTAGAGAAAACCGGCCCTTATAACCTGCTGTTTCTCTTTAg                                                                <  1:178726/56‑1 (MQ=255)
                      gaACATGGCATCGTTGTTGAGAAAACCGGTCCGTATAACCTGCTGTTcct                                                                      >  1:2411829/1‑50 (MQ=16)
                      gaACATGGCATCGTTGTTGAGAAAACCGGTCCGTATAACCTGCTGTTcct                                                                      >  1:1916936/1‑50 (MQ=16)
                         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:1480316/4‑70 (MQ=255)
                            ggcaTCGTTGTTGAGAAAACCGGTCCGTATAACCTGCTGTTCCTGTTCAGCATCGGTATCGATAAGACCa                                            >  1:3703444/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)
                                         agaAAACCGGCCCTTATAACCTGCTGTTTCTCttt                                                                  <  1:2328990/35‑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:1104122/57‑1 (MQ=255)
                                                     cTTATAACCTGCTGTTTCTCTTTAGTATTGGCATCGATAAAACCAAAGCAATGGGat                                <  1:528183/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: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:1081613/1‑71 (MQ=255)
                                                        ataACCTGCTGTTTCTCTTTAGTATTGGCATCGATAAAACCAAAGCAATGGGATTATTGCGTGGGTTGACg               >  1:494207/1‑71 (MQ=255)
                                                        ataACCTGCTGTTTCTCTTTAGTATTGGCATCGATAAAACCAAAGCAATGGGATTATTGCGTGGGTTGACg               >  1:3672948/1‑71 (MQ=255)
                                                        ataACCTGCTGTTTCTCTTTAGTATTGGCATCGATAAAACCAAAGCAATGGGATTATTGCGTGGGTTGACg               >  1:3735814/1‑71 (MQ=255)
                                                        ataACCTGCTGTTTCTCTTTAGTATTGGCATCGATAAAACCAAAGCAATGGGATTATTGCGTGGGTTGACg               >  1:634149/1‑71 (MQ=255)
                                                        ataACCTGCTGTTTCTCTTTAGTATTGGCATCGATAAAACCAAAGCAATGGGATTATTGCGTGGGTTGACg               >  1:1326514/1‑71 (MQ=255)
                                                        ataACCTGCTGTTTCTCTTTAGTATTGGCATCGATAAAACCAAAGCAATGGGATTATTGCGTGGGTTGACg               >  1:1787046/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:418897/35‑1 (MQ=255)
                                                                   tttCTCTTTAGTATTGGCATCGATAAAACCAAAGc                                        <  1:631527/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)
                                                                     |                                                                      
GCTGGTAGCAAAATTCCTCGACGAACGTGGGATCGTAGTAGAGAAAACCGGCCCTTATAACCTGCTGTTTCTCTTTAGTATTGGCATCGATAAAACCAAAGCAATGGGATTATTGCGTGGGTTGACGGAATTCAAACGCT  >  minE/190039‑190178

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: