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,4800CA90.0% 70.6 / ‑3.4 30intergenic (+57/‑107)trpS/zapEtryptophan‑‑tRNA ligase/nucleoside triphosphate hydrolase domain‑containing protein
Reads supporting (aligned to +/- strand):  ref base C (1/2);  new base A (12/15);  total (13/17)
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,4810AC93.3% 77.5 / ‑5.5 30intergenic (+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 (13/15);  total (13/17)
Fisher's exact test for biased strand distribution p-value = 4.92e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

GGGCCGCAGTTCGACAACGCTGAAGCCCGTGATGCGGCCATTCAGAGCCTGAACGAAGCCTTGGCACCACAGCAGGACTGAGGCTGCCAAACGGTGGGAGCGGGCGTGCCCGCGAAACAGGCGCCGCGGTGCATGGCCAC‑GGCTTTGCCGGTGT‑T‑TCGGG‑CA‑CG‑CCGCTCCCACAGACTGATTGCCATTAAGCGGGCCTATCGCTACAGTGACGGCCCGTTTTTTG‑TTGCCTCGCTAACGAATCATGACTCCCTTAGAACGCTATCAAGCAG  >  NC_002947/1499343‑1499614
                                                                                                                                         ||                                                                                                                                            
gggCCGCAGTTCGACAACGCTGAAGCCCGTGATGCGGCCATTCAGAGCCTGAACGAAGCCTTGGCACCACAGCAGGACTGAGGCTGCCAAACGGTGGGAGCGGGCGTGCCCGCGAAACAGGCGCCGCGGTGCATGGCACC‑GGCTTTg                                                                                                                                     <  6:248684/147‑1 (MQ=255)
    cGCAGTTCGACAACGCTGAAGCCCGTGATGCGGCCATTCAGAGCCTGAACGAAGCCTTGGCACCACAGCAGGACTGAGGCTGCCAAACGGTGGGAGCGGGCGTGCCCGCGAAACAGGCGCCGCGGTGCATGGCACC‑GGCTTTGCCgg                                                                                                                                 >  1:138228/1‑147 (MQ=255)
    cGCAGTTCGACAACGCTGAAGCCCGTGATGCGGCCATTCAGAGCCTGAACGAAGCCTTGGCACCACAGCAGGACTGAGGCTGCCAAACGGTGGGAGCGGGCGTGCCCGCGAAACAGGCGCCGCGGTGCATGGCACC‑GGCTTTGCCgg                                                                                                                                 <  2:138228/147‑1 (MQ=255)
             aCAACGCTGAAGCCCGTGATGCGGCCATTCAGAGCCTGAACGAAGCCTTGGCACCACAGCAGGACTGAGGCTGCCAAACGGTGGGAGCGGGCGTGCCCGCGAAACAGGCGCCGCGGTGCATGGCACC‑GGCTTTGCCGgt                                                                                                                                <  8:152545/139‑1 (MQ=255)
             aCAACGCTGAAGCCCGTGATGCGGCCATTCAGAGCCTGAACGAAGCCTTGGCACCACAGCAGGACTGAGGCTGCCAAACGGTGGGAGCGGGCGTGCCCGCGAAACAGGCGCCGCGGTGCATGGCACC‑GGCTTTGCCGgt                                                                                                                                >  7:152545/1‑139 (MQ=255)
                       aGCCCGTGATGCGGCCATTCAGAGCCTGAACGAAGCCTTGGCACCACAGCAGGACTGAGGCTGCCAAACGGTGGGAGCGGGCGTGCCCGCGAAACAGGCGCCGCGGTGCATGGCACC‑CGCTTTGCCGGTGT‑TCGCGGG‑CA‑cg‑cccgc                                                                                                          >  6:79543/1‑144 (MQ=25)
                           cGTGATGCGGCCATTCAGAGCCTGAACGAAGCCTTGGCACCACAGCAGGACTGAGGCTGCCAAACGGTGGGAGCGGGCGTGCCCGCGAAACAGGCGCCGCGGTGCATGGCACC‑GGCTTTGCCGGTGT‑TCGCGGG‑CA‑CGCCCGCTCCca                                                                                                      >  3:238458/1‑148 (MQ=25)
                                    gCCATTCAGCGCCGCAACGAAGCCTTGGCACCACAGCAGGACTGAGGCTGCCAAACGGTGGGAGCGGGCGGGCCCGCGAAACAGGCGCCGCGGTGCAGGGCACC‑GGCTTTGCCGGTGT‑TCGCGGG‑CA‑CGCCCGCTCCCACAGACTGAt                                                                                             <  7:39828/148‑1 (MQ=21)
                                      cATTCAGAGCCTGAACGAAGCCTTGGCACCACAGCAGGACTGAGGCTGCCAAACGGTGGGAGCGGGCGTGCCCGCGAAACAGGCGCCGCGGTGCATGGCACC‑GGCTTTGCCGGTGT‑TCTCGGG‑CA‑CGCCCGCTCCCACAGACTGATTg                                                                                           <  1:30471/148‑1 (MQ=25)
                                                cTGAACGAAGCCTTGGCACCACAGCAGGACTGAGGCTGCCAAACGGTGGGAGCGGGCGTGCCCGCGAAACAGGCGCCGCGGTGCATGGCACC‑GGCGTTGCCGGTGT‑TCGCGGG‑CA‑CGCCCGCTCCCCCAGACTGATTGCCATTAAGCg                                                                                 >  2:20244/1‑148 (MQ=21)
                                                 tGAACGAAGCCTTGGCACCACAGCAGGACTGAGGCTGCCAAACGGTGGGAGCGGGCGTGCCCGCGCAACCGGCGCCCCGGTGCCTGGCCCC‑CGCGTTGCCGGTGT‑TGGCGGG‑CC‑CGCCCGCGCCCACAGCCGGGGTGGCCTTAAGCgg                                                                                >  7:245829/1‑148 (MQ=16)
                                                 tGAACGAAGCCTTGGCACCACAGCAGGACTGAGGCTGCCAAACGGTGGGAGCGGGCGTGCCCGCGAAACAGGCGCCGCGGTGCATGGCACC‑GGCTTTGCCGGTGT‑TCGCGGG‑CA‑CGCCCGCTCCCACAGACTGATTGCCATTAAGCg                                                                                 <  5:101100/147‑1 (MQ=21)
                                                   aaCGAAGCCTTGGCACCACAGCAGGACTGAGGCTGCCAATCGGTGGGAGCGGGGGTGCCCGCGAAACAGGCGCCGCTGTGCATGGCACC‑GGCTTTGCCGGTGT‑TCGCGGG‑CA‑CGCCCGCTCCCACAGACTGATTGCCATTAAGCGGGc                                                                              <  7:74600/148‑1 (MQ=21)
                                                    aCGAAGCCTTGGCACCACAGCAGGACTGAGGCTGCCAAACGGTGGGAGCGGGCGTGCCCGCGAAACAGGCGCCGCGGTGCATGGCACC‑GGCTTTGCCGGTGT‑TCGCGGG‑CA‑CGCCCGCTCCCACAGACTGATTGCCATTAAGCgga                                                                               <  5:79543/146‑2 (MQ=21)
                                                            ttGGCACCACAGCAGGACTGAGGCTGCCAAACGGTGGGAGCGGGCGTGCCCGCGAAACAGGCGCCGCGGTGCATGGCACC‑GGCTTTGCCGGTGT‑TCGCGGG‑CA‑CGCCCGCTCCCACAGACTGATTGCCATTAAGCGGGCCTATCGCTa                                                                     >  2:217022/1‑148 (MQ=21)
                                                                  ccacAGCAGGACTGAGGCTGCCAAACGGTGGGAGCGGGCGTGCCCGCGAAACAGGCGCCCCGGTGCATGGCACC‑CGCATTGCCGGTGT‑TAGCGGGCCC‑CC‑CCCCGCCCCCAGACTGATTGCCAATAAGCCGGCCTATCGCGACAGTg                                                                >  5:285251/1‑147 (MQ=18)
                                                                    acaGCAGGACTGAGGCTGCCAAACGGTGGGAGCGGGCGTGCCCGCGAAACAGGCGCCGCGGTGCATGGCACC‑GGCTTTGCCGGTGT‑TAGCGGG‑CA‑CGCCCGCTCCCACAGACTGATTTCCAATAAGCCGGGCGATCGCGACAGTGACg                                                             >  2:192401/1‑148 (MQ=18)
                                                                          gcaCTGCGGCTCCCAACCGGGGGGAGCGGGCGGCCCCGCTAAACCGGCGCCGCGGGGCA‑GGCCACGGGCTTTGCCGGTGT‑TCGCGGG‑CA‑CGCCCGCTCCCACAGACTGATTGCCATTAAGCGGGCCTATCGCTACAGTGACGGCCCg                                                        <  6:15710/145‑1 (MQ=21)
                                                                               gAGGCTGCCAAACGGTGGGAGCGGGCGTGCCCCCGAAACAGGCGCCCCCGTGCATGGCACC‑CGCGTTGCCGGTGT‑TAGCCGG‑CACCC‑CCCCGCCCACAGACTGATTGCCATTAAACCGGCCGATCGCTACAGTGACGGCCCGttttt                                                   >  6:128226/1‑147 (MQ=255)
                                                                                aGGCTGCCAAACGGTGGGAGCGGGCGTGCCCGCGAAACAGGCGCCGCGGTGCATGGCACC‑GGCTTTGCCGGTGT‑TCGCGGG‑CA‑CGCCCGCTCCCACAGACTGATTGCCATTAAGCGGGCCTATCGCTACAGTGACGGCCCGttttt                                                   <  7:136130/146‑1 (MQ=21)
                                                                                   cTGCCAAACGGGGGGAGCGGGCGGGCCCGCGAAACAGGCGCCGCGGTGCATGGCACC‑GGCTTTGCCGGTGT‑TCGCGGG‑CA‑CGCCCGCTCCCACAGACTGATTGCCATTAAGCGGGCCTATCGCTACAGTGACGGCCCGTTTTttg‑ttg                                             <  3:293817/148‑1 (MQ=21)
                                                                                                      ggCGTGCCCGCGAAACAGGCGCCGCGGTGCATGGCACC‑GGCTTTGCCGGTGT‑TCGCGGG‑CA‑CGCCCGCTCCCACAGACTGATTGCCATTAAGCGGGCCTATCGCTACAGTGACGGCCCGTTTTTTG‑TTGCCTCGCTAACGAATCATGa                          <  4:238458/148‑1 (MQ=37)
                                                                                                       gCGTGCCCGCGAAACAGGCGCCGCGGTGCATGGCACC‑GGCTTTGCCGGTGT‑TCGCGGG‑CA‑CGCCCGCTCCCACAGACTGATTGCCATTAAGCGGGCCTATCGCTACAGTGACGGCCCGTTTTTTG‑TTGCCTCGCTAACGAATCATGAc                         >  7:148470/1‑148 (MQ=37)
                                                                                                        cGTGCCCGCGAAACAGGCGCCGCGGTGCATGGCACC‑GGCTTTGCCGGTGT‑TCGCGGG‑CA‑CGCCCGCTCCCACAGACTGATTGCCATTAAGCGGGCCTATCGCTACAGTGACGGCCCGTTTTTTG‑TTGCCTCGCTAACGAATCATGACt                        >  1:178838/1‑148 (MQ=25)
                                                                                                          tGCCCGCGAAACAGGCGCCGCGGTGCAGGGCACC‑GGCTTGGCCGGTGT‑TCGCGGG‑CA‑CGCCCGCTCCCACAGACTGATTGCCATTAAGCGGGCCTATCGCTACAGTGACGGCCCGTTTTTTG‑TTGCCTCGCTAACGAATCATGACTcc                      <  3:37322/148‑1 (MQ=25)
                                                                                                             ccGCGAAACAGGCGCCGCGGTGCATGGCACC‑GGCTTTGCCGGTGT‑TCGCGGG‑CA‑CGCCCGCTCCCACAGACTGATTGCCATTAAGCGGGCCTATCGCTACAGTGACGGCCCGTTTTTTGCTTGCCTCGCTAACGAATCATGACTCCCt                    <  1:192401/148‑1 (MQ=25)
                                                                                                                ataaaCCGGCGCCGCGGT‑CCTTGCCACGGGCTTGGCGGTTGTGC‑GCGGC‑CA‑CGCCCGCTCCCACAGACTGATTGCCATTAAGCGGGCCTATCGCTACAGTGACGGCCCGTTTTTTG‑TTGCCTCGCTAACGAATCATGACTCCCTTAGa                <  1:318712/146‑1 (MQ=255)
                                                                                                                         cgccgcGGTGCATGGCACC‑GGCTTTGCCGGTGT‑TCGCGGG‑CA‑CGCCCGCTCCCACAGACTGATTGCCATTAAGCGGGCCTATCGCTACAGTGACGGCCCGTTTTTTG‑TTGCCTCGCTAACGAATCATGACTCCCTTAGAACGCTATCa       <  2:52533/148‑1 (MQ=37)
                                                                                                                           ccgcGGTGCATGGCACC‑GGCTTTGCCGGTGT‑TCGCGGG‑CA‑CGCCCGCTCCCACAGACTGATTGCCATTAAGCGGGCCTATCGCTACAGTGACGGCCCGATTTTTG‑TTGCCTCGCTAACGAATAATGACTCCCTTAGAACGCTATCAAg     >  2:126630/1‑148 (MQ=37)
                                                                                                                              gggTGCAGGGCACC‑GGCTTTGCCGTTGT‑TCGCGGG‑CA‑CGCCCGCTCCCACAGACTGATTGCCATTAAGCGGGCCTATCGCTACAGTGACGGCCCGTTTTTTG‑TTGCCTCGCTAACGAATCATGTCTCCCTTAGAACGCTATCATGCAg  <  8:49950/147‑1 (MQ=37)
                                                                                                                                         ||                                                                                                                                            
GGGCCGCAGTTCGACAACGCTGAAGCCCGTGATGCGGCCATTCAGAGCCTGAACGAAGCCTTGGCACCACAGCAGGACTGAGGCTGCCAAACGGTGGGAGCGGGCGTGCCCGCGAAACAGGCGCCGCGGTGCATGGCCAC‑GGCTTTGCCGGTGT‑T‑TCGGG‑CA‑CG‑CCGCTCCCACAGACTGATTGCCATTAAGCGGGCCTATCGCTACAGTGACGGCCCGTTTTTTG‑TTGCCTCGCTAACGAATCATGACTCCCTTAGAACGCTATCAAGCAG  >  NC_002947/1499343‑1499614

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