Predicted mutation
evidence seq id position mutation annotation gene description
RA NC_000913 4,637,791 G→A A366T (GCG→ACG)  creC → sensory histidine kinase CreC

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009134,637,7910GA96.5% 75.9 / ‑4.0 29A366T (GCG→ACG) creCsensory histidine kinase CreC
Reads supporting (aligned to +/- strand):  ref base G (0/1);  new base A (12/16);  total (12/17)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 5.32e-01

GACGCAAAATGCGCGTATGCAGGCATTGGTAGAAACGTTACTACGCCAGGCAAGACTGGAGAATCGTCAGGAAGTCGTTCTGACTGCTGTTGATGTGGCGGCATTATTCCGCCGCGTCAGCGAAGCGCGCACCGTGCAGTTGGCAGAAAAAAAAATCACTCTGCATGTTACGCCCACCGAGGTTAACGTTGCTGCTGAACCGGCGTTACTGGAGCAGGCGCTGGGAAATTTACTGGATAACGCCATCGATTTTACCCCCGAGAGCGGTTGCATCACGCTAAGCGCCGAAGTGGATCAGGAACACGTCACGCTTAAGGTGCTGGATACCGGTAGTGGTATTCCTGACTACGCGCTTTCACGTATTTTTGAACGCTTTTACTCTTTGCCTCGTGCAAATGGGCAAAAAAGCAGCGGTCTGGGGTTGGCGTTCGTCAG  >  NC_000913/4637589‑4638023
                                                                                                                                                                                                          |                                                                                                                                                                                                                                        
gACGCAAAATGCGCGTATGCAGGCATTGGTAGAAACGTTACTACGCCAGGCAAGACTGGAGAATCGTCAGGAAGTCGTTCTGACTGCTGTTGATGTGGCGGCATTATTCCGCCGCGTCAGCGAAGCGCGCACCGTGCAGTTGGCAGAAAAAAAAATCACTCTGCATGTTACGCCCACCGAGGTTAACGTTGCTGCTGAACCGACGTTACTGGAGCAGGCGCTGGGAAATTTACTGGATAACGCCATCGAtt                                                                                                                                                                                          <  1:40968/251‑1 (MQ=255)
                              aGAAACGTTACTCCGCCAGGCAAGACTGGAGAATCGTCAGGAAGTCGTTCTGAATGCTGTTGATGTGGCGGCATTATTCCGCCGCGTCAGCGAAGCGCGCACCGTGCAGTTGGCAGAAAAAAAAATCACTCTGCATGTTACGCCCCCCGAGGTTAACGTTGCTGCTGAAACGACGTTACTGGAGCAGGAGCTGGGAAATTTACTGGATAACGCCATCGATTTTACCCCCGAGAGCGGTTGGATTACGCtc                                                                                                                                                             >  1:324224/1‑249 (MQ=255)
                                                     gACTGGAGAATCGTCAGGAAGTCGTTCTGACTGCTGTTGATGTGGCGGAATTATTCCGCCGCGTCAGCGAAGCGCGCACCGTGCAGTTGGCAGAAAAAAAAATCACTCTGCCTGTTACGCCCACCGAGGTTAACGTTGCTGCTGAACCGACGTTACTGGAGCAGGCGCTGGGAAATTTACTGGATAACGCCATCGATTTTACCCCCGAGAGCGGTTGCATCACGCTAAGCGCCGCAGTGGATCAGGAacac                                                                                                                                     >  1:235839/1‑251 (MQ=255)
                                                                           cGTTCTGACTGCTGTTGATGTGGCGGCATTATTCCGCCGCGTCAGCGAAGCGCGCACCGTGCAGTTGGCAGAAAAAAAAAGCACTCTGCATGGTACGCCCACCGAGGTTAACGTTGCTGCTGAACCGGCGTTACTGGAGCAGGCGCTGTGAGATGTACTGGATAACGCCATCGATTTTACCCCCGAGAGCGGTTGCATCACGc                                                                                                                                                               <  1:75511/203‑1 (MQ=255)
                                                                               cgtcctgctgTTGATGTGGCGGCATTATTCCGCCGCGTCAGCGAAGCGCGCACCGTGCAGGTGGCAGAAAAAAAAATCACTCGGCATGTTACTCCCACCGAGGTTAACGTTGCTGCTGAACCGACGTTACTGGAGCAGGCGCTGGGAAATTTTCTGAATATCGCAATCGATTTTTCCCGCTAGAGCgg                                                                                                                                                                          <  1:401631/184‑1 (MQ=255)
                                                                                                     cATTATTCCGCCGCGTCAGCGAAGCGCGCACCGTGCAGTTGGCAGAAAAAAAAATCACTCTGCATGTTACGCCCACCGAGGTTAACGTTGCTGCTGAACCGACGTTACTGGAg                                                                                                                                                                                                                               <  1:14813/113‑1 (MQ=255)
                                                                                                     cATTATTCCGCCGCGTCAGCGAAGCGCGCACCGTGCAGTTGGCAGAAAAAAAAATCACTCTGCATGTTACGCCCACCGAGGTTAACGTTGCTGCTGAACCGACGTTACTGGAGCAGGCGCTGGGAAATTTACTGGATAACGCCATCGATTTTACCCCCGagagn                                                                                                                                                                            <  1:58792/164‑2 (MQ=255)
                                                                                                                     cAGCGAAGCGCGCACCGTGCAGTTGGCTGAAAAAAAAATCACTCTGCATGTTACGCCCACCGAGGTTAACGTTGCTGCTGAACCGACGTTACTGGAGCAGGCGCTGGGAAATTTACTGGATAACGCCATCGATTTTACCCCCGAGAGCGGTTGCATCACGCTAAGCGCCGAAGTGGATCAGGAACACGTCACGCTTAAGGTGCTGGATACCGGTAGTGGTATTCCTGACTACGCGCTTTCACGTAtttttn                                                                     <  1:132426/251‑2 (MQ=255)
                                                                                                                                 cACCGTGCAGTTGGCAGAAAAAAAAATCACTNTGCATGTTACNCCCNNCGAGNTNANCGTTGCTGCTGAACCGACGTTACTGGAGCNGGCGCTGGGAAATTTACTGGATANNGCNNNCGNNTTTAc                                                                                                                                                                                      <  2:376547/126‑1 (MQ=255)
                                                                                                                                 cACCGTGCAGTTGGCAGAAAAAAAAATCACTCTGCATGTTACGCCCACCGAGGTTAACGTTGCTGCTGAACCGACGTTACTGGAGCAGGCGCTGGGAAATTTACTGGATAACGCCATCGATTTTAc                                                                                                                                                                                      >  1:376547/1‑126 (MQ=255)
                                                                                                                                     gTGCAGTTGGCAGAAAAAAAAATCACTCTGCATGTTACGCCCACCGAGGTTAACGTTGCTGCTGAACCGACGTTACTGGAGCAGGCGCTGGGAAATTTACTGGATAACGCCATCGATTTTACCCCCGAGAGCGGTTGCATCACGCTAAGCGcc                                                                                                                                                       >  1:376574/1‑153 (MQ=255)
                                                                                                                                     gTGCAGTNGNCAGAAAAAAAAATCACTCTGCATGTTACGCCCACCGAGGTTAACGTTGNTGCTGAACCGNCGTNNCTGGNGNANGCGCTGGGAAATTTACTGGATAACGCCATNGATTTTACCCCCGAGAGCGGTTGNNGCNNNCTNNGCGcc                                                                                                                                                       <  2:376574/153‑1 (MQ=255)
                                                                                                                                          gTTGGCAGAAAAAAAAATCACTCTGCATGTTACGCCCACCGAGGTTAACGTTGCTGATGAACCGACGTTACTGGAGCAGGAGCTGGGAAATTTACTGGATAACGCCATCGATTTTACCCCCGAGAGCGGTTGCATCACGCTAAGCGCCGAAGTGGATCAGGAACACGTCACGCTTAAGGTGCTGGATACCGGTAGTGGTATTCCTGACTACGCGCTTTCACGTATTTTTGAACGCTTTTACTCTTTGCCTc                                                >  1:504876/1‑251 (MQ=255)
                                                                                                                                               cAGAAAAAAAAATCACTCTGCATGTTACGCCCACCGAGGTTAACGTTGCTGCTGAACCGACGTTACTGGAGCAGGCGCTGGGAAATTTACTGGATAACGCCATCGATTTTAc                                                                                                                                                                                      >  1:30978/1‑112 (MQ=255)
                                                                                                                                                             aCTCTGCATGTTACGCCCACCGAGGTTAACGTTGCTGCTGAACCGACGTTACCGGAGCATGCGCTGGGAAATTTACTGTATAACGCCATCGATTTTACCCCCGAGATCGGTTGCATCACGCTAAGCGCCGAAGTGGATCAGGAACACGTCACGCTTAAGGTGCTGGa                                                                                                                 <  1:277501/167‑1 (MQ=255)
                                                                                                                                                              ntctGCANNTTACGCCCACCGAGGTTAACGTTGCTGCTGAACCGACGTTACTGGAGCAGGCGCTGGGAAATTTACTGGATAACGCCATCGATTTTAc                                                                                                                                                                                      <  2:30978/96‑1 (MQ=255)
                                                                                                                                                                ctGCATGTTACGCCCGCCGAGGTTAACGTTGCTGGTGATCCGACCTTTCTGGAGCAGGCGCTGGGAAATTTACTGGATAACGCCATCGATTTTACCCCCGAGAGCGGTTGCATCACGCTAAGCGCCGAAGTGGATCAGGAACGCGTCACGCTTAAGGTGCTGGATACCGGTAGTAGTATTCCTGACTACGCACTGTCACGTATTTTTGAACGCTTGTACTCTTTGCCTCGTCCAAATGGGCAAAAAagc                            <  1:239680/249‑1 (MQ=255)
                                                                                                                                                                  gCATGTTACGCCCACCGAGGTTAACGTTGCTGCTGAACCGACGTTACTGGAGCAGGCGCTGGGAAATTTACTGGATAACGCCATCGATTTTACCCCCGAGAGCGGTTGCATCACGCTAAGCGCCGAAGTGGATCAGGAACACGTCACGCTTAAGGTGCTGGATACCGGTAGTGGTATTCCTGACGACGCGCTTTCACGTATTTTTGAACGCTTTTACTCTTTGCCTCGTGCAAATGGGCAAAAAagcagc                         <  1:153854/250‑1 (MQ=255)
                                                                                                                                                                     tGTTACGCCCACCGAGGTTAACGTTGCTGCGGAACCGACGTTACTGGAGCAGGCGCTGGGAAATTTACTGGATAACGCCATCGATTTTACCCCCGAGAGCGGTTGCATCACGCTAAGCTCCGAAGTGGATCAGGAACACGTCACGCTTAAgg                                                                                                                        <  1:397716/152‑1 (MQ=255)
                                                                                                                                                                      gTTACGCCCACCGAGGTTAACGTTGCTGCTGAACCGACGTTACTGGAGCAGGCGCTGGGAAATTTACTGGATAACGCCATCGATTTTACCCCCGAGAGCGGTTGCATCACGCt                                                                                                                                                              <  1:180467/113‑1 (MQ=255)
                                                                                                                                                                               nccGAGGTTAACGTTGCTGCTGAACCGACGTTACTGGAGCAGGCGCTGGGAAATTTACTGGATAACGCCATCGATTTTACCCCCGAGAGCGGTTGCATCACGCTAAGCGCCGAAGTGGATCAGGAACACGTCACGCTTAAGGTGCTGGATACCGGTAGTGGTATTCCTGACTACGCGCTTTCAc                                                                              >  1:297329/2‑184 (MQ=255)
                                                                                                                                                                                    ggTTAACGTTGCTGCTGAACCGACGTTACTGGAGCAGGCGCTGGGAAATTTACTGGATAACGCCATCGATTTTACCCCCGAGAGCGGTTGCATCACGCTAAGCGCCGAAGTGGATCAGGAACACGTCACGCTTAAGGTGCTGGATACCTGTAGTGGTATTCCTGACTACGCGCTTTCACGTATTTTTGATCGCTTTTAcn                                                         <  1:97401/200‑2 (MQ=255)
                                                                                                                                                                                    ggTTAACGTTGCTGCTGAACCGACGTTACTGGAGCAGGCGCTGGGAAATTTACTGGATAACGCCATCGATTTTACCCCCGAGAGCGGTTGCATCACGCTAAGCGCCGAAGTGGATCAGGAACACGTCACGCTTAAGGTGCTGGATACCGGTAGTGGTATTCCTGACTACGCGCTTTCACGTATTTTTGAACGCTTTTACTCTTTGc                                                   >  1:2511/1‑206 (MQ=255)
                                                                                                                                                                                        aaCGTTGCTGCTGAACCGACGTTACTGGAGCAGGCGCTGGGAAATTTACTGGATAACGGCATCGATTTTACCCCCGAGAGCGGTTGCATCACGCTAAGCGCCGAAGTGGATCAGGAACACGTCACGCTTAAGGTGCTGGATACCGGTAGTGGTATTCCTGACTACGCGCTTTCACGTATTTTTGAACGCTTTTACTCTTTGCCTCGTGCAAATGGGCAAAAAAGCAGCGGTCTGGTGTTGGCGTTCGTCAg  <  1:179705/251‑1 (MQ=255)
                                                                                                                                                                                         acggtgctgcTGAACCGACGTTACTGNAGCAGGNNCTNGGAAATTTACTGGATAACGNCANNNANNNTNAACCCNAGNNCGGTTGCATNACGCTAAGCGCCGAAGTgg                                                                                                                                                >  2:114610/5‑108 (MQ=255)
                                                                                                                                                                                         aCGTTGCTGCTGAACCGACGTTACTGGAGCAGGCGCTGGGAAATTTACTGGATAACGCCATCGATTTTACCCCCGAGAGCGGTTGCATCACGCTAAGCGCCGAAGTgg                                                                                                                                                <  1:114610/108‑1 (MQ=255)
                                                                                                                                                                                               ntgctgAACCGACGTTACTGGAGCAGGCGATGGGAAATTTACTGGATAACGCCATCGATTTTACCCCCGAGAGCGGTTGCATCACGCTAAGCGCCGAAGTGGATCAGGAACACGTCACGCTTAAGGTGCTGGATACCGGTAGTGGTATTCCTGACTACGCGCTTTCACGTATTTTTGAACGCTTTTACTCTTTGCCTCGTGCAAATGGGCAAAAAAGCAGCGGTCTGGGGTTGGCGTTc       >  1:404114/2‑239 (MQ=255)
                                                                                                                                                                                               ctgctgAACCGACGTTACTGGAGCAGGCGCTGGGAAATTTACTGGATAACGCCATCGATTTTACCCCCGAGAGCGGTTGCATCACGCTAAgcgc                                                                                                                                                        >  1:21627/1‑94 (MQ=255)
                                                                                                                                                                                               ctgctgAACCGACGTTACTGGAGCAGGCGCTGGGAAATTTACTGGATAACGCCATCGATTTTACCCCCGAGAGCGGTTGCATCACGCTAAGCGCCGAAGTGGATCAGGAACACGTCACGCTTAAGGTGCTGGATACCGGTAGTGGTATTCCTGACTa                                                                                         >  1:171952/1‑157 (MQ=255)
                                                                                                                                                                                               ctgcnnAACCGACGTTACTGGAGCAGGCGCTGGGAAATTTACTGGATAACGCCATCGATTTTACCCCCGAGAGCGGTTGCATCACGCTAAgcgc                                                                                                                                                        <  2:21627/94‑1 (MQ=255)
                                                                                                                                                                                                          |                                                                                                                                                                                                                                        
GACGCAAAATGCGCGTATGCAGGCATTGGTAGAAACGTTACTACGCCAGGCAAGACTGGAGAATCGTCAGGAAGTCGTTCTGACTGCTGTTGATGTGGCGGCATTATTCCGCCGCGTCAGCGAAGCGCGCACCGTGCAGTTGGCAGAAAAAAAAATCACTCTGCATGTTACGCCCACCGAGGTTAACGTTGCTGCTGAACCGGCGTTACTGGAGCAGGCGCTGGGAAATTTACTGGATAACGCCATCGATTTTACCCCCGAGAGCGGTTGCATCACGCTAAGCGCCGAAGTGGATCAGGAACACGTCACGCTTAAGGTGCTGGATACCGGTAGTGGTATTCCTGACTACGCGCTTTCACGTATTTTTGAACGCTTTTACTCTTTGCCTCGTGCAAATGGGCAAAAAAGCAGCGGTCTGGGGTTGGCGTTCGTCAG  >  NC_000913/4637589‑4638023

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: