Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 4,505,657 T→C 9.1% L144P (CTG→CCG)  PP_3997 → glycosyl transferase

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0029474,505,6570TC 9.1% 114.7 / 3.2 44L144P (CTG→CCG) PP_3997glycosyl transferase
Reads supporting (aligned to +/- strand):  ref base T (21/19);  new base C (2/2);  total (23/21)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 7.48e-01

GCCCACGTACGCCGGCAAGAAGCGTCACCTGCCCTGGTACCTGCTGGCCGCCAAGTGCCTGGCCAACAATGGCGTGCGCATTCTGATGCACGGCGGCGGCGCACACA‑CTGCCGGGCGTATGTATACAGAGCAACTCCTCGACCTGCTGCAGATTCCGCTGTGCCGCGACTGGGCCGCCGTCGGCAACGCGCTGGACCAGCACCAGCTGGCATTCTTCCCGCTGCACGACTGGGCGCCGCAGCTGCAACGCATGATCGACCTGCGCAATACCCTGGGCCTGCGCTCGCCCAT  >  NC_002947/4505511‑4505801
                                                                                                                                                   |                                                                                                                                                
gCCCACGTACGCCGGCAAGAAGCGTCACCTGCCCTGGTACCTGCTGGCCGCCAAGTGCCTGGCCAACAATGGCGTGCGCATTCTGATGCACGGCGGCGGCGCACACA‑CTGTCGGGCGTATGTATACAGAGCAACTCCTCGACctgct                                                                                                                                                  <  7:28561/147‑1 (MQ=255)
    aCGTACGGCGGCAAGTAGCGTCACCTGCCCTGGTACCTGCTGGCCGCCAAGTGCCTGGCCAACAATGGCGTGCGGATTCTGATGCACGGCGGCGGCGCACACACCTGCCGGGCGTATGTATACAGAGCAACTCCTCGACCTGCTGCAGa                                                                                                                                             <  1:73542/149‑1 (MQ=255)
             ggCAAGAAGCGTCACCTGCCCTGGTACCTGCTGGCCGCCAAGTGCCTGGCCAACAATGGCGTGCGCATTCTGATGCACGGCGGCGGCGCACACA‑CTGCCGGGCGTATGTATACAGAGCAACTCCTCGACCTGCTGCAGATTCCGCtgtg                                                                                                                                   <  6:260181/149‑1 (MQ=255)
                        tCACCTGCCCTGGTACCTGCTGGCCGCCAAGTGCCTGGCCAACAATGGCGTGCGCATTCTGATGCACGGCGGCGGCGCACACA‑CTGCCGGGCGTATGTATACAGAGCAACTCCTCGACCTGCTGCAGATTCCGCTGTGCCGCGACTggg                                                                                                                        >  2:218731/1‑149 (MQ=255)
                          aCCTGCCCTGGTACCTGCTGGCCGCCAAGTGCCTGGCCAACAATGGCGTGCGCATTCTGATGCACGGCGGCGGCGCACACA‑CTGCCGGGCGTATGTATACAGAGCAACTCCTCGACCTGCTGCAGATTCCGc                                                                                                                                       <  8:134813/132‑1 (MQ=255)
                          aCCTGCCCTGGTACCTGCTGGCCGCCAAGTGCCTGGCCAACAATGGCGTGCGCATTCTGATGCACGGCGGCGGCGCACACA‑CTGCCGGGCGTATGTATACAGAGCAACTCCTCGACCTGCTGCAGATTCCGc                                                                                                                                       >  7:134813/1‑132 (MQ=255)
                            cTGCCCTGGTACCTGCTGGCCGCCAAGTGCCTGGCCAACAATGGCGTGCGCATTCTGATGCACGGCGGCGGCGCACACA‑CTGCCGGGCGTATGTATACAGAGCAACTCCTCGACCTGCTGCAGATTCCGCTGTGCCGCGACTGGgcc                                                                                                                      <  3:262928/147‑1 (MQ=255)
                            cTGCCCTGGTACCTGCTGGCCGCCAAGTGCCTGGCCAACAATGGCGTGCGCATTCTGATGCACGGCGGCGGCGCACACA‑CTGCCGGGCGTATGTATACAGAGCAACTCCTCGACCTCCTGCAGATTCCGCTGTGCCGCGACTcggac                                                                                                                      >  4:262928/1‑145 (MQ=255)
                                       ccTGCTGGCCGCCAAGTGCCTGGCCAACAATGGCGTGCGCATTCTGATGCACGGCGGCGGCGCACACA‑CTGCCGGGCGTATGTATACAGAGCAACTCCTCGTCCTGCTGCAGATTCCGCTGTGCCGCGACTGGGCCGCCGTCGGCAAcg                                                                                                         >  5:91105/1‑149 (MQ=255)
                                       ccTGCTGGCCGCCAAGTGCCTGGCCAACAATGGCGTGCGCATTCTGATGCACGGCGGCGGCGCACACA‑CTGCCGGGCGTATGTATACAGAGCAACTCCTCGACCTGCTGCAGATTCCGCTGTGCCGCGACTGGGCCGCCGTCGGCAAcg                                                                                                         >  5:86032/1‑149 (MQ=255)
                                        ctgctgGCCGCCAAGTGCCTGGCCAACAATGGCGTGCGCATTCTGATGCACGGCGGCGGCGCACACA‑CTGCCGGGCGTATGTATACAGAGCAACTCCTCGACCTGCTGCAGATTCCGCTGTGCCGCGACTGGGCCGCCGTCGGCa                                                                                                            >  5:229743/1‑145 (MQ=255)
                                        ctgctgGCCGCCAAGTGCCTGGCCAACAATGGCGTGCGCATTCTGATGCACGGCGGCGGCGCACACA‑CTGCCGGGCGTATGTATACAGAGCAACTCCTCGACCTGCTGCAGATTCCGCTGTGCCGCGACTGGGCCGCCGTCGGCa                                                                                                            <  6:229743/145‑1 (MQ=255)
                                           ctgGCCGCCAAGTGCCTGGCCAACAATGGCGTGCGCATTCTGATGCACGGCGGCGGCGCACACA‑CTGCCGGGCGTATGTATACAGAGCAACTCCTCGACCTGCCGCAGATTCCGCTGTGCCGCGa                                                                                                                             >  1:322114/1‑125 (MQ=255)
                                           ctgGCCGCCAAGTGCCTGGCCAACAATGGCGTGCGCATTCTGATGCACGGCGGCGGCGCACACA‑CTGCCGGGCGTATGTATACAGAGCAACTCCTCGACCTGCCGCAGATTCCGCTGTGCCGCGa                                                                                                                             <  2:322114/125‑1 (MQ=255)
                                                 gccAAGTGCCTGGCCAACAATGGCGTGCGCATTCTGATGCACGGCGGCGGCGCACACA‑CTGCCGGGCGTATGTATACAGAGCAACTCCTCGACCTGCTGCAGATTCCGCTGTGCCGCGACTGGGCCGCCGTCGGCAACGCGCTGGac                                                                                                 >  1:199429/1‑147 (MQ=255)
                                                             gCCAACAATGGCGTGCGCATTCTGATGCACGGCGGCGGCGCACACA‑CTGCCGGGCGTATGTATACAGAGCAACTCCTCGACCTGCTGCAGATTCCGCTGTGCCGCGACTGGGCCGCCGTCGGCAACGCGCTGGACCAGCACCAGCTGGc                                                                                   <  6:91105/149‑1 (MQ=255)
                                                             gCCAACAATGGCGTGCGCATTCTGATGCACGGCGGCGGCGCACACA‑CTGCCGGGCGTATGTATACAGAGCAACTCCTCGACCTGCTGCAGATTCCGCTGTGCCGCGACTGGGCCGCCGTCGGCAACGCGCTGGACCAGCACCAGCTGGc                                                                                   <  3:75905/149‑1 (MQ=255)
                                                             gCCAACAATGGCGTGCGCATTCTGATGCACGGCGGCGGCGCACACA‑CTGCCGGGCGTATGTATACAGAGCAACTCCTCGACCTGCTGCAGATTCCGCTGTGCCGCGACTGGGCCGCCGTCGGCAACGCGCTGGACCAGCACCAGCTGGc                                                                                   <  6:86032/149‑1 (MQ=255)
                                                                   aaTGGCGTGCGCATTCTGATGCACGGCGGCGGCGCACACA‑CTGCCGGGCGTATGTATACAGAGCAACTCCTCGACCTGCTGCAGATTCCGCTGTGCCGCGACTGGGCCGCCGTCGGCAACGCGCTGGac                                                                                                 <  3:86834/129‑1 (MQ=255)
                                                                   aaTGGCGTGCGCATTCTGATGCACGGCGGCGGCGCACACA‑CTGCCGGGCGTATGTATACAGAGCAACTCCTCGACCTGCTGCAGATTCCGCTGTGCCGCGACTGGGCCGCCGTCGGCAACGCGCTGGac                                                                                                 >  4:86834/1‑129 (MQ=255)
                                                                     tGGCGTGCGCATTCTGATGCACGGCGGCGGCGCACACA‑CTGCCGGGCGTATGTATACAGAGCAACTCCTCGACCTGCTGCAGATTCCGCTGTGCCGCGACTGGGCCGCCGTCGGCAACGCGCTGGaccagcaccag                                                                                        <  6:17435/136‑1 (MQ=255)
                                                                     tGGCGTGCGCATTCTGATGCACGGCGGCGGCGCACACA‑CTGCCGGGCGTATGTATACAGAGCAACTCCTCGACCTGCTGCAGATTCCGCTGTGCCGCGACTGGGCCGCCGTCGGCAACGCGCTGGaccagcaccag                                                                                        >  5:17435/1‑136 (MQ=255)
                                                                                         aCGGCGGCGGCGCACACA‑CTGCCGGGCGTATGTATACAGAGCAACTCCTCGACCTGCCGCAGATTCCGCTGTGCCGCGACTGGGCCGCCGTCGGCAACGCGCTGGACCAGCACCAGCTGGCATTCTTCCCGCTGCACGACTGGGcgccg                                                       >  2:177395/1‑149 (MQ=255)
                                                                                              ggcggcgCACACA‑CTGCCGGGCGTATGTATACAGAGCAACTCCTCGACCTGCCGCAGATTCCGCTGTGCCGCGACTGGGCCGCCGTCGGCAACGCGCTGGACCAGCACCAGCTGGCATTCTTCCCGCTGCACGACTGGGCGCCGCAGCt                                                  <  1:177395/149‑1 (MQ=255)
                                                                                                        aca‑cTGCCGGGCGTATGTATACAGAGCAAATCCTCGACCTGCTGAAGATTCCTCTGTGCCGCTACTGTGCAGCCGTCTGCAACTCGCTGGACCAGCACCATATGGCATTCTTCCCGCTGCACGACTGGGCGACGCAGCTGCAAAGCAt                                         >  7:189802/1‑148 (MQ=255)
                                                                                                          a‑cTGCCGGGCGTATGTATACAGAGCAACTCCTCGACCTGCTGCAGATTCCGCTGTGCCGCGACTGGGCCGCCGTCGGCAACGCGCTGGACCAGCACCAGCTGGCATTCTTCCCGCTGCACGACTGGGCGCCGCAGCTGCAACGCATGAt                                      >  4:78280/1‑149 (MQ=255)
                                                                                                          a‑cTGCCGGGCGTATGTATACAGAGCAACTCCTCGACCTGCTGCAGATTCCGCTGTGCCGCGACTGGGCCGCCGTCGGCAACGCGCTGGACCAGCACCAGCTGGCATTCTTCCCGCTGCACGACTGGGCGCCGCAGCTGCAACGCATGAt                                      >  3:106889/1‑149 (MQ=255)
                                                                                                               ccGGGCGTATGTATACAGAGCAACTCCTCGACCTGCTGCAGATTCCGCTGTGCCGCGACTGGGCCGCCGTCGGCAACGCGCTGGACCAGCACCAGCTGGCATTCTTCCCGCTGCACGACTGGGCGCCGCAGCTGCAACGCATGATCGAc                                  <  2:199429/149‑1 (MQ=255)
                                                                                                                   gCGTATGTATACAGAGCAACTCCTCGACCTGCTGCAGATTCCGCTGTGCCGCGACTGGGCCGCCGTCGGCAACGCGCTGGACCAGCACCAGCTGGCATTCTTCCCGCTGCACGACTGGGCGCCGCAGCTGCAACGCATGATCGACCt                                <  8:313355/147‑1 (MQ=255)
                                                                                                                   gCGTATGTATACAGAGCAACTCCTCGACCTGCTGCAGATTCCGCTGTGCCGCGACTGGGCCGCCGTCGGCAACGCGCTGGACCAGCACCAGCTGGCATTCTTCCCGCTGCACGACTGGGCGCCGCAGCTGCAACGCATGATCGACCt                                >  7:313355/1‑147 (MQ=255)
                                                                                                                   gCGTATGTATACAGAGCAACTCCTCGACCTGCTGCAGATTCCGCTGTGCCGCGACTGGGCCGCCGTCGGCAACGCGCTGGACCAGCACCAGCTGGCATTCTTCCCGCTGCACGACTGGGCGCCGCAGCTGCAACGCATGATCGACCTgc                              >  7:167465/1‑149 (MQ=255)
                                                                                                                    cGTATGTAAACAGCCCACCTCCTCGACCTGCTGCCGATGCCGCTGCGCCGCGCCTGGGCCGCCGTCGGCAACGCGCTGGACCCCCACCAGCTGGCATTCTCCCCGCTGCACGACGGGCCGCCGCAGCTCCACCCCATGATCGACCGgcg                             <  8:167465/149‑1 (MQ=255)
                                                                                                                       atgtatACAGAGCAACTCCTAGACCTGCTGCAGATTCCGCTGTGCCGCGACTGGGCCGCCGTCGGCAACGCGCTGGACCAGCACCAGCTGGCATTCTTCCCGCTGCACGACTGGGCGCCGCAGCTGCAACGCATGATCGACCTGCGCa                           <  8:189802/148‑1 (MQ=255)
                                                                                                                              cAGAGCAACTCCTCGACCTGCTGCAGATTCCGCTGTGCCGCGACTGGGCCGCCGTCGGCAACGCGCTGGACCAGCACCAGCTGGCATTCTTCCCGCTGCACGACTGGGCGCCGCAGCTGCAACGCATGATCGACCTGCGCAATACCCTg                   <  4:106889/149‑1 (MQ=255)
                                                                                                                               agagCAACTCCTCGACCTGCTGCAGATTCCGCTGTGCCGCGACTGGGCCGCCGTCGGCAACGCGCTGGACGAGCACCAGCTGGCATTCTTCCCGCTGCACGACTGGGCGCCGCAGCTGCAACGCATGATCGACCTGCGCAATACCCTgg                  <  3:78280/149‑1 (MQ=255)
                                                                                                                               agagCAACTCCTCGACCTGCTGCAGATTCCGCTGTGCCGCGACTGGGCCGCCGTCGGCAACGCGCTGGACCAGCACCAGCTGGCATTCTTCCCGCTGCACGACTGGGCGCCGCAGCTGCAACGCATGATCGACCTGCGCAATACCCTgg                  >  7:288243/1‑149 (MQ=255)
                                                                                                                                     aCTCCTCGACCTGCTGCAGATTCCGCTGTGCCGCGACTGGGCCGCCGTCGGCAACGCGCTGGACCAGCACCAGCTGGCATTCTTCCCGCTGCACGACTGGGCGCCGCAGCTGCAACGCATGATCGACCTGCGCAATACCCTg                   >  3:88980/1‑142 (MQ=255)
                                                                                                                                     aCTCCTCGACCTGCTGCAGATTCCGCTGTGCCGCGACTGGGCCGCCGTCGGCAACGCGCTGGACCAGCACCAGCTGGCATTCTTCCCGCTGCACGACTGGGCGCCGCAGCTGCAACGCATGATCGACCTGCGCAATACCCTg                   <  4:88980/142‑1 (MQ=255)
                                                                                                                                         ctcGACCTGCTGCAGATTCCGCTGTGCCGCGACTGGGCCGCCGTCGGCAACGCGCTGGACCAGCACCAGCTGGCATTCTTCCCGCTGCACGACTGGGCGCCGCAGCTGCAACGCATGa                                       <  2:140306/118‑1 (MQ=255)
                                                                                                                                         ctcGACCTGCTGCAGATTCCGCTGTGCCGCGACTGGGCCGCCGTCGGCAACGCGCTGGACCAGCACCAGCTGGCATTCTTCCCGCTGCACGACTGGGCGCCGCAGCTGCAACGCATGa                                       >  1:140306/1‑118 (MQ=255)
                                                                                                                                         ctcGACCTGCTGCAGATTCCGCTGTGCCGCGACTGGGCCGCCGTCGGCAACGCGCTGGACCAGCACCAGCTGGCATTCTTCCCGCTGCACGACTGGGCGCCGCAGCTGCAACGCATGATCGACCTGCGCAATACCCTGGGCCTGCGCTc        >  5:31514/1‑149 (MQ=255)
                                                                                                                                          tcGACCTGCTGCAGATTCCGCTGTGCCGCGACTGGGCCGCCGTCGGCAACGCGCTGGACCAGCACCAGCTGGCATTCTTCCCGCTGCACGACTGGGCGCCGCAGCTGCAACGCATGATCGACCTGCGCAATACCCTGGGCCTGCGCTCg       >  7:39969/1‑149 (MQ=255)
                                                                                                                                          tcGACCTGCTGCAGATTCCGCTGTGCCGCGACTGGGCCGCCGTCGGCAACGCGCTGGACCAGCACCAGCTGGCATTCTTCCCGCTGCACGACTGGGCGCCGCAGCTGCAACGCATGATCGACCTGCGCAATACCCTGGGCCTGCGCTCg       >  7:44809/1‑149 (MQ=255)
                                                                                                                                            gACCTGCTGCAGATTCCGCTGTGCCGCGACTGGGCCGCCGTCGGCAACGCGCTGGACCAGCACCAGCTGGCATTCTTCCCGCTGCACGACTGGGCGCCGCAGCTGCAACGCATGATCGACCTGCGCAATACCCTGGGCCTGCGCTCg       <  8:288243/147‑1 (MQ=255)
                                                                                                                                              ccTGCTGCAGATTCCGCTGTGCCGCGACTGGGCCGCCGTCGGCAACGCGCTGGACCAGCACCAGCTGGCATTCTTCCCGCTGCACGACTGGGCGCCGCAGCTGCAACGCATGATCGACCTGCGCAATACCCTGGGCCTGCGCTCGCCCa   >  1:54297/1‑149 (MQ=255)
                                                                                                                                               ctgctgCAGATTCCGCTGTGCCGCGACTGGGCCGCCGTCGGCAACGCGCTGGACCAGCACCAGCTGGCATTCTTCCCGCTGCACGACTGGGCGCCGCAGCTGCAACGCATGATCGACCTGCGCAATACCCTGGGCCTGCGCTCGCCCAt  <  6:31514/149‑1 (MQ=255)
                                                                                                                                                   |                                                                                                                                                
GCCCACGTACGCCGGCAAGAAGCGTCACCTGCCCTGGTACCTGCTGGCCGCCAAGTGCCTGGCCAACAATGGCGTGCGCATTCTGATGCACGGCGGCGGCGCACACA‑CTGCCGGGCGTATGTATACAGAGCAACTCCTCGACCTGCTGCAGATTCCGCTGTGCCGCGACTGGGCCGCCGTCGGCAACGCGCTGGACCAGCACCAGCTGGCATTCTTCCCGCTGCACGACTGGGCGCCGCAGCTGCAACGCATGATCGACCTGCGCAATACCCTGGGCCTGCGCTCGCCCAT  >  NC_002947/4505511‑4505801

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