Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 2,777,007 C→G 9.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,0070CG 9.1% 117.1 / 4.6 44W166C (TGG→TGC) PP_2429membrane protein
Reads supporting (aligned to +/- strand):  ref base C (24/16);  new base G (2/2);  total (26/18)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 2.69e-01

TCGCGAGGCGTTGACACGTTGAGCCGAGGCCGCTCAGGCCCCTTGTCGGTCAAGTAGAGGCTCAGGCCAATACTGCAGCCCATGCCGATATCAACAAAACCCCCGCCAACACTTGGTTGAACACCCGCAACCGCCCCGGTGCCTGCAGCCACCGGGCGCTCCCCACCCCCAACACCGCCCAGGCGGCCATGCTCGGCAGAGCGATCAGCAAGAACACCAGCGCCAGCTGCCACACTGGCAACGACGGCGCGGCAAATACCCCGACCACCGCCACGGCCATCAACCAGACCTT  >  NC_002947/2776859‑2777150
                                                                                                                                                    |                                                                                                                                               
tCGCGAGGCGTTGACACGTTGAGCCGAGGCCGCTCAGGCCCCTTGTCGGTCAAGTAGAGGCTCAGGCCAATACTGCAGCCCATGCCGATATCAACAAAACCCCCGCCAACACTTGGTTGAACACCCGCAACCGCCCCGGTGCCTGCAGc                                                                                                                                                 >  4:253083/1‑149 (MQ=255)
    gAGGCGTTGACACGTTGAGCCGAGGCCGCTCAGGCCCCTTGTCGGTCAAGTAGAGGCTCAGGCCAATACTGCAGCCCATGCCGATATCAACAAAACCCCCGCCAACACTTGGTTGAACACCCGCAACCGCCCCGGTGCCTGCAGCCAcc                                                                                                                                             >  7:274477/1‑149 (MQ=255)
        cGTTGACACGTTGAGCCGAGGCCGCTCAGGCCCCTTGTCGGTCAAGTAGAGGCTCAGGCCAATACTGCAGCCCATGCCGATATCAACAAAACCCCCGCCAACACTTGGTTGAACACCCGCAACCGCCCCGGTGCCTGCAGCCACCGGgc                                                                                                                                         <  5:192812/149‑1 (MQ=255)
         gTTGACACGTTGAGCCGAGGCCGCTCAGGCCCCTTGTCGGTCAAGTAGAGGCTCAGGCCAATACTGCAGCCCATGCCGATATCAACAAAACCCCCGCCAACACTTGGTTGAACACCCGCAACCGCCCCGGTGGCTGCAGGCACCggggc                                                                                                                                        <  4:162838/149‑3 (MQ=255)
               acGTTGAGCCGAGGCCGCTCAGGCCCCTTGTCGGTCAAGTAGAGGCTCAGGCCAATACTGCAGCCCATGCCGATATCAACAAAACCCCCGCCAACACTTGGTTGAACACCCGCAACCGCCCCGGTGCCTGCAGCCACCGGGCGCTcccc                                                                                                                                  >  2:383119/1‑149 (MQ=255)
                 gTTGAGCCGAGGCCGCTCAGGCCCCTTGTCGGTCAAGTAGAGGCTCAGGCCAATACTGCAGCCCATGCCGATATCAACAAAACCCCCGCCAACACTTGGTTGAACACCCGCAACCGCCCCGGTGCCTGCAGCCACCGGGCGCTCCCCAc                                                                                                                                >  4:59289/1‑149 (MQ=255)
                       ccGAGGCCGCTCAGGCCCCTTGTCGGTCAAGTAGAGGCTCAGGCCAATACTGCAGCCCATGCCGATATCAACAACACCCCCGCCCACACATGGTTGGCCACCCCCCAACGCCCCCGTGCCGGCACCCCCCGGGCCCTCCCCCCCCCCaa                                                                                                                          >  1:39916/1‑149 (MQ=255)
                              cGCTCAGGCCCCTTGTCGGTCAAGTAGAGGCTCAGGCCAATACTGCAGCCCATGCCGATATCAACAAAACCCCCGCCAACACTTGGTTGAACACCCGCAACCGCCCCGGTGCCTGCCGCCCCCGGGCCCGCCCCCCCCCCCACcccccc                                                                                                                   >  8:17837/1‑149 (MQ=255)
                              cGCTCAGGCCCCTTGTCGGTCAAGTAGAGGCTCAGGCCAATACTGCAGCCCATGCCGATATCAACAAAACCCCCGCCAACACTTGGTTGAACACCCGCAACCGCCCCGGTGCCTGCAGCCACCGGGCGCTCCCCACCCCCAACACCGcc                                                                                                                   <  5:388908/149‑1 (MQ=255)
                                       cccTTGTCGGTCAAGTAGAGGCTCAGGCCAATACTGCAGCCCATGCCGATATCAACAAAACCCCCGCCAACACTTGGTTGAACACCCGCAACCGCCCCGGTGCCTGCAGCCACCGGGCGCTCCCCACCCCCAACACCCCCCCggccgc                                                                                                           >  2:354266/1‑148 (MQ=255)
                                        ccTTGTCGGTCAAGTAGAGGCTCAGGCCAATACTGCAGCCCATGCCGATATCAACAAAACCCCCGCCAACACTTGGTTGAACACCCGCAACCGCCCCGGTGCCTGCCGCCACCGGGCGGTCCCCCCCCCCACCCCCGCCCCggcggc                                                                                                           >  4:51913/1‑147 (MQ=255)
                                        ccTTGTCGGTCAAGTAGAGGCTCAGGCCAATACTGCAGCCCATGCCGATATCAACAAAACCCCCGCCAACACTTGGTTGAACACCCGCAACCGCCCCGGTGCCTGCAGCCACCGGGCGCTCCCCACCCCCAACACCGCCCAGGCGGCCa                                                                                                         <  1:139633/149‑1 (MQ=255)
                                            gTCGGTCAAGTAGAGGCTCAGGCCAATACTGCAGCCCATGCCGATATCAACAAAACCCCCGCCAAAACTTGGTTGAACACCCGCAACCGCCCCCGTGCCAGCAGCCACCGGGCGCAGCCCACCCCCAACACCCCCCCCGCGGCCAttca                                                                                                     >  3:352302/1‑146 (MQ=255)
                                              cGGTCAAGTAGAGGCTCAGGCCAATACTGCAGCCCATGCCGATATCAACAAAACCCCCGCCAACACTTGGTTGAACACCCGCAACCGCCCCGGTGCCTGCAGCCACCGGGCGCTCCCCACCCCCAACACCGCCCAGGCGGCCATGCTCg                                                                                                   <  8:363279/149‑1 (MQ=255)
                                                   aaGTAGAGGCTCAGGCCAATACTGCAGCCCATGCCGATATCAACAAAACCCCCGCCAACACTTGGTTGAACACCCGCAACCGCCCCGGTGCCTGCAGCCACCGGGCGCTCCCCACCCCCAACACCGCCCAGGCGGCCATGCTCGGCaga                                                                                              <  1:383119/149‑1 (MQ=255)
                                                            cTCAGGCCAATACTGCAGCCCATGCCGATATCAACAAAACCCCCGCCAACACTTGGTTGAACACCCGCAACCGCCCCGGTGCCTGCAGCCACCGGGCGCTCCCCACCCCCAACACCGCCCAGGCGGCCATGCTCGGCAGAGCGATCAGc                                                                                     <  1:354266/149‑1 (MQ=255)
                                                               aGGCCAATACTGCAGCCCATGCCGATATCAACAAAACCCCCGCCAACACTTGGTTGAACACCCGCAACCGCCCCGGTGCCTGCAGCCACCGGGCGCTCCCCACCCCCAACACCGCCCAGGCGGCCATGCTCGGCAGAGCGATCAGCAAg                                                                                  <  2:70297/149‑1 (MQ=255)
                                                                   cAATACTGCAGCCCATGCCGATATCAACAAAACCCCCGCCAACACTTGGTTGAACACCCGCAACCGCCCCGGTCCCTGCAGCCACCGGGCGCTCCCCACACCCAACACCCCCCAGGCCGCCATCAACGGCAGAGCGATCAGCAAGAaca                                                                              >  8:325593/1‑149 (MQ=255)
                                                                    aaTACTGCAGCCCATGCCGATATCAACAAAACCCCCGCCAACACTTGGTTGAACACCCGCAACCGCCCCGGTGCCGGCAGCCACCGGGCGCTCCCCACCCCCCCCACCGCCCAGGCGGCCATGCTCGGCAGAGCGATGAGCAAGAccac                                                                             <  2:39916/149‑1 (MQ=255)
                                                                               ccATGCCGATATCAACAAAACCCCCGCCAACACTTGGTTGAACACCCGCAACCGCCCCGGTGCCTGCAGCCACCGGGCGCTCCCCACCCCCAACACCGCCCAGGCGGCCATGCTCGGCAGAGCGATCAGCAAGAACACCAGCGCCAGCt                                                                  <  3:253083/149‑1 (MQ=255)
                                                                                 aTGCCGATATCAACAAAACCCCCGCCAACACTTGGTTGAACACCCGCAACCGCCCCGGTGCCTGCACCCACCGGGCGCACCCCACCCCCAACACCGCCCCGGCGGCCATGCTCGGCAGAGCGATCAGCAAGAACACCAGCGCCAGCTGc                                                                >  1:79529/1‑149 (MQ=255)
                                                                                  tGCCGATATCAACAAAACCCCCGCCAACACTTGGTTGAACACCCGCAACCGCCCCGGTGCCTGCAGCCACCGGGCGCTCCCCACCCCCAACACCGCCCAGGCGGCCATGCTCGGCAGAGCGATCAGCAAGAACACCAGCGCCAGCTGcc                                                               >  2:77872/1‑149 (MQ=255)
                                                                                         atCAACAAAACCCCCGCCAACACTTGGTTGAACACCCGCAACCGCCCCGGTGCCTGCAGCCACCGGGCGCTCCCCACCCCCAACACCGCCCAGGCGGCCATGCTCGGCAGAGCGATCAGCAAGAACACCAGCGCCAGCTGCCACACTgg                                                        >  2:344994/1‑149 (MQ=255)
                                                                                          tCAACAAAACCCCCGCCAACACTTGGTTGAACACCCGCAACCGCCCCGGTGCCTGCAGCCACCGGGCGCTCCCCACCCCCAACACCGCCCAGGCGGCCATGCTCGGCAGAGCGATCAGCAAGAACACCAGCGCCAGCTGCCACACTgg                                                        <  8:274477/148‑1 (MQ=255)
                                                                                          tCAACAAAACCCCCGCCAACACTTGGTTGAACACCCGCAACCGCCCCGGTGCCTGCAGCCACCGGGCGCTCCCCACCCCCAACACCGCCCAGGCGGCCATGCTCGGCAGAGCGATCAGCAAGAACACCAGCGCCAGCTGCCACACTGGc                                                       <  6:224208/149‑1 (MQ=255)
                                                                                                 aaCCCCCGCTAACACTTGGTTGAACACCCGCAACCGCCCCGGTGCCTGCAGCCACCGGGCGCTCCCCACCCCCAACACCGCCCAGGCGGCCAt                                                                                                        <  1:225084/93‑1 (MQ=255)
                                                                                                 aaCCCCCGCTAACACTTGGTTGAACACCCGCAACCGCCCCGGTGCCTGCAGCCACCGGGCGCTCCCCACCCCCAACACCGCCCAGGCGGCCAt                                                                                                        >  2:225084/1‑93 (MQ=255)
                                                                                                       cGCCAACACTTGGTTGAACACCCGCAACCGCCCCGGTGCCTGCAGCCACCGGGCGCTCCCCACCCCCAACACCGCCCAGGCGGCCATGCTCGGCAGAGCGATCAGCAAGAACACCAGCGCCAGCTGCCACACTGGCAACGACGGCGCgg                                          <  1:344994/149‑1 (MQ=255)
                                                                                                           aaCACTTGGTTGAACACCCGCAACCGCCCCGGTGCCTGCAGCCAACGGGCGCTCCCCACCCCCAACACCGCCCAGGCGGCCATGCTCGGCAGAGCGATCAGCAAGAACACCAGCGCCAGCTGCCACACTGGCAACGACGGCGCGGCaaa                                      <  8:243674/149‑1 (MQ=255)
                                                                                                                ttGGTTGAACACCCGCAACCGCCCCGGTGCCGGCAGCCACCGGGCGCTCCCCACCCCCAACACCGCCCAGGCGGCCATGCTCGGCAGAGCGATCAGCAAGAACACCAGCGCCAGCTGCCACACTGGCAACGACGGCGCGGCAAATAcac                                 >  3:35643/1‑147 (MQ=255)
                                                                                                                  ggTTGAACACCCGCAACCGCCCCGGTGCCTGCAGCCACCGGGCGCTCCCCACCCCCAACACCGCCCAGGCGGCCATGCTCGGCAGAGCGATCAGCAAGAACACCAGCGCCAGCTGCCACACTGGCAACGACGGCGCGGCAAATAcccc                                >  4:90307/1‑148 (MQ=255)
                                                                                                                      gAACACCCGCAACCGCCCCGGTGCCTGCAGCCACCGGGCGCTCCCCACCCCCAACACCGCCCAGGCGGCCATGCTCGGCAGAGCGATCAGCAAGAACACCAGCGCCAGCTGCCACACTGGCAACGACGGCGCGGCAAATACCCCGacca                           >  4:139027/1‑149 (MQ=255)
                                                                                                                         cacCCGCAACCGCCCCGGTGCCTGCAGCCACCGGGCGCTCCCCACCCCCAACACCGCCCAGGCGGCCATGCTAGGCAGAGCGATCAGCAAGAACACCAGCGCCAGCTGCCACACTGGCAACGACGGCGCGGCAAATACCCCGACCACCg                        >  6:334211/1‑149 (MQ=255)
                                                                                                                            ccGCAACCGCCCCGGTGCCTGCAGCCACCGGGCGCTCCCCACCCCCAACACCGCCCAGGCGGCCATGCTCGGCAGAGCGATCAGCAAGAACACCAGCGCCAGCTGCCACACTGGCAACGACGGCGCGGCAAATACCCCGACCACCGCCa                     <  7:17837/149‑1 (MQ=255)
                                                                                                                            ccGCAACCGCCCCGGTGCCTGCAGCCACCGGGCGCTCCCCACCCCCAACACCGCCCAGGCGGCCATGCTCGGCAGAGCGATCAGCAAGAACACCAGCGCCAGCTGCCACACTGGCAACGACGGCGCGGCAAATACCCCGACCACCGCCa                     >  8:173424/1‑149 (MQ=255)
                                                                                                                            ccGCAACCGCCCCGGTGCCTGCAGCCACCGAGCGCTCCACACCCCCAACAACGCCCAGGCGGCCATGCTAGGCAGAGCGATCAGCAAGCACACCCGCGCCAGCTGCCACACTGGCAACGACGGCGCGGCAAATACCCCGACCACCGccc                     >  1:35545/1‑148 (MQ=255)
                                                                                                                              gCAACCGCCCCGGTGCCTGCAGCCACCGGGCGCTCCCCACCCCCAACACCGCCCAGGCGGCCATGCTCGGCAGCGCGATCAGCAAGAACACCAGCGCCCGCTGCCCCACTGGCAACGACGGCGCGGCAAATACCCCGACCACCGCCACg                   >  1:334789/1‑149 (MQ=255)
                                                                                                                               cAACCGCCCCGGTGCCTGCAGCCACCGGGCGCTCCCCACCCCCAACACCGCCCAGGCGGCCATGCTCGGCAGAGCGATCAGCAAGAACACCAGCGCCAGCTGCCACACTGGCAACGACGGCGCGGCAAATACCCCGACCACCGCCACgg                  <  3:59289/149‑1 (MQ=255)
                                                                                                                                    gCCCCGGTGCCTGCAGCCACCGGGCGCTCCCCACCCCCAACACCGCCCAGGCGGCCATGCTCGGCAGAGCGATCAGCAAGAACACCAGCGCCAGCTGCCACACTGGCAACGACGGCGCGGCAAATACCCCGCCCACCGCCACGGccct              >  8:129319/1‑146 (MQ=255)
                                                                                                                                       ccGGTGGCTGCAGGCACCGGGCGCTCCCCACCCCCAACACCGCCCAGGCGGCCATGCTCGGCAGAGCGATCAGCAAGAACACCAGCGCCAGCTGCCACACTGGCAACGACGGCGCGGCAAATACCCCGACCACCGCCACGGCCATCAAc          >  7:306650/1‑149 (MQ=255)
                                                                                                                                       ccGGTGGCTGCAGGCACCGGGCGCTCCCCACCCCCAACACCGCCCAGGCGGCCATGCTCGGCAGAGCGATCAGCAAGAACACCAGCGCCAGCTGCCACACTGGCAACGACGGCGCGGCAAATACCCCGACCACCGCCACGGCCATCAAc          >  6:123468/1‑149 (MQ=255)
                                                                                                                                           tggcTGCAGGCACCGGGCGCTCCCCACCCCCAACACCGCCCAGGCGGCCATGCTCGGCAGAGCGATCAGCAAGAACACCAGCGCCAGCTGCCACACTGGCAACGACGGCGCGGCAAATACCCCGACCACCGCCACGGCCATCAACCAGa      <  8:306650/146‑1 (MQ=255)
                                                                                                                                              cTGCAGCCACCGGGCGCTCCCCACCCCCAACACCGCCCAGGCGGCCATGCTCGGCAGAGCGATCAGCAAGAACACCAGCGCCAGCTGCCACACTGGCAACGACGGCGCGGCAAATACCCCGACCACCGCCACGGCCATCAACCAGACCt   >  6:144362/1‑149 (MQ=255)
                                                                                                                                               tGCAGCCCCCGCGCGCTCCCCACCCCCAACACCGCCCAGGCGGCCATGCTCGCCAGAGCGAACAGCAACAACACCAGCGCCAGATGCCACACTCGCAACGACCGCGCGGAAAATACCCCGACCACACCCACGGCCATCAACCAGACAtt  >  2:305431/1‑149 (MQ=255)
                                                                                                                                               tGCAGCCACCGGGCGCTCCCCACCCCCAACACCGCCCAGGCGGCCATGCTCGGCAGAGCGATCAGCAAGAACACCAGCGCCAGCTGCCACACTGGCAACGACGGCGCGGCAAATACCCCGACCACCGCCACGGCCATCAACCAGACCtt  >  3:313534/1‑149 (MQ=255)
                                                                                                                                                    |                                                                                                                                               
TCGCGAGGCGTTGACACGTTGAGCCGAGGCCGCTCAGGCCCCTTGTCGGTCAAGTAGAGGCTCAGGCCAATACTGCAGCCCATGCCGATATCAACAAAACCCCCGCCAACACTTGGTTGAACACCCGCAACCGCCCCGGTGCCTGCAGCCACCGGGCGCTCCCCACCCCCAACACCGCCCAGGCGGCCATGCTCGGCAGAGCGATCAGCAAGAACACCAGCGCCAGCTGCCACACTGGCAACGACGGCGCGGCAAATACCCCGACCACCGCCACGGCCATCAACCAGACCTT  >  NC_002947/2776859‑2777150

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