Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 1,488,066 T→A 17.6% Y86F (TAC→TTC)  degS ← serine endoprotease DegS

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0029471,488,0660TA17.6% 72.9 / 10.2 34Y86F (TAC→TTC) degSserine endoprotease DegS
Reads supporting (aligned to +/- strand):  ref base T (16/12);  new base A (3/3);  total (19/15)
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.34e-01

AGCAGGTAGCCCTCTGGGCTCATGATTACCGCCGAACCCAGGCTTGATTCCCAGCGGCGCTGTTTGGGCAGGTTATCGCCGAAGAAGCGGCGGAACTGCGGGTCTTCGAACAGCGGGTGGGCGCTTTTGTTCACCACCTTGGTGGTGTACAGGTTGGCCACCGCCGGGGCGGCCAGGGTCACGGCGTCGGCGTACGACACCGGGCCTTGCATGATCCGGGTAGTCTGTGGCGCTTGCTGCAGGTTGACGTCCTGGCTGGGCAGCCCGACCCATTGCGGAAAACGCT  >  NC_002947/1487919‑1488204
                                                                                                                                                   |                                                                                                                                          
agcagGTAGCCCTCTGGGCTCATGATTACCGCCGAACCCAGGCTTGATTCCCAGCGGCGCTGTTTGGGCAGGTTATCGCCGAAGAAGCGGCGGAACTGCGGGTCTTCGAACAGCGGGTGGGCGCTTTTGTTCACCACCTTGGTGgtgt                                                                                                                                            <  4:138791/148‑1 (MQ=255)
 gcagGTAGCCCTCTGGGCTCATGATTACCGCCGAACCCAGGCTTGATTCCCAGCGGCGCTGTTTGGGCAGGTTATCGCCGAAGAAGCGGCGGAACTGCGGGTCTTCGAACAGCGGGTGGGCGCTTTTGTACACCACCAAGGTGGTGAAc                                                                                                                                          <  6:3072/149‑1 (MQ=255)
                     aTGATTACCGCCGAACCCAGGCTTGATTCCCAGCGGCGCTGTTTGGGCAGGTTATCGCCGAAGAAGCGGCGGAACTGCGGGTCTTCGAACAGCGGGTGGGCGCTTTTGTTCACCACCTTGGTGGTGTACAGGTTGGCCACCGCCggg                                                                                                                        <  6:205677/147‑1 (MQ=255)
                             cgccgAACCCAGGCTTGATTCCCAGCGGCGCTGTTTGGGCAGGTTATCGCCGAAGAAGCGGCGGAACTGCGGGTCTTCGAACAGCGGGTGGGCGCTTTTGTTCACCACCTTGGTGGTGTACAGGTTGGCCACCGCCGGGGCGGCCAggg                                                                                                              >  3:11277/1‑149 (MQ=255)
                               ccgAACCCAGGCTTGATTCCCAGCGGCGCTGTTTGGGCAGGTTAGCGCCGAAGCAGCGGGGGAACGGCGGGGGGTCGAGCCGCGGGGGGGGGGTTTTGTTCACCCCCGGGGGGGGGTACGGGTTGGCCACCGCGGGGGGGGCCAGGGTc                                                                                                            >  8:73590/1‑149 (MQ=255)
                                       aGGCTTGATTCCCAGCGGCGCTGTTTGGGCAGGTTATCGCCGAAGAAGCGGCGGAACTGCGGGTCTTCGAACAGCGGGTGGGCGCTTTTGTTCACCACCTTGGTGGTGTACAGGTTGGCCACCGCCGGGGCGGCCAGGGTCAcggcgtc                                                                                                    >  6:94755/1‑149 (MQ=255)
                                         gCTTGATTCCCAGCGGCGCTGTTTGGGCAGGTTATCGCCGAAGAAGCGGCGGAACTGCGGGTCTTCGAACAGCGGGTGGGCGCTTTTGTTCACCACCTTGGTGGTGTACAGGTTGGCCACCGCCGGGGCGGCCAGGGTCAcggcgtcgg                                                                                                  <  8:2939/149‑1 (MQ=255)
                                                   cAGCGGCGCTGTTTGGGCAGGTTATCGCCGAAGAAGCGGCGGAACTGCGGGTCTTCGAACAGCGGGTGGGCGCTTTTGTTCACCACCTTGGTGGTGTACAGGTTGGCCACAGCCGGGGCGGCCAGGGTCACGGCGTCGGCGTACGacac                                                                                        >  3:193315/1‑149 (MQ=255)
                                                       ggcgCTGTTTGGGCAGGTTATCGCCGAAGAAGCGGCGGAACTGCGGGTCTTCGAACAGCGGGTGGGCGCTTTTGTTCACCACCTTGGTGGTGTACAGGTTGGCCACCGCCGGGGCGGCCAGGGTCACGGCGTCGGCGTACGACACCggg                                                                                    <  4:193315/149‑1 (MQ=255)
                                                                 gggCAGGTTATCGCCGAAGAAGCGGCGGAACTGCGGGTCTTCGAACAGCGGGTGGGCGCTTTTGTTCACCACCTTGGTGGTGTACAGGTTGGCCACCGCCGGGGCGGCCAGGGTCACGGCGTCGGCGTACGACACCGGGCCTTGCATGa                                                                          >  3:236063/1‑149 (MQ=255)
                                                                            cGCCGAAGAAGCGGCGGAACTGCGGGTCTTCGAACAGCGGGTGGGCGCTTTTGTTCACCACCTTGGTGGTGTACAGGTTGGCCACCGCCGGGGCGGCCAGGGTCACGGCGTCGGCGTACGACACCGGGCCTTGCATGATCCGGGTAGTc                                                               >  8:1052/1‑149 (MQ=255)
                                                                             gCCGAAGAAGCGGCGGAACTGCGGGTCTTCGAACAGCGGGTGGGCGCTTTTGTTCACCACCTTGGTGGTGTACAGGTTGGCCACCGCCGGGGCGGCCAGGGTCACGGCGTCGGCGTACGACACCGGGCCTTGCATGATCCGGGTAGTCt                                                              >  8:200623/1‑149 (MQ=255)
                                                                             gCCGAAGAAGCGGCGGAACTGCGGGTCTTCGAACAGCGGGTGGGCGCTTTTGTTCACCACCTTGGTGGTGTACAGGTTGGCCACCGCCGGGGCGGCCAGGGTCACGGCGTCGGCGTACGACACCGGGCCTCGCATGATCCGGGTAGTc                                                               <  8:248269/148‑1 (MQ=255)
                                                                              ccGAAGAAGCGGCGGAACTGCGGGTCTTCGAACAGCGGGTGGGCGCTTTTGTTCACCACCTTGGTGGTGTACa                                                                                                                                         <  6:298188/73‑1 (MQ=255)
                                                                              ccGAAGAAGCGGCGGAACTGCGGGTCTTCGAACAGCGGGTGGGCGCTTTTGTTCACCACCTTGGTGGTGTACa                                                                                                                                         >  5:298188/1‑73 (MQ=255)
                                                                              ccGAAGAAGCGGCGGAACTGCGGGTCTTCGAACAGCGGGTGGGCGCTTTTGTTCACCACCTTGGTGGTGTACAGGTTGGCCACCGCCGGGGCGGCCAGGGTCACGGCGTCGGCGTACGACACCGGGCCTTGCATGATCCGGGTAGTCtg                                                             <  7:200623/149‑1 (MQ=255)
                                                                                        ggcggAACTGCGGGTCTTCGAACAGCGGGTGGGCGCTTTTGTTCACCACCTTGGTGGTGTACAGGTTGGCCACCGCCGGGGCGGCCAGGGTCACGGCGTCGGCGTACGACACCGGGCCTTGCATGATCCGGGTAGTCTGTGGCGCTtgc                                                   <  4:236063/149‑1 (MQ=255)
                                                                                                  cGGGTCTTCGAACAGCGGGTGGGCGCTTTTGTTCACCACCTTGGTGGTGTACAGGTTGGCCACCGCCGGGGCGGCCAGGGTCACGGCGTCGGCGTACGACACCGGGCCTTGCATGATCCGGGTAGTCTGTGGCGCTTGCTGCAGGtt                                           <  7:1052/147‑1 (MQ=255)
                                                                                                   gggTCTTCGAACAGCGGGTGGGCGCTTTTGTTCACCACCTTGGTGGTGTACAGGTTGGCCACCGCCGGGGCGGCCAGGGTCACGGCGTCGGCGTACGACACCGGGCCTTGCATGATCCGGGTAGTCTGTGGCGCTTGCTGCAGGTTGAc                                        <  7:64846/149‑1 (MQ=255)
                                                                                                     gTCTTCGAACAGCGGGTGGGCGCTTTTGTTCACCACCTTGGTGGTGTACAGGTTGGCCACCGCCGGGGCGGCCAGGGTCACGGCGTCGGCGTACGACACCGGGCCTTGCATGATCCGGGTAGt                                                                <  4:85423/123‑1 (MQ=255)
                                                                                                     gTCTTCGAACAGCGGGTGGGCGCTTTTGTTCACCACCTTGGTGGTGTACAGGTTGGCCACAGCCGGGGCGGCCAGGGTCACGGCGTCGGCGTACGACACCGGGCCTTGCATGATCCGGGTAGt                                                                >  3:85423/1‑123 (MQ=255)
                                                                                                        ttCGAACAGCGGGTGGGCGCTTTTGTTCACCACCTTGGTGGTGTACAGGTTGGCCACCGCCGGGGCGGCCAGGGTCACGGCGTCGGCGTACGACACCGGGCCTTGCATGATCCGGGTAGTCTGTGGCGCTTGCTGCAGGTTGACGTCct                                   >  6:235200/1‑149 (MQ=255)
                                                                                                            aaCAGCGGGTGGGCGCTTTTGTTCACCACCTTGGTGGTGTACAGGTTGGCCACCGCCGGGGCGGCCAGGGTCACGGCGTCGGCGTACGACACCGGGCCTTGCATGATCCGGGTAGTCTGTGGCGCTTGCTGCAgg                                             <  3:78978/135‑1 (MQ=255)
                                                                                                            aaCAGCGGGTGGGCGCTTTTGTTCACCACCTTGGTGGTGTACAGGTTGGCCACCGCCGGGGCGGCCAGGGTCACGGCGTCGGCGTACGACACCGGGCCTTGCATGATCCGGGTAGTCTGTGGCGCTTGCTGCAgg                                             >  4:78978/1‑135 (MQ=255)
                                                                                                            aaCAGCGGGTGGGCGCTTTTGTTCACCACCTTGGTGGTGTACAGGTTGGCCACCGCCGGGGCGGCCAGGGTCACGGCGTCGGCGTACGACACCGGGCCTTGCATGATCCGGGTAGTCTGTGGCGCTTGCTGCAGGTTGACGTCctggct                               >  6:32105/1‑149 (MQ=255)
                                                                                                             aCAGCGGGTGGGCGCTTTTGTTCACCACCTTGGTGGTGTACAGGTTGGCCACCGCCGGGGCGGCCAGGGTCACGGCGTCGGCGTACGACACCGGGCCTTGCATGAGCCGGGTAGTCTGTGGCGCTTGCTGCTGGTTGACGTCctggctg                              >  7:273802/1‑149 (MQ=255)
                                                                                                              cAGCGGGTGGGCGCTTTTGTTCACCACCTTGGTGGTGTACAGGTTGGCCACCGCCGGGGCGGCCAGGGTCACGGCGTCGGCGTACGACACCGGGCCTTGCATGATCCGGGTAGTCTGTGGCGCTTGCTGCAGGTTGACGTCctggctgg                             <  3:188230/149‑1 (MQ=255)
                                                                                                                         ccaacctGTACACCACCAAGGTGGTGAACAGGTTGGCCACCGCCGGGGCGGCCAGGGTCACGGCGTCGGCGTACGACACCGGGCCTTGCATGATCCGGGTAGTCTGTGGCGCTTGCTGCAGGTTGACGTCCTGGCTggg                            >  6:584/7‑139 (MQ=255)
                                                                                                                         ccaacctGTACACCACCAAGGTGGTGAACAGGTTGGCCACCGCCGGGGCGGCCAGGGTCACGGCGTCGGCGTACGACACCGGGCCTTGCATGATCCGGGTAGTCTGTGGCGCTTGCTGCAGGTTGACGTCCTGGCTggg                            <  5:584/133‑1 (MQ=255)
                                                                                                                            ttttGTTCACCACCTTGGTGGTGTACAGGTTGGCCACCGCCGGGGTGGCCAGGGTCACGGCGTCGGCGTACGACACCGGGCCTTGCATGATCCGGGTAGTCTGTGGCGCTTGCTGCAGGTTGACGTCCTGGCTGGGCAGCCCGACCCAt               >  3:121483/1‑149 (MQ=255)
                                                                                                                             cctGTACACCACCAAGGTGGTGAACAGGTTGGCCACCGCCGGGGCGGCCAGGGTCACGGCGTCGGTGTACGACACCGGGCCTTGCATGATCCGGGTAGTCTGTGGCGCTTGCTGCAGGTTGACGTCCTGGCTGGGCAGCCCGACCCAtt              >  3:70992/3‑149 (MQ=255)
                                                                                                                              ctGTACACCACCAAGGTGGTGAACAGGTTGGCCACCGCCGGGGCGGCCAGGGTCACGGCGTCGGCGTACGACACCGGGCCTTGCATGATCCGGGTAGTCTGTCGCGCTTGCTCCAGGTTGACGTCCTGGCTGGGCAGCCCGACCCATTg             >  8:33538/2‑149 (MQ=255)
                                                                                                                                 ttCACCACCTTGGTGGTGTACAGGTTGGCCACCGCCGGGGCGGCCAGGGTCACGGCGTCGGCGTACGACACCGGGCCTTGCATGATCCGGGTAGTCTGTGGCGCTTGCTGCAGGTTGACGTCCTGGCTGGGCAGCCCGACCCATTGCg           >  5:64481/1‑148 (MQ=255)
                                                                                                                                      caccTTGGTGGTGTACAGGTTGGCCACCGCCGGGGCGGCCAGGGTCACGGCGTCGGCGTACGACACCGGGCCTTGCATGATCCGGGTAGTCTGTGGCGCTTGCTGCAGGTTGACGTCCTGGCTGGGCAGCCCGACCCATTGCGGAAAAc     >  1:48344/1‑149 (MQ=255)
                                                                                                                                         caaggtggtGAACAGGTTGGCCACCGCCGGGGCGGCCAGGGTCACGGCGTCGGCGTACGACACCGGGCCTTGCATGATCCGGGTAGTCTGTGGCGCTTGCTGCAGGTTGACGTCCTGGCTGGGCAGCCCGACCCATTGCGGAAAACGCt  <  7:33538/146‑1 (MQ=255)
                                                                                                                                                   |                                                                                                                                          
AGCAGGTAGCCCTCTGGGCTCATGATTACCGCCGAACCCAGGCTTGATTCCCAGCGGCGCTGTTTGGGCAGGTTATCGCCGAAGAAGCGGCGGAACTGCGGGTCTTCGAACAGCGGGTGGGCGCTTTTGTTCACCACCTTGGTGGTGTACAGGTTGGCCACCGCCGGGGCGGCCAGGGTCACGGCGTCGGCGTACGACACCGGGCCTTGCATGATCCGGGTAGTCTGTGGCGCTTGCTGCAGGTTGACGTCCTGGCTGGGCAGCCCGACCCATTGCGGAAAACGCT  >  NC_002947/1487919‑1488204

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: