Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 4,528,861 C→A 11.4% N139K (AAC→AAA PP_4017 → hypothetical protein

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0029474,528,8610CA11.4% 97.0 / 2.7 35N139K (AAC→AAAPP_4017hypothetical protein
Reads supporting (aligned to +/- strand):  ref base C (17/14);  new base A (2/2);  total (19/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.89e-01

TGGACCCTGCTGGTGCAGGCCGTTGACCAGTTCGTTCCGGAAGTGGCCGAACTGCTGGAAAACTTCCGCTTCCTGCCCAGCTGGCGTATCGATGACGTGATGATCAGCTTCGCCGCCCCGGGTGGCAGCGTCGGCCCGCACTTCGACAACTACGATGTGTTCCTGCTGCAAGGCCACGGCCAGCGCAACTGGAAGATCGGCCAGATGTGCAGCAGTGACAGTCCGCTGCTGGAGCACGCCGACCTGCGTATCCTGGCCGAGTTCGAGCAGAGCGGCGAGTGGACCCTGGAGCCCGGCGA  >  NC_002947/4528712‑4529010
                                                                                                                                                     |                                                                                                                                                     
tGGACCCTGCTGGTGCAGGCCGTTGACCAGTTCGTTCCGGAAGTGGCCGAACTGCTGGAAAACTTCCGCTTCCTGCCCAGCTGGCGTATCGATGACGTGATGATCAGCTTCGCCGCCCCGGGTGGCAGCGTCGGCCCGCACTTCGACAAc                                                                                                                                                       <  1:659898/150‑1 (MQ=255)
               cAGGCCGTTGACCAGTTCGTTCCGGAAGTGGCCGAACTGCTGGAAAACTTCCGCTTCCTGCCCAGCTGGCGTATCGATGACGTGATGATCAGCTTCGCCGCCCCGGGTGGCAGCGTCGGCCCGCACTTCGACAACTACGATGTGTTCctg                                                                                                                                        <  2:359605/150‑1 (MQ=255)
                  gCCGTTGACCAGTTCGTTCCGGAAGTGGCCGAACTGCTGGAAAACTTCCGCTTCCTGCCCAGCTGGCGTATCGATGACGTGATGATCAGCTTCGCCGCCCCGGGTGGCAGCGTCGGCCCGCACTTCGACAAATACGATGTGTTCctgctg                                                                                                                                     >  1:756903/1‑150 (MQ=255)
                    cGTTGACCAGTTCGTTCCGGAAGTGGCCGAACTGCTGGAAAACTTCCGCTTCCTGCCCAGCTGGCGTATCGATGACGTGATGATCAGCTTCGCCGCCCCGGGTGGCAGCGTCGGCCCGCACTTCGACAACTACGATGTGTTCCTGCTGCa                                                                                                                                   >  1:558159/1‑150 (MQ=255)
                                      ggAAGTGGCCGAACTGCTGGAAAACTTCCGCTTCCTGCCCAGCTGGCGTATCGATGACGTGATGATCAGCTTCGCCGCCCCGGGTGGCAGCGTCGGCCCGCACTTCGACAACTACGATGTGTTCCTGCTGCAAGGCCACGGCCAGCGCaa                                                                                                                 >  2:421594/1‑150 (MQ=255)
                                                      ctgGAAAACTTCCGCTTCCTGCCCAGCTGGCGTATCGATGACGTGATGATCAGCTTCGCCGCCCCGGGTGGCAGCGTCGGCCCGCACTTCGACAACTACGATGTGTTCCTGCTGCAAGGCCACGGCCAGCGCAACTGGAAGATCGGCCAg                                                                                                 >  1:250039/1‑150 (MQ=255)
                                                         gAAAACTTCCGCTTCCTGCCCAGCTGGCGTATCGATGACGTGATGATCAGCCTCGCCGCCCCGGGTGGCAGCGTCGGCGCGGACTTCGACAACTACCAGGTGCTCCTGCTCCAAGGCCACGGCCAGCGCAACTGGAAGGTCGACCAGAtg                                                                                              >  2:154070/1‑150 (MQ=255)
                                                               ttCCGCTTCCTGCCCAGCTGGCGTATCGATGACGTGATGATCAGCTTCGCCGCCCCGGGTGGCAGCGTCGGCCCGCACTTCGACAACTACGATGTGTTCctgc                                                                                                                                       <  1:631836/103‑1 (MQ=255)
                                                               ttCCGCTTCCTGCCCAGCTGGCGTATCGATGACGTGATGATCAGCTTCGCCGCCCCGGGTGGCAGCGTCGGCCCGCACTTCGACAACTACGATGTGTTCCTGCTGCAAGGCCACGGCCAGCGCAACTGGAAGATCGGCCAGATGTgcagc                                                                                        >  1:504240/1‑150 (MQ=255)
                                                               ttCCGCTTCCTGCACAGCTGGCGTATCGATGACGTGATGATCAGCTTCGCCGCCCCGGGTGGCAGCGTCGGCCCGCACTTCGACAACTACGATGTGTTCctgc                                                                                                                                       >  2:631836/1‑103 (MQ=255)
                                                                 ccGCTTCCTGCCCAGCTGGCGTATCGATGACGTGATGATCAGCTTCGCCGCCCCGGGTGGCAGCGTCGGCCCGCACTTCGACAACTACGATGTGTTCCTGCTGCAAGGCCACCGCCAGCGCAACTGGAAGATCGGCCAGATGTGcagcag                                                                                      <  1:646386/150‑1 (MQ=255)
                                                                      tCCTGCCCAGCTGGAGTATCGATGACGTGATGATCAGCTTCGCCGCCCCGGGTGGCAGCGTCGGCCCGCACTTCGACAACTACGATGTGTTCCTGCTGCAAGGCCACGGCCAGCGCAACTGGAAGATCGGCCAGATGTGCAGCAGTGACa                                                                                 <  1:281596/150‑1 (MQ=255)
                                                                        cTGCCCAGCTGGCGTATCGATGACGTGATGATCAGCTTCGCCGCCCCGGGTGGCAGCGTCGGCCCGCACTTCGACAACTACGATGTGTTCCTGCTGCAAGGCCACGGCCAGCGCAACTGGAAGATCGGCCAGATGTGCAGCAGTGACAGt                                                                               <  1:421594/150‑1 (MQ=255)
                                                                         tGCCCAGCTGGCGTATCGATGACGTGATGATCAGCTTCGCCGCCCCGGGTGGCAGCGTCGGCCCGCACTTCGACAACTACGATGTGTTCCTGCTGCAAGGCCACGGCCAGCGCAACTGGAAGATCg                                                                                                      <  2:497907/126‑1 (MQ=255)
                                                                         tGCCCAGCTGGCGTATCGATGACGTGATGATCAGCTTCGCCGCCCCGGGTGGCAGCGTCGGCCCGCACTTCGACAACTACGATGTGTTCCTGCTGCAAGGCCACGGCCAGCGCAACTGGAAGATCg                                                                                                      >  1:497907/1‑126 (MQ=255)
                                                                         tGCCCAGCTGGCGTATCGATGACGTGATGATCAGCTTCGCCGCCCCGGGTGGCAGCGTCGGCCCGCACTTCGACAACTACGATGTGTTCCTGCTGCAAGGCCACGGCCAGCGCAACTGGAAGATCGGCCAGATGTGCAGCAGTGACAGTc                                                                              >  1:777119/1‑150 (MQ=255)
                                                                           cccAGCTGGCGTATCGATGACGTGATGATCAGCTTCGCCGCCCCGGGTGGCAGCGTCGGCCCGCACTTCGACAACTACGATGTGTTCCTGCTGCAAGGCCACGGCCAGCGCAACTGGAAGATCGGCCAGATGTGCAGCAGTGACAGTCCg                                                                            >  2:414507/1‑150 (MQ=255)
                                                                           cccAGCTGGCGTATCGATGACGTGATGATCAGCTTCGCCGCCCCGGGTGGCAGCGTCGGCCCGCACTTCGACAAATACGATGTGTTCCTGCTGCAAGGCCACGGCCAGCGCAACTGGAAGATc                                                                                                       >  1:297949/1‑123 (MQ=255)
                                                                           cccAGCTGGCGTATCGATGACGTGATGATCAGCTTCGCCGCCCCGGGTGGCAGCGTCGGCCCGCACTTCGACAAATACGATGTGTTCCTGCTGCAAGGCCACGGCCAGCGCAACTGGAAGATc                                                                                                       <  2:297949/123‑1 (MQ=255)
                                                                                   gCGTATCGATGACGTGATGATCAGCTTCGCCGCCCCGGGTGGCAGCGTCGGCCCGCACTTCGACAACTACGATGTGTTCCTGCTGCAAGGCCACGGCCAGCGCAACTGGAAGATCGGCCAGATGTGCAGCAGTGACAGTCCGCTGCTGGa                                                                    <  2:282785/150‑1 (MQ=255)
                                                                                          gATGACGTGATGATCAGCTTCGCCGCCCCGGGTGGCAGCGTCGGCCCGCACTTCGACAACTACGATGTGTTCCTGCTGCAAGGCCACGGCCAGCGCAACTGGAAGATCGGCCAGATGTGCAGCAGTGACAGTCCGCTGCTGGAGCACGcc                                                             >  2:51047/1‑150 (MQ=255)
                                                                                          gATGACGTGATGATCAGCTTCGCCGCCCCGGGTGGCAGCGTCGGCCCGCACTTCGACAACTACGATGTGTTCCTGCTGCAAGGCCACGGCCAGCGCAACTGGAAGATCGGCCAGATGTGCAGCAGTGACAGTCCGCTGCTGGAGCACGcc                                                             >  2:50970/1‑150 (MQ=255)
                                                                                               cGTGATGATCAGCTTCGCCGCCCCGGGTGGCAGCGTCGGCCCGCACTTCGACAACTACGATGTGTTCCTGCTGCAAGGCCACGGCCAGCGCAACTGGAAGATCGGCCAGATGTGCAGCAGTGACAGTCCGCTGCTGGAGCACGCCGACCt                                                        <  2:250039/150‑1 (MQ=255)
                                                                                                gTGATGATCAGCTTCGCCGCCCCGGGTGGCAGCGTCGGCCCGCACTTCGACAACTACGATGTGTTCCTGCTGCAAGGCCACGGCCAGCGCAACTGGAAGATCGGCCAGATGTGCAGCAGTGACAGTCCGCTGCTGGAGCACGCCGACCTg                                                       >  2:593304/1‑150 (MQ=255)
                                                                                                gTGATGATCAGCTTCGCCGCCCCGGGTGGCAGCGTCGGCCCGCACTTCGACAACTACGATGTGTTCCTGCTGCAAGGCCACGGCCAGCGCAACTGGAAGATCGGCCAGATGTGCAGCAGTGACAGTCCGCTGCTGGAGCACGCCGACCTg                                                       >  2:277124/1‑150 (MQ=255)
                                                                                                 tgatgaTCAGCTTCGCCGCCCCGGGTGGCAGCGTCGGCCCGCACTTCGACAAATACGATGTGTTCCTGCTGCAAGGCCACGGCCAGCGCAACTGGAAGATCGGCCAGATGTGCAGCAGTGACAGTCCGCTGCTGGAGCACGCCGACCTGc                                                      <  2:756903/150‑1 (MQ=255)
                                                                                                     gatCAGCTTCGCCGCCCCGGGTGGCAGCGTCGGCCCGCACTTCGACAACTAc                                                                                                                                                    >  1:445819/1‑52 (MQ=255)
                                                                                                     gatCAGCTTCGCCGCCCCGGGTGGCAGCGTCGGCCCGCACTTCGACAACTAc                                                                                                                                                    <  2:445819/52‑1 (MQ=255)
                                                                                                          gCTTCGCCGCCCCGGGTGGCAGCGTCGGCCCGCACTTCGACAACTACGATGTGTTCCTGCTGCAAGGCCACGGCCAGCGCAACTGGAAGATCGGCCAGATGTGCAGCAGTGACAGTCCGCTGCTGGAGCACGCCGACCTGCGTATCCTgg                                             >  1:390500/1‑150 (MQ=255)
                                                                                                                                 gTCGGCCCGCACTTCGACAACTACGATGTGTTCCTGCTGCAAGGCCACGGCCAGCGCAACTGGAAGATCGGCCAGATGTGCAGCAGTGACAGTCCGCTGCTGGAGCACGCCGACCTGCGTATCCTGGCCGAGTTCGAGCAGAGCGGCGAg                      <  1:453139/150‑1 (MQ=255)
                                                                                                                                       ccGCACTTCGACAACTACGATGTGTTCCTGCTGCAAGGCCACGGCCAGCGCAACTGGAAGATCGGCCAGATGTGCAGCAGTGACAGTCCGCTGCTGGAGCACGCCGACCTGCGTATCCTGGCCGAGTTCGATCAGAGCGGCGAGt                     >  2:130218/1‑145 (MQ=255)
                                                                                                                                       ccGCACTTCGACAACTACGATGTGTTCCTGCTGCAAGGCCACGGCCAGCGCAACTGGAAGATCGGCCAGATGTGCAGCAGTGACAGTCCGCTGCTGGAGCACGCCGACCTGCGTATCCTGGCCGAGTTCGAGCAGAGCGGCGAGt                     <  1:130218/145‑1 (MQ=255)
                                                                                                                                          cACTTCGACAACTACTATGTGTTTCTGCTGTAAGGCCACGGCCAGCGCAACTTGAAGATCGGCCAGATGTGCAGCAGTGACAGTCCGCTGCTCGAGCACTTTGACCAGCGGTTTGTGGCGGGCTTCGTTCATAGCGGCGAGTGGACCCTg             <  1:100112/150‑1 (MQ=255)
                                                                                                                                           aCTTCGACAACTACGATGTGTTCCTGCTGCAAGGCCACGGCCAGCGCAACTGGAAGATCGGCCAGATGTGCAGCAGTGACAGTCCGctgctg                                                                      <  1:790964/92‑1 (MQ=255)
                                                                                                                                           aCTTCGACAACTACGATGTGTTCCTGCTGCAAGGCCACGGCCAGCGCAACTGGAAGATCGGCCAGATGTGCAGCAGTGACAGTCCGctgctg                                                                      >  2:790964/1‑92 (MQ=255)
                                                                                                                                            cTTCGACAACTACGATGTGTTCCTGCTGCAAGGCCACGGCCAGCGCAACTGGAAGATCGGCCAGATGTGCAGCAGTGACAGTCCGCTGCTGGAGCACGCCGACCTGCGTATCCTGGCCGAGTTCGAGCAGAGCGGCGAGTGGACCCTGGa           <  1:414507/150‑1 (MQ=255)
                                                                                                                                                     cTACGATGTGTTCCTGCTGCAAGGCCACGGCCAGCGCAACTGGAAGATCGGCCAGATGTGCAGCAGTGACAGTCCGCTGCTGGAGCACGCCGACCTGCGTATCCTGGCCGAGTTCGAGCAGAGCGGCGAGTGGACCCTGGAGCCCGGCGa  <  1:51047/150‑1 (MQ=255)
                                                                                                                                                     cTACGATGTGTTCCTGCTGCAAGGCCACGGCCAGCGCAACTGGAAGATCGGCCAGATGTGCAGCAGTGACAGTCCGCTGCTGGAGCACGCCGACCTGCGTATCCTGGCCGAGTTCGAGCAGAGCGGCGAGTGGACCCTGGAGCCCGGCGa  <  1:50970/150‑1 (MQ=255)
                                                                                                                                                     cTACGATGTGTTCCTGCTGCAAGGCCACGGCCAGCGCAACTGGAAGATCGGCCAGATGTGCAGCAGTGACAGTCCGCTGCTGGAGCACGCCGACCTGCGTATCCTGGCCGAGTTCGAGCAGAGCGGCGAGTGGACCCTGGAGCCCGGCGa  >  2:745118/1‑150 (MQ=255)
                                                                                                                                                     |                                                                                                                                                     
TGGACCCTGCTGGTGCAGGCCGTTGACCAGTTCGTTCCGGAAGTGGCCGAACTGCTGGAAAACTTCCGCTTCCTGCCCAGCTGGCGTATCGATGACGTGATGATCAGCTTCGCCGCCCCGGGTGGCAGCGTCGGCCCGCACTTCGACAACTACGATGTGTTCCTGCTGCAAGGCCACGGCCAGCGCAACTGGAAGATCGGCCAGATGTGCAGCAGTGACAGTCCGCTGCTGGAGCACGCCGACCTGCGTATCCTGGCCGAGTTCGAGCAGAGCGGCGAGTGGACCCTGGAGCCCGGCGA  >  NC_002947/4528712‑4529010

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