Predicted mutation
evidence seq id position mutation freq annotation gene description
RA minE 190,111 C→G 14.3% L536L (CTC→CTG ldcC → lysine decarboxylase 2, constitutive

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*minE190,1110CG14.3% 193.4 / 22.2 77L536L (CTC→CTGldcClysine decarboxylase 2, constitutive
Reads supporting (aligned to +/- strand):  ref base C (33/33);  new base G (6/5);  total (39/38)
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.62e-01

GTAGCAAAATTCCTCGACGAACGTGGGATCGTAGTAGAGAAAACCGGCCCTTATAACCTGCTGTTTCTCTTTAGTATTGGCATCGATAAAACCAAAGCAATGGGATTATTGCGTGGGTTGACGGAATTCAAACGCTCT  >  minE/190043‑190180
                                                                    |                                                                     
gTAGCAAAATTCCTCGACGAACGTGGGATCGTAGTAGAGAAAACCGGCCCTTATAACCTGCTGTTtctc                                                                       <  1:4049635/69‑1 (MQ=255)
gTAGCAAAATTCCTCGACGAACGTGGGATCGTAGTAGAGAAAACCGGCCCTTATAACCTGCTGTTTCTCtt                                                                     <  1:1680610/71‑1 (MQ=255)
gTAGCAAAATTCCTCGACGAACGTGGGATCGTAGTAGAGAAAACCGGCCCTTATAACCTGCTGTTTCTCtt                                                                     <  1:2393807/71‑1 (MQ=255)
gTAGCAAAATTCCTCGACGAACGTGGGATCGTAGTAGAGAAAACCGGCCCTTATAACCTGCTGTTTCTCtt                                                                     >  1:2480545/1‑71 (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)
           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)
                     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)
                                    gagaAAACCGGTCCGTATAACCTGCTGTTCCTGTTCAGCATCGGTATCGAt                                                     <  1:1978557/51‑1 (MQ=255)
                                     agaAAACCGGCCCTTATAACCTGCTGTTTCTCttt                                                                    <  1:2328990/35‑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: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: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:634149/1‑71 (MQ=255)
                                                    ataACCTGCTGTTTCTCTTTAGTATTGGCATCGATAAAACCAAAGCAATGGGATTATTGCGTGGGTTGACg                 >  1:3672948/1‑71 (MQ=255)
                                                    ataACCTGCTGTTTCTCTTTAGTATTGGCATCGATAAAACCAAAGCAATGGGATTATTGCGTGGGTTGACg                 >  1:3735814/1‑71 (MQ=255)
                                                    ataACCTGCTGTTTCTCTTTAGTATTGGCATCGATAAAACCAAAGCAATGGGATTATTGCGTGGGTTGACg                 >  1:494207/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: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)
                                                                    |                                                                     
GTAGCAAAATTCCTCGACGAACGTGGGATCGTAGTAGAGAAAACCGGCCCTTATAACCTGCTGTTTCTCTTTAGTATTGGCATCGATAAAACCAAAGCAATGGGATTATTGCGTGGGTTGACGGAATTCAAACGCTCT  >  minE/190043‑190180

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: