Predicted mutation
evidence seq id position mutation annotation gene description
RA NC_002947 1,499,480 2 bp→AC intergenic (+57/‑106) 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,4800CA84.0% 51.7 / 2.3 25intergenic (+57/‑107)trpS/zapEtryptophan‑‑tRNA ligase/nucleoside triphosphate hydrolase domain‑containing protein
Reads supporting (aligned to +/- strand):  ref base C (1/3);  new base A (8/13);  total (9/16)
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,4810AC92.0% 66.4 / ‑3.1 25intergenic (+58/‑106)trpS/zapEtryptophan‑‑tRNA ligase/nucleoside triphosphate hydrolase domain‑containing protein
Reads supporting (aligned to +/- strand):  ref base A (0/2);  new base C (9/14);  total (9/16)
Fisher's exact test for biased strand distribution p-value = 5.20e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

GACGGGCCGCAGTTCGACAACGCTGAAGCCCGTGATGCGGCCATTCAGAGCCTGAACGAAGCCTTGGCACCACAGCAGGACTGAGGCTGCCAAACGGTGGGAGCGGGCGTGCCCGCGAAACAGGCGCCGCGGTGCATGGCCA‑C‑GGCTTTGCCGGTGT‑T‑TCGGG‑CACG‑‑CCGCTCCCACAGACTGATTGCCATTAAGCGGGCCTATCGCTACAGTGACGG‑CCCGTTTTTTGTTGCCTCGCTAACGAATCATGACTCCCTTAGAACG  >  NC_002947/1499340‑1499603
                                                                                                                                            ||                                                                                                                                  
gACGGGCCGCAGTTCGACAACGCTGAAGCCCGTGATGCGGCCATTCAGAGCCTGAACGAAGCCTTGGCACCACAGCAGGACTGAGGCTGCCAAACGGTGGGAGCGGGCGTGCCCGCGAAACAGGCGCCGCGGTGCATGGCAC‑C‑GGCtt                                                                                                                            <  8:144655/148‑1 (MQ=255)
     gCCGCAGTTCGACAACGCTGAAGCCCGTGATGCGGCCATTCAGAGCCTGAACGAAGCCTTGGCACCACAGCAGGACTGAGGCTGCCAAACGGTGGGAGCGGGCGTGCCCGCGAAACAGGCGCCGCGGTGCATGGCAC‑C‑GGCTTTGCCg                                                                                                                       <  8:107476/148‑1 (MQ=255)
     gCCGCAGTTCGACAACGCTGAAGCCCGTGATGCGACCATTCAGAGCCTGAACGAAGCCTTGGCACCACAGCAGGACTGAGGCTGCCAAACGGTGGGAGCGGGCGTGCCCGCGAAACAGGCGCCGCGGTGCATGGCAC‑C‑GGCTTTGCCg                                                                                                                       <  4:169204/148‑1 (MQ=255)
               gACAACGCTGAAGCCCGTGATGCGGCCATTCAGAGCCTGAACGAAGCCTTGGCACCACAGCAGGACTGAGGCGGCCAAACGGGGGGAGCGGGCGTGCCCGCGAAACAGGCGCCGCGGTGCATGGCAC‑C‑GGCTTTGCCGGTGttcgcgg                                                                                                             >  1:94325‑M2/1‑141 (MQ=255)
               gACAACGCGGAAGCCCGTGATGCGGCCATTCAGAGCCTGAACGAAGCCTTGGCACCACAGCAGGACTGAGGCTGCCAAACGGTGGGAGCGGGCGTGCCCGCGAAACAGGCGCCGCGGTGCATGGCAC‑C‑GGCTTTGCCGGTGttcgcg                                                                                                              >  2:49439‑M2/1‑141 (MQ=255)
                    cGCTGAAGCCCGTGATGCGGCCATTCAGAGCCTGAACGAAGCCTTGGCACCACAGCAGGACTGAGGCTGCCAAACGGTGGGAGCGGGCGTGCCCGCGAAACAGGCGCCGCGGTGCATGGCAC‑C‑GGCTTTGCCGGTGttcgcgggcacg                                                                                                        >  1:119149‑M2/1‑136 (MQ=255)
                     gCTGAAGCCCGTGATGCGGCCATTCAGAGCCTGAACGAAGCCTTGGCACCACAGCAGGACTGAGGCTGCCAAACGGTGGGAGCGGGCGTGCCCGCGCAACAGGCGCCGCGGGGCCTGGCAC‑C‑CGCGTTGCCGGGGttagcgggcccgc                                                                                                       <  7:19764‑M2/148‑14 (MQ=255)
                       tGAAGCCCGTGATGCGGCCATTCAGAGCCTGAACGAAGCCTTGGCACCACAGCAGGACTGAGGCTGCCAAACGGTGGGAGCGGGCGTGCCCGCGCAACAGGCGCCGCGGTGCATGGCCC‑C‑GGCGTTGCCGGTGttcgcgggcccgccc                                                                                                     <  7:57613‑M2/148‑16 (MQ=255)
                              cGTGATGCGGCCATTCAGAGCCTGAACGAAGCCTTGGCACCACAGCAGGACTGAGGCTGCCAAACGGTGGGAGCGGGCGTGCCCGCGAAACAGGCGCCGCGGTGCATGGCAC‑C‑GGCTTTGCCGGTGttcgcgggcacgcccgctccc                                                                                               <  4:111746‑M2/147‑22 (MQ=255)
                                                   cTGAACTACGCCTTGGCACAACAGCAGGAGTTCGGCTGCCAAAGGGTGGGTGCGGGCTTGCCCGCGAAACTGGCGCCGGGGTGCAGGGCAC‑C‑GGCTTGGCCGGTGT‑T‑CGGGG‑GAAGGCCCGCTCCCACAGACTGATTGCCATTAAGCg                                                                      <  8:19764/148‑1 (MQ=255)
                                                      aaCGAAGCCTTGGCACCACAGCAGGACTGAGGCTGCCAAACGGTGGGAGCGGGCGTGGCCCCGCAACCCGCGCCCCCGTGGAAGGCCC‑C‑CGCGTTGGCGGTGT‑TAGCGGG‑CACGC‑CCGCGACCCCAGCCGGATTGCCAATAAGCGGGc                                                                   >  2:38057/1‑148 (MQ=255)
                                                       aCGAAGCCTTGGCACCACAGCTGGACTGAGGCTGCCAAACGGTGGGAGCGGGCGTGCCCGCGAAACAGGCGCCGCGGTGCATGGCAC‑C‑GGCTTTGCCGGTGT‑TCGCGGG‑CACGC‑CCGCTCCCACAGACTGATTGCCATTAAGCGGGc                                                                   <  2:26782/147‑1 (MQ=21)
                                                                           cagGAGTGAGGCGGCAAAAGGGGGGGCGCGGGCGGGCCCGCTAAACCGGCGCCGGGGGGCA‑GGCCA‑CGGGGTTTGCGGGTGT‑TCGGGGG‑CACGC‑CCGCTCCCACAGACTGATTGCCATTAAGCGGGCCTATCGCTACAGTGACGG‑ccc                                             <  2:28455/148‑1 (MQ=255)
                                                                                   aGGCTGCCAAACGGTGGGAGCGGGCGTGCCCCCGCAACCGGCGCCGCGGTGCATGGCAC‑C‑GGCTTTTCCGGTGT‑TAGCGGG‑CCCGC‑CCGCTCCCCCAGCCGGAGTGGCATTAAGCGGGCCTATAGCTACAGTGACGG‑CCCGgtttttt                                     >  8:60150/1‑147 (MQ=16)
                                                                                    ggCTGCCAAACGGTGGGAGCGGGCGTGGCCCCGAAACAGGCGCCGCGGTGCCTTGCAC‑C‑CGCTTTTCCGGTGT‑TAGCGGGCCCCC‑‑CCGCGCCCACAGACAGATTTCCATTAAACCGGCCGATAGCGACACTGACGGCCCCGttttttt                                     >  4:52705/1‑147 (MQ=255)
                                                                                               ggTGGGAGCGGGCGTGCCCGCTAAACAGGCGCCGCGGTGCATGGCAC‑C‑GGCTTTGCCGGTGT‑TCGCGGG‑CACGC‑CCGCTCCCACAGACTGATTGCCATTAAGCGGGCCTATCGCTACAGTGACGG‑CCCGTTTTTTGTTGCCTCGCTaa                         <  4:166966/148‑1 (MQ=21)
                                                                                                      gCGGGGGTGCCCGCGAAACAGGCGCCGCGGTG‑ATGGGCACC‑GGCTTGGCGGGTTT‑TCGCGGG‑CACGC‑CCGCTCCCACAGACTGATTGCCATTAAGCGGGCCTATCGCTACAGTGACGG‑CCCGTTTTTTGTTGCCTCGCTAACGAATCa                  <  6:17012/148‑1 (MQ=25)
                                                                                                         ggCGTGCCCGCGAAACAGGCGCCGCGGTGCATGGCAC‑C‑GGCTTTGCCGGTGT‑TCGCGGG‑CACGC‑CCGCTCCCACAGACTGATTGCCATTAAGCGGGCCTATCGCTACAGTGACGG‑CCCGTTTTTTGTTGCCTCGCTAACGAATCATGa               <  8:141977/148‑1 (MQ=37)
                                                                                                            gTGCCCGCGAAACCGGCGCCGCGGTGCATGGCAC‑C‑GGCTTTGCCGGGGTGC‑GCGGG‑CACGC‑CCGCTCCCACAGACTGATTGCCATTAAGCGGGCCTATCGCTACAGTGACGG‑CCCGTTTTTTGTTGCCTCGCTAACGAATCATGACTc            <  1:132138/148‑1 (MQ=25)
                                                                                                               cccGCGAAACAGGCGCCGCGGTGCATGGCAC‑C‑GGCTTTGCCGGTGT‑TCGCGGG‑CACGC‑CCGCTCCCACAGACTGATTGCCATTAAGCGGGCCTATCGCTACAGTGACGG‑CCCGTTTTTTGTTGCCTCGCTAACGAATCATGACTCCCt         <  3:111746/148‑1 (MQ=37)
                                                                                                                ccGCGAAACAGGCGCCGCGGTGCATGGCAC‑C‑GGCTTTGCCGGTGttcgcgggcacgcccgctcccaccgtttggcagcctcagtcctgctgtggtgccaaggcttcgttcaggctctgaatggccgcatcacgggcttcagcgttg              <  5:59607‑M2/146‑103 (MQ=255)
                                                                                                                 cgcgAAACAGGCGCCGCGGTGCATGGCAC‑C‑GGCTTTGCCGGTGttcgcgggcacgcccgctcccaccgtttggcagcctcagtcctgctgtggtgccaaggcttcgttcaggctctgaatggccgcatcacgggcttcagcgttgt             >  6:59607‑M2/1‑43 (MQ=255)
                                                                                                                 cgcgAAACAGGCGCCGCGGTGCATGGCAC‑C‑GGCTTTGCCGGTGT‑TCGCGGG‑CACGC‑CCGCTCCCACAGACTGATTGCCATTAAGCGGGCCTATCGCTACAGTGACGG‑CCCGTTTTTTGTTGCCTCGCTAACGAATCATGACTc            >  2:11582/1‑143 (MQ=37)
                                                                                                                  gcgAAACAGGCGCCGCGGTGCATGGCAC‑C‑GGCTTTGCCGGTGT‑TCGCGGG‑CACGC‑CCGCTCCCACAGACTGATTGCCATTAAGCGGGCCTATCGCTACAGTGACGG‑CCCGTTTTTTGTTGCCTCGCTAACGAATCATGACTcc           <  1:11582/143‑1 (MQ=37)
                                                                                                                       aCAGGCGCCGCGGTGCATGGCAC‑C‑GGCTTTGCCGGTGT‑TCGCGGG‑CACGC‑CCGCTCCCACAGACTGATTGCCATTAAGCGGGCCTATCGCTACAGTGACGG‑CCCGTTTTTTGTTGCCTCGCTAACGAATCATGACTCCCTTAGAACg  >  5:77372/1‑147 (MQ=37)
                                                                                                                                            ||                                                                                                                                  
GACGGGCCGCAGTTCGACAACGCTGAAGCCCGTGATGCGGCCATTCAGAGCCTGAACGAAGCCTTGGCACCACAGCAGGACTGAGGCTGCCAAACGGTGGGAGCGGGCGTGCCCGCGAAACAGGCGCCGCGGTGCATGGCCA‑C‑GGCTTTGCCGGTGT‑T‑TCGGG‑CACG‑‑CCGCTCCCACAGACTGATTGCCATTAAGCGGGCCTATCGCTACAGTGACGG‑CCCGTTTTTTGTTGCCTCGCTAACGAATCATGACTCCCTTAGAACG  >  NC_002947/1499340‑1499603

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: ‑