Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 2,777,050 C→G 8.1% W166C (TGG→TGC)  PP_2429 ← membrane protein

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0029472,777,0500CG 8.1% 148.7 / 3.7 49W166C (TGG→TGC) PP_2429membrane protein
Reads supporting (aligned to +/- strand):  ref base C (20/25);  new base G (2/2);  total (22/27)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.68e-01

GAGGCGTTGACACGTTGAGCCGAGGCCGCTCAGGCCCCTTGTCGGTCAAGTAGAGGCTCAGGCCAATACTGCAGCCCATGCCGATATCAACAAAACCCCCGCCAACACTTGGTTGAACACCCGCAACCGCCCCGGTGCCTGCAGCCACC‑GGGCGCTCCCCACCCCCAACACCGCCCAGGCGGCCATGCTCGGCAGAGCGATCAGCAAGAACACCAGCGCCAGCTGCCACACTGGCAACGACGGCGCGGCAAATACCCCGACCACCGCCACGGCCATCAACCAGACCTTGGGGT  >  NC_002947/2776906‑2777198
                                                                                                                                                |                                                                                                                                                     
gAGGCGTTGACACGTTGAGCCGAGGCCGCTCAGGCCCCTTGTCGGTCAAGTAGAGGCTCAGGCCAATACTGCAGCCCATGCCGATATCAACAAAACCCCCGCCAACACTTGGTTGAACACCCGCAACCGCCCCGGTGCCTGCAGCCACC‑g                                                                                                                                                 >  1:526824/1‑150 (MQ=255)
 aGGCGTTGACACGTTGAGCCGAGGCCGCTCAGGCCCCTTGTCGGTCAAGTAGAGGCTCAGGCCAATACTGCAGCCCATGCCGATATCAACAAAACCCCCGCCAACACTTGGTTGAACACCCGCAACCGCCCCGGTGCCTGCAGCCACC‑gg                                                                                                                                                <  2:146934/150‑1 (MQ=255)
        gACACGTTGAGCCGAGGCCGCTCAGGCCCCTTGTCGGTCAAGTAGAGGCTCAGGCCAATACTGCAGCCCATGCCGATATCAACAAAACCCCCGCCAACACTTGGTTGAACACCCGCAACCG‑CCCGGTGGCTGCAGGCACC‑ggggcggttg                                                                                                                                        <  1:572030/150‑8 (MQ=255)
           acGTTGAGCCGAGGCCGCTCAGGCCCCTTGTCGGTCAAGTAGAGGCTCAGGCCAATACTGCAGCCCATGCCGATATCAACAAAACCCCCGCCAACACTTGGTTGAACACCCGCAACCGCCCCGGTGCCTGCAGCCACC‑GGGCGCTCCCCa                                                                                                                                      >  1:286081/1‑150 (MQ=255)
           acGTTGAGCCGAGGCCGCTCAGGCCCCTTGTCGATATTTTAGAGGCTCAGGCCAATACTGCAGCCCATGCCAATATAAGAAAAACCCCCGCCAACACTTGGTTGAACACCCGCAACCGCCCCGGTGCCTGCAGCCACC‑GGGCGCTCCCCa                                                                                                                                      >  1:462028/1‑150 (MQ=255)
             gTTGAGCCGAGGCCGCTCAGGCCCCTTGTCGGTCAAGTAGAGGCTCAGGCCAATACTGCAGCCCATGCCGATATCAACAAAACCCCCGCCAACACTTGGTTGAACACCCGCAACCGCCCCGGTGCCTGCAGCCACC‑GGGCGCTCCCCAcc                                                                                                                                    >  2:509935/1‑150 (MQ=255)
                                ggCCCCTTGTCTGTCAAGTAGAGGCTCAGGCCAATACTGCAGCCCATGCCGATATCAACAAAACCCCCGCCAACACTTGGTTGAACACCCGCAACCGCCCCGGTGCCTGCAGCCACC‑GGGCGCTCCCCACCCCCAACACCGCCCAggcgg                                                                                                                 <  1:797297/150‑1 (MQ=255)
                                 gCCCCTTGTCGGTCAAGTAGAGGCTCAGGCCAATACTGCAGCCCATGCCGATATCAACAAAACCCCCGCCAACACTTGGTTGAACACCCGCAACCGCCCCGGTGCCTGCAGCCACC‑gg                                                                                                                                                >  2:92382/1‑118 (MQ=255)
                                 gCCCCTTGTCGGTCAAGTAGAGGCTCAGGCCAATACTGCAGCCCATGCCGATATCAACAAAACCCCCGCCAACACTTGGTTGAACACCCGCAACCGCCCCGGTGCCTGCAGCCACC‑gg                                                                                                                                                <  1:92382/118‑1 (MQ=255)
                                           ggTCAAGTAGAGGCTCAGGCCAATACTGCAGCCCATGCCGATATCAACAAAACCCCCGCCAACACTTGGTTGAACACCCGCAACCGCCCCGGTGCCTGCAGCCACC‑GGGCGCTCCCCACCCCCAACACCGCCCAGGCGGCCATGCTCGGc                                                                                                      >  1:92450/1‑150 (MQ=255)
                                               aaGTAGAGGCTCAGGCCAATACTGCAGCCCATGCCGATATCAACAAAACCCCCGCCAACACTTGGTTGAACACCAGCAACCGCCCCGGTGCCTGCAGCCACC‑GGGCGCTCCCCACCCCCAACACCGCCCAGGCGGCCATGCTCGGCagag                                                                                                  >  1:233892/1‑150 (MQ=255)
                                                  tAGAGGCTCAGGCCAATACTGCAGCCCATGCCGATATCAACAAAACCCCCGCCAACACTTGGTTGAACACCCGCAACCGCCCCGGTGCCTGCAGCCACC‑GGgcgc                                                                                                                                            >  1:162984/1‑105 (MQ=255)
                                                  tAGAGGCTCAGGCCAATACTGCAGCCCATGCCGATATCAACAAAACCCCCGCCAACACTTGGTTGAACACCCGCAACCGCCCCGGTGCCTGCAGCCACC‑GGgcgc                                                                                                                                            <  2:162984/105‑1 (MQ=255)
                                                   agagGCTCAGGCCAATACTGCAGCCCATGCCGATATCAACAAAACCCCCGCCAACACTTGGTTGAACACCCGCAACCGCCCCGGTGCCTGCAGCCACC‑GGGCGCTCCCCACCCCCAACACCGCCCAGGCGGCCATGCTCGGCAGAGCGAt                                                                                              <  1:197910/150‑1 (MQ=255)
                                                     agGCTCAGGCCAATACTGCAGCCCATGCCGATATCTACACATCCCCCGCCAACACTTGGTTGAACACCCGCAACCGCCCCGGTGCCTGCAGCCACC‑GGGCGCTCCCCACCCCCAACACCGCCCAGGCGGCCATGCTCGGCAGAGCGATCa                                                                                            <  2:462028/150‑1 (MQ=255)
                                                         tCAGGCCAATACTGCAGCCCATGCCGATATCAACAAAACCCCCGCCAACACTTGGTTGAACACCCGCAACCGCCCCGGTGCCTGCAGCCACC‑GGGCGCTCCCCACCCCCAACACCGCCCAGGCGGCCATGCTCGGCAGAGCGATCAGCaa                                                                                        <  2:278341/150‑1 (MQ=255)
                                                         tCAGGCCAATACTGCAGCCCATGCCGATATCAACAAAACCCCCGCCAACACTTGGTTGAACACCCGCAACCGCCCCGGTGCCTGCAGCCACC‑GGGCGCTCCCCACCCCAAACACCGCCCAGGCGGCCATGCTCGGCAGAGCGATCAGCaa                                                                                        <  2:286081/150‑1 (MQ=255)
                                                             gCCAATACTGCAGCCCATGCCGATATCAACAAAACCCCCGCCAACACTTGGTTGAACACCCGCAACCGCCCCGGTGCCTGCAGCCACC‑GGGCGCTCCCCACCCCCAACACCGCCCAGGCGGcc                                                                                                               <  2:565249/123‑1 (MQ=255)
                                                             gCCAATACTGCAGCCCATGCCGATATCAACAAAACCCCCGCCAACACTTGGTTGAACACCCGCAACCGCCCCGGTGCCTGCAGCCACC‑GGGCGCTCCCCACCCCCAACACCGCCCAGGCGGcc                                                                                                               >  1:565249/1‑123 (MQ=255)
                                                              ccAATACTGCAGCCCATGCCGATATCAACAAAACCCCCGCCAACACTTGGTTGAACACCCGCAACCGCCCCGGTGCCTGCAGCCACC‑GGGCGCTCCCCACCCCCAACACCGCCCAGGCGGCCATGCTCGGCAGAGCGATCAGCAAGAaca                                                                                   <  1:502123/150‑1 (MQ=255)
                                                                aaTACTGCAGCCCATGCCGATATCAACAAAACCCCCGCCAACACTTGGTTGAACACCCGCAACCGCCCCGGTGCCTGCAGCCACC‑GGGCGCTCCCCACCCCCAACACCGCCCAGGCGGCCATGCTCGGCAGAGCGATCAGCAAGAACAcc                                                                                 >  1:455924/1‑150 (MQ=255)
                                                                   aCTGCAGCCCATGCCGATATCAACAAAACCCCCGCCAACACTTGGTTGAACACCCGCAACCGCCCCGGTGCCTGCAGCCACC‑GGGCGCTCCCCACCCCCAACACCGCCCAGGCGGCCATGCTCGGCAGAGCGATCAGCAAGAACACCAgc                                                                              <  1:509935/150‑1 (MQ=255)
                                                                          cccATGCCGATATCAACAAAACCCCCGCCAACACTTGGTTGAACACCCGCAACCGCCCCGGTGCCTGCAGCCACC‑GGGCGCTCCCCACCCCCAACACCGCCCAGGCGGCCATGCTCGGCAGAGCGATCAGCAAGAACACCAGCGCCAGCt                                                                       >  2:798060/1‑150 (MQ=255)
                                                                              agCCGATATCAACAAAACCCCCGCCAACACTTGGTTGAACACCAGCAACCGCCCCGGTGCCTGCAGCCACC‑GGGCGCTCCCCACCCCCAACACCGCCCAGGCGGCCATGCTCGGCAGAGCGATCAGCAAGAACACCAGCGCCAGCTGCca                                                                   <  2:233892/149‑1 (MQ=255)
                                                                                 cGATATCAACAAAACCCCCGCCAACACTTGGTTGAACACCCGCAACCGCCCCGGTGCCTGCAGCCACC‑GGGCGCTCCCCACCCCCAACACCGCCCAGGCGGCCATGCTCGGCAGAGCGATCAGCAAGAACACCAGCGCCAGCTGCCacac                                                                >  2:531458/1‑150 (MQ=255)
                                                                                      tCAACAAAACCCCCGCCAACACTTGGTTGAACACCCGCAACCGCCCCGGTGCCTGCAGCCACC‑GGGCGCTCCCCACCCCCAACACCGCCCAggcgg                                                                                                                 <  2:194436/96‑1 (MQ=255)
                                                                                      tCAACAAAACCCCCGCCAACACTTGGTTGAACACCCGCAACCGCCCCGGTGCCTGCAGCCACC‑GGGCGCTCCCCACCCCCAACACCGCCCAggcgg                                                                                                                 >  1:194436/1‑96 (MQ=255)
                                                                                        aacaaAACCCCCGCCAACACTTGGTTGAACACCCGCAACCGCCCCGGTGCCTGCAGCCACC‑GGGCGCTCCCCACCCCCAACACCGCCCAGGCGGCCAt                                                                                                             <  1:54662/98‑1 (MQ=255)
                                                                                        aacaaAACCCCCGCCAACACTTGGTTGAACACCCGCAACCGCCCCGGTGCCTGCAGCCACC‑GGGCGCTCCCCACCCCCAACACCGCCCAGGCGGCCAt                                                                                                             >  2:54662/1‑98 (MQ=255)
                                                                                        aacaaAACCCCCGCCAACACTTGGTTGAACACCCGCAACCGCCCCGGTGCCTGCAGCCACC‑GGGCGCTCCCCACCCCCAACACCGCCCAGGCGGCCATGCTCTGCCGAGcatggc                                                                                            >  1:427714/1‑110 (MQ=255)
                                                                                        aacaaAACCCCCGCCAACACTTGGTTGAACACCCGCAACCGCCCCGGTGCCTGCAGCCACC‑GGGCGCTCCCCACCCCCAACACCGCCCAGGCGGCCATGCTCTGCCGAGcatggc                                                                                            <  2:427714/115‑6 (MQ=255)
                                                                                           aaaaCCCCCGCCAACACTTGGTTGAACACCCGCAACCGCCCCGGTGCCTGCAGCCACC‑GGGCGCTCCCCACCCCCAACACCGCCCAGGCGGCCATGCTCGGCAGAGCGATCAGCAAGAACACCAGCGCCAGCTGCCACACTGGCAacgac                                                      >  2:273131/1‑150 (MQ=255)
                                                                                                 cccGCCAACACTTGGTTGAACACCCGCAACCGCCCCGGTGCCTGCAGCCACC‑GGGCGCTCCCCACCCCCAACACCGCCCAGGCGGCCATGCTCGGCAGAGCGATCAGCAAGAACACCAGCGCCAGCTGCCACACTGGCAACGACGGcgcg                                                >  1:205593/1‑150 (MQ=255)
                                                                                                      cAACACTTGGTTGAACACCCGCAACCGCCCCGGTGCCTGCAGCCACC‑GGGCGCTCCCCACCCCCAACACCGCCCAGGCGGCCATGCTCGGCAGAGCGATCAGCAAGAACACCAGCGCCAGCTGCCACACTGGCAACGACGGCGCGGCaaa                                           <  2:205593/150‑1 (MQ=255)
                                                                                                           cTTGGTTGAACACCCGCAACCGCCCCGGTGCCTGCAGCCACC‑GGGCGCTCCCCACCCCCAACACCGCCCAgg                                                                                                                    >  2:341200/1‑72 (MQ=255)
                                                                                                           cTTGGTTGAACACCCGCAACCGCCCCGGTGCCTGCAGCCACC‑GGGCGCTCCCCACCCCCAACACCGCCCAgg                                                                                                                    <  1:341200/72‑1 (MQ=255)
                                                                                                            ttGGTTGAACACCCGCAACCGCCCCGGTGCCTGCAGCCACC‑GGGCGCTCCCCACCCCCaa                                                                                                                               <  1:724035/60‑1 (MQ=255)
                                                                                                            ttGGTTGAACACCCGCAACCGCCCCGGTGCCTGCAGCCACC‑GGGCGCTCCCCACCCCCaa                                                                                                                               >  2:724035/1‑60 (MQ=255)
                                                                                                            ttGGTTGAACACCCGCAACCGCCCCGGTGCCTGCAGCCACC‑GGGCGCTCCCCACCCCCAACACCGCCCAGGCGGCCATGCTCGGCAGAGCGATCAGCAAGAACACCAGCGCCAGCTGCCACACTGGCAacga                                                       <  1:88587/132‑1 (MQ=255)
                                                                                                            ttGGTTGAACACCCGCAACCGCCCCGGTGCCTGCAGCCACC‑GGGCGCTCCCCACCCCCAACACCGCCCAGGCGGCCATGCTCGGCAGAGCGATCAGCAAGAACACCAGCGCCAGCTGCCACACTGGCAacga                                                       >  2:88587/1‑132 (MQ=255)
                                                                                                               gTTGAACACCCGCAACCGCCCCGGTGCCTGCAGCCACC‑GGGCGCTCCCCACCCCCAACACCGCCCAGGCGGCCATGCTCGGCAGAGCGATCAGCAAGAACACCAGCGCCAGCTGCCACACTGTCAATGACGGCGCGGCAAATACCCCGac                                  <  1:419936/150‑1 (MQ=255)
                                                                                                               gTTGAACACCCGCAACCGCCCCGGTGCCTGCAGCCACC‑GGGCGCTCCCCACCCCCAACACCGCCCAGGCGGCCATGCTCGGCAGAGCGATCAGCAAGAACACCAGCGCCAGCTGCCACACTGGCAACGACGGCGCGGCAAATACCCCGac                                  <  2:526824/150‑1 (MQ=255)
                                                                                                                                 gcccGGTGGCTGCAGGCACCGGGGCGCTCCCCACCCCCAACACCGCCCAGGCGGCCATGCTCGGCAGAGCGATCAGCAAGAACACCAGCGCCAGCTGCCACACTGGCAACGACGGCGCGGCAAATACCCCGACCACCGCCACGGCCATCa                 >  2:143709/2‑150 (MQ=255)
                                                                                                                                 ccccGGTGCCTACAGCCACC‑GGGCGCTCCCCACCCCCAACACCGTCAAGGCGGCCATGCTCGGCAGAGCGATCAGCAAGAACACCAGCACCAGCTGCCACACTGGCAACGACGACGCGGCAAATACCCCCACCACCGCCACGGCCATCaa                <  2:77557/150‑1 (MQ=255)
                                                                                                                                  cccGGTGGCTGCAGGCACC‑GGGCGCTCCCCACCCCCAACACCGCCCAGGCGGCCATGCTCGGCAGAGCGATCAGCAAGAACACCAGCGCCAGCTGCCACACTGGCAACGACGGCGCGGc                                              <  2:351288/119‑1 (MQ=255)
                                                                                                                                  cccGGTGGCTGCAGGCACC‑GGGCGCTCCCCACCCCCAACACCGCCCAGGCGGCCATGCTCGGCAGAGCGATCAGCAAGAACACCAGCGCCAGCTGCCACACTGGCAACGACGGCGCGGc                                              >  1:351288/1‑119 (MQ=255)
                                                                                                                                  cccGGTGCCTGCAGCCACC‑GGGCGCGCCCCACCCCCAACACCGCCCAGGCGGCCATGCTCGGCAGAGCGATCAGCAAGAACACCAGCGCCAGCTGCCACACTGGCAACGACGGCGCGGCAAATACCCCGACCACCGCCACGGCCATCAAc               <  1:491792/150‑1 (MQ=255)
                                                                                                                                             cAGCCACC‑GGGCGCTCCCCACCCCCAACACCGCCCAGGCGGCCATGCTCGGCAGAGCGATCAGCAAGAACACCAGCGCCAGCTGCCACACTGGCAACGACGGCGCGGCAAATACCCCGACCACCGCCACGGCCATCAACCAGACCTTggg    <  2:92450/150‑1 (MQ=255)
                                                                                                                                               gCCACC‑GGGCGCTCCCCACCCCCAACACCGCCCAGGCGGCCATGCTCGGCAGAGCGATCAGCAAGAACACCAGCGCCAGCTGCCACACTGGCAACGACGGCGCGGCAAATACCCCGACCACCGCCACGGCCATCAACCAGACCTTGGGGt  <  2:455924/150‑1 (MQ=255)
                                                                                                                                                |                                                                                                                                                     
GAGGCGTTGACACGTTGAGCCGAGGCCGCTCAGGCCCCTTGTCGGTCAAGTAGAGGCTCAGGCCAATACTGCAGCCCATGCCGATATCAACAAAACCCCCGCCAACACTTGGTTGAACACCCGCAACCGCCCCGGTGCCTGCAGCCACC‑GGGCGCTCCCCACCCCCAACACCGCCCAGGCGGCCATGCTCGGCAGAGCGATCAGCAAGAACACCAGCGCCAGCTGCCACACTGGCAACGACGGCGCGGCAAATACCCCGACCACCGCCACGGCCATCAACCAGACCTTGGGGT  >  NC_002947/2776906‑2777198

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