Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 4,651,193 C→T 8.3% R317Q (CGG→CAG)  PP_4114 ← hypothetical protein

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0029474,651,1930CT 8.3% 133.7 / 2.5 48R317Q (CGG→CAG) PP_4114hypothetical protein
Reads supporting (aligned to +/- strand):  ref base C (19/25);  new base T (2/2);  total (21/27)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 8.31e-01

GCCTTCCATACCAGCTCCTTGTCGCGCAGGGCATCAAGGGCCTTCTGTACGTTTGCCGTGTTGACGTTGGCTACCGCCTGGGCCTGCTCCAGCTTGCGGCTGACCTCAACCAGGGTTTGCTCGGTAAACGGGGAAAACGGCTGCCGGCCCAGCGTGCGCTCGGCCATGACCTCTAGCAC‑GGCCTGCTGCAACGGCGTAAGTT‑CGTTCCAAGCGCTTTCGTAGTCACTCCACACCCCAGCTTGATGGTTCAGCGCCCCATTGCGCAGCAGTTCACCCAGGTTGCCCGCCTCACC  >  NC_002947/4651049‑4651341
                                                                                                                                                |                                                                                                                                                      
gCCTTCCATACCAGCTCCTTGTCGCGCAGGGCATCAAGGGCCTTCTGTACGTTTGCCGTGTTGACGTTGGCTACCGCCTGGGCCTGCTCCAGCTTGCGGCTGACCTCAACCAGGGTTTGCTCGGTAAACGGGGAAAACGGCTGCCGGcc                                                                                                                                                    <  8:273436/149‑1 (MQ=255)
  cTTCCATACCAGCTCCTTGTCGCGCAGGGCATCAAGGGCCTTCTGTACGTTTGCCGTGTTGACGTTGGCTACCGCCTGGGCCTGCTCCAGCTTGCGGCTGACCTCAACCAGGGTTTGCTCGGTAAACGGGGAAAACGGCTGCCGGcc                                                                                                                                                    >  5:316568/1‑147 (MQ=255)
     ccATACCAGCTCCTTGTCGCGCAGGGCATCAAGGGCCTTCTGTACGTTTGCCGTGTTGACGTTGGCTACCGCCTGGGCCTGCTCCAGCTTGCGGCTGACCTCAACCAGGGTTTGCTCGGTAAACGGGGAAAACGGCTGCCGGCCCAGc                                                                                                                                                <  1:56651/148‑1 (MQ=255)
     ccATACCAGCTCCTTGTCGCGCAGGGCATCAAGGGCCTTCTGTACGTTTGCCGTGTTGACGTTGGCTACCGCCTGGGCCTGCTCCAGCTTGCGGCTGACCTCAACCAGGGTTTGCTCGGTAAACGGGGAAAACGGCTGCCGGCCCAGc                                                                                                                                                <  1:58394/148‑1 (MQ=255)
      cATACCAGCTCCTTGTCGCGCAGGGCATCAAGGGCCTTCTGTACGTTTGCCGTGTTGACGTTGGCTACCGCCTGGGCCTGCTCCAGCTTGCGGCTGACCTCAACCAGCGTTTGCTCGGTAAACGGGGAAAACGGCTGCCGGCCCAGCGt                                                                                                                                              >  5:281921/1‑149 (MQ=255)
                  ttGTCGCGCAGGGCATCAAGGGCCTTCTGTACGTTTGCCGTGTTGACGTTGGCTACCGCCTGGGCCTGCTCCAGCTTGCGGCTGACCTCAACCAGGGTTTGCTCGGTAAACGGGGAAAACGGCTGCCGGCCCAGCGTGCGCTCGGCCAt                                                                                                                                  >  5:249766/1‑149 (MQ=255)
                     tCGCGCAGGGCATCAAGGGCCTTCTGTACGTTTGCCGTGTTGACGTTGGCTACCGCCTGGGCCTGCTCCAGCTTGCGGCTGACCTCAACCAGGGTTTGCTCGGTAAACGGGGAAAACGGCTGCCGGCCCAGCGTGCGCTCGGCCATGAc                                                                                                                               <  6:316568/149‑1 (MQ=255)
                       gcgcAGGGCATCAAGGGCCTTCTGTACGTTTGCCGTGTTGACGGTGGCTACCGCCTGGGCCTGCTCCAGCTTGCGGCTGACCTCAACCAGGGTTTGCTCGGTAAACGGGGAAAACGGCTGCCGGCCCAGCGTGCGCTCGGCCATGACct                                                                                                                             <  1:240071/149‑1 (MQ=255)
                                aTCAAGGGCCTTCTGTACGTTTGCCGTGTTGACGTTGGCTACCGCCTGGGCCTGCTCCAGCTTGCGGCTGACCTCAACCAGGGTTTGCTCGGTAAACGGGGAAAACGGCTGCCGGCCCAGCGTGCGCTCGGCCATGACCTCTAGCAC‑gg                                                                                                                   >  1:148845/1‑149 (MQ=255)
                                              gTACGTTTGCCGTGTTGACGTTGGCTACCGCCTGGGCCTGCTCCAGCTTGCGGCTGACCTCAACCAGGGTTTGCTCGGTAAACGGGGAAAACGGCTGCCGGCCCAGCGTGCGCTCGGCCATGACCTCTAGCAC‑GGCCTGCTGCAACgg                                                                                                      <  3:165785/148‑1 (MQ=255)
                                              gTACGTTTGCCGTGTTGACGTTGGCTACAGCCTGCGCCTGCTACAGCTTGCGGCTGACCTCAACCAGGGTTAGCTCGGTAAACGGTGAAAACGCCTGCCGGCCCAGCGTGCGCTCGGCCATGACCTCTATCAC‑TCCCTCCTGGAACgg                                                                                                      >  4:165785/1‑148 (MQ=255)
                                               tACGTTTGCCGTGTTGACGTTGGCTACCGCCTGGGCCTGCTCCAGCTTGCGGCTGACCTCAACCAGGGTTTGCTCGGTAAACGGGGAAACCGGCTGCCGGCCCAGCGTGCGCTCGGCCATGACCTCTAg                                                                                                                         <  4:66376/129‑1 (MQ=255)
                                               tACGTTTGCCGTGTTGACGTTGGCTACCGCCTGGGCCTGCTCCAGCTTGCGGCTGACCTCAACCAGGGTTTGCTCGGTAAACGGGGAAACCGGCTGCCGGCCCAGCGTGCGCTCGGCCATGACCTCTAg                                                                                                                         >  3:66376/1‑129 (MQ=255)
                                               tACGTTTGCCGTGTTGACGTTGGCTACCGCCTGGGCCTGCTCCAGCTTGCGGCTGACCTCAACCAGGGTTTGCTCGGTAAACGGGGAAAACGGCTGCTGGCCCAGCGTGCGCTCGGCCATGACCTCTAGCAC‑GGCCTGCTGCAACGGc                                                                                                     >  5:70173/1‑148 (MQ=255)
                                               tACGTTTGCCGTGTTGACGTTGGCTACCGCCTGGGCCTGCTCCAGCTTGCGGCTGACCTCAACCAGGGTTTGCTCGGTAAACGGGGAAAACGGCTGCCGGCCCAGCGTGCGCTCGGCCATGACCTCTAGCAC‑GGCCTGCTGCAACg                                                                                                       <  5:249085/146‑1 (MQ=255)
                                               tACGTTTGCCGTGTTGACGTTGGCTACCGCCTGGGCCTGCTCCAGCTTGCGGCTGACCTCAACCAGGGTTTGCTCGGTAAACGGGGAAAACGGCTGCCGGCCCAGCGTGCGCTCGGCCATGACCTCTAGCAC‑GGCCTGCTGCAACGGc                                                                                                     <  2:18687/148‑1 (MQ=255)
                                               tACGTTTGCCGTGTTGACGTTGGCTACCGCCTGGGCCTGCTCCAGCTTGCGGCTGACCTCAACCAGGGTTTGATCGGTAAACGGGGAAAACGGCTGCCGGCCCACCGTGCGAACGGCCATGACCGCTAGCAC‑GCCCTGCTGCAACggg                                                                                                     >  7:204725/1‑147 (MQ=255)
                                                     tGCCGTGTTGACGTTGGCTACCGCCTGGGCCTGCTCCAGCTTGCGGCTGACCTCAACCAGGGTTTGCTCGGTAAACGGGGAAAACGGCTGCCGGCCCAGCGTGCGCTCGGCCATGACCTCTAGCAC‑GGCCTGCTGCAACGGCGTAAgt                                                                                               <  7:40119/148‑1 (MQ=255)
                                                            ctGACGTTGGCTAACGCCTGGGCCTGCTCCATCTTGCGGCTGAACTCAACCAGGGTTTGCTCGGTAAAAGGGGAAAACGGCTGCCGGCCCAGCGTGCGCTCGGCCATGACCTCTAGGACAGGCCTGCTGCAACGGCGTAAgtt‑cgttc                                                                                        <  1:187910/147‑1 (MQ=255)
                                                              gACGTTGGCTACCGCCTGGGCCTGCTCCAGCTTGCGGCTGACCTCAACCAGGGTTTGCTCGGTAAACGGGGAAAACGGCTGCCGGCCCAGCGTGCGCTCGGCCATGACCTCTAGCAC‑GGCCTACTGCAACGGCGTAAGTT‑CGTTCCaa                                                                                     <  6:249766/148‑1 (MQ=255)
                                                               aCGTTGGCTACCGCCTGGGCCTGCTCCAGCTTGCGGCTGACCTCAACCAGGGTTTGCTCGGTAAACGGGGAAAACGGCTGCCGGCCCAGCGTGCGCTCGGCCATGACCTCTAGCAC‑GGCCTGCTGCAACGGCGTAAgtt‑cgttc                                                                                        >  2:187910/1‑144 (MQ=255)
                                                                 gTTGGCTACCGCCTGGGCCTGCTCCAGCTTGCGGCTGACCTCAACCAGGGTTTGCTCGGTAAACGGGGAAAACGGCTGCCGGCCCAGCGTGCGCTCGGCCATGAc                                                                                                                               <  7:217211/105‑1 (MQ=255)
                                                                 gTTGGCTACCGCCTGGGCCTGCTCCAGCTTGCGGCTGACCTCAACCAGGGTTTGCTCGGTAAACGGGGAAAACGGCTGCCGGCCCAGCGTGCGCTCGGCCATGAc                                                                                                                               >  8:217211/1‑105 (MQ=255)
                                                                    ggCTACCGCCTGGGCCTGCTCCAGCTTGCGGCTGACCTCAACCAGGGTTTGCTCGGTAAACGGGGAAAACGGCTGCCGGCCCAGCGTGCGCTCGGCCATGACCTCTAGCAC‑GGCCTGCTGCAACGGCGTAAGTT‑CGTTCCAAGCGCtt                                                                               <  6:174404/148‑1 (MQ=255)
                                                                      cTACCGCCTGGGCCTGCTCCAGCTTGCGGCTGACCTCAATAAGGGTTTGCTCGGTAAACGGGGAAAACGGCTGCCGGCCCAGCGTGCGCTCGGCCATGACCTCTAGCAC‑GGCCTGCTGCAACGGCGTAAGTT‑CGTTCCAAGCGCTTTCg                                                                            <  1:110362/149‑1 (MQ=255)
                                                                       tACTGCCTGGGCCTGCTCCAGCTTGCGGCTGACCTCAACCAGGCTTTGCTCGGTAAACGGGGAAAACGGCTGCCGGCCCAGCGTGCGCTCGGCCATGACCTCTAGCAC‑GGCCTGCTGCAACGGCGTAAGTT‑CGTTCCCAGCGCTTTCg                                                                            >  1:148410/1‑148 (MQ=255)
                                                                                      cTCCAGCTTGCGGCTGACCTCAACCAGGGTTTGCTCGGTAAACGGGGAAAACGGCTGCCGGCCCAGCGTGCGCTCGGCCATGACCTCTAGCAC‑GGCCTGCTGCAACGGCGTAAGTTCCGTTCCAAGCGCTTTCGTAGTCACTCcaca                                                               <  1:24521/147‑1 (MQ=255)
                                                                                      cTCCAGCTTGCGGCTGACCTCAACCAGGGTTTGCTCGGTAAACGGGGAAAACGGCTGCCGGCCCAGCGTGCGCTCGGCCATGACCTCTAGCAC‑GGCCTGCTGCAACGGCGTAAGTT‑CGTTCCAAGCGCTTTCGTAGTCACTCCacac                                                              <  8:98720/147‑1 (MQ=255)
                                                                                       tCCAGCTTGCGGCTGACCTCAACCAGGGTTTGCTCGGTAAACGGGGAAAACGGCTGCTGGCCCAGCGTGCGCTCGGCCATGACCTCTAGCAC‑GGCCTGCTGCAACGGCGTAAGTT‑CGTTCCAAGCGCTTTCGTAGTCACTCCACAccc                                                            >  1:272610/1‑148 (MQ=255)
                                                                                         cAGCTTGCGGCTGACCTCAACCAGGGTTTGCTCGGTAAACGGGGAAAACGGCTGCCGGCCCAGCGTGCGCTCGGCCATGACCTCTAGCAC‑GGCCTGCTGCAACGGCGTAAGTC‑CGTTCCAAGCGCTTTCGTAGTCACTCCACACCCCAg                                                         >  2:138716/1‑149 (MQ=255)
                                                                                                 ggCTGACCTCAACCAGGGTTTGCTCGGTAAACGGGGAAAACGGCTGCCGGCCCAGCGTGCGCTCGGCCATGACCTCTAGCAC‑GGCCTGCTGCAACGGCGTAAGTT‑CGTTCCAAGCGCTTTCGTAGTCACTCCACACCCCAGCTTGATgg                                                 >  7:77570/1‑149 (MQ=255)
                                                                                                        cTCAACCAGGGTTTGCTCGGTAAACGGGGAAAACGGCTGCCGGCCCAGCGTGCGCTCGGCCATGACCTCTAGCAC‑GGCCTGCTGCAACGGCGTAAGTT‑CGTTCCAAGCGCTTTCGTAGTCACTCCACACCCCAGCTTGATGGTTCAgcg                                          >  1:241686/1‑149 (MQ=255)
                                                                                                        cTCAACCAGGGTTTGCTCGGT‑AACGGGGAAAACGGCTGCCGGCCCAGCGTGCGCTCGGCCATGACCTCTAGCAC‑GGCCTGCTGCAACGGCGTAAGTT‑CGTTCCAAGCGCTTTCGTAGTCACTCCACACCCCAGCTTGATGGTTCAgcgc                                         >  4:28297/1‑149 (MQ=255)
                                                                                                              cAGGGTTTGCTCGGTAAACGGGGAAAACGGCTGCCGGCCCAGCGTGCGCTCGGCCATGACCTCTAGCAC‑GGCCTGCTGCAACGGCGTAAGTT‑CGTTCCAAGCGCTTTCGTAGTCACTCCACACCCCAGCTTGATGGTTCAGCGCCCCAt                                    >  2:45280/1‑149 (MQ=255)
                                                                                                                 ggTTTGCTCGGTAAACGGGGAAAACGGCTGCCGGCCCAGCGTGCGCTCGGCCATGACCTCTAGCAC‑GGCCTGCTGCAACGGCGTAAGTT‑CGTTCCAAGCGCTTTCGTAGTCACTCCACACCCCAGCTTGATGGTTCAGCGCCCCATTgc                                 <  4:318762/149‑1 (MQ=255)
                                                                                                                  gTTTGCTCGGTAAACGGGGAAAACGGCTGCTGGCCCAGCGTGCGCTCGGCCATGACCTCTAGCAC‑GGCCTGCTGCAACGGCGTAAGTT‑CGTTCCAAGCGCTTTCGTAGTCACTCCACACCCCAGCTTGATGGTTCAGCGCCCCATTgcg                                <  2:272610/149‑1 (MQ=255)
                                                                                                                     tGCTCGGTAAACGGGGAAAACGGCTGCCGGCCCAGCGTGCGCTCGGCCATGACCTCTAGCAC‑GGCCTGCTGCAACGGCGTAAGTT‑CGTTCCAAGCGCTTTCGTAGTCACTCCACTCCCCAGCTTGATGGTTCAGCGCCCCATTGCgca                              <  1:45280/148‑1 (MQ=255)
                                                                                                                      gCTCGGTAAACGGGGAAAACGGCTGCCGGCCCAGCGTGCGCTCGGCCATGACCTCTAGCAC‑GGCCTGCTGCAACGGCGTAAGTT‑CGTTCCAAGCGCTTTCGTAGTCACTCCACACCCCAGCTTGATGGTTCAGCGCCCCATTGCgcag                             >  3:276375/1‑148 (MQ=255)
                                                                                                                       cTCGGTAAACGGGGAAAACGGCTGCCGGCCCAGCGTGCGCTCGGCCATGACCTCTAGCAC‑GGCCTGCTGCAACGGCGTAAGTT‑CGTTCCAAGCGCTTTCGTAGTCACTCCACACCCCAGCTTGATGGTTCAGCGCCCCATTGCgcagca                           <  8:77570/149‑1 (MQ=255)
                                                                                                                       cTCGGTAAACGGGGAAAACGGCTGCCGGCCCAGCGTGCGCTCGGCCATGACCTCTAGCAC‑GGCCTGCTGCAACGGCGTAAGTT‑CGTTCCAAGCGCTTTCGTAGTCACTCCACACCCCAGCTTGATGGTTCAGCGCCCCATTGCgcagca                           <  2:173146/149‑1 (MQ=255)
                                                                                                                          ggTAAACGGGGAAAACGGCTGCTGGCCCAGCGTGCGCTCGGCCATGACCTCTAGCAC‑GGCCTGCTGCAACGGCGTAAGTT‑CGTTCCAAGCGCTTTCGTAGTCACTCCACACCCCAGCTTGATGGTTCAGCGCCCGATTGCGCAGCAGtt                        <  6:70173/149‑1 (MQ=255)
                                                                                                                          ggTAAACGGGGAAAACGGCTGCCGGCCCAGCGTGCGCTCGGCCATGACCTCTAGCAC‑GGCCTGCTGCAACGGCGTAAGTT‑CGTTCCAAGCGCTTTCGTAGTCACTCCACACCCCAGCTTGATGGTTCAGCGCCCCATTGCgcagca                           >  6:86214/1‑146 (MQ=255)
                                                                                                                          ggTAAACGGGGAAAACGGCTGCCGGCCCAGCGTGCGCTCGGCCATGACCTCTAGCAC‑GGCCTGCTGCAACGGCGTAAGTT‑CGTTCCAAGCGCTTTCGTAGTCACTCCACACCCCAGCTTGATGGTTCAGCGCCCCATTGCgcagca                           <  5:86214/146‑1 (MQ=255)
                                                                                                                          ggTAAACGGGGAAAACGGCTGCCGGCCCAGCGTGCGCTCGGCCATGACCTCTAGCAC‑GGCCTGCTGCAACGGCGTAAGTC‑CGTTCCAAGCGCTTTCGTAGTCACTCCACACCCCAGCTTGATGGTTCAGCGCCCCATTGCGCAGCAGt                         <  1:138716/148‑1 (MQ=255)
                                                                                                                               aCGGGGAAAACGGCTGCCGGCCCAGCGTGCGCTCGGCCATGACCTCTAGCAC‑GGCCTGCTGCAACGGCGTAAGTT‑CGTTCCAAGCGCTTTCGTAGTCACTCCACACCCCAGCTTGa                                                    >  7:71421/1‑116 (MQ=255)
                                                                                                                               aCGGGGAAAACGGCTGCCGGCCCAGCGTGCGCTCGGCCATGACCTCTAGCAC‑GGCCTGCTGCAACGGCGTAAGTT‑CGTTCCAAGCGCTTTCGTAGTCACTCCACACCCCAGCTTGa                                                    <  8:71421/116‑1 (MQ=255)
                                                                                                                                           gCTGCCGGCCCAGCGTGCGCTCGGCCATGACCTCTAGCAC‑GGCCTGCTGCAACGGCGTAAGTT‑CGTTCCAAGCGCTTTCGTAGTCACTCCACACCCCAGCTTGATGGTTCAGCGCCCCAtt                                   >  5:173933/1‑121 (MQ=255)
                                                                                                                                           gCTGCCGGCCCAGCGTGCGCTCGGCCATGACCTCTAGCAC‑GGCCTGCTGCAACGGCGTAAGTT‑CGTTCCAAGCGCTTTCGTAGTCACTCCACACCCCAGCTTGATGGTTCAGCGCCCCAtt                                   <  6:173933/121‑1 (MQ=255)
                                                                                                                                               ccGGCCCAGCGTGCGCTCGGCCATGACCTCTAGCAC‑GGCCTGCTGCAACGGCGTAAGTT‑CGTTCCAAGCGCTTTCGTAGTCACTCCACACCCCTGCTTGATGGTTCAGCGCCCCATTGCGCAGCAGTTCACCCAGGTTGCCCGCCTCAc   <  2:148410/149‑1 (MQ=255)
                                                                                                                                                cGGCCCAGCGTGCGCTCGGCCATGACCTCTAGCAC‑GGCCTGCTGCAACGGCGTAAGTT‑CGTTCCAAGCGCTTTCGTAGTCACTCCACACCCCAGCTTGATGGTTCAGCGCCCCATTGCGCAGCAGTTCACCCAGGTTGCCCGCCTCAcc  <  6:281921/149‑1 (MQ=255)
                                                                                                                                                |                                                                                                                                                      
GCCTTCCATACCAGCTCCTTGTCGCGCAGGGCATCAAGGGCCTTCTGTACGTTTGCCGTGTTGACGTTGGCTACCGCCTGGGCCTGCTCCAGCTTGCGGCTGACCTCAACCAGGGTTTGCTCGGTAAACGGGGAAAACGGCTGCCGGCCCAGCGTGCGCTCGGCCATGACCTCTAGCAC‑GGCCTGCTGCAACGGCGTAAGTT‑CGTTCCAAGCGCTTTCGTAGTCACTCCACACCCCAGCTTGATGGTTCAGCGCCCCATTGCGCAGCAGTTCACCCAGGTTGCCCGCCTCACC  >  NC_002947/4651049‑4651341

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