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

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*minE190,1140TC14.3% 195.5 / 22.9 77F537F (TTT→TTCldcClysine decarboxylase 2, constitutive
Reads supporting (aligned to +/- strand):  ref base T (37/29);  new base C (6/5);  total (43/34)
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

GCAAAATTCCTCGACGAACGTGGGATCGTAGTAGAGAAAACCGGCCCTTATAACCTGCTGTTTCTCTTTAGTATTGGCATCGATAAAACCAAAGCAATGGGATTATTGCGTGGGTTGACGGAATTCAAACGCTCT  >  minE/190046‑190180
                                                                    |                                                                  
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: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: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:3672948/1‑71 (MQ=255)
                                                 ataACCTGCTGTTTCTCTTTAGTATTGGCATCGATAAAACCAAAGCAATGGGATTATTGCGTGGGTTGACg                 >  1:634149/1‑71 (MQ=255)
                                                 ataACCTGCTGTTTCTCTTTAGTATTGGCATCGATAAAACCAAAGCAATGGGATTATTGCGTGGGTTGACg                 >  1:494207/1‑71 (MQ=255)
                                                 ataACCTGCTGTTTCTCTTTAGTATTGGCATCGATAAAACCAAAGCAATGGGATTATTGCGTGGGTTGACg                 >  1:3735814/1‑71 (MQ=255)
                                                 ataACCTGCTGTTTCTCTTTAGTATTGGCATCGATAAAACCAAAGCAATGGGATTATTGCGTGGGTTGACg                 >  1:1081613/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)
                                                                  tttAGTATTGGCATCGATAAAACCAAAGCAATGGGATTATTGCGTGGGTTGACGGAATTCAAACGctct  <  1:1242821/69‑1 (MQ=255)
                                                                    tAGTATTGGCATCGATAAAACCAAAGCAATGGGATTATTGCGTGGGTTg                    <  1:3907176/49‑1 (MQ=255)
                                                                    |                                                                  
GCAAAATTCCTCGACGAACGTGGGATCGTAGTAGAGAAAACCGGCCCTTATAACCTGCTGTTTCTCTTTAGTATTGGCATCGATAAAACCAAAGCAATGGGATTATTGCGTGGGTTGACGGAATTCAAACGCTCT  >  minE/190046‑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: