Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 5,538,335 G→T 12.9% W122L (TGG→TTG)  PP_4871 → lysine decarboxylase family protein

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0029475,538,3350GT12.9% 77.8 / 3.6 31W122L (TGG→TTG) PP_4871lysine decarboxylase family protein
Reads supporting (aligned to +/- strand):  ref base G (13/14);  new base T (2/2);  total (15/16)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 8.58e-01

CGCTGAAGTCGGCCACAAGGGCCTGACCCGCCTGGAAGTGGTGGATGGCATGCATGCGCGCAAGGCCCGCATGGCCGAGCTGAGCGATGCATTCATTGCCTTGCCGGGTGGGCTGGGTACGCTGGAGGAACTGTTCGAAGTGTGGACCTGGGGGCAGCTGGGCTATCACGCCAAGCCGCTCGGGCTACTGGATGTGAACGGTTTCTACGACAAGCTGGGCGGGTTCCTCGACCATATCGTCGAAGAAGGTTTCGTGCGCCCGCAACACCGGGCGATGTTGCTGCTGGGGCAGCA  >  NC_002947/5538186‑5538479
                                                                                                                                                     |                                                                                                                                                 
cGCTGAAGTCGGCCACAAGGGCCTGACCCGCCTGGAAGTGGTGGATGGCATGCATGCGCGCAAGGCCCGCATGGCCGAGCTGAGCGATGCATTCATTGCCTTGCCGGGTGGGCTGGGTACGCTGGAGGAACTGTTCGAAGTGTGGACCTg                                                                                                                                                   <  2:26234/150‑1 (MQ=255)
            ccgcaAGGGCCTGACCCGCCTGGAAGTGGTGGATGGCATGCATGCGCGCAAGGCCCGCATGGCCGAGCTGAGCGATGCCTTCATTGCCTTGCCGGGTGGGCTGGGTACGCTGGAGGAACTGTTCGAAGTGTGGACCTGGGGGCAGCTggg                                                                                                                                       <  2:591176/147‑1 (MQ=255)
            ccACAAGGGCCTGACCCGCCTGGAAGTGGGGGATGGCAGGGATGCGCGCAAGGCCCGCATGGCCGAGCTGAGCGATGCATTCATTGCCTTGCCGGGTGGACTCGGTACGCTGGACGAAGGGGGCGAAGTGGTGACCTGGGGGCAGCTggg                                                                                                                                       >  2:375614/1‑150 (MQ=255)
                   ggCCTGACCCGCCTGGAAGTGGTGGATGGCATGCATGCGCGCAAGGCCCGCATGGCCGAGCTGAGCGATGCATTCATTGCCTTGCCGGGTGGGCTGGGTACGCTGGAGGAACTGTTCGAAGTGTGGACCTGGGGGCAGCTGGGCTATCAc                                                                                                                                <  1:585336/150‑1 (MQ=255)
                          cccGCCTGGAAGTGGTGGATGGCATGCATGCGCGCAAGGCCCGCATGGCCGAGCTGAGCGATGCATTCATTGCCTTGCCGGGTGGGCTGGGTACGCTGGAGGAACTGTTCGAAGTGTGGACCTGGGGGCAGCTGGGCTATCACGCCAAGc                                                                                                                         <  1:531683/150‑1 (MQ=255)
                           ccGCCTGGAAGTGGTGGATGGCATGCATGCGCGCACGGCGCGCATGATAGCGCTGAGCGATGCATTCATTTTCTTGCCGGGTGGGATGGGTACGCTGGGGGAACGGTTCGAAGTGTGGACCTGGGGGCAGCTGGGCTATCACGCCAAGcc                                                                                                                        >  1:474487/1‑150 (MQ=255)
                                tGGAAGTGGTGGATGGCATGCATGCGCGCAAGGCCCGCATGGCCGAGCTGAGCGATGCATTCATTGCCTTGCCGGGTGGGCTGGGTACGCTGGAGGAACTGTTCGAAGTGTGGACCTGGGGGCAGCTGGGCTATCACGCCAAGCCGCTCg                                                                                                                   >  2:524340/1‑150 (MQ=255)
                                        gtggATGGCATGCATGCGCGCAAGGCCCGCATGGCCGAGCTGAGCGATGCATTCATTGCCTTGCCGGGTGGGCTGGGTACGCTGGAGGAACTGTTCGAAGTGTGGACCTTGGGGCAGCTGGGCTATCACGCCAAGCCGCTCGGGCTACTg                                                                                                           >  2:252262/1‑150 (MQ=255)
                                        gtggATGGCATGCATGCGCGCAAGGCCCGCATGGCCGAGCTGAGCGATGCATTCATTGCCTTGCCGGGTGGGCTGGGTACGCTGGAGGAACTGTTCGAAGTGTGGACCTTGGGGCAGCTGGGCTATCACGCCAAGCCGCTCGGGCTACTg                                                                                                           >  2:252299/1‑150 (MQ=255)
                                           gATGGCATGCATGCGCGCAAGGCCCGCATGGCCGAGCTGAGCGATGCATTCATTGCCTTGCCGGGTGGGCTGGGTACGCTGGAGGAACTGTTCGAAGTGTGGACCTGGGGGCAGCTGGGCTATCACGCCAAGCCGCTCGGGCTACTGGAt                                                                                                        <  2:474487/150‑1 (MQ=255)
                                                         gcgcAAGGCCCGCATGGCCGAGCTGAGCGATGCATTCATTGCCTTGCCGGGTGGGCTGGGTACGCTGGAGGAACTGTTCGAAGTGTGGACCTGGGGGCAGCTGGGCTATCACGCCAAGCCGCTCGGGCTACTGGATGTGAACGGTTTCTa                                                                                         <  1:555915/150‑1 (MQ=255)
                                                                          ccGAGCTGAGCGATGCATTCATTGCCTTGCCGGGTGGGCTGGGTACGCTGGAGGAACTGTTCGAAGTGTGGACCTGGGGGCAGCTGGGCTATCACGCCAAGCCGCTCGGGCTACTGGATGTGAACGGTTTCTACGACAAGCTGGGCGGGt                                                                        >  2:35160/1‑150 (MQ=255)
                                                                           cGAGCTGAGCGATGCATTCATTGCCTTGCCGGGTGGGCTGGGTACGCTGGAGGAACTGTTCGAAGTGTGGACCTGGGGGCAGCTGGGCTATCACGCCAAGCCGCTCGGGCTACTGGATGTGAACGGTTTCTACGACAAGCTGGGCGGGtt                                                                       <  1:427433/150‑1 (MQ=255)
                                                                             aGCTGAGCGATGCATACATTGGCTTGCCGGGTGGGCTGGGTACGCTGGAGGAACTGTTCGAAGTGTGGACCTGGGGGCAGCTGGGCTATCACGCCAAGCCGCTCGGGCTACTGGATGTGAACGGTTTCTACGACAAGCTGGGCGGGTTcc                                                                     <  1:217090/150‑1 (MQ=255)
                                                                              gCTGAGCGATGCATTCATTGCCTTGCCGGGTGGGCTGGGTACGCTGGAGGAACTGTTCGAAGTGTGGACCTGGGGGCAGCTGGGCTATCACGCCAAGCCGCTCGGGCTACTGGAtgt                                                                                                      >  1:10399/1‑117 (MQ=255)
                                                                              gCTGAGCGATGCATTCATTGCCTTGCCGGGTGGGCTGGGTACGCTGGAGGAACTGTTCGAAGTGTGGACCTGGGGGCAGCTGGGCTATCACGCCAAGCCGCTCGGGCTACTGGAtgt                                                                                                      <  2:10399/117‑1 (MQ=255)
                                                                                tGAGCGATGCATTCATTGCCTTGCCGGGTGGGCTGGGTACGCTGGAGGAACTGTTCGAAGTGTGGACCTGGGGGCAGCTGGGCTATCACGCCAAGCCGCTCGGGCTACTGGa                                                                                                         >  1:584555/1‑112 (MQ=255)
                                                                                tGAGCGATGCATTCATTGCCTTGCCGGGTGGGCTGGGTACGCTGGAGGAACTGTTCGAAGTGTGGACCTGGGGGCAGCTGGGCTATCACGCCAAGCCGCTCGGGCTACTGGa                                                                                                         <  2:584555/112‑1 (MQ=255)
                                                                                 gAGCGATGCATTCATTGCCTTGCCGGGTGGGCTGGGTACGCTGGGGGAACTGTTCGAAGTGTGGACCTGGGGGCAGCTGGGCTATCACGCCAAGCCGCTCGGGCTACTGGATGTGAACGGTTTCTACGACAAGCTGGGCGGGTTCCt                                                                    >  2:140203/1‑147 (MQ=255)
                                                                                 gAGCGATGCATTCATTGCCTTGCCGGGTGGGCTGGGTACGCTGGAGGAACTGTTCGAAGTGTGGACCTGGGGGCAGCTGGGCTATCACGCCAAGCCGCTCGGGCTACTGGATGTGAACGGTTTCTACGACAAGCTGGGCGGGTTCCt                                                                    <  1:140203/147‑1 (MQ=255)
                                                                                          attcattGCCTTGCCGGGTGGGCTGGGTACGCTGGAGGAACTGTTCGAAGTGTGGACCTGGGGGCAGCTGCGCTTTCATGCCAAGCCGCTAGGGCTAATGAGGTGGAAGGGTTTCGACGACAAGCTGGGCGGGTTCCTCGACCATAtcgt                                                        >  1:139343/1‑150 (MQ=255)
                                                                                           ttcattGCCTTGCCGGGTGGGCTGGGTACGCTGGAGGAACTGTTCGAAGTGTGGACCTTGGGGCAGCTGGGCTATCACGCCAAGCCGCTCGGGCTACTGGATGTGAACGGTTTCTACGACAAGCTGGGCGGGTTCCTCGACCATAtcgtc                                                       <  1:252299/150‑1 (MQ=255)
                                                                                           ttcattGCCTTGCCGGGTGGGCTGGGTACGCTGGAGGAACTGTTCGAAGTGTGGACCTTGGGGCAGCTGGGCTATCACGCCAAGCCGCTCGGGCTACTGGATGTGAACGGTTTCTACGACAAGCTGGGCGGGTTCCTCGACCATAtcgtc                                                       <  1:252262/150‑1 (MQ=255)
                                                                                                   cttgccGGGTGGGCTGGGTACGCTGGAGGAACTGTTCGAAGTGTGGACCTGGGGGCAGCTGGGCTATCACGCCAAGCCGCTCGGGCTACTGGATGTGAACGGTTTCTAc                                                                                        >  2:145079/1‑109 (MQ=255)
                                                                                                   cttgccGGGTGGGCTGGGTACGCTGGAGGAACTGTTCGAAGTGTGGACCTGGGGGCAGCTGGGCTATCACGCCAAGCCGCTCGGGCTACTGGATGTGAACGGTTTCTAc                                                                                        <  1:145079/109‑1 (MQ=255)
                                                                                                    ttgccGGGTGGGCTGGTTACGCTGGAGGAACTGTTCGAAGTGTGGACCTGGGGGCAGCTGGGCTATCACGCCAAGCCGCTCGGGCTACTGGATGTGAACGGTTTCTACGACAAGCTGGGCGGGTTCCTCGACCATATCGTCGAAGAAGGt                                              <  1:524340/150‑1 (MQ=255)
                                                                                                              ggCTGGGTACGCTGGAGGAACTGTTCGAAGTGTGGACCTGGGGGCAGCTGGGATATCCCGAACCGCCGCTCGCGATCCGGGATGTGAACGGTTCCTACGACAAGCTGTGCTGGTTCCTCGACCATATCGTCGAAGAAGGTTTCGTGCGcc                                    <  2:378187/150‑1 (MQ=255)
                                                                                                                            gaggaACTGTTCGAAGTGTGGACCTGGGGGCAGCTGGGCTATCACGCCAAGCCGCTCGGGCTACTGGATGTGAACGGTTTCTACGACAAGCTGGGCGGGTTCCTCGACCATATCGTCGAAGAAGGTTTCGTGCGCCCGCAACACCGGGCg                      >  1:497002/1‑150 (MQ=255)
                                                                                                                                   tGTTCGAAGTGTGGACCTGGGGGCAGCTGGGCTATCACGCCAAGCCGCTCGGGCTACTGGATGTGAACGGTTTCTACGACAAGCTGGGCGGGTTCCTCGACCATATCGTCGAAGACAGTTTCGTGCGCACGCCACACCGGGCGATGTtgt               >  2:139147/1‑149 (MQ=255)
                                                                                                                                            tgtgGACCTGGGGGCAGCTGGGCTATCACGCCAAGCCGCTCGGGCTACTGGATGTGAACGGTTTCTACGACAAGCTGGGCGGGTTCCTCGACCATATCGTCGAAGAAGGTTTCGTGCGCCCGCAACACCGGGCGATGTTGCTGCTGGGgc      >  2:391053/1‑150 (MQ=255)
                                                                                                                                              tgGACCTGGGGGCAGCTGGGCTATCACGCCAAGCCGCTCGGGCTACTGGATGTGAACGGTTTCTACGACAAGCTGGGCGGGTTCCTCGACCATATCGTCGAAGAAGGTTTCGTGCGCCCGCAACACCGGGCGATGTTGCTGCTGGGgcag    <  1:35160/150‑1 (MQ=255)
                                                                                                                                                gACCTGGGGGCAGCTGGGCTATCACGCCAAGCCGCTCGGGCTACTGGATGTGAACGGTTTCTACGACAAGCTGGGCGGGTTCCTCGACCATATCGTCGAAGAAGGTTTCGTGCGCCCGCAACACCGGGCGATGTTGCTGCTGGGgcagca  >  1:206857/1‑150 (MQ=255)
                                                                                                                                                     |                                                                                                                                                 
CGCTGAAGTCGGCCACAAGGGCCTGACCCGCCTGGAAGTGGTGGATGGCATGCATGCGCGCAAGGCCCGCATGGCCGAGCTGAGCGATGCATTCATTGCCTTGCCGGGTGGGCTGGGTACGCTGGAGGAACTGTTCGAAGTGTGGACCTGGGGGCAGCTGGGCTATCACGCCAAGCCGCTCGGGCTACTGGATGTGAACGGTTTCTACGACAAGCTGGGCGGGTTCCTCGACCATATCGTCGAAGAAGGTTTCGTGCGCCCGCAACACCGGGCGATGTTGCTGCTGGGGCAGCA  >  NC_002947/5538186‑5538479

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