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,3100CT95.5% 61.1 / ‑4.6 22G321G (GGG→GGA) basSsensor histidine kinase BasS
Reads supporting (aligned to +/- strand):  ref base C (0/1);  new base T (20/1);  total (20/2)
Fisher's exact test for biased strand distribution p-value = 9.09e-02
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 8.56e-01

CTCTGACGCGCATACTCTCCTCCAGGTTAACGGAGGAGAGTGCAATGAAAAACCGTGTTTATGAAAGTTTAACTACCGTGTTCAGCGTGCTGGTGGTCAGCAGCTTTCTTTATATCTGGTTTGCCACGTACTGATCTTTCTTCAGCCGTACCCAGGCCCGCGTGCCGGAAGTCTCTTGCCGGTTTTGCAGGAAAAACTGCCCGTGATGCAACTGTGTAATGCGGCTGACAATACTTAACCCCAGACCAATCCCGCCATAACGGCTGTCCATACG‑TACAAACGCTTTACTCAACTCCCCGCATTTACTCTCATCAATACCTGGTCCTTCATCTTCAACTGCCATGACCGCTCCGTCATCTTCTTGCAGCTTAATCATAATGTTGCTGCCTTGCGGGCTGTAACGATGGGCGTTTTCTACCAGGTTTCGCAATAACATCCGCAGCAGGGTTGCATCACCCTGAACGGTGATGTCGGCGGCGC  >  NC_000913/4332072‑4332551
                                                                                                                                                                                                                                              |                                                                                                                                                                                                                                                  
ctctGACGCGCATACTCTACTCCAGGTTAACGGAGGAGAGTGCAATGAAAAACCGTGTTTATGAAAGTTTAACTACCGTGTTCAGCGTGCTGGTGGTCAGCAGCTTTCTTTATATCTGGTTTGCCACGTACTGATCTTTCTTCAGCCGTACCCAGTCCCGCGTGCCGGAAGTCTCTTGCCGGTTTTGCAGGAAAAACTGCCCGTGATGCAACTGTGTAATGCGGCTGACAATACTTAATCCCAGACCAATc                                                                                                                                                                                                                                        >  1:348296/1‑251 (MQ=255)
    gACGCGCATACTCTCCTCCAGGTTAACGGAGGAGAGTGCAATGAAAAACCGTGTTTATGAAAGTTTAACTACCGTGTTCAGCGTGCTGGTGGTCAGCAGCTTTCTTTATATCTGGTTTGCCACGTACTGATCTTTCTTCAGCCGTACCCAGGCCCGCGTGCCGGAAGTCTCTTGCCGGTTTTGCAGGAAAAACTGCCCGTGATGCAACTGTGTAATGCGGCTGACAATACTTAATCCCAGACCAATCCCg                                                                                                                                                                                                                                     <  2:181223/250‑1 (MQ=255)
                   ctccAGGTTAACGGAGGAGAGTGCAATGAAAAACCGTGTTTATGAAAGTTTAACTACCGTGTTCAGCGTGCTGGTGGTCAGCAGCTTTCTTTATATCTGGTTTGCCACGTACTGATCTTTCTTCAGCCGTACCCAGGCCCGCGTGCCGGAAGTCTCTTGCCGGTTTTGCAGGACAAACTGCCCGTGATGCAACTGTGTAATGCGGCTGACAATACTTAATCCCAGACCAATCCCGCCATAACGGCTGTcc                                                                                                                                                                                                                      >  1:379067/1‑250 (MQ=255)
                      cAGGTTAACGGAGGAGAGTGCAATGAAAAACCGTGTTTATGAAAGCTTAACTACCGTGTTCAGCGTGCTGGTCGTCAGCAGATTTCTTTATATCTGGTTTGCCACGTACTGATCTTTCTTCAGCAGTACCCAGGCCCGCGTGCCGGAAGTCTCTTGCCGGTTTTGCAGGAAAAACTGCCCGCGATGCAACTGTGTAATGCGGCTGACAATACTTAATCCCAGACCAATCCCGCCATAACGGCt                                                                                                                                                                                                                          >  1:436011/1‑243 (MQ=255)
                           tAACGGAGGAGAGTGCAATGAAAAACCGTGTTTATGAAAGTTTAACTACCGTGTTCAGCGTGCTGGTGGTCAGCAGCTTTCTTTATATCTGGTTTGCCACGTACTGATCTTTCTTCAGCCGTACCCAGGCCCGCGTGCCGGAAGTCTCTTGCCGGTTTTGCAGGAAAAACTGCCCGTGATGCAACTGTGTAATGCGGCTGACAATACTTAATCCCAGACCAATCCCGCCATAACGGCTGTCCATACG‑TAc                                                                                                                                                                                                             >  1:183859/1‑250 (MQ=255)
                               ggaggaGAGTGCAATGAAAACCCGTGTTTATGAAAGTTTAACTACCGTGTTCAGCGTGCTGGTGGTCAGCAGCTTTCTTTATATCTGGTTTGCCACGTACTGATCTTTCTTCAGCCGTACCCAGGCCCGCGTGCCGGAAGTCTCTTGCCGGTTTTGCAGGAAAAACTGCCCGTGATGCAACTGTGTAATGCGGCTGACAATACTTAATCCCAGACCAATCCCGCCATAACGGCTGTCCATACG‑TACAAACg                                                                                                                                                                                                        >  2:207262/1‑251 (MQ=255)
                                      agTGCAATGAAAAACCGAGTTTATGAAAGTTTAACTACCGTGTTCAGCGTGCTGGTGGTCAGCAGCTTTCTTTATATCTGGTTTGTCACGTACTGATCTTTCTTCAGCCGTACCCAGGCCCGCGTGCCGGAAGTCTCTTGCCGGTTTTGCAGGAAAAACTGCCCGTGATGCAACTGTGTCATGCGGCTGACAATACTTAATCCCAGACCAATCCCGCCATAACGGCTGTCCATACG‑TACAAACGCTTTACt                                                                                                                                                                                                 >  2:308847/1‑251 (MQ=255)
                                              gAAAAACCGTGTTTATGAAAGTTTAACTACCGTGTTCAGCGTGCTGGTGGTCAGCAGCTTTCTTTATATCTGGTTTGCCACGTACTGATCTTTCTTCAGCCGTACCCAGGCCCGCGTGCCGGAAGTCTCTTGCCGGTTTTGCAGGAAAAACTGCCCGTGATGCAACTGTGTAATGCGGCTGACAATACTTAATCCCAGACCAATCCCGCCATAACGGCTGTCCATACG‑TACAAACGCTTTACTCAACTccc                                                                                                                                                                                         >  1:278543/1‑251 (MQ=255)
                                                                aaGTTTCACTACAGTGTTCAGCGTGCTGGTGGTCAGCAGCGTTCTTTATATCTGGTCTGCCACGTACTGATCTTTCTTCAGCCGTACCCAGGCCCGCGTGCCGGAAGTCTCTTGCCGGTTTTGCAGGAAACACTGCCCGTGATGCAGCTGTGTAATGCGGCTGACAATACTTAATCCCACACCATTACAGCCATAACGGCGGTCCATACGATA‑AACCTCTTTACTCACCTCCACGCagttcatctactca                                                                                                                                                                        >  2:145989/1‑237 (MQ=255)
                                                                  gTTTAACTCCCGAGTTCAGCGTGCTGGTGGTCAGCAGCTTTCTTTATATCTGGTTTGCCACGTACTGATCTTTCTTCAGCCGTACCCAGGCCAGCGTGCCGGAAGTCTCTTGCCGGTTTTGCAGGAAAACCTGCCCGTGATGCAACTGTGTAATGCGGCTGACACTCCTTAATCCCCGACCAATCCCGCCATAACGGCTGTCCATACG‑TACAAACGCTATCCTCCACTCCCCACATTTACTCtcatcatt                                                                                                                                                                      >  2:431765/1‑248 (MQ=255)
                                                                     tAACTACCGTGTTCAGCGTGCTGGTGGTCAGCAGCTTTCTTTATATCTGGTTTGCCACGTACTGATCTTTCTTCAGCCGTACCCAGGCCCGCGTGCCGGACGTCTCTTGCAGGTTTTGCAGGAAAAACTGCCCGTGATGCAACTGTGTAATGCGGCTGACAATACTTAATCCCAGACCAATCCCGCCATAACGGCTGTCCATACG‑TACAAACGCTTTACTCAACTCCCCGCATTTACTCTCATCAATACCt                                                                                                                                                                  >  2:291624/1‑251 (MQ=255)
                                                                        cTACCGTGTTCAGCGTGCTGGTGGTCAGCAGCTTTCTTTATATCTGGTTTGCCACGTACTGATCTTTCTTCAGCCGTACCCAGGCCCGCGTGCCGGAAGTCTCTTGCCGGTTTTGCAGGAAAAACTGCCCGTGATGCAACTGTGTAATGCGGCTGACAATACTTAATCCCAGACCAATCCCGCCATAACGGCTGTCCATACG‑TACAAACGCTTTACTCAACTCCCCGCATTTAc                                                                                                                                                                                >  1:285229/1‑234 (MQ=255)
                                                                                     cGTGCTAGCGGTCCGCAGCGTTCTTTACATCTGGTTTGCCACGTACTGATCTTTCTTCAGCCGTCCCCAGGCCCGCGTGCCGCAAGTCTCTTGCCGGTTTTGCAGGAAAAACTGCCCGTGATGCAACTGTGTAATGCGGCTGACAATACTTAATCCCAGACCCATCCCGCCATACCGGCTGTCCATACG‑TACAAACGCTTTACTCAACTCCCCGCATTTACTCTCATCAATCCCTGGTCCTTCATCTTCa                                                                                                                                                   >  2:406055/1‑250 (MQ=255)
                                                                                                                   cTGGTTTGCCACGTACTGCTCTTTCTTCAGCCGTACCCAGGCCCGCGTGCCGGAAGTCTCTTGCCGGTTTTGCAGGAAAAACTGCCCGTGATGCAACTGTGTAATGCGGCTGACAATACTTAATCCCAGACCAATCCCGCCATAACGGCTGTCCATACG‑TACAAACGCTTTACTCAACTCCCCGCATTTACTCTCATCAATACCTGGTCCTTCATCTTCAACTGCCATGACCGCTCCGTCATCTTCTTGc                                                                                                                     >  1:361712/1‑250 (MQ=255)
                                                                                                                               gTACTGATCTTTCTTCAGCCGTACCCAGGCCCGCGTGCCGGAAGTCTCTTGCCGGTTTTGCAGGAAAAACTGCCCGTGATGCAACTGTGTAATGCGGCTGACAATACTTAATCCCAGACCAATCCCGCCATAACGGCTGTCCATACG‑TACAAACGCTTTACTCAACTCCCCGCATTTACTCTCATCAATACCTGGTCCTTCATCTTCAACTGCCATGACCGCTCCGTCATCTTCTTGCAGCTTAATCATa                                                                                                         >  2:30027/1‑250 (MQ=255)
                                                                                                                                                 ccGTACCCAGGCCCGCGTGCCGGAAGTCTCTTGCCGGTTTTGCAGGAAAAACTGCCCGTGATGCACCTGTGTAATGCGGCTGACAATACTTAATCCCAGCCCAATCCCGCCATAACGGCTGTCCATACG‑TACAAACGCTTTACTCAACTCCCCGCCTTTACTCTCCTCAATCCCTGGTCCTTCATCTTCAACTGCCATGCCCGCTCCGTCATCTTCTTGCAGCTTAATCATAATGTTGCTGCCTTGCGGGc                                                                                      >  2:194213/1‑251 (MQ=255)
                                                                                                                                                         aGGCCCGCGTGCCGGAAGTCTCTTGCCGGTTTTGCAGGACCAACTGCCCGTGATGCGACTGTGTAATGCGGCTGACAATACTTAATCCCAGACCAATCCCGCCATAACGGCTGTCCATACG‑TACACACGCTTTACTCAACTCCCCGCAt                                                                                                                                                                                    >  2:36012/1‑149 (MQ=255)
                                                                                                                                                                    ccGGAAGTCTCGTCCCGGATTTGAAGGCAAAACTGCCCGTGAAGCCACTGTGTAATGCGGCTGACCATACTTACTCCCAGACCAATCCCGCCATAACGGCTGTCCATACG‑TACAAACGCTTTCCTCAACTCCCCCCATTTACTCTCATCAATAACAGGTCCTTCATCTTCACCTGCCATGACCGCTCCGTCATCTTCTTGCAGCTTCATCCTAATGTTGCTCCCTTGCGGGCTGTAACGATGGGCGtttt                                                                    >  2:413206/1‑250 (MQ=255)
                                                                                                                                                                                                        ctgTGATGCAACTGTGTAATCCGGCTGACACTACTTAATCCCAGACCAATCCAGCCATAACGGCTGTCCATACG‑TACAAACGCTTTACTCACAACCCCGCATTTACTCTCATCAATACCTGGTCCTGCATCTTCAACTGCAATGACCCCCCAGTCATCTTCTTGCAGCTTACTCATAATGTTGCTGCCTTGCGGGCTGTAACGATGGGCGTTTTCTACCAGGTTTCGCAATAACATGCGAAGCAgggatac                               >  2:118769/3‑247 (MQ=255)
                                                                                                                                                                                                                        tAATGCGGCTGACAATACTTAATCCCAGACCAATCCCGCCATAACGGCTGTCCATACG‑TACAAACGCTTTACTCAACTCCCCGCATTTACTCTCATCAATACCTGGTCCTTCATCTTCAACTGCCATGACCGCTCCGTCATCTTCTTGCAGCTTAATCATAATGTTGCTGCCTTGCGGGCTGTAACGATGGGCGTTTTCTACCAGGTTTCGCAATAACATCCGCAGCAGGGTTGCATCACCCTGAACGGTg               >  2:352074/1‑251 (MQ=255)
                                                                                                                                                                                                                              ggtgtccattctctacccccACACCAACCCCGCCATAACGGCTGTCCATACG‑TACAAACGTTTTACTCAACTCCCCGCATTTACTCTCATCAATACCTGGTCCTTCATCTTCAACTGGCATGACCGCTCCGTCATCTTCTTGCAGCTTAAGCATAATGTTGCTGCCTTGCGGGCTGTAACGATGGGCGTTTTCTACCAGGTTTCGCAATAACATCCGCAGCAGGGTTGCATCACCCTGAACGGTGATGTcg         <  1:78231/235‑1 (MQ=255)
                                                                                                                                                                                                                                 ggacataacgcaaCCCCACACCATCCCCGCCATAACGGCTGTCCATACG‑TACAAACGCTTTACTCAACTCCCCGCATTTACTCTCATCAATACCTGGTCCTTCATCTTCAACTGCCATGACCGCTCCGTCATCTTCTTGCAGCTTAATCATAATGTTGCTGCCTTGCGGGCTGTAACGATGGGCGTTTTCTACCAGGTTTCGCAAGAACATCCGCAGCAGGGTTGCATCACCCTGAACGGTGATGTcggcg      <  1:194213/240‑1 (MQ=255)
                                                                                                                                                                                                                                     aaTACTTAATCCCAGACCAATCCCGCCATAACGGCTGTCCATACG‑TACAAACGCTTTACTCAACTCCCCGCATTTACTCTCATCAATACCTGGTCCTTCATCTTCAACTGCCATGACCGCTCCGTCATCTTCTTGCAGCTTAATCATAATGTTGCTGCCTTGCGGGCTGTAACGATGGGCGTTTTCTACCAGGTTTCGCAATAACATCCGCAGCAGGGTTGCATCACCCTGAACGGTGATGTCGGCGgcgc  >  1:132436/1‑251 (MQ=255)
                                                                                                                                                                                                                                              |                                                                                                                                                                                                                                                  
CTCTGACGCGCATACTCTCCTCCAGGTTAACGGAGGAGAGTGCAATGAAAAACCGTGTTTATGAAAGTTTAACTACCGTGTTCAGCGTGCTGGTGGTCAGCAGCTTTCTTTATATCTGGTTTGCCACGTACTGATCTTTCTTCAGCCGTACCCAGGCCCGCGTGCCGGAAGTCTCTTGCCGGTTTTGCAGGAAAAACTGCCCGTGATGCAACTGTGTAATGCGGCTGACAATACTTAACCCCAGACCAATCCCGCCATAACGGCTGTCCATACG‑TACAAACGCTTTACTCAACTCCCCGCATTTACTCTCATCAATACCTGGTCCTTCATCTTCAACTGCCATGACCGCTCCGTCATCTTCTTGCAGCTTAATCATAATGTTGCTGCCTTGCGGGCTGTAACGATGGGCGTTTTCTACCAGGTTTCGCAATAACATCCGCAGCAGGGTTGCATCACCCTGAACGGTGATGTCGGCGGCGC  >  NC_000913/4332072‑4332551

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