Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 4,477,206 A→G 8.5% P262P (CCT→CCC)  PP_3968 ← sensor histidine kinase

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0029474,477,2060AG 8.5% 127.2 / 2.1 47P262P (CCT→CCC) PP_3968sensor histidine kinase
Reads supporting (aligned to +/- strand):  ref base A (20/23);  new base G (2/2);  total (22/25)
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

CATCGCTGGGAATCGGCTCTCCGGAATTGTGCACGCTGAAGCGCGGAACTCCGGCTTCGCTGAAGAAAGTGACCTGGATCGGTAGCGTCGAGTCACCATGGAGCAGGGCATTGCCAATCAAATTGGAGAAGACCTGAGCCATGCGATCAGGATCAAACTGACCGGTGATCTGTACGGATTCGTTCAAAATAATGTTCGCCCAAGGGGAACCCGCGCGCAGTTCCTCGATGACCGAGCGACAGATCGGGTTCAGCTCCACCTGCTGCAGGTGCACCGGGATACCTGTTCC  >  NC_002947/4477058‑4477346
                                                                                                                                                    |                                                                                                                                            
caTCGCTGGGAATCGGCTCTCCGGAATTGTGCACGCTGAAGCGCGGAACTCCGGCTTCGCTGAAGAAAGTGACCTGGATCGGTAGCGTCGAGTCACCATGGAGCAGGGCATTGCCAATCAAATTGGAGAAGACCTGAGCCATGCGATCa                                                                                                                                              <  1:147327/149‑1 (MQ=255)
       gggAATCGGCTCTCCGGAATTGTGCACGCTGAAGCGCGGAACTCCGGCTTCGCTGAAGAAAGTGACCTGGATCGGTAGCGTCGAGTCACCATGGAGCAGGGCATTGCCAATCAAATTGGAGAAGACCTGAGCCATGCGATCAGGATCa                                                                                                                                        >  4:302138/1‑148 (MQ=255)
        ggAATCGGCTCTCCGGAATTGTGCACGCTGAAGCGCGGAACTCCGGCTTCGCTGAAGAAAGTGACCTGGATCGGTAGCGTCGAGTCACCATGGAGCAGGGCATTGCCAATCAAATTGGAGAAGACCTGAGCCATGCGATCAGGATCaaa                                                                                                                                      <  3:23372/149‑1 (MQ=255)
         gattCGGCTCTCCGGAATTGTCCACGCTGAAGCGCGGAACTCCGGCATCGCTGAAGAAAGTGACCTGGATCGGTAGCGTCGAGTCACCATGGAGCAGGGCATTGCCAATCAAATTGGAGAAGACCTGAGCCATGCGATCAGGATCaaa                                                                                                                                      >  5:56703/4‑148 (MQ=255)
          aaTCGGCTCTCCGGTATTGTGCACGCTGAAGCGCGGAACTCCGGCTTCGCTGAAGAAAGTGACCTGGATCGGTAGCGTCGAGTCACCATGGAGCAGGGCATTGCCAATCAAATTGGAGAAGACCTGAGCCATGCGATCAGGa                                                                                                                                           >  6:148084/1‑142 (MQ=255)
          aaTCGGCTCTCCGGTATTGTGCACGCTGAAGCGCGGAACTCCGGCTTCGCTGAAGAAAGTGACCTGGATCGGTAGCGTCGAGTCACCATGGAGCAGGGCATTGCCAATCAAATTGGAGAAGACCTGAGCCATGCGATCAGGa                                                                                                                                           <  5:148084/142‑1 (MQ=255)
              ggCTCTCCGGAATTGTGCACGCTGAAGCGCGGAACTCCGGCTTCGCTGAAGAAAGTGACCTGGATCGGTAGCGTCGAGTCACCATGGAGCAGGGCATTGCCAATCAAATTGGAGAAGACCTGAGCCATGCGATCAGGATCAAACTGAc                                                                                                                                 <  2:124798/148‑1 (MQ=255)
                   tcCGGAATTGTGCACGCTGAAGCGCGGAACTCCGGCTTCGCTGAAGAAAGTGACCTGGATCGGTAGCGTCGAGTCACCATGGAGCAGGGCATTGCCAATCAAATTGGAGAAGACCTGAGCCATGCGATCAGGATCAAACTGACCGGTGa                                                                                                                           >  1:75446/1‑149 (MQ=255)
                      ggAATTGTGCACGCTGAAGCGCGGAACTCCGGCTTCGCTGAAGAAAGTGACCTGGATCGGTAGCGTCGAGTCACCATGGAGCAGGGCATTGCCAATCAAATTGGAGAAGACCTGAGCCATGCGATCAGGATCAAACTGACCGGTGATCt                                                                                                                        <  6:221884/149‑1 (MQ=255)
                         aTTGTGCACGCTGAAGCGCGGAACTCCGGCTTCGCTGAAGAAAGTGACCTGGATCGGTAGCGTCGAGTCACCATGGAGCAGGGCATTGCCAATCAAATTGGAGAAGACCTGAGCCATGCGATCAGGATCAAACTGACCGGTGATCTGTa                                                                                                                     <  2:289567/149‑1 (MQ=255)
                          ttGTGCACGCTGAAGCGCGGAACTCCGGCTTCGCTGAAGAAAGTGACCTGGATCGGTAGCGTCGAGTCACCATGGAGCAGGGCATTGCCAATCAAATTGGAGAAGACCTGAGCCATGCGATCAGGATCAAACTGACCGGTGATCTGttc                                                                                                                    <  4:311392/149‑3 (MQ=255)
                                aCGCTGAAGCGCGGAACTCCGGCTTCGCTGAAGAAAGTGACCTGGATCGGTAGCGTCGAGTCACCATGGAGCAGGGCATTGCCAATCAAATTGGAGAAGACCTGAGCCATGCGATCAGGATCAAACTGACCGGTGATCTGTACGGATTc                                                                                                              <  2:75446/149‑1 (MQ=255)
                                       aGCGCGGAACTCCGGCTTCGCTGAAGAAAGTGACCTGGATCGGTAGCGTCGAGTCACCATGGAGCAGGGCATTGCCAATCAAATTGGAGAAGACCTGAGCCATGCGATCAGGATCAAACTGACCGGTGATCTGTACGTATTCGTTCaaa                                                                                                       <  4:166201/149‑1 (MQ=255)
                                                   cGGCTTCGCTGAAGAAAGTGACCTGGATCGGTAGCGTCGAGTCACCATGGAGCAGGGCATTGCCAATCAAATTGGAGAAGACCTGAGCCATGCGATCAGGATCAAACTGACCGGTGATCTGTACGGATTCGTTCAAAATAATGTTCGcc                                                                                           >  4:205728/1‑149 (MQ=255)
                                                     gCTTCGCTGAAGAAAGTGACCTGGATCGGTAGCGTCGAGTCACCATGGAGCAGGGCATTGCCAATCAAATTGGAGAAGACCTGAGCCATGCGATCAGGATCAAACTGACCGGTGATCTGTACGGATTCGTTCAAAATAATGTTc                                                                                              <  4:48799/144‑1 (MQ=255)
                                                     gCTTCGCTGAAGAAAGTGACCTGGATCGGTAGCGTCGAGTCACCATGGAGCAGGGCATTGCCAATCAAATTGGAGAAGACCTGAGCCATGCGATCAGGATCAAACTGACCGGTGATCTGTACGGATTCGTTCAAAATAATGTTc                                                                                              >  3:48799/1‑144 (MQ=255)
                                                      ctacGCTGAAGAAAGTGACCTGGATCGGTAGCGTCGAGTCACCATGGAGCAGGGCATTGCAAAACAAATTGGAGAAGCCCTGAGCCATGCGAACAGCATCAAACTGACCGGTGATCTGTACGGATTCGTTAAAAATAATGTTCGcca                                                                                          >  7:151228/4‑146 (MQ=255)
                                                            tGAAGAAAGTGACCTGGATCGGTAGCGTCGAGTCACCATGGAGCAGGGCATTGCCAATCAAATTGGAGAAGACCTGAGCCATGCGATCAGGATCAAACTGACCGGTGATCTGTACGGATTCGTTCAAAATAATGTTc                                                                                              >  2:144095/1‑137 (MQ=255)
                                                            tGAAGAAAGTGACCTGGATCGGTAGCGTCGAGTCACCATGGAGCAGGGCATTGCCAATCAAATTGGAGAAGACCTGAGCCATGCGATCAGGATCAAACTGACCGGTGATCTGTACGGATTCGTTCAAAATAATGTTc                                                                                              <  1:144095/137‑1 (MQ=255)
                                                              aagaaAGTGACCTGGATCGGTAGCGTCGAGTCACCATGGAGCAGGGCATTGCCAATCAAATTGGAGAAGACCTGAGCCATGCGATCAGGATCaa                                                                                                                                       >  3:293720/1‑94 (MQ=255)
                                                              aagaaAGTGACCTGGATCGGTAGCGTCGAGTCACCATGGAGCAGGGCATTGCCAATCAAATTGGAGAAGACCTGAGCCATGCGATCAGGATCAAACTGACCGGTGATCTGTACGGATTCGTTCAAAATAATGTTCGCCCAAGGGGaa                                                                                  <  3:205728/147‑1 (MQ=255)
                                                              aagaaAGTGACCTGGATCGGTAGCGTCGAGTCACCATGGAGCAGGGCATTGCCAATCAAATTGGAGAAGACCTGAGCCATGCGATCAGGATCAAACTGACCGGTGATCTGTACGGATTCGTTCAAAATAATGTTCGCCCAAGGGGAAc                                                                                 <  7:299701/148‑1 (MQ=255)
                                                              aagaaAGTGACCTGGATCGGTAGCGTCGAGTCAACATGGAGCAGGGCATTGCCAATCAAATTGGAGAAGACCTGAGCCATGCGATCAGGATCaa                                                                                                                                       <  4:293720/94‑1 (MQ=255)
                                                                gaaAGTGACCTGGATCGGTAGCGTCGAGTCACCATGGAGCAGGGCATTGCCAATCAAATTGGAGAAGACCTGAGCCATGCGATCAGGATCAAACTGACCGGTGATCTGTACGGATTc                                                                                                              >  3:226912/1‑117 (MQ=255)
                                                                gaaAGTGACCTGGATCGGTAGCGTCGAGTCACCATGGAGCAGGGCATTGCCAATCAAATTGGAGAAGACCTGAGCCATGCGATCAGGATCAAACTGACCGGTGATCTGTACGGATTc                                                                                                              <  4:226912/117‑1 (MQ=255)
                                                                             aTCGGTAGCGTCGAGTCACCATGGAGCAGGGCATTGCCAATCAAATTGGAGAAGACCTGAGCCATGCGATCAGGATCAAACTGACCGGTGATCTGTACGGATTCGTTCAAAATAATGTTCGCCCAAGGGGAACCCGCGCGCAGTTCCt                                                                  <  3:302138/148‑1 (MQ=255)
                                                                              tCGGTAGCGTCGAGTCACCATGGATCAGGGCATTGCCAATCAAATTGGAGAACACCTGAGCCATGCGATCAGCATCAAACTGACCGGTGATCTGTACGGATTCGTTCAAAATAATGTTCGCCCAAGGCGAACCCGCGCGCAGTTccgc                                                                 >  4:11908/1‑146 (MQ=255)
                                                                                gctAGCGACGAGTCACCATGCAGCAGGTCATTGCCAATCAAATTGGAGAAGACCTGAGCCATGCGATCAGCATCAAACTGACCGGTGATCTGTACGGATTCGTTCAAAATAATGTTCGCCCAAGGGGAACCCGCGCGCAGTTCCTCGa                                                               <  8:158399/146‑1 (MQ=255)
                                                                                 gTAGCGTCGAGTCACCATGGAGCAGGGCATTGCCAATCAAATTGGAGAAGACCTGAGCCATGCGATCAGGATCAAACTGACCGGTGATCTGTACGGATTCGTTCAAAATAATGTTCGCCCAAGGGGAACCCGCGCGCAGTTCCTCg                                                                >  1:309375/1‑146 (MQ=255)
                                                                                 gTAGCGTCGAGTCACCATGGAGCAGGGCATTGCCAATCAAATTGGAGAAGACCTGAGCCATGCGATCAGGATCAAACTGACCGGTGATCTGTACGGATTCGTTCAAAATAATGTTCGCCCAAGGGGAACCCGCGCGCAGTTCCTCg                                                                <  2:309375/146‑1 (MQ=255)
                                                                                    gCGTCGAGTCACCATGGAGCAGGGCATTGCCAATCAATTTGGAGAAGACCTGAGCCATGCGATGAGGATCAAACTGACCGGTGATCTGTACGGATTCGTTCAAAATAATGTTCGCCCAAGGGGAACCCGGGCGCAGTTCCTCGATGAc                                                           <  3:11908/148‑1 (MQ=255)
                                                                                        cGAGTCACCATGGAGCAGGGCATTGCCAATCAAATTGGAGAAGACCTGAGCCATGCGATCAGGATCAAACTGACCGGTGATCTGTACGGATTCGTTCAAAATAATGTTCGCCCAAGGGGAACCCGCGCGCAGTTCCTCGATGACCGAGc                                                      >  3:170522/1‑149 (MQ=255)
                                                                                                          ggCATTGCCAATCAAATTGGAGAAGACCTGAGCCATGCGATCAGGATCAAACTGAGCGGTGATCTGTACGGATTCGTTCAAAATAATGTTCGCCCAAGGGGAACCCGCGCGCAGTTCCTCGATGACCGa                                                        <  6:240606/129‑1 (MQ=255)
                                                                                                          ggCATTGCCAATCAAATTGGAGAAGACCTGAGCCATGCGATCAGGATCAAACTGACCGGTGATCTGTACGGATTCGTTCAAAATAATGTTCGCCCAAGGGGAACCCGCGCGCAGTTCCTCGATGACCGa                                                        >  5:240606/1‑129 (MQ=255)
                                                                                                            cATTGCCAATCAAATTGGAGAAGACCTGAGCCATGCGATCAGGATCAAACTGACCGGTGATCTGTACGGATTCGTTCAAAATAATGTTCGCCCAAGGGGAACCCGCGCGCAGTTCCTCGATGACCGAGCGACAGATCGGGTTCAGCTc                                   >  1:330409/1‑148 (MQ=255)
                                                                                                                  cAATCAAATTGGAGAAGACCTGAGCCATGCGATCAGGATCAAACTGACCGGTGATCTGTACGGATTCGTTCAAAATAATGTTCGCCCAAGGGGAACCCGCGCGCAGTTCCTCGATGACCg                                                         >  3:43754/1‑120 (MQ=255)
                                                                                                                  cAATCAAATTGGAGAAGACCTGAGCCATGCGATCAGGATCAAACTGACCGGTGATCTGTACGGATTCGTTCAAAATAATGTTCGCCCAAGGGGAACCCGCGCGCAGTTCCTCGATGACCg                                                         <  4:43754/120‑1 (MQ=255)
                                                                                                                        aaTTGGAGTAGATCTGAGCCATGCGATCACGATCAAACTGACCGGTGATCTGTACGGATTCGTTCAAAATAATGTTCGCGCTATGGGAACCCGCGCGCAGTTCCTCGATGACCGAGCGACAGAGCGGGTTCAGATCCACCTGCTGGAg                       <  8:151228/148‑1 (MQ=255)
                                                                                                                         aTTGGAGAAGACCTGAGCCATGCGATCAGGATCAAACTGACCAGTGATCTGTACGGATTCGTTCAAAATAATGTTCGCCCAAGGGGAACCCGCGCGCAGTTCCTCGATGACCGAGCGACAGATCGGGTTCAGCTCCACCTGCTGCAGGt                     >  1:68846/1‑149 (MQ=255)
                                                                                                                             gataagaCATGAGCCATGCGATCGTGATCAAACTGACCGGTGATCTGTACGGATTCGTTCAAAATAATGTTCGCCCAAGGGGAACCCGCGCGCAGTTCCTCGATGACCGAGCGACAGATCGGGTTCAg                                      <  6:255378/125‑1 (MQ=255)
                                                                                                                             gagaAGACCTGAGCCATGCGATCGGGATCAAACTGACCGGTGATCTGTACGGATTCGTTCAAAATAATGTTCGCCCAAGGGGAACCCGCGCGCAGTTCCTCGATGACCGAGCGACAGATCGGGTTCAg                                      >  5:255378/1‑128 (MQ=255)
                                                                                                                             gagaAGACCTGAGCCATGCGATCGGGATCAAACTGACCGGTGATCTGTACGGATTCGTTCAAAATAATGTTCGCCCAAGGGGAACCCGCGCGCAGTTCCTCGATGACCGAGCGACAGATCGGGTTCAg                                      >  5:260026/1‑128 (MQ=255)
                                                                                                                             gagaAGACCTGAGCCATGCGATCGGGATCAAACTGACCGGTGATCTGTACGGATTCGTTCAAAATAATGTTCGCCCAAGGGGAACCCGCGCGCAGTTCCTCGATGACCGAGCGACAGATCGGGTTCAg                                      <  6:260026/128‑1 (MQ=255)
                                                                                                                             gagaAGACCTGAGCCATGCGATCAGGATCAAACTGACCGGTGATCTGTACGGATTCGTTCAAAATAATGTTCGCCCAAGGGGAACCCGCGCGCAGTTCCTCGATGACCGAGCGACAGATCGGGTTCAGCTCCACCTGCTGCAGGTGCAc                 >  3:186367/1‑149 (MQ=255)
                                                                                                                             gagaAGACCTGAGCCATGCGATCAGGATCAAACTGACCGGTGATCTGTACGGATTCGTTCAAAATAATGTTCGCCCAAGGGGAACCCGCGCGCAGTTCCTCGATGACCGAGCGACAGATCGGGTTCAGCTCCACCTGCAGCAGGTggag                 <  4:170522/149‑4 (MQ=255)
                                                                                                                              agaagaCCTGAGCCATGCGATCAGGATCAAACTGACCGGTGATCTGTACGGATTCGTTCAAAATAATGTTCGCCCAAGGGGAACCCGCGCGCAGTTCCTCGATGACCGAGCGACAGATCGGGTTCAGCTCCACCTGCTGCAGGTGCAcc                >  1:191136/1‑149 (MQ=255)
                                                                                                                                        aGCCATGCGATCAGGATCAAACTGACCGGTGATCTGTACGGATTCGTTCAAAATAATGTTCGCCCAAGGGGAACCCGCGCGCAGTTCCTCGATGACCGAGCGACAGATCGGGTTCAGCTCCACCTGCTGCAGGTGCACCGGGATACCt       <  2:330409/148‑1 (MQ=255)
                                                                                                                                            aTGCGATCAGGATCAAACTGACCGGTGATCTGTACGGATTCGTTCAAAATAATGTTCGCCCAAGGGGAACCCGCGCGCAGTTCCTCGATGACCGACCGACAGATCGGGTTCAGCTCCACCTGCTGCAGGTGCACCGGGATACATGTTcc  >  3:281066/1‑149 (MQ=255)
                                                                                                                                                    |                                                                                                                                            
CATCGCTGGGAATCGGCTCTCCGGAATTGTGCACGCTGAAGCGCGGAACTCCGGCTTCGCTGAAGAAAGTGACCTGGATCGGTAGCGTCGAGTCACCATGGAGCAGGGCATTGCCAATCAAATTGGAGAAGACCTGAGCCATGCGATCAGGATCAAACTGACCGGTGATCTGTACGGATTCGTTCAAAATAATGTTCGCCCAAGGGGAACCCGCGCGCAGTTCCTCGATGACCGAGCGACAGATCGGGTTCAGCTCCACCTGCTGCAGGTGCACCGGGATACCTGTTCC  >  NC_002947/4477058‑4477346

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: ‑