Predicted mutation
evidence seq id position mutation freq annotation gene description
RA minE 2,463,067 A→G 31.7% E118E (GAA→GAG gidB → methyltransferase, glucose‑inhibited cell‑division protein

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*minE2,463,0670AG31.7% 66.9 / 42.9 85E118E (GAA→GAGgidBmethyltransferase, glucose‑inhibited cell‑division protein
Reads supporting (aligned to +/- strand):  ref base A (21/37);  new base G (13/14);  total (34/51)
Fisher's exact test for biased strand distribution p-value = 3.46e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 8.93e-02

GCCTTGGTAAACGCGTGCGTTTCCTTCGTCAGGTGCAACATGAGCTTAAACTGGAGAATATTGAACCAGTACAGAGCAGGGTAGAAGAGTTTCCTTCAGAGCCGCCATTTGATGGCGTAATT  >  minE/2463003‑2463124
                                                                |                                                         
gCCTTGGTAAACGCGTGCGTTTCCTTCGTCAGGTGCAACATGAGCTTAAACTGGAGAATATTGAGCCAg                                                       <  1:681144/69‑1 (MQ=255)
gCCTTGGTAAACGCGTGCGTTTCCTTCGTCAGGTGCAACATGAGCTTAAACTGGAGAATATTGAGCCAg                                                       <  1:2908637/69‑1 (MQ=255)
gCCTTGGTAAACGCGTGCGTTTCCTTCGTCAGGTGCAACATGAGCTTAAACTGGAGAATATTGAGCCAg                                                       <  1:2452462/69‑1 (MQ=255)
gCCTTGGTAAACGCGTGCGTTTCCTTCGTCAGGTGCAACATGAGCTTAAACTGGAGAATATTGAGCCAg                                                       <  1:587109/69‑1 (MQ=255)
gCCTTGGTAAACGCGTGCGTTTCCTTCGTCAGGTGCAACATGAGCTTAAACTGGAGAATATTGAGCCAg                                                       <  1:269245/69‑1 (MQ=255)
gCCTTGGTAAACGCGTGCGTTTCCTTCGTCAGGTGCAACATGAGCTTAAACTGGAGAATATTGAGCCAg                                                       <  1:1441426/69‑1 (MQ=255)
gCCTTGGTAAACGCGTGCGTTTCCTTCGTCAGGTGCAACATGAGCTTAAACTGGAGAATATTGAGCCAg                                                       <  1:1488395/69‑1 (MQ=255)
gCCTTGGTAAACGCGTGCGTTTCCTTCGTCAGGTGCAACATGAGCTTAAACTGGAGAATATTGAGCCAg                                                       <  1:2051605/69‑1 (MQ=255)
gCCTTGGTAAACGCGTGCGTTTCCTTCGTCAGGTGCAACATGAGCTTAAACTGGAGAATATTGAGCCAg                                                       <  1:2014761/69‑1 (MQ=255)
gCCTTGGTAAACGCGTGCGTTTCCTTCGTCAGGTGCAACATGAGCTTAAACTGGAGAATATTGAGCCAg                                                       <  1:272296/69‑1 (MQ=255)
gCCTTGGTAAACGCGTGCGTTTCCTTCGTCAGGTGCAACATGAGCTTAAACTGGAGAATATTGAGCCAg                                                       <  1:3166193/69‑1 (MQ=255)
gCCTTGGTAAACGCGTGCGTTTCCTTCGTCAGGTGCAACATGAGCTTAAACTGGAGAATATTGAGCCAg                                                       <  1:1778419/69‑1 (MQ=255)
 ccTTGGTAAACGCGTGCGTTTCCTTCGTCAGGTGCAACATGAGCTTAAACTGGAGAATATTGAGCCAg                                                       <  1:2733229/68‑1 (MQ=255)
 ccTTGGTAAACGCGTGCGTTTCCTTCGTCAGGTGCAACATGAGCTTAAACTGGAGAATATTGAGCCAg                                                       <  1:3024290/68‑1 (MQ=255)
       tAAACGCGTGCGTTTCCTTCGTCAGGTGCAACATGAGCTTAAACTGGAGAATATTGAACCAGTACag                                                  <  1:2713675/67‑1 (MQ=255)
       tAAACGCGTGCGTTTCCTTCGTCAGGTGCAACATGAGCTTAAACTGGAGAATATTGAACCAGTACag                                                  <  1:3331788/67‑1 (MQ=255)
       tAAACGCGTGCGTTTCCTTCGTCAGGTGCAACATGAGCTTAAACTGGAGAATATTGAACCAGTACag                                                  <  1:2694228/67‑1 (MQ=255)
               tgcgtTTCCTTCGTCAGGTGCAACATGAGCTTAAACTTGAGAATATTGAACCAGTACAGAGCAGGGt                                          >  1:1828543/1‑67 (MQ=255)
               tgcgtTTCCTTCGTCAGGTGCAACATGAGCTTAAACTGGAGAATATTGAACCAGTACAGAGCAGGGt                                          >  1:425871/1‑67 (MQ=255)
                 cgtTTCCTTCGTCAGGTGCAACATGAGCTTAATCTGGAGAATATTGAACCAGTACAGAGCAGGGTagaa                                      <  1:430407/69‑1 (MQ=255)
                 cgtTTCCTTCGTCAGGTGCAACATGAGCTTAAACTGGAGAATATTGAACCAGTACAGAGCAGGGTagaa                                      <  1:2997247/69‑1 (MQ=255)
                 cgtTTCCTTCGTCAGGTGCAACATGAGCTTAAACTGGAGAATATTGAACCAGTACAGAGCAGGGTagaa                                      <  1:2215538/69‑1 (MQ=255)
                 cgtTTCCTTCGTCAGGTGCAACATGAGCTTAAACTGGAGAATATTGAACCAGTACAGAGCAGGGTagaa                                      <  1:3037330/69‑1 (MQ=255)
                 cgtTTCCTTCGTCAGGTGCAACATGAGCTTAAACTGGAGAATATTGAACCAGTACAGAGCAGGGTagaa                                      <  1:1697254/69‑1 (MQ=255)
                 cgtTTCCTTCGTCAGGTGCAACATGAGCTTAAACTGGAGAATATTGAACCAGTACAGAGCAGGGTagaa                                      <  1:654394/69‑1 (MQ=255)
                 cgtTTCCTTCGTCAGGTGCAACATGAGCTTAAACTGGAGAATATTGAACCAGTACAGAGCAGGGTagaa                                      <  1:2690220/69‑1 (MQ=255)
                     tCCTTCGTCAGGTGCAACATGAGCTTAAACTGGAGAATATTGAACCAGTACAGAGCAGGGTAGAAGAGt                                  <  1:2557014/69‑1 (MQ=255)
                     tCCTTCGTCAGGTGCAACATGAGCTTAAACTGGAGAATATTGAACCAGTACAGAGCAGGGTAGAAGAGt                                  >  1:3134417/1‑69 (MQ=255)
                     tCCTTCGTCAGGTGCAACATGAGCTTAAACTGGAGAATATTGAACCAGTACAGAGCAGGGTAGAAGAGt                                  <  1:374755/69‑1 (MQ=255)
                     tCCTTCGTCAGGTGCAACATGAGCTTAAACTGGAGAATATTGAACCAGTACAGAGCAGGGTAGAAGAGt                                  >  1:1725900/1‑69 (MQ=255)
                     tCCTTCGTCAGGTGCAACATGAGCTTAAACTGGAGAATATTGAACCAGTACAGAGCAGGGTAGAAGAGt                                  <  1:562182/69‑1 (MQ=255)
                       cTTCGTCAGGTGCAACATTAGCTTAAACTGGAGAATATTGAACCAGTACAGAGCAGGGTAGAAGAGtt                                 <  1:604101/68‑1 (MQ=255)
                       cTTCGTCAGGTGCAACATGAGCTTAAACTGGAGAATATTGAACCAGTACAGAGCAGGGTAGAAGAGtt                                 <  1:796520/68‑1 (MQ=255)
                       cTTCGTCAGGTGCAACATGAGCTTAAACTGGAGAATATTGAACCAGTACAGAGCAGGGTAGAAGAGtt                                 <  1:109096/68‑1 (MQ=255)
                       cTTCGTCAGGTGCAACATGAGCTTAAACTGGAGAATATTGAACCAGTACAGAGCAGGGTAGAAGAGtt                                 <  1:1903988/68‑1 (MQ=255)
                         tCGTCAGGTGCAACATGAGCTTAAACTGGAGAATATTGAACCAGTACAGAGCAGGGTAGAAGAGttta                               <  1:1378129/68‑2 (MQ=255)
                         tCGTCAGGTGCAACATGAGCTTAAACTGGAGAATATTGAACCAGTACAGAGCAGGGTAGAAGAGTTTc                               <  1:828490/68‑1 (MQ=255)
                         tCGTCAGGTGCAACATGAGCTTAAACTGGAGAATATTGAACCAGTACAGAGCAGGGTAGAAGAGTTTc                               <  1:1191164/68‑1 (MQ=255)
                         tCGTCAGGTGCAACATGAGCTTAAACTGGAGAATATTGAACCAGTACAGAGCAGGGTAGAAGAGTTTc                               <  1:474738/68‑1 (MQ=255)
                         tCGTCAGGTGCAACATGAGCTTAAACTGGAGAATATTGAACCAGTACAGAGCAGGGTAGAAGAGTTTc                               <  1:301678/68‑1 (MQ=255)
                         tCGTCAGGTGCAACATGAGCTTAAACTGGAGAATATTGAACCAGTACAGAGCAGGGTAGAAGAGTTTc                               <  1:2722703/68‑1 (MQ=255)
                          cGTCAGGTGCAACATGAGCTTAAACTGGAGAATATTGAACCAGTACAGAGCAGGGTAGAAGAGTTTc                               <  1:913854/67‑1 (MQ=255)
                                    aaCATGAGCTTTTACTGGAGAATATTGAACCAGTACAGAGCAGGGTAGAAGAGTTTc                               <  1:3256134/57‑1 (MQ=255)
                                       aTGAGCTTAAACTGGAGAATATTGAACCAGTACAGAGCAGGGTAGAAGAGTTTCCTTCAGAGCCTCCAt                <  1:1204602/69‑1 (MQ=255)
                                         gAGCTTAAACTGGAGAATATTGAACCAGTACAGAGCAGGGTAGAAGAGTTTCCTTCAGAGCCGCCAttt              <  1:1952438/69‑1 (MQ=255)
                                          aGCTTAAACTGGAGAATATTGAACCAGTACAGAGCAGGGTAGAAGAGTTTCCTTCAGAGCTGCCAttt              <  1:1996856/68‑1 (MQ=255)
                                          aGCTTAAACTGGAGAATATTGAACCAGTACAGAGCAGGGTAGAAGAGTTTCCTTCAGAGCCGCCAttt              <  1:677255/68‑1 (MQ=255)
                                          aGCTTAAACTGGAGAATATTGAACCAGTACAGAGCAGGGTAGAAGAGTTTCCTTCAGAGCCGCCAttt              <  1:3144474/68‑1 (MQ=255)
                                          aGCTTAAACTGGAGAATATTGAACCAGTACAGAGCAGGGTAGAAGAGTTTCCTTCAGAGCCGCCATTTg             >  1:2456657/1‑69 (MQ=255)
                                          aGCTTAAACTGGAGAATATTGAACCAGTACAGAGCAGGGTAGAAGAGTTTCCTTCAGAGCCGCCATTTg             >  1:211249/1‑69 (MQ=255)
                                                     gagaATATTGAACCAGTACAGAGCAGGGTAGAAGTGTTTCCTTCAGAGCCGCCATTTGATGGCGTAAtt  >  1:3150063/1‑69 (MQ=255)
                                                     gagaATATTGAACCAGTACAGAGCAGGGTAGAAGAGTTTCCTTCAGAGCCGCCATTTGATGGCGTAAtt  >  1:338224/1‑69 (MQ=255)
                                                     gagaATATTGAACCAGTACAGAGCAGGGTAGAAGAGTTTCCTTCAGAGCCGCCATTTGATGGCGTAAtt  >  1:2585538/1‑69 (MQ=255)
                                                     gagaATATTGAACCAGTACAGAGCAGGGTAGAAGAGTTTCCTTCAGAGCCGCCATTTGATGGCGTAAtt  >  1:3205182/1‑69 (MQ=255)
                                                     gagaATATTGAACCAGTACAGAGCAGGGTAGAAGAGTTTCCTTCAGAGCCGCCATTTGATGGCGTAAtt  >  1:225189/1‑69 (MQ=255)
                                                     gagaATATTGAACCAGTACAGAGCAGGGTAGAAGAGTTTCCTTCAGAGCCGCCATTTGATGGCGTAAtt  >  1:251739/1‑69 (MQ=255)
                                                     gagaATATTGAACCAGTACAGAGCAGGGTAGAAGAGTTTCCTTCAGAGCCGCCATTTGATGGCGTAAt   >  1:2537760/1‑68 (MQ=255)
                                                     gagaATATTGAACCAGTACAGAGCAGGGTAGAAGAGTTTCCTTCAGAGCCGCCATTTGATGGCGTAAt   >  1:2792058/1‑68 (MQ=255)
                                                          tatTGAACCAGTACAGAGCAGGGTAGAAGAGTTTCCTTCAGAGCCGCCATTTGATGGCGt      <  1:3092872/60‑1 (MQ=255)
                                                          tatTGAACCAGTACAGAGCAGGGTAGAAGAGTTTCCTTCAGAGCCGCCATTTGATGGCGt      <  1:2064552/60‑1 (MQ=255)
                                                          tatTGAACCAGTACAGAGCAGGGTAGAAGAGTTTCCTTCAGAGCCGCCATTTGATGGCGt      <  1:3191663/60‑1 (MQ=255)
                                                          tatTGAACCAGTACAGAGCAGGGTAGAAGAGTTTCCTTCAGAGCCGCCATTTGATGGCGt      <  1:2175334/60‑1 (MQ=255)
                                                          tatTGAACCAGTACAGAGCAGGGTAGAAGAGTTTCCTTCAGAGCCGCCATTTGATGGCGt      <  1:858497/60‑1 (MQ=255)
                                                          tatTGAACCAGTACAGAGCAGGGTAGAAGAGTTTCCTTCAGAGCCGCCATTTGATGGCGt      <  1:263317/60‑1 (MQ=255)
                                                          tatTGAACCAGTACAGAGCAGGGTAGAAGAGTTTCCTTCAGAGCCGCCATTTGATGGCGt      <  1:164842/60‑1 (MQ=255)
                                                            ttGAGCCAGTACATAGTAGGGTAGAAGAGTTTCCGTCAGAGCCGCCATTTGa            >  1:3262625/1‑52 (MQ=255)
                                                            ttGAGCCAGTACAGAGTAGGGTAGAAGAGTTTCCGTCAGAGCCGCCATTTGa            >  1:1136190/1‑52 (MQ=255)
                                                            ttGAGCCAGTACAGAGTAGGGTAGAAGAGTTTCCGTCAGAGCCGCCATTTGa            >  1:2515282/1‑52 (MQ=255)
                                                            ttGAGCCAGTACAGAGTAGGGTAGAAGAGTTTCCGTCAGAGCCGCCATTTGa            >  1:564840/1‑52 (MQ=255)
                                                            ttGAGCCAGTACAGAGTAGGGTAGAAGAGTTTCCGTCAGAGCCGCCATTTGa            >  1:1415401/1‑52 (MQ=255)
                                                            ttGAGCCAGTACAGAGTAGGGTAGAAGAGTTTCCGTCAGAGCCGCCATTTGa            >  1:2448611/1‑52 (MQ=255)
                                                            ttGAGCCAGTACAGAGTAGGGTAGAAGAGTTTCCGTCAGAGCCGCCATTTGa            >  1:434724/1‑52 (MQ=255)
                                                            ttGAGCCAGTACAGAGTAGGGTAGAAGAGTTTCCGTCAGAGCCGCCATTTGa            >  1:1694010/1‑52 (MQ=255)
                                                            ttGAGCCAGTACAGAGTAGGGTAGAAGAGTTTCCGTCAGAGCCGCCATTTGa            >  1:38053/1‑52 (MQ=255)
                                                            ttGAGCCAGTACAGAGTAGGGTAGAAGAGTTTCCGTCAGAGCCGCCATTTGa            >  1:1850322/1‑52 (MQ=255)
                                                            ttGAGCCAGTACAGAGTAGGGTAGAAGAGTTTCCGTCAGAGCCGCCATTTGa            >  1:2243143/1‑52 (MQ=255)
                                                            ttGAGCCAGTACAGAGTAGGGTAGAAGAGTTTCCGTCAGAGCCGCCATTTGa            >  1:3038029/1‑52 (MQ=255)
                                                            ttGAGCCAGTACAGAGTAGGGTAGAAGAGTTTCCGTCAGAGCCGCCATTTGa            >  1:2176773/1‑52 (MQ=255)
                                                            ttGAACCAGTACAGAGCAGGGTAGAAGAGTTTCCtt                            >  1:3376277/1‑36 (MQ=255)
                                                            ttGAACCAGTACAGAGCAGGGTAGAAGAGTTTCCtt                            >  1:2658492/1‑36 (MQ=255)
                                                            ttGAACCAGTACAGAGCAGGGTAGAAGAGTTTCCtt                            >  1:329701/1‑36 (MQ=255)
                                                            ttGAACCAGTACAGAGCAGGGTAGAAGAGTTTCCtt                            >  1:2382346/1‑36 (MQ=255)
                                                            ttGAACCAGTACAGAGCAGGGTAGAAGAGTTTCCtt                            >  1:1014411/1‑36 (MQ=255)
                                                            ttGAACCAGTACAGAGCAGGGTAGAAGAGTTTCCtt                            >  1:2845344/1‑36 (MQ=255)
                                                            ttGAACCAGTACAGAGCAGGGTAGAAGAGTTTCCtt                            >  1:2182294/1‑36 (MQ=255)
                                                                |                                                         
GCCTTGGTAAACGCGTGCGTTTCCTTCGTCAGGTGCAACATGAGCTTAAACTGGAGAATATTGAACCAGTACAGAGCAGGGTAGAAGAGTTTCCTTCAGAGCCGCCATTTGATGGCGTAATT  >  minE/2463003‑2463124

Alignment Legend
Aligned base mismatch/match (shaded by quality score): ATCG/ATCG < 3 ≤ ATCG/ATCG < 21 ≤ ATCG/ATCG < 32 ≤ ATCG/ATCG < 36 ≤ ATCG/ATCG
Unaligned base: atcg    Masked matching base: atcg    Alignment gap:     Deleted base: