Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 889,441 A→C 10.4% intergenic (+178/‑16) PP_0769 → / → PP_0770 sensor histidine kinase/PemI‑like protein

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_002947889,4410AC10.4% 134.1 / 7.1 48intergenic (+178/‑16)PP_0769/PP_0770sensor histidine kinase/PemI‑like protein
Reads supporting (aligned to +/- strand):  ref base A (23/20);  new base C (2/3);  total (25/23)
Fisher's exact test for biased strand distribution p-value = 6.60e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.25e-01

TGGCTTCCCAAGGTGATGTAGGGCATTTCCTTCTTGAGGTTGCAGGGCTTTTGTTCGGTCAGGGCAGATGACAGTGTAGATACGACGTTTCTACATTTCGCTTGCAATCATTGGGAGGGTTTCTACTATGTATCTACCTTGGGAGATACTCATGCAGATCAAGATTCAACAGTGGGGCAACAGCGCCGCGATCCGCTTGCCCGCCGCAGTACTCAAGCAGATGCGCCTCGGTGTCGGCTCCACCCTGAGCCTTGACACAACGGGTGAGACGATGG  >  NC_002947/889306‑889580
                                                                                                                                        |                                                                                                                                           
tGGCTTCCCAAGGTGATGTAGGGCATTTCCTTCTTGAGGTTGCAGGGCTTTTGTTCGGTCAGGGCAGATGACAGTGTAGATACGACGTTTCTACATTTCGCTTGCAATCATTGGGAGGGTTTCTACTATGTATCTACCTTGGGAGATa                                                                                                                                 <  6:21665/148‑1 (MQ=255)
   cTTCCCAAGGTGATGTAGGGCATTTCCTTCTTGAGGTTGCAGGGCTTTTGTTCGGTCAGGGCAGATGACAGTGTAGATACGACGTTTCTACATTTCGCTTGCAATCATTGGGAGGGTTTCTACTATGTATCTCCCAAGGTAGATAcata                                                                                                                             <  5:141572/149‑4 (MQ=255)
   cTTCCCAAGGTGATGTAGGGCATTTCCTTCTTGAGGTTGCAGGGCTTTTGTTCGGTCAGGGCAGATGACAGTGTAGATACGACGTTTCTACATTTCGCTTGCAATCATTGGGAGGGTTTCTACTATGTATCTCCCAAGGTAGATAcata                                                                                                                             <  1:263225/149‑4 (MQ=255)
   cTTCCCAAGGTGATGTAGGGCATTTCCTTCTTGAGGTTGCAGGGCTTTTGTTCGGTCAGGGCAGATGACAGTGTAGATACGACGTTTCTACATTTCGCTTGCAATCATTGGGAGGGTTTCTACTATGTATCTACCTTGGGAGATACTCa                                                                                                                             <  3:278807/149‑1 (MQ=255)
                 gTAGGGCATTTCCTTCTTGAGGTTGCAGGGGCTTTTGTTCGGTCAGGGCAGATGACAGTGTAGATACGACGTTTCTACATTTCGCTTGCAATCATTGGGAGGGTTTCTACTATGTATCTACCTTggg                                                                                                                                      >  1:13917/1‑127 (MQ=255)
                 gTAGGGCATTTCCTTCTTGAGGTTGCAGGGCTTTTGTTCGGTCAGGGCAGATGACAGTGTAGATACGACGTTTCTACATTTCGCTTGCAATCATTGGGAGGGTTTCTACTATGTATCTACCTTggg                                                                                                                                      <  2:13917/126‑1 (MQ=255)
                      gCATTTCCTTCTTGAGGTTGCAGGGCTTTTGTTCGGTCAGGGCAGATGACAGTGTAGATACGACTTTTCTACATTTCGCTTGCAATCATTGGGAGGGTTTCTACTATGTATCTACCTTGGGAGATACTCATGCAGATCAAGATTCaa                                                                                                            >  4:129235/1‑147 (MQ=255)
                             cttgttTAGGTTGCAGGGCTTTTGTTCTGTCAGGGCAGATGACAGTGTAGATACGACGTTTCTACATTTCGCTTGCAATCATTGGGAGGGTTTCTACTATGTATCTACCTTGGGAGATACTCATGCAGATCAAGATTCAACAGTGGGGc                                                                                                   <  8:219468/149‑1 (MQ=255)
                                   gAGGTTGCAGGGCTTTTGTTCGGTCAGGGCAGATGACAGTGTAGATACGACGTTTCTACATTTCGCTTGCAATCATTGGGAGGGTTTCTACTATGTATCTACCTTGGGAGATACTCATGCAGATCAAGATTCAACAGTGGGGCAAc                                                                                                <  3:129235/146‑1 (MQ=255)
                                             aggCTTTTGTTCGGGCAGGGCAGATGACAGTTTAGATACGAAGTTTCTACATTTCGCTTGCAATCATTTGGAGGGTTTCTACTATGTATCTACCTTGGGAGATACTCATGCAGATCAAGATTCAACAGTGTGGCAACAGCGCCGCGATc                                                                                    <  2:284708/148‑1 (MQ=255)
                                                 ttttGTTCGGTCAGGGCAGATGACAGTGTAGATACGACGTTTCTACATTTCGCTTGCAATCATTGGGAGGGTTTCTACTATGTATCTACCTTGGGAGATACTCATGCAGATCAAGATTCAACAGTGGGGCAACAGCGCCGCGATCCGCt                                                                                >  3:174238/1‑149 (MQ=255)
                                                      ttCGGTCAGGGCAGATGACAGTGTAGATACGACGTTTCTACATTTCGCTTGCAATCATTGGGAGGGTTTCTACTATGTATCTACCTTGGGAGATACTCATGCAGATCAAGATTCAACAGTGGGGCAACAGCGCCGCGATCCGCTTGccc                                                                           >  1:195917/1‑149 (MQ=255)
                                                          gTCAGGGCAGATGACAGTGTAGATACGACGTTTCTACATTTCGCTTGCAATCATTGGGAGGGTTTCTACTATGTATCTACCTTGGGAGATACTCATGCAGATCAAGATTCAACAGTGGGGCAACAGCGCCGCGATCCGc                                                                                 <  1:47579/139‑1 (MQ=255)
                                                          gTCAGGGCAGATGACAGTGTAGATACGACGTTTCTACATTTCGCTTGCAATCATTGGGAGGGTTTCTACTATGTATCTACCTTGGGAGATACTCATGCAGATCAAGATTCAACAGTGGGGCAACAGCGCCGCGATCCGc                                                                                 >  2:47579/1‑139 (MQ=255)
                                                               ggCAGATGACAGTGTAGATACGACGTTTCTACATTTCGCTTGCAATCATTGGGAGGGTTTCTACTATGTATCTACCTTGGGAGATACTCATGCAGATCAAGATTCAACAGTGGGGCAACAGcgccgc                                                                                        <  8:43663/127‑1 (MQ=255)
                                                               ggCAGATGACAGTGTAGATACGACGTTTCTACATTTCGCTTGCAATCATTGGGAGGGTTTCTACTATGTATCTACCTTGGGAGATACTCATGCAGATCAAGATTCAACAGTGGGGCAACAGcgccgc                                                                                        >  5:291031/1‑127 (MQ=255)
                                                               ggCAGATGACAGTGTAGATACGACGTTTCTACATTTCGCTTGCAATCATTGGGAGGGTTTCTACTATGTATCTACCTTGGGAGATACTCATGCAGATCAAGATTCAACAGTGGGGCAACAGcgccgc                                                                                        <  6:291031/127‑1 (MQ=255)
                                                               ggCAGATGACAGTGTAGATACGACGTTTCTACATTTCGCTTGCAATCATTGGGAGGGTTTCTACTATGTATCTACCTTGGGAGATACTCATGCAGATCAAGATTCAACAGTGGGGCAACAGcgccgc                                                                                        >  7:43663/1‑127 (MQ=255)
                                                                  aGATGACAGTGTAGATACGACGTTTCTACATTTCGCTTGCAATCATTGGGAGGGTTTCTACTATGTATCTACCTTGGGAGATACTCATGCAGATCAAGATTCAACAGTGGGGCAACAGCGCCGCGATCCGCTTGCCCGCCGCAGTACTc                                                               <  7:115155/149‑1 (MQ=255)
                                                                     tGACAGTGTAGATACGACGTTTCTACATTTCGCTTGCAATCATTGGGAGTGTTTCTACTATGTATCTACCTTGGGAGATACTCATGCAGATCAAGATTCAACAGTGGGGCAACAGCGCCGCGATCCGCTTGCCCGCCGCAGTACTCaa                                                             <  4:174238/148‑1 (MQ=255)
                                                                      gACAGTGTAGATACGACGTTTCTACATTTCCCTTGCAATCATTGGGAGGGTTTCTACTATGTATCTACCTTGGGAGATACTCATGCAGATCAAGATTCAACAGTGGGGCAACAGCGCCGCGATCCGCTTGCCCGCCGCAGTACTCAAg                                                            >  2:180009/1‑148 (MQ=255)
                                                                                  acgacgTTTCTACATTTCGCTTGCAATCATTGGGAGGGTTTCTACTATGTATCTACCTTGGGAGATACTCATGCAGATCAAGATTCAACAGTGGGGCAACAGCGCCGCGATCCGCTTGCCCGCCGCAGTACTCAAGCAGATGCGCCTCg                                               >  1:102898/1‑149 (MQ=255)
                                                                                    gacgTTTCTACATTTCGCTTGCAATCATTGGGAGGGTTTCTACTATGTATCTACCTTGGGAGATACTCATGCAGATCAAGATTCAACAGTGGGGCAACAGCGCCGCGATCCGCTTGCCCGCCGCAGTACTCAAGCAGATGCGCCTCGgt                                             >  1:44389/1‑149 (MQ=255)
                                                                                            tACATTTCGCTTGCAATCATTGGGAGGGTTTCTACTATGTATCTACCTTGGGAGATACTCATGCAGATCAAGATTCAACAGTGGGGCAACAGCGCCGCGATCCGCTTGCCCGCCGCAGTACTCAAGCAGATGCGCCTCGGTGTCGGCTc                                     >  5:196454/1‑149 (MQ=255)
                                                                                             aCATTTCGCTTGCAATCATTGGGAGGGTTTCTACTATGTATCTACCTTGGGAGATACTCATGCTGATCAAGATTCAACAGTGGGGCAACAGCGCCGCGATCCGCTTGCCCGCCGCAGTACTCAAGCAGATGCGCCTCGGTGTcaaggc                                     >  4:185031/1‑143 (MQ=255)
                                                                                             aCATTTCGCTTGCAATCATTGGGAGGGTTTCTACTATGTATCTACCTTGGGAGATACTCATGCAGATCAAGATTCAACAGTGCGGCAACAGCGCCGCGATCCGCTTGCCCGCCGCAGTACTCAAGCAGATGCGCCTCGGTGTCGGCTc                                     >  2:208129/1‑148 (MQ=255)
                                                                                               aTTTCGCTTGCAATCATTGGGAGGGTTTCTACTATGTATCTACCTTGGGAGATACTCATGCAGATCAAGATTCAACAGTGGGGCAACAGCGCCGCGATCCGCTTGCCCGCCGCAGTACTCAAGCAGATGCGCCTCGGTGTCGGCTCCAc                                  >  2:115381/1‑149 (MQ=255)
                                                                                                    gCTTGCAATCATTGGGAGGGTTTCTACTATGTATCTACCTTGGGAGATACTCATGCAGATCAAGATTCAACAGTGGGGCAACAGCGCCGCGATCCGCTTGCCCGCCGCAGTACTCAAGCAGATGCGCCTCGGTGTCGGCTCCACCCTGa                             <  6:196454/149‑1 (MQ=255)
                                                                                                        gCAATCATTGGGAGGGTTTCTACTATGTATCTACCTTGGGAGATACTCATGCAGATCAAGATTCAACAGTGGGGCAACAGCGCCGCGATCCGCTTGCCCGCCGCAGTACTCAAGCAGAAGCGCCTCGGTGTCGGCTCCACCCTGAGc                           <  1:180009/147‑1 (MQ=255)
                                                                                                          aaTCATTGGGAGGGTTTCTACTATGTATCTACCTTGGGAGATACTCATGCAGATCAAGATTCAACAGTGGGGCAACAGCGCCGCGATCCGCTTGCCCGCCGCAGTACTCAAGCAGATGCGCCTCg                                               <  6:19146/125‑1 (MQ=255)
                                                                                                          aaTCATTGGGAGGGTTTCTACTATGTATCTACCTTGGGAGATACTCATGCAGATCAAGATTCAACAGTGGGGCAACAGCGCCGCGATCCGCTTGCCCGCCGCAGTACTCAAGCAGATGCGCCTCg                                               >  5:19146/1‑125 (MQ=255)
                                                                                                          aaTCATTGGGAGGGTTTCTACTATGTATCTACCTTGGGAGATACTCATGCAGATCAAGATTCAACAGTGGGGCAACAGCGCCGCGATCCGCTTGCCCGCCGCAGTACTCAAGCAGATGCGCCTCGGTGTCGGCTCCACCCTGAGCCtt                        >  7:153834/1‑148 (MQ=255)
                                                                                                          aaTCATTGGGAGGGTTTCTACTATGTATCTACCTTGGGAGATACTCATGCAGATCAAGATTCAACAGTGGGGCAACAGCGCCGCGATCCGCTTGCCCGCCGCAGTACTCAAGCAGATGCGCCTCGGTGTCGGCTCCACCCTGAGCCtt                        >  4:264803/1‑148 (MQ=255)
                                                                                                           aTCATTGGGAGGGTTTCTACTATGTATCTACCTTGGGAGATACTCATGCAGATCAAGATTCAACAGTGGGGCAACAGCGCCGCGATCCGCTTGCCCGCCGCAGTACTCAAGCAGATGCGCCTCGGTGTCGGCTCCACCCTGAGCCTTGa                      >  3:237388/1‑149 (MQ=255)
                                                                                                           aTCATTGGGAGGGTTTCTACTATGTATCTACCTTGGGAGATACTCATGCAGATCAAGATTCAACAGTGGGGCAACAGCGCCGCGATCCGCTTGCCCGCCGCAGTACTCAAGCAGATGCGCCTCGGTGTCGGCTCCACCCTGAGCCTTGa                      <  1:115381/149‑1 (MQ=255)
                                                                                                                  ggAGGGTTTCTACTATGTATCTACCTTGGGAGATACTCATGCAGATCAAGATTCAACAGTGGGGCAACAGCGCCGCGATCCGCTTGCCCGCCGCAGTACTCAAGCAGATGCGCCTCGGTGTCGGCTCCACCCTGAGCCTTGACACAACg               <  2:44389/149‑1 (MQ=255)
                                                                                                                   gAGGGTTTCTACTATGTATCTACCTTGGGAGATACTCATGCAGATCAAGATTCAACAGTGGGGCAACAGCGCCGCGATCCGCTTGCCCGCCGCAGTACTCaa                                                             <  2:12520/102‑1 (MQ=255)
                                                                                                                   gAGGGTTTCTACTATGTATCTACCTTGGGAGATACTCATGCAGATCAAGATTCAACAGTGGGGCAACAGCGCCGCGATCCGCTTGCCCGCCGCAGTACTCaa                                                             >  1:12520/1‑102 (MQ=255)
                                                                                                                   gAGGGTTTCTACTATGTATCTACCTTGGGAGATACTCATGCAGATCAAGATTCAACAGTGGGGCAACAGCGCCGCGATCCGCTTGCCCGCCGCAGTACTCAAGCAGATGCGCCTCGGTGTCGGCTCCACCCTGAGCCTTGACACAACgg              >  2:55854/1‑149 (MQ=255)
                                                                                                                    aGGGTTTCTACTATGTATCTACCTTGGGAGATACTCATGCAGATCAAGa                                                                                                                 <  4:175170/49‑1 (MQ=255)
                                                                                                                    aGGGTTTCTACTATGTATCTACCTTGGGAGATACTCATGCAGATCAAGa                                                                                                                 >  3:175170/1‑49 (MQ=255)
                                                                                                                           tgcatgagTATCTCCCAAGGTAGATACTCATGCAGATCAAGATTCAACAGTGGGGCAACAGCGCCGCGATCCGCTTGCCCGCCGCAGTACTCAAGCAGATGCGCCTCGGTGTCGGCTCCACCCTGAGCCTTGACACAACGGGTGAGAc       >  1:293429/8‑148 (MQ=255)
                                                                                                                              ctaTGTATCTACCTTGGGAGATACTCATGCAGATCAAGATTCAACAGTGGGGCAACAGCGCCGCGATCCGCTTGCCCGCCGCAGTACTCAAGCAGATGCGCCTCGGTGTCGGCTCCACCCTGAGCCTTGACACAACGGGTGAGACGATg   >  2:113014/1‑149 (MQ=255)
                                                                                                                               tgagTATCTCCCAAGGTAGATACTCATGCAGATCAAGATTCAACAGTGGGGCAACAGCGCCGCGATCCGCTTGCCCGCCGCAGTACTCAAGCAGATGCGCCTCGGTGTCGGCTCCACCCTGAGCCTTGACACAACGGGTGAGAc       >  7:10298/4‑144 (MQ=255)
                                                                                                                               tgagTATCTCCCAAGGTAGATACTCATGCAGATCAAGATTCAACAGTGGGGCAACAGCGCCGCGATCCGCTTGCCCGCCGCAGTACTCAAGCAGATGCGCCTCGGTGTCGGCTCCACCCTGAGCCTTGACACAACGGGTGAGAc       <  8:10298/141‑1 (MQ=255)
                                                                                                                                 tGTATCTACCTTGGGAGATACTCATGCAGATCAAGATTCAACAGTGGGGCAACAGCGCCGCGATCCGCTTGCCCGCCGCAGTACTCAAGCAGATGCGCCTCGGTGTCGGCTCCACCCTGAGCCTTGACACAACGGGTGAGACGATgg  <  2:102898/147‑1 (MQ=255)
                                                                                                                                    aTCTACCTTGGGAGATACTCATGCAGATCAAGATTCAACAGTGGGGCAACAGCGCCGCGATCCGCTTGCCCGCCGCAGTACTCAAGCAGATGCGCCTCGGTGTCGGCTCCACCCTGAg                            >  7:190987/1‑118 (MQ=255)
                                                                                                                                    aTCTACCTTGGGAGATACTCATGCAGATCAAGATTCAACAGTGGGGCAACAGCGCCGCGATCCGCTTGCCCGCCGCAGTACTCAAGCAGATGCGCCTCGGTGTCGGCTCCACCCTGAg                            <  8:190987/118‑1 (MQ=255)
                                                                                                                                        |                                                                                                                                           
TGGCTTCCCAAGGTGATGTAGGGCATTTCCTTCTTGAGGTTGCAGGGCTTTTGTTCGGTCAGGGCAGATGACAGTGTAGATACGACGTTTCTACATTTCGCTTGCAATCATTGGGAGGGTTTCTACTATGTATCTACCTTGGGAGATACTCATGCAGATCAAGATTCAACAGTGGGGCAACAGCGCCGCGATCCGCTTGCCCGCCGCAGTACTCAAGCAGATGCGCCTCGGTGTCGGCTCCACCCTGAGCCTTGACACAACGGGTGAGACGATGG  >  NC_002947/889306‑889580

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