Predicted mutation
evidence seq id position mutation annotation gene description
RA NC_000913 4,636,925 G→C R77P (CGC→CCC)  creC → sensory histidine kinase CreC

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009134,636,9250GC95.0% 56.5 / NA 21R77P (CGC→CCC) creCsensory histidine kinase CreC
Reads supporting (aligned to +/- strand):  ref base G (1/0);  major base C (12/7);  minor base A (1/0);  total (14/7)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

GGGCTATTTTTTACTGGTGGCGGTGGCAGCCTGGTTCGTACTGGCCATTTTTGTCAAAGAAGTTAAACCGGGCGTGCGAAGAGCAACGGAGGGGACGTTGATCGACACCGCAACGTTGCTGGCGGAGCTGGCGCGTCCCGATTTGCTCTCTGGGGACCCAACGCATGGGCAACTGGCGCAGGCGTTTAATCAGCTACAACATCGCCCGTTTCGCGCCAATATCGGTGGCATTAACAAAGTGCGCAATGAATATCATGTCTATATGACCGATGCGCAGGGCAAAGTATTGTTCGATTCGGCAAATAAAGCCGTTGGACAGGATTATTCGCGCTGGAATGACGTCTGGCTAACGTTGCGTGGTCAGTATGGTGCGCGCAGCACGTTGCAAAATCCTGCCGATCCCGAAAGTTCTGTGATGTATGTTGCCGCACC  >  NC_000913/4636722‑4637153
                                                                                                                                                                                                           |                                                                                                                                                                                                                                    
gggCTATTTTTTACTGGTGGCGGTGGCAGCCTGGTTCGTACTGGCCATTTTTGTCAAAGAAGTTAAACCGGGCGTGCGAAGAGCAACGGAGGGGACGTTGATCGACACCGCAACGTTGCTGGCGGAGCTGGCGCGTCCCGATTTGCTCTCTGGGGACCCAACGCATGGGCAACTGGCGCAGGCGTTTAATCAGCTACAACATCCCCCGTTTCGCGCCAATATCGGTGGCATTAACAAAGTGCGCAATGAat                                                                                                                                                                                       >  1:2425/1‑251 (MQ=255)
                                            ccATTTTAGTCACAGAAGTTCAACCGGGCGTGAGCAGAGAAACGGAGGGGACGTTGATCGACAACGCAACGTTGCTGGTGCAGCTGACGCGCCCCGATTTGATCTCTCGGGACCCAACGCACGGGCAACTGGAGACGGCGTTTAATCAGCTACAACATCCCCAGTTTCGCGCAAATATCCGTGGCATTAAACAAGGGCACAATGAATATCACGTCTATATGAACGATGCGCCGGGCACAGTATTGTTCGa                                                                                                                                            >  1:420064/1‑250 (MQ=255)
                                                               tAAACCGGGCGTGCGAAGAGCAACGGAGGGGACGTGGAGCGACACCGCAACGTTGCTGGCGGAGCTGGCGCGTCCCGATTTGCTCTCTGGGGACCCAACGCATGGGGAACTGGCGCAGGCGTTTAATCAGCTACAACATACCCCGGTTCGCGACAATGTCGGCGGCATTAAAAAAGTGCGCAATGAATATCATGGCTATATGACCGAGGCGCCGGGCAAAGTATTGTTCGATTCGGCAAATAAAGCCGTTg                                                                                                                        >  1:221044/1‑251 (MQ=255)
                                                                      ggCGTGCGAAGAGCAACGGAGGGGACGTTGATCGACACCGCCACGTTGCTGGCGGAGCTGGCGCGTCCCGATTTGCTCTCTGGGGACCCAACGCATGGGCAACTGGCGCAGGCGTTTAATCAGCTACCACATCCCCCGTTTCGCGCCAATATCGGTgg                                                                                                                                                                                                              >  1:198740/1‑158 (MQ=255)
                                                                            cGAAGAGCAACGGAGGGGACGTTGATCGACACCGCAACGTTGCTGGCGGAGCTGGCGCGTCCCGATTTGCTCTCTGGGGACCCAACGCATGGGCAACTGGCGCAGGCGTTTAATCAGCTACAACATCCCCCGTTTCGCGCCAATATCGGTGGCa                                                                                                                                                                                                            >  1:391902/1‑154 (MQ=255)
                                                                            cGAAGAGCAACGGAGGGGACGTTGATCGACACCGCAACGTTGCTGGCGGAGCTGGCGCGTCCCGATTTGCTCTCTGGGGACCCAACGCATGGGCAACTGGCGCAGGCGTTTAATCAGCTACAACATCCCCCGTTTCGCGCCAATATCGGTGGCATTAACAAAGTGCGCAATGAATATCATGTCTATATGACCGATGCGCAGGGCAACGTATTGTTCGATTCGGCAAATAAAGCCGTTGGAc                                                                                                                     >  1:177543/1‑241 (MQ=255)
                                                                                    aaCGGAGGGGACGTTGATCGACACCGCAACGTTGCTGGCGGCGCTGGCGCGTCCCGATTTGCTCTCTGGGGACCCAACGCATGGGCAACTGGCGCAGGCGTTTAATCAGCTACAACATCCCCCGTTTCGCGCCAATATCGGTGGCATTAACAAAGTGAGCAATGAATATCATGTCTATATGACCGATGCGCAGGGCAAAGTATTGTTCGATTCGGc                                                                                                                                      >  1:265417/1‑216 (MQ=255)
                                                                                                          acCGCAACGTTGCTGGCGGAGCTGGCGCGTCCCGATTTGCTCTCTGGGGACCCAACGCATGGGCAACTGGCGCAGGCGTTTAATCAGCTACAACATCCCCCGTTTcgc                                                                                                                                                                                                                            <  1:81817/108‑1 (MQ=255)
                                                                                                                                 ggCGCGTCCCGAGTTGCTCTCTGGGGACCCAACGCATGGGCAACTGGCGCAGGCGTTTAATCAGCTACAACATCCCCCGTTTCGCGCCAATATCGGTGGCATTAACAAAGTGCGCAATGAATATCATGTCTATATGACCGATGCGCAGGGCAAAGTATTGTACGATTCGGCAAATAAAGCCGTTGGACAGGAATATTCGCGCTGGAATGACGTCTGGCTAACGTTGCGTGGTCAGTATGGTGCGCgcag                                                        <  1:73573/249‑1 (MQ=255)
                                                                                                                                     ngTCCCGATTTGCTCTCTGGGGACCCAACGCATGGGCAACTGGCGCAGGCGTTTAATCAGCTACAACATCGCCAGTTTCGCGCCAATATCGGTGGCATTAACAAAGTGCGCAATGCATATCATGTCTATATGACCGATGCGCAGGGCAAAGTATTGTTCGATTCGGCAAATAAAGCCGTTGGACAGGATTATTCGCGCTGGAATGACGTCTGGCTAACGTTGCGTGGTCAGTATGGTGCGCGCAGCACGtt                                                  >  1:45138/2‑251 (MQ=255)
                                                                                                                                                   ctctGGGGACCCAACGCATGGGCAACTGGCGCAGGCGTTTAATCAGCTCCAACATCCCCCGTTTCGCGCAAATCTCGGTGGCATTAACAAAGTGCGCAATGAATATCATGTCTATATGAACGATGCGCAGGGCAACGTATTGTTCGATTCGGCAAATAAAGCCGTTGGACAGGATTATTCGCGCTGGAATGACGTCTGGCTAACGTTGCGTGGTCAGTATGGTGCGCGCAGCACGTTGCAAAATCCTGcc                                     >  1:320601/1‑250 (MQ=255)
                                                                                                                                                   ctctGGGGACCCAACGCATGGGCAACTGGCGCAGGCGTTTAATCAGCTACAACATCCCCCGTTTCGCGCCAATATCGGTGGCATTAACAAAGTGCGCAATGAATATCATGTCTATATGACCGATGCGCAGGGCAAAGTATTGTTCGATTCGGCAAATAAAGCCGTTGGACAGGATTATTCGCGCTGGAATGACGTCTGGCTAACGTTGCGTGGTCAGTATGGTGCGCGCAGCACGTTGCAAAATCCtn                                       <  1:62050/248‑2 (MQ=255)
                                                                                                                                                       ggggACCCNANNNANNGGCAACTGGCGCAGGCGTTTAATCAGCTACAACATCCCCCGTTTNGCGCCAATATCGGTGNCANTAACNAAGNGCGCAATGAATATCATGTCTATATGACCGATGCGCAGGGCAAAGTATTGTTCGATTCGGCAAATAACGCCGTTGGACAGGATTATTCGCGCTGGCATGACGTCTGNCTANCGttt                                                                               >  2:356830/1‑203 (MQ=255)
                                                                                                                                                       ggggACCCAACGCATGGGCAACTGGCGCAGGCGTTTAATCAGCTACAACATCCCCCGTTTCGCGCCAATATCGGTGGCATTAACAAATTGCGCAATGAATATCATGTCTATATTACCGATGCGCAGGGCAAAGGATTGTTCTATTCGGCAAATAAAGCCGTTGGACAGGATTATTCGCGCTGGAATGACGTCTGGCTAACGTTg                                                                               <  1:356830/204‑1 (MQ=255)
                                                                                                                                                            nccAACGCATGGGCAACTGGCGCAGGCGTTTAATCAGCTACAACATCCCCCGTTTCGAGccc                                                                                                                                                                                                                        >  1:59046/2‑61 (MQ=255)
                                                                                                                                                            cccAACGCATGGGCAACTGGCGCAGGCGTTTAATCAGCTACAACATCCCCCGTTTCGCGCCAATATCGGTGGCATTAACAAAGTGCGCAATGAATATCATGTCTATATGACCGATGCGCAGGGCAAAGTATTGTTCGATTCGGCAAATAAAGCCGTTGGACAGGATTATTCGCGCTGGAATGACGTCTGGCTAACGTTGCGTGGTCAGTATGGTGCGCGCAGCACGTTGCAAAATCCTGCCGATCCCGaa                            >  1:140303/1‑250 (MQ=255)
                                                                                                                                                                 ngNATNGGCAANTGGCNCAGGCGTTTAATCAGNTACAACATCCCCCGTTTCGNNCNAGTATNGGTgg                                                                                                                                                                                                              <  2:198740/66‑1 (MQ=255)
                                                                                                                                                                   cATGGGCAACTGGCGCAGGCGTTTAATCAGCTACAACATCCCCCGTTTCGCGCCAATATCGGTGGCATTAACAAAGTGCGCAATGAATATCATGTCTATATGACCGATGCGCAGGGCAAAGTATTGTTCGATTCGGCAAATAAAGCCGTTGGACAGGATTATTCGCGCTGGAATGACGTCTGGCTAACGTTGCGTGGTCAGTATGGTGCGCGCAGCACGTTGCAAAATCCTGcc                                     >  1:136347/1‑234 (MQ=255)
                                                                                                                                                                                             tCAGCTACAACATCCCCCGTTTCGCGCCAATATCGGTGGCATTAACAAAGTGCGCAATGAATATCATGTCTATATGACCGATGCGCAGGGCAAAGTATTGTTCGATTCGGCAAATAAAGCCGTTGGACAGGATTATTCGCGCTGTAATGACGTCTGGCTATCGTTGCGTGGTCAGTATGGTGCGCGCAGCACGtt                                                  <  1:189215/195‑1 (MQ=255)
                                                                                                                                                                                             tCAGCTACAACATCCCCCGTTTCGCGCAAATATCGGGGGCATTAACAAAGTGCGCAGTGAATATCATGTCTATATGACCCATGCGCAGTGCAAAGTATTGTTCGATTCGGc                                                                                                                                      <  1:424926/111‑1 (MQ=255)
                                                                                                                                                                                             tCAGCNACAACNTNNGCCGTTTAGCGCCAATATCNGTGGCATTAACAAAGNGCGCNATGAANATNa                                                                                                                                                                                   >  2:189215/1‑66 (MQ=255)
                                                                                                                                                                                                gCTACAACATCACCCGTTTCGCGCCAATATCGGTGGCATTAACAAAGTGCGCAATGAATATCATGTCTATATGACCGATGCGCAGGGCAAAGTATTGTTCGATTCGGCAAATAAAGCCGTTGGACAGGATTATTCGCGCTGGAATGACGTCTGGCTAACGTTGCGTGGTCAGTATGGTGCGCGCAGCACGTTGCAAAATCCTGCCGATCCCGAAAGTTCTGTGATGTATGTTGCCGCAcc  >  1:425910/1‑240 (MQ=255)
                                                                                                                                                                                                           |                                                                                                                                                                                                                                    
GGGCTATTTTTTACTGGTGGCGGTGGCAGCCTGGTTCGTACTGGCCATTTTTGTCAAAGAAGTTAAACCGGGCGTGCGAAGAGCAACGGAGGGGACGTTGATCGACACCGCAACGTTGCTGGCGGAGCTGGCGCGTCCCGATTTGCTCTCTGGGGACCCAACGCATGGGCAACTGGCGCAGGCGTTTAATCAGCTACAACATCGCCCGTTTCGCGCCAATATCGGTGGCATTAACAAAGTGCGCAATGAATATCATGTCTATATGACCGATGCGCAGGGCAAAGTATTGTTCGATTCGGCAAATAAAGCCGTTGGACAGGATTATTCGCGCTGGAATGACGTCTGGCTAACGTTGCGTGGTCAGTATGGTGCGCGCAGCACGTTGCAAAATCCTGCCGATCCCGAAAGTTCTGTGATGTATGTTGCCGCACC  >  NC_000913/4636722‑4637153

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