Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 436,899 C→A 7.7% W130L (TGG→TTG)  PP_0358 ← membrane protein

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_002947436,8990CA 7.7% 162.6 / 4.2 52W130L (TGG→TTG) PP_0358membrane protein
Reads supporting (aligned to +/- strand):  ref base C (25/23);  new base A (2/2);  total (27/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 = 3.54e-01

GCCGGGCAGCGCAGCCACCGGCGCCTGCTTACGGATCAGCCCATAGAAGCCGAAGCTCAACGCCAGCACCAGTGATACCCACGGCAGGCTGCCCACTTGCCACACCTGCTGCGCCACACCCACGGCAGCCATGCCCACCGCCAGCCATTGCAGGCGGCGCAGGCGCTCACCAAGAATCAGCATGCCCAGCAGTACATTGATCAGCGGGTTGATGTAATAGCCCAGGCTGGCCTCGAGCATACGGCCGTTGTTCACCGACCACACGTAGGTCAGCCAGTTGCCGGCGATCAGTGC  >  NC_002947/436754‑437047
                                                                                                                                                 |                                                                                                                                                    
gCCGGGCAGCGCAGCCACCGGCGCCTGCTTACGGATCAGCCCATAGAAGCCGAAGCTCAACGCCAGCACCAGTGATACCCACGGCAGGCTGCCCACTTGCCACACCTGCTGCGCCACACCCACGGCAGCCATGCCCACCGCCAGCCAtt                                                                                                                                                   >  4:162555/1‑149 (MQ=255)
       agcgcagcCACCGGCGCCTGCTTACGGATCAGCCCATAGAAGCCGAAGCTCAACGCCAGCACCAGTGATACCCACGGCAGGCTGCCCACTTGCCACACCTGCTGCGCCACACCCACGGCAGCCATGCCCACCGCCAGCCATTGCAggcg                                                                                                                                            <  3:162555/149‑1 (MQ=255)
                           cTTACGGATCAGCCCATAGAAGCCGAAGCTCAACGCCAGCACCAGTGATACCCACGGCAGGCTGCCCACTTGCCACACCTGCTGCGCCACACCCACGGCAGCCATGCCCACCGCCAGCCATTGCAGGCGGCGCAGGCGCTCACCAAGa                                                                                                                         >  2:23107/1‑148 (MQ=255)
                            ttACGGATCAGCCCATAGAAGCCGAAGCTCAACGCCAGCACCAGTGATACCCACGGCAGGCTGCCCACTTGCCACACCTGCTGCGCCACACCCACGGCAGCCATGCCCACCGCCTGCAATGGCTggcggtg                                                                                                                                         >  4:320467/1‑129 (MQ=255)
                            ttACGGATCAGCCCATAGAAGCCGAAGCTCAACGCCAGCACCAGTGATACCCACGGCAGGCTGCCCACTTGCCACACCTGCTGCGCCACACCCACGGCAGCCATGCCCACCGCCTGCAATGGCTggcggtg                                                                                                                                         <  3:320467/131‑3 (MQ=255)
                            ttACGGATCAGCCCATAGAAGCCGAAGCTCAACGCCAGCACCAGTGATACCCACGGCAGGCTGCCCACTTGCCACACCTGCTGCGCCACACCCACGGCAGCCATGCCCACCGCCAGCCATTGCAGGCGGCGCAGGCGCTCACCAAGaa                                                                                                                        >  1:49484/1‑148 (MQ=255)
                                ggATCAGCCCATAGAAGCCGAAGCTCAACGCCAGCACCAGTGATACCCACGGCAGGCTGCCCACTTGCCACACCTGCTGCGCCACACCCACGGCAGCCATGCCCACCGCCAACCATTGCAGGCGGCGCAGGCGCTCACCAAGa                                                                                                                         >  6:333557/1‑143 (MQ=255)
                                ggATCAGCCCATAGAAGCCGAAGCTCAACGCCAGCACCAGTGATACCCACGGCAGGCTGCCCACTTGCCACACCTGCTGCGCCACACCCACGGCAGCCATGCCCACCGCCAACCATTGCAGGCGGCGCAGGCGCTCACCAAGa                                                                                                                         <  5:333557/143‑1 (MQ=255)
                                      gCCCATAGAAGCCGAAGCTCAACGCCAGCACCAGTGATACCCACGGCAGGCTGCCCACTTGCCACACCTGCTGCGCCACACCCACGGCAGCCATGCCCACCGCCAGCCATTGCAGGCGGCGCAGGCGCTCACCAAGa                                                                                                                         >  6:108494/1‑137 (MQ=255)
                                      gCCCATAGAAGCCGAAGCTCAACGCCAGCACCAGTGATACCCACGGCAGGCTGCCCACTTGCCACACCTGCTGCGCCACACCCACGGCAGCCATGCCCACCGCCAGCCATTGCAGGCGGCGCAGGCGCTCACCAAGa                                                                                                                         <  5:108494/137‑1 (MQ=255)
                                        ccATAGAAGCCGAAGCTCAACGCCAGCACCAGTGATACCCACGGCAGGCTGCCCACTTGCCACACCTGCTGCGCCACACCCACGGCAGCCATGCCCACCGCCAGCCATTGCAGGCGGCGCAGGCGCTCACCAAGAATCAGCATGCCcag                                                                                                           >  5:1381/1‑149 (MQ=255)
                                          aTAGAAGCCGAAGCTCAACGCCAGCACCAGTGATACCCACGGCAGGCTGCCCACTTGCCACACCTGCTGCGCCACACCCACGGCAGCCATGCCCACCGCCAGCCATTGCAGGCGGCGCAGGCGCTCACCAAGAATCAGCATGCCcag                                                                                                           >  5:96768/1‑147 (MQ=255)
                                            aGAAGCCGAAGCTCAACGCCAGCACCAGTGATACCCACGGCAGGCTGCCCACTTGCCACACCTGCCGCGCCACACCCACGGCAGCCATGCCCACCGCCAGCCATTGCAGGCGGCGCAGGCGCTCACCAAGAAt                                                                                                                       >  5:67111/1‑133 (MQ=255)
                                            aGAAGCCGAAGCTCAACGCCAGCACCAGTGATACCCACGGCAGGCTGCCCACTTGCCACACCTGCCGCGCCACACCCACGGCAGCCATGCCCACCGCCAGCCATTGCAGGCGGCGCAGGCGCTCACCAAGAAt                                                                                                                       <  6:67111/133‑1 (MQ=255)
                                              aaGCCGAAGCTCAACGCCAGCACCAGTGATACCCACGGCAGGCTGCCCACTTGCCACACCTGCTGCGCCACACCCACGGCAGCCATGCCCACCGCCAGCCATTGCAGGCGGCGCAGGCGCTCACCAAGAATCAGCATGCCcag                                                                                                           <  6:340993/143‑1 (MQ=255)
                                              aaGCCGAAGCTCAACGCCAGCACCAGTGATACCCACGGCAGGCTGCCCACTTGCCACACCTGCTGCGCCACACCCACGGCAGCCATGCCCACCGCCAGCCATTGCAGGCGGCGCAGGCGCTCACCAAGAATCAGCATGCCcag                                                                                                           >  5:340993/1‑143 (MQ=255)
                                                gCCGAAGCTCAACGCCAGCACCAGTGATACCCACGGCAGGCTGCCCACTTGCCACACCTGCTGCGCCTCACCCACGGCAGCCATGCCCACCGCCAGCCATTGCAGGCGGCGCAGGCGCTCACCAAGAATCAGCATGCCCAGCAGTACAt                                                                                                   <  6:1381/149‑1 (MQ=255)
                                                  cGAAGCTCAACGCCAGCACCAGTGATACCCACGGCAGGCTGCCCACTTGCCACACCTGCTGCGCCACACCCACGGCAGCCATGCCCACCGCCAGCCATTGCAGGCGGCGCAGGCGCTCACCAAGAATCAGCATGCCCAGCAGTACAtt                                                                                                  >  3:36745/1‑148 (MQ=255)
                                                    aaGCTCAACGCCAGCACCAGTGATACCCACGGCAGGCTGCCCACTTGCCACACCTGCTGCGCCACCCCCACGGCAGCCATGCCCACCGCCAGCCATTGCAGGCGGCGCAGGCGCTCACCAAGAATCAGCATGCCCAGCAGTACATTGAt                                                                                               >  7:85949/1‑149 (MQ=255)
                                                        tCAACGCCAGCACCAGTGATACCCACGGCAGGCTGCCCACTTGCCACACCTGCTGCGCCACACCCACGGCAGCCATGCCCACCGCCAGCCATTGCAGGCGGCGGAGGCGCTCACCAAGAATCa                                                                                                                     <  8:261664/123‑1 (MQ=255)
                                                        tCAACGCCAGCACCAGTGATACCCACGGCAGGCTGCCCACTTGCCACACCTGCTGCGCCACACCCACGGCAGCCATGCCCACCGCCAGCCATTGCAGGCGGCGCAGGCGCTCACCAAGAATCAGCATGCCCAGCAGTACATTGATCAGc                                                                                           >  4:207337/1‑149 (MQ=255)
                                                          aaCGCCAGCACCAGTGATACCCACGGCAGGCTGCCCACTTGCCACACCTGCTGCGCCACACCCACGGCAGCCATGCCCACCgccagcc                                                                                                                                                      >  1:151322/1‑88 (MQ=255)
                                                          aaCGCCAGCACCAGTGATACCCACGGCAGGCTGCCCACTTGCCACACCTGCTGCGCCACACCCACGGCAGCCATGCCCACCgccagcc                                                                                                                                                      <  2:151322/88‑1 (MQ=255)
                                                           aCGCCAGCACCAGTGATACCCACGGCAGGCTGCCCACTTGCCACACCTGCTGCGCCACACCCACGGCAGCCATGCCCACCGCCTGCAATGGCTggcggt                                                                                                                                          <  8:72675/99‑2 (MQ=255)
                                                           aCGCCAGCACCAGTGATACCCACGGCAGGCTGCCCACTTGCCACACCTGCTGCGCCACACCCACGGCAGCCATGCCCACCGCCTGCAATGGCTggcggt                                                                                                                                          >  7:72675/1‑98 (MQ=255)
                                                                aGCACCAGTGATACCCACGGCAGGCTGCCCACTTGCCACACCTGCTGCGCCACACCCACGGCAGCCATGCCCACCGCCAGCCATTGCAGGCGGCGCAGGCGCTCACCAAGAATCAGCATGCCCAGCAGTACATTGATCAGCGGGTTGAt                                                                                   <  4:36745/149‑1 (MQ=255)
                                                                 gCACCAGTGATACCCACGGCAGGCTGCCCACTTGCCACACCTGCTGCGCCACACCCACGGCAGCCATGCCCACCGCCAGCCATTGCAGGCGGCGCAGGCGCTCACCAAGAATCAGCATGCCCAGCAGTACATTGATCAGCGGGTTGATg                                                                                  <  1:23107/149‑1 (MQ=255)
                                                                       gTGATACCCACGGCAGGCTGCCCACTTGCCACACCTGCTGCGCCACACCCACGGCAGCCATGCCCACCGCCAGCCATTGCAGGCGGCGCAGGCGCTCACCAAGAATCAGCATGCCCAGCAGTACATTGATCAGCg                                                                                          >  5:334023/1‑135 (MQ=255)
                                                                       gTGATACCCACGGCAGGCTGCCCACTTGCCACACCTGCTGCGCCACACCCACGGCAGCCATGCCCACCGCCAGCCATTGCAGGCGGCGCAGGCGCTCACCAAGAATCAGCATGCCCAGCAGTACATTGATCAGCg                                                                                          <  6:334023/135‑1 (MQ=255)
                                                                                  ggCAGGCTGCCCACTTGCCACACCTGCTGCGCCACACCCACGGCAGCCATGCCCACCGCCAGCCATTGCAGGCGGCGCAGGCGCTCACCAAGAATCAGCATGCCCAGCAGTACATTGATCAGCGGGTTGATGTAATAGCCCAGGCTGGc                                                                 <  3:207337/149‑1 (MQ=255)
                                                                                       gCTGCCCACTTGCCACACCTGCTGCGCCACACCCACGGCAGCCATGCCCACCGCCAGCCATTGCAGGCGGCGCAGGCGCTCACCAAGAATCAGCATGCCCAGCAGTACATTGATCAGCGGGTTGATGTAATAGCCCAGGCTGGCCTCGa                                                            <  2:49484/149‑1 (MQ=255)
                                                                                       gCTGCCCACTTGCCACACCTGCTGCGCCACACCCACGGCAGCCATGCCCACCGCCAGCCATTGCAGGCGGCGCAGGCGCTCACCAAGAATCAGCATGCCCAGCAGTACATTGATCAGCGGGTTGATGTAATAGCCCAGGCTGGCCTCGa                                                            <  5:339208/149‑1 (MQ=255)
                                                                                        cTGCCCACTTGCCACACCTGCTGCGCCACACCCACGGGAGCCATGCCCACCGCCAGCCATTGCAGGCGGCGCAGTCGCTCACCAAGAATCAGCATGCCCAGCAGTACATTGATCAGCGGGTTGATGTAATAGCCCAGGCTGGCCTGGat                                                           <  8:19590/149‑2 (MQ=255)
                                                                                            ccACTTGCCACACCTGCTGCGCCACACCCACGGCAGCCATGCCCACCGCCAGCCATTGCAGGCGGCGCAGGCGCTCACCAAGAATCAGCATGCCCAGCAGTACATTGATCAGCGGGTTGATGTAATAGCCCAGGCTGGCCTCGAGCATa                                                       >  3:163722/1‑149 (MQ=255)
                                                                                             cACTTGCCACACCTGCTGCGCCACACCCACGGCAGCCATGCCCACCGCCAGCCATTGCAGGCGGCGCAgg                                                                                                                                     >  7:23550/1‑70 (MQ=255)
                                                                                             cACTTGCCACACCTGCTGCGCCACACCCACGGCAGCCATGCCCACCGCCAGCCATTGCAGGCGGCGCAgg                                                                                                                                     <  8:23550/70‑1 (MQ=255)
                                                                                             cACTTGCCACACCTGCTGCGCCACACCCACGGCAGCCATGCCCACCGCCAGCCATTGCAGGCGGCGCAGGCGCTCACCAAGAATCAGCATGCCCAGCAGTACATTGATCAGCGGGTTGATGTAATAGCCCAGGCTGGc                                                                 <  8:151567/138‑1 (MQ=255)
                                                                                             cACTTGCCACACCTGCTGCGCCACACCCACGGCAGCCATGCCCACCGCCAGCCATTGCAGGCGGCGCAGGCGCTCACCAAGAATCAGCATGCCCAGCAGTACATTGATCAGCGGGTTGATGTAATAGCCCAGGCTGGc                                                                 >  7:151567/1‑138 (MQ=255)
                                                                                                  gCCACACCTGCTGCGCCACACCCACGGCAGCCATGCCCACCGCCAGCCATTGCAGGCGGCGCAGGCGCTCACCAAGAATCAGCATGCCCAGCAGTACATTGATCAGCGGGTTGATGTAATAGCCCAGGCTGGCCTCGAGCATACGGc                                                   <  6:96768/147‑1 (MQ=255)
                                                                                                    cacaCCTGCTGCGCCACACCCACGGCAGCCATGCCCACCGCCAGCCATTGCAGGCGGCGCAGGCGCTCACCAAGAATCAg                                                                                                                    >  4:114446/1‑80 (MQ=255)
                                                                                                    cacaCCTGCTGCGCCACACCCACGGCAGCCATGCCCACCGCCAGCCATTGCAGGCGGCGCAGGCGCTCACCAAGAATCAGCATGCCCAGCAGTACATTGATCAGCGGGTTGATGTAATAGCCCAGGCTGGCCTCGAGCATACGGCCgtt                                               <  8:85949/149‑1 (MQ=255)
                                                                                                     acacCTGCTGCGCCACACCCACGGCAGCCATGCCCACCGCCAGCCATTGCAGGCGGCGCAGGCGCTCACCAAGAATCAGCATGCCCAGCAGTACATTGATCAGCGGGTTGATGTAATAGCCCAGGCTGGCCTCGAGCATACGGCCgttg                                              >  5:64383/1‑149 (MQ=255)
                                                                                                                  cacaCCCACGGCAGCCATGCCCACCGCCAGCCATTGCAGGCGGCGCAGGCGCTCACCAAGAATCAGCATGCCCAGCAGTACATTGATCAGCGGGTTGATGTAATAGCCCAGGCTGGCCTCGAGCATACGGCCGTTGTTCACCGAc                                     >  4:164344/1‑145 (MQ=255)
                                                                                                                  cacaCCCACGGCAGCCATGCCCACCGCCAGCCATTGCAGGCGGCGCAGGCGCTCACCAAGAATCAGCATGCCCAGCAGTACATTGATCAGCGGGTTGATGTAATAGCCCAGGCTGGCCTCGAGCATACGGCCGTTGTTCACCGAc                                     <  3:164344/145‑1 (MQ=255)
                                                                                                                       ccACGGCAGCCATGCCCACCGCCAGCCATTGCAGGCGGCGCAGGCGCTCACCAAGAATCAGCATGCCCAGCAGTACATTGATCAGCGGGTTGATGTAATAGCCCAGGCTGGCCTCGAGCATACGGCCGTTGTTCAc                                         >  1:159840/1‑136 (MQ=255)
                                                                                                                       ccACGGCAGCCATGCCCACCGCCAGCCATTGCAGGCGGCGCAGGCGCTCACCAAGAATCAGCATGCCCAGCAGTACATTGATCAGCGGGTTGATGTAATAGCCCAGGCTGGCCTCGAGCATACGGCCGTTGTTCAc                                         <  2:159840/136‑1 (MQ=255)
                                                                                                                       ccACGGCAGCCATGCCCACCGCCAGCCATTGCAGGCGGCGCAGGCGCTCACCAAGAATCAGCATGCCCAGCAGTACATTGATCAGCGGGTTGATGTAATAGCCCAGGCTGGCCTCGAGCATACGGCCGTTGTTCACCGACCACACGTAg                            >  3:143512/1‑149 (MQ=255)
                                                                                                                        cACGGCAGCCATGCCCACCGCCAGCCATTGCAGGCGGCGCAGGCGCTCACCAAGAATCAGCATGCCCAGCAGTACATTGATCAGCGGGTTGAGGTAATAGCCCAGGCTGGCCTCGAGCATACGGCCGTTGTTCACCGACCACACGTAgg                           >  5:268134/1‑149 (MQ=255)
                                                                                                                               gCCATGCCCACCGCCAGCCATTGCAGGCGGCGCAGGCGCTCACCAAGAATCAGCATGCCCAGCAGTACATTGATCAGCGGGTTGATGTAATAGCCCAGGCTGGCCTCGAGCATACGGCCGTTGTTCACCGACCACACGTAGGTCAGc                      >  1:131619/1‑147 (MQ=255)
                                                                                                                               gCCATGCCCACCGCCAGCCATTGCAGGCGGCGCAGGCGCTCACCAAGAATCAGCATGCCCAGCAGTACATTGATCAGCGGGTTGATGTAATAGCCCAGGCTGGCCTCGAGCATACGGCCGTTGTTCACCGACCACACGTAGGTCAGc                      <  2:131619/147‑1 (MQ=255)
                                                                                                                                     cccACCGCCAGCCATTGCAGGCGGCGCAGGCGCTCACCAAGAATCAGCATGCCCAGCAGTACATTGATCAGCGGGTTGATGTAATAGCCCAGGCTGGCCTCTAGCATACGGCCGTTGTTCACAGACCACACGTAGGTCAGCCAGTTGcc              >  8:122502/1‑149 (MQ=255)
                                                                                                                                       cACCGCCAGCCATTGCAGGCGGCGCAGGCGCTCACCAAGAATCAGCATGCCCAGCAGTACATTGATCAGCGGGTTGATGTAATAGCCCAGGCTGGCCTCGAGCATACGGCCGTTGTTCACCGACCACACGTAGGTCAGCCAGTTGCCgg            >  6:38366/1‑149 (MQ=255)
                                                                                                                                       cACCGCCAGCCATTGCAGGCGGCGCAGGCGCTCACCAAGAATCAGCATGCCCAGCAGTACATTGATCAGCGGGTTGATGTAATAGCCCAGGCTGGCCTCGAGCATACGGCCGTTGTTCACCGACCACACGTAGGTCAGCCAGTTGCCgg            <  6:64383/149‑1 (MQ=255)
                                                                                                                                          cGCCAGCCATTGCAGGCGGCGCAGGCGCTCACCAAGAATCAGCATGCCCAGCAGTACATTGATCAGCGGGTTGATGTAATAGCCCAGGCTGGCCTCGAGCATACGGCCGTTGTTCACCGACCACACGTAGGTCAGCCAGTTGCCGGCGa         <  4:143512/149‑1 (MQ=255)
                                                                                                                                            ccagccaTTGCAGGCGGCGCAGGCGCTCACCAAGAATCAGCATGCCCAGCAGTACATTGATCAGCGGGTTGATGTAATAGCCCAGGCTGGCCTCGAGCATACGGCCGTTGTTCACCGACCACACGTAGGTCAGCCAGTTGCCGGCGAt        >  5:264694/1‑148 (MQ=255)
                                                                                                                                             cagccaTTGCAGGCGGCGCAGGCGCTCACCAAGAATCAGCATGTCCAGCAGTACATTGATCAGCGGGTTGATGTAATAGCCCAGGCTGGCCTCGAGCATACGGCCGTTGTTCACCGACCACACGTAGGTCAGCCAGTTGCCGGCGAt        >  4:275184/1‑147 (MQ=255)
                                                                                                                                             cagccaTTGCAGGCGGCGCAGGCGCTCACCAAGAATCAGCATGCCCAGCAGTACATTGATCAGCGGGTTGATGTAATAGCCCAGGCTGGCCTCGAGCATACGGCCGTTGTTCACAGACCACACGTAGGTCAGCCAGTTGCCGGCGATCa      <  7:122502/149‑1 (MQ=255)
                                                                                                                                              agccaTTGCAGGCGGCGCAGGCGCTCACCAAGAATCAGCATGCCCAGCAGTACATTGATCAGCGGGTTGATGTAATAGCCCAGGCTGGCCTCGAGCATACGGCCGTTGTTCACCGACCACACGTAGGTCAGCCAGTTGCCGGCGATCAg     <  5:38366/149‑1 (MQ=255)
                                                                                                                                                 caTTGCAGGCGGCGCAGGCGCTCACCAAGAATCAGCATGCCCAGCAGTACATTGATCAGCGGGTTGATGTAATAGCCCAGGCTGGCCTCGAGCATACGGCCGTTGTTCACCGACCACACGTAGGTCAGCCAGTTGCCGGCGATCAGTGc  <  2:280563/149‑1 (MQ=255)
                                                                                                                                                 caTTGCAGGCGGCGCAGGCGCTCACCAAGAATCAGCATGCCCAGCAGTACATTGATCAGCGGGTTGATGTAATAGCCCAGGCTGGCCTCGAGCATACGGCCGTTGTTCACCGACCACACGTAGGTCAGCCAGTTGCCGGCGATCAGTGc  >  1:280563/1‑149 (MQ=255)
                                                                                                                                                 |                                                                                                                                                    
GCCGGGCAGCGCAGCCACCGGCGCCTGCTTACGGATCAGCCCATAGAAGCCGAAGCTCAACGCCAGCACCAGTGATACCCACGGCAGGCTGCCCACTTGCCACACCTGCTGCGCCACACCCACGGCAGCCATGCCCACCGCCAGCCATTGCAGGCGGCGCAGGCGCTCACCAAGAATCAGCATGCCCAGCAGTACATTGATCAGCGGGTTGATGTAATAGCCCAGGCTGGCCTCGAGCATACGGCCGTTGTTCACCGACCACACGTAGGTCAGCCAGTTGCCGGCGATCAGTGC  >  NC_002947/436754‑437047

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