Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 3,839,252 G→A 9.1% G59S (GGC→AGC)  PP_3389 → membrane protein

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0029473,839,2520GA 9.1% 120.3 / 2.7 44G59S (GGC→AGC) PP_3389membrane protein
Reads supporting (aligned to +/- strand):  ref base G (19/21);  new base A (2/2);  total (21/23)
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.60e-01

TGGCAACTCAGTCAGAAGCTGCACGTAGAACGTACGAAGCAGGGATGACACCAGCAGAAAAGCGCGCAAATAGAGATTTTTATGAAGGCGACCTGAATGAAGAGATTCACTTTTCCCGGGTAAGTGAAACACCCCTCATACGCGGCCCTGTTTATGCATTCAATGAAACCCACATGGAAATGCGCTGTGGCGGCTCAGAAGAGAAGCGCGGCCTCATTACGCTAGTCACGATCGGCGGGGTGCTGCCTGGGATAATTTTCGGAACCGGCATAGCTGTCGGTTCCATCTGGCT  >  NC_002947/3839109‑3839400
                                                                                                                                               |                                                                                                                                                    
ttgCAACTCAGTCAGAAGCTGCAAGTAGAACGTACGAAGCAGGGATGAAACCAGCAGAAAAGCGCGCAAATAGAGATTTTTATGAAGGCGACCTGAATGAAGAGATTCACTTTTCCCGGGTAAGTGAAACACCCCTCATACGCGGccc                                                                                                                                                  <  7:230774/146‑1 (MQ=255)
      cTCAGTCAGAAGCTGCACGTAGAACGTACGAAGCAGGGATGACACCAGCAGAAAAGCGCGCAAATAGAGATTTTTATGAAGGCGACCTGAATGAAGAGATTCACTTTTCCCGGGTAAGTGAAACACCCCTCATACGCGGCCCTGTTTAt                                                                                                                                           <  2:93984/149‑1 (MQ=255)
      cTCAGTCAGAAGCTGCACGTAGAACGTACGAAGCAGGGATGACACCAGCAGAAAAGCGCGCAAATAGAGATTTTTATGAAGGCGACCTGAATGAAGAGATTCACTTTTCCCGGGTAAGTGAAACACCCCTCATACGCGGCCCTGTTTAt                                                                                                                                           <  1:138170/149‑1 (MQ=255)
      cTCAGTCAGAAGCTGCACGTAGAACGTACGAAGCAGGGATGACACCAGCAGAAAAGCGCGCAAATAGAGATTTTTATGAAGGCGACCTGAATGAAGAGATTCACTTTTCCCGGGTAAGTGAAACACCCCTCATACGCGGCCCTGTTTAt                                                                                                                                           <  3:45550/149‑1 (MQ=255)
         agtcagAAGCTGCACGTAGAACGTACGAAGCAGGGATGACACCAGCAGAAAAGCGCGCAAATAGAGATTTTTATGAAGGCGACCTGAATGAAGAGATTCACTTTTCCCGGGTAAGTGAAACACCCCTCATACGCGGCCCTGTTTATGc                                                                                                                                         <  4:97007/148‑1 (MQ=255)
               aaGCTGCACGTAGAACGTACGAAGCAGGGATGACACCAGCAGAAAAGCGCGCAAATAGAGATTTTTATGAAGGCGACCTGAATGAAGAGATTCACTTTTCCCGGGTAAGTGAAACACCCCTCATACGCGGCCCTGTTTATGCATTCAAt                                                                                                                                  >  7:178843/1‑149 (MQ=255)
                  cTGCAAGTAGAACGTACGAAGCAGGGATGACACCAGCAGAAAGGCGCGCAAATAGAGATTTTTATGAAGGCGACCTGAATGAAGAGATTCACTTTTCCCGGGTAAGTGAAACACCCCTCATACGCGGCCCTGTTTATGCATTCAATGa                                                                                                                                >  7:242734/1‑148 (MQ=255)
                       cGTAGAACGTACGAAGCAGGGATGACACCAGCAGAAAAGCGCGCAAATAGAGATTTTTATGAAGGCGACCTGAATGAAGAGATTCACTTTTCCCGGGTAAGTGAAACACCCCTCATACGCGGCCCTGTTTATGCATTCAATGAAACCca                                                                                                                          <  1:41677/149‑1 (MQ=255)
                         tAGAACGTACGAAGCAGGGATGACACCAGCAGAAAAGCGCGCAAATAGAGATTTTTATGAAGGCGACCTGAATGAAGAGATTCACTTTTCCCGGGTAAGTGAAACACCCCTCATACGCGGCCCTGTTTATGCATTCAATGAAACCcaca                                                                                                                        >  5:214001/1‑149 (MQ=255)
                         tAGAACGTACGAAGCAGGGATGACACCAGCAGAAAAGCGCGCAAATAGAGATTTTTATGAAGGCGACCTGAATGAAGAGATTCACTTTTCCCGGGTAAGTGAAACACCCCTCATACGCGGCCCTGTTTATGCATTCAATGAAACCcaca                                                                                                                        >  6:133396/1‑149 (MQ=255)
                              cGTACGAAGCAGGGATGACACCAGCAGAAAAGCGCGCAAATAGAGATTTTTATGAAGGCGACCTGAATGAAGAGATTCACTTTTCCCGGGTAAGTGAAACACCCCTCATACGCGGCCCTGTTTATGCATTCAATGAAACCCACATGGa                                                                                                                    >  3:85407/1‑148 (MQ=255)
                                         gggATGACACCAGCAGAAAAGCGCGCATATAGAGATTTTTATGAAGTCGACCTTAATGAAGAGATTCACTTTTACCTGGTAAGTGAAACACCCCTCATACGCGGCCCTGTTTATGCATTCAATTAAACCCACATGGAAATGCTCTGTgg                                                                                                        <  6:121006/149‑1 (MQ=255)
                                          ggATGACACCAGCAGAAATGCGCGCAAATAGAGATTTTTATGAAGGCGACCTGAATGAAGAGATTCACTTTTCCCGGGTAAGTGAAACACCCCTCATACGCGGCCCTGTTTATGCATTCAATGAAACCCACATGGAAATGCGCTGTggc                                                                                                       <  8:242734/149‑1 (MQ=255)
                                          ggATGACACCAGCAGAAAAGCGCGCAAATAGAGATTTTTATGAAGGCGACCTGAATGAAGAGATTCACTTTTCCCGGGTAAGTGAAACACCCCTCATACGCAGCCCTGTTTATGCATTCAATGAAACCCACATGGAAATGCGCTGTggc                                                                                                       >  2:72301/1‑149 (MQ=255)
                                          ggATGACACCAGCAGAAAAGCGCGCAAATAGAGATTTTTATGAAGGCGACCTGAATGAAGAGATTCACTTTTCCCGGGTAAGTGAAACACCCCTCATACGCAGCCCTGTTTATGCATTCAATGAAACCCACATGGAAATGCGCTGTggc                                                                                                       >  2:68809/1‑149 (MQ=255)
                                           gATGACACCAGCAGAAAAGCGCGCAAATAGAGATTTTTATGAAGGCGACCTGAATGAAGAGATTCACTTTTCCCGGGTAAGTGAAACACCCCTCATACGCGGCCCTGTTTATGCATTCAATGAAACCCACATGGAAATGCGCTGTggcg                                                                                                      <  5:223436/149‑1 (MQ=255)
                                               acacCAGCAGAAAAGCGCGCAAATAGAGATTTTTATGAAGGCGACCTGAATGAAGAGATTCACTTTTCCCGGGTAAGTGAAACACCCCTCATACGCAGCCCTGTTTATGCATTCAATGAAACCCACATGGAAATGCGCTGTGGCGGCTc                                                                                                  <  1:68809/149‑1 (MQ=255)
                                                cacCAGCAGAAAAGCGCGCAAATAGAGATTTTTATGAAGGCGACCTGAATGAAGAGATTCACTTTTCCCGGGTAAGTGAAACACCCCTCATACGCGGCCCTGTTTATGCATTCAATGAAACCCACATGGAAATGCGCTGTggcgg                                                                                                     >  3:40210/1‑145 (MQ=255)
                                                cacCAGCAGAAAAGCGCGCAAATAGAGATTTTTATGAAGGCGACCTGAATGAAGAGATTCACTTTTCCCGGGTAAGTGAAACACCCCTCATACGCGGCCCTGTTTATGCATTCAATGAAACCCACATGGAAATGCGCTGTggcgg                                                                                                     <  4:40210/145‑1 (MQ=255)
                                                cacCAGCAGAAAAGCGCGCAAATAGAGATTTTTATGAAGGCGACCTGAATGAAGAGATTCACTTTTCCCGGGTAAGTGAAACACCCCTCATACGCGGCCCTGTTTATGCATTCAATGAAACCCACATGGAAATGCGCTGTGGCGGCTCa                                                                                                 <  6:214001/149‑1 (MQ=255)
                                                cacCAGCAGAAAAGCGCGCAAATAGAGATTTTTATGAAGGCGACCTGAATGAAGAGATTCACTTTTCCCGGGTAAGTGAAACACCCCTCATACGCAGCCCTGTTTATGCATTCAATGAAACCCACATGGAAATGCGCTGTGGCGGCTc                                                                                                  <  1:72301/148‑1 (MQ=255)
                                                  ccAGCAGAAAAGCGCGCAAATAGAGATTTTTATGAAGGCGACCTGAATGAAGAGATTCACTTTTCCCGGGTAAGTGAAACACCCCTCATACGCGGCCCTGTTTATGCATTCAATGAAACCCACATGGAAATGCGCTGTGGCGGCTCaga                                                                                               <  4:22543/149‑1 (MQ=255)
                                                    agcatAAAATCGCGCAAATAGAGATTTTTATGAAGGCGACCTTAATGAAGAGATTCACTTTTCCCGGGTAAGTGAAACACCGCTCATACGCGGCCCTGTTTATGCATTCAATGAAACCCACATGGAAATGCGCTGTGGCGGCTCagaag                                                                                             <  4:85407/149‑1 (MQ=255)
                                                                   aaaTAGAGATTTTTATGAAGGCGACCTGAATGAAGAGATTCACTTTTCCCGGGTAAGTGAAACACCCCTCATACGCGGCCCTGTTTATGCATTCAATGAAACCCACATGGAAATGCGCTGTGGCGGCTCAGAAGAGAAGCGCGGCCTCa                                                                              >  6:54867/1‑149 (MQ=255)
                                                                                     aGGCGACCTGAATGAAGAGATTCACTTTTCCCGGGTAAGTGAAACACCCCTCATACGCGGCCCTGTTTATGCATTCAATGAAACCCACATGGAAATGCGCTGTGGCGGCTCAGAAGAGAAGCGCGGCCTCAtt                                                                            <  3:159656/133‑1 (MQ=255)
                                                                                     aGGCGACCTGAATGAAGAGATTCACTTTTCCCGGGTAAGTGAAACACCCCTCATACGCGGCCCTGTTTATGCATTCAATGAAACCCAAATGGAAATGCGCTGTGGCGGCTCAGAAGACAAGCGCGGCCTCaat                                                                            >  4:159656/1‑131 (MQ=255)
                                                                                        cGACCTGACTGAAGTGATTCACTTTTCCCGGGTTAGTGAAACACCCCTCATACGCGGCCCTGTTTATGCATTCAATGAAACCCACATGGAAATGCGCTGTGGCGGCTCAGAAGAGAAGCGCGGCCTCATTACGCTAGTCACGATcggcg                                                         <  8:222931/149‑1 (MQ=255)
                                                                                            cTGAATGAAGAGATTCACTTTTCCCGGGTAAGTGAAACACCCCTCATACGCGGCCCTGTTTATGCATTCAATGAAACCCACATGGAAATGCGCTGTGGCGGCTCAGAAGAGAAGCGCGGCCTCATTACGCTAGTCACGATCGGCGGGGt                                                     <  5:133396/149‑1 (MQ=255)
                                                                                              gaatgaaGAGATTCACTTTTCCCGGGTAAGTGAAACACCCCTCATACGCGGCCCTGTTTATGCATTCAATGAAACCCACATGGAAATGCGCTGTGGCGGCTCAGAAGAGAAGCGCGGCCTCATTACGCTAGTCACGATCGGCGGGGtgc                                                   >  1:43440/1‑149 (MQ=255)
                                                                                              gaatgaaGAGATTCACTTTTCCCGGGTAAGTGAAACACCCCTCATACGCGGCCCTGTTTATGCATTCAATGAAACCCACATGGAAATGCGCTGTGGCGGCTCAGAAGAGAAGCGCGGCCTCATTACGCTAGTCACGATCGGCGGGGtgc                                                   >  6:80080/1‑149 (MQ=255)
                                                                                                  gaaGAGATTCACTTTTCCCGGGTAAGTGAAACACCCCTCATACGCGGCCCTGTTTATGCATTCAATGAAACCCACATGGAAATGCGCTGTGGCGGCTCAGAAGAGAAGCGCGGCCTCATTACGCTAGTCACGATCGGCGGGGTGCTGcc                                               >  5:169223/1‑149 (MQ=255)
                                                                                                    agagATTCACTTTTCCCGTGTAAGTGAAACACCCCTCATACGCGGCCCTGTTTATGCATTCAATGAAACCCACATGGAAATGCGCTGTGGCGGCTCAGAAGAGAAGCGCGGCCTCATTACGCTAGTCACGATCGGCGGGGTGCTGCCTg                                             <  8:178843/149‑1 (MQ=255)
                                                                                                        aTTCACTTTTCCCGGGTAAGTGAAACACCCCTCATACGCGGCCCTGTTTATGCATTCAATGAAACCCACATGGAAATGCGCTGTGGCGGCTCAGAAGAGAAGCGCGGCCTCATTACGCTAGTCACGATcggcgg                                                        >  5:113471/1‑134 (MQ=255)
                                                                                                        aTTCACTTTTCCCGGGTAAGTGAAACACCCCTCATACGCGGCCCTGTTTATGCATTCAATGAAACCCACATGGAAATGCGCTGTGGCGGCTCAGAAGAGAAGCGCGGCCTCATTACGCTAGTCACGATcggcgg                                                        <  6:113471/134‑1 (MQ=255)
                                                                                                        aTTCACTTTTCCCGGGTAAGTGAAACACCCCTCATACGCGGCCCTGTTTATGCATTCAATGAAACCCACATGGAAATGCGCTGTGGCGGCTCAGAAGAGAAGCGCGGCCTCATTACGCTAGTCACGATCGGCGGGGTGCTGCCTGGGAt                                         <  3:13768/149‑1 (MQ=255)
                                                                                                         ttCACTTTTCCCGGGTAAGTGAAACACCCCTCATACGCGGCCCTGTTTATGCATTCAATGAAACCCACATGGAAATGCGCTGTGGCGGCTCAGAAGAGAAGCGCGGCCTCATTACGCTAGTCACGATCGGCGGGGTGCTGCCTGGGAt                                         >  7:49351/1‑148 (MQ=255)
                                                                                                          tCACTTTTCCCGGGTAAGTGAAACACCCCTCATACGCGGCCCTGTTTATGCATTCAATGAAACCCACATGGAAATGCGCTGTGGCGGCTCAGAAGAGAAGCGCGGCCTCATTACGCTAGTCACGATCGGCGGGGTGCTGCCTGGGATaa                                       >  1:264476/1‑149 (MQ=255)
                                                                                                             cTTTTCCCGGGTAAGTGAAACACCCCTCATACGCGGCCCTGTTTATGCATTCAATGAAACCCACATGGAAATGCGCTGTGGCGGCTCAGAAGAGAAGCGCGGCCTCATTACGCTAGTCACGATCGGCGGGGTGCTGCCTGGGATAAttt                                    >  7:26074/1‑149 (MQ=255)
                                                                                                              ttttCCCGGGTAAGTGAAACACCCCTCATACGCGGCCCTGTTTATGCATTCAATGAAACCCACATGGAAATGCGCTGTGGCGGCTCAGAAGAGAAGCGCGGCCTCATTACGCTAGTCACGATCGGCGGGGTGCTGCCTGGGATAAtttt                                   >  7:208141/1‑149 (MQ=255)
                                                                                                              ttttCCCGGGTAAGTGAAACACCCCTCATACGCGGCCCTGTTTATGCATTCAATGAAACCCACATGGAAATGCGCTGTGGCGGCTCAGAAGAGAAGCGCGGCCTCATTACGCTAGTCACGATCGGCGGGGTGCTGCCTGGGATAAtttt                                   >  1:262768/1‑149 (MQ=255)
                                                                                                                           gTGAAACACCCCTCATACGCGGCCCTGTTTATGCATTCAATGAAACCCACATGGAAATGCGCTGTGGCGGCTCAGAAGAGAAGCGCGGCCTCATTACGCTAGTCACGATCGGCGGGGTGCTGCCTGGGATAATTTTCg                                 <  3:30608/138‑1 (MQ=255)
                                                                                                                           gTGAAACACCCCTCATACGCGGCCCTGTTTATGCATTCAATGAAACCCACATGGAAATGCGCTGTGGCGGCTCAGAAGAGAAGCGCGGCCTCATTACGCTAGTCACGATCGGCGGGGTGCTGCCTGGGATAATTTTCg                                 >  4:30608/1‑138 (MQ=255)
                                                                                                                                  acCCCTCATACGCGGCCCTGTTTATGCATTCAATGAAACCCACATGGAAATGCGCTGTGGCGGCTCAGAAGAGAAGCGCGGCCTCATTACGCTAGTCACGATCGGCGGGGTGCTGCCTGGGATAATTTTCGGAACCGGCATAGCTGTCg               <  2:200240/149‑1 (MQ=255)
                                                                                                                                   ccccTCATACGCGGCCCTGTTTATGCATTCAATGAAACCCACATGGAAATGCGCTGTGGCGGCTCAGAAGAGAAGCGCGGCCTCATTACGCTAGTCACGATCGGCGGGGTGCTGCCTGGGATAATTTTCGGAACCGGCATAGCTGTCgg              >  1:35094/1‑149 (MQ=255)
                                                                                                                                               ggCCCTGTTTATGCATTCAATGAAACCCACATGGAAATGCGCTGTGGCGGCTCAGAAGAGAAGCGCGGCCTCATTACGCTAGTCACGATCGGCGGGGTGCTGCCTGGGATAATTTTCGGAACCGGCATAGCTGTCGGTTCCATCTGGCt  <  8:208141/149‑1 (MQ=255)
                                                                                                                                               |                                                                                                                                                    
TGGCAACTCAGTCAGAAGCTGCACGTAGAACGTACGAAGCAGGGATGACACCAGCAGAAAAGCGCGCAAATAGAGATTTTTATGAAGGCGACCTGAATGAAGAGATTCACTTTTCCCGGGTAAGTGAAACACCCCTCATACGCGGCCCTGTTTATGCATTCAATGAAACCCACATGGAAATGCGCTGTGGCGGCTCAGAAGAGAAGCGCGGCCTCATTACGCTAGTCACGATCGGCGGGGTGCTGCCTGGGATAATTTTCGGAACCGGCATAGCTGTCGGTTCCATCTGGCT  >  NC_002947/3839109‑3839400

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: