Predicted mutation
evidence seq id position mutation annotation gene description
RA NC_000913 4,638,120 G→A *475* (TAG→TAA creC → sensory histidine kinase in two‑component regulatory system with CreB or PhoB

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009134,638,1200GA97.8% 143.8 / NA 46*475* (TAG→TAAcreCsensory histidine kinase in two‑component regulatory system with CreB or PhoB
Reads supporting (aligned to +/- strand):  ref base G (0/0);  major base A (19/26);  minor base C (0/1);  total (19/27)
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

CAAAAAAGCAGCGGTCTGGGGTTGGCGTTCGTCAGTGAGGTCGCCCGTTTGTTTAACGGCGAAGTCACGCTGCGCAACGTGCAGGAAGGTGGCGTGCTGGCCTCGCTTCGACTTCACCGTCACTTCACATAGCTTCAAATTCTTCCCACATAGTCTTCGTATCCTGCTGCCATTGCAAAGGAGAAGACTATGTTGAAATCCCCCCTGTTCTGGAAAATGACTAGCCTGTTTGGTGCAGTATTGCTGTTGTTGATTCCGA  >  NC_000913/4637989‑4638247
                                                                                                                                   |                                                                                                                               
cAAAAAAGCAGCGGTCTGGGGTTGGCGTTCGTCAGTGAGGTCGCCCGTTTGTTTAACGGCGAAGTCACGCTGCGCAACGTGCAGGAAGGTGGCGTGCTGGCCTCGCTTCGACTTCACCGTCACTTCACATAACTTCaaa                                                                                                                          <  1:195796/139‑1 (MQ=255)
     aaGCAGCGGTCTGGGGTTGGCGTTCGTCAGTGAGGTCGCCCGTTTGTTTAACGGCGAAGTCACGCTGCGCAACGTGCAGGAAGGTGGCGTGCTGGCCTCGCTTCGACTTCACCGTCACTTCACATAACTTCAAAttctt                                                                                                                     <  1:196817/139‑1 (MQ=255)
          gcGGTCTGGGGTTGGCGTTCGTCAGTGAGGTCGCCCGTTTGTTTAACGGCGAAGTCACGCTGCGCAACGTGCAGGAAGGTGGCGTGCTGGCCTCGCTTCGACTTCACCGTCACTTCACATAACTTCAAATTCTTCCcac                                                                                                                >  1:982975/1‑139 (MQ=255)
           cGGTCTGGGGTTGGCGTTCGTCAGTGAGGTCGCCCGTTTGTTTAACGGCGAAGTCACGCTGCGCAACGTGCAGGAAGGTGGCGTGCTGGCCTCGCTTCGACTTCACCGTCACTTCACATAACTTCAAATTCTTCCcaca                                                                                                               <  1:698436/139‑1 (MQ=255)
           cGGTCTGGGGTTGGCGTTCGTCAGTGAGGTCGCCCGTTTGTTTAACGGCGAAGTCACGCTGCGCAACGTGCAGGAAGGTGGCGTGCTGGCCTCGCTTCGACTTCACCGTCACTTCACATAACTTCAAATTCTTCCcaca                                                                                                               <  2:1165606/139‑1 (MQ=255)
              tCTGGGGTTGGCGTTCGTCAGTGAGGTCGCCCGTTTGTTTAACGGCGAAGTCACGCTGCGCAACGTGCAGGAAGGTGGCGTGCTGGCCTCGCTTCGACTTCACCGTCACTTCACATAACTTCAAATTCTTCCCACATAg                                                                                                            <  1:900309/139‑1 (MQ=255)
                       ggCGTTCGTCAGTGAGGTCGCCCGTTTGTTTAACGGCGAAGTCACGCTGCGCAACGTGCAGGAAGGTGGCGTGCTGGCCTCGCTTCGACTTCACCGTCACTTCACATAACTTCAAATTCTTCCCACATAGTCTTCGTAt                                                                                                   >  1:1239336/1‑139 (MQ=255)
                        gCGTTCGTCAGTGAGGTCGCCCGTTTGTTTAACGGCGAAGTCACGCTGCGCAACGTGCAGGAAGGTGGCGTGCTGGCCTCGCTTCGACTTCACCGTCACTTCACATAACTTCAAATTCTTCCCACATAg                                                                                                            >  1:173204/1‑129 (MQ=255)
                        gCGTTCGTCAGTGAGGTCGCCCGTTTGTTTAACGGCGAAGTCACGCTGCGCAACGTGCAGGAAGGTGGCGTGCTGGCCTCGCTTCGACTTCACCGTCACTTCACATAACTTCAAATTCTTCCCACATAg                                                                                                            <  2:173204/129‑1 (MQ=255)
                        gCGTTCGTCAGTGAGGTCGCCCGTTTGTTTAACGGCGAAGTCACGCTGCGCAACGTGCAGGAAGGTGGCGTGCTGGCCTCGCTTCGACTTCACCGTCACTTCACATAACTTCAAATTCTTCCCACATAGTCTTCGTATc                                                                                                  >  1:1192293/1‑139 (MQ=255)
                               tCAGTGAGGTCGCCCGTTTGTTTAACGGCGAAGTCACGCTGCGCAACGTGCAGGAATGTGGCGTGCTGGCCTCGCTTCGACTTCACCGTCACTTCACATAACTTCAAATTCTTCCCACATAGTCTTCGTATCCtgctgc                                                                                           <  2:791372/139‑1 (MQ=255)
                                    gAGGTCGCCCGTTTGTTTAACGGCGAAGTCACGCTGCGCAACGTGCAGGAAGGTGGCGTGCTGGCCTCGCTTCGACTTCACCGTCACTTCACATAACTTCAAATTCTTCCCACATAGTCTTCGTATCCTGCTGCCATTg                                                                                      >  2:156615/1‑139 (MQ=255)
                                       gTCGCCCGTTTGTTTAACGGCGAAGTCACGCTGCGCAACGTGCAGGAAGGTGGCGTGCTGGCCTCGCTTCGACTTCACCGTCACTTCACATAACTTCAAATTCTTCCCACATAGTCTTCGTATCCTGCTGCCATTGCaa                                                                                   <  2:291121/139‑1 (MQ=255)
                                          gCCCGTTTGTTTAACGGCGAAGTCACGCTGCGCAACGTGCAGGAAGGTGGCGTGCTGGCCTCGCTTCGACTTCACCGTCACTTCACATAACTTCAAATTCTTCCcaca                                                                                                               <  2:882367/108‑1 (MQ=255)
                                          gCCCGTTTGTTTAACGGCGAAGTCACGCTGCGCAACGTGCAGGAAGGTGGCGTGCTGGCCTCGCTTCGACTTCACCGTCACTTCACATAACTTCAAATTCTTCCcaca                                                                                                               >  1:882367/1‑108 (MQ=255)
                                          gCCCGTTTGTTTAACGGCGAAGTCACGCTGCGCAACGTGCAGGAAGGTGGCGTGCTGGCCTCGCTTCGACTTCACCGTCACTTCACATAACTTCAAATTCTTCCCACATAGTCTTCGTATCCTGCTGCCATTGCAAAgg                                                                                <  1:321079/139‑1 (MQ=255)
                                           cccGTTTGTTTAACGGCGAAGTCACGCTGCGCAACGTGCAGGAAGGTGGCGTGCTGGCCTCGCTTCGACTTCACCGTCACTTCACATAACTTCAAATTCTTCCcaca                                                                                                               >  1:462558/1‑107 (MQ=255)
                                           cccGTTTGTTTAACGGCGAAGTCACGCTGCGCAACGTGCAGGAAGGTGGCGTGCTGGCCTCGCTTCGACTTCACCGTCACTTCACATAACTTCAAATTCTTCCcaca                                                                                                               <  2:462558/107‑1 (MQ=255)
                                           cccGTTTGTTTAACGGCGAAGTCACGCTGCGCAACGTGCAGGAAGGTGGCGTGCTGGCCTCGCTTCGACTTCACCGTCACTTCACATAACTTCAAATTCTTCCCACATAGTCTTCGTATCCTGCTGCCATTGCa                                                                                    >  2:206080/1‑134 (MQ=255)
                                           cccGTTTGTTTAACGGCGAAGTCACGCTGCGCAACGTGCAGGAAGGTGGCGTGCTGGCCTCGCTTCGACTTCACCGTCACTTCACATAACTTCAAATTCTTCCCACATAGTCTTCGTATCCTGCTGCCATTGCa                                                                                    <  1:206080/134‑1 (MQ=255)
                                           cccGTTTGTTTAACGGCGAAGTCACGCTGCGCAACGTGCAGGAAGGTGGCGTGCTGGCCTCGCTTCGACTTCACCGTCACTTCACATAACTTCAAATTCTTCCCACATAGTCTTCGTATCCTGCTGCCATTGCAAAGga                                                                               <  2:982975/139‑1 (MQ=255)
                                                ttgtttAACGGCGAAGTCACGCTGCGCAACGTGCAGGAAGGTGGCGTGCTGGCCTCGCTTCGACTTCACCGTCACTTCACATAACTTCAAATTCTTCCCACATAGTCTTCGTATCCTGCTGCCATTGCAAAGGagaaga                                                                          <  2:11431/139‑1 (MQ=255)
                                                             aaGTCACGCTGCGCAACGTGCAGGAAGGTGGCGTGCTGGCCTCGCTTCGACTTCACCGTCACTTCACATAACTTCAAATTCTTCCCACAt                                                                                                              <  1:1185992/90‑1 (MQ=255)
                                                             aaGTCACGCTGCGCAACGTGCAGGAAGGTGGCGTGCTGGCCTCGCTTCGACTTCACCGTCACTTCACATAACTTCAAATTCTTCCCACAt                                                                                                              >  2:1185992/1‑90 (MQ=255)
                                                                     cTGCGCAACGTGCAGGAAGGTGGCGTGCTGGCCTCGCTTCGACTTCACCGTCACTTCACATAACTTCAAATTCTTCCCACATAGTCTTCGTATCCTGCTGCCATTGCa                                                                                    >  2:48400/1‑108 (MQ=255)
                                                                     cTGCGCAACGTGCAGGAAGGTGGCGTGCTGGCCTCGCTTCGACTTCACCGTCACTTCACATAACTTCAAATTCTTCCCACATAGTCTTCGTATCCTGCTGCCATTGCa                                                                                    <  1:48400/108‑1 (MQ=255)
                                                                        cgcAACGTGCAGGAAGGTGGCGTGCTGGCCTCGCTTCGACTTCACCGTCACTTCACATAACTTCAAATTCTTCCCACATAGTCTTCGTATCCTGCTGCCATTGCAAAGga                                                                               <  2:407061/110‑1 (MQ=255)
                                                                        cgcAACGTGCAGGAAGGTGGCGTGCTGGCCTCGCTTCGACTTCACCGTCACTTCACATAACTTCAAATTCTTCCCACATAGTCTTCGTATCCTGCTGCCATTGCAAAGga                                                                               >  1:407061/1‑110 (MQ=255)
                                                                             cGTGCAGGAAGGTGGCGTGCTGGCCTCGCTTCGACTTCACCGTCACTTCACATAACTTCAAATTCTTCCCACATAGTCTTCGTATCCTGCTGCCATTGCAAAGga                                                                               <  1:744035/105‑1 (MQ=255)
                                                                             cGTGCAGGAAGGTGGCGTGCTGGCCTCGCTTCGACTTCACCGTCACTTCACATAACTTCAAATTCTTCCCACATAGTCTTCGTATCCTGCTGCCATTGCAAAGga                                                                               >  2:744035/1‑105 (MQ=255)
                                                                                         tGGCGTGCTGGCCTCGCTTCGACTTCACCGTCACTTCACATAACTTCAAATTCTTCCCACATAGTCTTCGTATCCTGCTGCCATTGCAAAGGAGAAGACTATGTTGAAATCCCCCCTGTTCTGGAAAATGACTAGCCTg                                 >  1:944867/1‑139 (MQ=255)
                                                                                           gCGTGCTGGCCTCGCTTCGACTTCACCGTCACTTCACATAACTTCAAATTCTTCCCACATAGTCTTCGTATCCTGCTGCCATTGCAAAGGAGAAGACTATGTTGAAATCCCCCCTGTTCTGGAAAATGACTAGCCTGtt                               <  2:944867/139‑1 (MQ=255)
                                                                                                cTGGCCTCGCTTCGACTTCACCGTCACTTCACATACCTTCATATTCTTCCCACATAGTGTTCGTGTCCTACTGCCATTGCAAAGGAGAAGACTATGTTGAAATCCCCCCTGTTCTGGAAAATGACTAGCCTGGTTGGTg                          <  2:1134758/139‑1 (MQ=255)
                                                                                                cTGGCCTCGCTTCGACTTCACCGTCACTTCACATAACTTCAAATTCTTCCca                                                                                                                 >  2:1243642/1‑52 (MQ=255)
                                                                                                cTGGCCTCGCTTCGACTTCACCGTCACTTCACATAACTTCAAATTCTTCCca                                                                                                                 <  1:1243642/52‑1 (MQ=255)
                                                                                                cTGGCCTCGCTTCGACTTCACCGTCACTTCACATAACTTCAAATTCTTCCCACATAGTCTTCGTATCCTGCTGCCATTGCAAAGGAGAAGACTATGTTGAAATccc                                                           >  2:1166404/1‑106 (MQ=255)
                                                                                                cTGGCCTCGCTTCGACTTCACCGTCACTTCACATAACTTCAAATTCTTCCCACATAGTCTTCGTATCCTGCTGCCATTGCAAAGGAGAAGACTATGTTGAAATccc                                                           <  1:1166404/106‑1 (MQ=255)
                                                                                                  ggCCTCGCTTCGACTTCACCGTCACTTCACATAACTTCAAATTCTTCCCACATAGTCTTCGTATCCTGCTGCCATTGCAAAGGAGAAGACTATGTTGAAATCCCCCCTGTTCTGGAAAATGACTAGCCTGTTTGGTGCa                        >  1:497189/1‑139 (MQ=255)
                                                                                                        gCTTCGACTTCACCGTCACTTCACATAACTTCAAATTCTTCCCACATAGTCTTCGTATCCTGCTGCCATTGCa                                                                                    <  1:902383/73‑1 (MQ=255)
                                                                                                        gCTTCGACTTCACCGTCACTTCACATAACTTCAAATTCTTCCCACATAGTCTTCGTATCCTGCTGCCATTGCa                                                                                    >  2:902383/1‑73 (MQ=255)
                                                                                                        gCTTCGACTTCACCGTCACTTCACATAACTTCAAATTCTTCCCACATAGTCTTCGTATCCTGCTGCCATTGCAAAGGAGAAGACTATGTTGAAATCCCCCCTGTTCTGGAAAATGACTAGCCTGTTTGGTGCAGTATTg                  <  2:1192293/139‑1 (MQ=255)
                                                                                                                  cACCGTCACTTCACATAACTTCAAATTCTTCCCACATAGTCTTCGTATCCTGCTGCCATTGCAAAGGagaaga                                                                          <  1:366945/73‑1 (MQ=255)
                                                                                                                  cACCGTCACTTCACATAACTTCAAATTCTTCCCACATAGTCTTCGTATCCTGCTGCCATTGCAAAGGagaaga                                                                          >  2:366945/1‑73 (MQ=255)
                                                                                                                        cACTTCACATAACTTCAAATTCTTCCCACATAGTCTTCGTATCCTGCTGCCATTGCAAAGGAGAAGACTATGTTGAAATCCCCCCTGTTCTGGAAAATGACTAGCCTGTTTGGTGCAGTATTGCTGTTGTTGATTCCGa  <  1:156615/139‑1 (MQ=255)
                                                                                                                         aCTTCACATAACTTCAAATTCTTCCCACATAGTCTTCGTATCCTGCTGCCATTGCAAAGGagaa                                                                            >  1:582189/1‑64 (MQ=255)
                                                                                                                         aCTTCACATAACTTCAAATTCTTCCCACATAGTCTTCGTATCCTGCTGCCATTGCAAAGGagaa                                                                            <  2:582189/64‑1 (MQ=255)
                                                                                                                                   |                                                                                                                               
CAAAAAAGCAGCGGTCTGGGGTTGGCGTTCGTCAGTGAGGTCGCCCGTTTGTTTAACGGCGAAGTCACGCTGCGCAACGTGCAGGAAGGTGGCGTGCTGGCCTCGCTTCGACTTCACCGTCACTTCACATAGCTTCAAATTCTTCCCACATAGTCTTCGTATCCTGCTGCCATTGCAAAGGAGAAGACTATGTTGAAATCCCCCCTGTTCTGGAAAATGACTAGCCTGTTTGGTGCAGTATTGCTGTTGTTGATTCCGA  >  NC_000913/4637989‑4638247

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