Predicted mutation
evidence seq id position mutation freq annotation gene description
RA CP000730 306,351 C→A 14.5% A119E (GCG→GAG)  lytS → sensor histidine kinase LytS

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*CP000730306,3510CA14.5% 90.5 / 6.6 29A119E (GCG→GAG) lytSsensor histidine kinase LytS
Reads supporting (aligned to +/- strand):  ref base C (11/13);  new base A (2/2);  total (13/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 = 8.98e-01

GTTTAGATGATGATGTATCTTTAGCTAACACACGTGTATTAACGATAGGTGTCGCAGGATTAGTTGGTGGCCCTTTTGTAGGTCTATTTGTTGGCGTTATTTCAGGTATTTTCAGAGTGTATATGGGTGGGGCGGATGCACAAGTTTATCTTATCTCATCTATATTTATTGGTATAATTGCTGGTTATTTTGGCTTACAAGCTCAAAGACGCAAGCGTTACCCGAGTATTGCGAAAAGTGCCATGATTGGAATTGTTATGGAAATGATTC  >  CP000730/306219‑306488
                                                                                                                                    |                                                                                                                                         
gTTTAGATGATGATGTATCTTTAGCTAACACACGTGTATTAACGATAGGTGTCGCAGGATTAGTTGGTGGCCCTTTTGTAGGTCTATTTGTTGGCGTTATTTCAGGTATTTTCAGAGTGTATATGGGTGGGGCGGATGcac                                                                                                                                   <  1:157755/141‑1 (MQ=255)
               tATCTTTAGCTAACACACGTGTATTAACGATAGGTGTCGCAGGATTAGTTGGTGGCCCTTTTGTAGGTCTATTTGTTGGCGTTATTTCAGGTATTTTCAGAGTGTATATGGGTGGGGCGGATGCACAAGTTtatcttatct                                                                                                                    <  2:50148/141‑1 (MQ=255)
                aTCTTTAGCTAACACACGTGTATTAACGATAGGTGTCGCAGGATTAGTTGGTGGCCCTTTTGTAGGTCTATTTGTTGGCGTTATTTCAGGTATTTTCAGAGTGTATATGGGTGGGGAGGATGCACAAGTTTATCTTAtctc                                                                                                                   >  2:230158/1‑141 (MQ=255)
                aTCTTTAGCTAACACACGTGTATTAACGATAGGTGTCGCAGGATTAGTTGGTGGCCCTTTTGTAGGTCTATTTGTTGGCGTTATTTCAGGTATTTTCAGAGTGTATATGGGTGGGGAGGATGCACAAGTTTATCTTAtctc                                                                                                                   >  2:185398/1‑141 (MQ=255)
                         tAACACACGTGTATTAACGATAGGTGTCGCAGGATTAGTTGGTGGCCCTTTTGTAGGTCTATTTGTTGGCGTTATTTCAGGTATTTTCAGAGTGTATATGGGTGGGGCGGATGCACAAGTTTATCTTATCTCATCTAtatt                                                                                                          <  1:16176/141‑1 (MQ=255)
                                   gtATTAACGATAGGTGTCGCAGGCTTAGTTGGTGGCCCTTTTGTAGGTCTATTTGTTGGCGTTATTTCAGGTATTTTCAGAGTGTATATGGGTGGGGCGGATGCACAAGTTTATCTTATCTCATCTATATTTATTGGtata                                                                                                >  2:21999/1‑141 (MQ=255)
                                     aTTAACGATAGGTGTCGCAGGATTAGTTGGTGGCCCTTTTGTAGGTCTATTTGTTGGCGTTATTTCAGGTATTTTCAGAGTGTATATGGGTGGGGCGGATGCACAAGTTTATCTTATCTCATCTATATTTATTGGTATAAt                                                                                              <  1:261001/141‑1 (MQ=255)
                                     aTTAACGATAGGTGTCGCAGGATTAGTTGGTGGCCCTTTTGTAGGTCTATTTGTTGGCGTTATTTCAGGTATTTTCAGAGTGTATATGGGTGGGGCGGATGCACAAGTTTATCTTATCTCATCTATATTTATTGGTATAAt                                                                                              <  2:5855/141‑1 (MQ=255)
                                                gtgtCGCAGGATTAGTTGGTGGCCCTATTGTAGGTCTATTTGTTGGCGTTATTTCAGGTATTTTCAGAGTGTATATGGGTGGGGCGGATGCACAAGTTTATCTTATCTCATCTATATTTATTGGTATAATTGCTGGTTAtt                                                                                   >  1:104415/1‑141 (MQ=255)
                                                            tAGTTGGTGGCCCTTTTGTAGGTCTATTTGTTGGCGTTATTTCAGGTATTTTCAGAGTGTATATGGGTGGGGCGGATGCACAAGTTTTTCTTATCTCATCTATATTTATTGGTATAATTGCTGGTTATTTTGGCTTACaat                                                                       >  1:91120/1‑140 (MQ=255)
                                                            tAGTTGGTGGCCCTTTTGTAGGTCTATTTGTTGGCGTTATTTCAGGTATTTTCAGAGTGTATATGGGTGGGGCGGATGCACAAGTTTATCTTATCTCATCTATATTTATTGGTATAATTGCTGGTTATTTTGGCTTACAAg                                                                       >  1:61553/1‑141 (MQ=255)
                                                               ttGGTGGCCCTTTTGTAGGTCTATTTGTTGGCGTTATTTCAGGTATTTTCAGAGTGTATATGGGTGGGGCGGATGCACAAGTTTATCTTCTCTCATCTATATTTATTGGTATAATTGCTGGTTATTTTGGCTTACAAGCTc                                                                    >  2:10569/1‑141 (MQ=255)
                                                               ttGGTGGCCCTTTTGTAGGTCTATTTGTTGGCGTTATTTCAGGTATTTTCAGAGTGTATATGGGTGGGGCGGATGCACAAGTTTATCTTATCTCATCTATATTTATTGGTATAATTGCTGGTTATTTTGGCTTACAAGCTc                                                                    <  1:21999/141‑1 (MQ=255)
                                                                        cTTTTGTAGGTCTATTTGTAGGCGTTATTTCAGGTATTTTCAGAGTGTATATGGGTGGGGCGGATGCACAAGTCTATCTTATCTCATCTATATTTATTGGTATAATTACTGGTTATTTTGGCTTACAAGCTCAAAGACGCa                                                           >  2:132729/1‑141 (MQ=255)
                                                                              tAGGTCTATTTGTTGGCGTTATCTCAGGTATTTTCAGAGTGTATATGGGAGCGCGGGACGAAAATGCATACCTTATCTCATCTATATTTCTTGGTATAATTGCTGGTTATTTTGGCTTACAAGCTCAAAGACGCAAGCGtt                                                     <  1:224354/141‑1 (MQ=255)
                                                                                     aTTTGTTGGCGTTATTTCAGGTATTTTCAGAGTGTATATGGGTGGGGCGGATACACAAGTTTATCTTATCTCATCTATATTTATTGGTATAATTGCTGGTTATTTTGGCTTACAAGCTCAAAGACGCAAGCGTTACCCGAg                                              >  2:42864/1‑141 (MQ=255)
                                                                                         gttgGCGTTATTTCAGGTGTTTTCAGAGTGTATATGGGTGGGGCGGATGCACAAGTTTATCTTATCTCATCTATATTTATTGGTATAATTGCTGGTTATTTTGGCTTACAAGCTCAAAGACGCAAGCGTTACCCGAGTAtt                                          <  1:188851/141‑1 (MQ=255)
                                                                                         gttgGCGTTATTTCAGGTATTTTCAGAGTGTATATGGGTGGGGCGGATGCACAAGTTTATCTTATCTCATCTATATTTATTGGTATAATTGCTGGTTATTTTGGCTTACAAGCTCAAAGACGCAAGCGTTACCCGAGTAtt                                          <  2:254425/141‑1 (MQ=255)
                                                                                         gttgGCGTTATTTCAGGTATTTTCAGAGTGTATATGGGTGGGGCGGATGCACAAGTTTATCTTATCTCATCTATATTTATTGGTATAATTGCTGGTTATTTTGGCTTACAAGCTCAAAGACGCAAGCGTTACCCGAGTAtt                                          <  2:254414/141‑1 (MQ=255)
                                                                                         gttgGCGTTATTTCAGGTATTTTCAGAGTGTATATGGGTGGGGCGGATGCACAAGTTTATCTTATCTCATCTATATTTATTGGTATAATTGCTGGTTATTTTGGCTTACAAGCTCAAAGACGCAAGCGTTACCCGAGTAtt                                          <  2:1668/141‑1 (MQ=255)
                                                                                            ggCGTTATTTCAGGTATTTTCAGAGTGTATATGGGTGGGGCGGATGCACAAGTTTATCTTATCTCATCTATATTTATTGGTATAATTGCTGGTTATTTTGGCTTACAAGCTCAAAGACGCAAGCGTTACCCGAGTATTGCg                                       >  2:82644/1‑141 (MQ=255)
                                                                                                  aTTTCAGGTATTTTCAGAGTGTATATGGGTGGGGAGGATGCACAAGTTTATCTTATCTCATCTATATTTATTGGTATAATTGCTGGTTATTTTGGCTTACAAGCTCAAAGACGCAAGCGTTACCCGAGTATTGCGAAAAGt                                 <  1:230158/141‑1 (MQ=255)
                                                                                                  aTTTCAGGTATTTTCAGAGTGTATATGGGTGGGGAGGATGCACAAGTTTATCTTATCTCACCTCTATTTTTTGGTATAATTGCTGGTTATTTTGGCTTACAAGCTCAAAGACGCAAGCGTTACCCGAGTATTGCGAAAAGt                                 <  1:185398/141‑1 (MQ=255)
                                                                                                      cAGGTATTTTCAGAGTGTATATGGGTGGGGCGGATGCACAAGTTTATCTTATCTCATCTATATTTATTGGTATAATTGCTGGTTATTTTGGCTTACAAGCTCAAAGACGCAAGCGTTACCCGAGTATTGCGAAAAGTGCCa                             <  2:231800/141‑1 (MQ=255)
                                                                                                        ggTATTTTCAGAGTGTATATGGGTGGGGCGGATGCACAAGTTTATCTTATCTCATCTATATTTATTGGTATAATTGCTGGTTATTTTGGCTTACAAGCTCAAAGACGCAAGCGTTACCCGAGTATTGCGAAAAGTGCCATg                           >  1:111259/1‑141 (MQ=255)
                                                                                                              ttCAGAGTGTATATGGGTGGGGCGGATGCACAAGTTTATCTTATCTCATCTATATTTATTGGTATAATTGCTGGTTATTTTGGCTTACAAGCTCAAAGACGCAAGCGTTACCCGAGTATTGCGAAAAGTGCCATGATTGGa                     >  2:171293/1‑141 (MQ=255)
                                                                                                                 agagTGTATATGGGTGGGGCGGATGCACAAGTTTATCTTATCTCATCTATATTTATTGGTATAATTGCTGGTTATTTTGGCTTACAAGCTCAAAGACGCAAGCGTTACCCGAGTATTGCGAAAAGTGCCATGATTGGAAtt                  <  1:21792/141‑1 (MQ=255)
                                                                                                                                ggggCGGATGCACAAGTTTATCTTATCTCATCTATATTTATTGGTATAATTGCTGGTTATTTTGGCTTACAAGCTCAAAGACGCAAGCGTTACCCGAGTATTGCGAAAAGTGCCATGATTGGAATTGTTATGGAAATGAtt   <  1:182742/141‑1 (MQ=255)
                                                                                                                                 gggCGGATGCACAAGTTTATCTTATCTCATCTATATTTGTTGGTATAATTGCTGGTTATTTTGGCTTACAAGCTCAAAGACGCAAGCGTTACCCGAGTATTGCGAAAAGTGCCATGATTGGAATTGTTATGGAAATGATTc  >  2:256290/1‑141 (MQ=255)
                                                                                                                                    |                                                                                                                                         
GTTTAGATGATGATGTATCTTTAGCTAACACACGTGTATTAACGATAGGTGTCGCAGGATTAGTTGGTGGCCCTTTTGTAGGTCTATTTGTTGGCGTTATTTCAGGTATTTTCAGAGTGTATATGGGTGGGGCGGATGCACAAGTTTATCTTATCTCATCTATATTTATTGGTATAATTGCTGGTTATTTTGGCTTACAAGCTCAAAGACGCAAGCGTTACCCGAGTATTGCGAAAAGTGCCATGATTGGAATTGTTATGGAAATGATTC  >  CP000730/306219‑306488

Alignment Legend
Aligned base mismatch/match (shaded by quality score): ATCG/ATCG < 3 ≤ ATCG/ATCG < 12 ≤ ATCG/ATCG < 22 ≤ ATCG/ATCG < 37 ≤ ATCG/ATCG < 41 ≤ ATCG/ATCG
Unaligned base: atcg    Masked matching base: atcg    Alignment gap:     Deleted base: