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.C83.2% 20.2 / ‑3.4 19intergenic (+74/‑90)trpS/zapEtryptophan‑‑tRNA ligase/nucleoside triphosphate hydrolase domain‑containing protein
Reads supporting (aligned to +/- strand):  ref base . (3/0);  major base C (3/10);  minor base A (3/0);  total (9/10)
Fisher's exact test for biased strand distribution p-value = 3.57e-02
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: Variant not supported by required number of reads on each strand.
*NC_0029471,499,4980TG84.2% 24.1 / 2.7 19intergenic (+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 (6/10);  minor base C (3/0);  total (9/10)
Fisher's exact test for biased strand distribution p-value = 8.67e-02
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 8.69e-01
Rejected as polymorphism: Variant not supported by required number of reads on each strand.

CCCGTGATGCGGCCATTCAGAGCCTGAACGAAGCCTTGGCACCACAGCAGGACTGAGGCTGCCAAACGGTGGGAGCGGGCGTGCCCGCGAAACAGGCGCCGCGGTGCATGGCCACGGCTTTGCCGGTGTTTCGGGCACGCCGCTCCCACAGACTGATTGCCATTAAGCGGGCCTATCGCTACAGTGACGGCCCGTTTTTTGTTGCCTCGCTAACGAATCATGACTCCCT  >  NC_002947/1499368‑1499596
                                                                                                                                   |||                                                                                                    
cccGTGATGCGGCCATTCAGAGCCTGAACGAAGCCTTGGCACCACAGCAGGACTGAGGCTGCCAAACGGTGGGAGCGGGCGTGCCCGCGCAACAGGCGCCGCGGTGCATGGCACCGGCGTTGCCGGTGTTAGCGGGCCcgcccgcgccc                                                                                    >  8:77785/1‑142 (MQ=25)
  cGTGATGCGGCCATTCAGAGCCTGAACGAAGCCTTGGCACCACAGCAGGACTGCGGCTGCCAAACGGTGGGAGCGGGCGTGCCCCCGAAACACGCGCCCCGGTGCAAGGCCCCCGCGTTGCCCGTGTTAGCGGGCACGcccccccccc                                                                                  >  3:132613/1‑146 (MQ=255)
                     gCCTGAACGAAGCCTTGGCACCACAGCAGGACTGAGGCTGCCAAACGGTGGGAGCGGGCGTGCCCGCGAAACAGGCGCCGCGGTGCATGGCACCGGCTTTGCCGGTGTTCGCGGGCACGCCCGCTCCCACAGACTGATTGCCATTAAGc                                                               <  3:49821/149‑1 (MQ=21)
                          aaCGAAGCCTTGGCACCACAGCAGGACTGAGGCTGCCAAACGGTGGGAGCGGGCGTGCCCGCGAAACAGGCGCCGCGGTGCATGGCACCGGCTTTGCCGGTGTTCGCGGGCACGCCCGCTCCCACAGACTGATTGCCATTAAGCGGGc                                                           <  7:77785/148‑1 (MQ=21)
                          aaCGAAGCCTTGCCACCACAGCAGGACTGAGGCTGCCAAACGGTGGGAGCGGGCGTGCCCGCGAAACAGGCGCCGCGGTGCATGGCACCGGCTTTGCCGGTGTTCGCGGGCACGCCCGCTCCCACAGACTGATTGCCATTAAGCggg                                                            <  3:216054/147‑1 (MQ=21)
                                gCCTTGGCACCACAGCAGGACTGAGGCTGCCAAACGGTGGGAGCGGGCGTGCCCGGGAAACAGGCGCCCCGGGGGATGGCACCGGCTTTGCCGGTGTTAGCGGGCCCGCCCGCGCCCACAGAATGCATGCCATTAAACGGGCCAAtagg                                                    >  6:262141/1‑146 (MQ=18)
                                     ggCACCACAGCAGGAGTGGGGCTGCCAAAGGGTGGGAGGGGGGGTGCCCGCTAAACTGGCGCGGGGGTGCATGGCACGGGCTTTGCGGGTGTTCGCGGGCACGCCCGCTCCCACAGACTGATTGCCATTAAGCGGGCCTATCGCTACa                                                <  4:132613/148‑1 (MQ=21)
                                     ggCACCACAGCAGGACTGTGGCGGCCAAACGGTGGGAGCGGGCGTGCCCGCGAAACAGGCGCCGCGGTGCAGGGCACCGGCTTTGCCGGTGTTCGCGGGCACGCCCGCTCCCACAGACTGATTGCCATTAAGCGGGCCTATCGCTACa                                                <  3:63644/148‑1 (MQ=21)
                                     ggCACCACAGCAGGACTGAGGCTGCCAAACGGTGGGAGCGGGCGTGCCCGCGAAACAGGCGCCGCGGTGCATGGCACCGGCTTTGCCGGTGTTCGCGGGCACCCCCGCTCCCACAGACTGATTGCCATTAAGCGGGCCTATCGCTACa                                                <  1:242070/148‑1 (MQ=21)
                                      gCACCACAGCAGGACTGAGGCTGCCAAACGGTGGGAGCGGGCGTGCCCGCGAAACAGGCGCCGCGGTGCATGGCACCGGCTTTGCCGGTGTTCGCGGGCACGCCCGCTCCCACAGACTGATTGCCATTAAGCGGGCCTATCGCTACAGt                                              <  3:112071/149‑1 (MQ=21)
                                             agcagGACTGAGGCTGCCAAACGGTGGGAGCGGGCGTGCCCGCGAAACAGGCGCCGCGGTGCATGGCACCGGCTTTGCCGGTGTTCGCGGGCCCGCCCGCTCCCACAGACTGATTGCCATTAAGCGGGCCTATCGCTACAGTGACGGcc                                       >  7:235036/1‑149 (MQ=21)
                                                       aGGCTGCCAAACGGTGGGAGCGGGCGTGCCCGCGAAACAGGCGCCGCGGTGCATGGCACCGGCTTTGCCGGTGTTCGCGGGCACGCCCGCTCCCACAGACTGATTGCCATTAAGCGGGCCTATCGCTACAGTGACGGCCCGtttttt                               <  1:228925/147‑1 (MQ=21)
                                                             ccAAACGGTGGGAGCGGGCGTGCACGCGAAACAGGCGCCGCGGTGCATGGCACCGGCTTTGCCGGTGTTCGCGGGCACGCCCGCTCCCACAGACTGATTGCCa                                                                     >  7:310694/1‑103 (MQ=14)
                                                             ccAAACGGTGGGAGCGGGCGTGCACGCGAAACAGGCGCCGCGGTGCATGGCACCGGCTTTGCCGGTGTTCGCGGGCACGCCCGCTCCCACAGACTGATTGCCa                                                                     <  8:310694/103‑1 (MQ=14)
                                                                        tagggggCGTTCCCGCGAAACAGGCGCCGCGGTGCAGGGCACGGGCTTTGCCGGTGTTCGGGGGAACGCCCGCTCCCACAGACTGATTGCCATTAAGCGGGCCTATCGCTACAGTGACGGCCCGTTTTTTGTTGCCTCGCTAACGAATc            <  4:147555/145‑1 (MQ=255)
                                                                        GAGCGGGCGTGCCCGCGAAACAGGCGCCGCGGTGCATGGCACCGGCTTTGCCGGTGTTCGCGGGCACGCCCGCTCCCACAGACTGATTGCCATTAAGCGGGCCTATCGCTACAGTGACGGCCCGTTTTTTGTTGCCTCGCTAACGAATc            >  8:255317‑M1/1‑149 (MQ=255)
                                                                          GCGGGCGTGCCCGCGAAACAGGCGCCGCGGTGCATGGCACCGGCTTTGCCGGTGTTCGCGGGCACGCCCGCTCCCACAGACTGATTGCCATTAAGCGGGCCTATCGCTACAGTGACGGCCCGTTTTTTGTTGCCTCGCTAACGAATCAt          >  7:151225‑M1/1‑149 (MQ=255)
                                                                           CGGGCGTGCCCGCGAAACAGGCGCCGCGGTGCATGGCACCGGCTTTGCCGGTGTTCGCGGGCACGCCCGCTCCCACAGACTGATTGCCATTAAGCGGGCCTATCGCTACAGTGACGGCCCGTTTTTTGTTGCCTCGCTAACGAATCATg         >  3:220638‑M1/1‑149 (MQ=255)
                                                                                    ccGCGAAACAGGCGCCGCGGTGCATGGCACCGGCTTTGCCGGTGTTCGCGGGCACGCCCGCTCCCACAGACTGATTGCCATTAAGCGGGCCTATCGCTACAGTGACGGCCCGTTTTTTGTTGCCTCGCTAACGAATCATGACTCCCt  >  1:253442/1‑147 (MQ=37)
                                                                                                                                   |||                                                                                                    
CCCGTGATGCGGCCATTCAGAGCCTGAACGAAGCCTTGGCACCACAGCAGGACTGAGGCTGCCAAACGGTGGGAGCGGGCGTGCCCGCGAAACAGGCGCCGCGGTGCATGGCCACGGCTTTGCCGGTGTTTCGGGCACGCCGCTCCCACAGACTGATTGCCATTAAGCGGGCCTATCGCTACAGTGACGGCCCGTTTTTTGTTGCCTCGCTAACGAATCATGACTCCCT  >  NC_002947/1499368‑1499596

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: