Predicted mutation
evidence seq id position mutation annotation gene description
RA NC_000913 1,054,270 C→T A637T (GCC→ACC)  torS ← sensory histidine kinase TorS

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009131,054,2700CT96.4% 79.6 / NA 28A637T (GCC→ACC) torSsensory histidine kinase TorS
Reads supporting (aligned to +/- strand):  ref base C (0/0);  major base T (11/16);  minor base A (1/0);  total (12/16)
Fisher's exact test for biased strand distribution p-value = 4.29e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

AACCGTCAAGACGCACCGCCTGATTGACTGTTTTTGGCACGGGTGCCGTGGCAACACGTAACGGCAAGCGTAAACAAAAACAGCTTCCAACCTCCGGCGTGCTGGTGGCGCTCAGTTCGCCGCCCATCGCCTGGGCCAGACGGCTACTGATAGTCAGTCCCAGCCCGGTGCCGCCGCGTTTGCCGCTTACCTGCACAAATGGCTGGAAGATTTCTGCCAGTTTCGCGGGATCAATACCGCAGCCGCTGTCTTCCACTTCGACCAGCCAT  >  NC_000913/1054135‑1054403
                                                                                                                                       |                                                                                                                                     
aaCCGTCAAGACGCACCGCCTGATTGACTGTTTTTGGCACGGGTGCCGTGGCAACACGTAACGGCAAGCGTAAACAAAAACAGCTTCCAACCTCCGGCGTGCTGGTGGCGCTCAGTTCGCCGCCCATCGCCTGGGTCAGa                                                                                                                                   >  2:604654/1‑140 (MQ=255)
     tCAAGACGCACCGCCTGATTGACTGTTTTTGGCACGGGTGCCGTGGCAACACGTAACGGCAAGCGTAAACAAAAACAGCTTCCAACCTCCGGCGTGCTGGTGGCGCTCAGTTCGCCGCCCATCGCCTGGGTCAGACGGCt                                                                                                                              >  1:115408/1‑140 (MQ=255)
             cACCGCCTGATTGACTGTTTTTGGCACGGGTGCCGTGGCAACACGTAACGGCAAGCGTAAACAAAAACAGCTTCCAACCTCCGGCGTGCTGGTGGCGCTCAGTTCGCCGCCCATCGCCTGGGTCAGACGGCTACTGATag                                                                                                                      >  1:735052/1‑140 (MQ=255)
              aCCGCCTGATTGACTGTTTTTGGCACGGGTGCCGTGGCAACACGTAACGGCAAGCGTAAACAAAAACAGCTTCCAACCTCCGGCGTGCTGGTGGCGCTCAGTTCGCCGCCCATCGCCTGGGTCAGACGGCTACTGATagt                                                                                                                     <  2:172069/140‑1 (MQ=255)
               ccGCCTGATTGACTGTTTTTGGAACGGGTTCCGTGGCCACACGTAACGGCAAGCGTAAACAAAAACAGCTTCCAACCTCCGGCGTGCTGGTGGCGCTCAGTTCGCCGCCCATCGCCTGGGTCAGACGGCTACTGATagtc                                                                                                                    <  1:378840/140‑1 (MQ=255)
                cGCCTGATTGACTGTTTTTGGCACGGGTGCCGTGGCAACACGTAACGGCAAGCGTAAACAAAAACAGCTTCCAACCTCCGGCGTGCTGGTGGCGCTCAGTTCGCCGCCCAACGCTTGTGTCAGATGGCTACTGATagtca                                                                                                                   >  1:286779/1‑140 (MQ=255)
                  ccTGATTGACTGTTTTTGGCACGGGTGCCGTGGCAACACGTAACGGCAAGCGTAAACAAAAACAGCTTCCAACCTCCGGCGTGCTGGTGGCCCTCAGTTCGCCGACCACCGCCTGGGACAGAGGGCTACTGATagtcagt                                                                                                                 >  2:332331/1‑140 (MQ=255)
                   cTGATTGACTGTTTTTGGCACGGGTGCCGTGGCAACACGTAACGGCAAGCGTAAACAAAAACAGCTTCCAACCTCCGGCGTGCTGGTGGCGCTCAGTTCGCCGCCCATCGACTGAGTGAGACGGGTACTGATagtcagtc                                                                                                                >  1:229907/1‑140 (MQ=255)
                            tGTTTTTGGCACGGGATCCGTGGCAACACGTAACGGCAAGCGCAAACAAAAACAGCTTCCAACCTCCGGCGTGCTGGTGGCGATCAGTTCGACGCCCATCGCCTGGGTTAGACGGCTGCTGAGAGTCAGTCCCAGCCCgg                                                                                                       >  2:246885/1‑140 (MQ=255)
                             gTTTTTGGCACGGGTGCCGTGGCAACACGTAACGGCAAGCGTAAACAAAAACAGCTTCCAACCTCCGGCGTGCTGGTGGCGCTCAGTTCGCCGCCCATCGCCTGGGTCAGACGGCTACTGATAGTCAGTCCCAGCCCGGt                                                                                                      >  1:168471/1‑140 (MQ=255)
                                        gggTGCCGTGGCAACACGTAACGGCAAGCGTAAATAAAAACAGCTTCCAACCTCCGGCGTGCTGGTGGCGCTCAGTTCGCCGCCCATCGCCTGGGTCAGACGTCTACTGAAAGTAAGGCCCAGCACGGTGCCGCCGCGtt                                                                                           >  1:397262/1‑140 (MQ=255)
                                        gggTGCCGTGGCAACACGTAACGGCAAGCGTAAACAAAAACAGCTTCCAACCTCCGGCGTGCTGGTGGCGCTCAGTTCGCCGCCCATCGCCTGGGTCAGACGGCTACTGATAGTCAGTCCCAGCCCGGTGCCGCCGCGtt                                                                                           <  1:568899/140‑1 (MQ=255)
                                             ccGTGGCAACACGTACCGGCAAACGCAAACAAAAACACCTTCCAACCTCCGGCGTGCTGGTGGCGCTCAGTTCGCCGCCCATCGCCTGGGTCAGACGGCTACTGATAGTCAGTCCCAGCCCGGTGCCGCCGCGTTTGCCg                                                                                      <  1:246885/140‑1 (MQ=255)
                                             ccGTGGCAACACGTAACGGCAAGCGTAAACAAAAACAGCTTCCAACCTCCGGCGTGCTGGTGGCGCTCAGTTCGCCGCCCATCGCCTGGGTCAGACGGCTACTGATAGTCAGTCCCAGCCCGGTGCCGCCGCGTTTGCCg                                                                                      <  1:604654/140‑1 (MQ=255)
                                                  gCAACACGTAACGGCCAGCGTATACTAAAACAGCTTCTAACCTCCGGCGTGCTGGTGGCGCTCAGTTCGCCGCCCATCGCCTGGGTCAGACGGCTACTGATAGTCAGTCCCAGCCCGGTGCCGCCGCGTTTGCCGCTTAc                                                                                 <  1:329951/140‑1 (MQ=255)
                                                  gCAACACGTAACGGCCAGCGTACACAAAACCAGCTTCAAACCTCCGGCGTGCTGGTGGCGCTCAGTTCGCCGCCCATCGCCTGGGTCAGACGGCTACTGATAGTCAGTCCCAGCCCGGTGCCGCCGCGTTTGCCGCTTAc                                                                                 <  1:263567/140‑1 (MQ=255)
                                                  gCAACACGTAACGGCAAGCGTAAACAAAAACAGCTTCCAACCTCCGGCGTGCTGGTGGCGCTCAGTTCGCCGCCCATCGCCTGGGTCAGACGGCTACTGATAGTCAGTCCCAGCCCGGTGCCGCCGCGTTTGCCGCTTAc                                                                                 <  2:168471/140‑1 (MQ=255)
                                                     acacGTAACGGCAAGTGTAAACGAATACAGCTTCCAAGCTCCGGCGTGCTGGTGGCGCTCAGTTCGCCGCCCATCGCCTGGGTCAGACGGCTACTGATAGTCAGTCCCAGCCCGGTGCCGCCGCGTTTGCCGCTTACCTg                                                                              <  1:332331/140‑1 (MQ=255)
                                                          tAACGGCAAGCGTAAACAAAAACAGCTTCCAACCTCCGGCGTGCTGGTGGCGCTCAGTTCGCCGCCCATCGCCTGGGTCAGACGGCTACTGATAGTCAGTCCCAGCCCGGTGCCGCCGCGTTTGCCGCTTACCTGCACaa                                                                         <  2:115408/140‑1 (MQ=255)
                                                               gCAAGCGTAAATAAAACCAGCTCCCAACCTCCGGCCTGCCGGTGGCGCACAGTTCGCCGCCCATCGCCTGGGTCAGACGGCTACTGATAGTCAGTCCCAGCCCGGTGCCGCCGCGTTTGCCGCTTACCTGCACAAATGGc                                                                    <  2:397262/140‑1 (MQ=255)
                                                                   gCGTAAACAAAAACAGCTTCCAACCTCCGGCGTGCTGGTGGCGCTCAGTTCGCCGCCCATCGCCTGGGTCAGACGGCTACTGATAGTCAGTCCCAGCCCGGTGCCGCCGCGTTTGCCGCTTACCTGCACAAATGGCTGGa                                                                <  1:430904/140‑1 (MQ=255)
                                                                            aaaaCAGCTTCCAACCTCCGGCGTGCTGGTGGCGCTCAGTTCGCCGCCCATCGCCTGGGTCAGACGGCTACTGATAGTCAGTCCCAGCCCGGTGCCGCCGCGTTTGCCGCTTACCTGCACAAATGGCTGGAAGATTTCTg                                                       >  1:102611/1‑140 (MQ=255)
                                                                             aaaCAGCTTCCAACCTCCGGCGTGCTGGTGGCGCTCAGTTCGCCGCCCATCGCCTGGGTCAGACGGCTACTGATAGTCAGTCCCAGCCCGGTGCCGCCGCGTTTGCCGCTTACCTGCACAAATGGCTGGAAGATTTCTGc                                                      >  2:425785/1‑140 (MQ=255)
                                                                                      ccAACCTCCGGCGTGCTGGTGGCGCTCAGTTCGCCGCCCATCGCCTGGGTCAGACGGCTACTGATAGTCAGTCCCAGCCCGGTGCCGCCGCGTTTGCCGCTTACCTGCACAAATGGCTGGAAGATTTCTGCCAGTTTcgc                                             <  1:73131/140‑1 (MQ=255)
                                                                                      ccAACCTCCGGCGTGCTGGTGGCGCTCAGTTCGCCGCCCATCGCCTGGGTCAGACGGCTACTGATAGTCAGTCCCAGCCCGGTGCCGCCGCGTTTGCCGCTTACCTGCACAAATGGCTGGAAGATTTCTGCCAGTTTcgc                                             <  1:75354/140‑1 (MQ=255)
                                                                                              cGGCGTGCTGGTGGCGCTCAGTTCGCCGCCCATCGCCTGGGTCAGACGGCTACTGATAGTCAGTCCCAGCCCGGTGCCGCCGCGTTTGCCGCTTACCTGCACAAATGGCTGGAAGATTTCTGCCAGTTTCGCGGGATCaa                                     >  2:183120/1‑140 (MQ=255)
                                                                                                         ccaggcTCAGTTCGCCACCCATCGGCTGGGTCAGCCGGCTACTGATAGTCAGTCCCAGCGCGGTGCCCACGCGGTTGCCGCTTACCTGCACAAATGGCTGGAAGATTTCTGCCAGTTTCGCGGGATCAATACCGCAGCCg                          <  2:811808/136‑1 (MQ=255)
                                                                                                                        cgccCATCGCCTGGGTCAGACGGCTACTGATAGTCAGTCCCAGCCCGGTGCCGCCGCGTTTGCCGCTTACCTGCACAAATGGCTGGAAGATTTCTGCCAGTTTCGCGGGATCAATACCGCAGCCGCTGTCTTCCAcctgt           >  1:811808/1‑136 (MQ=255)
                                                                                                                              tCGCCTGGGTCAGACGGCTACTGATAGTCAGTCCCAGCCCGGTGCCGCCGCGTTTGCCGCTTACCTGCACAAATGGCTGGAAGATTTCTGCCAGTTTCGCGGGATCAATACCGCAGCCGCTGTCTTCCACTTCGACCAGc     <  1:2280/140‑1 (MQ=255)
                                                                                                                                 ccTGGGTCAGACGGCTACTGATAGTCAGTCCCAGCCCGGTGCCGCCGCGTTTGCCGCTTACCTGCACAAATGGCTGGAAGATTTCTGCCAGTTTCGCGGGATCAATACCGCAGCCGCTGTCTTCCACTTCGACCAGCCAt  <  1:183120/140‑1 (MQ=255)
                                                                                                                                       |                                                                                                                                     
AACCGTCAAGACGCACCGCCTGATTGACTGTTTTTGGCACGGGTGCCGTGGCAACACGTAACGGCAAGCGTAAACAAAAACAGCTTCCAACCTCCGGCGTGCTGGTGGCGCTCAGTTCGCCGCCCATCGCCTGGGCCAGACGGCTACTGATAGTCAGTCCCAGCCCGGTGCCGCCGCGTTTGCCGCTTACCTGCACAAATGGCTGGAAGATTTCTGCCAGTTTCGCGGGATCAATACCGCAGCCGCTGTCTTCCACTTCGACCAGCCAT  >  NC_000913/1054135‑1054403

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