Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 1,499,506:1 +C 100% intergenic (+83/‑81) trpS → / → zapE tryptophan‑‑tRNA ligase/nucleoside triphosphate hydrolase domain‑containing protein

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0029471,499,5061.C95.8% 61.2 / ‑2.9 24intergenic (+83/‑81)trpS/zapEtryptophan‑‑tRNA ligase/nucleoside triphosphate hydrolase domain‑containing protein
Reads supporting (aligned to +/- strand):  ref base . (1/0);  new base C (8/15);  total (9/15)
Fisher's exact test for biased strand distribution p-value = 3.75e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 8.37e-01
Rejected as polymorphism: E-value score below prediction cutoff.
Rejected as polymorphism: Frequency below/above cutoff threshold.
Rejected as polymorphism: Variant not supported by required number of reads on each strand.

TGAAGCCCGTGATGCGGCCATTCAGAGCCTGAACGAAGCCTTGGCACCACAGCAGGACTGAGGCTGCCAAACGGTGGGAGCGGGCGTGCCCGCGAAACAGGCGCCGCGGTGCATGGCCACGGCTTTGCCGGTGTTTCGGGCACGCCGCTCCCACAGACTGATTGCCATTAAGCGGGCCTATCGCTACAGTGACGGCCCGTTTTTTGTTGCCTCGCTAACGAATCATGACTCCCTTAGAACGCTATCAAGCAGATC  >  NC_002947/1499363‑1499617
                                                                                                                                                  |                                                                                                                
tGAAGCCCGTGATGCGGCCATTCAGAGCCTGAACGAAGCCTTGGCACCACAGCAGGACTGAGGCTGCCAAACGGTGGGAGCGGGCGTGCCCGCGAAACACGCGCCGCGGTGCATGGCACCGGCTTTGCCGGTGTTAGCGGGCCcgcccg                                                                                                              >  5:164114/1‑147 (MQ=25)
                gCCTTTCAGAGCCTGAACGAAGCCTTGGCACCAGAGCAGGACTGAGGCTGCCAAACGGTGGGAGCGGGCGGGCCCGCGAAACAGGCGCCGCGGTGCATGGCACCGGGTTTGCCGGTGTTCGCGGGCACGCCCGCTCCCACAGACTGAtt                                                                                              <  8:356912/149‑1 (MQ=21)
                gCCATTCAGAGCCTGAACGAAGCCTTGGCACCACAGCAGGACTGAGGCTGCCAAACGGTGGGAGCGGGCGTGCCCGCGAAACAGGCGCCGCGGTGCATGGCACCGGGTTTGCCGGTGTTAGCGGGCCCGCCCGCTCCCCCAGACTGAtt                                                                                              >  6:302248/1‑149 (MQ=25)
                gCCATTCAGAGCCTGAACGAAGCCTTGGCACCACAGCAGGACTGAGGCTGCCAAACGGTGGGAGCGGGCGTGCCCGCGAAACAGGCGCCGCGGTGCATGACACCGGCGTTGCCGGTGTTCGCGGGCACGCCCGCTCCCACAGACTGAtt                                                                                              <  2:77009/149‑1 (MQ=21)
                    ttCAGAGCCTGAACGAAGCCTTGGCACCACAGCAGGACTGAGGCTGCCAAACGGTGGGAGCGGGCGTGCCCGCGAAACAGGCGCCGCGGTGCATGGCACCGGCTTTGCCGGTGTTCGCGGGCACGCCCGCTCCCACAGACTGATTGCCa                                                                                          <  3:283052/149‑1 (MQ=21)
                      cAGAGCCTGAACGAAGCCTTGGCACCACAGCAGGACTGAGGCTGCCAAACGGTGGGAGCGGGCGTGCCCGCGAAACAGGCGCCGCGGTGCATGGCACCGGCTTTGCCGGTGTTCGCGGGCACGCCCGCTCCCACAGACTGATTGCCAt                                                                                         <  2:31751/148‑1 (MQ=21)
                         agCCTGAACGAAGCCTTGGCACCACAGCAGGACTGAGGCTGCCAAACGGTGGGAGCGGGCGTGCCCGCGAAACAGGCGCCGCGGTGCATGGCACCGGCTTTGCCGGTGTTCGCGGGCACGCCCGCTCCCACAGACTGATTGCCATTAAg                                                                                     <  1:240597/149‑1 (MQ=21)
                             tGAACGAAGCCTTGGCACCACAGCAGGACTGAGGCTGCCAAACGGTGGGAGCGGGCGTGCCCGCGAAACAGGCGCCGCGGTGCATGGCACCGGCTTTGCCGGTGTTCGCGGGCACGCCCGCTCCCCCAGACTGATTGCCATTAAGCggg                                                                                 >  7:281275/1‑149 (MQ=21)
                                        ttGGCACCACAGCAGGACTGAGGCTGCCAAACGGTGGGAGCGGGCGTGCCCGCGAAACCGGCGCCGCGGTGCATGGCACCCGCGTTGCCGGTGTTAGCGGGCACGCCCGCTCCCACAGACTGATTGCCATTAAGCGGGCCGATCGCTa                                                                       >  2:329290/1‑148 (MQ=21)
                                        ttGGCACCACAGCAGGACTGAGGCTGCCAAACGGTGGGAGCGGGCGTGCCCGCGAAACAGGCGCCGCGGTGCATGGCACCGGCTTTGCCGGTGTTCGCGGGCCCGCCCGCTCCCCCAGACTGATTTCCAATAAGCGGGCCTATCGCTa                                                                       >  3:79953/1‑148 (MQ=21)
                                              ccacAGCAGGACTGAGGCTGCCAAACGGTGGGAGCGGGCGTGCCCGCGAAACAGGCGCCGCGGTGCATGGCACCGGCTTTGCCGGTGTTCGCGGGCACGCCCGCTCCCACAGACTGATTGCCATTAAGCGGGCCTATCGCTACAGTGAc                                                                <  4:245245/149‑1 (MQ=21)
                                                             ggCTGCCAAACGGTGGGAGCGGGCGTGCCCGCGAAACATGGTCCGCTGTGCATGGCACCGGCTTTGCCGGTGTTCGCGGGCACGCCCGCTCCCACAGACTGATTGCCATTAAGCGGGCCTATCGCTACAGTTACGGCCCGTTTTttgtt                                                 <  6:194891/149‑1 (MQ=21)
                                                                  ccAAACGGTGGGAGCGGGCGTGGCCGCGGAACAGGCGCCGCGGTGGATGGCACCGGCGTTGGCGGGGTTAGCGGGGCCGCCCGCGCCCCCAGCCGGGTTGCCATTAAGCGGGGCGATCGGGACAGTGACGGGCCGGTGTTTGTTGccgc                                            >  6:108785/1‑147 (MQ=255)
                                                                    aaaCGGGGGGAGCGGGCGTGCCCGCGAAACAGGCGCCGCGGTGCATGGCACCGGCTTTGCCGGTGTTCGCGGGCACGCCCGCTCCCACAGACTGATTGCCATTAAGCGGGCCTATCGCTACAGTGACGGCCCGTTTTTTGTTGCCTCGc                                          <  7:190922/149‑1 (MQ=21)
                                                                     aaCGGTGGGAGCGGGCGTGCCCGCGAAACAGGCGCCGCGGTGCATGGCACCGGCTTTGCCGGTGTTCGCGGGCACGCCCGCTCCCACAGACTGATTGCCATTAAGCGGGCCTATCGCTACAGTGACGGCCCGTTTTTTGTTGCCTCGc                                          >  1:269393/1‑148 (MQ=21)
                                                                        ggTGGGAGGGGGCGGGCCCGCTAAACAGGCGCCGCGGGGCATGGCACGGGCTTTGCCGGTGTTCGCGGGCACGCCCGCTCCCACAGACTGATTGCCATTAAGCGGGCCTATCGCTACAGTGACGGCCCGTTTTTTGTTGCCTCGCTAAc                                      <  3:381698/149‑1 (MQ=25)
                                                                            ggAGCGGGCGTGCCCGCGAAACAGGCGCCGCGGTGCATGGCACCGGCTTTGCCGGTGTTCGCGGGCACGCCCGCTCCCACAGACTGATTGCCATTAAGCGGGCCTATCGCTACAGTGACGGCCCGTTTTTTGTTGCCTCGCTAACGAAt                                  <  5:302248/149‑1 (MQ=25)
                                                                             gcgCGGGCGTGCCCGCTAAACCGGCGCGGCGGCGCAGGGCACCGGCTTTGCCGGTGTTCGCGGGCACGCCCGCTCCCACAGACTGATTGCCATTAAGCGGGCCTATCGCTACAGTGACGGCCCGTTTTTTGTTGCCTCGCTAACGAATc                                 <  8:289551/147‑1 (MQ=25)
                                                                                 gggCGTGCCCGCGAAACAGGCGCCGCGGTGCATGGCACCTGCTTTGCCGGTGTTCGCGGGCACGCCCGCTCCCACAGACTGATTGCCATTAAGCGGGCCTATCGCTACAGTGACGGCCCGTTTTTTGTTGCCTCGCTAACGAATCATGa                             <  7:297403/149‑1 (MQ=37)
                                                                                           gcgAAACAGGCGCCGCGGTGCATGGCACCGGCTTTGCCGGTGTTCGCGGGCACGCCGGCTCCCACAGACTGATTGCCATTAAGCGGGCCTATCGCTACAGTGACGGCACCTTTTTTTTTTCCTCGCTAACTAATCATGACTCCCTTAg                    >  5:328470/1‑148 (MQ=25)
                                                                                              aaaCAGGCGCCGCGGTGCAGGCCACGGGGTTTGCCGTTGTTCGCGGGCACGCCCGCTCCCACAGACTGATTGCCATTAAGCGGGCCTATCGCTACAGTGACGGCCCGTTTTTTGTTGCCTCGCTAACGAATCATGACTCCCTTAGAACg                <  5:15275/149‑1 (MQ=255)
                                                                                                aCAGGCGCCGCGGTGCATGGCACCGGCTTTGCCGGTGTTCGCGGGCACGCCCGCTCCCACAGACTGATTGCCATTAAGCGGGCCTATCGCTACAGTGACGGCCCGTTTTTTGTTGCCTCGCTAACGAATCATGACTCCCTTAGAACGCt              <  6:164114/149‑1 (MQ=37)
                                                                                                  aGGCGCCGCGGTGCATGGCACCGGCTTTGCCGGTGTTCGCGGGCACGCCCGCTCCCACAGACTGATTGCCATTAAGCGGGCCTATCGCTACAGTGACGGCCCGTTTTTTGTTGCCTCGCTAACGAATCATGACTCCCTTAGAACGCTAt            >  3:43523/1‑149 (MQ=37)
                                                                                                      gccgcGGTGCATGGCACCGGCTTTGCCGGTGTTCGCGGGCACGCCCGCTCCCACAGACTGATTGCCATTAAGCGGGCCTATCGCTACAGTGACGGCCCGTTTTTTGTTGCCTCGCTAACGAATCATGACTCCCTTAGAACGCTATCAAg        <  7:41217/149‑1 (MQ=37)
                                                                                                            gTGCATGGCACCGGCTTTGCCGGTGTTCGCGGGCACGCCCGCTCCCACAGACTGATTGCCATTAAGCGGGCCTATCGCTACAGTGACGGCCCGTTTTTTGTTGCCTCGCTAACGAATCATGACTCCCTTAGAACGCTATCAACCAGATc  >  1:350054/1‑149 (MQ=37)
                                                                                                                                                  |                                                                                                                
TGAAGCCCGTGATGCGGCCATTCAGAGCCTGAACGAAGCCTTGGCACCACAGCAGGACTGAGGCTGCCAAACGGTGGGAGCGGGCGTGCCCGCGAAACAGGCGCCGCGGTGCATGGCCACGGCTTTGCCGGTGTTTCGGGCACGCCGCTCCCACAGACTGATTGCCATTAAGCGGGCCTATCGCTACAGTGACGGCCCGTTTTTTGTTGCCTCGCTAACGAATCATGACTCCCTTAGAACGCTATCAAGCAGATC  >  NC_002947/1499363‑1499617

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