Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 4,158,879 A→G 8.9% S297G (AGC→GGC)  PP_3661 → membrane protein

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0029474,158,8790AG 8.9% 152.9 / 5.0 56S297G (AGC→GGC) PP_3661membrane protein
Reads supporting (aligned to +/- strand):  ref base A (30/21);  new base G (2/3);  total (32/24)
Fisher's exact test for biased strand distribution p-value = 6.42e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.46e-01

GCCTGTGGATCAGCCGTTCGCGCAAGGCAGGCCAGGCGCAGGGCAACGGGCGCATTGCCATGCCCTGGTTCGCTTTCGGCTTCCTGGCCCTGGTCCTGGTGAATTCGATGCAAGTGCTGCCTGGCAGCGTGACCCAGGCCGTGAACAGCCTGGACACCTTTGCCTTGACCATGGCCATGACGGCGTTGGGCATGGAGACCCGCTTCAGCCAGATTCGCCAGGCTGGGCCACGGGCGCTGGCCACCGGTGCGATCCTCAACCTGTGGCTGGTTGGCGGGGGGCTGGCGAT  >  NC_002947/4158733‑4159021
                                                                                                                                                  |                                                                                                                                              
gCCTGTGGATCAGCCGTTCGCGCAAGGCAGGCCAGGCGCAGGGCAACGGGCGCATTGCCATGCCCTGGTTCGCTTTCGGCTTCCTGGCCCTGGTCCTGGTGAATTCGATGCAAGTGCTGCCTGGCAGCGTGACCCAGGCCGTGAACAGc                                                                                                                                              >  4:253293/1‑149 (MQ=255)
gCCTGTGGATCAGCCGTTCGCGCAAGGCAGGCCAGGCGCAGGGCAACGGGCGCATTGCCATGCCCTGGTTCGCTTTCGGCTTCCTGGCCCTGGTCCTGGTGAATTCGATGCAAGTGCTGCCTGGCAGCGTGACCCAGGCCGTGAACAGc                                                                                                                                              <  6:266847/149‑1 (MQ=255)
 ccTGTGGATCAGCCGTTCGCGCAAGGCAGGCCAGGCGCAGGGCAACGGGCGCATTGCCATGCCCTGGTTCGCTTTCGGCTTCCTGGCCCTGGTCCTGGTGAATTCGATGCAAGTGCTGCCTGGCAGCGTGACCCAGGCCGTGAACAGcc                                                                                                                                             >  3:199467/1‑149 (MQ=255)
     tgGATCAGCCGTTCGCGCAAGGCAGGCCAGGCGCAGGGCAACGGGCGCATTGCCATGCCCTGGTTCGCTTGCGGCTTCCTGGCCCTGGTCCTGGTGAATTCGATGCAAGTGCTGCCTGGCAGCGTGACCCAGGCTGTTCACGGCCTgg                                                                                                                                          <  3:66225/148‑1 (MQ=255)
      ggATCAGCCGTTCGCGCAAGGCAGGCCAGGCGCAGGGCAACGGGCGCATTGCCATGCCCTGGTTCGCTTTCGGCTTCCTGGCCCTGGTCCTGGTGAATTCGATGCAAGTGCTGCCTGGCAGCGTGACCCAGGCCGTGAACAGCCTGGa                                                                                                                                         >  2:140507/1‑148 (MQ=255)
       gATCAGCCGTTCGCGCAAGGCAGGCCAGGCGCAGGGCAACGGGCGCATTGCCATGCCCTGGTTCGCTTTCGGCTTCCTGGCCCTGGTCCTGGTGAATTCGATGCAAGTGCTGCCTGGCAGCGTGACCCAGGCCGTGAACAGCCTGGaca                                                                                                                                       >  2:139363/1‑149 (MQ=255)
        aTCAGCCGTTCGCGCAAGGCAGGCCAGGCGCAGGGCAACGGGCGCATTGCCATGCCCTGGTTCGCTTTCGGCTTCCTGGCCCTGGTCCTGGTGAATTCGATGCAAGTGCTGCCTGGCAGCGTGACCCAGGCCGTGAACAGCCTGGacac                                                                                                                                      <  1:140507/149‑1 (MQ=255)
           aGCCGTTCGCGCAAGGCAGGCCAGGCGCAGGGCAACGGGCGCATTGCCATGCCCTGGTTCGCTTTCGGCTTCCTGGCCCTGGTCCTGGTGAATTCGATGCAAGTGCTGCCTGGCAGCGTGACCCAGGCCGTGAACAGCCTGGACACCt                                                                                                                                    <  1:274734/148‑1 (MQ=255)
            gCCGTTCGCGCAAGGCAGGCCAGGCGCAGGGCAACGGGCGCATTGCCATGCCCTGGTTCGCTTTCGGCTTCCTGGCCCTGGTCCTGGTGAATTCGATGCAAGTGCTGCCTGGCAGCGTGACCCAGGCCGTGAACAGCCTGGACACCttt                                                                                                                                  <  6:233846/149‑1 (MQ=255)
                  cgcgCAAGGCAGGCCAGGCGCAGGGCAACGGGCGCATTGCCATGCCCTGGTTCGCTTTCGGCTTCCTGGCCCTGGTCCTTGTGAATTCGATGCAAGTGCTGCCTGGCAGCGTGACCCAGGCCGTGAACAGCCTGGACACCTTTGCCtt                                                                                                                             <  5:214783/148‑1 (MQ=255)
                     gcAAGGCAGGCCAGGCGCAGGGCAACGGGCGCATTGCCATGCCCTGGTTCGCTTTCGGCTTCCTGGCCCTGGTCCTGGTGAATTCGATGCAAGTGCTGCCTGGCAGCGTGACCCAGGCCGTGAACAGCCTGGACACCTTTGCCTTGAc                                                                                                                          >  3:259602/1‑148 (MQ=255)
                          gcaggcCAGGCGCAGGGCAACGGGCGCATTGCCATGCCCTGGTTCGCTTTCGGCTTCCTGGCCCTGGTCCTGGTGAATTCGATGCAAGTGCTGCCTGGCAGCGTGACCCAGGCCGTGAACAGCCTGGACACCTTTGCCTTGACCATGGc                                                                                                                    <  3:253293/149‑1 (MQ=255)
                                            aaCGGGCGCATTGCCATGCCCTGGTTCGCTTTCGGCTTCCTGGCCCTGGTCCTGGTGAATTCGATGCAAGTGCTGCCTGGCAGCGTGACCCAGGCCGTGAACAGCCTGGACACCTTTGCCTTGACCATGGCCATGACGGCGTTGGGCAt                                                                                                  >  3:296780/1‑149 (MQ=255)
                                               gggCGCATTGCCATGCCCTGGTTCGCTTTCGGCTTCCTGGCCCTGGTCCTGGTGAATTCGATGCAAGTGCTGCCTGGCAGCGTGACCCAGGCCGTGAACAGCCTGGACACCTTTGCCTTGACCATGGCCATGACGGCGTTGGGCATGga                                                                                               <  4:296780/149‑1 (MQ=255)
                                                   gcATTGCCATGCCCTGGTTCGCTTTCGGCTTCCTGGCCCTGGTCCTGGTGAATTCGATGCAAGTGCTGCCTGGCAGCGTGACCCAGGCCGTGAACAGCCTGGACACCTTTGCCTTGACCATGGCCATGACGGCGTTGGGCATGGAGAcc                                                                                           >  4:16359/1‑149 (MQ=255)
                                                    cATTGCCATGCCCTGGTTCGCTTTCGGCTTCCTGGCCCTGGTCCTGGTGAATTCGATGCAAGTGCTGCCTGGCAGCGTGACCCAGGCCGTGAACAGCCTGGACACCTTTGCCTTGACCATGGCCATGACGGCGTTGGGCATGGAGAccc                                                                                          >  2:150424/1‑149 (MQ=255)
                                                     aTTGCCATGCCCTGGTTCGCTTTCGGCTTCCTGGCCCTGGTCCTGGTGAATTCGATGCAAGTGCTGCCTGGCAGCGTGACCCAGGCCGTGAACAGCCTGGACACCTTTGCCTTGACCATGGCCATGACGGCGTTGGGCATGGAGACCCg                                                                                         >  2:38574/1‑149 (MQ=255)
                                                     aTTGCCATGCCCTGGTTCGCTTTCGGCTTCCTGGCCCTGGTCCTGGTGAATTCGATGCAAGTGCTGCCTGGCAGCGTGACCCAGGCCGTGAACAGCCTGGACACCTTTGCCTTGACCATGGCCATGACGGCGTTGGGCATGGAGACCCg                                                                                         >  2:181827/1‑149 (MQ=255)
                                                          cATGCCCTGGTTCGCTTTCGGCTTCCTGGCCCTGGTCCTGGTGAATTCGATGCAAGTGCTGCCTGGCAGCGTGACCCAGGCCGTGAACAGCCTGGac                                                                                                                                        <  6:235754/97‑1 (MQ=255)
                                                          cATGCCCTGGTTCGCTTTCGGCTTCCTGGCCCTGGTCCTGGTGAATTCGATGCAAGTGCTGCATGGCAGCGTGACCCAGGCCGTGAACAGCCTGGac                                                                                                                                        >  5:235754/1‑97 (MQ=255)
                                                             gCCCTGGTTCGCTTTCGGCTTCCTGGCCCTGGTCCTGGTGAATTCGATGCAAGTGCTGCCTGGCAGCGTGACCCAGGCCGTGAACAGCCTGGACACCTTTGCCTTGACCATGGCCATGACGGCGTTGGGCATGGAGACACGCTTCAGc                                                                                  >  8:259199/1‑148 (MQ=255)
                                                              cccTGGTTCGCTTTCGGCTTCCTGGCCCTGGTCCTGGTGAATTCGATGCAAGTGCTGCCTGGCAGCGTGACCCAGGCCGTGAACAGCCTGGACACCTTTGCCTTGACCATGGCCATGACGGCGTTGGGCATGGAGACCCGCTTCAGcc                                                                                 <  4:259602/148‑1 (MQ=255)
                                                               ccTGGTTCGCTTTCGGCTTCCTGGCCCTGGTCCTGGTGAATTCGATGCAAGTGCTGCCTGGCAGCGTGACCCAGGCCGTGAACAGCCTGGACACCTTTGCCTTGACCATGGCCATGACGGCGTTGGGCATGGAGACCCGCTTCAGc                                                                                  >  8:183666/1‑146 (MQ=255)
                                                                    ttCGCTTTCGGCTTCCTGGCCCTGGTCCTGGTGAATTCGATGCAAGTGCTGCCTGGCAGCGTGACCCAGGCCGTGAACAGCCTGGACACCTTTGCCTTGACCATGGCCATGACGGCGTTGGGCATGGAGACCCGCTTCAGCCAGATTCg                                                                          >  4:59993/1‑149 (MQ=255)
                                                                       gCTTTCGGCTTCCTGGCCCTGGTCCTGGTGAATTCGATGCAAGTGCTGCCTGGCAGCGTGACCCAGGCCGTGAACAGCCTGGACACCTTTGCCTTGACCATGGCCATGACGGCGTTGGGCATGGAGACCCGCTTCAGCCAGATTCGcc                                                                        <  4:199467/148‑1 (MQ=255)
                                                                         tttCGGCTTCCTGGCCCTGGTCCTGGTGAATTCGATGCAAGTGCTGCCTGGCAGCGTGACCCAGGCCGTGAACAGCCTGGACACCTTTGCCTTGACCATGGCCATGACGGCGTTGGGCATGGAGACCCGCTTCAGCCAGATTCGCCAgg                                                                     >  3:70170/1‑149 (MQ=255)
                                                                               cTTCCTGGCCCTGGTCCTGGTGAATGCGATGCAAGTGCTGCCTGGCAGCGTGACCCAGGCCGTGAACAGCCTGGACACCTTTGCCTTGACCATGGCCATGACGGCGTTGGGCATGGAGACCCGCTTCAGCCAGATTCGCCAGGCTGGGc                                                               <  8:157900/149‑1 (MQ=255)
                                                                                    tGGCCCTGGTCCTGGTGAATTCGATGCAAGTGCTGCCTGGCAGCGTGACCCAGGCCGTGAACAGCCTGGACACCTTTGCCTTGACCATGGCCATGACGGCGTTGGGCATGGAGACCCGCTTCAGCCAGATTCGCCAGGCTGGGCCACgg                                                          >  4:180113/1‑149 (MQ=255)
                                                                                        cctggtcctggtGAATTCGATGCAAGTGCTGCCTGGCAGCGTGACCCAGGCCGTGAACAGCCTGGACACCTTTGCCTTGACCATGGCCATGACGGCGTTGGGCATGGAGACCCGCTTCAGCCAGATTCGCCAGGCTGGGCCACGGgcgc                                                      >  7:12075/1‑149 (MQ=255)
                                                                                                tggtGAATTCGATGCAAGTGCTGCCTGGCAGCGTGACCCAGGCCGTGAACAGCCTGGACACCTTTGCCTTGACCATGGCCATGACGGCGTTGGGCATGGAGACCCGCTTCAGCCAGATTCGCCAGGCTGGGCCACGGGCGCTGGCCAcc                                              >  1:245834/1‑149 (MQ=255)
                                                                                                ccgcccaTTCGATGCAAGTGCTGCCTGGCAGCGTGACCCAGGCCGTGAACAGCCTGGACACCTTTGCCTTGACCATGGCCATGACGGCGTTGGGCATGGAGACCCGCTTCAGCCAGATTCGCCAGGCTGGGCCACGGGCGCTGGCCAc                                               <  2:329069/142‑1 (MQ=255)
                                                                                                 ggtGAATTCGATGCAAGTGCTGCCTGGCAGCGTGACCCAGGCCGTGAACAGCCTGGACACCTTTGCCTTGACCATGGCCATGACGGCGTTGGGCATGGAGACCCGCTTCAGCCAGATTCGCCAGGCTGGGCCACGGGCGCTGGCCACCg                                             <  1:139363/149‑1 (MQ=255)
                                                                                                            tGCAAGTGCTGCCTGGCAGCGTGACCCAGGCCGTGAACAGCCTGGACACCTTTGCCTTGACCATGGCCATGACGGCGTTGGGCATGGAGACCCGCTTCAGCCAGATTCGCCAGGCTGGGCCACGGGCGCTGGCCACCGGTGCGATCCTc                                  <  7:183666/149‑1 (MQ=255)
                                                                                                               aaGTGCTGCCTGGCAGCGTGACCCAGGCCGTGAACAGCCTGGACACCTTTGCCTTGACCATGGCCATGACGGTGTTGGGCATGGAGACCCGCTTCAGCCAGATTCGCCAGGCTGGGCCACGGGCGCTGGCCACCGGTGCGATCCTCAAc                               >  4:177337/1‑149 (MQ=255)
                                                                                                               aaGTGCTGCCTGGCAGCGTGACCCAGGCCGTGAACAGCCTGGACACCTTTGCCTTGACCATGGCCATGACGGCGTTGGGCATGGAGACCCGCTTCAGCCAGATTCGCCAGGCTGGGCCACGGGCGCTGGCCACCGGTGCGATCCTCAAc                               >  5:103484/1‑149 (MQ=255)
                                                                                                               aaGTGCTGCCTGGCAGCGTGACCCAGGCCGTGAACAGCCTGGACACCTTTGCCTTGACCATGGCCATGACGGCGTTGGGCATGGAGACCCGCTTCAGCCAGATTCGCCAGGCTGGGCCACGGGCGCTGGCCACCGGTGCGATCCTCAAc                               >  2:179937/1‑149 (MQ=255)
                                                                                                                 gTGCTGCCTGGCAGCGTGACCCTGGCCGTGAACAGCCTGGACACCTTTGCCTTGACCATGGCCATGACGGCGTTGGGCATGGAGACCCGCTTCAGCCAGATTCGCCAGGCTGGGCCACGGGCGCTGGCCACCGGTGCGATCCTCAACCt                             >  5:197452/1‑149 (MQ=255)
                                                                                                                    ctgcCTGGCAGCGTGACCCAGGCCGTGAACAGCCTGGACACCTTTGCCTTGACCATGGCCATGACGGCGTTGGGCATGGAGACCCGCTTCAGCCAGATTCGCCAGGCTGGGCCACGGGCGCTGGCCACCGGTGCGATCCTCAACCtgtg                          >  4:104144/1‑149 (MQ=255)
                                                                                                                    ctgcCTGGCAGCGTGACCCAGGCCGTGAACAGCCTGGACACCTTTGCCTTGACCATGGCCATGACGGCGTTGGGCATGGAGACCCGCTTCAGCCAGATTCGCCAGGCTGGGCCACGGGCGCTGGCCACCGGTGCGATCCTCAACCtgt                           <  8:12075/148‑1 (MQ=255)
                                                                                                                    ctgcCTGGCAGCGTGACCCAGGCCGTGAACAGCCTGGACACCTTTGCCTTGACCATGGCCATGACGGCGTTGGGCATGGAGACCCGCTTCAGCCAGATTCGCCAGGCTGGGCCACGGGCGCTGGCCACCGGTGCGATCCTCAACCtgt                           <  1:38574/148‑1 (MQ=255)
                                                                                                                    ctgcCTGGCAGCGTGACCCAGGCCGTGAACAGCCTGGACACCTTTGCCTTGACCATGGCCATGACGGCGTTGGGCATGGAGACCCGCTTCAGCCAGATTCGCCAGGCTGGGCCACGGGCGCCGGCCACCGGTGCGATCCTCAACCtgt                           <  1:181827/148‑1 (MQ=255)
                                                                                                                     tgcCTGGCAGCGTGACCCAGGCCGTGAACAGCCTGGACACCTTTGCCTTGACCATGGCCATGACGGCGTTGGGCATGGAGACCCGCTTCAGCCAGATTc                                                                           <  5:256808/99‑1 (MQ=255)
                                                                                                                     tgcCTGGCAGCGTGACCCAGGCCGTGAACAGCCTGGACACCTTTGCCTTGACCATGGCCATGACGGCGTTGGGCATGGAGACCCGCTTCAGCCAGATTc                                                                           >  6:256808/1‑99 (MQ=255)
                                                                                                                       ccTGGCAGCGTGACCCAGGCCGTGAACAGCCTGGACACCTTTGCCTTGACCATGGCCATGACGGCGTTGGGCATGGAGACCCGCTTCAGCCAGATTCGCCAGGCTGGGCCACGGGCGCTGGCCACCGGTGCGATCCTCAACCTGTGGCt                       >  8:208010/1‑149 (MQ=255)
                                                                                                                           gCAGCGTGACCCTGGCCGTGAACAGCCTGGACACCTTTGCCTTGACCATGGCCATGACGGCGTTGGGCATGGAGACCCGCTTCAGCCAGATTCGCCAGGCTGGGCCACGGGCGCTGGCCACCGGTGCGATCCTCAACCTGTGGCTGGtt                   <  6:197452/149‑1 (MQ=255)
                                                                                                                           gCAGCGTGACCCAGGCCGTGAACAGCCTGGACACCTTTGCCTTGACCATGGCCATGACGGCGTTGGGCATGGAGACCCGCTTCAGCCAGATTCGCCAGGCTGGGCCACGGGCGCTGGCCAGCGGTGCGATCCTCAACCTGTGGCtcgat                   >  4:215285/1‑145 (MQ=255)
                                                                                                                               aggtgtccAGGCTGTTCACGGCCTGGGCACCTTTGCCTTGACCATGGCCATGACGGCGTTGGGCATGGAGACCCGCTTCAGCCAGATTCGCCAGGCTGGGCCACGGGCGCTGGCCACCGGt                                           <  6:236884/115‑1 (MQ=255)
                                                                                                                               aggtgtccAGGCTGTTCACGGCCTGGGCACCTTTGCCTTGACCATGGCCATGACGGCGTTGGGCATGGAGACCCGCTTCAGCCAGATTCGCCAGGCTGGGCCACGGGCGCTGGCCACCGGt                                           >  5:236884/7‑121 (MQ=255)
                                                                                                                                gTGACCCAGGCCGTGAACAGCCTGGACACCTTTGCCTTGACCATGGCCATGACGGCGTTGGGCATGGAGACCCGCTTCAGCCAGATTCGCCAGGCTGGGCCACGGGCGCTGGCCACCGGTGCGATCCTCAACCTGTGGCTGGTTGGCgg              <  3:180113/149‑1 (MQ=255)
                                                                                                                                    tccAGGCTGTTCACGGCCTGGACACCTTTGCCTTGACCATGGCCATGACGGCGTTGGGCATGGAGACCCGCTTCAGCCAGATTCGCCAGGCTGGGGCACGGGCGCTGGCCACCGGTGCGATCCTCAACCtgtg                          >  1:279084/2‑133 (MQ=255)
                                                                                                                                    tccAGGCTGTTCACGGCCTGGACACCTTTGCCTTGACCATGGCCATGACGGCGTTGGGCATGGAGACCCGCTTCAGCCAGATTCGCCAGGCTGGGCCACGGGCGCTGGCCACCGGTGCGATCCTCAACCtgtg                          <  2:279084/132‑1 (MQ=255)
                                                                                                                                    cccAGGCCGTGAACAGCCTGGACACCTTTGCCTTGACCATGGCCATGACGGCGTTGGGCATGGAGACCCGCTTCAGCCAGATTCGCCAGGCTGGGCCACGGGCGCTGGCCACCGGTGCGATCCTCAACCTGTgg                         <  4:223071/134‑1 (MQ=255)
                                                                                                                                    cccAGGCCGTGAACAGCCTGGACACCTTTGCCTAGACCATGGCCATGACGGCGTTGGGCATGGAGACCCGCTTCAGCCAGATTCGCCAGGCTGGGCCACGGGCGCTGGCCACCGGTGCGATCCTCAACCTGTgg                         >  3:223071/1‑134 (MQ=255)
                                                                                                                                      cAGGCCGTGAACAGCCTGGACACCTTTGCCTTGACCATGGCCATGACGGCGTTGGGCATGGAGACCCGCTTCAGCCAGATTCGCCAGGCTGGGCCACGGGCGCTGGCCACCGGTGCGATCCTCAACCTGTGGCTGGTTGGCGGGGGGCt        >  6:224823/1‑149 (MQ=255)
                                                                                                                                       aGGCCGTGAACAGCCTGGACACCTTTGCCTTGACCATGGCCATGACGGCGTTGGGCATGGAGACCCGCTTCAGCCAGATTCGCCAGGCTGGGCCACGGGCGCTGGCCACCGGTGCGATCCTCAACCTGTGGCTGGTTGGCGGGGGGCTg       >  3:27840/1‑149 (MQ=255)
                                                                                                                                            gTGAACAGCCTGGACACCTTTGCCTTGACCATGGCCATGACGGCGTTGGGCATGGAGACCCGCTTCAGCCAGATTCGCCAGGCTGGGCCACGGGCGCTGGCCACCGGTGCGATCCTCAACCTGTGGCTGGTTGGCGGGGGGCTGGCGAt  >  8:65032/1‑149 (MQ=255)
                                                                                                                                                  |                                                                                                                                              
GCCTGTGGATCAGCCGTTCGCGCAAGGCAGGCCAGGCGCAGGGCAACGGGCGCATTGCCATGCCCTGGTTCGCTTTCGGCTTCCTGGCCCTGGTCCTGGTGAATTCGATGCAAGTGCTGCCTGGCAGCGTGACCCAGGCCGTGAACAGCCTGGACACCTTTGCCTTGACCATGGCCATGACGGCGTTGGGCATGGAGACCCGCTTCAGCCAGATTCGCCAGGCTGGGCCACGGGCGCTGGCCACCGGTGCGATCCTCAACCTGTGGCTGGTTGGCGGGGGGCTGGCGAT  >  NC_002947/4158733‑4159021

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