Predicted mutation
evidence seq id position mutation annotation gene description
RA NC_002947 1,499,498 1 bp→CG 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.C97.4% 109.3 / ‑3.6 38intergenic (+74/‑90)trpS/zapEtryptophan‑‑tRNA ligase/nucleoside triphosphate hydrolase domain‑containing protein
Reads supporting (aligned to +/- strand):  ref base . (0/1);  new base C (14/23);  total (14/24)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00
*NC_0029471,499,4980TG100.0% 112.8 / NA 38intergenic (+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 (14/24);  total (14/24)

CGCTGAAGCCCGTGATGCGGCCATTCAGAGCCTGAACGAAGCCTTGGC‑‑‑ACCACAGCAGGACTGAGGCTGCCAAACGGTGGGAGCGGGCGTGCCCGCGAAACAGGCGCCGCGGTGCATGGCCACGGCTTTGCCGGTGTTTCGGG‑‑‑CACGCCGCTCCCACAGACTGATTGCCATTAAGCGGGCCTATCGCTACAGTGACGGCCCGTTTTTTGTTGCCTCGCTAACGAATCATGACTCCCTTAGAACGCTATCAAGCAGATCTGAAACGTCCCGACTTCTTCCATG  >  NC_002947/1499360‑1499641
                                                                                                                                              ||                                                                                                                                                  
cGCTGAAGCCCGTGATGCGGCCATTCAGAGCCTGAACGAAGCCTTGGC‑‑‑ACCACAGCAGGACTGAGGCTGCCAAACGGTGGGAGCGGGCGTGCCCGCGAAACAGGCGCCGCGGTGCATGGCACCGGCTTTGCCGGTGTTCGCGGG‑‑‑CAcgc                                                                                                                                        <  2:37615/149‑1 (MQ=37)
 gCTGAAGCCCGTGATGCGGCCATTCAGAGCCTGAACGAAGCCTTGGCCTTAACTCAGCAGGACTGAGGCTGCCAAACGGTGGGAGCGGGCGTGCCCGCGAAACAGGCGCCGCGGTGCATGGCACCGGCTTTGCCGGTGTTCGCGGG‑‑‑CAc                                                                                                                                           <  1:137929/149‑1 (MQ=21)
  cTGAAGCCCGTGATGCGGCCATTCAGAGCCTGAACGAAGCCTTGGC‑‑‑ACCACAGCAGGACTGAGGCTGCCAAACGGTGGGAGCGGGCGTGCCCGCGAAACAGGCGCCGCGGTGCATGGCACCGGCTTTGCCGGTGTTCGCGGG‑‑‑CAcgccc                                                                                                                                      <  3:143505/149‑2 (MQ=37)
       gCCCGTGATGCGGCCATTCAGAGCCTGAACGAAGCCTTGGC‑‑‑ACCACAGCAGGACTGAGGCTGCCAAACGGTGGGAGCGGGCGTGCCCGCGAAACAGGCGCCGCGGTGCATGGCACCGGCTTTGCCGGTGTTCGCGGG‑‑‑CACGCCCGCTcc                                                                                                                                  >  7:342386/1‑149 (MQ=25)
                 cGGCCATTCAGAGCCTGAACGAAGCCTTGGC‑‑‑ACCACAGCAGGACTGAGGCTGCCAAACGGTGGGAGCGGGCGTGCCCGCGAAACAGGCGCCGCGGTGCATGGCACCGGCTTTGCCGGTGTTCGCGGG‑‑‑CAcgcccg                                                                                                                                     >  5:399210/1‑133 (MQ=25)
                 cGGCCATTCAGAGCCTGAACGAAGCCTTGGC‑‑‑ACCACAGCAGGACTGAGGCTGCCAAACGGTGGGAGCGGGCGTGCCCGCGAAACAGGCGCCGCGGTGCATGGCACCGGCTTTGCCGGTGTTCGCGGG‑‑‑CAcgcccg                                                                                                                                     <  6:399210/135‑3 (MQ=25)
                  ggCCATTCAGAGCCTGAACGAAGCCTTGGC‑‑‑ACCACAGCAGGACTGAGGCTGCCAAACGGTGGGAGCGGGCGTGCCCGCGAAACAGGCGCCGCGGTGCATGGCACCGGCTTTGCCGGTGTTCGCGGG‑‑‑CAcgcccgctc                                                                                                                                  <  7:376523/137‑6 (MQ=21)
                   gCCATTCAGAGCCTGAACGAAGCCTTGGC‑‑‑ACCACAGCAGGACTGAGGCTGCCAAACGGTGGGAGCGGGCGTGCCCGCGAAACAGGCGCCGCGGTGCATGGCACCGGCTTTGCCGGTGTTCGCGGG‑‑‑CACGCCCGCTCCCACAGACTGAtt                                                                                                                      <  5:314660/149‑1 (MQ=21)
                     cATTCATAGCCTGAAGGAAGCCTTGGC‑‑‑ACCACAGCTGGACTGAGGCTGCCAAACGGTGGGAGCGGGCGTGCCCGCGAAACAGGCGCCGCGGTGCATGGCACCGGCTTTGCCGGTGTTCGCGGG‑‑‑CAGGCCCGCTCCCACAGACTGATTGc                                                                                                                    <  6:115743/149‑1 (MQ=18)
                          agagCCTGAACGAAGCCTTGGC‑‑‑ACCACAGCAGGACTGAGGCTGCCAAACGGTGGGAGCGGGCGTGCCCGCGAAACAGGCGCAGCGGTGCATGGCACCGGCTTTGCCGGTGTTCGCGGG‑‑‑CAcgcccgct                                                                                                                                   <  1:311837/128‑5 (MQ=21)
                               cTGAACGAAGCCTTGGC‑‑‑ACCACAGCAGGACTGAGGCTGCCAAACGGTGGGAGCGGGCGTGCCCGCGAAACAGGCGCCGCGGTGCATGGCACCGGCTTTGCCGGTGTTCGCGGG‑‑‑CACGCCCGCTCCCACAGACTGATTGCCATTAAGCgg                                                                                                          >  7:321540/1‑149 (MQ=21)
                                  aaCGAAGCCTTGGC‑‑‑ACCACACCAGGACTGAGGCTGCCAAACGGTGGGAGCGGGCGTGCCCGCGAAACAGGCGCCGCGGTGCATGGCACCGGCTTTGCCGGTGTTCGCGGG‑‑‑CACGCCCGCTCCCACAGACTGATTGCCATTAAGCGTGc                                                                                                        <  2:63656/148‑1 (MQ=21)
                                              gC‑‑‑ACCACAGCAGGACTGAGGCTGCCAACCGGTGGGAGCGGGCGTGCCCGCTAAACAGGCGCCGCGGTGCATGGCACCGGCTTTGCCGGTGTTCGCGGG‑‑‑CACGCCCGCTCCCACAGACTGATTGCCATTAAGCGGGCCTATCGCTACAGt                                                                                           <  3:133593/149‑1 (MQ=21)
                                                    ccacAGCAGGACTGAGGCTGCAAAAGGGGGGCTGCGGGCGTGCCCGCTAAACTGGCGCCGCGGTGCAGGGCACGGGCTTTGCCGGTGTTCGCGGG‑‑‑CACGCCCGCTCCCACAGACTGATTGCCATTAAGCGGGCCTATCGCTACAGTGAc                                                                                        <  4:414339/149‑1 (MQ=21)
                                                      acaGCAGGACTGAGGCTGCCAAACGGTGGGAGCGGGCGTGCCCGCGAAACAGGCGCCGCGGTGCATGGCACCGGCTTTGCCGGTGTTCGCGGG‑‑‑CACGCCCGCTCCCACAGACTGATTGCCATTAAGCGGGCCTATCGCTACAGTGACgg                                                                                      >  5:427971/1‑149 (MQ=21)
                                                                  aGGCGGCCAAACGGGGGGAGCGGGCGGCCCCGCGAAACAGGCGCCGCGGGGCATGGCACCGGCTTTGCCGGTGTTCGCGGG‑‑‑CACGCCCGCTCCCACAGACTGATTGCCATTAAGCGGGCCTATCGCTACAGTGACGGCCCGtttttt                                                                            <  3:179131/147‑1 (MQ=21)
                                                                    gCTGCCAAACGGTGGGAGCGGGCGTGCCCGCGAAACAGGCGCCGCGGTGCATGGCACCGGCTTTGCCGGTGTTCGCGGG‑‑‑CACGCCCGCTCCCACAGACTGATTGCCATTAAGCGGGCCTATCGCTACAGTGACGGCCCGTTTTttgtt                                                                         <  2:256361/148‑1 (MQ=21)
                                                                     cTGCCAAACGGTGGGAGCGGGCGTGCCCGCGAAACAGGCGCCGCGGTGCATGGCACCGGCTTTGCCGGTGTTCGCGGG‑‑‑CACGCCCGCTCCCACAGACTGATTGCCATTAAGCGGGCCTATCGCTACAGTGACGGCCCGTTTTTTGTTGc                                                                       <  8:342386/149‑1 (MQ=21)
                                                                                 gggAGCGGGCGTGCCCGCGAAACAGGCGCCGCGGTGCATGGCACCGGCTTTGCCGGTGTTCGCGGG‑‑‑CACGCCCGCTCCCACAGACTGATTGCCATTAAGCGGGCCTATCGCTACAGTGACGGCCCGTTTTTTGTTGCCTCGCTAACGaa                                                           <  3:310004/149‑1 (MQ=25)
                                                                                             gCCCGCCAAACCGGCGGCGGGGTGCCCGGCCCGGGGTTTGGGGGTGTGGCGGGGCCCCCCCCGCTCCCACAGACTGATTGCCATTAAGCGGGCCTATCGCTACAGTGACGGCCCGTTTTTTGTTGCCTCGCTAACGAATCATGAcccct                                                <  2:194897/149‑2 (MQ=255)
                                                                                                      aCAGGCGCCGCGGTGCATGGCACCGGCTTTGCCGGTGTTCGCGGG‑‑‑CACGCCCGCTCCCACAGACTGATTGCCATTAAGCGGGCCTATCGCTACAGTGACGGCCCGTTTTTTGTTGCCTCGCTAACGAATCATGACTCCCTTAGAACGCt                                      <  6:427971/149‑1 (MQ=37)
                                                                                                      aCAGGCGCCGCGGTGCATGGCACCGGCTTTGCCGGTGTTCGCGGG‑‑‑CACGCCCGCTCCCACAGACTGATTGCCATTAAGCGGGCCTATCGCTACAGTGACGGCCCGTTTTTTGTTGCCTCGCTAACGAATCATGACTCCCTTAGAACGCt                                      <  8:321540/149‑1 (MQ=37)
                                                                                                       cAGGCGCCGCGGTGCATGGCACCGGCTTTGCCGGTGTTCGCGGG‑‑‑CACGCCCGCTCCCACAGACTGATTGCCATTAAGCGGGCCTATCGCTACAGTGACGGCCCGTTTTTTGTTGCCTCGCTAACGAATCATGACTCCCTTAGAACGCTa                                     >  3:346028/1‑149 (MQ=37)
                                                                                                          gcgcCGCGGTGCATGGCACCGGCTTTGCCGGTGTTCGCGGG‑‑‑CACGCCCGCTCCCACAGACTGATTGCCATTAAGCGGGCCTATCGCTACa                                                                                             >  6:292970/1‑90 (MQ=17)
                                                                                                          gcgcCGCGGTGCATGGCACCGGCTTTGCCGGTGTTCGCGGG‑‑‑CACGCCCGCTCCCACAGACTGATTGCCATTAAGCGGGCCTATCGCTACa                                                                                             <  5:292970/90‑1 (MQ=17)
                                                                                                            gccgcGGTGCATGGCACGGGCTTTGCCGGTGTTCGGGGG‑‑‑AACGCCCGCTCCCACAGACTGATTGCCATTAAGCGGGCCTATCGCTACAGTGACGGCCCGTTTTTTGTTGCCTCGCTAACGAATCATGACTCCCTTAGAACGCTATCAAg                                <  1:399416/149‑1 (MQ=255)
                                                                                                             ccgcGGTGCATGGCACCGGCTTTGCCGGTGTTCGCGGG‑‑‑CACGCCCGCTCCCACAGACTGATTGCCATTAAGCGGGCCTATCGCTACAGTGACGGCCCGTTTTTTGTTGCCTCGCTAACGAATCATGACTCCCTTAGAACGCTATCAAg                                >  5:392424/1‑148 (MQ=37)
                                                                                                                         gcaccGGCTTTGCCGGTGTTCGCGGG‑‑‑CACGCCCGCTCCCACAGACTGATTGCCATTAAGCGGGCCTATCGCTACAGTGACGGCCCGTTTTTTGTTGCCTCGCTAACGAATCATGACTCCCTTAGAACGCTATCAAGCAGATCTGAAACg                   >  8:382001/5‑149 (MQ=38)
                                                                                                                           accGGCTTTGCCGGTGTTCGCGGG‑‑‑CACGCCCGCTCCCACAGACTGATTGCCATTAAGCGGGCCTATCGCTACAGTGACGGCCCGTTTTTTGTTGCCTCGCTAACGAATCATGACTCCCTTAGAACg                                        >  2:339501/3‑126 (MQ=38)
                                                                                                                           accGGCTTTGCCGGTGTTCGCGGG‑‑‑CACGCCCGCTCCCACAGACTGATTGCCATTAAGCGGGCCTATCGCTACAGTGACGGCCCGTTTTTTGTTGCCTCGCTAACGAATCATGACTCCCTTAGAACg                                        <  1:339501/124‑1 (MQ=38)
                                                                                                                           accGGCTTTGCCGGTGTTCGCGGG‑‑‑CACGCCCGCTCCCACAGACTGATTGCCATTAAGCGGGCCTATCGCTACAGTGACGGCCCGTTTTTTGTTGCCTCGCTAACGAATCATGACTCCCTTAGAACGCTATCAAGCAGATCTGAAACGTc                 <  4:346028/147‑1 (MQ=38)
                                                                                                                           accGGCTTTGCCGGTGTTCGCGGG‑‑‑CACGCCCGCTCCCACAGACTGATTGCCATTAAGCGGGCCTATCGCTACAGTGACGGCCCGTTTTTTGTTGCCTCGCTAACGAATCATGACTCCCTTAGAACGCTATCAAGCAGATCTGAAACGTc                 >  1:207440/3‑149 (MQ=38)
                                                                                                                            ccGGCTTTGCCGGTGTTCGCGGG‑‑‑CACGCCCGCTCCCACAGACTGATTGCCATTAAGCGGGCCTATCGCTACAGTGACGGCCCGTTTTTTGTTGCCTCGCTAACGAATCATGACTCCCTTAGAACGCTATCAAGCAGATCTGAAACGTcc                >  3:82604/2‑149 (MQ=39)
                                                                                                                             cGGCTTTGCCGGTGTTCGCGGG‑‑‑CACGCCCGCTCCCACAGACTGATTGCCATTAAGCGGGCCTATCGCTACAGTGACGGCCCGTTTTTTGTTGCCTCGCTAACGAATCATGACTCCCTTAGAACGCTATCAAGCAGATCTGAAACGTccc               <  6:392424/149‑1 (MQ=39)
                                                                                                                                 tttGCCGGTGTTCGCGGG‑‑‑CACGCCCGCTCCCACAGACTGATTGCCATTAAGCGGGCCTATCGCTACAGTGACGGCCCGTTTTTTGTTGCCTCGc                                                                  >  5:429184/1‑94 (MQ=37)
                                                                                                                                 tttGCCGGTGTTCGCGGG‑‑‑CACGCCCGCTCCCACAGACTGATTGCCATTAAGCGGGCCTATCGCTACAGTGACGGCCCGTTTTTTGTTGCCTCGc                                                                  <  6:429184/94‑1 (MQ=37)
                                                                                                                                  ttGCCGGTGTTCGCGGG‑‑‑CACGCCCGCTCCCACAGACTGATTGCCATTAAGCGGGCCTATCGCTACAGTGACGGCCCGTTTTTTGTTGCCTCGCTAACGAATCATGACTCCCTTAGAACGCTATCAAGCAGATCTGAAACGTCCCGActt          >  1:212845/1‑149 (MQ=255)
                                                                                                                                     ccGGTGTTCGCGGG‑‑‑CACGCCCGCTCCCACAGACTGATTGCCATTAAGCGGGCCTATCGCTACAGTGACGGCCCGTTTTTTGTTGCCTCGCTAACGAATCATGACTCCCTTAGAACGCTATCAAGCAGATCTGAAACGTCCCGActtctt       >  2:67010/1‑149 (MQ=255)
                                                                                                                                         ttttcgcGGG‑‑‑CACGCCCGCTCCCACAGACTGATTGCCATTAAGCGGGCCTATCGCTACAGTGACGGCCCGTTTTTTGTTGCCTCGCTAACGAATCATGACTCCCTTAGAACGCTATCAAGCAGATCTGAAACGTCCCGACTTCTTCCAt   <  1:369690/143‑1 (MQ=255)
                                                                                                                                         tgttcgcGGG‑‑‑CACGCCCGCTCCCACAGACTGATTGCCATTAAGCGGGCCTATCGCTACAGTGACGGCCCGTTTTTTGTTGCCTCGCTAACGAATCATGACTCCCTTAGAACGCTATCAAGCAGATCTGAAACGTCCCGACTTCTTCCAt   <  7:159640/143‑1 (MQ=255)
                                                                                                                                           ttcgcGGG‑‑‑CACGCCCGCTCCCACAGACTGATTGCCATTAAGCGGGCCTATCGCTACAGTGACGGCCCGTTTTTTGTTGCCTCGCTAACGAATCATGACTCCCTTAGAACGCTATCAAGCAGATCTGAAACGTCCCGACTTCTTCCATg  <  4:82604/144‑1 (MQ=255)
                                                                                                                                              ||                                                                                                                                                  
CGCTGAAGCCCGTGATGCGGCCATTCAGAGCCTGAACGAAGCCTTGGC‑‑‑ACCACAGCAGGACTGAGGCTGCCAAACGGTGGGAGCGGGCGTGCCCGCGAAACAGGCGCCGCGGTGCATGGCCACGGCTTTGCCGGTGTTTCGGG‑‑‑CACGCCGCTCCCACAGACTGATTGCCATTAAGCGGGCCTATCGCTACAGTGACGGCCCGTTTTTTGTTGCCTCGCTAACGAATCATGACTCCCTTAGAACGCTATCAAGCAGATCTGAAACGTCCCGACTTCTTCCATG  >  NC_002947/1499360‑1499641

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: