Predicted mutation
evidence seq id position mutation annotation gene description
RA NC_000913 3,849,980 G→A L53L (CTG→TTG)  uhpB ← sensory histidine kinase UhpB

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009133,849,9800GA96.0% 77.2 / ‑4.5 25L53L (CTG→TTG) uhpBsensory histidine kinase UhpB
Reads supporting (aligned to +/- strand):  ref base G (1/0);  new base A (20/4);  total (21/4)
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

TACGGGCAGTAACGTCAGTAAACTACCGATCATCAATAACGGAAAATGGGTTAAACCGACCGCCTGCGTTAGCCAGTAAATCAGCAGCCACTCCGCGCCCAGCAATACGGGCCAGTATCCGCGCGGGCATTGCAGCATTAGCCCCAGACGCAGACCAAACGGAAATAACAGCACCGCCATATCAGGGCGCTCAACCAGATGCAGGCTGATACTCCACAGGCAAAACCATGCGGCAGAGAAGATAAAAAAGCAGGCAATA  >  NC_000913/3849846‑3850104
                                                                                                                                      |                                                                                                                            
tACGGGCAGTAACGTCAGTAAACTACCGATCATCAATAACGGAAAATGGGTTAAACCGACCGCCTGCGTTAGCCAGTAAATCAGCAGCCACTCCGCGCCCAGCAATACGGGCCAGTATCCGCGCGGGCATTGCAACATTa                                                                                                                         <  1:637523/140‑1 (MQ=255)
           aCGTCAGTAAACTACCGATCATCAATAACGGAAAATGGGTTAAACCGACCGCCTGCGTTAGCCAGTAAATCAGCAGCCACTCCGCGCCCAGCAATACGGGCCAGTATCCGCGCGGGCATTGCAACATTAGCCCCAGACGc                                                                                                              <  1:699536/140‑1 (MQ=255)
             gTCAGTAAACTACCGATCATCAATAACGGAAAATGGGTTAAACCGACCGCCTGCGTTAGCCAGTAAATCAGCAGCCACTCCGCGCCCAGCAATACGGGCCAGTATCCGCGCGGGCATTGCAACATTAGCCCCAGACGCAg                                                                                                            >  1:477320/1‑140 (MQ=255)
                aGTAAACTACCGATCATCAATAACGGAAAATGGGTTAAACCGACCGCCTGCGTTAGCCAGTAAATCAGCAGCCACTCCGCGCCCAGCAATACGGGCCAGTATCCGCGCGGGCATTGCAACATTAGCCCCAGACGCAGAcc                                                                                                         >  1:780351/1‑140 (MQ=255)
                 gTAAACTACCGATCATCAATAACGGAAAATGGGTTAAACCGACCGCCTGCGTTAGCCAGTAAATCAGCAGCCACTCCGCGCCCAGCAATACGGGCCAGTATCCGCGCGGGCATTGCAACATTAGCCCCAGACGCAGACCa                                                                                                        >  2:64942/1‑140 (MQ=255)
                  tataCTCCCGATCACCAATAACGGAAAATGGGTAAAACGGGCCGCCTGCGTTAGCCAGTAAATCAGCAGCCACTCCGCGCCCAGCAATACGGGCCAGTATCCGCGCGGGCATTGCAACATTAGCCCCAGACGCAGACCaa                                                                                                       <  2:938120/137‑1 (MQ=255)
                       tACCGATCATAAATAACGGAAAATGGGTTAAACCGACCGCCTGCGTTAGCCAGTAAATCAGCAGCCACTCCGCGCCCAGCAATACGGGCCAGTATCCGCGCGGGCATTGCAACATTAGCCCCAGACGCAGACCAAACGGa                                                                                                  <  1:968436/140‑1 (MQ=255)
                         ccGATCATCAATAACGGAAAATGGGTTAAACCGACCGCCTGCGTTAGCCAGTAAATCAGCAGCCACTCCGCGCCCAGCAATACGGGCCAGTATCCGCGCGGGCATTGCAACATTAGCCCCAGACGCAGACCAAACGGaaa                                                                                                >  1:49591/1‑140 (MQ=255)
                           gATCATCAATAACGGAAAATGGGTTAAACCGACCGCCTGCGTTAGCCAGTAAATCAGCAGCCACTCCGCGCCCAGCAATACGGGCCAGTATCCGCGCGGGCATTGCAACATTAGCCCCAGACGCAGACCAAACGGAAATa                                                                                              >  1:727028/1‑140 (MQ=255)
                              catcaATAACGGAAAATGGGTTAAACCGACCGCCTGCGTTAGCCAGTAAATCAGCAGCCACTCCGCGCCCAGCAATACGGGCCAGTATCCGCGCGGGCATTGCAACATTAGCCCCAGACGCAGACCAAACGGAAATAACa                                                                                           >  1:700118/1‑140 (MQ=255)
                              catcaATAACGGAAAATGGGTTAAACCGACCGCCTGCGTTAGCCAGTAAATCAGCAGCCACTCCGCGCCCAGCAATACGGGCCAGTATCCGCGCGGGCATTGCAACATTAGCCCCAGACGCAGACCAAACGGAAATAACa                                                                                           >  1:780254/1‑140 (MQ=255)
                                   aTAACGGAAAATGGGTTAAACCGACCGCCTGCGTTAGCCAGTAAATCAGCAGCCACTCCGCGCCCAGCAATACGGGCCAGTATCCGCGCGGGCATTGCAACATTAGCCCCAGACGCAGACCAAACGGAAATAACAGCAcc                                                                                      >  2:328739/1‑140 (MQ=255)
                                   aTAACGGAAAATGGGTTAAACCGACCGCCTGCGTTAGCCAGTAAATCAGCAGCCACTCCGCGCCCAGCAATACGGGCCAGTATCCGCGCGGGCATTGCAACATTAGCCCCAGACGCAGACCAAACGGAAATAACAGCAcc                                                                                      >  2:202908/1‑140 (MQ=255)
                                           aaaTGGGTTAAACCGACCGCCTGCGTTAGCCAGTAAATCAGCAGCCACTCCGCGCCCAGCAATACGGGCCAGTATCCGCGCGGGCATTGCAACATTAGCCCCAGACGCAGACCAAACGGGAATAACAGCACCGCCATATc                                                                              >  1:860210/1‑140 (MQ=255)
                                            aaTGGGTTAAACCGACCGCCTGCGTTAGCCAGTAAATCAGCAGCCACTCCGCGCCCAGCAATACGGGCCAGTATCCGCGCGGGCATTGCAACATTAGCCCCAGACGCAGACCAAACGGAAATAACAGCACCGCCATATCa                                                                             >  1:186673/1‑140 (MQ=255)
                                                 gTTAAACCGACCGCCTGCGTTAGCCAGTAAATCAGCAGCCACTCCGCGCCCAGCAATACGGGCCAGTATCCGCGCGGGCATTGCAACATTAGCCCCAGACGCAGACCAAACGGAAATAACAGCACCGCCATATCAGGgcg                                                                        >  1:685471/1‑140 (MQ=255)
                                                          accgCCTGCGTTAGCCAGTAAATCAGCAGCCACTCCGCGCCCAGCAATACGGGCCAGTATCCGCGCGGGCATTGCAACATTAGCCCCAGACGCAGACCAAACGGAAATAACAGCACCGCCATATCAGGGCGCTCAACCAg                                                               >  1:304373/1‑140 (MQ=255)
                                                               cTGCGTTAGCCAGTAAATCAGCAGCCACTCCGCGCCCAGCAATACGGGCCAGTATCCGCGCGGGCATTGCAGCATTAGCCCCAGACGCAGACCAAACGGAAATAACAGCACCGCCATATCAGGGCGCTCAACCAGATGCa                                                          >  1:1021273/1‑140 (MQ=255)
                                                                            aaaaTCAGCAGCCCCCCCGCCCCCAGCAATACGGGCCAGTAGCCGCGGGGGCATTGCAACATTAGCCCCGGACGCAGACAAAACGGAAATAACGGCGCCGCCATATCAGGGCGCTCAACCAGATGCAGGCTGATACTCca                                             <  2:633612/139‑1 (MQ=255)
                                                                                              gcgcCCAGCAATACGGGCCAGTATCCGCGCGGGCATTGCAACATTAGCCCCAGACGCAGACCAAACGGAAATAACAGCACCGCCATATCAGGGCGCTCAACCAGATGCAGGCTGATACTCCACAGGCAAAACCATGCGGc                           >  1:1014557/1‑140 (MQ=255)
                                                                                                    aGCAATACGGGCCAGTATCCGCGCGGGCATTGCAACATTAGCCCCAGACGCAGACCAAACGGAAATAACAGCACCGCCATATCAGGGCGCTCAACCAGATGCAGGCTGATACTCCACAGGCAAAACCATGCGGCAGagaa                     >  1:984341/1‑140 (MQ=255)
                                                                                                     gCAATACGGGCCAGTATCCGCGCGGGCATTGCAACATTAGCCCCAGACGCAGACCAAACGGAAATAACAGCACCGCCATATCAGGGCGCTCAACCAGATGCAGGCTGATACTCCACAGGCAAAACCATGCGGCAGagaag                    >  1:517196/1‑140 (MQ=255)
                                                                                                      cAATACGGGCCAGTATCCGCGCGGGCATTGCAACATTAGCCCCAGACGCAGACCAAACGGAAATAACAGCACCGCCATATCAGGGCGCTCAACCAGATGCAGGCTGATACTCCACAGGCAAAACCATGCGGCAGagaaga                   >  1:918799/1‑140 (MQ=255)
                                                                                                              gCCAGTATCCGCGCGGGCATTGCAACATTAGCCCCAGACGCAGACCAAACGGAAATAACAGCACCGCCATATCAGGGCGCTCAACCAGATGCAGGCTGATACTCCACAGGCAAAACCATGCGGCAGAGAAGATAAAAAAg           >  1:996598/1‑140 (MQ=255)
                                                                                                              gCCAGTATCCGCGCGGGCATTGCAACATTAGCCCCAGACGCAGACCAAACGGAAATAACAGCACCGCCATATCAGGGCGCTCAACCAGATGCAGGCTGATACTCCACAGGCAAAACCATGCGGCAGAGAAGATAAAAAAg           >  1:121158/1‑140 (MQ=255)
                                                                                                                       cgcgcgGGCATTGCAACATTAGCCCCAGACGCAGACCAAACGGAAATAACAGCACCGCCATATCAGGGCGCTCAACCAGATGCAGGCTGATACTCCACAGGCAAAACCATGCGGCAGAGAAGATAAAAAAGCAGGCAATa  >  1:25164/1‑140 (MQ=255)
                                                                                                                                      |                                                                                                                            
TACGGGCAGTAACGTCAGTAAACTACCGATCATCAATAACGGAAAATGGGTTAAACCGACCGCCTGCGTTAGCCAGTAAATCAGCAGCCACTCCGCGCCCAGCAATACGGGCCAGTATCCGCGCGGGCATTGCAGCATTAGCCCCAGACGCAGACCAAACGGAAATAACAGCACCGCCATATCAGGGCGCTCAACCAGATGCAGGCTGATACTCCACAGGCAAAACCATGCGGCAGAGAAGATAAAAAAGCAGGCAATA  >  NC_000913/3849846‑3850104

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