Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 2,427,982 G→A 8.1% intergenic (‑255/‑94) PP_5503 ← / → PP_2128 CheA signal transduction histidine kinase/CheV‑like chemotaxis protein

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0029472,427,9820GA 8.1% 135.2 / 2.7 49intergenic (‑255/‑94)PP_5503/PP_2128CheA signal transduction histidine kinase/CheV‑like chemotaxis protein
Reads supporting (aligned to +/- strand):  ref base G (26/19);  new base A (2/2);  total (28/21)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.44e-01

AAAGGATATCGCGAACATCAGGTTGCACCCGATTGCTGCTAAGGGTTTCCCTGATCAAGGGCTAGGGTAGGGAACTTTTGTCGTGCCCCTGATTTGCCCCGGCGTGTCGGCAGTAAAGCATGATGGTAAAGTGATATCAATCATCCGCGCGGCATTTTCCCGCAAGGGTCAAGGTGTGCGTGAAGCCGCCGATACAGGGTAGATGATTTCACCTTTTACGACTCAAGCCTGGGGGCTTTCCAATGGCTGGCATTCTCGACACGGTAGACCAACGCACGCAACTGGTGGGTGAGAA  >  NC_002947/2427834‑2428128
                                                                                                                                                    |                                                                                                                                                  
aaaGGATATCGCGAACATCAGGTTGCACCCGATTGCTGCTAAGGGTTTCCCTGATCAAGGGCTAGGGTAGGGAACTTTTGTCGTGCCCCTGATTTGCCCCGGCGTGTCGGCAGTAAAGCATGATGGTAAAGTGATATCAATCATCcgcg                                                                                                                                                    >  4:178388/1‑149 (MQ=255)
          gcgAACATCAGGTTGCACCCGATTGCTGCTAAGGGTTTCCCTGATCAAGGGCTAGGGTAGGGAACTTTTGTCGTGCCCCTGATTTGCCCCGGCGTGTCGGCAGTAAAGCATGATGGTAAAGTGATATCAATCATCCGCACGGCATTTTc                                                                                                                                          >  5:342560/1‑149 (MQ=255)
           cgAACATCAGGTTGCACCCGATTGCTGCTAAGGGTTTCCCTGATCAAGGGCTAGGGTAGGGAACTTTTGTCGTGCCCCTGATTTGCCCCGGCGTGTCGGCAGTAAAGCATGATGGTAAAGTGATATCAATCATCCGCACGGCATTTTcc                                                                                                                                         <  6:342560/149‑1 (MQ=255)
             aaCATCAGGTTGCACCCGATTGCTGCTAAGGGTTTCCCTGATCAAGGGCTAGGGTAGGGAACTTTTGTCGTGCCCCTGATTTGCCCCGGCGTGTCGGCAGTAAAGCATGATGGTAAAGTGATATCAATCATCCGCGCGGCATTTTCCCg                                                                                                                                       <  8:80524/149‑1 (MQ=255)
                aTCAGGTTGCACCCGATTGCTGCTAAGGGTTTCCCTGATCAAGGGCTAGGGTAGGGAACTTTTGTCGTGCCCCTGATTTGCCCCGGCGTGTCGGCAGTAAAGCATGATGGTAAAGTGATATCAATCATCCGCGCGGCATTTTCCCGc                                                                                                                                      >  6:94449/1‑147 (MQ=255)
                      ttGCACCCGATTGCTGCTAAGGGTTTCCCTGATCAAGGGCTAGGGTAGGGAACTTTTGTCGTGCCCCTGATTTGCCCCGGCGTGTCGGCAGTAAAGCATGATGGTAAAGTGATATCAATCATACGCGCGGCATTTTCCCGCAAGGGTc                                                                                                                               >  1:332057/1‑148 (MQ=255)
                        gCACCCGATTGCTGCTAAGGGTTTCCCTGATCAAGGGCTAGGGTAGGGAACTTTTGTCGTGCCCCTGATTTGCCCCGGCGTGTCGGCAGTAAAGCATGATGGTAAAGTGATATCAATCATACGCGCGGCATTTTCCCGCAAGGGTCaa                                                                                                                             <  2:332057/148‑1 (MQ=255)
                         cACCCGATTGCTGCTAATGGTTTCCCTGATCAAGGGCTAGGGTAGGGAACTTTTGTCGTGCCCCTGATTTTCCCCGGCGTGTCGGCAGTAAAGCATGATGGTAAAGTGATATCAATCATCCGCGCGGCATTTTCCCGCAAGGGTCAAg                                                                                                                            <  7:158770/148‑1 (MQ=255)
                          aCCCGATTGCTGCTAAGGGTTTCCCTGATCAAGGGCTAGGGTAGGGAACTTTTGTCGTGCCCCTGATTTGCCAAGGCGTGTCGGCAGTAAAGCATGATGGTAAAGTGATATCAATCATCCGCGCGGCATTTTCACGCAAGGGTCAAGgt                                                                                                                          >  1:305857/1‑149 (MQ=255)
                                ttGCTGCTAAGGGTTTCCCTGATCAAGGGCTAGGGTAGGGAACTTTTGTCGTGCCCCTGATTTGCCCCGGCGTGTCGGCAGTAAAGCATGATGGTAAAGTGATATCAATCATCCGCGCGGCATTTTCCCGCAAGGGTCAAGGTGTGCGt                                                                                                                    <  5:114213/149‑1 (MQ=255)
                                     gctAAGGGTTTCCCTGATCAAGGGCTAGGGTAGGGAACTTTTGTCGTGCCCCTGATTTGCCCCGGCGTGTCGGCAGTAAAGCATGATGGTAAAGTGATATCAATCATCCGCGCGGCATTTTACCGCAAGGGTCAAGGTGTGCGTGAAgc                                                                                                               >  8:311796/1‑149 (MQ=255)
                                                 ccTGATCAAGGGCTAGGGTAGGGAACTTTTGTCGTGCCCCTGATTTGCCCCGGCGTGTCGGCAGTAAAGCATGATGGTAAAGTGATATCAATCATCCGCGCGGCATTTTCCCGCAAGGGTCAAGGTGTGCGTGTAGCCGCCGATACAgg                                                                                                   >  1:319525/1‑149 (MQ=255)
                                                  cTGATCAAGGGCTAGGGTAGGGAACTTTTGTCGTGCCCCTGATTTGCCCCGGCGTGTCGGCAGTAAAGCATGATGGTAAAGTGATATCAATCATCCGCGCGGCATTTTCCCGCaa                                                                                                                                    <  7:210171/115‑1 (MQ=255)
                                                  cTGATCAAGGGCTAGGGTAGGGAACTTTTGTCGTGCCCCTGATTTGCCCCGGCGTGTCGGCAGTAAAGCATGATGGTAAAGTGATATCAATCATCCGCGCGGCATTTTCCCGCaa                                                                                                                                    >  8:210171/1‑115 (MQ=255)
                                                  cTGATCAAGGGCTAGGGTAGGGAACTTTTGTCGTGCCCCTGATTAGCCCCGGCGTGTCGGCAGTAAAGCATGATGGTAAAGTGATATCAATCATCCGCGCGGCATTTTCCCGCAAGGGTCAAGGTGTGCGTGAAGCCGCCGATACAggg                                                                                                  >  7:2602/1‑149 (MQ=255)
                                                  cTGATCAAGGGCTAGGGTAGGGAAATTTTGTCGTGCCAATGATTTGACACGGCGTGTCGGCAGTAAAGAATTATTTTAAAGTGATATCAATCATCCGCGCGGCATTTTCCCGCAAGGGTCAAGGTGTGCGTGAAGCGGCCGATACAggg                                                                                                  >  6:298722/1‑149 (MQ=255)
                                                   tGATCAAGGGCTAGGGTAGGGAACTTTTGTCGTGCCCCTGATTTGCCCCGGCGTGTCGGCAGTAAAGCATGATGGTAAAGTGATATCAATCATCCGCGCGGCATTTTCCCGCAAGGGTCAAGGTGTGCGTGAAGCCGCCGATACAggg                                                                                                  >  7:86115/1‑148 (MQ=255)
                                                     aTCAAGGGCTAGGGTAGGGAACTTTTGTCGTGCCCCTGATTTGCCCCGGCGTGTCGGCAGTAAAGCATGATGGTAAAGTGATATCAATCATCCGCGCGGCATTTTCCCGCAAGGGTCAAGGTGTGCGTGAAGCCGCCGATACAGGGt                                                                                                 <  8:86115/147‑1 (MQ=255)
                                                      tCAAGGGCTAGGGTAGGGAACTTTTGTCGTGCCCCTGATTTGCCCCGGCGTGTCGGCAGTAAAGCATGATGGTAAAGTGATATCAATCATCCGCGCGGCATTTTCCCGCAAGGGTCAAGGTGTGCGTGAAGCCGCCGATACAGGGTAga                                                                                              <  5:94449/149‑1 (MQ=255)
                                                        aaGGGCTAGGGTAGGGAACTTTTGTCGTGCCCCTGATTTGCCCCGGCGTGTCGGCAGTAAAGCATGATGGTAAAGTGATATCAATCATCCGCGCGGCATTTTCCCGCAAGGGTCAAGGTGTGCGTGAAGCCGCCGATACAGGGTAgat                                                                                             >  6:169643/1‑148 (MQ=255)
                                                                           ttttGTCGTGCCCCTGATTTGCCCCGGCGTGTCGGCAGTAAAGCATGATGGTAAAGTGATATCAATCATCCGCGCGGCATTTTCCCGCAAGGGTCAAGGTGTGCGTGAAGCCGCCGATACAGGGTAGATGATTTCACCTTTTactact                                                                          >  5:354682/1‑148 (MQ=255)
                                                                           ttttGTCGTGCCCCTGATTTGCCCCGGCGTGTCGGCAGTAAAGCATGATGGTAAAGTGATATCAATCATCCGCGCGGCATTTTCCCGCAAGGGTCAAGGTGTGCGTGAAGCCGCCGATACAGGGTAGATGATTTCACCTTTTACGACt                                                                          >  6:69439/1‑148 (MQ=255)
                                                                                  gTGCCCCTGATTTGCCCCGGCGTGTCGGCAGTAAAGCATGATGGTAAAGTGATATCAATCATCCGCGCGGCATTTTCCCGCAAGGGTCAAGGTGTGCGTGAAGCCGCCGATACAGGGTAGATGATTTCACCTTTTACGACTCAAGCCTg                                                                  <  5:169643/149‑1 (MQ=255)
                                                                                       ccTGATTTGCCCCGGCGTGTCGGCAGTAAAGCATGATGGTAAAGTGATATCAATCATCCGCGCGGCATTTTCCCGCAAGGGTCAAGGTGTGCGTGAAGCCGCCGATACAGGGTAGATGATTTCACCTTTTACGACTCAAGCCTGGGGGc                                                             >  3:96513/1‑149 (MQ=255)
                                                                                                ccccGGCGTGTCGGCAGTAAAGCATGATGGTAAAGTGATATCAATCATCCGCGCGGCATTTTCCCGCAAGGGTCAAGGTGTGCGTGAAGCCGCCGATACAGGGTAGATGATTTCACCTTTTACGACTCAAGCCTGGGGGCTTTCCAAtg                                                    <  4:96513/149‑1 (MQ=255)
                                                                                                ccccGGCGTGTCGGCAGTAAAGCATGATGGTAAAGTGATATCAATCATCCGCACGGCATTTTCCCGCAAGGGTCAAGGTGTGCGTGAAGCCGCCGATACAGGGTAGATGATTTCAcc                                                                                    >  4:372641/1‑117 (MQ=255)
                                                                                                ccccGGCGTGTCGGCAGTAAAGCATGATGGTAAAGTGATATCAATCATCCGCACGGCATTTTCCCGCAAGGGTCAAGGTGTGCGTGAAGCCGCCGATACAGGGTAGATGATTTCAcc                                                                                    <  3:372641/117‑1 (MQ=255)
                                                                                                 cccGGCGTGTCGGCAGTAAAGCATGATGGTAAATTGATATCAATCATCCGCGCGGCATTTTCCCGCAAGGGTCAAGGTGTGCGTGAAGCCGCCGATACAGGGTAGATGAtt                                                                                         <  7:240342/111‑1 (MQ=255)
                                                                                                 cccGGCGTGTCGGCAGTAAAGCATGATGGTAAAGTGATATCAATCATCCGCGCGGCATTTTCCCGCAAGGGTCAAGGTGTGCGTGAAGCCGCCGATACAGGGTAGATGAtt                                                                                         >  8:240342/1‑111 (MQ=255)
                                                                                                  ccGGCGTGTCGGCAGTAAAGCATGATGGTAAAGTGATATCAATCATCCGCGCGGCATTTTCCCGCAAGGGTCAAGGTGTGCGTGAAGCCGCCGATACAGGGTAGATGATTTCACCTTTTACGACTCAAGCCTGGGGGCTTTCCAAtggc                                                  >  8:45002/1‑149 (MQ=255)
                                                                                                    ggCGTGTCGGCAGTAAAGCATGATGGTAAAGTGATATCAATCATCCGCGCGGCATTTTCCCGCAAGGGTCAAGGTGTGCGTGAAGCCGCCGATACAGGGTAGATGATTTCACCTTTTACGACTCAAGCCTGGGGGCTTTCCAAtggctg                                                <  2:319525/149‑1 (MQ=255)
                                                                                                        tgtCGGCAGTAAAGCATGATGGTAAAGTGATATCAATCATCCGCGCGGCATTTTCCCGCAAGGGTCAAGGTGTGCGTGAAGCCGCCGATACAGGGTAGATGATTTCACCTTTTCCGACTCAAGCCTGGGGGCTTTCCAATGGCTGGCAt                                            >  8:69242/1‑149 (MQ=255)
                                                                                                            ggCAGTAAAGCATGATGGTAAAGTGATATCAATCATCCGCGCGGCATTTTCCCGCAAGGGTCAAGGTGTGCGTGAAGCCGCCGATACAGGGTAGATGATTTCACCTTTTACGACTCAAGCCTGGGGGCTTTCCAAtggctggc                                              <  1:194276/143‑1 (MQ=255)
                                                                                                            ggCAGTAAAGCATGATGGTAAAGTGATATCAATCATCCGCGCGGCATTTTCCCGCAAGGGTCAAGGTGTGCGTGAAGCCGCCGATACAGGGTAGATGATTTCACCTTTTACGACTCAAGCCTGGGGGCTTTCCAAtggctggc                                              >  2:194276/1‑143 (MQ=255)
                                                                                                            ggCAGTAAAGCATGATGGTAAAGTGATATCAATCATCCGCGCGGCATTTTCCCGCAAGGGTCAAGGTGTGCGTGAAGCCGCCGATACAGGGTAGATGATTTCACCTTTTACGACTCAAGCCTGGGGGCTTTCCAATGGCTGGCATtctc                                        <  8:2602/149‑1 (MQ=255)
                                                                                                             gCAGTAAAGCATGATGGTAAAGTGATATCAATCATCCGCGCGGCATTTTCCCGCAAGGGTCAAGGTGTGCGTGAAGCCGCCGATACAGGGTAGATGATTTCACCTTTTACGACTCAAGCCTGGGGGCTTTCCAATGGCTGGCATtctc                                        >  8:154467/1‑148 (MQ=255)
                                                                                                               aGTAAAGCATGATGGTAAAGTGATATCAATCATCCGCGCGGCATTTTCCCGCAAGGGTCAAGGTGTGCGTGAAGCCGCCGATACAGGGTAGATGATTTCACCTTTTACGACTCAAGCCTGGGGGCTTTCCAATGGCTGGCATTCTCGa                                      >  2:318183/1‑148 (MQ=255)
                                                                                                                  aaaGCATGATGGTAAAGTGATATCAATCATCCGCGCGGCATTTTCCCGCAAGGGTCAAGGTGTGCGTGAAGCCGCCGATACAGGGTAGATGATTTGACCTTTTACGACTCAAGCCTGGGGGCTTTCCAATGGCTGGCATTCTCGACACg                                  >  4:137215/1‑149 (MQ=255)
                                                                                                                   aaGCATGATGGTAAAGTGATATCAATCATCCGCGCGGCATTTTCCCGCAAGGGTCAAGGTGTGCGTGAAGCCGCCGATACAGGGTAGATGATTTCACCTTTTACGACTCAAGCCTGGGGGCTTTCCAATGGCTGGCATTCTCGACACgg                                 <  5:69439/149‑1 (MQ=255)
                                                                                                                     gCATGATGGTAAAGTGATATCAATCATCCGCGCGGCATTTTCCCGCAAGGGTCAAGGTGTGCGTGAAGCCGCCGATACAGGGTAGATGATTTCACCTTTTACGACTCAAGCCTGGGGGCTTTCCAATGGCTGGCATTCTCGACACGGTa                               >  8:22584/1‑149 (MQ=255)
                                                                                                                      cATGATGGTAAAGTGATATCAATCATCCGCGCGGCATTTTCCCGCAAGGGTCAAGGTGTGCGTGaa                                                                                                                 >  2:318562/1‑66 (MQ=255)
                                                                                                                      cATGATGGTAAAGTGATATCAATCATCCGCGCGGCATTTTCCCGCAAGGGTCAAGGTGTGCGTGaa                                                                                                                 <  1:318562/66‑1 (MQ=255)
                                                                                                                        tgatgGTAAAGTGTTATCATTCATCCGCGCGGCATTTTCCCGCAAGGGTCAAGGTGTGCGTGATGCCGCCGATACAGGGTAGATGATTTCTCCTTTTACGACTCAAGCCTGGGGGCTTTCCAATGGCTGGCATTCTCGACACGGTAGAc                            <  2:305857/149‑1 (MQ=255)
                                                                                                                                 aGTGATATCAATCATCCGCGCGGCATTTTCCCGCAAGGGTCAAGGTGTGCGTGAAGCCGCCGATACAGGGTAGATGATTTCACCTTTTACGACTCAAGCCTGGGGGCTTTCCAATGGCTGGCATTCTCGACACGGTAGACCAacgcacg                   >  5:359687/1‑149 (MQ=255)
                                                                                                                                  gTGATATCAATCATCCGCGCGGCATTTTCCCGCAAGGGTCAAGGTGTGCGTGAAGCCGCCGATACAGGGTAGATGATTTCACCTTTTACGACTCAAGCCTGGGGGCTTTCCAATGGCTGGCATTCTCGACACGGTAGACCAacgcacg                   >  6:286704/1‑148 (MQ=255)
                                                                                                                                       atcaatcaTCCGCGCGACATTTTCACGCAAGGGTCAAGGTGTGCGTGAAGCCGCCGATACAGGGTAGATGATTTCACCTTTTACGACTCAAGCCTGGGGGCTTTCCAATGGCTGGCATTCTCGACACGGTAGACCAACGCACGCAACtg             >  2:252765/1‑149 (MQ=255)
                                                                                                                                          aatcaTCCGCGCGGCATTTTCCCGCAAGGGTCAAGGTGTGCGTGAAGCCGCCGATACAGGGTAGATGATTTCACCTTTTACGACTCAAGCCTGGGGGCTTTCCAATGGCTGGCATTCTCGACACGGTAGACCAACGCACGCAACtggtg          >  8:324706/1‑149 (MQ=255)
                                                                                                                                          aatcaTCCGCGCGGCATTTTCCCGCAAGGGTCAAGGTGTGCGTGAAGCCGCCGATACAGGGTAGATGATTTCACCTTTTACGACTCAAGCCTGGGGGCTTTCCAATGGCTGGCATTCTCGACACGGTAGACCAACGCACGCAACtggtg          <  7:22584/149‑1 (MQ=255)
                                                                                                                                           atcatcCGCGCGGCATTTTCCCGCAAGGGTCAAGGTGTGCGTGAAGCCGCCGATACAGGGTAGATGATTTCACCTTTTACGACTCAAGCCTGGGGGCTTTCCAATGGCTGGCATTCTCGACACGGTAGACCAACGCACGCAACtggtgg         <  7:311796/149‑1 (MQ=255)
                                                                                                                                              atcCGCGCGGCATTTTCCCGCAAGGGTCAAGGTGTGCGTGAAGCCGCCGATACAGGGTAGATGATTTCACCTTTTACGACTCAAGCCTGGGGGCTTTCCAATGGCTGGCATTCTCGACACGGTAGACCAACGCACGCAACTggtgggt       <  6:359687/148‑1 (MQ=255)
                                                                                                                                                  gcgcgGCATTTTCCCGCAAGGGTCAAGGTGTGCGTGAAGCCGCCGATACAGGGTAGATGATTTCACCTTTTACGACTCAAGCCTGGGGGCTTTCCAATGGCTGGCATTCTCGACACGGTAGACCAACGCACGCAACTGGTGGGTGAGaa  <  3:178388/149‑1 (MQ=255)
                                                                                                                                                    |                                                                                                                                                  
AAAGGATATCGCGAACATCAGGTTGCACCCGATTGCTGCTAAGGGTTTCCCTGATCAAGGGCTAGGGTAGGGAACTTTTGTCGTGCCCCTGATTTGCCCCGGCGTGTCGGCAGTAAAGCATGATGGTAAAGTGATATCAATCATCCGCGCGGCATTTTCCCGCAAGGGTCAAGGTGTGCGTGAAGCCGCCGATACAGGGTAGATGATTTCACCTTTTACGACTCAAGCCTGGGGGCTTTCCAATGGCTGGCATTCTCGACACGGTAGACCAACGCACGCAACTGGTGGGTGAGAA  >  NC_002947/2427834‑2428128

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: