Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 3,630,365 T→C 9.2% W315R (TGG→CGG) ‡ PP_3200 → hypothetical protein

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0029473,630,3650TC 9.2% 148.4 / 5.1 54W315R (TGG→CGG) ‡PP_3200hypothetical protein
Reads supporting (aligned to +/- strand):  ref base T (26/23);  new base C (3/2);  total (29/25)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 7.36e-01

ATTTCTCCACAGACGACTTCTTTGGCGATGGCGCCGAGCGCATGTTCGTCGGCGCGCCGGTGTTCCCCAGTGGCCAACCTCTGGCGTTCTTCCACGGCAACGACAAGCAAGCCAAGGACTCAGGGCAGGGCGGCATGCAGGTGCGCTGGCGCACCGAGACGGTCGACTGGGGCCTGTACGCCATCCGTTTCCACGACAAGTCGCCGCAACTCAATGTACGCCCGAACTTCGCCAACCTCGACCCGGTCAGCGGCAAGGCCGGCGAGTTCTACTGGGTGTATCCCGAGG  >  NC_002947/3630219‑3630506
                                                                                                                                                  |                                                                                                                                             
aTTTCTCCACAGACGACTTCTTTGGCGATGGCGCCGAGTGCATGTTCGTCGGCGCGCCGGTGTTCCCCAGTGGCCAACCGCTGGCGTTCTTCCACGGCAACGACAAGCAAGCCAAGGACTCAGGGCAGGGCGGCATGCAGGTGCGCTGgc                                                                                                                                            <  1:239648/150‑1 (MQ=255)
      ccACAGACGACTTCTTTGGCGATGGCGCCGAGCGCATGTTCGTCGGCGCGCCGGTGTTCCCCAGTGGCCAACCTCTGGCGTTCTTCCACGGCAACGACAAGCAAGCCAAGGACTCAGGGCAGGGCGGCATGCAGGTGCGCTGGCGCAcct                                                                                                                                      <  2:669000/150‑2 (MQ=255)
           gacgacTTCTTTGGCGATGGCGCCGAGCGCATGTTCGTCGGCGCGCCGGTGTTCCCCAGTGGCCAACCTCTGGCGTTCTTCCACGGCAACGACAAGCAAGCCAAGGACTCAGGGCAGGGCGGCATGCAGGTGCGCTGGCGCACCGAGACg                                                                                                                                 <  1:981923/150‑1 (MQ=255)
                 ttcttTGGCGATGGCGCCGAGCGCATGTTCGTCGGCGCGCCGGTGTTCCCCAGTGGCCAACCTCTGGCGTTCTTCCACGGCAACGACAAGCAAGCCAAGGACTCAGGGCAGGGCGGCATGCAGGTGCGCTGGCGCACCGAGACGGTCGAc                                                                                                                           >  1:1225390/1‑150 (MQ=255)
                          gATGGCGCCGAGCGCATGTTCGTCGGCGCGCCGGTGTTCCCCAGTGGCCAACCTCTGGCGTTCTTCCACGGCAACGACAAGCAAGCCAAGGACTCAGGGCAGGGCGGCATGCAGGTGCGCTGGCGCACCGAGACGGTCGACTGGGGCCTg                                                                                                                  >  1:1065367/1‑150 (MQ=255)
                           aTGGCGCCGAGCGCATGTTCGTCGGCGCGCCGGTGTTCCCCAGTGGCCAACCTCTGGCGTTCTTCCACGGCAACGACAAGCAAGCCAAGGACTCAGGGCAGGGCGGCATGCAGGTGCGCTGGCGCACCGAGACGGTCGACTGGGGCCTGt                                                                                                                 >  2:991138/1‑150 (MQ=255)
                              gcgcCGAGCGCATGTTCGTCGGCGCGCCGGTGTTCCCCAGTGGCCAACCTCTGGCGTTCTTCCACGGCAACGACAAGCAAGCCAAGGACTCAGGGCAGGGCGGCATGCAGGTGCGCTGGCGCACCGAGAAGGTCGACTGGGGCCTGTACg                                                                                                              >  2:1227908/1‑150 (MQ=255)
                                gcCGAGCGCATGTTCGTCGGCGCGCCGGTGTTCCCCAGTGGCCAACCTCTGGCGTTCTTCCACGGCAACGACAAGCAAGCCAAGGACTCAGGGCAGGGCGGCATGCAGGTGCGCTGGCGCACCGAGACGGTCGACTGGGGCCTGTACGcc                                                                                                            <  1:121892/150‑1 (MQ=255)
                                gcCGAGCGCATGTTCGTCGGCGCGCCGGTGTTCCCCAGTGGCCAACCTCTGGCGTTCTTCCACGGCAACGACAAGCAAGCCAAGGACTCAGGGCAGGGCGGCATGCAGGTGCGCCAGCGCAcctgcat                                                                                                                                  <  1:836753/128‑6 (MQ=255)
                                gcCGAGCGCATGTTCGTCGGCGCGCCGGTGTTCCCCAGTGGCCAACCTCTGGCGTGCTTCCACGGCAACGACAAGCAAGCCAAGGACTCAGGGCAGGGCGGCATGCAGGTGCGCCAGCGCAcctgcat                                                                                                                                  >  2:836753/1‑123 (MQ=255)
                                     gcgcATGTTCGTCGGCGCGCCGGTGTTCCCCAGTGGCCAACCTCTGGCGTTCTTCCACGGCAACGACAAGCAAGCCAAGGACTCAGGGCAGGGCGGCATGCAGGTGCGCTGGCGCACCGAGACGGTCGACTGGGGCCTGTACGCCATCCg                                                                                                       >  2:1330040/1‑150 (MQ=255)
                                       gcATGTTCGTCGGCGCGCCGGTGTTCCCCAGTGGCCAACCTCTGGCGTTCTTCCACGGCAACGACAAGCAAGCCAAGGACTCAGGGCAGGGCGGCATGCAGGTGCGCTGGCGCACCGAGACGGTCGACTGGGGCCTGTACGCCATCCGtt                                                                                                     <  1:652308/150‑1 (MQ=255)
                                         aTGTTCGTCGGCGCGCCGGTGTTCCCCAGTGGCCAACCTCTGGCGTTCTTCCACGGCAACGACAAGCAAGCCAAGGACTCAGGGCAGGGCGGCATGCAGGTGCGCTGGTGCACCGAGACGGTCGACTGGGGCCTGTACGCCATCCGTTTc                                                                                                   <  1:329693/150‑1 (MQ=255)
                                         aTGTTCGTCGGCGCGCCGGTGTTCCCCAGTGGCCAACCTCTGGCGTTCTTCCACGGCAACGACAAGCAAGCCAAGGACTCAGGGCAGGGCGGCATGCAGGTGCGCTGGCGCACCGAGACGGTCGACTgggg                                                                                                                      <  2:719259/131‑1 (MQ=255)
                                         aTGTTCGTCGGCGCGCCGGTGTTCCCCAGTGGCCAACCTCTGGCGTTCTTCCACGGCAACGACAAGCAAGCCAAGGACTCAGGGCAGGGCGGCATGCAGGTGCGCTGGCGCACCGAGACGGTCGACTgggg                                                                                                                      >  1:719259/1‑131 (MQ=255)
                                              cgtcgGCGCGCCGGTGTTCCCCAGTGGCCAACCTCTGGCGTTCTTCCACGGCAACGACAAGCAAGCCAAGGACTCAGGGCAGGGCGGCATGCAGGTGCGCTGGCGCACCGAGACGGTCGACTGGGGCCTGTACGCCATCCGTTTCCAcca                                                                                              <  2:1065367/150‑3 (MQ=255)
                                                tcgGCGCGCCGGTGTTCCCCAGTGGCCAACCTCTGGCGTTCTTCCACGGCAACGACAAGCAAGCCAAGGACTCAGGGCAGGGCGGCATGCAGGTGCGCTGGCGCACCGAGACGGTCGACTGGGGCCTGTACGCCATCCGTTTCCACGACa                                                                                            >  1:191593/1‑150 (MQ=255)
                                                    cgcgcCGGTGTTCCCCAGTGGCCAACCTCTGGCGTTCTTCCACGGCAACGACAAGCAAGCCAAGGACTCAGGGCAGGGCGGCATGCAGGTGCGCTGGCGCACCGAGACGGTCGACTGGGGCCTGTACGCCATCCGTTTCCACGACAAGTc                                                                                        <  1:1330040/150‑1 (MQ=255)
                                                             gtTCCCCAGTGGCCAACCTCTGGCGTTCTTCCACGGCAACGACAAGCAAGCCAAGGACTCAGGGCAGGGCGGCATGCAGGTGCGCTGGCGCACCGAGACGGTCGACTGGGGCCTGTACGCCATCCGTTTCCACGACAAGTCGCCGCAACt                                                                               <  2:1167481/150‑1 (MQ=255)
                                                                ccccAGTGGCCAACCTCTGGCGTTCTTCCACGGCAACGACAAGCAAGCCAAGGACTCAGGGCAGGGCGGCATGCAGGTGCGCTGGCGCACCGAGACGGTCGACTGGGGCCTGTACGCCa                                                                                                           >  2:1064259/1‑119 (MQ=255)
                                                                ccccAGTGGCCAACCTCTGGCGTTCTTCCACGGCAACGACAAGCAAGCCAAGGACTCAGGGCAGGGCGGCATGCAGGTGCGCTGGCGCACCGAGACGGTCGACTGGGGCCTGTACGCCa                                                                                                           <  1:1064259/119‑1 (MQ=255)
                                                                 accAGTGGCCAACCTCTGGCGTTCTTCCACGGCAACGACAAGCAAGCCAAGGACTCAGGGCAGGGCGGCATGCAGGTGCGCTGGCGCACCGAGACGGTCGACTGGGGCCTGTACGCCATCCGTTTCCACGACAAGTCGCCGCAACTCAAt                                                                           >  1:275455/2‑150 (MQ=255)
                                                                           aaCCTCTGGCGTTCTTCCACGGCAACGACAAGCAAGCCAAGGACTCAGGGCAGGGCGGCATGCAGGTGCGCTGGCGCACCGAGACGGTCGACTGGGGCCTGTACGCCATCCGTTTCCACGACAAGTCGCCGCAACTCAATGTACGCCCGa                                                                 >  2:816499/1‑150 (MQ=255)
                                                                             ccTCTGGCGTTCTTCCACGGCAACGACAAGCAAGCCAAGGACTCAGGGCAGGGCGGCATGCAGGTGCGCTGGCGCACCGAGACGGTCGACTGGGGCCTGTACGCCATCCGTTTCCACGACAAGTCGCCGCAACTCAATGTACGCCCGAAc                                                               >  2:1177412/1‑150 (MQ=255)
                                                                               tctGGCGTTCTTCCACGGCAACGACAAGCAAGCCAAGGACTCAGGGCAGGGCGGCATGCAGGTGCGCTGGCGCACCGAGACGGTCGACTGGGGCCTGTACGCCATCCGTTTCCACGACAAGTCGCCGCAACTCAATGTACGCCCGAACtt                                                             >  2:404127/1‑150 (MQ=255)
                                                                                 tGGCGTTCTTCCACGGCAACGACAAGCAAGCCAAGGACTCAGGGCAGGGCGGCATGCAGGTGCGCTGGCGCACCGAGACGGTCGACTGGGGCCTGTACGCCATCCGTTTCCACGACAAGTCGCCGCAACTCAATGTACGCCCGAACTTCg                                                           <  2:1169682/150‑1 (MQ=255)
                                                                                         ttcCACGGCAACGACAAGCAAGCCAAGGACTCAGGGCAGGGCGGCATGCAGGTGCGCTGGCGCACCGAGACTGTAGACGGTGGCCTGTACGCCATCCGTTTCCACGACAAGTCGCCGCAACTCAATGTACGCCCGAACTTCGCCAACCTc                                                   >  2:1118666/1‑150 (MQ=255)
                                                                                         ttcCACGGCAACGACAAGCAAGCCAAGGACTCAGGGCAGGGCGGCATGCAGGTGCGCTGGCGCACCGAGACGGTCGACTGGGGCCTGTACGCCATCCGTTTCCACGACAAGTCGCCGCAACTCAATGTACGCCCGAACTTCGCCAAcgcc                                                   >  1:1266674/1‑147 (MQ=255)
                                                                                                 cAACGACAAGCAAGCCAAGGACTCAGGGCAGGGCGGCATGCAGGTGCGCTGGCGCACCGAGACGGTCGACTGGGGCCTGTACGCCATCCGTTTCCACGACAAGTCGCCGCAACTCAATGTACGCCCGaa                                                                >  2:817568/1‑129 (MQ=255)
                                                                                                 cAACGACAAGCAAGCCAAGGACTCAGGGCAGGGCGGCATGCAGGTGCGCTGGCGCACCGAGACGGTCGACTGGGGCCTGTACGCCATCCGTTTCCACGACAAGTCGCCGCAACTCAATGTACGCCCGaa                                                                <  1:817568/129‑1 (MQ=255)
                                                                                                  aaCGACAAGCAAGCCAAGGACTCAGGGCAGGGCGGCATGCAGGTGCGCTGGCGCACCGAGACGGTCGAc                                                                                                                           <  2:447686/69‑1 (MQ=255)
                                                                                                  aaCGACAAGCAAGCCAAGGACTCAGGGCAGGGCGGCATGCAGGTGCGCTGGCGCACCGAGACGGTCGAc                                                                                                                           >  1:447686/1‑69 (MQ=255)
                                                                                                      aCAAGCAAGCCAAGGACTCAGGGCAGGGCGGCATGCAGGTGCGCTGGCGCACCGAGACGGTCGACTGGGGCCTGTACGCCATCCGTTTCCACGACAAGTCGCCGCAACTCAATGTACGCCCGAACTTCGCCAACCTCGACCCGGTCAGCg                                      <  1:934114/150‑1 (MQ=255)
                                                                                                        aagcaagcCAAGGACTCAGGGCAGGGCGGCATGCAGGTGCGCTGGCGCACCGAGACGGTCGACTGGGGCCTGTACGCCATCCGTTTCCACGACAAGTCGCCGCAACTCAAt                                                                           <  2:1359037/111‑1 (MQ=255)
                                                                                                        aagcaagcCAAGGACTCAGGGCAGGGCGGCATGCAGGTGCGCTGGCGCACCGAGACGGTCGACTGGGGCCTGTACGCCATCCGTTTCCACGACAAGTCGCCGCAACTCAATGTACGCCCGAACTTCGCCAACCTCGACCCGGTCAGCGGc                                    >  1:405068/1‑150 (MQ=255)
                                                                                                        aagcaagcCAAGGACTCAGGGCAGGGCGGCATGCAGGTGCGCTGGCGCACCGAGACGGTCGACTGGGGCCTGTACGCCATCCGTTTCCACGACAAGCCGCCGCAACTCAAt                                                                           >  1:1359037/1‑111 (MQ=255)
                                                                                                              gcCAAGGACTCAGGGCAGGGCGGCATGCAGGTGCGCTGGCGCACCGAGACGGTCGACTGGGGCCTGTACGCCATCCGTTTCCACGACAAGTCGCCGCAACTCAATGTACGCCCGAACTTCGCCAACCTCGACCCGGTCAGCGGCAAGGcc                              >  1:1052919/1‑150 (MQ=255)
                                                                                                               ccAAGGACTCAGGGCAGGGCGGCATGCAGGTGCGCTGGCGCACCGAGACGGTCGACTGGGGCCTGTACGCCATCCGTTTCCACGACAAGTCGCCGCAACTCAATGTACGCCCGAACTTCGCCAACCTCGACCCGGTCAGCGGCAAGGCCg                             >  1:704142/1‑150 (MQ=255)
                                                                                                                cAAGGACTCAGGGCAGGGCGGCATGCAGGTGCGCTGGCGCACCGAGACGGTCGACTGGGGCCTGTACGCCATCCGTTTCCACGACAAGTCGCCGCAACTCAATGTACGCCCGAACTTCGCCAACCTCGACCCGGTCAGCGGCAAGGCCgg                            >  2:369976/1‑150 (MQ=255)
                                                                                                                 aaGGACTCAGGGCAGGGCGGCATGCAGGTGCGCTGGCGCACCGAGACGGTCGACTGGGGCCTGTACGCCATCCGTTTGCACGACAAGTCGCCGCAACTCAATGTACGCCCGAACTTCGCCAACCTCGACCCGGTCAGCGGCAcggccggc                           >  2:391698/1‑150 (MQ=255)
                                                                                                                 aaGGACTCAGGGCAGGGCGGCATGCAGGTGCGCTGGCGCACCGAGACGGTCGACTGGGGCCTGTACGCCATCCGTTTCCACGACAAGTCGCCGCAACTCAATGTACGCCCGAACTTCGCCAACCTCGACCCGGTCAGCGGCAAGGCCGGc                           <  1:965296/150‑1 (MQ=255)
                                                                                                                   ggACTCAGGGCAGGGCGGCATGCAGGTGCGCTGGCGCACCGAGACGGTCGACTGGGGCCTGTACGCCATCCGTTTCCACGACAAGTCGCCGCAGCGCAATGTACGCCCGAACTTCGCCAACCTCGACCCGGTCAGCGGCAAGGCCGGCGa                         <  2:275455/150‑1 (MQ=255)
                                                                                                                         agggcagggcGGCATGCAGGTGCGCTGGCGCACCGAGACGGTCGACTGGGGCCTGTACGCCATCCGTTTCCACGACAAGTCGCCGCAACTCAATGTACGCCCGAACTTCGCCAACCTCGACCCGGTCAGCGGCAAGGCCGGCGAGTTCTa                   >  2:1108732/1‑150 (MQ=255)
                                                                                                                          gggcagggcGGCATGCAGGTGCGCTGGCGCACCGAGACGGTCGACTGGGGCCTGTACGCCATCCGTTTCCACGACAAGTCGCCGCAACTCAATGTACGCCCGAACTTCGCCAACCTCGACCCGGTCAGCGGCAAGGCCGGCGAGTTCTAc                  >  1:1050738/1‑150 (MQ=255)
                                                                                                                              agggcGGCATGCAGGTGCGCTGGCGCACCGAGACGGTCGACTGGGGCCTGTACGCCATCCGTTTCCACGACAAGTCGCCGCAACTCAATGTACGCCCGAACTTCGCCAACCTCGACCCGGTCAGCGGCAAGGCCGGCGAGTTCTACTggg              <  1:404127/150‑1 (MQ=255)
                                                                                                                                ggcggcATGCAGGTGCGCTGGCGCACCGAGACGGTCGACTGGGGCCTGTACGCCATCCGTTTCCACGACAAGTCGCCGCAACTCAATGTACGCCCGAACTTCGCCAACCTCGACCCGGTCAGCGGCAAg                                 >  1:1094900/1‑129 (MQ=255)
                                                                                                                                ggcggcATGCAGGTGCGCTGGCGCACCGAGACGGTCGACTGGGGCCTGTACGCCATCCGTTTCCACGACAAGTCGCCGCAACTCAATGTACGCCCGAACTTCGCCAACCTCGACCCGGTCAGCGGCAAg                                 <  2:1094900/129‑1 (MQ=255)
                                                                                                                                   ggcATGCAGGTGCGCTGGCGCACCGAGANNGTCGACTGGGGCCTNTNCNCCATCCGTTTCCCCGACAAGTCGCCGCAACTCAATGTACGCCCGAACTTCGACAACCTCGACCCGGTCAGCGGCAAGGCCGGCGAGTTCTACTGGGTGTAt         >  1:835586/1‑150 (MQ=255)
                                                                                                                                     cgtctcggTGCGCCAGCGCACCTAGACGGTCGACTGGGGCCTGTACGCCATCCGTTTCCACGACAAGTCGCCGCAACTCAATGTACGCCCGAACTTCGCCAACCTCGACCCCGTCAACGACAGGGCCGGCGAGTGCTACTAGGTGTAGcc       >  2:180268/7‑150 (MQ=255)
                                                                                                                                      aTGCAGGTGCGCTGGCGCACCGAGACGGTCGACTGGGGCCTGTACGCCATCCGTTTCCACGACAAGTCGCCGCAACTCAATGTACGCCCGAACTTCGCCAACCTCGACCCGGTCAGCGGCAAGGCCGGCGAGTTCTACTGGGTGTCTccc      <  1:816499/150‑1 (MQ=255)
                                                                                                                                         tcggTGCGCCAGCGCACCGAGACGGTCGACTGGGGCCTGTACGCCATCCGTTTCCACGACAAGTCGCCGCAACTCAATGTACGCCCGAACTTCGCCAACCTCGACCCGGTCAGCGGCAAg                                 >  2:977142/3‑120 (MQ=255)
                                                                                                                                         tcggTGCGCCAGCGCACCGAGACGGTCGACTGGGGCCTGTACGCCATCCGTTTCCACGACAAGTCGCCGCAACTCAATGTACGCCCGAACTTCGCCAACCTCGACCCGGTCAGCGGCAAg                                 <  1:977142/118‑1 (MQ=255)
                                                                                                                                         cAGGTGCGCTGGCGCACCGAGACGGTCGACTGGGGCCTGTACGCCATCCGTTTCCACGACAAGTCGCCGCAACTCAATGTACGCCCGAACTTCGCCAACCTCGACCCGGTCAGCGGCAAGGCCGGCGAGTTCTACTGGGTGTATCCCGAg   <  2:1052919/150‑1 (MQ=255)
                                                                                                                                          aGGTGCGCTGGCGCACCGAGACGGTCGACTGGGGCCTGTACGCCATCCGTTTCCACGACAAGTCGCCGCAACTCAATGTACGCCCGAACTTCGCCAACCTCGACCCGGTCAGCGGCAAGGCCGGCGAGTTCTACTGGGTGTATCCCGAgg  <  2:191593/150‑1 (MQ=255)
                                                                                                                                                  |                                                                                                                                             
ATTTCTCCACAGACGACTTCTTTGGCGATGGCGCCGAGCGCATGTTCGTCGGCGCGCCGGTGTTCCCCAGTGGCCAACCTCTGGCGTTCTTCCACGGCAACGACAAGCAAGCCAAGGACTCAGGGCAGGGCGGCATGCAGGTGCGCTGGCGCACCGAGACGGTCGACTGGGGCCTGTACGCCATCCGTTTCCACGACAAGTCGCCGCAACTCAATGTACGCCCGAACTTCGCCAACCTCGACCCGGTCAGCGGCAAGGCCGGCGAGTTCTACTGGGTGTATCCCGAGG  >  NC_002947/3630219‑3630506

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