Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 867,872 G→A 20.0% E94E (GAG→GAA PP_0750 → hypothetical protein

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_002947867,8720GA20.0% 39.6 / 5.0 20E94E (GAG→GAAPP_0750hypothetical protein
Reads supporting (aligned to +/- strand):  ref base G (5/11);  new base A (2/2);  total (7/13)
Fisher's exact test for biased strand distribution p-value = 5.87e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 4.07e-01

TGCGCGTGGACGACTTCGAGCGCTTTGTAGGCTTCTTCGTCGAAGCTGGCCTGGACCTGGACGGCAAGGACGAGCACGGCAAGACCTTCGTCGAGCAGATTGCCGACCAGCGCAATGCCCCTGAGTACATCGAGATCATCGATAAAGCCCGTGGCTGATCTGTAGCCAGTGCCGGCCTCTTCGCGGGCGCGCCCGCTCCCACAGGAGTCGCGCCGCACCTGCGGGCATTACCGAACCTGTG  >  NC_002947/867748‑867988
                                                                                                                            |                                                                                                                    
tGCGCGTGGACGACTTCGAGCGCTTTGTAGGCTTCTTCGTCGAAGCTGGCCTGGACCTGGACGGCAAGGACGAGCACGGCAAGACCTTCGTCGAGCAGATTGCCGACCAGCGCAATGCCCCTGAGTACATCgagtt                                                                                                           >  1:486828/1‑134 (MQ=255)
     gTGGACGACTTCGAGCGCTTTGTAGGCTTCTTCGTCGAAGCTGGCCTGGACCTGGACGGCAAGGACGAGCACGGCAAGACCTTCGTCGAGCAGATTGCCGACCAGCGCAATGCCCCTGAGTACATCGAGATCATCg                                                                                                      <  2:196202/136‑1 (MQ=255)
     gTGGACGACTTCGAGCGCTTTGTAGGCTTCTTCGTCGAAGCTGGCCTGGACCTGGACGGCAAGGACGAGCACGGCAAGACCTTCGTCGAGCAGATTGCCGACCAGCGCAATGCCCCTGAGTACATCGAGATCATCg                                                                                                      >  2:203084/1‑136 (MQ=255)
                 gAGCGCTTTGTAGGCTTCTTCGTCGAAGCTGGCCTGGACCTGGACGGCAAGGACGAGCACGGCAAGACCTTCGTCGAGCAGATTGCCGACCAGCGCAATGCCCCTGAGTACATCGAGATCATCGATAAAGCCCGTg                                                                                          <  1:350146/136‑1 (MQ=255)
                 gAGCGCTTTGTAGGCTTCTTCGTCGAAGCTGGCCTGGACCTGGACGGCAAGGACGAGCACGGCAAGACCTTCGTCGAGCAGATTGCCGACCAGCGCAATGCCCCTGAGTACATCGAGATCATCGATAAAGCCCGTg                                                                                          <  1:37816/136‑1 (MQ=255)
                                    tcgtcgAAGCTGGCCTGGACCTGGACGGCAAGGACGAGCACGGCAAGACCTTCGTCGAGCAGATTGCCGACCAGCGCAATGCCCCTGAGTACATCGAGATCATCGATAAAGCCCGTGGCTGATCTGTAGCCAGTGc                                                                       >  1:358093/1‑136 (MQ=255)
                                    tcgtcgAAGCTGGCCTGGACCTGGACGGCAAGGACGAGCACGGCAAGACCTTCGTCGAGCAGATTGCCGACCAGCGCAATGCCCCTGAGTACATCGAGATCATCGATAAAGCCCGTGGCTGATCTGTAGCCAGTGc                                                                       >  2:103202/1‑136 (MQ=255)
                                     cgtcgAAGCTGGCCTGGACCTGGACGGCAAGGACGAGCACGGCAAGACCTTCGTCGAGCAGATTGCCGACCAGCGCAATGCCCCTGAGTACATCGAGATCATCGATAAAGCCCGTGGCTGATCTGTAGCCAGTGcc                                                                      <  2:358093/136‑1 (MQ=255)
                                                                   ggACGAGCACGGCAAGACCTTCGTCGAGCAGATTGCCGACCAGCGCAATGCCCCTGAATACATCgaga                                                                                                            <  1:126883/68‑1 (MQ=255)
                                                                   ggACGAGCACGGCAAGACCTTCGTCGAGCAGATTGCCGACCAGCGCAATGCCCCTGAATACATCgaga                                                                                                            <  1:129611/68‑1 (MQ=255)
                                                                   ggACGAGCACGGCAAGACCTTCGTCGAGCAGATTGCCGACCAGCGCAATGCCCCTGAATACATCgaga                                                                                                            >  2:126883/1‑68 (MQ=255)
                                                                   ggACGAGCACGGCAAGACCTTCGTCGAGCAGATTGCCGACCAGCGCAATGCCCCTGAATACATCgaga                                                                                                            >  2:129611/1‑68 (MQ=255)
                                                                     aCGAGCACGGCAAGACCTTCGTCGAGCAGATTGCCGACCAGCGCAATGCCCCTGAGTACATCGAGATCATCGATAAAGCCCGTGGCTGATCTGTAGCCAGTGCCGGCCTCTTCGCGGGCGCGCCCGCTCCCACAgg                                      <  2:26077/136‑1 (MQ=255)
                                                                           aCGGCAAGACCTTCGTCGAGCAGATTGCCGACCAGCGCAATGCCCCTGAGTACATCGAGATCATCGATAAAGCCCGTGGCTGATCTGTAGCCAGTGCCGGCCTCTTCGCGGGCGCGCCCGCTCCCACAGGAGTcgc                                >  2:111423/1‑136 (MQ=255)
                                                                                             aGCAGATTGCCGACCAGCGCAATGCCCCTGAGTACATCGAGATCATCGATAAAGCCCGTGGCTGATCTGTAGCCAGTGCCGGCCTCTTCGCGGGCGCGCCCGCTCCCACAGGAGTCGCGCCGCACCTGCGGGCAtt              <  1:111423/136‑1 (MQ=255)
                                                                                             aGCAGATTGCCGACCAGCGCAATGCCCCTGAGTACATCGAGATCATCGATAAAGCCCGTGGCTGATCTGTAGCCAGTGCCGGCCTCTTCGCGGGCGCGCCCGCTCCCACAGGAGTCGCGCCGCACCTGCGGGCAtt              <  2:290392/136‑1 (MQ=255)
                                                                                               cAGATTGCCGACCAGCGCAATGCCCCTGAGTACATCGAGTTCATCTATAAAGCCCGTGGCTGATCTGTAGCCAGTGCCGGCCTCTTCGCGGGCGCGCCCGCTCCCACAGGAGTCGCGCCGCACCTGCGGGCATTAc            <  2:486828/136‑1 (MQ=255)
                                                                                               cAGATTGCCGACCAGCGCAATGCCCCTGAGTACATCGAGATCATCGATAAAGCCCGTGGCTGATCTGTAGCCAGTGCCGGCCTCTTCGCGGGCGCGCCCGCTCCCACAGGAGTCGCGCCGCACCTGCGGGCATTAc            <  1:103202/136‑1 (MQ=255)
                                                                                                 ttttGCCGACCAGCGCTATGCCCCTGAGTACTTCGATATCATCGATAAAGCCCGTGGCTGATCTGTTGCCAGTGCCGGCCTCTTCGCGGGCGGGCAGGCTCCCACAGGATTCGCGCCGCACCTGCGGGAATTACCg          <  1:203084/134‑1 (MQ=255)
                                                                                                         cccAGCGCAATGCCCCTGAGTACATCGAGATCATCGATAAAGCCCGTGGCTGATCTGTAGCCAGTGCCGGCCTCTTCGCGGGCGCGCCCGCTCCCACAGGAGTCGCGCCGCACCTGCGGGCATTACCGAACCtgtg  <  1:164329/135‑1 (MQ=255)
                                                                                                                            |                                                                                                                    
TGCGCGTGGACGACTTCGAGCGCTTTGTAGGCTTCTTCGTCGAAGCTGGCCTGGACCTGGACGGCAAGGACGAGCACGGCAAGACCTTCGTCGAGCAGATTGCCGACCAGCGCAATGCCCCTGAGTACATCGAGATCATCGATAAAGCCCGTGGCTGATCTGTAGCCAGTGCCGGCCTCTTCGCGGGCGCGCCCGCTCCCACAGGAGTCGCGCCGCACCTGCGGGCATTACCGAACCTGTG  >  NC_002947/867748‑867988

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