Predicted mutation
evidence seq id position mutation freq annotation gene description
RA minE 48,236 C→T 46.8% S222N (AGC→AAC)  apaH ← diadenosine tetraphosphatase

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*minE48,2360CT46.8% 10.2 / 65.1 79S222N (AGC→AAC) apaHdiadenosine tetraphosphatase
Reads supporting (aligned to +/- strand):  ref base C (17/25);  new base T (22/15);  total (39/40)
Fisher's exact test for biased strand distribution p-value = 1.17e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 4.64e-01

CAGCGCGTATATACCTTCCGGCGTACCTTTGCCCTCCAGCGATGCCCAGTGACCAAAGGCGATGCTGTATTCTTCAGCGACAGGGCCAGGAATCGCAAACCACGGTTTCAGTGGGGCAGGGGCCTCTTCCG  >  minE/48172‑48302
                                                                |                                                                  
cAGCGCGTATATACCTTCCGGCGTACCTTTGCCCTCCAGCGATGCCCAGTGACCAAAGGCGATGCTGt                                                                 >  1:1135835/1‑68 (MQ=255)
cAGCGCGTATATACCTTCCGGCGTACCTTTGCCCTCCAGCGATGCCCAGTGACCAAAGGCGATGCTGTa                                                                >  1:1918757/1‑69 (MQ=255)
cAGCGCGTATATACCTTCCGGCGTACCTTTGCCCTCCAGCGATGCCCAGTGACCAAAGGCGATGCTGTa                                                                >  1:2042071/1‑69 (MQ=255)
       tatataCCTTCCGGCGTACCTTTGCCCTCCAGCGATGCCCAGTGACCAAAGGCGATGCTGTATTCTTCa                                                         <  1:966704/69‑1 (MQ=255)
       tatataCCTTCCGGCGTACCTTTGCCCTCCAGCGATGCCCAGTGACCAAAGGCGATGCTGTATTCTTCa                                                         <  1:896332/69‑1 (MQ=255)
       tatataCCTTCCGGCGTACCTTTGCCCTCCAGCGATGCCCAGTGACCAAAGGCGATGCTGTATTCTTCa                                                         <  1:3406277/69‑1 (MQ=255)
       tatataCCTTCCGGCGTACCTTTGCCCTCCAGCGATGCCCAGTGACCAAAGGCGATGCTGTATTCTTCa                                                         <  1:2959823/69‑1 (MQ=255)
         tataCCTTCCGGCGTACCTTTGCCCTCCAGCGATGCCCAGTGACCAAAGGCGATGCTg                                                                  <  1:3243518/58‑1 (MQ=255)
         tataCCTTCCGGCGTACCTTTGCCCTCCAGCGATGCCCAGTGACCAAAGGCGATGCTg                                                                  <  1:2020744/58‑1 (MQ=255)
         tataCCTTCCGGCGTACCTTTGCCCTCCAGCGATGCCCAGTGACCAAAGGCGATGCTg                                                                  <  1:2444032/58‑1 (MQ=255)
         tataCCTTCCGGCGTACCTTTGCCCTCCAGCGATGCCCAGTGACCAAAGGCGATGCTg                                                                  <  1:2526908/58‑1 (MQ=255)
         tataCCTTCCGGCGTACCTTTGCCCTCCAGCGATGCCCAGTGACCAAAGGCGATGCTg                                                                  <  1:2600932/58‑1 (MQ=255)
         tataCCTTCCGGCGTACCTTTGCCCTCCAGCGATGCCCAGTGACCAAAGGCGATGCTg                                                                  <  1:2853037/58‑1 (MQ=255)
         tataCCTTCCGGCGTACCTTTGCCCTCCAGCGATGCCCAGTGACCAAAGGCGATGCTg                                                                  <  1:3355499/58‑1 (MQ=255)
         tataCCTTCCGGCGTACCTTTGCCCTCCAGCGATGCCCAGTGACCAAAGGCGATGCTg                                                                  <  1:1430586/58‑1 (MQ=255)
         tataCCTTCCGGCGTACCTTTGCCCTCCAGCGATGCCCAGTGACCAAAGGCGATGCTg                                                                  <  1:370181/58‑1 (MQ=255)
         tataCCTTCCGGCGTACCTTTGCCCTCCAGCGATGCCCAGTGACCAAAGGCGATGCTg                                                                  <  1:1255120/58‑1 (MQ=255)
         tataCCTTCCGGCGTACCTTTGCCCTCCAGCGATGCCCAGTGACCAAAGGCGATGCTg                                                                  <  1:794952/58‑1 (MQ=255)
          ataCCTTCCGGCGTACCTTTGCCTTCCAGCGATGCCCAGTGACCAAAGGCGATGTTGTActctt                                                           <  1:234302/64‑1 (MQ=255)
          ataCCTTCCGGCGTACCTTTGCCTTCCAGCGATGCCCAGTGACCAAAGGCGATGTTGTActctt                                                           <  1:3328017/64‑1 (MQ=255)
          ataCCTTCCGGCGTACCTTTGCCTTCCAGCGATGCCCAGTGACCAAAGGCGATGTTGTActctt                                                           <  1:2804061/64‑1 (MQ=255)
          ataCCTTCCGGCGTACCTTTGCCTTCCAGCGATGCCCAGTGACCAAAGGCGATGTTGTActctt                                                           <  1:2771941/64‑1 (MQ=255)
          ataCCTTCCGGCGTACCTTTGCCTTCCAGCGATGCCCAGTGACCAAAGGCGATGTTGTActctt                                                           <  1:2383250/64‑1 (MQ=255)
          ataCCTTCCGGCGTACCTTTGCCTTCCAGCGATGCCCAGTGACCAAAGGCGATGTTGTActctt                                                           <  1:1007013/64‑1 (MQ=255)
          ataCCTTCCGGCGTACCTTTGCCTTCCAGCGATGCCCAGTGACCAAAGGCGATGTTGTActctt                                                           <  1:829272/64‑1 (MQ=255)
          ataCCTTCCGGCGTACCTTTGCCTTCCAGCGATGCCCAGTGACCAAAGGCGATGTTGTActctt                                                           <  1:904914/64‑1 (MQ=255)
          ataCCTTCCGGCGTACCTTTGCCTTCCAGCGATGCCCAGTGACCAAAGGCGATGTTGTActctt                                                           <  1:1836437/64‑1 (MQ=255)
          ataCCTTCCGGCGTACCTTTGCCTTCCAGCGATGCCCAGTGACCAAAGGCGATGTTGTActctt                                                           <  1:1788736/64‑1 (MQ=255)
          ataCCTTCCGGCGTACCTTTGCCTTCCAGCGATGCCCAGTGACCAAAGGCGATGTTGTActctt                                                           <  1:1596732/64‑1 (MQ=255)
          ataCCTTCCGGCGTACCTTTGCCTTCCAGCGATGCCCAGTGACCAAAGGCGATGTTGTActctt                                                           <  1:1447196/64‑1 (MQ=255)
          ataCCTTCCGGCGTACCTTTGCCTTCCAGCGATGCCCAGTGACCAAAGGCGATGTTGTActctt                                                           <  1:1457016/64‑1 (MQ=255)
           taCCTTCCGGCGTACCTTTGCCTTCCAGCGATGCCCAGTGACCAAAGGCGATGTTGTActctt                                                           <  1:86557/63‑1 (MQ=255)
           taCCTTCCGGCGTACCTTTGCCTTCCAGCGATGCCCAGTGACCAAAGGCGATGTTGTActctt                                                           <  1:1499381/63‑1 (MQ=255)
           taCCTTACGGCGTACCTTTGCCCTCCAGCGATGCCCAGTGACCAAAGGCGATGCTGTATTCTTCAGc                                                       >  1:1965055/1‑67 (MQ=255)
                                  tCCAGCGATGGCCAGTGACCAAAGGCGATGTTGTACTCTTCAGCGACgggg                                                >  1:944285/1‑51 (MQ=255)
                                  tCCAGCGATGCCCAGTGACCAAAGGCGATGTTGTACTCTTCAGCGACgggg                                                >  1:1631680/1‑51 (MQ=38)
                                  tCCAGCGATGCCCAGTGACCAAAGGCGATGTTGTACTCTTCAGCGACgggg                                                >  1:1748688/1‑51 (MQ=38)
                                  tCCAGCGATGCCCAGTGACCAAAGGCGATGTTGTACTCTTCAGCGACgggg                                                >  1:2714162/1‑51 (MQ=38)
                                  tCCAGCGATGCCCAGTGACCAAAGGCGATGTTGTACTCTTCAGCGACgggg                                                >  1:2790819/1‑51 (MQ=38)
                                  tCCAGCGATGCCCAGTGACCAAAGGCGATGTTGTACTCTTCAGCGACgggg                                                >  1:1567823/1‑51 (MQ=38)
                                  tCCAGCGATGCCCAGTGACCAAAGGCGATGTTGTACTCTTCAGCGACgggg                                                >  1:2811699/1‑51 (MQ=38)
                                  tCCAGCGATGCCCAGTGACCAAAGGCGATGTTGTACTCTTCAGCGACgggg                                                >  1:3036721/1‑51 (MQ=38)
                                  tCCAGCGATGCCCAGTGACCAAAGGCGATGTTGTACTCTTCAGCGACgggg                                                >  1:1458945/1‑51 (MQ=38)
                                  tCCAGCGATGCCCAGTGACCAAAGGCGATGTTGTACTCTTCAGCGACgggg                                                >  1:3344896/1‑51 (MQ=38)
                                  tCCAGCGATGCCCAGTGACCAAAGGCGATGTTGTACTCTTCAGCGACgggg                                                >  1:1302562/1‑51 (MQ=38)
                                  tCCAGCGATGCCCAGTGACCAAAGGCGATGTTGTACTCTTCAGCGACgggg                                                >  1:42269/1‑51 (MQ=38)
                                  tCCAGCGATGCCCAGTGACCAAAGGCGATGTTGTACTCTTCAGCGACgggg                                                >  1:650285/1‑51 (MQ=38)
                                  tCCAGCGATGCCCAGTGACCAAAGGCGATGTTGTACTCTTCAGCGACgggg                                                >  1:886419/1‑51 (MQ=38)
                                  tCCAGCGATGCCCAGTGACCAAAGGCGATGTTGTACTCTTCAGCGACgggg                                                >  1:2562890/1‑51 (MQ=38)
                                  tCCAGCGATGCCCAGTGACCAAAGGCGATGTTGTACTCTTCAGCGACgggg                                                >  1:1639453/1‑51 (MQ=38)
                                  tCCAGCGATGCCCAGTGACCAAAGGCGATGTTGTACTCTTCAGCGACgggg                                                >  1:2422598/1‑51 (MQ=38)
                                  tCCAGCGATGCCCAGTGACCAAAGGCGATGTTGTACTCTTCAGCGACgggg                                                >  1:1704622/1‑51 (MQ=38)
                                  tCCAGCGATGCCCAGTGACCAAAGGCGATGTTGTACTCTTCAGCGACgggg                                                >  1:1677291/1‑51 (MQ=38)
                                  tCCAGCGATGCCCAGTGACCAAAGGCGATGTTGTACTCTTCAGCGACgggg                                                >  1:2115234/1‑51 (MQ=38)
                                  tCCAGCGATGCCCAGTGACCAAAGGCGATGTTGTACTCTTCAGCGACgggg                                                >  1:2580052/1‑51 (MQ=38)
                                  tCCAGCGATGCCCAGTGACCAAAGGCGATGTTGTACTCTTCAGCGACgggg                                                >  1:2283186/1‑51 (MQ=38)
                                   ccAGCGCTGCCCAGTGACCAAAGGCGATGCTGTATTCTTCAGCGACAGGGCCAGGAATCGCAAACCACg                             >  1:1058/1‑69 (MQ=255)
                                   ccAGCGATGCCCAGTGACCAAAGGCGATGCTGTATTCTTCAGCGACAGGGCCAGGAATCGCAAACCACg                             >  1:2180961/1‑69 (MQ=255)
                                   ccAGCGATGCCCAGTGACCAAAGGCGATGCTGTATTCTTCAGCGACAGGGCCAGGAATCGCAAACCACg                             >  1:1938307/1‑69 (MQ=255)
                                   ccAGCGATGCCCAGTGACCAAAGGCGATGCTGTATTCTTCAGCGACAGGGCCAGGAATCGCAAACCACg                             >  1:1171706/1‑69 (MQ=255)
                                   ccAGCGATGCCCAGTGACCAAAGGCGATGCTGTATTCTTCAGCGACAGGGCCAGGAATCGCAAACCACg                             >  1:516598/1‑69 (MQ=255)
                                   ccAGCGATGCCCAGTGACCAAAGGCGATGCTGTATTCTTCAGCGACAGGGCCAGGAATCGCAAACCACg                             >  1:1648728/1‑69 (MQ=255)
                                   ccAGCGATGCCCAGTGACCAAAGGCGATGCTGTATTCTTCAGCGACAGGGCCAGGAATCGCAAACCACg                             >  1:2398717/1‑69 (MQ=255)
                                        gATGCCCAGTGACCAAAGGCGATGCTGTATTCTTCAGCGACAGGGCCAGGAATCGCAAACCACGGTTTc                        <  1:3336682/69‑1 (MQ=255)
                                        gATGCCCAGTGACCAAAGGCGATGCTGTATTCTTCAGCGACAGGGCCAGGAATCGCAAACCACGGTTTc                        <  1:1635084/69‑1 (MQ=255)
                                        gATGCCCAGTGACCAAAGGCGATGCTGTATTCTTCAGCGACAGGGCCAGGAATCGCAAACCACGGTTTc                        <  1:1393088/69‑1 (MQ=255)
                                        gATGCCCAGTGACCAAAGGCGATGCTGTATTCTTCAGCGACAGGGCCAGGAATCGCAAACCACGGTTTc                        <  1:307255/69‑1 (MQ=255)
                                        gATGCCCAGTGACCAAAGGCGATGCTGTATTCTTCAGCGACAGGGCCAGGAATCGCAAACCACGGTTTc                        <  1:390271/69‑1 (MQ=255)
                                        gATGCCCAGTGACCAAAGGCGATGCTGTATTCTTCAGCGACAGGGCCAGGAATCGCAAACCACGGTTTc                        <  1:913657/69‑1 (MQ=255)
                                        gATGCCCAGTGACCAAAGGCGATGCTGTATTCTTCAGCGACAGGGCCAGGAATCGCAAACCACGGTTTc                        <  1:1475502/69‑1 (MQ=255)
                                        gATGCCCAGTGACCAAAGGCGATGCTGTATTCTTCAGCGACAGGGCCAGGAATCGCAAACCACGGTTTc                        <  1:2919386/69‑1 (MQ=255)
                                        gATGCCCAGTGACCAAAGGCGATGCTGTATTCTTCAGCGACAGGGCCAGGAATCGCAAACCACGGTTTc                        <  1:1987867/69‑1 (MQ=255)
                                                   aCCAAAGGCGATGCTGTATTCTTCAGCGACAGGGCCAGGAATCGCAAACCACGGTTTCAGTGGGGCAgg             >  1:1015809/1‑69 (MQ=255)
                                                   aCCAAAGGCGATGCTGTATTCTTCAGCGACAGGGCCAGGAATCGCAAACCACGGTTTCAGTGGGGCAgg             >  1:164949/1‑69 (MQ=255)
                                                         ggCGATGCTGTATTCTTCAGCGACAGGGCCAGGAATCGCAAACCACGGTTTCAGTGGGGCAGGGGc          <  1:2350760/66‑1 (MQ=255)
                                                         ggCGATGCTGTATTCTTCAGCGACAGGGCCAGGAATCGCAAACCACGGTTTCAGTGGGGCAGGGGc          >  1:2339667/1‑66 (MQ=255)
                                                               gCTGTATTCTTCAGCGACAGGGCCAGGAATCGCAAACCACGGTTTCAGTGGGGCAGGGGCCTCTTCCg  >  1:1946255/1‑68 (MQ=255)
                                                               gCTGTATTCTTCAGCGACAGGGCCAGGAATCGCAAACCACGGTTTCAGTGGGGCAGGGGCCTCTTCCg  >  1:1022228/1‑68 (MQ=255)
                                                               gCTGTATTCTTCAGCGACAGGGCCAGGAATCGCAAACCACGGTTTCAGTGGGGCAGGGGCCTCTTCCg  >  1:1018326/1‑68 (MQ=255)
                                                                |                                                                  
CAGCGCGTATATACCTTCCGGCGTACCTTTGCCCTCCAGCGATGCCCAGTGACCAAAGGCGATGCTGTATTCTTCAGCGACAGGGCCAGGAATCGCAAACCACGGTTTCAGTGGGGCAGGGGCCTCTTCCG  >  minE/48172‑48302

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