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→GGAbasSsensor 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

CTCTGACGCGCATACTCTCCTCCAGGTTAACGGAGGAGAGTGCAATGAAAAACCGTGTTTATGAAAGTTTAACTACCGTGTTCAGCGTGCTGGTGGTCAGCAGCTTTCTTTATATCTGGTTTGCCACGTACTGATCTTTCTTCAGCCGTACCCAGGCCCGCGTGCCGGAAGTCTCTTGCCGGTTTTGCAGGAAAAACTGCCCGTGATGCAACTGTGTAATGCGGCTGACAATACTTAACCCCAGACCAATCCCGCCATAACGGCTGTCCATACGTACAAACGCTTTACTCAACTCCCCGCATTTACTCTCATCAATACCTGGTCCTTCATCTTCAACTGCCATGACCGCTCCGTCATCTTCTTGCAGCTTAATCATAATGTTGCTGCCTTGCGGGCTGTAACGATGGGCGTTTTCTACCAGGTTTCGCAATAACATCCGCAGCAGGGTTGCATCACCCTGAACGGTGATGTCGGCGGCGC  >  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)
                           tAACGGAGGAGAGTGCAATGAAAAACCGTGTTTATGAAAGTTTAACTACCGTGTTCAGCGTGCTGGTGGTCAGCAGCTTTCTTTATATCTGGTTTGCCACGTACTGATCTTTCTTCAGCCGTACCCAGGCCCGCGTGCCGGAAGTCTCTTGCCGGTTTTGCAGGAAAAACTGCCCGTGATGCAACTGTGTAATGCGGCTGACAATACTTAATCCCAGACCAATCCCGCCATAACGGCTGTCCATACGTAc                                                                                                                                                                                                             >  1:183859/1‑250 (MQ=255)
                               ggaggaGAGTGCAATGAAAACCCGTGTTTATGAAAGTTTAACTACCGTGTTCAGCGTGCTGGTGGTCAGCAGCTTTCTTTATATCTGGTTTGCCACGTACTGATCTTTCTTCAGCCGTACCCAGGCCCGCGTGCCGGAAGTCTCTTGCCGGTTTTGCAGGAAAAACTGCCCGTGATGCAACTGTGTAATGCGGCTGACAATACTTAATCCCAGACCAATCCCGCCATAACGGCTGTCCATACGTACAAACg                                                                                                                                                                                                        >  2:207262/1‑251 (MQ=255)
                                      agTGCAATGAAAAACCGAGTTTATGAAAGTTTAACTACCGTGTTCAGCGTGCTGGTGGTCAGCAGCTTTCTTTATATCTGGTTTGTCACGTACTGATCTTTCTTCAGCCGTACCCAGGCCCGCGTGCCGGAAGTCTCTTGCCGGTTTTGCAGGAAAAACTGCCCGTGATGCAACTGTGTCATGCGGCTGACAATACTTAATCCCAGACCAATCCCGCCATAACGGCTGTCCATACGTACAAACGCTTTACt                                                                                                                                                                                                 >  2:308847/1‑251 (MQ=255)
                                              gAAAAACCGTGTTTATGAAAGTTTAACTACCGTGTTCAGCGTGCTGGTGGTCAGCAGCTTTCTTTATATCTGGTTTGCCACGTACTGATCTTTCTTCAGCCGTACCCAGGCCCGCGTGCCGGAAGTCTCTTGCCGGTTTTGCAGGAAAAACTGCCCGTGATGCAACTGTGTAATGCGGCTGACAATACTTAATCCCAGACCAATCCCGCCATAACGGCTGTCCATACGTACAAACGCTTTACTCAACTccc                                                                                                                                                                                         >  1:278543/1‑251 (MQ=255)
                                                                aaGTTTCACTACAGTGTTCAGCGTGCTGGTGGTCAGCAGCGTTCTTTATATCTGGTCTGCCACGTACTGATCTTTCTTCAGCCGTACCCAGGCCCGCGTGCCGGAAGTCTCTTGCCGGTTTTGCAGGAAACACTGCCCGTGATGCAGCTGTGTAATGCGGCTGACAATACTTAATCCCACACCATTACAGCCATAACGGCGGTCCATACGATAAACCTCTTTACTCACCTCCACGCagttcatctactca                                                                                                                                                                        >  2:145989/1‑237 (MQ=255)
                                                                  gTTTAACTCCCGAGTTCAGCGTGCTGGTGGTCAGCAGCTTTCTTTATATCTGGTTTGCCACGTACTGATCTTTCTTCAGCCGTACCCAGGCCAGCGTGCCGGAAGTCTCTTGCCGGTTTTGCAGGAAAACCTGCCCGTGATGCAACTGTGTAATGCGGCTGACACTCCTTAATCCCCGACCAATCCCGCCATAACGGCTGTCCATACGTACAAACGCTATCCTCCACTCCCCACATTTACTCtcatcatt                                                                                                                                                                      >  2:431765/1‑248 (MQ=255)
                                                                     tAACTACCGTGTTCAGCGTGCTGGTGGTCAGCAGCTTTCTTTATATCTGGTTTGCCACGTACTGATCTTTCTTCAGCCGTACCCAGGCCCGCGTGCCGGACGTCTCTTGCAGGTTTTGCAGGAAAAACTGCCCGTGATGCAACTGTGTAATGCGGCTGACAATACTTAATCCCAGACCAATCCCGCCATAACGGCTGTCCATACGTACAAACGCTTTACTCAACTCCCCGCATTTACTCTCATCAATACCt                                                                                                                                                                  >  2:291624/1‑251 (MQ=255)
                                                                        cTACCGTGTTCAGCGTGCTGGTGGTCAGCAGCTTTCTTTATATCTGGTTTGCCACGTACTGATCTTTCTTCAGCCGTACCCAGGCCCGCGTGCCGGAAGTCTCTTGCCGGTTTTGCAGGAAAAACTGCCCGTGATGCAACTGTGTAATGCGGCTGACAATACTTAATCCCAGACCAATCCCGCCATAACGGCTGTCCATACGTACAAACGCTTTACTCAACTCCCCGCATTTAc                                                                                                                                                                                >  1:285229/1‑234 (MQ=255)
                                                                                     cGTGCTAGCGGTCCGCAGCGTTCTTTACATCTGGTTTGCCACGTACTGATCTTTCTTCAGCCGTCCCCAGGCCCGCGTGCCGCAAGTCTCTTGCCGGTTTTGCAGGAAAAACTGCCCGTGATGCAACTGTGTAATGCGGCTGACAATACTTAATCCCAGACCCATCCCGCCATACCGGCTGTCCATACGTACAAACGCTTTACTCAACTCCCCGCATTTACTCTCATCAATCCCTGGTCCTTCATCTTCa                                                                                                                                                   >  2:406055/1‑250 (MQ=255)
                                                                                                                   cTGGTTTGCCACGTACTGCTCTTTCTTCAGCCGTACCCAGGCCCGCGTGCCGGAAGTCTCTTGCCGGTTTTGCAGGAAAAACTGCCCGTGATGCAACTGTGTAATGCGGCTGACAATACTTAATCCCAGACCAATCCCGCCATAACGGCTGTCCATACGTACAAACGCTTTACTCAACTCCCCGCATTTACTCTCATCAATACCTGGTCCTTCATCTTCAACTGCCATGACCGCTCCGTCATCTTCTTGc                                                                                                                     >  1:361712/1‑250 (MQ=255)
                                                                                                                               gTACTGATCTTTCTTCAGCCGTACCCAGGCCCGCGTGCCGGAAGTCTCTTGCCGGTTTTGCAGGAAAAACTGCCCGTGATGCAACTGTGTAATGCGGCTGACAATACTTAATCCCAGACCAATCCCGCCATAACGGCTGTCCATACGTACAAACGCTTTACTCAACTCCCCGCATTTACTCTCATCAATACCTGGTCCTTCATCTTCAACTGCCATGACCGCTCCGTCATCTTCTTGCAGCTTAATCATa                                                                                                         >  2:30027/1‑250 (MQ=255)
                                                                                                                                                 ccGTACCCAGGCCCGCGTGCCGGAAGTCTCTTGCCGGTTTTGCAGGAAAAACTGCCCGTGATGCACCTGTGTAATGCGGCTGACAATACTTAATCCCAGCCCAATCCCGCCATAACGGCTGTCCATACGTACAAACGCTTTACTCAACTCCCCGCCTTTACTCTCCTCAATCCCTGGTCCTTCATCTTCAACTGCCATGCCCGCTCCGTCATCTTCTTGCAGCTTAATCATAATGTTGCTGCCTTGCGGGc                                                                                      >  2:194213/1‑251 (MQ=255)
                                                                                                                                                         aGGCCCGCGTGCCGGAAGTCTCTTGCCGGTTTTGCAGGACCAACTGCCCGTGATGCGACTGTGTAATGCGGCTGACAATACTTAATCCCAGACCAATCCCGCCATAACGGCTGTCCATACGTACACACGCTTTACTCAACTCCCCGCAt                                                                                                                                                                                    >  2:36012/1‑149 (MQ=255)
                                                                                                                                                                    ccGGAAGTCTCGTCCCGGATTTGAAGGCAAAACTGCCCGTGAAGCCACTGTGTAATGCGGCTGACCATACTTACTCCCAGACCAATCCCGCCATAACGGCTGTCCATACGTACAAACGCTTTCCTCAACTCCCCCCATTTACTCTCATCAATAACAGGTCCTTCATCTTCACCTGCCATGACCGCTCCGTCATCTTCTTGCAGCTTCATCCTAATGTTGCTCCCTTGCGGGCTGTAACGATGGGCGtttt                                                                    >  2:413206/1‑250 (MQ=255)
                                                                                                                                                                                                        ctgTGATGCAACTGTGTAATCCGGCTGACACTACTTAATCCCAGACCAATCCAGCCATAACGGCTGTCCATACGTACAAACGCTTTACTCACAACCCCGCATTTACTCTCATCAATACCTGGTCCTGCATCTTCAACTGCAATGACCCCCCAGTCATCTTCTTGCAGCTTACTCATAATGTTGCTGCCTTGCGGGCTGTAACGATGGGCGTTTTCTACCAGGTTTCGCAATAACATGCGAAGCAgggatac                               >  2:118769/3‑247 (MQ=255)
                                                                                                                                                                                                                        tAATGCGGCTGACAATACTTAATCCCAGACCAATCCCGCCATAACGGCTGTCCATACGTACAAACGCTTTACTCAACTCCCCGCATTTACTCTCATCAATACCTGGTCCTTCATCTTCAACTGCCATGACCGCTCCGTCATCTTCTTGCAGCTTAATCATAATGTTGCTGCCTTGCGGGCTGTAACGATGGGCGTTTTCTACCAGGTTTCGCAATAACATCCGCAGCAGGGTTGCATCACCCTGAACGGTg               >  2:352074/1‑251 (MQ=255)
                                                                                                                                                                                                                              ggtgtccattctctacccccACACCAACCCCGCCATAACGGCTGTCCATACGTACAAACGTTTTACTCAACTCCCCGCATTTACTCTCATCAATACCTGGTCCTTCATCTTCAACTGGCATGACCGCTCCGTCATCTTCTTGCAGCTTAAGCATAATGTTGCTGCCTTGCGGGCTGTAACGATGGGCGTTTTCTACCAGGTTTCGCAATAACATCCGCAGCAGGGTTGCATCACCCTGAACGGTGATGTcg         <  1:78231/235‑1 (MQ=255)
                                                                                                                                                                                                                                 ggacataacgcaaCCCCACACCATCCCCGCCATAACGGCTGTCCATACGTACAAACGCTTTACTCAACTCCCCGCATTTACTCTCATCAATACCTGGTCCTTCATCTTCAACTGCCATGACCGCTCCGTCATCTTCTTGCAGCTTAATCATAATGTTGCTGCCTTGCGGGCTGTAACGATGGGCGTTTTCTACCAGGTTTCGCAAGAACATCCGCAGCAGGGTTGCATCACCCTGAACGGTGATGTcggcg      <  1:194213/240‑1 (MQ=255)
                                                                                                                                                                                                                                     aaTACTTAATCCCAGACCAATCCCGCCATAACGGCTGTCCATACGTACAAACGCTTTACTCAACTCCCCGCATTTACTCTCATCAATACCTGGTCCTTCATCTTCAACTGCCATGACCGCTCCGTCATCTTCTTGCAGCTTAATCATAATGTTGCTGCCTTGCGGGCTGTAACGATGGGCGTTTTCTACCAGGTTTCGCAATAACATCCGCAGCAGGGTTGCATCACCCTGAACGGTGATGTCGGCGgcgc  >  1:132436/1‑251 (MQ=255)
                                                                                                                                                                                                                                              |                                                                                                                                                                                                                                                  
CTCTGACGCGCATACTCTCCTCCAGGTTAACGGAGGAGAGTGCAATGAAAAACCGTGTTTATGAAAGTTTAACTACCGTGTTCAGCGTGCTGGTGGTCAGCAGCTTTCTTTATATCTGGTTTGCCACGTACTGATCTTTCTTCAGCCGTACCCAGGCCCGCGTGCCGGAAGTCTCTTGCCGGTTTTGCAGGAAAAACTGCCCGTGATGCAACTGTGTAATGCGGCTGACAATACTTAACCCCAGACCAATCCCGCCATAACGGCTGTCCATACGTACAAACGCTTTACTCAACTCCCCGCATTTACTCTCATCAATACCTGGTCCTTCATCTTCAACTGCCATGACCGCTCCGTCATCTTCTTGCAGCTTAATCATAATGTTGCTGCCTTGCGGGCTGTAACGATGGGCGTTTTCTACCAGGTTTCGCAATAACATCCGCAGCAGGGTTGCATCACCCTGAACGGTGATGTCGGCGGCGC  >  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: