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.C95.7% 54.3 / ‑4.5 24intergenic (+74/‑90)trpS/zapEtryptophan‑‑tRNA ligase/nucleoside triphosphate hydrolase domain‑containing protein
Reads supporting (aligned to +/- strand):  ref base . (0/1);  major base C (10/12);  minor base A (1/0);  total (11/13)
Fisher's exact test for biased strand distribution p-value = 4.78e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.06e-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.
*NC_0029471,499,4980TG100.0% 68.2 / NA 24intergenic (+75/‑89)trpS/zapEtryptophan‑‑tRNA ligase/nucleoside triphosphate hydrolase domain‑containing protein
Reads supporting (aligned to +/- strand):  ref base T (0/0);  new base G (11/13);  total (11/13)
Rejected as polymorphism: Frequency below/above cutoff threshold.
Rejected as polymorphism: Variant not supported by required number of reads on each strand.

GCTGAAGCCCGTGATGCGGCCATTCAGAGCCTGAACGAAGCCTTGGCACCACAGCAGGACTGAGGCTGCCAAACGGTGGGAGCGGGCGTGCCCGCGAAACAGGCGCCGCGGTGCATGGCCACGGCTTTGCCGGTGT‑T‑TCGGGCACG‑CCGCTCCCACAGACTGATTGCCATTAAGCGGGCCTATCGCTACAGTGACGGCCCGTTTTTTGTTGCCTCGCTAACGAATCATGACTCCCTTAGAACGCTATCAAGCAGATCTGAAACGTCCCGACTTCT  >  NC_002947/1499361‑1499635
                                                                                                                                           ||                                                                                                                                         
gCTGAAGCCCGTGATGCGGCCATTCAGAGCCTGAACGAAGCCTTGGCACCACAGCAGGACTGAGGCTGCCAAACGGTGGGAGCGGGCGTGCCCCCGAAACAGGCGCCGCGGTGCATGGCACCGGCTTTGCCGGTGT‑TCGCGGGCAcg‑cc                                                                                                                                 >  8:204593/1‑149 (MQ=37)
    aaGCCCGTGATGCGGCCATTCAGAGCCTGAACGAAGCCTTGGCACCACAGCAGGACTGAGGCTGCCAAACGGTGGGAGCGGGCGTGCCCGCGAAACAGGCGCCGCGGTGCATGGCACCGGCTTTGCCGGTGT‑TCGCGGGCAcg‑cccgct                                                                                                                             <  6:123081/149‑5 (MQ=25)
       cccGTGATGCGGCCATTCAGAGCCTGAACGAAGCCTTGGCACCACAGCAGGACTGAGGCTGCCAAACGGTGGGAGCGGGCGTGCCCGCGAAACAGGCGCCGCGGTGCATGGCACCGGCTTTGCCGGTGT‑TCGCGGGCACGCCCGCTccc                                                                                                                           <  8:99398/149‑1 (MQ=25)
                              cTGAACGAAGCCTTGGCACCACAGCAGGACTGAGGCTGCCAAACGGTGGGAGCGGGCGTGCCCGCGAAACAGGCGCCGCGGTGCAGGGCACCGGCTTTGCCGGTGT‑TCGCGGGCACGCCCGCTCCCACAGACTGATTGCCATTAAGCgg                                                                                                    <  3:285497/149‑1 (MQ=21)
                                 aaCGAAGCCTGGCCACCACCCCTGCACTGCGGCTGCAAAACGGGGGGTGGGGGCGTGCCCGCTAAACCGGCGCCGCGGTGAAGGCCCCGGGCTTTGCGGGTGT‑TCGCGGGAACGCCCGCTCCCACAGACTGATTGCCATTAAGCGGGcc                                                                                                 <  4:66640/149‑1 (MQ=255)
                                  aCGAAGCCTTGGCACCACAGCAGGACTGAGGCTGCCAAACGGTGGGAGCGGGCGTGCCCGCGCAACAGGCGCCGCGGTGGATGGCACCGGGTTTGCCGGTGT‑TAGGGGGCACGCCCGCGGCCCCAGACTGCTTGGCATTAAGCggggc                                                                                                 >  3:66640/1‑146 (MQ=255)
                                                   cagcagGCCTCAGGCTGCCAAACGGGGGGAGCGGGCGTGCCCGCGAAACCGGCGCCGCGGTGCATGGCACCGGCTTTGCCGGTGT‑TCGCGGGCACGCCCGCTCCCACAGACTGATTGCCATTAAGCGGGCCTATCGCTACAGTGACGGc                                                                               <  8:874/149‑1 (MQ=21)
                                                   cagcagGACTGAGGCTGCCAAACGGTGGGAGCGGGCGTGCCCGCGAAACAGGCGCCGCGGTGCATGGCACCGGCTTTGCCGGTGT‑TCGCGGGCACGCCCGCTCCCACAGACTGATTGCCATTAAGCGGGCCTATCGCTACAGTGACGGc                                                                               <  7:204593/149‑1 (MQ=21)
                                                              aGGCTGCCAAACGGTGGGAGCGGGCGTGCCCGCGAAACAGGCGCCGCGGTGCATGGCACCGGGTTTGCCGGTGT‑TCGCGGGCCCGCCCGCGCCCACAGACGGCTTGcc                                                                                                             >  3:42407/1‑108 (MQ=11)
                                                                                   gggcgggcCCGCTAAACAGGCGCCGGGGTGCATGGCACCGGCTTTGCCGGTGTGC‑GCGGGAAGGCCCGCTCCCACAGACTGATTGCCATTAAGCGGGCCTATCGCTACAGTGACGGCCCGTTTTTTGTTGCCTCGCTAACGAATCATGa                                               <  1:112602/149‑1 (MQ=25)
                                                                                       gTGCCCGCGAAACAGGCGCCGCGGTGCATGGCACCGGCTTTGCCGGTGT‑TCGCGGGCACGCCCGCTCCCACAGACTGATTGCCATTAAGCGGGCCTATCGCTACAGTGACGGCCCGTTTTTTGTTGCCTCGCTAACGAATCATGACTcc                                           >  2:110265/1‑149 (MQ=25)
                                                                                          cccGCGAAACAGGCGCCGCGGTGCATGGCACCGGCTTTGCCGGTGT‑TCGCGGGCACGCCCGCTCCCACAGACTGATTGCCATTAAGCGGGCCTATCGCTACAGTGACGGCCCGTTTTTTGTTGCCTCGCTAACGAATCATGACTCCCtt                                        <  1:110265/149‑1 (MQ=37)
                                                                                           ccGCGAAACAGGCGCCGCGGTGCATGGCACCGGCTTTGCCGGTGT‑TCGCGGGCACGCCCGCTCCCACAGACTGATTGCCATTAAGCGGGCCTATCGCTACAGTGACGGCCCGTTTTttg                                                                     <  8:70661/119‑1 (MQ=21)
                                                                                           ccGCGAAACAGGCGCCGCGGTGCATGGCACCGGCTTTGCCGGTGT‑TCGCGGGCACGCCCGCTCCCACAGACTGATTGCCATTAAGCGGGCCTATCGCTACAGTGACGGCCCGTTTTttg                                                                     >  7:70661/1‑119 (MQ=21)
                                                                                            cgcgAAACAGGCGCCGCGGTGCATGGCACCGGCTTTGCCGGTGT‑TCGCGGGCACGCCCGCTCCCACAGACTGATTGCCATTAAGCGGGCCTATCGCTACAGTGACGGCCCGTTTTTTGTTGCCTCGCTAACGAATCAt                                                 >  6:189471/1‑138 (MQ=25)
                                                                                                 aaCAGGCGCCGCGGTGCATGGCACCGGCTTTGCCGGTGT‑TCGCGGGCACGCCCGCTCCCACAGACTGATTGCCATTAAGCGGGCCTATCGCTACAGTGACGGCCCGTTTTttgtt                                                                   <  7:116858/115‑1 (MQ=25)
                                                                                                 aaCAGGCGCCGCGGTGCATGGCACCGGCTTTGCCGGTGT‑TCGCGGGCACGCCCGCTCCCACAGACTGATTGCCATTAAGCGGGCCTATCGCTACAGTGACGGCCCGTTTTttgtt                                                                   >  8:116858/1‑115 (MQ=25)
                                                                                                             ggTGCATTGCTCCTGCTTTGCCGGTTT‑TCGCGGGCACGCCCGCTCCCACAGACTGATTGCCATTAAGCGGTCCTATCTCTACAGTGACTGCCCGTTTTTTGTTGCCTCGCTAACTAATCATGACTCCCTTAGAACGCTATCAAGCAGAt                     <  3:107684/149‑1 (MQ=37)
                                                                                                                   tggcaccGGCTTTGCCGGTGT‑TCGCGGGCACGCCCGCTCCCACAGACTGATTGCCATTAAGCGGGCCTATCGCTACAGTGACGGCCCGTTTTTTGTTGCCTCGCTAACGAATCATGACTCCCTTAGAACGCTATCAAGCAGATCTGaaa               <  4:208596/143‑1 (MQ=38)
                                                                                                                     gcaccGGCTTTGCCGGTGT‑TCGCGGGCACGCCCGCTCCCACAGACTGATTGCCATTAAGCGGGCCTATCGCTACAGTGACGGCCCGGTTTTTGTTGCCTCGCTAACGAATCATGACTCCCTTAGAACGCTATCAAGCAGATCTGAAACg             >  5:74717/5‑149 (MQ=38)
                                                                                                                       accGGCTTTGCCGGTGT‑TCGCGGGCACGCCCGCTCCCACAGACTGATTGCCATTAAGCGGGCCTATCGCTACAGTGACGGCCCGTTTTTTGTTGCCTCGCTAACGAATCATGACTCCCTTAGAACGCTATCAAGCAGATCTGAAACGTc           >  3:113249/3‑149 (MQ=38)
                                                                                                                           gCTTTGCCGGTTT‑TCGCGGGCACGCCCGCTCCCACAGACTGATTGCCATTAAGCGGGCCTATCGCTACAGTGACGGCCCGTTTTTTGTTGCCTCGCTAACGAATCATGACTCCCTTAGAACGCTATCAAGCAGATCTGAAACGTCCCGa       <  3:135307/149‑1 (MQ=255)
                                                                                                                                gCCGGTGT‑TCGCGGGCACGCCCGCTCCCACAGACTGATTGCCATTAAGCGGGCCTATCGCTACAGTGACGGCCCGTTTTTTGTTGCCTCGCTAACGAATCATGACTCCCTTAGAACGCTATCAAGCAGATCTGAAACGTCCCGActtct  >  5:234763/1‑149 (MQ=255)
                                                                                                                                 ccGGTGT‑TCGCGGGCACGCCCGCTCCCACAGACTGATTGCCATTAAGCGGGCCTATCGCTACAGTGACGGCCCGTTTTTTGTTGCCTCGCTAACGAATCATGACTCCCTTAGAACGCTATCAAGCAGATCTGAAACGTCCCGActtct  >  4:72748/1‑148 (MQ=255)
                                                                                                                                          cgcGGGCACGCCCGCTCCCACAGACTGATTGCCATTAAGCGGGCCTATCGCTACAGTGACGGCCCGTTTTTTGTTGCCTCGCTAACGAATCATGACTCCCTTAGAACGCTATCAAGCAGATCTGAAACGTcc          >  2:5788/3‑132 (MQ=255)
                                                                                                                                          cgcGGGCACGCCCGCTCCCACAGACTGATTGCCATTAAGCGGGCCTATCGCTACAGTGACGGCCCGTTTTTTGTTGCCTCGCTAACGAATCATGACTCCCTTAGAACGCTATCAAGCAGATCTGAAACGTcc          <  1:5788/130‑1 (MQ=255)
                                                                                                                                           gcGGGCACGCCCGCTCCCACAGACTGATTGCCATTAAGCGGGCCTATCGCTACAGTGACGGCCCGTTTTTTGTTGCCTCGCTAACGAATCATGACTCCCTTa                                       <  6:163412/101‑1 (MQ=255)
                                                                                                                                           gcGGGCACGCCCGCTCCCACAGACTGATTGCCATTAAGCGGGCCTATCGCTACAGTGACGGCCCGTTTTTTGTTGCCTCGCTAACGAATCATGACTCCCTTa                                       >  5:163412/2‑102 (MQ=255)
                                                                                                                                           ||                                                                                                                                         
GCTGAAGCCCGTGATGCGGCCATTCAGAGCCTGAACGAAGCCTTGGCACCACAGCAGGACTGAGGCTGCCAAACGGTGGGAGCGGGCGTGCCCGCGAAACAGGCGCCGCGGTGCATGGCCACGGCTTTGCCGGTGT‑T‑TCGGGCACG‑CCGCTCCCACAGACTGATTGCCATTAAGCGGGCCTATCGCTACAGTGACGGCCCGTTTTTTGTTGCCTCGCTAACGAATCATGACTCCCTTAGAACGCTATCAAGCAGATCTGAAACGTCCCGACTTCT  >  NC_002947/1499361‑1499635

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