Predicted mutation
evidence seq id position mutation annotation gene description
RA NC_000913 4,332,310 C→T G321G (GGG→GGA basS ← sensor histidine kinase BasS

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009134,332,3100CT85.0% 47.1 / ‑2.6 20G321G (GGG→GGAbasSsensor histidine kinase BasS
Reads supporting (aligned to +/- strand):  ref base C (0/3);  new base T (15/2);  total (15/5)
Fisher's exact test for biased strand distribution p-value = 8.77e-03
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

TAATCTCTGACGCGCATACTCTCCTCCAGGTTAACGGAGGAGAGTGCAATGAAAAACCGTGTTTATGAAAGTTTAACTACCGTGTTCAGCGTGCTGGTGGTCAGCAGCTTTCTTTATATCTGGTTTGCCACGTACTGATCTTTCTTCAGCCGTACCCAGGCCCGCGTGCCGGAAGTCTCTTGCCGGTTTTGCAGGAAAAACTGCCCGTGATGCAACTGTGTAATGCGGCTGACAATACTTAACCCCAGACCAATCCCGCCATAACGGCTGTCCATACGTACAAACGCTTTACTCAACTCCCCGCATTTACTCTCATCAATACCTGGTCCTTCATCTTCAACTGCCATGACCGCTCCGTCATCTTCTTGCAGCTTAATCATAATGTTGCTGCCTTGCGGGCTGTAACGATGGGCGTTTTCTACCAGGTTTCGCAATAACATCCGCAGCAGGGTTGCATCACCCTGAACGGTGATGTC  >  NC_000913/4332068‑4332543
                                                                                                                                                                                                                                                  |                                                                                                                                                                                                                                         
tAATCTCTGACGCGCATACTCTCCTCCAGGTTAACGGAGGAGAGTGCAATGAAAAACCGTGTTTATGAAAGTTTAACTACCGTGTTCAGCGTGCTGGTGGTCAGCAGCTTTCTTTATATCTGGTTTGCCACGTACTGATCTTTCTTCAGCCGTACCCAGGCCCGCGTGCCGGAAGTCTCTTGCCGGTTTTGCAGGAAAAACTGCCCGTGATGCAACTGTGTAATGCGGCTGACAATACTTAATCCCAGAc                                                                                                                                                                                                                                    <  1:130111/250‑1 (MQ=255)
        gACGAGCAGACTCTCCTCCAGGTTAACGGAGGAGAGTGCAATGAAAAACCGTGTTTATGAAAGTTTAACGACCTTGTTCAGCGTGCTGGTGGTCAGCAGCTTTCTTTATATCTGGTTTGCCACGTACTGATCTTTCTTCAGCCGTACCCAGGCCCGCGGGCCGGAAGTCTCGTGCCGTTTTTGCAGGAAAAACTGCCCGTGATGCAACTGTGTAATTCGGCTGACAATAGGTAATCCCATACCAAGCCCGc                                                                                                                                                                                                                           <  2:163882/251‑1 (MQ=255)
             gcATACTCTCCTCCAGGTTAACGGAGGAGAGTGCAATGAAAAACCGTGTTTATGAAAGTTTAACTACCGTGTTCAGCGTGCTGGTGGTCAGCAGCTTTCTTTTTATCTGGTTTGCCACGTACTGATCTTTCTTCAGCCGTACCCAGGCCCGCGTGCCGGAAGTCTCTTGCCGGTTTTGCAGGAAAAACTGCCCGTGATGCAACTGTGTAATGCGGCTGACAATACTTAATCCCAGACCAATCCCGCCATa                                                                                                                                                                                                                       >  1:32222/1‑250 (MQ=255)
                 actatcctccAGGTTAACGGAGGAGAGTGCAATGAAAAACCGTGTTTATGAAAGTTTAACTACCGTGTTCAGCGTGCTGGTGGTCAGCAGCTTTCTTTATATCTGGTTTGCCACGTACTGATCTTTCTTCAGCCGTACCCAGGCCCGCGTGCCGGAAGTCTCTTGCCGGTTTTGCAGGAAAAACTGCCCGTGATGCAACTGTGTAATGCGGCTGACAATACTTAATCCCAGACCAATCCCGCCATAACGGc                                                                                                                                                                                                                  >  1:370405/5‑251 (MQ=255)
                             gTTAACGGAGGAGAGTGCAATGAAAAACCGTGTTTATGAAAGTTTAACTACCGTGTTCAGCGTGCTGGTGGTCAGCAGCTTTCTTTATATCTGGTTTGCCACGTACTGATCTTTCTTCAGCCGTACCCAGGCCCGCGTGCCGGAAGTCTCTTGCCGGTTTTGCAGGAAAAACTGCCCGTGATGCAACTGTGTAATGCGGCTGACAATACTTAATCCCAGACCAATCCCGCCATAACGGCTGTCCATACGt                                                                                                                                                                                                       >  1:80868/1‑250 (MQ=255)
                                   ggaggaGAGTGCAATGAAAAACCGTGTTTATGAAAGTTTAACTACCGTGTTCAGCGTGCTGGTGGTCAGCAGCTTTCTTTATATCTGGTTTGCCACGTAATGATATTTCTTCAACAGTAACAAGGCACAAGAGAAGAAAGTCTATTGCAGGTTTTGAAAAAAAAACTACACGTGATGCACCTGTGTAATGCGGCGGACAATACTTAATCCCAGCCCACTCCCGCCATAACGGCTGTCCATACGTACAAACg                                                                                                                                                                                                >  1:283795/1‑251 (MQ=255)
                                          agTGCAATGAAAAACCGTGTTTATGAAAGTTTAACTACCGTGTTCAGCGTGCTGGTGGTCAGCAGCTTTCTTTATATCTGGTTTGCCACGTACTGATCTTTCTTCAGCCGTACCCAGGCCCGCGTGCCGGAAGTCTCTTGCCGGTTTTGCAGGAAAAACTGCCCGTGATGCAACTGTGTAATGCGGCTGACAATACTTAATCCCAGACCAATCCCGCCATAACGGCTGTCCATACGTACAAACGCTTTACt                                                                                                                                                                                         >  1:64447/1‑251 (MQ=255)
                                                                                cgttttCAGGGTGTTGGTGGTCAGAAGTTTTCTTTATATTTGGTTTGCCACGTCCTGATTTTTCTTCAGCCGTACCCAGGCCCGCGTGCCGGAATTCTCTTGCGGTTTTTGCAGAAAAACCTGCCGGTGATCCATCTGGGAAAGGCGGCTGAAAATCTTTAACCCCAACCCAATCCCGCCATAACGGCTGTCCATACGTACAAACGCTTTACTCAACTCCCCGCATTTACTCTCATCAATACCTGGTCCtt                                                                                                                                                   <  1:134303/247‑1 (MQ=255)
                                                                                                                   atatCTGGTTTGCCACGTACTGATCTTTCTTCAGCCGTACCCAGGCCCGCGTGCCGGAAGTCTCTTGCCGGTTTTGCAGGAAAAACTGCCCGTGATGCAACTGTGTAATGCGGCTGACAATACTTAATCCCAGACCAATCCCGCCATAACGGCTGTCCATACGTACAAACGCTTTACTCAACTCCCCGCATTTACTCTCATCAATACCTGGTCCTTCATCTTCAACTGCCATGACCGCTCCGTCAtcttct                                                                                                                >  2:136583/1‑251 (MQ=255)
                                                                                                                                     aCTGATCTTTCTTCAGCCGTACCCAGGCCCGCGTGCCGGAAGTCTCTTGCCGGTTTTGCAGGAAAAACTGCCCGTGATGCAACTGTGTAATGCGGCTGACAATACTTAATCCCAGACCAATCCCGCCATAACGGCTGTCCATACGTACAAACGCTTTACTCAACTCCCCGCATTTACTCTCATCAATACCTGGTCCTTCATCTTCAACTGCCATGACCGCTCCGTCATCTTCTTGCAGCTTAATCATAATg                                                                                              >  1:274578/1‑251 (MQ=255)
                                                                                                                                     aCTGATCTTTCTTCAGCCGTACCCAGGCCCGCGTGCCGGAAGTCTCTTGCCGGTTTTGCAGGAAAAACTGCCCGTGATGCAACTGTGTAATGCGGCTGACAATACTTAATCCCAGACCAATCCCGCCATAACGGCTGTCCATACGTACAAACGCTTTACTCAACTCCCCGCATTTACTCTCATCAATACCTGGTCCTTCATCTTCAACTGCCATGACCGCTCCGTCATCTTCTTGCAGCTTAATCATAATg                                                                                              >  1:4964/1‑251 (MQ=255)
                                                                                                                                              tcttcAGCCGTACCCAAGCCCGCGTGCCGGAAGTCTCTTGCCGGTTTTGCAGGCAAAACTGCCCGTGATGCAAATGCGTAATGCGGCTGACAATACTTAATCCCACACAAATCCCGCCATAACGGCCGTCCATACGTACAAAAGCTTTACTCAACTCCCCGCATTTCCTCTCATCAATACCTGGTCCTTCATCTTCAACTGCCATGACCGCTCCGTCATCTTCTTGCAGCTTAATCATAATGTTGCTGCCt                                                                                     >  2:189646/1‑251 (MQ=255)
                                                                                                                                               cttcAGCCGGACCCAGTCCCGAGTGCCGGAAGTATCTCGTCGGTTTTGCCGGACAAACTGCCCGTGATGCAACTGTGTAATGCGGCTGACAATACTTAATCCCAGACCAATCCCGCAATAACGGCTGTCCATACGCACAAACGCTTTACTCAACTCCCCGCACTTACTATCATCACTACATGGTACTTCATCTTCACTTGCCATGACCTCTCCGTCATCTTCTTGCAGCTTACTAAGAATGTagctgcatt                                                                                    >  2:83459/1‑248 (MQ=255)
                                                                                                                                                                     gTGCCGGAAGTCTCTTGCCGGTTTTGCAGGAAAAACTGCCCGTGATGCAACTGTGTAATGCGGCTGACAATACTTAATCCCAGACCAATCCCGCCATAACGGCTGTCCATACGTACAAACGCTTTACTCAACTCCCCGCATTTACTCTCATCAATACCTGGTCCTTCATCTTCAACTGCCATGACCGCTCCGTCATCTTCTTGCAGCTTAATCATAATGTTGCTGCCTTGCGGGCTGTAACGATGGGCGtt                                                              >  1:169968/1‑251 (MQ=255)
                                                                                                                                                                          ggAAGTCTCTTGCCGGTTTTGCAGGAAAAACTGCCCGTGATGCAACTGTGTAATGCGGCTGACAATACTTAATCCCAGACCAATCCCGCCATACCGGCTGTACATACGTACAAACGCTTTACTCAACTCACCGCATTTACTCTCATCAATACCTGGTCCTTAATCTTCAACTGCCATAACCGCTCCGTCATCTTCTTGCAGCTTAATCATAATGTTGCTGCCTTGCGGGCTGTAACGATGGGCGTTTTCt                                                          >  1:284480/1‑250 (MQ=255)
                                                                                                                                                                               tctcTTGCCGGTTTTGCAGGAAAAACTGCCCGTGATGCAACTGTGTAATGCGGCTGACAATACTTAATCCCAGACCAATCCCGCCATAACGGCTGTCCATACGTACAAACGCTTTACTCAACTCCCCGCATTTACTCTCATCAATACCTGGTCCTTCATCTTCAACTGCCATGACCGCTCCGTCATCTTCTTGCAGCTTAATCATAATGTTGCTGCCTTGCGGGCTGTAACGATGGGCGTTTTCTACCAgg                                                    >  1:110765/1‑251 (MQ=255)
                                                                                                                                                                                                                   acaacttgaaaagggggcaaaCAATTCTAAACCCCACACCAACCCCGCCATAACGGCTGTCCATACGTACAAACGCTTTACTCAACTCCCCGCATTTACTCTCATCAATACCTGGTCCTTCATCTTCAACTGCCATGACCGCTCCGTCATCTTCTTGCAGCTTAATCATAATGTTGCTGCCTTGCGGGCTGTAACGATGGGCGTTTTCTACCAGGTTTCGCAATAACATCCGCAGCAGGGTTGCATCAccc                <  1:44846/231‑1 (MQ=255)
                                                                                                                                                                                                                      aCTGTGTAATGCGGCTGACAATACTTAATCCCAGACCAATCCCGCCATAACGGCTGTCCATACGTACAAACGCTTTACTCAACTCCCCGCATTTACTCTCATCAATACCTGGTCCTTCATCTTCAACTGCCATGACCGCTCCGTCATCTTCTTGCAGCTTAATCATAATGTTGCTGCCTTGCGGGCTGTAACGATGGGCGTTTTCTACCAGGTTTCGCAATAACATCCGCCGCAGGGTTGCATCACCCTg              >  2:201869/1‑250 (MQ=255)
                                                                                                                                                                                                                      aCTGTGTAATGCGGCTGACAATAATTAATCCCAGAACAATCCCGCCATAACGGCTGTCCATACGTACAAACGCTTTACTCAACTCCCCGCATTTACTCTCATCAATACCTGGTCCTTCATCTTCAACTGCCATGACCGCTCCGTCATCTTCTTGCAGCTTAATCATAATGTTGCTGCCTTGCGGGCTGTAACGATGGGCGTTTTCTACCAGGTTTCGCAATAACATCCGCAGCAGGGTTGCATCACCCTg              >  1:116748/1‑250 (MQ=255)
                                                                                                                                                                                                                                 cgtctaaCAATACTTATCCCCATCCCAATCCCGCCATAACGGCTGTCCATACGTACAAACGCTTTACTCAACTCCCCGCATTTACTCTCATCAATACCTGGTCCTTCATCTTCAACTGCCATGACCGCTCCGTCATCTTCTTGCAGCTTAATCATAATGTTGCTGCCTTGCGGGCTGTAACGATGGGCGTTTTCTACCAGGTTTCGCAATAACATCCGCAGCAGGGTTGCATCACCCTGAACGGTGATGTc  <  1:300760/245‑1 (MQ=255)
                                                                                                                                                                                                                                                  |                                                                                                                                                                                                                                         
TAATCTCTGACGCGCATACTCTCCTCCAGGTTAACGGAGGAGAGTGCAATGAAAAACCGTGTTTATGAAAGTTTAACTACCGTGTTCAGCGTGCTGGTGGTCAGCAGCTTTCTTTATATCTGGTTTGCCACGTACTGATCTTTCTTCAGCCGTACCCAGGCCCGCGTGCCGGAAGTCTCTTGCCGGTTTTGCAGGAAAAACTGCCCGTGATGCAACTGTGTAATGCGGCTGACAATACTTAACCCCAGACCAATCCCGCCATAACGGCTGTCCATACGTACAAACGCTTTACTCAACTCCCCGCATTTACTCTCATCAATACCTGGTCCTTCATCTTCAACTGCCATGACCGCTCCGTCATCTTCTTGCAGCTTAATCATAATGTTGCTGCCTTGCGGGCTGTAACGATGGGCGTTTTCTACCAGGTTTCGCAATAACATCCGCAGCAGGGTTGCATCACCCTGAACGGTGATGTC  >  NC_000913/4332068‑4332543

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