Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 1,499,498 1 bp→CG 100% intergenic (+75/‑89) 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,4971.C94.4% 47.5 / ‑2.9 18intergenic (+74/‑90)trpS/zapEtryptophan‑‑tRNA ligase/nucleoside triphosphate hydrolase domain‑containing protein
Reads supporting (aligned to +/- strand):  ref base . (0/1);  new base C (13/4);  total (13/5)
Fisher's exact test for biased strand distribution p-value = 2.78e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.01e-01
Rejected as polymorphism: E-value score below prediction cutoff.
Rejected as polymorphism: Variant not supported by required number of reads on each strand.
*NC_0029471,499,4980TG95.4% 47.6 / ‑5.9 18intergenic (+75/‑89)trpS/zapEtryptophan‑‑tRNA ligase/nucleoside triphosphate hydrolase domain‑containing protein
Reads supporting (aligned to +/- strand):  ref base T (0/0);  major base G (13/4);  minor base C (0/1);  total (13/5)
Fisher's exact test for biased strand distribution p-value = 2.78e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00
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.

CCGCGAAACAGGCGCCGCGGTGCATGGCCACGGCTTTGCCGGTGTT‑TCGG‑GCACG‑CCGCTCCCACAGACTGATTGCCATTAAGCGGGCCTATCGCTACAGTGACGGCCCGTTTTTTGTTGCCTCGCTAACGAATCATGACTCCCTTAGAACGCTATCAAGCAGATCTGAAACGTCCCGACTTCTTCCATGA  >  NC_002947/1499452‑1499642
                                               ||                                                                                                                                                 
ccGGGAAACTGGCGCGGCGGTGCTTGGCATGGGGTTTGCCGGTGTT‑CGGGTGAAGGCCCGCTCCCACAGACTGATTGCCATTAAGCGGGCCTATCGCTACAGTGACGGCCCGTTTTTTGTTGCCTCGCTAACGAATCATGACTCCCtt                                               <  2:128851/148‑1 (MQ=255)
ccGCGAAACAGGCGCCGCGGTGCATGGCACCGGCTTTGCCGGTGTTCGCGG‑GCACGCCCGCTCCCACAGACTGATTGCCATTAAGCGGGc                                                                                                         <  3:86407/90‑1 (MQ=14)
ccGCGAAACAGGCGCCGCGGTGCATGGCACCGGCTTTGCCGGTGTTCGCGG‑GCACGCCCGCTCCCACAGACTGATTGCCATTAAGCGGGc                                                                                                         >  4:86407/1‑90 (MQ=14)
ccGCGAAACAGGCGCCGCGGTGCATGGCACCGGCTTTGCCGGTGTTCGCGG‑GCACGCCCGCTCCCACAGACTGATTGCCATTAAGCGGGCCTATCGCTACAGTGa                                                                                          >  3:55756/1‑105 (MQ=18)
ccGCGAAACAGGCGCCGCGGTGCATGGCACCGGCTTTGCCGGTGTTCGCGG‑GCACGCCCGCTCCCACAGACTGATTGCCATTAAGCGGGCCTATCGCTACAGTGa                                                                                          <  4:55756/105‑1 (MQ=18)
ccGCGAAACAGGCGCCGCGGTGCATGGCACCGGCTTTGCCGGTGTTCGCGG‑GCACGCCCGCTCCCACAGACTGATTGCCATTAAGCGGGCCTATCGCTACAGTGACGGCCCGTTTTTTGTTGCCTCGCTAACGAATCATGACTCCCt                                                >  2:83045/1‑147 (MQ=37)
ccGCGAAACAGGCGCCGCGGTGCATGGCACCGGCTTTGCCGGTGTTCGCGG‑GCACGCCCGCTCCCACAGACTGATTGCCATTAAGCGGGCCTATCGCTACAGTGACGGCCCGTTTTTTGTTGCCTCGCTAACGAATCATGACTCCCt                                                >  1:230826/1‑147 (MQ=37)
ccGCGAAACAGGCGCCGCGGTGCATGGCACCGGCTTTGCCGGTGTTCGCGG‑GCACGCCCGCTCCCACAGACTGATTGCCATTAAGCGGGCCTATCGCTACAGTGACGGCCCGTTTTTTGTTGCATCGCTAACGAATCATGAATCCCt                                                >  6:162995/1‑147 (MQ=25)
        cAGGCGCCGCGGTGCATGGCACCGGCTTTGCCGGTGTTCGCGG‑GCACGCCCGCTCCCACAGACTGATTGCCATTAAGCGGGCCTATCGCTACAGTGACGGCCCGTTTTTTGTTGCCTCGCTAACGAATCATGACTCCCTTAGAACGCTa                                      >  8:176516/1‑149 (MQ=37)
          ggCGCCGCGGTGCATGGCACCGGCTTTGCCGGTGTTCGCGG‑GCACGCCCGCTCCCACAGACTGATTGCCATTAAGCGGGCCTATCGCTACAGTGACGGCCCGTTTTTTGTTGCCTCGCTAACGAATCATGACTCCCTTAGAACGCTAt                                     >  5:93592/1‑148 (MQ=37)
                     gCATGGCACCGGCTTTGCCGGTGTTCGCGG‑GCACGCCCGCTCCCACAGACTGATTGCCATTAAGCGGGCCTATCGCTACAGTGACGGCCCGTTTTTTGTTGCCTCGCTAACGAATCATGACTCCCTTAGAACGCTATCAAGAAGAACt                          >  3:989/1‑148 (MQ=37)
                          gcaccGGCTTTGCCGGTGTTCGCGG‑GCACGCCCGCTCCCACAGACTGATTGCCATTAAGCGGGCCTATCGCTACAGTGACGGCCCGTTTTTTGTTGCCTCGCTAACGAATCATGACTCCCTTAGAACGCTATCAAGCAGAt                            >  1:256585/5‑141 (MQ=38)
                          gcaccGGCTTTGCCGGTGTTCGCGG‑GCACGCCCGCTCCCACAGACTGATTGCCATTAAGCGGGCCTATCGCTACAGTGACGGCCCGTTTTTTGTTGCCTCGCTAACGAATCATGACTCCCTTAGAACGCTATCAAGCAGAt                            <  2:256585/137‑1 (MQ=38)
                          gcaccGGCTTTGCCGGTGTTCGCGG‑GCACGCCCGCTCCCACAGACTGATTGCCATTAAGCGGGCCTATCGCTACAGTGACGGCCCGTTTTTTGTTGCCTCGCTAACGAATCATGACTCCCTTAGAACGCTATCAAGCAGATCTGAAACg                    >  4:129805/5‑149 (MQ=38)
                            accGGCTTTGCCGGTGTTCGCGG‑GCACGCCCGCTCCCACAGACTGATTGCCATTAAGCGGGCCTATCGCTACAGTGACGGCCCGTTTTTTGTTGCCTCGCTAACGAATCATGACTCCCTTAGAACGCTATCAAGCAGATCTGAAACGTc                  >  4:250292/3‑149 (MQ=38)
                            accGGCTTTGCCGGTGTTCGCGG‑GCACGCCCGCTCCCACAGACTGATTGCCATTAAGCGGGCCTATCGCTACAGTGACGGCCCGTTTTTTGTTGCCTCGCTAACGAATCATGACTCCCTTAGAACGCTATCAAGCAGATCTGAAACGTc                  >  4:4431/3‑149 (MQ=38)
                             ccGGCTTTGCCGGTGTTCGCGG‑GCACGCCCGCTCCCACAGACTGATTGCCATTAAGCGGGCCTATCGCTACAGTGACGGCCCGTTTTTTGTTGCCTCGCTAACGAATCATGACTCCCTTAGAACGCTATCAAGCAGATCTGAAACGTcc                 <  3:4431/148‑1 (MQ=39)
                                  tttGCCGGTGTTCGCGG‑GCACGCCCGCTCCCACAGACTGATTGCCATTAAGCGGGCCTATCGCTACAGTGACGGCCCGTTTTTTGTTGCCTCGCTAACGAATCATGACTCCCTTAGAACGCTATCAAGCAGATCTGAAACGTCCCGAct            >  5:197514/1‑149 (MQ=39)
                                          tgttcgcGG‑GCACGCCCGCTCCCACAGACTGATTGCCATTAAGCGGGCCTATCGCTACAGTGACGGCCCGTTTTTTGTTGCCTCGCTAACGAATCATGACTCCCTTAGAACGCTATCAAGCAGATCTGAAACGTCCCGACTTCTTCCAt    >  5:166374/7‑149 (MQ=255)
                                            ttcgcGG‑GCACGCCCGCTCCCACAGACTGATTGCCATTAAGCGGGCCTATCGCTACAGTGACGGCCCGTTTTTTGTTGCCTCGCTAACGAATCATGACTCCCTTAGAACGCTATCAAGCAGATCTGAAACGTCCCGACTTCTTCCATGa  <  6:93592/145‑1 (MQ=255)
                                            ttcgcGG‑GCACGCCCGCTCCCACAGACTGATTGCCATTAAGCGGGCCTATCGCTACAGTGACGGCCCGTTTTTTGTTGCCTCGCTAACGAATCATGACTCCCTTAGAACGCTATAAAGCAGAACTCAAACGTCCCGACATCTTACATGa  >  2:100216/5‑149 (MQ=255)
                                             tcgcGG‑GCACGCCCGCTCCCACAGACTGATTGCCATTAAGCGGGCCTATCGCTACAGTGACGGCCCGTTTTTTGTTGCCTCGCTAACGAATCATGACTCCCTTAGAACGCTATCAAGCAGATCTGAAACGTCCCGACttcttc       <  7:297020/140‑1 (MQ=255)
                                             tcgcGG‑GCACGCCCGCTCCCACAGACTGATTGCCATTAAGCGGGCCTATCGCTACAGTGACGGCCCGTTTTTTGTTGCCTCGCTAACGAATCATGACTCCCTTAGAACGCTATCAAGCAGATCTGAAACGTCCCGACttcttc       >  8:297020/4‑143 (MQ=255)
                                               ||                                                                                                                                                 
CCGCGAAACAGGCGCCGCGGTGCATGGCCACGGCTTTGCCGGTGTT‑TCGG‑GCACG‑CCGCTCCCACAGACTGATTGCCATTAAGCGGGCCTATCGCTACAGTGACGGCCCGTTTTTTGTTGCCTCGCTAACGAATCATGACTCCCTTAGAACGCTATCAAGCAGATCTGAAACGTCCCGACTTCTTCCATGA  >  NC_002947/1499452‑1499642

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