Predicted mutation
evidence seq id position mutation freq annotation gene description
RA minE 1,449,160 C→T 49.5% R522R (CGC→CGT atoS → sensory histidine kinase in two‑component regulatory system with AtoC

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*minE1,449,1600CT49.5% 0.1 / 101.8 111R522R (CGC→CGTatoSsensory histidine kinase in two‑component regulatory system with AtoC
Reads supporting (aligned to +/- strand):  ref base C (22/34);  new base T (18/37);  total (40/71)
Fisher's exact test for biased strand distribution p-value = 5.54e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.34e-01

Only 100 of 111 total aligned reads displayed.

CAAACAGGTACTACTGAATATCCTGATCAATGCCGTCCAGGCTATCAGCGCACGAGGGAAAATTCGCATTCAAACCTGGCAATACAGCGACTCACAACAGGCCATTTCGATAGAGGACAACGGCTGTGG  >  minE/1449094‑1449222
                                                                  |                                                              
cAAACAGGTACTACTGAATATCCTGATCAATGCCGTCCAGGCTATCAGTGCACGAGGGAAAATTCGTAt                                                              >  1:3013087/1‑69 (MQ=255)
cAAACAGGTACTACTGAATATCCTGATCAATGCCGTCCAGGCTATCAGTGCACGAGGGAAAATTCGTAt                                                              >  1:1533919/1‑69 (MQ=255)
cAAACAGGTACTACTGAATATCCTGATCAATGCCGTCCAGGCTATCAGTGCACGAGGGAAAATTCGTAt                                                              >  1:2534678/1‑69 (MQ=255)
cAAACAGGTACTACTGAATATCCTGATCAATGCCGTCCAGGCTATCAGTGCACGAGGGAAAATTCGTAt                                                              >  1:1168308/1‑69 (MQ=255)
          ctactGAATATCCTGATCAATGCCGTCCAGGCTATCAGCGCACGAGGGAAAATTCGCATTCAAACCTg                                                     >  1:1153204/1‑68 (MQ=255)
            actTAATATCCTGATCAATGCCGTCCAGGCTATCAGCGCACGAGGGAAAATTCGCATTCa                                                           >  1:2904150/1‑60 (MQ=255)
            actGAATATCCTGATCAATGCCGTCCAGGCTATCAGCGCACGAGGGAAAATTCGCATTCa                                                           >  1:2408954/1‑60 (MQ=255)
            actGAATATCCTGATCAATGCCGTCCAGGCTATCAGCGCACGAGGGAAAATTCGCATTCa                                                           <  1:1582214/60‑1 (MQ=255)
            actGAATATCCTGATCAATGCCGTCCAGGCTATCAGCGCACGAGGGAAAATTCGCATTCa                                                           <  1:1610935/60‑1 (MQ=255)
            actGAATATCCTGATCAATGCCGTCCAGGCTATCAGCGCACGAGGGAAAATTCGCATTCa                                                           <  1:1712866/60‑1 (MQ=255)
            actGAATATCCTGATCAATGCCGTCCAGGCTATCAGCGCACGAGGGAAAATTCGCATTCa                                                           >  1:1469585/1‑60 (MQ=255)
            actGAATATCCTGATCAATGCCGTCCAGGCTATCAGCGCACGAGGGAAAATTCGCATTCa                                                           <  1:1946323/60‑1 (MQ=255)
            actGAATATCCTGATCAATGCCGTCCAGGCTATCAGCGCACGAGGGAAAATTCGCATTCa                                                           >  1:2173607/1‑60 (MQ=255)
            actGAATATCCTGATCAATGCCGTCCAGGCTATCAGCGCACGAGGGAAAATTCGCATTCa                                                           >  1:2390532/1‑60 (MQ=255)
            actGAATATCCTGATCAATGCCGTCCAGGCTATCAGCGCACGAGGGAAAATTCGCATTCa                                                           >  1:2410253/1‑60 (MQ=255)
            actGAATATCCTGATCAATGCCGTCCAGGCTATCAGCGCACGAGGGAAAATTCGCATTCa                                                           >  1:299007/1‑60 (MQ=255)
            actGAATATCCTGATCAATGCCGTCCAGGCTATCAGCGCACGAGGGAAAATTCGCATTCa                                                           >  1:2551468/1‑60 (MQ=255)
            actGAATATCCTGATCAATGCCGTCCAGGCTATCAGCGCACGAGGGAAAATTCGCATTCa                                                           <  1:510765/60‑1 (MQ=255)
            actGAATATCCTGATCAATGCCGTCCAGGCTATCAGCGCACGAGGGAAAATTCGCATTCa                                                           <  1:505401/60‑1 (MQ=255)
            actGAATATCCTGATCAATGCCGTCCAGGCTATCAGCGCACGAGGGAAAATTCGCATTCa                                                           >  1:2858956/1‑60 (MQ=255)
            actGAATATCCTGATCAATGCCGTCCAGGCTATCAGCGCACGAGGGAAAATTCGCATTCAAACCTGGc                                                   <  1:2886770/68‑1 (MQ=255)
            actGAATATCCTGATCAATGACGTCCAGGCTATCAGCGCACGAGGGAAAATTCGCATTCa                                                           >  1:101861/1‑60 (MQ=255)
               gAATATCCTGATCAATGCCGTCCAGGCTATCAGCGCACGAGGGAAAATTCGCATTCAAACCTGGCaa                                                 <  1:1663734/67‑1 (MQ=255)
               gAATATCCTGATCAATGCCGTCCAGGCTATCAGCGCACGAGGGAAAATTCGCATTCAAACCTGGCaa                                                 <  1:1495167/67‑1 (MQ=255)
               gAATATCCTGATCAATGCCGTCCAGGCTATCAGCGCACGAGGGAAAATTCGCATTCAAACCTGGCaa                                                 <  1:761387/67‑1 (MQ=255)
                  tatCCTGATCAATGCCGTCCAGGCTATCAGTGCACGAGGGAAAATTCGTATTCGAACCTGGCAGTata                                             >  1:12699/1‑66 (MQ=255)
                  tatCCTGATCAATGCCGTCCAGGCTATCAGTGCACGAGGGAAAATTCGTATTCGAACCTGGCAGTata                                             >  1:2469579/1‑66 (MQ=255)
                  tatCCTGATCAATGCCGTCCAGGCTATCAGTGCACGAGGGAAAATTCGTATTCGAACCTGGCAGTata                                             >  1:2590608/1‑66 (MQ=255)
                  tatCCTGATCAATGCCGTCCAGGCTATCAGTGCACGAGGGAAAATTCGTATTCGAACCTGGCAGTata                                             >  1:510310/1‑66 (MQ=255)
                  tatCCTGATCAATGCCGTCCAGGCTATCAGTGCACGAGGGAAAATTCGTATTCGAACCTGGCAGTata                                             >  1:59150/1‑66 (MQ=255)
                  tatCCTGATCAATGCCGTCCAGGCTATCAGTGCACGAGGGAAAATTCGTATTCGAACCTGGCAGTata                                             >  1:2236183/1‑66 (MQ=255)
                  tatCCTGATCAATGCCGTCCAGGCTATCAGTGCACGAGGGAAAATTCGTATTCGAACCTGGCAGTata                                             >  1:595410/1‑66 (MQ=255)
                  tatCCTGATCAATGCCGTCCAGGCTATCAGTGCACGAGGGAAAATTCGTATTCGAACCTGGCAGTata                                             >  1:2796989/1‑66 (MQ=255)
                  tatCCTGATCAATGCCGTCCAGGCTATCAGTGCACGAGGGAAAATTCGTATTCGAACCTGGCAGTata                                             >  1:2083292/1‑66 (MQ=255)
                  tatCCTGATCAATGCCGTCCAGGCTATCAGTGCACGAGGGAAAATTCGTATTCGAACCTGGCAGTata                                             >  1:1110414/1‑66 (MQ=255)
                  tatCCTGATCAATGCCGTCCAGGCTATCAGTGCACGAGGGAAAATTCGTATTCGAACCTGGCAGTata                                             >  1:195473/1‑66 (MQ=255)
                  tatCCTGATCAATGCCGTCCAGGCTATCAGTGCACGAGGGAAAATTCGTATTCGAACCTGGCAGTata                                             >  1:1936701/1‑66 (MQ=255)
                  tatCCTGATCAATGCCGTCCAGGCTATCAGTGCACGAGGGAAAATTCGTATTCGAACCTGGCAGTata                                             >  1:3152604/1‑66 (MQ=255)
                  tatCCTGATCAATGCCGTCCAGGCTATCAGTGCACGAGGGAAAATTCGTATTCGAACCTGGCAGTat                                              <  1:2844645/67‑2 (MQ=255)
                  tatCCTGATCAATGCCGTCCAGGCTATCAGTGCACGAGGGAAAATTCGTATTCGAACCTGGCAGTat                                              <  1:1071942/67‑2 (MQ=255)
                  tatCCTGATCAATGCCGTCCAGGCTATCAGTGCACGAGGGAAAATTCGTATTCGAACCTGGCAGTat                                              <  1:2941292/67‑2 (MQ=255)
                  tatCCTGATCAATGCCGTCCAGGCTATCAGTGCACGAGGGAAAATTCGTATTCGAACCTGGCAGTat                                              <  1:3285189/67‑2 (MQ=255)
                  tatCCTGATCAATGCCGTCCAGGCTATCAGTGCACGAGGGAAAATTCGTATTCGAACCTGGCAGTat                                              <  1:1194441/67‑2 (MQ=255)
                  tatCCTGATCAATGCCGTCCAGGCTATCAGTGCACGAGGGAAAATTCGTATTCGAACCTGGCAGTat                                              <  1:2978200/67‑2 (MQ=255)
                  tatCCTGATCAATGCCGTCCAGGCTATCAGTGCACGAGGGAAAATTCGTATTCGAACCTGGCAGTat                                              <  1:102179/67‑2 (MQ=255)
                  tatCCTGATCAATGCCGTCCAGGCTATCAGTGCACGAGGGAAAATTCGTATTCGAACCTGGCAGTat                                              <  1:1230833/67‑2 (MQ=255)
                  tatCCTGATCAATGCCGTCCAGGCTATCAGTGCACGAGGGAAAATTCGTATTCGAACCTGGCAGTat                                              <  1:3265581/67‑2 (MQ=255)
                  tatCCTGATCAATGCCGTCCAGGCTATCAGTGCACGAGGGAAAATTCGTATTCGAACCTGGCAGTat                                              <  1:1012612/67‑2 (MQ=255)
                  tatCCTGATCAATGCCGTCCAGGCTATCAGTGCACGAGGGAAAATTCGTATTCGAACCTGGCAGTat                                              <  1:3250917/67‑2 (MQ=255)
                  tatCCTGATCAATGCCGTCCAGGCTATCAGTGCACGAGGGAAAATTCGTATTCGAACCTGGCAGTat                                              <  1:1448491/67‑2 (MQ=255)
                  tatCCTGATCAATGCCGTCCAGGCTATCAGTGCACGAGGGAAAATTCGTATTCGAACCTGGCAGTat                                              <  1:1957157/67‑2 (MQ=255)
                  tatCCTGATCAATGCCGTCCAGGCTATCAGTGCACGAGGGAAAATTCGTATTCGAACCTGGCAGTat                                              <  1:1963489/67‑2 (MQ=255)
                  tatCCTGATCAATGCCGTCCAGGCTATCAGTGCACGAGGGAAAATTCGTATTCGAACCTGGCAGTat                                              <  1:1983580/67‑2 (MQ=255)
                  tatCCTGATCAATGCCGTCCAGGCTATCAGTGCACGAGGGAAAATTCGTATTCGAACCTGGCAGTat                                              <  1:2372602/67‑2 (MQ=255)
                  tatCCTGATCAATGCCGTCCAGGCTATCAGTGCACGAGGGAAAATTCGTATTCGAACCTGGCAGTat                                              <  1:235193/67‑2 (MQ=255)
                  tatCCTGATCAATGCCGTCCAGGCTATCAGTGCACGAGGGAAAATTCGTATTCGAACCTGGCAGTat                                              <  1:2342260/67‑2 (MQ=255)
                  tatCCTGATCAATGCCGTCCAGGCTATCAGTGCACGAGGGAAAATTCGTATTCGAACCTGGCAGTat                                              <  1:2071841/67‑2 (MQ=255)
                  tatCCTGATCAATGCCGTCCAGGCTATCAGTGCACGAGGGAAAATTCGTATTCGAACCTGGCAGTat                                              <  1:848917/67‑2 (MQ=255)
                  tatCCTGATCAATGCCGTCCAGGCTATCAGTGCACGAGGGAAAATTCGTATTCGAACCTGGCAGTat                                              <  1:611154/67‑2 (MQ=255)
                  tatCCTGATCAATGCCGTCCAGGCTATCAGTGCACGAGGGAAAATTCGTATTCGAACCTGCCAGTat                                              <  1:3077125/67‑2 (MQ=255)
                   atCCTGATCAATGCCGTCCAGGCTATCAGCGCACGAGGGAAAATTCGCATTCAAACCTGGCAATACAGc                                           <  1:1281718/69‑1 (MQ=255)
                    tCCTGATCAATGCCGTCCAGGCTATCAGCGCACGAGGGAAAATTCGCATTCAAACCTGGCAATACAGc                                           >  1:55848/1‑68 (MQ=255)
                    tCCTGATCAATGCCGTCCAGGCTATCAGCGCACGAGGGAAAATTCGCATTCAAACCTGGCAATACAGCg                                          >  1:658808/1‑69 (MQ=255)
                    tCCTGATCAATGCCGTCCAGGCTATCAGCGCACGAGGGAAAATTCGCATTCAAACCTGGCAATACAGCg                                          >  1:84124/1‑69 (MQ=255)
                    tCCTGATCAATGCCGTCCAGGCTATCAGCGCACGAGGGAAAATTCGCATTCAAACCTGGCAATACAGCg                                          >  1:1240268/1‑69 (MQ=255)
                    tCCTGATCAATGCCGTCCAGGCTATCAGCGCACGAGGGAAAATTCGCATTCAAACCTGGCAATACAGCg                                          >  1:137343/1‑69 (MQ=255)
                      cTGATCAATGCCGTCCAGGCTATCAGCGCACGAGGGAAAATTCGCATTCAAACCTGGCAATACAGCGAc                                        <  1:2676514/69‑1 (MQ=255)
                      cTGATCAATGCCGTCCAGGCTATCAGCGCACGAGGGAAAATTCGCATTCAAACCTGGCAATACAGCGAc                                        <  1:184743/69‑1 (MQ=255)
                        gATCAATGCCGTCCAGGCTATCAGCGCACGAGGGAAAATTCGCATTCAAACCTGGCAATACAGCGACt                                       <  1:2590189/68‑1 (MQ=255)
                           cAATGCCGTCCAGGCTATCAGCGCACGAGGGAAAATTCGCATTCAAACCTGGCAATACAGCGACTcac                                    >  1:943828/1‑68 (MQ=255)
                              tGCCGTCCAGGCTATCAGTGCACGAGTGAAAATTCGTATTCGAACCTGGCAGTat                                              <  1:1761224/55‑2 (MQ=25)
                                ccGTCCAGGCTATCAGCGCACGAGGGAAAATTCGCATTCAAACCTGGCaa                                                 <  1:2755652/50‑1 (MQ=255)
                                ccGTCCAGGCTATCAGCGCACGAGGGAAAATTCGCATTCAAACCTGGCaa                                                 <  1:1134164/50‑1 (MQ=255)
                                ccGTCCAGGCTATCAGCGCACGAGGGAAAATTCGCATTCAAACCTGGCaa                                                 <  1:2886098/50‑1 (MQ=255)
                                ccGTCCAGGCTATCAGCGCACGAGGGAAAATTCGCATTCAAACCTGGCaa                                                 <  1:2488233/50‑1 (MQ=255)
                                ccGTCCAGGCTATCAGCGCACGAGGGAAAATTCGCATTCAAACCTGGCaa                                                 <  1:3103783/50‑1 (MQ=255)
                                ccGTCCAGGCTATCAGCGCACGAGGGAAAATTCGCATTCAAACCTGGCaa                                                 <  1:738711/50‑1 (MQ=255)
                                ccGTCCAGGCTATCAGCGCACGAGGGAAAATTCGCATTCAAACCTGGCaa                                                 <  1:2403168/50‑1 (MQ=255)
                                ccGTCCAGGCTATCAGCGCACGAGGGAAAATTCGCATTCAAACCTGGCAATACAGCGACTCACAACAg                               >  1:2420907/1‑68 (MQ=255)
                                ccGTCCAGGCTATCAGCGCACGAGGGAAAATTCGCATTCAAACCTGGCAATACAGCGACTCACAACAg                               >  1:2135882/1‑68 (MQ=255)
                                         cTATCAGTGCACTAGGGAAAATTCGTATTCGAACCTGGCAGTATAGCGACTCGCAACAGGCCATTTCGa                     <  1:263251/69‑1 (MQ=25)
                                         cTATCAGTGCACGAGGGAAAATTCGTATTCGAACCTGGCAGTATAGCGACTCGCAACAGGCCATTTCGa                     <  1:2784673/69‑1 (MQ=37)
                                         cTATCAGTGCACGAGGGAAAATTCGTATTCGAACCTGGCAGTATAGCGACTCGCAACAGGCCATTTCGa                     <  1:2760562/69‑1 (MQ=37)
                                         cTATCAGTGCACGAGGGAAAATTCGTATTCGAACCTGGCAGTATAGCGACTCGCAACAGGCCATTTCGa                     <  1:2646710/69‑1 (MQ=37)
                                         cTATCAGTGCACGAGGGAAAATTCGTATTCGAACCTGGCAGTATAGCGACTCGCAACAGGCCATTTCGa                     <  1:2618777/69‑1 (MQ=37)
                                         cTATCAGTGCACGAGGGAAAATTCGTATTCGAACCTGGCAGTATAGCGACTCGCAACAGGCCATTTCGa                     <  1:2198203/69‑1 (MQ=37)
                                         cTATCAGTGCACGAGGGAAAATTCGTATTCGAACCTGGCAGTATAGCGACTCGCAACAGGCCATTTCGa                     <  1:2143471/69‑1 (MQ=37)
                                         cTATCAGTGCACGAGGGAAAATTCGTATTCGAACCTGGCAGTATAGCGACTCGCAACAGGCCATTTCGa                     <  1:1753491/69‑1 (MQ=37)
                                         cTATCAGTGCACGAGGGAAAATTCGTATTCGAACCTGGCAGTATAGCGACTCGCAACAGGCCATTTCGa                     <  1:882404/69‑1 (MQ=37)
                                                  cACGAGGGAAAATTCGCATTCAAACCTGGCAATACAGCGACTCACAACAGGCCATTTCGATAGAGGACa            <  1:1252678/69‑1 (MQ=255)
                                                  cACGAGGGAAAATTCGCATTCAAACCTGGCAATACAGCGACTCACAACAGGCCATTTCGATAGAGGACa            <  1:650196/69‑1 (MQ=255)
                                                            aaTTCGCATTCAAACCTGGCAATACAGCGACTCACAACAGGCCATTTCGATAGAGGACAACGGCTGTgg  <  1:958232/69‑1 (MQ=255)
                                                            aaTTCGCATTCAAACCTGGCAATACAGCGACTCACAACAGGCCATTTCGATAGAGGACAACGGCTGTgg  <  1:321168/69‑1 (MQ=255)
                                                            aaTTCGCATTCAAACCTGGCAATACAGCGACTCACAACAGGCCATTTCGATAGAGGACAACGGCTGTgg  <  1:1425343/69‑1 (MQ=255)
                                                            aaTTCGCATTCAAACCTGGCAATACAGCGACTCACAACAGGCCATTTCGATAGAGGACAACGGCTGTgg  <  1:1553342/69‑1 (MQ=255)
                                                            aaTTCGCATTCAAACCTGGCAATACAGCGACTCACAACAGGCCATTTCGATAGAGGACAACGGCTGTgg  <  1:819696/69‑1 (MQ=255)
                                                            aaTTCGCATTCAAACCTGGCAATACAGCGACTCACAACAGGCCATTTCGATAGAGGACAACGGCTGTgg  <  1:1744162/69‑1 (MQ=255)
                                                            aaTTCGCATTCAAACCTGGCAATACAGCGACTCACAACAGGCCATTTCGATAGAGGACAACGGCTGTgg  <  1:759604/69‑1 (MQ=255)
                                                            aaTTCGCATTCAAACCTGGCAATACAGCGACTCACAACAGGCCATTTCGATAGAGGACAACGGCTGTgg  <  1:226194/69‑1 (MQ=255)
                                                            aaTTCGCATTCAAACCTGGCAATACAGCGACTCACAACAGGCCATTTCGATAGAGGACAACGGCTGTgg  <  1:3197472/69‑1 (MQ=255)
                                                                  |                                                              
CAAACAGGTACTACTGAATATCCTGATCAATGCCGTCCAGGCTATCAGCGCACGAGGGAAAATTCGCATTCAAACCTGGCAATACAGCGACTCACAACAGGCCATTTCGATAGAGGACAACGGCTGTGG  >  minE/1449094‑1449222

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: