Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 3,630,366 G→A 9.2% W315* (TGG→TAG) ‡ PP_3200 → hypothetical protein

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0029473,630,3660GA 9.2% 143.0 / 5.4 54W315* (TGG→TAG) ‡PP_3200hypothetical protein
Reads supporting (aligned to +/- strand):  ref base G (26/23);  new base A (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

ATTTCTCCACAGACGACTTCTTTGGCGATGGCGCCGAGCGCATGTTCGTCGGCGCGCCGGTGTTCCCCAGTGGCCAACCTCTGGCGTTCTTCCACGGCAACGACAAGCAAGCCAAGGACTCAGGGCAGGGCGGCATGCAGGTGCGCTGGCGCACCGAGACGGTCGACTGGGGCCTGTACGCCATCCGTTTCCACGACAAGTCGCCGCAACTCAATGTACGCCCGAACTTCGCCAACCTCGACCCGGTCAGCGGCAAGGCCGGCGAGTTCTACTGGGTGTATCCCGAGGGTATCGAGG  >  NC_002947/3630219‑3630515
                                                                                                                                                   |                                                                                                                                                     
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                                                                                                                    <  1:1064259/119‑1 (MQ=255)
                                                                ccccAGTGGCCAACCTCTGGCGTTCTTCCACGGCAACGACAAGCAAGCCAAGGACTCAGGGCAGGGCGGCATGCAGGTGCGCTGGCGCACCGAGACGGTCGACTGGGGCCTGTACGCCa                                                                                                                    >  2:1064259/1‑119 (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                                                                         <  1:817568/129‑1 (MQ=255)
                                                                                                 cAACGACAAGCAAGCCAAGGACTCAGGGCAGGGCGGCATGCAGGTGCGCTGGCGCACCGAGACGGTCGACTGGGGCCTGTACGCCATCCGTTTCCACGACAAGTCGCCGCAACTCAATGTACGCCCGaa                                                                         >  2:817568/1‑129 (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                                          <  1:977142/118‑1 (MQ=255)
                                                                                                                                         tcggTGCGCCAGCGCACCGAGACGGTCGACTGGGGCCTGTACGCCATCCGTTTCCACGACAAGTCGCCGCAACTCAATGTACGCCCGAACTTCGCCAACCTCGACCCGGTCAGCGGCAAg                                          >  2:977142/3‑120 (MQ=255)
                                                                                                                                         cAGGTGCGCTGGCGCACCGAGACGGTCGACTGGGGCCTGTACGCCATCCGTTTCCACGACAAGTCGCCGCAACTCAATGTACGCCCGAACTTCGCCAACCTCGACCCGGTCAGCGGCAAGGCCGGCGAGTTCTACTGGGTGTATCCCGAg            <  2:1052919/150‑1 (MQ=255)
                                                                                                                                          aGGTGCGCTGGCGCACCGAGACGGTCGACTGGGGCCTGTACGCCATCCGTTTCCACGACAAGTCGCCGCAACTCAATGTACGCCCGAACTTCGCCAACCTCGACCCGGTCAGCGGCAAGGCCGGCGAGTTCTACTGGGTGTATCCCGAgg           <  2:191593/150‑1 (MQ=255)
                                                                                                                                                   ggCGCACCGAGACGGTCGACTGGGGCCTGTACGCCATCCGTTTCCACGACAAGTCGCCGCAACTCAATGTACGCCCGAACTTCGCCAACCTCGACCCGGTCAGCGGCAAGGCCGGCGAGTTCTACTGGGTGTATCCCGAGGGTATCGAgg  <  1:369976/150‑1 (MQ=255)
                                                                                                                                                   |                                                                                                                                                     
ATTTCTCCACAGACGACTTCTTTGGCGATGGCGCCGAGCGCATGTTCGTCGGCGCGCCGGTGTTCCCCAGTGGCCAACCTCTGGCGTTCTTCCACGGCAACGACAAGCAAGCCAAGGACTCAGGGCAGGGCGGCATGCAGGTGCGCTGGCGCACCGAGACGGTCGACTGGGGCCTGTACGCCATCCGTTTCCACGACAAGTCGCCGCAACTCAATGTACGCCCGAACTTCGCCAACCTCGACCCGGTCAGCGGCAAGGCCGGCGAGTTCTACTGGGTGTATCCCGAGGGTATCGAGG  >  NC_002947/3630219‑3630515

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: ‑