Predicted mutation
evidence seq id position mutation annotation gene description
RA NC_000913 2,585,779 (T)7→8 coding (49/1701 nt) narQ → sensory histidine kinase NarQ

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009132,585,7721.T86.4% 133.3 / 13.1 44coding (42/1701 nt)narQsensory histidine kinase NarQ
Reads supporting (aligned to +/- strand):  ref base . (3/3);  new base T (20/18);  total (23/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 = 6.32e-01

TTCGCTATAAACTGTGGCAGATCAAATAATCCCATCACTGACTAAATTGCGTTTCAGCGAACTGGAACATTAATGATTTTTTGTGGAGAAGACGCGTGTGATTGTTAAACGACCCGTCTCGGCCAGTCTGGCCCGGGCCTTTTTTTACATTGTGCTGCTGTCGATTCTTTCCACGGGTATCGCTCTGCTAACTCTGGCGAGCAGTTTGCGCGACGCTGAGGCTATCAATATTGCCGGATCGCTGCGTATGCAGAGTTACCGCCTG  >  NC_000913/2585634‑2585898
                                                                                                                                            |                                                                                                                               
ttCGCTATAAACTGTGGCAGATCAAATAATCCCATCACTGACTAAATTGCGTTTCAGCGAACTGGAACATTAATGATTTTTTGTGGAGAAGACGCGTGTGATTGTTAAACGACCCGTCTCGGCCAGTCTGGCCCGGGCCt                                                                                                                                <  1:490076/140‑1 (MQ=255)
ttCGCTATAAACTGTGGCAGATCAAATAATCCCATCACTGACTAAATTGCGTTTCAGCGAACTGGAACATTAATGATTTTTTGTGGAGAAGACGCGTGTGATTGTTAAACGACCCGTCTCGGCCAGTCTGGCCCGGGCCt                                                                                                                                <  1:377923/140‑1 (MQ=255)
   gCTATAAACTGTGGCAGATCAAATAATCCCATCACTGACTAAATTGCGTTTCAGCGAACTGGAACATTAATGATTTTTTGTGGAGAAGACGCGTGTGATTGTTAAACGACCCGTCTCGGCCAGTCTGGCCCGGGCCtttt                                                                                                                            >  2:52006/1‑140 (MQ=255)
        aaaCTGTGGTAGATCAAATAATCCCATCACTGACTAAATTGCGTTTCAGCGAACTGGAACATTAATGATTTTTTGTGGAGAAGACGCGTGTGATTGTTAAACGACCCGTCTCGGCCAGTCTGGCCCGGGCCtttttttta                                                                                                                       <  2:38274/140‑3 (MQ=255)
              tgGCAGATCAAATAATCCCATCACTGACTAAATTGCGTTTCAGCGAACTGGAACATTAATGATTTTTTGTGGAGAAGACGCGTGTGATTGTTAAACGACCCGTCTCGGCCAGTCTGGCCCGGGCCTTTTTTTTACATtgt                                                                                                                  >  1:374831/1‑140 (MQ=255)
              tgGCAGATCAAATAATCCCATCACTGACTAAATTGCGTTTCAGCGAACTGGAACATTAATGATTTTTTGTGGAGAAGACGCGTGTGATTGTTAAACGACCCGTCTCGGCCAGTCTGGCCCGGGCCTTTTTTTTACATtgt                                                                                                                  >  1:694033/1‑140 (MQ=255)
               ggCAGATCAAATAATCCCATCACTGACTAAATTGCGTTTCAGCGAACTGGAACATTAATGATTTTTTGTGGAGAAGACGCGTGTGATTGTTAAACGACCCGTCTCGGCCAGTCTGGCCCGGGCCttttttttaca                                                                                                                     <  1:542879/135‑5 (MQ=255)
               ggCAGATCAAATAATCCCATCACTGACTAAATTGCGTTTCAGCGAACTGGAACATTAATGATTTTTTGTGGAGAAGACGCGTGTGATTGTTAAACGACCCGTCTCGGCCAGTCTGGCCCGGGCCttttttttaca                                                                                                                     >  2:542879/1‑131 (MQ=255)
               ggCAGATCAAATAATCCCATCACTGACTAAATTGCGTTTCAGCGAACTGGAACATTAATGATTTTTTGTGGAGAAGACGCGTGTGATTGTTAAACGACCCGTCTCGGCCAGTCTGGCCCGGGCCttttttttaca                                                                                                                     >  1:230100/1‑131 (MQ=255)
               ggCAGATCAAATAATCCCATCACTGACTAAATTGCGTTTCAGCGAACTGGAACATTAATGATTTTTTGTGGAGAAGACGCGTGTGATTGTTAAACGACCCGTCTCGGCCAGTCTGGCCCGGGCCttttttttaca                                                                                                                     <  2:230100/135‑5 (MQ=255)
                 cAGATCAAATAATCCCATCACTGACTAAATTGCGTTTCAGCGAACTGGAACATTAATGATTTTTTGTGGAGAAGACGCGTGTGATTGTTAAACGACCCGTCTCGGCCAGTCTGGCCCGGGCCCTTTTTTTACATTGtgct                                                                                                               >  1:577414/1‑140 (MQ=255)
                   gATCAAATAATCCCATCACTGACTAAATTGCGTTTCAGCGAACTGGAACATTAATGATTTTTTGTGGAGAAGACGCGTGTGATTGTTAAACGACCCGTCTCGGCCAGTCTGGCCCGGGCCTTTTTTTTACATTGtgctgc                                                                                                             <  2:443810/140‑1 (MQ=255)
                    aTCAAATAATCCCATCACTGCCTAAATTGCGTTTCAGCGAACTGGAACATTAATGATTTTTTGTGGAGAAGACGCGTGTGATTGTTAAACGACCCGTCTCGGCCAGTCTGGCCCGGGCCtttttt                                                                                                                          <  2:593529/125‑1 (MQ=255)
                    aTCAAATAATCCCATCACTGACTAAATTGCGTTTCAGCGAACTGGAACATTAATGATTTTTTGTGGAGAAGACGCGTGTGATTGTTAAACGACCCGTCTCGGCCAGTCTGGCCCGGGCCTTTTTTTTACATTGTgctgct                                                                                                            >  1:287831/1‑140 (MQ=255)
                    aTCAAATAATCCCATCACTGACTAAATTGCGTTTCAGCGAACTGGAACATTAATGATTTTTTGTGGAGAAGACGCGTGTGATTGTTAAACGACCCGTCTCGGCCAGTCTGGCCCGGGCCtttttt                                                                                                                          >  1:593529/1‑125 (MQ=255)
                    aTCAAATAATCCCATCACTGACTAAATTGCGTTTCAGCGAACTGGAACATTAATGATTTTTTGTGGAGAAGACGCGTGTGATTGTTAAACGACCCGTCTCGGCCAGTCTGGCCCGGGCCTTTCTTTTACATTGTgctgct                                                                                                            >  1:287830/1‑140 (MQ=255)
                              cccATCACTGACTAAATTGCGTTTCAGCGAACTGGAACATTAATGATTTTTTGTGGAGAAGACGCGTGTGATTGTTAAACGACCCGTCTCGGCCAGTCTGGCCCGGGCCTTTTTTTTACATTGTGCTGCTGTCGAtt                                                                                                     >  2:584109/1‑137 (MQ=255)
                              cccATCACTGACTAAATTGCGTTTCAGCGAACTGGAACATTAATGATTTTTTGTGGAGAAGACGCGTGTGATTGTTAAACGACCCGTCTCGGCCAGTCTGGCCCGGGCCTTTTTTTTACATTGTGCTGCTGTCGAtt                                                                                                     <  1:584109/137‑1 (MQ=255)
                                 aTCACTGACTAAATTGCGTTTCAGCGAACTGGAACATTAATGATTTTTTGTGGAGAAGACGCGTGTGATTGTTAAACGACCCGTCTCGGCCAGTCTGGCCCGGGCCTTTTTTTTACATTGTGCTGCTGTCGATTCTTTcc                                                                                               >  1:347019/1‑140 (MQ=255)
                                  tCACTGACTAAATTGCGTTTCAGCGAACTGGAACATTAATGATTTTTTGTGGAGAAGACGCGTGTGATTGTTAAACGACCCGTCTCGGCCAGTCTGGCCCGGGCCTTTTTTTTACATTGTGCTGCTGTCGATTCTTTCCa                                                                                              >  2:35436/1‑140 (MQ=255)
                                  tCACTGACTAAATTGCGTTTCAGCGAACTGGAACATTAATGATTTTTTGTGGAGAAGACGCGTGTGATTGTTAAACGACCCGTCTCGGCCAGTCTGGCCCGGGCCTTTTTTTTACATTGTGCTGCTGTCGATTCTTTCCa                                                                                              <  1:205268/140‑1 (MQ=255)
                                   cACTGACTAAATTGCGTTTCAGCGAACTGGAACATTAATGATTTTTTGTGGAGAAGACGCGTGTGATTGTTAAACGACCCGTCTCGGCCAGTCTGGCCCGGGCCTTTTTTTACATTGTGCTGCTGTCGATTCTTTCCACg                                                                                            >  1:369154/1‑140 (MQ=255)
                                        aCTAAATTGCGTTTCAGCGAACTGGAACATTAATGATTTTTTGTGGAGAAGACGCGTGTGATTGTTAAACGACCCGTCTCGGCCAGTCTGGCCCGGGCCTTTTTTTTACATTGTGCTGCTGTCGATTCTTTCCACGGGTa                                                                                        >  1:205993/1‑140 (MQ=255)
                                        aCTAAATTGCGTTTCAGCGAACTGGAACATTAATGATTTTTTGTGGAGAAGACGCGTGTGACTGTTAAACGACCCGTCTCGGCCAGTCTGGCCCGGGCCTTTTTTTTACATTGTGCTGCTGTCGATTCTTTCCACGGGTa                                                                                        >  2:287719/1‑140 (MQ=255)
                                           aaaTTGCGTTTCAGCGAACTGGAACATTAATGATTTTTTGTGGAGAAGACGCGTGTGATTGTTAAACGACCCGTCTCGGCCAGTCTGGCCCGGGCCTTTTTTTTACATTGTGCTGCTGTCGATTCTTTCCACGGGTATCg                                                                                     <  2:622006/140‑1 (MQ=255)
                                            aaTTGCGTTTCAGCGAACTGGAACATTAATGATTTTTTGTGGAGAAGACGCGTGTGATTGTTAAACGACCCGTCTCGGCCAGTCTGGCCCGGGCCTTTTTTTTACATTGTGCTGCTGTCGAtt                                                                                                     <  2:265322/123‑1 (MQ=255)
                                            aaTTGCGTTTCAGCGAACTGGAACATTAATGATTTTTTGTGGAGAAGACGCGTGTGATTGTTAAACGACCCGTCTCGGCCAGTCTGGCCCGGGCCTTTTTTTTACATTGTGCTGCTGTCGAtt                                                                                                     >  1:265322/1‑123 (MQ=255)
                                            aaTTGCGTTTCAGCGAACTGGAACATTAATGATTTTTTGTGGAGAAGACGCGTGTGATTGTTAAACGACCCGTCTCGGCCAGTCTGGCCCGGGCCTTTTTTTTACATTGTGCTGCTGTCGATTCTTTCCACGGGTATCGc                                                                                    <  1:35436/140‑1 (MQ=255)
                                                  gTTTCAGCGAACTGGAACATTAATGATTTTTTGTGGAGAAGACGCGTGTGATTGTTAAACGACCCGTCTCGGCCAGTCTGGCCCGGGCCTTTTTTTTACATTGTGCTGCTGTCGATTCTTTCCACGGGTATCGCTCTGCt                                                                              >  1:273788/1‑140 (MQ=255)
                                                              tGGAACATTAATGATTTTTTGTGGAGAAGACGCGTGTGATTGTTAAACGACCCGTCTCGGCCAGTCTGGCCCGGGCCTTTTTTTTACATTGTGCTGCTGTCGATTCTTTCCACGGGTATCGCTCTGCTAACTCTGGCGAg                                                                  <  1:130303/140‑1 (MQ=255)
                                                                      tAATGATTTTTTGTGGAGAAGACGCGTGTGATTGTTAAACGACCCGTCTCGGCCAGTCTGGCCCGGGCCTTTTTTTTACATTGTGCTGCTGTCGATTCTTTCCACGGGTATCGCTCTGCTAACTCTGGCGAGCAGTTTgc                                                          >  1:418679/1‑140 (MQ=255)
                                                                       aaTGATTTTTTGTGGAGAAGACGCGTGTGATTGTTAAACGACCCGTCTCGGCCAGTCTGGCCCGGGCCTTTTTTTTACATTGTGCTGCTGTCGATTCTTTCCACGGGTATCGCTCTGCTAACTCTGGCGAGCAGTTTgcg                                                         <  2:418679/140‑1 (MQ=255)
                                                                        aTGATTTTTTGTGGAGAAGACGCGTGTGATTGTTAAACGACCCGTCTCGGCCAGTCTGGCCCGGGCCTTTTTTTTACATTGTGCTGCTGTCGATTCTTTCCACGGGTATCGCTCTGCTAACTCTGGCGAGCAGTTTgcgc                                                        >  1:650234/1‑140 (MQ=255)
                                                                        aTGATTTTTTGTGGAGAAGACGCGTGTGATTGTTAAACGACCCGTCTCGGCCAGTCTGGCCCGGGCCTTTTTTTTACATTGTGCTGCTGTCGATTCTTTCCACGGGTATCGCTCTGCTAACTCTGGCGAGCAGTTTgcgc                                                        >  2:319414/1‑140 (MQ=255)
                                                                              ttttGTGGAGAAGACGCGTGTGATTGTTAAACGACCCGTCTCGGCCAGTCTGGCCCGGGCCTTTTTTTTACATTGTGCTGCTGTCGATTCTTTCCACGGGTATCGCTCTGCTAACTCTGGCGAGCAGTTTGCGCGACGCt                                                  <  2:273788/140‑1 (MQ=255)
                                                                                 tgtgGAGAAGACGCGTGTGATTGTTAAACGACCCGTCTCGGCCAGTCTGGCCCGGGCCTTTTTTTTACATTGTGCTGCTGTCGATTCTTTCCACGGGTATTGCTCTGCTAACTCTGGCGAGCAGTTTGCGCGACGCTGAg                                               <  1:241799/140‑1 (MQ=255)
                                                                                 tgtgGAGAAGACGCGTGTGATTGTTAAACGACCCGTCTCGGCCAGTCTGGCCCGGGCCTTTTTTTTACATTGTGCTGCTGTCGATTCTTTCCACGGGTATCGCTCTGCTAACTCTGGCGAGCAGTTTGCGCGACGCTGAg                                               <  2:694033/140‑1 (MQ=255)
                                                                                    ggAGAAGACGCGTGTGATTGTTAAACGACCCGTCTCGGCCAGTCTGGCCCGGGCCCTTTTTTTACATTGTGCTGCTGTCGATTCTTTCCACGGGTATCGCTCTGCTAACTCTGGCGAGCAGTTTGCGCGACGCTGAGGCt                                            <  2:577414/140‑1 (MQ=255)
                                                                                     gagaAGACGCGTGTGATTGTTAAACGACCCGTCTCGGCCAGTCTGGCCCGGGCCTTTTTTTTACATTGTGCTGCTGTCGATTCTTTCCACGGGTATCGCTCTGCTAACTCTGGCGAGCAGTTTGCGCGACGCTGAGGCTa                                           <  2:287831/140‑1 (MQ=255)
                                                                                     gagaAGACGCGTGTGATTGTTAAACGACCCGTCTCGGCCAGTCTGGCCCGGGCCTTTTTTTTACATTGTGCTGCTGTCGATTCTTTCCACGGGTATCGCTCTGCTAACTCTGGCGAGCAGTTTGCGCGACGCTGAGGCTa                                           >  2:268654/1‑140 (MQ=255)
                                                                                     gagaAGACGCGTGTGATTGTTAAACGACCCGTCTCGGCCAGTCTGGCCCGGGCCTTTCTTTTACATTGTGCTGCTGTCGATTCTTTCCACGGGTATCGCTCTGCTAACTCTGGCGAGCAGTTTGCGCGACGCTGAGGCTa                                           <  2:287830/140‑1 (MQ=255)
                                                                                              cgTGTGATTGTTAAACGACCCGTCTCGGCCAGTCTGGCCCGGGCCTTTTTTTTACATTGTGCTGCTGTCGATTCTTTCCACGGGTATCGCTCTGCTAACTCTGGCGAGCAGTTTGCGCGACGCTGAGGCTATCAATATTg                                  <  2:347019/140‑1 (MQ=255)
                                                                                               gtgtGATTGTTAAACGACCCGTCTCGGCCAGTCTGGCCCGGGCCTTTTTTTTACATTGTGCTGCTGTCGATTCTTTCCACGGGTATCGCTCTGCTAACTCTGGCGAGCAGTTTGCGCGACGCTGAGGCTATCAATATTGc                                 >  2:179573/1‑140 (MQ=255)
                                                                                               gtgtGATTGTTAAACGACCCGTCTCGGCCAGTCTGGCCCGGGCCTTTTTTTTACATTGTGCTGCTGTCGATTCTTTCCACGGGTATCGCTCTGCTAACTCTGGCGAGCAGTTTGCGCGACGCTGAGGCTATCAATATTGc                                 >  2:124154/1‑140 (MQ=255)
                                                                                                tgtgATTGTTAAACGACCCGTCTCGGCCAGTCTGGCCCGGGCCTTTTTTTTACATTGTGCTGCTGTCGATTCTTTCCACGGGTATCGCTCTGCTAACTCTGGCGAGCAGTTTGCGCGACGCTGAGGCt                                            >  1:103924/1‑128 (MQ=255)
                                                                                                tgtgATTGTTAAACGACCCGTCTCGGCCAGTCTGGCCCGGGCCTTTTTTTTACATTGTGCTGCTGTCGATTCTTTCCACGGGTATCGCTCTGCTAACTCTGGCGAGCAGTTTGCGCGACGCTGAGGCt                                            <  2:103924/128‑1 (MQ=255)
                                                                                                      tGTTAAACGACCCGTCTCGGCCAGTCTGGCCCGGGCCTTTTTTTTACATTGTGCTGCTGTCGATTCTTTCCACGGGTATCGCTCTGCTAACTCTGGCGAGCAGTTTGCGCGACGCTGAGGCTATCAATATTGCCGGATCg                          <  2:537131/140‑1 (MQ=255)
                                                                                                         tAAACGACCCGTCTCGGCCAGTCTGGCCCGGGCCTTTTTTTTACATTGTGCTGCTGTCGATTCTTTCCACGGGTATCGCTCTGCTAACTCTGGCGAGCAGTTTGCGCGACGCTGAGGCTATCAATATTGCCGGATCGCTg                       <  1:52006/140‑1 (MQ=255)
                                                                                                               aCCCGTCTCGGCCAGTCTGGCCCGGGCCTTTTTTTTACATTGTGCTGCTGTCGATTCTTTCCACGGGTATCGCTCTGCTAACTCTGGCGAGCAGTTTGCGCGACGCTGAGGCTATCAATATTGCCGGATCGCTGCGTATg                 >  1:860359/1‑140 (MQ=255)
                                                                                                               aCCCGTCTCGGCCAGTCTGGCCCGGGCCTTTTTTTTACATTGTGCTGCTGTCGATTCTTTCCACGGGTATCGCTCTGCTAACTCTGGCGAGCAGTTTGCGCGACGCTGAGGCTATCAATATTGCCGGATCGCTGCGTATg                 <  1:287719/140‑1 (MQ=255)
                                                                                                               aCCCGTCTCGGCCAGTCTGGCCCGGGCCTTTTTTTTACATTGTGCTGCTGTCGATTCTTTCCACGGGTATCGCTCTGCTAACTCTGGCGAGCAGTTTGCGCGACGCTGAGGCTATCAATATTGCCGGATCGCTGCGTATg                 <  1:268654/140‑1 (MQ=255)
                                                                                                                  cGTCTCGGCCAGTCTGGCCCGGGCCTTTTTTTTACATTGTGCTGCTGTCGATTCTTTCCACGGGTATCGCTCTGCTAACTCTGGCGAGCAGTTTGCGCGACGCTGAGGCTATCAATATTGCCGGATCGCTGCGTATGCag              >  1:13264/1‑140 (MQ=255)
                                                                                                                       cGGCCAGTCTGGCCCGGGCCTTTTTTTACATTGTGCTGCTGTCGATTCTTTCCACGGGTATCGCTCTGCTAACTCTGGCGAGCAGTTTGCGCGACGCTGAGGCTATCAATATTGCCGGATCGCTGCGTATGCAGAGTTAc        <  2:369154/140‑1 (MQ=255)
                                                                                                                              tCTGGCCCGGGCCTTTTTTTTACATTGTGCTGCTGTCGATTCTTTCCACGGGTATCGCTCTGCTAACTCTGGCGAGCAGTTTGCGCGACGCTGAGGCTATCAATATTGCCGGATCGCTGCGTATGCAGAGTTACCGCCTg  >  2:342681/1‑140 (MQ=255)
                                                                                                                                            |                                                                                                                               
TTCGCTATAAACTGTGGCAGATCAAATAATCCCATCACTGACTAAATTGCGTTTCAGCGAACTGGAACATTAATGATTTTTTGTGGAGAAGACGCGTGTGATTGTTAAACGACCCGTCTCGGCCAGTCTGGCCCGGGCCTTTTTTTACATTGTGCTGCTGTCGATTCTTTCCACGGGTATCGCTCTGCTAACTCTGGCGAGCAGTTTGCGCGACGCTGAGGCTATCAATATTGCCGGATCGCTGCGTATGCAGAGTTACCGCCTG  >  NC_000913/2585634‑2585898

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