Predicted mutation
evidence seq id position mutation annotation gene description
RA NC_000913 2,484,827 C→T L152F (CTT→TTT)  evgS → sensor histidine kinase EvgS

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009132,484,8270CT95.5% 60.9 / ‑4.5 23L152F (CTT→TTT) evgSsensor histidine kinase EvgS
Reads supporting (aligned to +/- strand):  ref base C (0/1);  major base T (7/14);  minor base G (0/1);  total (7/16)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 6.49e-01

GTTGCATACCGATTCGCAGCAACGGGTTCGTGGTATTAATGCTGATTATTTAAATCTTTTAAAAAGAGCGTTAAATATCAAATTAACACTCCGGGAATACGCAGATCATCAAAAAGCAATGGACGCGCTTGCAGAAGGTGAAGTCGATATAGTGTTATCACATTTAGTTACTTCGCCGCCTCTTAATAATGACATTGCTGCAACCAAACCATTGATAATTACCTTTCCGGCGCTGGTAACCACCCTTCACGACTCAATGCGACCGCTTACCTCACCAAAACCAGTAAATATTGCTCGGGTAGCAAATTACCCCCCAGACGAGGTAATTCATCAATCATTTCCAAAAGCAACAATTATCTCTTTTACAAATTTATATCAGGCATTAGCATCCGTCTCAGCTGGGCACAATGATTACTTTATTGGTAGTAACATCATTACCAGCAGTATGATTTCCCGCTATTTCACTCACT  >  NC_000913/2484583‑2485052
                                                                                                                                                                                                                                                    |                                                                                                                                                                                                                                  
gtttctcatcGATTTGCAGAAACGTGTTCGTTGTCTTAATGCTGATTATTTAAATCTTTTTAAAATAGCGTTAAATATCACATTAACAATTCGGGAATACGCTGATCATCAAAAAGCAATGGACGCGCTTGTAGAAGGTGAAGTTGATATAGTGTTATCACATTTAGTTACTTCGCCGCCTCTTAATACTGACATTGCTGCAACCATATCATCGATAAGTACTGCTCCGGCGCAGGTATCGACCTTTCAc                                                                                                                                                                                                                              <  2:579334/241‑1 (MQ=255)
                  gcaACGGGTTCGTGATATTAATGCTGATTATTTAAATCTTTTAAAAAGAGCGTTAAATATCAAATTAACACTCCGGGAATACGCAGATCATCAAAAAGCAATGGACGCGCTTGCAGAAGGTGAAGTCGATATAGTGTTATCACATTTAGTTACTTCGCCGCCTCTTAATAATGACATTGCTGCAACCAAACCATTGATAATTTCCTTTCCGGCGCTGGTAACCACCTTTCACGACTCAATGCGACCGCTTa                                                                                                                                                                                                           <  2:231075/251‑1 (MQ=255)
                      cGGGTTCGTGGTATTAATGCTGATTATTTAAATCTTTTAAAAAGAGCGTTCAATATCAAATTAACCCTCCGGGAATACGCAGATCATCAAAAAGCAATGGACGCGCTTGCAGAAGGTGACGTCGATATAGTGTTATCACATTTAGTTACTTCGCCGCCTCTTAATAATGACATTGCTGCAACCAAACCATTGATAATTCCCTTTCCGGCGCTGGTAACCACCTTTCACGCCTCAATGCGACCGCTTCCCTc                                                                                                                                                                                                       >  2:250020/1‑251 (MQ=255)
                                           gATTATTTAAATCTTTTAAAAAGAGCGTTAAATATCAAATTAACACTCCGGGAATACGCAGATCATCAAAAAGCAATGGACGCGCTTGCAGAGGGCGAAGTCGATATAGTGTTATCACATTTAGTTTCTTCGCTGCCTCTTAATAATGCGAGTGCGGCAACCAAACCATGGATAATTTCCTATACGGCGCTGGTAACGAGCGTTCGCGACTCAATGCGACCGCTTACCTCACCAAAACCAGTAAATATTg                                                                                                                                                                                   <  1:274543/250‑1 (MQ=255)
                                                 ttAAATCTTTTAAAAAGAGCGTTAAATATCAAATTAACACTCCGGGAATACGCAGATCATCAAAAAGCAATGGACGCGCTTGCAGAAGGTGAAGTCGATATAGTGTTATCACATTTAGTTACTTCGCCGCCTCTTAATAATGACATTGCTGCAACCAAACCATTGATAATTACCTTTCCGGCGCTGGTAACCACCTTTCACGACTCAATGCGACCGCTTACATCACCAAAACAAGTAAATATTGCTCGGGt                                                                                                                                                                            >  2:489567/1‑251 (MQ=255)
                                                      tCTTTTAAAAAGAGCGTTAAATATCAAATTAACACTCCGGGAATACGAAGCTCATCAAAAAGCAATGGAAGCGCTTGCAGAAGGTGAAGTCGATATAGTGTTATCACATTTAGTTACTTCGCCGCCTCTTAATAATGACATTGCTGCAACCAAACCATTGATATTTACCTTTCCGGCGCTGGTAACCACCTTTCACGACTCAATGCGACCGCTTACCTCACCAAAACCAGTAAATATTGATCGGGTAGcc                                                                                                                                                                        >  2:536160/1‑249 (MQ=255)
                                                      tCTTTTAAAAAGAGCGTTAAAAATCAAATTAAAACTCCGGGAATACGCAGATCATCAAAAAGCAATGGACGCGCTTACAGAAGGTGAAGTCGATATAGTGTTATCACATTTAGTTACTTCGCCGCCTCTTAATAATGACATTGATGCAAACAAACCATTGATAATTACCTTTCCGGCGCTGGTAACCACCTTTCACGACTCAATGCTTCAGCTTACCACACCAAAACTAGTAAATATTGCTCGGGTAGCa                                                                                                                                                                        >  2:316381/1‑250 (MQ=255)
                                                            aaaaaGAGCGGTAAATATCAAATTAACACTCCGGGAATACGCAGATCATCAAAAAGCAATGTACGCGCTTGCAGAAGGTGAAGTCGATATAGTGTTATCACATTTAGTTTCTTCGCCGCCTCTTAATAATGACATTGCTGCAACCAAACCATTGATAATTACCTGTCCGGCGCTGGTAACCCCCTGTCTCGACTCAATGCGACCTCTTTCCTCACCAAAACCAGTGAATATTGCTCGTTTAGCAAATTAcc                                                                                                                                                                 <  2:573086/251‑1 (MQ=255)
                                                              aaaGAGCGTTAAATATCAAATTAACACTCCGGGAATACGCAGATCATCAAAAAGCAATGGACGCGCTTGCAGAAGGTGAAGTCGATATAGTGTTATCACATTTAGTTACTTCGCCGCCTCTTAATAATGACATTGCTGCAACCAAACCATTGATAATTACCTTGCCGGCGCTGGTAACCACCTTTCACGACTCAATGCGACCGCTTACCTCACCAAAACCAGTAAATATTGCTCGGGTAGCAAATTAcccc                                                                                                                                                               <  1:536160/251‑1 (MQ=255)
                                                              aaaGAGCGGTAAACATCAAATTAACACTCAAGGAATCCTCAGATCATCAAAAAGCAATGGAAGCGCTTGCAGAAGGTGAAGACGATATAGTGTTATCACATTTAGTCTCGTCCCCGCCTCTTAATAATGATATTGCTGCAACCAAACCATTGATAATTACCTTTCCGGCTCTGGGAACCACCTTTCTCGACTCAATGCTACCGCTTACCTCACCAAAACCAGTAATTATTGCTCGCGTAGCAAAtttcacc                                                                                                                                                               <  2:439421/251‑6 (MQ=255)
                                                                        aaaTATCAAATTAACACTCCGGGAATACGCAGATCATCAAAAAGCAATCGACGCGCTTGCAGAAGGTGAAGTCGATATAGTGTTATCACATTTAGTTACTTCGCCGCCTCTTAATAATGACATTGCTGCAACCAAACCATTGATAATTACCTTTCCGGCGCTGGTAACCACCTTTCACGACTCAATGCGACCGCTTACCTCACCAAAACCAGTAAATATTGCTCGGGTAGCAAATTACCCCCCAGACGAg                                                                                                                                                      <  1:402662/250‑1 (MQ=255)
                                                                            atCAAATTAACACTCCGGGAATACGCAGAGCATCAAAAAGCAATGGACGCGCTTGCAGAAGGTGAAGTCGATATAGTGTTATCACATTTAGTTACTTCGCCGCCTCTTAATAATGACATTGCTGCAACCAAACCATTGATAATTACCTTTCCGGCGCTGGTAACCACCTTTCACGACTCAATGCGACCGCTTACCTCACCAAAACCAGTAAATATTGCTCGGGTAGCAAATTACCCCCCAGACGATGTAAt                                                                                                                                                 <  2:551847/251‑1 (MQ=255)
                                                                                aaTTAACACTCCGGGAATACGCAGATCATCAAAAAGCAATGGACGCGCTTGCAGAAGGTGAAGTCGATATAGTGTTATCACATTTAGTTAATTCGCCGCCTCTTAATAATGACATTGCTGCAACCAAACCATTGATAATTACCTTTCCGGCGCTGGTAACCACCTTTCACGACTCAATGCGACCGCTTACCTCACCAAAACCAGTAAATATTGCTCGGGTAGCAAATTACCCCCCATACGAGGTAATtcat                                                                                                                                             <  2:97389/251‑1 (MQ=255)
                                                                                 aTTAACACTCCGGGAATACGCAGATCATCAAAAAGCAATGGACGCGCTTGCAGAAGGTGAAGTCGATATAGTGTTATCACATTTAGTTACTTCGCCGCCTCTTAATAATGACATTGCTGCAACCAAACCATTGATAATTACCTTTCCGGCGCTGGTAACCACCTTTCACGACTCAATGCGACCGCTTACCTCACCAAAACCAGTAAATATTGCTCGGGTAGCAAATTACCCCCCAGACGAGGTAATtcatc                                                                                                                                            >  1:302031/1‑251 (MQ=255)
                                                                                                       gATCATCAAAAAGCAATGGACGCGCTTGCAGAAGGTGAAGTCGATATAGTGTTATCACATTTAGTTACTTCGCCGCCTCTTAATAATGACATTGCTGCAACCAAACCATTGATAATTACCTTTCCGGCGCTGGTAACCACCTTTCACGACTCAATGCGACCGCTTACCTCACCAAAACCAGTAAATATTGCTCGGGTAGCAAATTACCCCCCAGACGATGTAATTCATCAATCATTTCCAAAAGCAACAAt                                                                                                                      <  1:489567/251‑1 (MQ=255)
                                                                                                       gATCATCAAAAAGCAATGGACGCGCTTGCAGAAGGTGAAGTCGATATAGTGTTATCACATTTAGTTACTTCGCCGCCTCTTAATAATGACATTGCTGCAACCAAACCATTGATAATTACCTTTCCGGCGCTGGTAACCACCTTTCACGACTCAATGCGACCGCTTACCTCACCAAAACCAGTAAATATTGCTCGGGTAGCAAATTACCCCCCAGACGAGGTAATTCATCAATCATTTCCAAAAGCAACAAt                                                                                                                      >  1:580566/1‑251 (MQ=255)
                                                                                                                  aGCAATGGACGCGCTTGCAGAAGGTGAAGTCGATATAGTGTTATCACATTTAGTTACTTCGCCGCCTCTTAATAATGACATTGCTGCAACCAAACCATTGATAATTACCTTTCCGGCGCTGGTAACCACCTTTCACGACTCAATGCGACCGCTTACCTCACCAAAACCAGTAAATATTGCTCGGGTAGCAAATTACCCCCCAGACGAGGTAATTCATCAATCATTTCCAAAAGCAACAATTATCTCTTTTa                                                                                                           <  1:66197/251‑1 (MQ=255)
                                                                                                                          aCGCGCTTGCAGAAGGTGAAGTCGATATAGTGTTATCACATTTAGTTACTTCGCCGCCTCTTAATAATGACATTGCTGCAACCAAACCATTGATAATTACCTTTCCGGCGCTGGTAACCACCTTTCACGACTCAATGCGACCGCTTACCTCACCAAAACCAGTAAATATTGCTCGGGTAGCAAATTACCCCCCAGACGAGGTAATTCATCAATCATTTCCAAAAGCAACAATTATCTCTTTTACAAATTta                                                                                                   <  1:8430/251‑1 (MQ=255)
                                                                                                                                                                                                 aTTGCTGCAACCAAGCCATTGATAATTACCTTTCCGGCGCTGGTAACCACCTTTCACGACTCAATGCGACCGCTTACCTCACCAAAACCAGTAAATATTGCTCGGGTAGCAAATTACCCCCCAGACGAGGTAATTCATCAATCATTTCCAAAAGCAACAATTATCTCTTTTGCAAATTTATATCAGGCATTAGCATCCGGCTCAGCTGGGCACAATGATTACTTTATTGGTAGTAACATCATTACcagcag                            <  1:35140/251‑1 (MQ=255)
                                                                                                                                                                                                        cAACCAAACCAATGATAATTTCCTTTCCGGCGCTGGTAACCACCTGTCACGACTCAATGCGACCGCCTACCGCACCAAAACCAGTAAATATTGCTCGGGTAGCAAATTACCCCCCAGACGAGGTAATTCATCAATCATTTCCAAAAGCAACAATTATCTCTTTTACAAATTTATATCAGGCATTAGCATCCGTCTCAGCTGGGCACAATGATTACTTTATTGGTAGTAACATCATTACCAGCAGTATGAtt                     <  2:580566/251‑1 (MQ=255)
                                                                                                                                                                                                              aaccaTTGATAATTACCTTTCCGGCGCTGGTAACCACCTTTCACGACTCAATGCGACCGCTTACCTCACCAAAACCAGTAAATATTGCTCGGGTAGCAAATTACCCCCCAGACGAGGTAATTCATCAATCATTTCCAAAAGCAACAATTATCTCTTTTACAAATTTATATCAGGCATTAGCATCCGTCTCAGCTGGGCACAATGATTACTTTATTGGTAGTAACATCATTACCAGCAGTATGATTTCCCGc               >  1:128705/1‑251 (MQ=255)
                                                                                                                                                                                                                        aaTTACCTTTCCGGCGCTGGTAACCACCTTTCACGACTCAATGCGACCGCTTACCTCACCAAAACCAGTAAATATTGCTCGGGTAGCAAATTACCCCCCAGACGAGGTAATTCATCAATCATTTCCAAAAGCAACAATTATCTCTTTTACAAATGTATATCAGGCATTAGCATCCGTCTCAGCTGGGCACAATGATTACTTTATTGGTAGTAACATCATTACCAGCAGTATGATTTCCCGCTATTtcactc     <  1:343885/251‑1 (MQ=255)
                                                                                                                                                                                                                            cccTTTCCGGCGCTGGTAACAATCTTTCACGAATCAATGCGACCGCTTACCTCAACAAAACCAGTAAATAGTGCTCGGGTAGCAAATTACCCCACAGATGAGGTAATTCATCAATCATTTCCAAAAGCAACAAGTAACGCTTTTACAAATTTATATCAGGCATTAGCATCCGTCTCAGCTGGGCACAATGATTACTTTAGTGGTAGTAACATCATTACCAGCAGTATGATTTCCCGCTATTTcactcact  <  2:181799/249‑1 (MQ=255)
                                                                                                                                                                                                                                                    |                                                                                                                                                                                                                                  
GTTGCATACCGATTCGCAGCAACGGGTTCGTGGTATTAATGCTGATTATTTAAATCTTTTAAAAAGAGCGTTAAATATCAAATTAACACTCCGGGAATACGCAGATCATCAAAAAGCAATGGACGCGCTTGCAGAAGGTGAAGTCGATATAGTGTTATCACATTTAGTTACTTCGCCGCCTCTTAATAATGACATTGCTGCAACCAAACCATTGATAATTACCTTTCCGGCGCTGGTAACCACCCTTCACGACTCAATGCGACCGCTTACCTCACCAAAACCAGTAAATATTGCTCGGGTAGCAAATTACCCCCCAGACGAGGTAATTCATCAATCATTTCCAAAAGCAACAATTATCTCTTTTACAAATTTATATCAGGCATTAGCATCCGTCTCAGCTGGGCACAATGATTACTTTATTGGTAGTAACATCATTACCAGCAGTATGATTTCCCGCTATTTCACTCACT  >  NC_000913/2484583‑2485052

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