Predicted mutation
evidence seq id position mutation freq annotation gene description
RA minE 1,449,142 C→T 54.9% S516S (AGC→AGT 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,1420CT54.9% 15.2 / 85.0 102S516S (AGC→AGTatoSsensory histidine kinase in two‑component regulatory system with AtoC
Reads supporting (aligned to +/- strand):  ref base C (23/23);  new base T (19/37);  total (42/60)
Fisher's exact test for biased strand distribution p-value = 1.11e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.84e-01

Only 100 of 102 total aligned reads displayed.

GCTCAAACAGGTACTACTGAATATCCTGATCAATGCCGTCCAGGCTATCAGCGCACGAGGGAAAATTCGCATTCAAACCTGGCAATACAGCGACTCACAACAGGCCATTTCGA  >  minE/1449091‑1449203
                                                   |                                                             
gCTCAAACAGGTACTACTGAATATCCTGATCAATGCCGTCCAGGCTATCAGCGCACGAGGGAAAATTc                                               >  1:388295/1‑68 (MQ=255)
   cAAACAGGTACTACTGAATATCCTGATCAATGCCGTCCAGGCTATCAGTGCACGAGGGAAAAt                                                 >  1:2117517/1‑63 (MQ=255)
   cAAACAGGTACTACTGAATATCCTGATCAATGCCGTCCAGGCTATCAGTGCACGAGGGAAAATTCGTAt                                           >  1:2534678/1‑69 (MQ=255)
   cAAACAGGTACTACTGAATATCCTGATCAATGCCGTCCAGGCTATCAGTGCACGAGGGAAAATTCGTAt                                           >  1:1772107/1‑69 (MQ=255)
   cAAACAGGTACTACTGAATATCCTGATCAATGCCGTCCAGGCTATCAGTGCACGAGGGAAAATTCGTAt                                           >  1:1533919/1‑69 (MQ=255)
   cAAACAGGTACTACTGAATATCCTGATCAATGCCGTCCAGGCTATCAGTGCACGAGGGAAAATTCGTAt                                           >  1:1168308/1‑69 (MQ=255)
   cAAACAGGTACTACTGAATATCCTGATCAATGCCGTCCAGGCTATCAGTGCACGAGGGAAAATTCGTAt                                           >  1:3013087/1‑69 (MQ=255)
             ctactGAATATCCTGATCAATGCCGTCCAGGCTATCAGCGCACGAGGGAAAATTCGCATTCAAACCTg                                  >  1:1153204/1‑68 (MQ=255)
               actTAATATCCTGATCAATGCCGTCCAGGCTATCAGCGCACGAGGGAAAATTCGCATTCa                                        >  1:2904150/1‑60 (MQ=255)
               actGAATATCCTGATCAATGCCGTCCAGGCTATCAGCGCACGAGGGAAAATTCGCATTCa                                        >  1:2390532/1‑60 (MQ=255)
               actGAATATCCTGATCAATGCCGTCCAGGCTATCAGCGCACGAGGGAAAATTCGCATTCa                                        <  1:1712866/60‑1 (MQ=255)
               actGAATATCCTGATCAATGCCGTCCAGGCTATCAGCGCACGAGGGAAAATTCGCATTCa                                        <  1:1946323/60‑1 (MQ=255)
               actGAATATCCTGATCAATGCCGTCCAGGCTATCAGCGCACGAGGGAAAATTCGCATTCa                                        <  1:1610935/60‑1 (MQ=255)
               actGAATATCCTGATCAATGCCGTCCAGGCTATCAGCGCACGAGGGAAAATTCGCATTCa                                        <  1:1582214/60‑1 (MQ=255)
               actGAATATCCTGATCAATGCCGTCCAGGCTATCAGCGCACGAGGGAAAATTCGCATTCa                                        >  1:299007/1‑60 (MQ=255)
               actGAATATCCTGATCAATGCCGTCCAGGCTATCAGCGCACGAGGGAAAATTCGCATTCa                                        >  1:2173607/1‑60 (MQ=255)
               actGAATATCCTGATCAATGCCGTCCAGGCTATCAGCGCACGAGGGAAAATTCGCATTCa                                        >  1:1469585/1‑60 (MQ=255)
               actGAATATCCTGATCAATGCCGTCCAGGCTATCAGCGCACGAGGGAAAATTCGCATTCa                                        >  1:2408954/1‑60 (MQ=255)
               actGAATATCCTGATCAATGCCGTCCAGGCTATCAGCGCACGAGGGAAAATTCGCATTCa                                        <  1:510765/60‑1 (MQ=255)
               actGAATATCCTGATCAATGCCGTCCAGGCTATCAGCGCACGAGGGAAAATTCGCATTCa                                        >  1:2551468/1‑60 (MQ=255)
               actGAATATCCTGATCAATGCCGTCCAGGCTATCAGCGCACGAGGGAAAATTCGCATTCa                                        >  1:257513/1‑60 (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:761387/67‑1 (MQ=255)
                  gAATATCCTGATCAATGCCGTCCAGGCTATCAGCGCACGAGGGAAAATTCGCATTCAAACCTGGCaa                              <  1:1663734/67‑1 (MQ=255)
                  gAATATCCTGATCAATGCCGTCCAGGCTATCAGCGCACGAGGGAAAATTCGCATTCAAACCTGGCaa                              <  1:1495167/67‑1 (MQ=255)
                     tatCCTGATCAATGCCGTCCAGGCTATCAGTGCACGAGGGAAAATTCGTATTCGAACCTGGCAGTata                          >  1:2469579/1‑66 (MQ=255)
                     tatCCTGATCAATGCCGTCCAGGCTATCAGTGCACGAGGGAAAATTCGTATTCGAACCTGGCAGTata                          >  1:1110414/1‑66 (MQ=255)
                     tatCCTGATCAATGCCGTCCAGGCTATCAGTGCACGAGGGAAAATTCGTATTCGAACCTGGCAGTata                          >  1:3152604/1‑66 (MQ=255)
                     tatCCTGATCAATGCCGTCCAGGCTATCAGTGCACGAGGGAAAATTCGTATTCGAACCTGGCAGTata                          >  1:2796989/1‑66 (MQ=255)
                     tatCCTGATCAATGCCGTCCAGGCTATCAGTGCACGAGGGAAAATTCGTATTCGAACCTGGCAGTata                          >  1:510310/1‑66 (MQ=255)
                     tatCCTGATCAATGCCGTCCAGGCTATCAGTGCACGAGGGAAAATTCGTATTCGAACCTGGCAGTata                          >  1:2590608/1‑66 (MQ=255)
                     tatCCTGATCAATGCCGTCCAGGCTATCAGTGCACGAGGGAAAATTCGTATTCGAACCTGGCAGTata                          >  1:12699/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:2083292/1‑66 (MQ=255)
                     tatCCTGATCAATGCCGTCCAGGCTATCAGTGCACGAGGGAAAATTCGTATTCGAACCTGGCAGTata                          >  1:1936701/1‑66 (MQ=255)
                     tatCCTGATCAATGCCGTCCAGGCTATCAGTGCACGAGGGAAAATTCGTATTCGAACCTGGCAGTata                          >  1:195473/1‑66 (MQ=255)
                     tatCCTGATCAATGCCGTCCAGGCTATCAGTGCACGAGGGAAAATTCGTATTCGAACCTGGCAGTat                           <  1:2844645/67‑2 (MQ=255)
                     tatCCTGATCAATGCCGTCCAGGCTATCAGTGCACGAGGGAAAATTCGTATTCGAACCTGGCAGTat                           <  1:1963489/67‑2 (MQ=255)
                     tatCCTGATCAATGCCGTCCAGGCTATCAGTGCACGAGGGAAAATTCGTATTCGAACCTGGCAGTat                           <  1:1957157/67‑2 (MQ=255)
                     tatCCTGATCAATGCCGTCCAGGCTATCAGTGCACGAGGGAAAATTCGTATTCGAACCTGGCAGTat                           <  1:1230833/67‑2 (MQ=255)
                     tatCCTGATCAATGCCGTCCAGGCTATCAGTGCACGAGGGAAAATTCGTATTCGAACCTGGCAGTat                           <  1:1194441/67‑2 (MQ=255)
                     tatCCTGATCAATGCCGTCCAGGCTATCAGTGCACGAGGGAAAATTCGTATTCGAACCTGGCAGTat                           <  1:3285189/67‑2 (MQ=255)
                     tatCCTGATCAATGCCGTCCAGGCTATCAGTGCACGAGGGAAAATTCGTATTCGAACCTGGCAGTat                           <  1:3265581/67‑2 (MQ=255)
                     tatCCTGATCAATGCCGTCCAGGCTATCAGTGCACGAGGGAAAATTCGTATTCGAACCTGGCAGTat                           <  1:3250917/67‑2 (MQ=255)
                     tatCCTGATCAATGCCGTCCAGGCTATCAGTGCACGAGGGAAAATTCGTATTCGAACCTGGCAGTat                           <  1:611154/67‑2 (MQ=255)
                     tatCCTGATCAATGCCGTCCAGGCTATCAGTGCACGAGGGAAAATTCGTATTCGAACCTGGCAGTat                           <  1:1012612/67‑2 (MQ=255)
                     tatCCTGATCAATGCCGTCCAGGCTATCAGTGCACGAGGGAAAATTCGTATTCGAACCTGGCAGTat                           <  1:1071942/67‑2 (MQ=255)
                     tatCCTGATCAATGCCGTCCAGGCTATCAGTGCACGAGGGAAAATTCGTATTCGAACCTGGCAGTat                           <  1:2941292/67‑2 (MQ=255)
                     tatCCTGATCAATGCCGTCCAGGCTATCAGTGCACGAGGGAAAATTCGTATTCGAACCTGGCAGTat                           <  1:2978200/67‑2 (MQ=255)
                     tatCCTGATCAATGCCGTCCAGGCTATCAGTGCACGAGGGAAAATTCGTATTCGAACCTGGCAGTat                           <  1:102179/67‑2 (MQ=255)
                     tatCCTGATCAATGCCGTCCAGGCTATCAGTGCACGAGGGAAAATTCGTATTCGAACCTGGCAGTat                           <  1:848917/67‑2 (MQ=255)
                     tatCCTGATCAATGCCGTCCAGGCTATCAGTGCACGAGGGAAAATTCGTATTCGAACCTGGCAGTat                           <  1:2120408/67‑2 (MQ=255)
                     tatCCTGATCAATGCCGTCCAGGCTATCAGTGCACGAGGGAAAATTCGTATTCGAACCTGGCAGTat                           <  1:2342260/67‑2 (MQ=255)
                     tatCCTGATCAATGCCGTCCAGGCTATCAGTGCACGAGGGAAAATTCGTATTCGAACCTGGCAGTat                           <  1:235193/67‑2 (MQ=255)
                     tatCCTGATCAATGCCGTCCAGGCTATCAGTGCACGAGGGAAAATTCGTATTCGAACCTGGCAGTat                           <  1:2372602/67‑2 (MQ=255)
                     tatCCTGATCAATGCCGTCCAGGCTATCAGTGCACGAGGGAAAATTCGTATTCGAACCTGGCAGTat                           <  1:1448491/67‑2 (MQ=255)
                     tatCCTGATCAATGCCGTCCAGGCTATCAGTGCACGAGGGAAAATTCGTATTCGAACCTGGCAGTat                           <  1:2071841/67‑2 (MQ=255)
                     tatCCTGATCAATGCCGTCCAGGCTATCAGTGCACGAGGGAAAATTCGTATTCGAACCTGGCAGTat                           <  1:687068/67‑2 (MQ=255)
                     tatCCTGATCAATGCCGTCCAGGCTATCAGTGCACGAGGGAAAATTCGTATTCGAACCTGGCAGTat                           <  1:1983580/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:137343/1‑69 (MQ=255)
                       tCCTGATCAATGCCGTCCAGGCTATCAGCGCACGAGGGAAAATTCGCATTCAAACCTGGCAATACAGCg                       >  1:2378509/1‑69 (MQ=255)
                       tCCTGATCAATGCCGTCCAGGCTATCAGCGCACGAGGGAAAATTCGCATTCAAACCTGGCAATACAGCg                       >  1:1240268/1‑69 (MQ=255)
                       tCCTGATCAATGCCGTCCAGGCTATCAGCGCACGAGGGAAAATTCGCATTCAAACCTGGCAATACAGCg                       >  1:84124/1‑69 (MQ=255)
                       tCCTGATCAATGCCGTCCAGGCTATCAGCGCACGAGGGAAAATTCGCATTCAAACCTGGCAATACAGCg                       >  1:658808/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)
                           gATCAATGCCGTCCAGGCTATCAGCGCACGAGGGAAAATTCGCATTCAAACCTGGCAATACAGCGACt                    <  1:3293919/68‑1 (MQ=255)
                              cAATGCCGTCCAGGCTATCAGCGCACGAGGGAAAATTCGCATTCAAACCTGGCAATACAGCGACTcac                 >  1:943828/1‑68 (MQ=255)
                                 tGCCGTCCAGGCTATCAGTGCACGAGTGAAAATTCGTATTCGAACCTGGCAGTat                           <  1:1761224/55‑2 (MQ=25)
                                   ccGTCCAGGCTATCAGCGCACGAGGGAAAATTCGCATTCAAACCTGGCaa                              <  1:3103783/50‑1 (MQ=255)
                                   ccGTCCAGGCTATCAGCGCACGAGGGAAAATTCGCATTCAAACCTGGCaa                              <  1:1134164/50‑1 (MQ=255)
                                   ccGTCCAGGCTATCAGCGCACGAGGGAAAATTCGCATTCAAACCTGGCaa                              <  1:1496253/50‑1 (MQ=255)
                                   ccGTCCAGGCTATCAGCGCACGAGGGAAAATTCGCATTCAAACCTGGCaa                              <  1:738711/50‑1 (MQ=255)
                                   ccGTCCAGGCTATCAGCGCACGAGGGAAAATTCGCATTCAAACCTGGCaa                              <  1:2403168/50‑1 (MQ=255)
                                   ccGTCCAGGCTATCAGCGCACGAGGGAAAATTCGCATTCAAACCTGGCaa                              <  1:2488233/50‑1 (MQ=255)
                                   ccGTCCAGGCTATCAGCGCACGAGGGAAAATTCGCATTCAAACCTGGCaa                              <  1:2755652/50‑1 (MQ=255)
                                   ccGTCCAGGCTATCAGCGCACGAGGGAAAATTCGCATTCAAACCTGGCaa                              <  1:2886098/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:2646710/69‑1 (MQ=37)
                                            cTATCAGTGCACGAGGGAAAATTCGTATTCGAACCTGGCAGTATAGCGACTCGCAACAGGCCATTTCGa  <  1:2618777/69‑1 (MQ=37)
                                            cTATCAGTGCACGAGGGAAAATTCGTATTCGAACCTGGCAGTATAGCGACTCGCAACAGGCCATTTCGa  <  1:2760562/69‑1 (MQ=37)
                                            cTATCAGTGCACGAGGGAAAATTCGTATTCGAACCTGGCAGTATAGCGACTCGCAACAGGCCATTTCGa  <  1:2198203/69‑1 (MQ=37)
                                            cTATCAGTGCACGAGGGAAAATTCGTATTCGAACCTGGCAGTATAGCGACTCGCAACAGGCCATTTCGa  <  1:2143471/69‑1 (MQ=37)
                                            cTATCAGTGCACGAGGGAAAATTCGTATTCGAACCTGGCAGTATAGCGACTCGCAACAGGCCATTTCGa  <  1:2784673/69‑1 (MQ=37)
                                            cTATCAGTGCACGAGGGAAAATTCGTATTCGAACCTGGCAGTATAGCGACTCGCAACAGGCCATTTCGa  <  1:2793104/69‑1 (MQ=37)
                                            cTATCAGTGCACGAGGGAAAATTCGTATTCGAACCTGGCAGTATAGCGACTCGCAACAGGCCATTTCGa  <  1:1753491/69‑1 (MQ=37)
                                            cTATCAGTGCACGAGGGAAAATTCGTATTCGAACCTGGCAGTATAGCGACTCGCAACAGGCCATTTCGa  <  1:1232717/69‑1 (MQ=37)
                                            cTATCAGTGCACGAGGGAAAATTCGTATTCGAACCTGGCAGTATAGCGACTCGCAACAGGCCATTTCGa  <  1:882404/69‑1 (MQ=37)
                                            cTATCAGCGCACGAGGGAAAATTCGCATTCAAACCTGGCAATACAGCGACTCACAACAGGCCATTTc    >  1:3004473/1‑67 (MQ=255)
                                                   |                                                             
GCTCAAACAGGTACTACTGAATATCCTGATCAATGCCGTCCAGGCTATCAGCGCACGAGGGAAAATTCGCATTCAAACCTGGCAATACAGCGACTCACAACAGGCCATTTCGA  >  minE/1449091‑1449203

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: