Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 5,677,577 C→A 23.5% intergenic (+38/‑181) PP_4979 → / → rhlE‑II substrate‑binding protein/ATP‑dependent RNA helicase

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0029475,677,5770CA23.5% 34.0 / 5.1 17intergenic (+38/‑181)PP_4979/rhlE‑IIsubstrate‑binding protein/ATP‑dependent RNA helicase
Reads supporting (aligned to +/- strand):  ref base C (8/5);  new base A (2/2);  total (10/7)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 5.60e-01

GCGATCAGCAACCGTTATACAACCGACACTGGAAACATTCTCAACGCCGTTGAGCAGCAGGAAAGCAGCACAGCGCGATAGCTCTGTTATACTCGGGCCTTCCCGCCCGGCTCACGCCCGGACGCTCGGCCTCGCAACAGGCATCCCGATCGGTACCGACGCCCCTTCGCGTCCAGCCTCCGTTTCCCGGATGTGCAGTGAAAGCCCAGCTGGACCGGACGCGATCGCATTTCACCGATGC  >  NC_002947/5677459‑5677699
                                                                                                                      |                                                                                                                          
gCGATCAGCAACCGTTATACAACCGACACTGGAAACATTCTCAACGCCGTTGAGCAGCAGGAAAGCAGCACAGCGCGATAGCTCTGTTATACTCGGGCCTTCCCGCCCGGCTCACGCCCGGACGCTCGGCCTCGCa                                                                                                           <  1:301633/136‑1 (MQ=255)
         aaCCGTTATACAACCGACACTGGAAACATTCTCAACGCCGTTGAGCAGCAGGAAAGCAGCACAGCGCGATAGCTCTGTTATACTCGGGCCTTCCCGCCCGGCTCACGCCCGGACGCTCGGCCTCGCAACAGGCATc                                                                                                  >  1:137454/1‑136 (MQ=255)
         aaCCGTTATACAACCGACACTGGAAACATTCTCAACGCCGTTGAGCAGCAGGAAAGCAGCACAGCGCGATAGCTCTGTTATACTCGGGCCTTCCCGCCCGGCTCACGCCCGGACGCTCGGCCTCGCAACAGGCATc                                                                                                  >  2:190466/1‑136 (MQ=255)
           ccGTTATACAACCGACACTGGAAACATTCTCAACGCCGTTGAGCAGCAGGAAAGCAGCACAGCGCGATAGCTCTGTTATACTCGGGCCTTCCCGCCCGGCTCACGCCCGGACGCTCGGCCTCGCAACAGGCATccc                                                                                                <  1:394936/136‑1 (MQ=255)
                         acacTGGAAACATTCTCAACGCCGTTGAGCAGCAGGAAAGCAGCACAGCGCGATAGCTCTGTTATACTCGGGCCTTCCCGCCCGGCTCACGCCCGGACGCTCGGCCTCGCAACAGGCATCCCGATCGGTACCGACg                                                                                  >  1:149590/1‑136 (MQ=255)
                         acacTGGAAACATTCTCAACGCCGTTGAGCAGCAGGAAAGCAGCACAGCGCGATAGCTCTGTTATACTCGGGCCTTCCCGCCCGGCTCACGCCCGGACGCTCGGCCTCGCAACAGGCATCCCGATCGGTACCGACg                                                                                  >  1:24939/1‑136 (MQ=255)
                                 aaCATTCTCAACGCCGTTGAGCAGCAGGAAAGCAGCACAGCGCGATAGCTCTGTTATACTCGGGCCTTCCCGCCCGGCTCACGCCCGGACGCTCGGCCTCGCAACAGGCATCCCGATCGGTACCGACGCCCCTTcg                                                                          >  2:239049/1‑136 (MQ=255)
                                            cGCCGTTGAGCAGCAGGAAAGCAGCACAGCGCGATAGCTCTGTTATACTCGGGCCTTCCCGCCCGGCTCACGCCAGGACGCTCGGCCTCGCAACAGGCATCCCGATCGGTACCGACGCCCCTTCGCGTCCAGCCTc                                                               >  2:37875/1‑136 (MQ=255)
                                            cGCCGTTGAGCAGCAGGAAAGCAGCACAGCGCGATAGCTCTGTTATACTCGGGCCTTCCCGCCCGGCTCACGCCAGGACGCTCGGCCTCGCAACAGGCATCCCGATCGGTACCGACGCCCCTTCGCGTCCAGCCTc                                                               >  2:339834/1‑136 (MQ=255)
                                                   gAGCAGCAGGAAAGCAGCACAGCGCGATAGCTCTGTTATACTCGGGCCTTCCCGCCCGGCTCACGCCCGGACGCTCGGCCTCGCAACAGGCATCCCGATCGGTACCGACGCCCCTTCGCGTCCAGCCTCCGTTTcc                                                        <  2:149590/136‑1 (MQ=255)
                                                            gAAAGCAGCACAGCGCGATAGCTCTGTTATACTCGGGCCTTCCCGCCCGGCTCACGCCCGGACGCTCGGCCTCGCAACAGGCATCCCGATCGGTACCGACGCCCCTTCGCGTCCAGCCTCCGTTTCCCGGATGTGc                                               <  1:248224/136‑1 (MQ=255)
                                                                                                 ccTTCCCGCCCGGCTCACGCCCGGACGCTCGGCCTCGCAACAGGCATCCCGATCGGTACCGACGCCCCTTCGCGTCCAGCCTCCGTTTCCCGGATGTGCAGTGAAAGCCCAGCTGGACCGGACGCGATCGCATTTc          >  2:279910/1‑136 (MQ=255)
                                                                                                 ccTTCCCGCCCGGCTCACGCCCGGACGCTCGGCCTCGCAACAGGCATCCCGATCGGTACCGACGCCCCTTCGCGTCCAGCCTCCGTTTCCCGGATGTGCAGTGAAAGCCCAGCTGGACCGGACGCGATCGCATTTc          >  1:128828/1‑136 (MQ=255)
                                                                                                   ttCCCGCCCGGCTCACGCCCGGACGCTCGGCCTCGCAACAGGCATCCCGATCGGTACCGACGCCCCTTCGCGTCCAGCCTCCGTTTCCCGGATGTGCAGTGAAAGCCCAGCTGGACCGGACGCGATCGCATTTCAc        >  1:180156/1‑136 (MQ=255)
                                                                                                         cccgGCTCACGCCCGGACGCTCGGCCTCGCAACAGGCATCCCGATCGGTACCGACGCCCCTTCGCGTCCAGCCTCCGTTTCCCGGATGTGCAGTGAAAGCCCAGCTGGACCGGACGCGATCGCATTTCACCGATGc  <  1:173292/136‑1 (MQ=255)
                                                                                                         cccgGCTCACGCCAGGACGCTCGGCCTCGCAACAGGCATCCCGATCGGTACCGACGCCCCTTCGCGTCCAGCCTCCGTTTCCCGGATGTGCAGTGAAAGCCCAGCTGGACCGGACGCGATCGCATTTCACCGATGc  <  1:37875/136‑1 (MQ=255)
                                                                                                         cccgGCTCACGCCAGGACGCTCGGCCTCGCAACAGGCATCCCGATCGGTACCGACGCCCCTTCGCGTCCAGCCTCCGTTTCCCGGATGTGCAGTGAAAGCCCAGCTGGACCGGACGCGATCGCATTTCACCGATGc  <  1:339834/136‑1 (MQ=255)
                                                                                                                      |                                                                                                                          
GCGATCAGCAACCGTTATACAACCGACACTGGAAACATTCTCAACGCCGTTGAGCAGCAGGAAAGCAGCACAGCGCGATAGCTCTGTTATACTCGGGCCTTCCCGCCCGGCTCACGCCCGGACGCTCGGCCTCGCAACAGGCATCCCGATCGGTACCGACGCCCCTTCGCGTCCAGCCTCCGTTTCCCGGATGTGCAGTGAAAGCCCAGCTGGACCGGACGCGATCGCATTTCACCGATGC  >  NC_002947/5677459‑5677699

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