Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 2,321,822 T→C 12.3% F594L (TTC→CTC)  atoS → sensory histidine kinase AtoS

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009132,321,8220TC12.3% 76.6 / 4.5 45F594L (TTC→CTC) atoSsensory histidine kinase AtoS
Reads supporting (aligned to +/- strand):  ref base T (18/21);  new base C (3/3);  total (21/24)
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.46e-01

CGCTGCAAAAAAAGATCTTCGATCCCTTTTTCACCACCAAAGCCTCAGGAACCGGGCTTGGTCTGGCGTTAAGTCAACGCATCATTAATGCCCATCAGGGTGATATTCGCGTCGCCAGTTTGCCGGGCTACGGCGCAACCTTCACGCTTATTTTACCGATCAACCCGCAGGGAAATCAGACTGTATGACTGCTATTAATCGCATCCTTATTGTGGATGATGAAGATAATGTTCGCCGTATGCTGAGCACCGCTTTTGCACTACAAGGATTCGAAACACATTGT  >  NC_000913/2321682‑2321964
                                                                                                                                            |                                                                                                                                              
cGCTGCAAAAAAAGATCTTCGATCCCTTTTTCACCACCAAAGCCTCAGGAACCGGGATTGGGCTGGGGTTAAGTCAACGCATCATTAATGCGCATCGGGGTGATGTTCGCGTCGCAAGTTTGCGGGGCGACGGCGCAGCCtgt                                                                                                                                              >  1:219840/1‑141 (MQ=255)
    gCAAAAAAAGATCTTCGATCCCTTTTTCACCACCAAAGCCTCAGGAACCGGGCTTGGTCTGGCGTTAAGTCAACGCATCATTAATGCCCATCAGGGTGATATTCGCGTCGCCAGCTTGGCGGGCTACGGCGGAGCCTTCACGc                                                                                                                                          >  2:236265/1‑143 (MQ=255)
                 ttCGATCCCTTTTTCACCACCAAAGCCTCAGGAACCGGGCTTGGTCTGGCGTTAAGTCAACGCATCATTAATGCCCATCAGGGTGATATTCGCGTCGCCAGTTTGCCGGGCTACGGCGCAACCTTCACGCTTATTTTACCGAt                                                                                                                             >  2:36764/1‑143 (MQ=255)
                 gtggATCCCTCTGTCAGCACCAAAGGCTCAGGAACCGGGCTTGGTCTGGCGTTAAGTCAACGCATCATTAATGCCCATCAGGGTGATATTCGCGTCGCCAGTTTGCCGGGCTACGGCGCAACCTTCACGCTTATTTTACCGAt                                                                                                                             <  1:504528/140‑1 (MQ=255)
                 gtCGATCCCTTTTTCACCCCCAAGTCCTCAGGAACCGGGCTTGGTCTGGCGTTAAGTCAACGCATCATTAATGCCCATCAGGGTGATATTCTCGTCGCCAGTTTGCCGGGCTACGGCGCAACCTTCACGCTTATTTTACCGAt                                                                                                                             <  2:396076/142‑1 (MQ=255)
                     aTCCCTTTTTCACCACCAAAGCCTCAGGAACCGGGCTTGGTCTGGCGTTAAGTCAACGCATCATTAATGCCCATCAGGGTGATATTCGCGTCGCCAGTTTGCCGGGCTACGGCGCAACCCTCACGCTTATTTTACCGATCAAc                                                                                                                         >  1:666499/1‑143 (MQ=255)
                     aTCCCTTTTTCACCACCAAAGCCTCAGGAACCGGGCTTGGTCTGGCGTTAAGTCAACGCATCATTAATGCCCATCAGGGTGATATTCGCGTCGCCAGTTTGCCGGGCTACGGCGCAACCCTCACGCTTATTTTACCGATCAAc                                                                                                                         >  1:666497/1‑143 (MQ=255)
                     aTCCCTTTTGCACCACCAAAGCCTCAGGAACCGGGCTTGGTCTGGCGTTAAGTCAACGCATCATTAATGCCCATCAGGGTGATATTCGCGTCGCCAGTTTGCCGGGCTACGGCGCAACCTTCACGCTTATTTTACCGATCAAc                                                                                                                         <  2:246430/143‑1 (MQ=255)
                           ttttCACCACCAAAGCCTCAGGAACCGGGCTTGGTCTGGCGTTAAGTCAACGCATCCTTAATGCCCATCAGGGTGATATTTGCGTCGCCCCTCTCCCGTGCCACGGCGCAACCTTCACGCTTATTTTTCCGCTCAACCCGcgg                                                                                                                   >  1:598878/1‑141 (MQ=255)
                           tcttCACCCACAATTCCTCAGGATCCGGGCTTGGTCTGGCGTTAAGTCCACGCATGATTAATGCCGATCTGGGTGATATTCGCTTCGCCAGTTTGCCGGGCTACGGCGCAACCTTCACGCTTATTTTACCGATCAACCCGCAg                                                                                                                   <  2:74175/141‑1 (MQ=255)
                              tCACCTCCATGGCCGCCGTAATCGGGCTTTGTCTGGCGTTCAGTGAACGCATGATTAAGGCCGATCAGGGTGTTATTCGCGTCGCCAGTTTGCCGGTCTACGGCGCAACCTTCACTCTTATTTTACCGATCAACCCGCAGGGa                                                                                                                <  2:262709/143‑1 (MQ=255)
                              tCACCACCAAAGCCTCAGGAACCGGGCTTGGTCTGGCGTTAAGTCAACGCATCATTAATGCCCATCAGGGTGATATTCGCGTCGCCAGTTTGCCGGGCTACGGCGCAACCCTCACGCTTATTTTACCGATCAACCCGCAGGGa                                                                                                                <  2:666499/143‑1 (MQ=255)
                              tCACCACCAAAGCCTCAGGAACCGGGCTTGGTCTGGCGTTAAGTCAACGCATCATTAATGCCCATCAGGGTGATATTCGCGTCGCCAGTTTGCCGGGCTACGGCGCAACCCTCACGCTTATTTTACCGATCAACCCGCAGGGa                                                                                                                <  2:666497/143‑1 (MQ=255)
                                    ccaAAGCCTCAGGAACCGGGCTTGGTCTGGCGTTAAGTCAACGCATCATTAATGCCCATCAGGGTGATATTCGCGTCGCCAGTTTGCCGGGCTACGGCGCAACCTTCACGCTTATTTTACCGATCAACCCGCAGGGCAATCAg                                                                                                          >  2:176485/1‑143 (MQ=255)
                                              aGGAACCGGGCTTGGTCTGGCATTAAGTCAACGCATCATTAATGCCCATCAGGGTGATATTAGCGTCGCCAATTTGCCGGGCTACGGCGCAACCTTCACGCCTCTGTTACCGATCAACCCGCAACGAAATCAGACCGTACGAc                                                                                                >  2:663857/1‑143 (MQ=255)
                                                 aaCCGGGCTTGGTCTGGCGTTAAGTCAACGCATCATTAATGCCCATCAGGGTGATATTCGCGTCGCCAGTTTGCCGGGCTACGGCGCAACCTTCACGCTTATTTTACCGATCAACCGGCAGGGAAATCGGCCTGGATGCCTGc                                                                                             >  1:389290/1‑143 (MQ=255)
                                                    cGGGCTTGGTCTGGCGTTAAGTCAACGCATCATTAATGCCCATCAGGGTGATTTTACTGCTGCCAGGAGGCCGGGCTACGGCGTAACCTTCACGCTTATTTTACCGATCAACCCGCAGGGAAATCAGACTGTATGACTGCTAt                                                                                          >  1:141150/1‑143 (MQ=255)
                                                    cGGGCTTGGTCTGGCGTTAAGTCAACGCATCATTAATGCCCATCAGGGTGATATTCGCGTCGCCAGTTTGCCGGGCTACGGCGCAACATTCACGCTTATTTTACCGATCAACCCGCAGGGAAATCAGACTGTATGACTGCTAt                                                                                          >  2:364784/1‑143 (MQ=255)
                                                    cGGGCTTGGTCTGGCGTTAAGTCAACGCATCATTAATGCCCATCAGGGTGATATTCGCGTCGCCAGTTTGCCGGGCTAAGGCGCAACCTTCACGCTTACGTTACCGATCAACCCGCAGGGAAATCAGACGGTATGACTGCTAt                                                                                          >  2:619846/1‑143 (MQ=255)
                                                    cGGGCTTGGTCTGGCATTAAGTCAACGCATCATTAATGCCCATCAGGGTGATATTCGCGTGACCAGTTTGCCGCGCTACGGCGCAACCTTCACGCATATGTTACCGATCAACCCGAAGAGAAAGCAGACTGTAAGCCTGcaaa                                                                                          >  2:222159/1‑140 (MQ=255)
                                                      ggCTTGGTCTGGCGTTAAGTCAACGCATCAGTAATGCCCAGCAGGGTGATATTCGCGTCGCCAGTTTGCCGGGCTACGGCGCAACCTTCACGCTTATTTTACCGATCAACCCGCAGGGAAATCAGACTGTATGACTGCTATTa                                                                                        <  1:734483/143‑1 (MQ=255)
                                                      ggCTTGGTCTGCCGTTAAGTCAACGCCTCATTAATGCCCATCAGGGTGATGTTAGGTACGCCATTTTGCCGGGCTACGGCGCAACCTTCACGCTTATTTTACCGATCAACCCGCAGGGAAATCAGACTGTATGACTGCTATTa                                                                                        <  1:172154/143‑1 (MQ=255)
                                                        cTTGGGCTGGCGTTAAGTCAACGCATCATTAATGCCCATCAGGGTGATATTCGCGTCGCCAGTTTGCCGTGCTACGGCGGAATCTTTTTGCTTATTTTGCCGATCAACCCGCAGGGAAATCAGACTGTATGACTGCTATTAAt                                                                                      <  2:125901/143‑1 (MQ=255)
                                                         ttGGTCTGGCGTTAAGTCAACGCATCATTAATGCCCATCAGGGTGATATTCGCGTCGCCAGTTTGCCGGGCTACGGCGCAACCTTCACGCTTATTTTACCGATCAACCCGCAGGGAAATCAGAATGTATGACTGCTATTAATc                                                                                     >  1:418065/1‑143 (MQ=255)
                                                                         tCAACGCATCATTAATGCCCATCAGGGTGATATTCGCGTCGCCAGTTTGACGGGCTACGGCGCAACCCTCACGCATATTTTACCGATCAACCCACAGGGAAATCAGACTGTATGACTGCCATTAATCACATACTAATAGTGga                                                                     >  2:443321/1‑143 (MQ=255)
                                                                            aCGCATCATTAAGGCCCATCAGGGTGATATTCGCGTCGCCAGTTTGCCGGGCTACGGCGCAACCTTCACGCTTATTTTACCGATCAACCCGCAGGGAAATCAGACTGTATGACTGCTATTAATCGCATCCTTATTGTGgatga                                                                  <  1:551881/143‑1 (MQ=255)
                                                                                 tcatTAATGCCCATCAGGGTGATATTCGCGTCGCCAGTTTGCCGGGCTACGGCGCAACCTTCAAGCATAGGGTTCCGATCAACCCGCAGGGAAATCAGACTGTCTGCCTGCTATTACTTGCTTCCTTATTGTTGATGATGAAg                                                             >  1:480425/1‑143 (MQ=255)
                                                                                   atTAATGCCCATCAGGGTGATATTCGCGTCGCCAGTTATCCAGGCTACACCGCAACCTTCACGCTTATGTTACCGATCAACCCGCAGAGAAATCAGACTGTATGACTCCTATTAACCGCAACCTTATCGTGCGCGATGAAgct                                                           >  2:506142/1‑141 (MQ=255)
                                                                                         gggcATCAGGGGCATTCTAGCGTCGGCAGTTTTCCGGGGTACGGCGGAACTTTCACGCTTATTTTACCGATCAACCCGCAAGGAAATCAGACTGTATGACTGCTATTAATCGCATCCTTATTGTGGATGATGAAGATAATGtt                                                     <  2:768747/140‑1 (MQ=255)
                                                                                              tCAGGGTGATATTCGCGTCGCCAGTTTGCCGGGCTACGGCGCAACCTTCACGCTTATTTTACCGATCAACCCGCAGGGAAATCAGACTGTATGACTGCTATTAATCGCATCCTTACTGTGGATCATGAAGATAATGTTAGCCg                                                >  2:656977/1‑143 (MQ=255)
                                                                                                      atatTCGCGTCGCCAGTTTGCCGGGCTACGGCGCAACCTTCACGCTTATTTTACCGATCAACCCGCAGGGAAATCAGACTGTATGACTGCTATTAATCGCATCCTTATTGTGGATGATGAAGATAATGTTCGCCGTATGCTGa                                        <  1:36764/143‑1 (MQ=255)
                                                                                                          tCGCGTCGCCAGTTTGCCGGGCTACGGCGCAACCTTCACGCTTATTTTACCGATCAACCCGCAGGGAAATCAGACTGTATGACCGCTATTAATCGCATCCTTCTTGTGGAAGCTGAAGATAATGTTCGCCGTAAGCTcaccac                                    >  2:21094/1‑143 (MQ=255)
                                                                                                            gcgTCGCTAGTTTGCCGAGCTACGGCGCAACCTTCACGCTTATTTTACCGATCAACCCGCAGGGAAATCAGACTGTATGACTGCTATTAATCGCATCCTTATTGTGGATGATGAAGATAATGTTCGCCGTATGCTGAGCACCg                                  <  1:614279/143‑1 (MQ=255)
                                                                                                             cgTCGCCAGTTTGCCGGGCTACGGCGCAACCTTCACGCTTATTTTACCGATCAACCCGCAGGGAAATCAGACTGTATGACTGCTATTAATCGCATCCTTATTGTGGATGATGAAGATACTGTTCGCCGTATGCTGAGCACCGc                                 >  1:750614/1‑143 (MQ=255)
                                                                                                              gTCGCCAGTTTGCCGTGCTACGGCGTAACCTTCTCGTTTATTTTACCTATCATTCCGCAGGGTAATCAGATTGTTTGACTGCTATTAATCGCATCCTTATTGTGGATGATGAAGATAATGTGCGCCGTATGCTGAGCACCGCt                                <  1:397873/143‑1 (MQ=255)
                                                                                                                 gCCAGTTTGCCGGGCTACGGCGCAACCTTCACGCTTATTTTACCGATCAACCCGCAGGGAAATCAGACTGTATGACTGCTATTAATCGCATCCTTATCGTGGATGATGAAGATAATGTTCGCCGTATGCTGAGCACCGCtttt                             >  2:47570/1‑143 (MQ=255)
                                                                                                                    aGTTTGCCGGGCTACGGCGCTACCTTCACGATTATTTTACCGATCAGCCCGCAGAGAAAGAAGACTGTAGGACGGCTATTAATCGCATCCTTATTGTGGATGATGAAGATAATGTTCGCCGTATGCTGAGCACCGCTTTTGCa                          <  1:189844/143‑1 (MQ=255)
                                                                                                                     gcatGCCGGGCAACGGCGCCATCCGCACGCTTATTGTACCGATCAACCCGCAGGGAAATCAGACTGTATTACTGCTATTAATCGCATCCTTATTGTGGATGATGAAGATAATGTTCGCCGTATGCTGAGCACCGCTTTTGCAc                         <  1:236265/140‑1 (MQ=255)
                                                                                                                     gTTTGCCGGGCTACGTCGCAACCTTCACGCTTATTTTACCGATCAACCCGCAGGGAAATCAGACTGTATGACTGCTATTAATCGCATCCTTATTGTGGATGATGAAGATAATGTTCGCCGTATGCTGAGCACCGCTTTTGCAc                         <  2:306382/143‑1 (MQ=255)
                                                                                                                      ttggcggggCTAGGGCGCAACCTTCTCGCTTATTTTACCGGTCAAGCCGCAGGGAAATCAGACTGTATGACTGCTAGTAATCGCATCCTTATTGTGGATGATGAAGATAATGTTCGCCGTATGCTGAGCACCGCTTTTGCACt                        <  2:157363/137‑1 (MQ=255)
                                                                                                                          ccGGGCTACGGCGCATCCTTCACGCTTATTTTACCGATCAACCCGCAGGGAAATCAGACTGTATGACTGCTATTAATCGCATCCTTATTGTGGATGATGAAGATAATGTTCGCCGTATGCTGAGCACCGCTTTTGCACTACaa                    <  1:364784/143‑1 (MQ=255)
                                                                                                                            gggCTACGGCGTAACCTTCACGTTTATTTTAATGATCAACACGAAGGGAAATAAGACTGTATGACTGCTATGAATCGCATCCTTATTGTGGATGATGAAGATAATGTTCGCCGTATGCTGAGCACCGCTTTTGCACTACAAgg                  <  1:271694/143‑1 (MQ=255)
                                                                                                                             ggCTACGGCGCAACCTTCACGCTTATTTTACCGATCATCCCGCAGGGAAATCAGACTGTATGACTGCTATTAATCGCATCCTTATTGTGGATGATGAAGATAATGTTCGCCGTATGCTGAGCACCGCTTTTGCACTACAAGGa                 <  1:579153/143‑1 (MQ=255)
                                                                                                                                tACGGCGCAACCTTCACGCTTATTTTACCGATCAACCCGCAGGGAAATCAGACTGTATGACTGCTATTAATCGCATCCTTATTGTGGATGATGAAGATAATGTTCGCCGTATGCTGAACACCGCTTTTGCACTACAAGGACTc              >  2:272687/1‑143 (MQ=255)
                                                                                                                                   tccacAACCTTCACGCTTATTTTACCGATCAACCCGCAGGGAAATCAGACTGTATGACTGCTATTAATCGCATCCTTATTGTGGATGATGAAGATAATGTTCGCCGTATGCTGAGCACCGCTTTTGCACTACAAGGATTCGaa           <  1:96179/139‑1 (MQ=255)
                                                                                                                                    gcgcAACCTTCACGCTTATTTTACCGATCAACCCGCAGGGAAATCAGACTGTATGACTGCTATTAATCGCATCCTTATTGTGGATGATGAAGATAATGTTCGCCGTATGCTGAGCACCGCTTTTGCACTACAAGGATTCGaaa          <  2:480425/143‑1 (MQ=255)
                                                                                                                                            ttCACGCTTATTTTACCGATCAACCCGCAGGGAAATCAGACTGTATGACTGCTATTAATCGCATCCTTATTGTGGATGATGAAGATAACGTTCGCCGGATGCTGAGCACCGCTTTCGCACTCCAAGGGTTCGacactcctgct  >  1:432741/1‑136 (MQ=255)
                                                                                                                                            ttCACGCTTATTTTACCGATCAACCCGCAGGGAAATCAGACTGTATGACTGCTATTAATCGCATCCTAATTGTGGATGATGAAGATAATGTTCGCCGTATGCTGAGCACCGCCTTTGCACTACAAGGATTCGAAACACATtgt  >  2:9102/1‑143 (MQ=255)
                                                                                                                                            |                                                                                                                                              
CGCTGCAAAAAAAGATCTTCGATCCCTTTTTCACCACCAAAGCCTCAGGAACCGGGCTTGGTCTGGCGTTAAGTCAACGCATCATTAATGCCCATCAGGGTGATATTCGCGTCGCCAGTTTGCCGGGCTACGGCGCAACCTTCACGCTTATTTTACCGATCAACCCGCAGGGAAATCAGACTGTATGACTGCTATTAATCGCATCCTTATTGTGGATGATGAAGATAATGTTCGCCGTATGCTGAGCACCGCTTTTGCACTACAAGGATTCGAAACACATTGT  >  NC_000913/2321682‑2321964

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