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,4800CA97.6% 121.3 / NA 42intergenic (+57/‑107)trpS/zapEtryptophan‑‑tRNA ligase/nucleoside triphosphate hydrolase domain‑containing protein
Reads supporting (aligned to +/- strand):  ref base C (0/1);  new base A (20/21);  total (20/22)
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,4810AC97.6% 122.7 / ‑5.0 42intergenic (+58/‑106)trpS/zapEtryptophan‑‑tRNA ligase/nucleoside triphosphate hydrolase domain‑containing protein
Reads supporting (aligned to +/- strand):  ref base A (0/1);  new base C (20/21);  total (20/22)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.99e-01

GCCGACGGGCCGCAGTTCGACAACGCTGAAGCCCGTGATGCGGCCATTCAGAGCCTGAACGAAGCCTTGGCACCACAGCAGGACTGAGGCTGCCAAACGGTGGGAGCGGGCGTGCCCGCGAAACAGGCGCCGCGGTGCATGGCCAC‑GGCTTTGCCGGTGT‑T‑TCGGG‑CA‑CG‑CCGCTCCCACAGACTGATTGCCATTAAGCGGGCCTATCGCTACAGTGACGGCCCGTTTTTTGTTGCCTCGCTAACGAATCATGACTCCCTTAGAACGCTATCAAGCAGAT  >  NC_002947/1499337‑1499616
                                                                                                                                               ||                                                                                                                                             
gCCGACGGGCCGCAGTTCGACAACGCTGAAGCCCGTGATGCGGCCATTCAGAGCCTGAACGAAGCCTTGGCACCACAGCAGGACTGAGGCTGCCAAACGGTGGGAGCGGGCGTGCCCGCGAAACAGGCGCCGCGGTGCATGGCACC‑GGc                                                                                                                                          >  7:138617/1‑149 (MQ=255)
       ggCCGCAGTTCGACAACGCTGAAGCCCGTGATGCGGCCATTCAGAGCCTGAACGAAGCCTTGGCACCACAGCAGGACTGAGGCTGCCAAACGGTGGGAGCGGGCGTGCCCGCGAAACAGGCGCCGCGGTGCATGGCACC‑GGCTTTGCCg                                                                                                                                   <  1:303365/149‑1 (MQ=255)
       ggCCGCAGTTCGACAACGCTGAAGCCCGTGATGCGGCCATTCAGAGCCTGAACGAAGCCTTGGCACCACAGCAGGACTGAGGCTGCCAAACGGTGGGAGCGGGCGTGCCCGCGAAACAGGCGCCGCGGTGCATGGCACC‑GGCTTTGCCg                                                                                                                                   <  8:42916/149‑1 (MQ=255)
        gCCGCAGTTCGACAACGCTGAAGCCCGTGATGCGGCCATTCAGAGCCTGAACGAAGCCTTGGCACCACAGCAGGACTGAGGCTGCCAAACGGTGGGAGCGGGCGTGCCCGCGAAACAGGCGCCGCGGTGCATGGCACC‑GGCTTTGCCgg                                                                                                                                  <  7:77153/149‑1 (MQ=255)
         ccGCAGTTCGACAACGCTGAAGCCCGTGATGCGGCCATTCAGAGCCTGAACGAAGCCGTGGCACCACAGCAGGACTGAGGCTGCCAAACGGTGGGAGCGGGCGTGCCCGCGAAACAGGCGCCGCGGTGCATGGCACC‑GGCTTTGCCGgt                                                                                                                                 <  1:90803/149‑1 (MQ=255)
           gCAGTTCGACAACGCTGAAGCCCGTGATGCGGCCATTCAGAGCCTGAACGAAGCCTTGGCACCACAGCAGGACTGAGGCTGCCAAACGGTGGGAGCGGGCGTGCCCGCGAAACAGGCGCCGCGGTGCATGGCACC‑GGCTTTGCCGgtgt                                                                                                                               <  1:117682/149‑1 (MQ=255)
            catttCGACATCGGTGATGCCCGTGTTGCGGCCATTCTGAGCCTGATGGTTGGCTTGGCACCACTGCAGGACTGAGGCTGCCAATCGGTGGGAGCGGGCGTGCCCGCGAAACAGGCGCCGCGGTGCATGGCACC‑GGCTTTGCCGgtgt                                                                                                                               <  2:201743/145‑1 (MQ=21)
                 cGACAACGCTGAAGCCCGTGATGCGGCCATTCAGAGCCTGAACGAAGCCTTGGCACCACAGCAGGACTGAGGCTGCCAAACGGTGGGAGCGGGCGTGCCCGCGAAACAGGCGCCGCGGTGCATGGCACC‑GGCTTTGCCg                                                                                                                                   <  7:281927/139‑1 (MQ=255)
                 cGACAACGCTGAAGCCCGTGATGCGGCCATTCAGAGCCTGAACGAAGCCTTGGCACCACAGCAGGACTGAGGCTGCCAAACGGTGGGAGCGGGCGTGCCCGCGAAACAGGCGCCGCGGTGCATGGCACC‑GGCTTTGCCg                                                                                                                                   >  8:281927/1‑139 (MQ=255)
                  gACAACGCTGAAGCCCGTGATGCGGCCATTCAGAGCCTGAACGAAGCCTTGGCACCACAGCAGGACTGAGGCTGCCAAACGGTGGGAGCGGGCGTGCCCGCGAAACAGGCGCCGCGGTGCATGGCACC‑GGCTTTGCCGGTGt‑tcgcgg                                                                                                                        <  8:138617/148‑6 (MQ=255)
                   aCAACGCTGAAGCCCGTGATGCGGCCATTCAGAGCCTGAACGAAGCCTTGGCACCACAGCAGGACTGAGGCTGCCAAACGGTGGGAGCGGGCGTGCCCGCGAAACAGGCGCCGCGGTGCATGGCACC‑GGCTTTGCCGGTGt‑tcgcgg                                                                                                                        <  6:215732/147‑6 (MQ=255)
                   aCAACGCTGAAGCCCGTGATGCGGCCATTCAGAGCCTGAACGAAGCCTTGGCACCACAGCAGGACAGAGGCTGCCAAACGGTGGGAGCGGGCGTGCCCGCGAAACAGGCGCCGCGGTGCATGGCACC‑CGCTTTGCCGGTGt‑tcgcgggc                                                                                                                      >  2:101025/1‑142 (MQ=37)
                        gCTGAAGCCCGTGATGCGGCCATTCAGAGCCTGAACGAAGCCTTGGCACCACAGCAGGACTGAGGCTGCCAAACGGTGGGAGCGGGCGTGCCCGCGAAACAGGCGCCGCGGTGCATGGCACC‑GGCTTTGCCGGTGT‑TCGCGGG‑CA‑cg‑cc                                                                                                              <  1:380406/149‑1 (MQ=37)
                         cTGAAGCCCGTGATGCGGCCATTCAGAGCCTGAACGAAGCCTTGGCACCACAGCAGGACTGAGGCTGCCAAACGGTGGGAGCGGGCGTGCCCGCGAAACAGGCGCCGCGGTGCATGGCACC‑GGCTTTGCCGGTGT‑TCGCGGG‑CA‑cg‑ccc                                                                                                             <  8:89265/149‑2 (MQ=37)
                                                    gCCTGAACGAAGCCTTGGCACCACAGCAGGACTGAGGCTGCCAAACGGTGGGAGCGGGCGTGCCCGCGAAACAGGCGCCGCGGTGCATGGCACC‑GGCTTTGCCGGTGT‑TCGCGGG‑CA‑CGCCCGCTCCCACAGACTGATTGCCATTAAg                                                                                    <  4:107731/148‑1 (MQ=21)
                                                     ccTGAACGAAGCCTTGGCACCACAGCAGGACTGAGGCTGCCAAACGGTGGGAGCGGGCGTGCCCCCCAAACACGCCCCGCGGTGCATGGCACC‑GGCATTGGCGGTGT‑TCGCGGG‑CA‑cg‑cccgcg                                                                                                          >  7:202250/1‑120 (MQ=21)
                                                         aaCGAAGCCTTGGCACCACAGCAGGACTGAGGCTGCCAAACGGTGGGAGCGGGCGTGCCCGCGAAACAGGCGCCGCGGTGCATGGCACC‑GGCTTTGCCGGTGT‑TCGCGGG‑CA‑CGCCCGCTCCCACAGACTGATTGc                                                                                           <  3:265367/136‑1 (MQ=21)
                                                         aaCGAAGCCTTGGCACCACAGCAGGACTGAGGCTGCCAAACGGTGGGAGCGGGCGTGCCCGCCAAACAGGCGCCGCGGTGCATGGCACC‑GGCTTTGCCGGTGT‑TAGCGGG‑CACCC‑CCGCTCCCACACACAGATTGc                                                                                           >  4:265367/1‑136 (MQ=18)
                                                            gAAGCCTTGGCACCACAGCAGGACTGAGGCTGCCAAACGGTGGGAGCGGGCGTGCCCGCGAAACAAGCGCCGCCGTGCATTGCACC‑CGCATTTCCCGTGT‑TAGCGGGCCA‑GC‑CCGCGCCCCCACACTGATTTGCATTAAGCggggca                                                                             >  6:138420/1‑144 (MQ=12)
                                                                                      aGGCTGCCAAACGGTGGGAGCGGGCGTGCCCGCGAAACAGGCGCCGCGGTGCATGGCACC‑GGCTTTGCCGGTGT‑TCGCGGG‑CA‑CGCCCGCTCCCACAGACTGATTGCCATTAAGCGGGCCTATAGCTACACTGACGGCCCGttttttt                                                  >  6:201827/1‑147 (MQ=21)
                                                                                               aaCGGGGGGAGCGGGCGTGCCCGCGAAACAGGCGCCGCGGTGCATGGCACC‑GGCTTTGCCGGTGT‑TCGCGGG‑CA‑CGCCCGCTCCCACAGACTGATTGCCATTATGCGGGCCTATCGCTACAGTGACGGCCCGTTTTTTGTTGCCTCGCt                                        <  1:279747/149‑1 (MQ=21)
                                                                                                aCGGTGGGAGCGGGCGTGCCCGCGAAACAGGCGCCGCGGTGCATGGCACC‑GGCTTTGCCGGTGT‑TCGCGGG‑CC‑CGCCCGCTCCCACAGACTGATTGGCATTAAGCGGGCCTATCGCGACAGTGACGGCCCGTTTTTTGTTGGCTCggt                                        >  3:68545/1‑146 (MQ=18)
                                                                                                            ggCGTGACCGCTAAACCGCCCCCGGGGT‑CCTGGCCACGGGCTTTGCCGGTGTGC‑GCGGG‑AA‑CGCCCGCTCCCACAGACTGATTGCCATTAAGCGGGCCTATCGCTACAGTGACGGCCCGTTTTTTGTTGCCTCGCTAACGAATCATGAc                           <  7:145113/149‑1 (MQ=255)
                                                                                                             gggtccccGCTAAACAGGCGCGGGGGTGATTGGCACG‑GGGTTGGCGGGTGTGC‑GCGGG‑CA‑GGCCCGCTCCCACAGACTGATTGCCATTAAGCGGGCCTATCGCTACAGTGACGGCCCGTTTTTTGTTGCCTCGCTAACGAATCATGAc                           <  6:31236/143‑1 (MQ=255)
                                                                                                             gCGTGCCCGCGAAACAGGCGCCGCGGTGCATGGCACC‑GGCTTTGCCGGTGT‑TCGCGGG‑CA‑CGCCCGCTCCCACAGACTGATTGCCATTAAGCGGGCCTATCGCTACAGTGACGGCCCGTTTTTTGTTGCCTCGCTAACGAATCATGACt                          >  1:243453/1‑149 (MQ=37)
                                                                                                                   ccGCGAAACAGGCGCCGCGGTGCATGGCACC‑GGCTTTGCCGGTGT‑TCGCGGG‑CA‑CGCCCGCTCCCACAGACTGATTGCCATTAAGCGGGCCTATCGCTACAGTGACGGCCCGTTTTTTGTTGCCTCGCTAACGAATCATGACTCCCt                      >  2:206154/1‑147 (MQ=37)
                                                                                                                   ccGCGAAACAGGCGCCGCGGTGCATGGCACC‑GGCTTTGCCGGTGT‑TCGCGGG‑CA‑CGCCCGCTCCCACAGACTGATTGCCATTAAGCGGGCCTATCGCTACAGTGACGGCCCGTTTTTTGTTGCCTCGCTAACGAATCATGACTCCCt                      >  3:351585/1‑147 (MQ=37)
                                                                                                                      ctaaaCAGGCGCCGCGGTGCAGGGCACC‑GGCTTTGCCGGTGT‑TCGCGGG‑CA‑CGCCCGCTCCCACACACTGATTGCCATTAAGCGGGCCTATCGCTACAGTGACGGCCCGTTTTTTGTTGCCTCGCTAACGAATCATGACTCCCTTAGaa                 <  1:347395/147‑1 (MQ=37)
                                                                                                                      cgAAACAGGCGCCGCGGTGCATGGCACC‑GGCTTTGCCGGTGT‑TCGCGGG‑CA‑CGCCCGCTCCCACAGACTGATTGCCATTAAGCGGGCCTATCGCTACAGTGACGGCCCG‑TTTTTGTTGCCTCGCTAACGAATCATGACTCCCTTAGaa                 <  5:201827/148‑1 (MQ=37)
                                                                                                                          aCAGGCGCCGCGGTGCATGGCACC‑GGCTTTGCCGGTGT‑TCGCGGG‑CA‑CGCCCGCTCCCACAGACTGATTGCCATTAAGCGGGCCTATCGCTACAGTGACGGCCCGTTTTTTGTTGCCTCGCTAACGAATCATGACTCCCTTAGAACg               >  6:281204/1‑147 (MQ=37)
                                                                                                                            aGGCGCCGCGGTGCATGGCACC‑GGCTTTGCCGGTGT‑TCGCGGG‑CA‑CGCCCGCTCCCACAGACTGATTGCCATTAAGCGGGCCTATCGCTACAGTGACGGCCCGTTTTTTGTTGCCTCGCTAACGAATCATGACTCCCTTAGAACGCTAt           >  8:160803/1‑149 (MQ=37)
                                                                                                                             ggCGCCGCGGTGCATGGCACC‑GGCTTTGCCGGTGT‑TCGCGGG‑CA‑CGCCCGCTCCCACAGACTGATTGCCATTAAGCGGGCCTATCGCTACAGTGACGGCCCGTTTTTTGTTGCCTCGCTAACGAATCATGACTCCCTTAGAACGCTAt           >  2:287447/1‑148 (MQ=37)
                                                                                                                              gcgcCGCGGTGCATGGCACC‑GGCTTTGCCGGTGT‑TCGCGGG‑CA‑CGCCCGCTCCCACAGACTGATTGCCATTAAGCGGGCCTAt                                                                           >  4:174217/1‑83 (MQ=14)
                                                                                                                              gcgcCGCGGTGCATGGCACC‑GGCTTTGCCGGTGT‑TCGCGGG‑CA‑CGCCCGCTCCCACAGACTGATTGCCATTAAGCGGGCCTAt                                                                           <  3:174217/83‑1 (MQ=14)
                                                                                                                              gcgcCGCGGTGCATGGCACC‑GGCTTTGCCGGTGT‑TCGCGGG‑CA‑CGCCCGCTCCCACAGACTGATTGCCATTAAGCGGGCCTAt                                                                           <  7:85736/83‑1 (MQ=14)
                                                                                                                              gcgcCGCGGTGCATGGCACC‑GGCTTTGCCGGTGT‑TCGCGGG‑CA‑CGCCCGCTCCCACAGACTGATTGCCATTAAGCGGGCCTAt                                                                           >  8:85736/1‑83 (MQ=14)
                                                                                                                              gcgcCGCGGTGCATGGCACC‑GGCTTTGCCGGTGT‑TCGCGGG‑CA‑CGCCCGCTCCCACAGACTGATTGCCATTAAGCGGGCCTATCGCTACAGTGACGGCCCGTTTTTTGTTGCCTCGCTAACGAATCATGACTCCCTTAGaa                 <  3:244519/141‑1 (MQ=37)
                                                                                                                              gcgcCGCGGTGCATGGCACC‑GGCTTTGCCGGTGT‑TCGCGGG‑CA‑CGCCCGCTCCCACAGACTGATTGCCATTAAGCGGGCCTATCGCTACAGTGACGGCCCGTTTTTTGTTGCCTCGCTAACGAATCATGACTCCCTTAGaa                 >  4:244519/1‑141 (MQ=37)
                                                                                                                              gcgcCGCGGTGCATGGCACC‑GGCTTTGCCGGTGT‑TCGCGGG‑CA‑CGCCCGCTCCCACAGACTGATTGCCATTAAGCGGGCCTATCGCTACAGTGACGGCCCGTGTTTTGTTGCCTCGCTAACGAATCATGACTCCCTTAGAACGCTATc          >  7:286241/1‑148 (MQ=37)
                                                                                                                                 ccgcGGTGCATGGCACC‑GGCTTTGCCGGTGT‑TCGCGGG‑CA‑CGCCCGCTCCCACAGACTGATTGCCATTAAGCGGGCCTATCGCTACAGTGACGGCCCGTTTTTTGTTGCCTCGCTAACGAATCATGACTCCCTTAGAACGCTATCAAGc      <  2:243453/149‑1 (MQ=37)
                                                                                                                                   gcgGTGCATGGCACC‑GGCTTTGCCGGTGT‑TCGCGGG‑CA‑CGCCCGCTCCCACAGACTGATTGCCATTAAGCGGGCCTATCGCTACAGTGACGGCCCGTTTTTTGTTGCCTCGCTAACGAATCATGACTCCCTTAGAACGCTATCAAg       >  8:215462/1‑146 (MQ=37)
                                                                                                                                      gTGCATGGCACC‑GGCTTTGCCGGTGT‑TCGCGGG‑CA‑CGCCCGCTCCCACAGACTGATTGCCATTAAGCGGGCCTATCGCTACAGTGACGGCCCGTTTTTTGTTGCCTCGCTAACGAATCATGACTCCCTTAGAACGCTATCAAGCAGAt  >  6:242737/1‑148 (MQ=38)
                                                                                                                                               ||                                                                                                                                             
GCCGACGGGCCGCAGTTCGACAACGCTGAAGCCCGTGATGCGGCCATTCAGAGCCTGAACGAAGCCTTGGCACCACAGCAGGACTGAGGCTGCCAAACGGTGGGAGCGGGCGTGCCCGCGAAACAGGCGCCGCGGTGCATGGCCAC‑GGCTTTGCCGGTGT‑T‑TCGGG‑CA‑CG‑CCGCTCCCACAGACTGATTGCCATTAAGCGGGCCTATCGCTACAGTGACGGCCCGTTTTTTGTTGCCTCGCTAACGAATCATGACTCCCTTAGAACGCTATCAAGCAGAT  >  NC_002947/1499337‑1499616

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