Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 1,973,431 C→T 7.4% A632T (GCA→ACA)  cheA ← fused chemotactic sensory histidine kinase in two‑component regulatory system with CheB and CheY: sensory histidine kinase/signal sensing protein

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009131,973,4310CT 7.4% 169.5 / 3.6 54A632T (GCA→ACA) cheAfused chemotactic sensory histidine kinase in two‑component regulatory system with CheB and CheY: sensory histidine kinase/signal sensing protein
Reads supporting (aligned to +/- strand):  ref base C (29/21);  new base T (2/2);  total (31/23)
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.29e-01

ATTCCTGGCCTGACGGCTCGCTGGCCAGCTTTGTTACATTCGTCATACCGGTCATATTGTTACCTTTTTACTCATTCAGGCGGCGGTGTTCGCCATACGTTGTTCGCGGTTTATCGCCTGCAAGGCGGAGACATCAACAATCAGTGCCACGCTGCCGTCGCCAAGAATGGTCGCAGCAGAAATGCCGGGGACTTTGCGATAGTTACTTTCAAGGTTTTTAACCACAACCTGGTGTTGACCAATTAATTGATCCACCAGCAAGGCGTAGCGGCGACCGCCACTTTGTAAGATCA  >  NC_000913/1973285‑1973577
                                                                                                                                                  |                                                                                                                                                  
aTTCCTGGCCTGACGGCTCGCTGGCCAGCTTTGTTACATTCGTCATACCGGTCATATTGTTACCTTTTTACTCATTCAGGCGGCGGTGTTCGCCATACGTTGTTCGCGGTTTATCGCCTGCAAGGCGGAGACATCAACAATCAGTGcc                                                                                                                                                   <  7:131218/148‑1 (MQ=255)
 ttCCTGGCCTGACGGCTCGCTGGCCAGCTTTGTTACATTCGTCATACCGGTCATATTGTTACCTTTTTACTCATTCAGGCGGCGGTGTTCGCCATACGTTGTTCGCGGTTTATCGCCTGCAAGGCGGAGACATCAACAATCAGTGCCa                                                                                                                                                  <  6:31817/148‑1 (MQ=255)
          tGACGGCTCGCTGGCCAGCTTTGTTACATTCGTCATACCGGTCATATTGTTACCTTTTTACTCATTCAGGCGGCGGTGTTCGCCATACGTTGTTCGCGGTTTATCGCCTGCAAGGCGGAGACATCAACAATCAGTGCCACGCTGCCGt                                                                                                                                         >  7:249372/1‑148 (MQ=255)
          tGACGGCTCGCTGGCCAGCTTTGTTACATTCGTCATACCGGTCATATTGTTACCTTTTTACTCATTCAGGCGGCGGTGTTCGCCATACGTTGTTCGCGGTTTATCGCCTGCAAGGCGGAGACATCAACAATCAGTGCCACGCTGCCGt                                                                                                                                         >  7:162272/1‑148 (MQ=255)
           gACGGCTCGCTGGCCAGCTTTGTTACATTCGTCATACCGGTCATATTGTTACCTTTTTACTCATTCAGGCGGCGGTGTTCGCCATACGTTGTTCGCGGTTTATCGCCTGCAAGGCGGAGACATCAACAATCAGTGCCACGCTGCCGTc                                                                                                                                        >  4:66908/1‑148 (MQ=255)
                  cGCTGGCCAGCTTTGTTACATTCGTCATACCGGTCATATTGTTACCTTTTTACTCATTCAGGCGGCGGTGTTCGCCATACGTTGTTCGCGGTTTATCGCCTGCAAGGCGGAGACATCAACAATCAGTGCCACGCTGCCGTCGCCAAGa                                                                                                                                 <  8:285452/148‑1 (MQ=255)
                       gCCAGCTTTGTTACATTCGTCATACCGGTCATATTGTTACCTTTTTACTCATTCAGGCGGCGGTGTTCGCCATACGTTGTTCGCGGTTTATCGCCTGCAAGGCGGAGACATCAACAATCAGTGCCACGCTGCCGTCGCCAAGAATgg                                                                                                                             <  8:162272/147‑1 (MQ=255)
                          aGCTTTGTTACATTCGTCATACCGGTCATATTGTTACCTTTTTACTCATTCAGGCGGCGGTGTTCGCCATACGTTGTTCGCGGTTTATCGCCTGCAAGGCGGAGACATCAACAATCAGTGCCACGCTGCCGTCGCCAAGAATgg                                                                                                                             <  6:159849/144‑1 (MQ=255)
                          aGCTTTGTTACATTCGTCATACCGGTCATATTGTTACCTTTTTACTCATTCAGGCGGCGGTGTTCGCCATACGTTGTTCGCGGTTTATCGCCTGCAAGGCGGAGACATCAACAATCAGTGCCACGCTGCCGTCGCCAAGAATgg                                                                                                                             >  5:159849/1‑144 (MQ=255)
                          aGCTTTGTTACATTCGTCATACCGGTCATATTGTTACCTTTTTACTCATTCAGGCGGCGGTGTTCGCCATACGTTGTTCGCGGTTTATCGCCTGCAAGGCGGAGACATCAACAATCAGTGCCACGCTGCCGTCGCCAAGAATGGTCg                                                                                                                          >  3:41736/1‑147 (MQ=255)
                                 ttACATTCGTCATACCGGTCATATTGTTACCTTTTTACTCATTCAGGCGGCGGTGTTCGCCATACGTTGTTCGCGGTTTATCGCCTGCAAGGCGGAGACATCAACAATCAGTGCCACGCTGCCGTCGCCAAGAATGGTCGCAGCAGa                                                                                                                   <  6:29014/147‑1 (MQ=255)
                                      ttCGTCATACCGGTCATATTGTTACCTTTTTACTCATTCAGGCGGCGGTGTTCGCCATACGTTGTTCGCGGTTTATCGCCTGCAAGGCGGAGACATCAACAATCAGTGCCACGCTGCCGTCGCCAAGAATGGTCGCAGCAGAAATGcc                                                                                                             <  2:98296/148‑1 (MQ=255)
                                      ttCGTCATACCGGTCATATTGTTACCTTTTTACTCATTCAGGCGGCGGTGTTCGCCATACGTTGTTCGCGGTTTATCGCCTGCAAGGCGGAGACATCAACAATCAGTGCCACGCTGCCGTCGCCAAGAATGGTCGCAGCAGAAATGcc                                                                                                             >  8:16722/1‑148 (MQ=255)
                                        cGTCATATCGGTCATATTGTTACCTTTTTACTCATTCAGGCGGCGGTGTTCGCCATACGTTGTTCGCGGTTTATCGCCTGCAAGGCGGAGACATCAACAATCAGTGCCACGCTGCCGTCGCCAAGAATGGTCGCAGCAGAAATGCCg                                                                                                            <  7:254209/147‑1 (MQ=255)
                                          tCATACCGGTCATATTGTTACCTTTTTACTCATTCAGGCGGCGGTGTTCGCCATACGTTGTTAGCGGTTTATCGCCTGCAAGGCAGAGCCATCAACAATCAGTGCCACGCTGCCGTCGCCAAGAATGGTCGCAGCAGAAATgaccag                                                                                                          >  1:230490/1‑142 (MQ=255)
                                           cATACCGGTCATATTGTTACCTTTTTACTCATTCAGGCGGCGGTGTTCGCCATACGTTGTTCGCGGTTTATCGCCTGCAAGGCGGAGACATCAACAATCAGTGCCACGCTGCCGTCGCCAAGAATGGTCGCAGCAGAAATGCCGGGGa                                                                                                        >  5:280605/1‑148 (MQ=255)
                                               ccGGTCATATTGTTACCTTTTTACTCATTCAGGCGGCGGTGTTCGCCATACGTTGTTCGCGGTTTATCGCCTGCAAGGCGGAGACATCAACAATCAGTGCCACGCTGCCGTCGCCAAGAATGGTCGCAGCAGAAATGCCGGGGACttt                                                                                                    >  2:312766/1‑148 (MQ=255)
                                                    cATATTGTTACCTTTTTACTCATTCAGGCGGCGGTGTTCGCCATACGTTGTTCGCGGTTTATCGCCTGCAAGGCGGAGACATCAACAATCAGTGCCACGCTGCCGTCGCCAAGAATGGTCGCAGCAGAAATGCCGGGGACTTTGCGa                                                                                                >  7:153745/1‑147 (MQ=255)
                                                        ttGTTACCTTTTTACTCATTCAGGCGGCGGTGTTCGCCATACGTTGTTCGCGGTTTATCGCCTGCAAGGCGGAGACATCAACAATCAGTGCCACGCTGCCGTCGCCAAGAATGGTCGCAGCAGAAATGCCGGGGACTTTGCGATAGtt                                                                                           <  7:305479/148‑1 (MQ=255)
                                                        ttGTTACCTTTTTACTCATTCAGGCGGCGGTGTTCGCCATACGTTGTTCGCGGTTTATCGCCTGCAAGGCGGAGACATCAACAATCAGTGCCACGCTGCCGTCGCCAAGAATGGTCGCAGCAGAAATGCCGGGGACTTTGCGATAGt                                                                                            >  8:305479/1‑147 (MQ=255)
                                                            tACCTTTTTACTCATTCAGGCGGCGGTGTTCGCCATACGTTGTTCGCGGTTTATCGCCTGCAAGGCGGAGACATCAACAATCAGTGCCACGCTGCCGTCGCCAAGAATGGTCGCAGCAGAAATGCCGGGGACTTTGCGATAGTTACtt                                                                                       >  8:152550/1‑148 (MQ=255)
                                                                 ttttACTCATTCAGGCGGCGGTGTTCGCCATACGTTGTTCGCGGTTTATCGCCTGCAAGGCGGAGACATCAACAATCAGTGCCACGCTGCCGTCGCCAAGAATGGTCGCAGCAGAAATGCCGGGGACTTTGCGATAGTTACTTTCAAg                                                                                  >  4:292863/1‑148 (MQ=255)
                                                                 ttttACTCATTCAGGCGGCGGTGTTCGCCATACGTTGTTCGCGGTTTATCGCCTGCAAGGCGGAGACATCAACAATCAGTGCCACGCTGCCGTCGCCAAGAATGGTCGCAGCAGAAATGCCGGGGACTTTGCGATAGTTACTTTCAAg                                                                                  >  7:71116/1‑148 (MQ=255)
                                                                     aCTCATTCAGGGGGCGGTGTTAGCCATACGTTGTTCGCGGTTTATCGCCTGCAAGGCGGAGACATCAACAATCAGTGCCACGCTGCCGTCGCCAAGAATGGTCGCAGCAGAAATGCCGGGGACTTTGCGATAGTTACTTTCAAGGtt                                                                               <  2:230490/147‑1 (MQ=255)
                                                                        cATTCAGGCGGCGGTGTTCGCCATACGTTGTTCGCGGTTTATCGCCTGCAAGGCGGAGACATCAACAATCAGTGCCACGCTGCCGTCGCCAAGAATGGTCGCAGCAGAAATGCCGGGGACTTTGCGATAGTTACTTTCAAGGttttt                                                                            >  7:193098/1‑147 (MQ=255)
                                                                         aTTCAGGCGGCGGTGTTCGCCATACGTTGTTCGCGGTTTATCGCCTGCAAGGCGGAGACATCAACAATCAGTGCCACGCTGCCGTCGCCAAGAATGGTCGCAGCAGAAATGCCGGGGACTTTGCGATAGTTACTTTCAAGGTTTTTa                                                                           <  7:16722/147‑1 (MQ=255)
                                                                         aTTCAGGCGGCGGTGTTCGCCATACGTTGTTCGCGGTTTATCGCCTGCAAGGCGGAGACATCAACAATCAGTGCCACGCTGCCGTCGCCAAGAATGGTCGCAGCAGAAATGCCGGGGACTTTGCGATAGTTACTTTCAAGGTTTTTa                                                                           <  1:312766/147‑1 (MQ=255)
                                                                                cggcggTGTTCGCCATACGTTGTTCGCGGTTTATCGCCTGCAAGGCGGAGACATCAACAATCAGTGCCACGCTGCCGTCGCCAAGAATGGTCGCAGCAGAAATGCCGGGGACTTTGCGATAGTTACTTTCAAGGTTTTTAACCACAAc                                                                   >  6:278364/1‑148 (MQ=255)
                                                                                  gcggTGTTCGCCATACGTTGTTCGCGGTTTATCGCCTGCAAGGCGGAGACATCAACAATCAGTGCCACGCTGCCGTCGCCAAGAATGGTCGCAGCAGAAATGCCGGGGACTTTGCGATAGTTACTTTCAAGGTTTTTAACCACAACCt                                                                 <  2:213363/148‑1 (MQ=255)
                                                                                   cggTGTTCGCCATACGTTGTTCGCGGTTTATCGCCTGCAAGGCGGAGACATCAACAATCAGTGCCACGCTGCCGTCGCCAAGAATGGTCGCAGCAGAAATGCCGGGGACTTTGCGATAGTTACTTTCAAGGTTTTTAACCACAACCTg                                                                <  8:153745/148‑1 (MQ=255)
                                                                                         tCGCCATACGTTGTTCGCGGTTTATCGCCTGCAAGGCGGAGACATCAACAATCAGTGCCACGCTGCCGTCGCCAAGAATGGTCGCAGCAGAAATGCCGGGGACTTTGCGATAGTTACTTTCAAGGTTTTTAACCACAACCTGGTGtt                                                           >  7:300877/1‑147 (MQ=255)
                                                                                             cATACGTTGTTCGCGGTTTATCGCCTGCAAGGCGGAGACATCAACAATCAGTGCCACGCTGCCGTCGCCAAGAATGGTCGCAGCAGAAATGCCGGGGACTTTGCGATAGTTACTTTCAAGGTTTTTAACCACAACCTGGTGTTGACCa                                                      <  3:66908/148‑1 (MQ=255)
                                                                                                aCGTTGTTCGCGGTTTATCGCCTGCAAGGCGGAGACATCAACAATCAGTGTCACGCTGCCGTCGCCAAGAATGGTCGCAGCAGAAATGCCGGGGACTTTGCGATAGTTACTTTCAAGGTTTTTAACCACAACCt                                                                 >  1:126902/1‑134 (MQ=255)
                                                                                                aCGTTGTTCGCGGTTTATCGCCTGCAAGGCGGAGACATCAACAATCAGTGTCACGCTGCCGTCGCCAAGAATGGTCGCAGCAGAAATGCCGGGGACTTTGCGATAGTTACTTTCAAGGTTTTTAACCACAACCt                                                                 >  1:122505/1‑134 (MQ=255)
                                                                                                aCGTTGTTCGCGGTTTATCGCCTGCAAGGCGGAGACATCAACAATCAGTGTCACGCTGCCGTCGCCAAGAATGGTCGCAGCAGAAATGCCGGGGACTTTGCGATAGTTACTTTCAAGGTTTTTAACCACAACCt                                                                 <  2:122505/134‑1 (MQ=255)
                                                                                                aCGTTGTTCGCGGTTTATCGCCTGCAAGGCGGAGACATCAACAATCAGTGTCACGCTGCCGTCGCCAAGAATGGTCGCAGCAGAAATGCCGGGGACTTTGCGATAGTTACTTTCAAGGTTTTTAACCACAACCt                                                                 <  2:126902/134‑1 (MQ=255)
                                                                                                   ttgttCGCGGTTTATCGCCTGCAAGGCGGAGACATCAACAATCAGTGCCACGCTGCCGTCGCCAAGAATGGTCGCAGCAGAAATGCCGGGGACTTTGCGATAGTTACTTTCAAGGTTTTTAACCACAACCTGGTGTTGACCaattaat                                                >  7:15069/1‑148 (MQ=255)
                                                                                                        cgcgGTTTATCGCCTGCAAGGCGGAGACATCAACAATCAGTGCCACGCTGCCGTCGCCAAGAATGGTCGCAGCAGAAATGCCGGGGACTTTGCGATAGTTACTTTCAAGGTTTTTAACCACAACCTGGTGTTGACCAATTAATTGAt                                            <  7:73930/147‑1 (MQ=255)
                                                                                                         gcgGTTTATCGCCTGCAAGGCGGAGACATCAACAATCAGTGCCACGCTGCCGTCGCCAAGAATGGTCGCAGCAGAAATGCCGGGGACTTTGCGATAGTTACTTTCAAGGTTTTTAACCACAACCTGGTGTTGACCAATTAATTGATc                                           <  4:41736/147‑1 (MQ=255)
                                                                                                          cgGTTTATCGCCTGCAAGGCGGAGACATCAACAATCAGTGCCACGCTGCCGTCGCCAAGAATGGTCGCAGCAGAAATGCCGGGGACTTTGCGATAGTTACTTTCAAGGTTTTTAACCACAACCTGGTGTTGACCAATTAATTGATcca                                         >  4:277077/1‑148 (MQ=255)
                                                                                                           ggTTTATCGCCTGCAAGGCGGAGACATCAACAATCAGTGCCACGCTGCCGTCGCCAAGAATGGTCGCAGCAGAAATGCCGGGGACTTTGCGATAGTTACTTTCAAGGTTTTTAACCACAACCTGGTGTTGACCAATTAATTGATccac                                        <  8:300877/148‑1 (MQ=255)
                                                                                                            gTTTATCGCCTGCAAGGCGGAGACATCAACAATCAGTGCCACGCTGCCGTCTCCAAGAATGGTCGCAGCAGAAATTCCGGGGACTTTGCTATAGTTACTTTCAAGGTTTTTAACCACAACCTGGTGTTGACCAATTAATTGATccac                                        >  3:233063/1‑147 (MQ=255)
                                                                                                                  cGCCTGCAAGGCGGAGACATCAACAATCAGTGCCACGCTGCCGTCGCCAAGAATGGTCGCAGCAGAAATGCCGGGGACTTTGCGATAGTTACTTTCAAGGTTTTTAACCACAACCTGGTGTTGACCaattaatt                                               >  1:199086/1‑134 (MQ=255)
                                                                                                                  cGCCTGCAAGGCGGAGACATCAACAATCAGTGCCACGCTGCCGTCGCCAAGAATGGTCGCAGCAGAAATGCCGGGGACTTTGCGATAGTTACTTTCAAGGTTTTTAACCACAACCTGGTGTTGACCaattaatt                                               <  2:199086/134‑1 (MQ=255)
                                                                                                                  cGCCTGCAAGGCGGAGACATCAACAATCAGTGCCACGCTGCCGTCGCCAAGAATGGTCGCAGCAGAAATGCCGGGGACTTTGCGATAGTTACTTTCAAGGTTTTTAACCACAACCTGGTGTTGACCAATTAATTGATCCACCAGCAAg                                 >  5:250460/1‑148 (MQ=255)
                                                                                                                   gCCTGCAAGGCGGAGACATCAACAATCAGTGCCACGCTGCCGTCGCCAAGAATGGTCGCAGCAGAAATGCCGGGGACTTTGCGATAGTTACTTTCAAGGTTTTTAACCACAACCTGGTGTTGACCAATTAATTGATCCACCAGCAAgg                                <  6:280605/148‑1 (MQ=255)
                                                                                                                      tGCAAGGCGGAGACATCAACAATCAGTGCCACGCTGCCGTCGCCAAGAATGGTCGCAGCAGAAATGCCGGGGACTTTGCGATAGTTACTTTCAAGGTTTTTAACCACAACCTGGTGTTGACCAATTAATTGATCCACCAGCAAGGCg                              >  5:50980/1‑147 (MQ=255)
                                                                                                                       gCAAGGCGGAGACATCAACAATCAGTGCCACGCTGCCGTCGCCAAGAATGGTCGCAGCAGAAATGCCGGGGACTTTGCGATAGTTACTTTCAAGGTTTTTAACCACAACCTGGTGTTGACCAATTAATTGATCCACCAGCAAGGCg                              >  8:11993/1‑146 (MQ=255)
                                                                                                                               gagaCATCAACAATCAGTGCCACGCTGCCGTCGCCAAGAATGGTCGCAGCAGAAATGCCGGGGACTTTGCGATAGTTACTTTCAAGGTTTTTAACCACAACCTGGTGTTGACCAATTAATTGATCCACCAGCAAGGCGTAGCGGCGAc                    >  6:9834/1‑148 (MQ=255)
                                                                                                                                      caacaaTCAGTGCCACGCTGCCGTCGCCAAGAATGGTCGCAGCAGAAATGCCGGGGACTTTGCGATAGTTACTTTCAAGGTTTTTAACCACAACCTGGTGTTGACCAATTAATTGATCCACCAGCAAGGCGTAGCGGCGACCGCCACt             >  6:187564/1‑148 (MQ=255)
                                                                                                                                      caacaaTCAGTGCCACGCTGCCGTCGCCAAGAATGGTCGCAGCAGAAATGCCGGGGACTTTGCGATAGTTACTTTCAAGGTTTTTAACCACAACCTGGTGTTGACCAATTAATTGATCCACCAGCAAGGCGTAGCGGCGACCGCCACt             >  6:191059/1‑148 (MQ=255)
                                                                                                                                       aacaaTCAGTGCCACGCTGCCGTCGCCAAGAATGGTCGCAGCAGAAATGCCGGGGACTTTGCGATAGTTACTTTCAAGGTTTTTAACCACAACCTGGTGTTGACCAATTAATTGATCCACCAGCAAGGCGTAGCGGCGACCGCCACtt            >  3:8761/1‑148 (MQ=255)
                                                                                                                                             cAGTGCCACGCTGCCGTCGCCAAGAATGGTCGCAGCAGAAATGCCGGGGACTTTGCGATAGTTACTTTCAAGGTTTTTAACCACAACCTGGTGTTGACCAATTAATTGATCCACCAGCAAGGCGTAGCGGCGACCGCCACTTTGt         >  7:144037/1‑145 (MQ=255)
                                                                                                                                              aGTGCCACGCTGCCGTCGCCAAGAATGGTCGCAGCAGAAATGCCGGGGACTTTGCGATAGTTACTTTCAAGGTTTTTAACCACAACCTGGTGTTGACCAATTAATTGATCCACCAGCAAGGCGTAGCGGCGACCGCCACTTTGTa        <  8:144037/145‑1 (MQ=255)
                                                                                                                                               gTGCCACGCTGCCGTCGCCAAGAATGGTCGCAGCAGAAATGCCGGGGACTTTGCGATAGTTACTTTCAAGGTTTTTAACCACAACCTGGTGTTGACCAATTAATTGATCCACCAGCAAGGCGTAGCGGCGACCGCCACTTTGTAAGAt    <  4:8761/148‑1 (MQ=255)
                                                                                                                                                  ccACGCTGCCTTCGCCAAGAATGGTCGCAGCAGAAATGCCTGTGACTTTGCGATAGTTACTTTCAAGGTTTTTAACCACAACCTGGTGTTGACCAATTAATTGATCCACCAGCAAGGCGTAGCGGCGACCGCCACTTTGTAAGa     <  5:252386/144‑1 (MQ=255)
                                                                                                                                                  ccACGCTGCCGTCGCCAAGAATGGTCGCAGCAGAAATGCCGGGGACTTTGCGATAGTTACTTTCAAGGTTTTTAACCACAACCTGGTGTTGACCAATTAATTGATCCACCAGCAAGGCGTAGCGGCGACCGCCACTTTGTAAGATca  <  5:278364/147‑1 (MQ=255)
                                                                                                                                                  |                                                                                                                                                  
ATTCCTGGCCTGACGGCTCGCTGGCCAGCTTTGTTACATTCGTCATACCGGTCATATTGTTACCTTTTTACTCATTCAGGCGGCGGTGTTCGCCATACGTTGTTCGCGGTTTATCGCCTGCAAGGCGGAGACATCAACAATCAGTGCCACGCTGCCGTCGCCAAGAATGGTCGCAGCAGAAATGCCGGGGACTTTGCGATAGTTACTTTCAAGGTTTTTAACCACAACCTGGTGTTGACCAATTAATTGATCCACCAGCAAGGCGTAGCGGCGACCGCCACTTTGTAAGATCA  >  NC_000913/1973285‑1973577

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