Predicted mutation
evidence seq id position mutation annotation gene description
RA NC_002947 1,499,508 +C intergenic (+85/‑79) 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,5061.C92.9% 72.7 / 1.1 28intergenic (+83/‑81)trpS/zapEtryptophan‑‑tRNA ligase/nucleoside triphosphate hydrolase domain‑containing protein
Reads supporting (aligned to +/- strand):  ref base . (2/0);  new base C (12/14);  total (14/14)
Fisher's exact test for biased strand distribution p-value = 4.81e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.16e-01

AGCCCGTGATGCGGCCATTCAGAGCCTGAACGAAGCCTTGGCACCACAGCAGGACTGAGGCTGCCAAACGGTGGGAGCGGGCGTGCCCGCGAAACAGGCGCCGCGGTGCATGGCCACGGCTTTGCCGGTGTTTCGGGCACGCCGCTCCCACAGACTGATTGCCATTAAGCGGGCCTATCGCTACAGTGACGGCCCGTTTTTTGTTGCCTCGCTAACGAATCATGACTCCCTTAGAACGCTATCAAGCAGATCTGAAACGTCCCGACTTC  >  NC_002947/1499366‑1499634
                                                                                                                                                  |                                                                                                                                 
aGCCCGTGATGCGGCCATTCAGAGCCTGAACGAAGCCTTGGCACCACAGCAGGACTGAGGCTGCCAAACGGTGGGAGCGGGCGTGCCCGCGAAACAGGCGCCGCGGTGCATGGCACCCGCTTTGCCGGTGTTCGCGGGCAcgcccgc                                                                                                                              >  6:79543/1‑144 (MQ=25)
    cGTGATGCGGCCATTCAGAGCCTGAACGAAGCCTTGGCACCACAGCAGGACTGAGGCTGCCAAACGGTGGGAGCGGGCGTGCCCGCGAAACAGGCGCCGCGGTGCATGGCACCGGCTTTGCCGGTGTTCGCGGGCACGCCCGCTCCca                                                                                                                          >  3:238458/1‑148 (MQ=25)
             gCCATTCAGCGCCGCAACGAAGCCTTGGCACCACAGCAGGACTGAGGCTGCCAAACGGTGGGAGCGGGCGGGCCCGCGAAACAGGCGCCGCGGTGCAGGGCACCGGCTTTGCCGGTGTTCGCGGGCACGCCCGCTCCCACAGACTGAt                                                                                                                 <  7:39828/148‑1 (MQ=21)
               cATTCAGAGCCTGAACGAAGCCTTGGCACCACAGCAGGACTGAGGCTGCCAAACGGTGGGAGCGGGCGTGCCCGCGAAACAGGCGCCGCGGTGCATGGCACCGGCTTTGCCGGTGTTCTCGGGCACGCCCGCTCCCACAGACTGATTg                                                                                                               <  1:30471/148‑1 (MQ=25)
                         cTGAACGAAGCCTTGGCACCACAGCAGGACTGAGGCTGCCAAACGGTGGGAGCGGGCGTGCCCGCGAAACAGGCGCCGCGGTGCATGGCACCGGCGTTGCCGGTGTTCGCGGGCACGCCCGCTCCCCCAGACTGATTGCCATTAAGCg                                                                                                     >  2:20244/1‑148 (MQ=21)
                          tGAACGAAGCCTTGGCACCACAGCAGGACTGAGGCTGCCAAACGGTGGGAGCGGGCGTGCCCGCGCAACCGGCGCCCCGGTGCCTGGCCCCCGCGTTGCCGGTGTTGGCGGGCCCGCCCGCGCCCACAGCCGGGGTGGCCTTAAGCgg                                                                                                    >  7:245829/1‑148 (MQ=16)
                          tGAACGAAGCCTTGGCACCACAGCAGGACTGAGGCTGCCAAACGGTGGGAGCGGGCGTGCCCGCGAAACAGGCGCCGCGGTGCATGGCACCGGCTTTGCCGGTGTTCGCGGGCACGCCCGCTCCCACAGACTGATTGCCATTAAGCg                                                                                                     <  5:101100/147‑1 (MQ=21)
                            aaCGAAGCCTTGGCACCACAGCAGGACTGAGGCTGCCAATCGGTGGGAGCGGGGGTGCCCGCGAAACAGGCGCCGCTGTGCATGGCACCGGCTTTGCCGGTGTTCGCGGGCACGCCCGCTCCCACAGACTGATTGCCATTAAGCGGGc                                                                                                  <  7:74600/148‑1 (MQ=21)
                             aCGAAGCCTTGGCACCACAGCAGGACTGAGGCTGCCAAACGGTGGGAGCGGGCGTGCCCGCGAAACAGGCGCCGCGGTGCATGGCACCGGCTTTGCCGGTGTTCGCGGGCACGCCCGCTCCCACAGACTGATTGCCATTAAGCgga                                                                                                   <  5:79543/146‑2 (MQ=21)
                                     ttGGCACCACAGCAGGACTGAGGCTGCCAAACGGTGGGAGCGGGCGTGCCCGCGAAACAGGCGCCGCGGTGCATGGCACCGGCTTTGCCGGTGTTCGCGGGCACGCCCGCTCCCACAGACTGATTGCCATTAAGCGGGCCTATCGCTa                                                                                         >  2:217022/1‑148 (MQ=21)
                                           ccacAGCAGGACTGAGGCTGCCAAACGGTGGGAGCGGGCGTGCCCGCGAAACAGGCGCCCCGGTGCATGGCACCCGCATTGCCGGTGTTAGCGGGCCCCCCCCCGCCCCCAGACTGATTGCCAATAAGCCGGCCTATCGCGACAGTg                                                                                    >  5:285251/1‑147 (MQ=18)
                                             acaGCAGGACTGAGGCTGCCAAACGGTGGGAGCGGGCGTGCCCGCGAAACAGGCGCCGCGGTGCATGGCACCGGCTTTGCCGGTGTTAGCGGGCACGCCCGCTCCCACAGACTGATTTCCAATAAGCCGGGCGATCGCGACAGTGACg                                                                                 >  2:192401/1‑148 (MQ=18)
                                                   gcaCTGCGGCTCCCAACCGGGGGGAGCGGGCGGCCCCGCTAAACCGGCGCCGCGGGGCAGGCCACGGGCTTTGCCGGTGTTCGCGGGCACGCCCGCTCCCACAGACTGATTGCCATTAAGCGGGCCTATCGCTACAGTGACGGCCCg                                                                            <  6:15710/145‑1 (MQ=21)
                                                        gAGGCTGCCAAACGGTGGGAGCGGGCGTGCCCCCGAAACAGGCGCCCCCGTGCATGGCACCCGCGTTGCCGGTGTTAGCCGGCACCCCCCCGCCCACAGACTGATTGCCATTAAACCGGCCGATCGCTACAGTGACGGCCCGttttt                                                                       >  6:128226/1‑147 (MQ=255)
                                                         aGGCTGCCAAACGGTGGGAGCGGGCGTGCCCGCGAAACAGGCGCCGCGGTGCATGGCACCGGCTTTGCCGGTGTTCGCGGGCACGCCCGCTCCCACAGACTGATTGCCATTAAGCGGGCCTATCGCTACAGTGACGGCCCGttttt                                                                       <  7:136130/146‑1 (MQ=21)
                                                            cTGCCAAACGGGGGGAGCGGGCGGGCCCGCGAAACAGGCGCCGCGGTGCATGGCACCGGCTTTGCCGGTGTTCGCGGGCACGCCCGCTCCCACAGACTGATTGCCATTAAGCGGGCCTATCGCTACAGTGACGGCCCGTTTTttgttg                                                                 <  3:293817/148‑1 (MQ=21)
                                                                               ggCGTGCCCGCGAAACAGGCGCCGCGGTGCATGGCACCGGCTTTGCCGGTGTTCGCGGGCACGCCCGCTCCCACAGACTGATTGCCATTAAGCGGGCCTATCGCTACAGTGACGGCCCGTTTTTTGTTGCCTCGCTAACGAATCATGa                                              <  4:238458/148‑1 (MQ=37)
                                                                                gCGTGCCCGCGAAACAGGCGCCGCGGTGCATGGCACCGGCTTTGCCGGTGTTCGCGGGCACGCCCGCTCCCACAGACTGATTGCCATTAAGCGGGCCTATCGCTACAGTGACGGCCCGTTTTTTGTTGCCTCGCTAACGAATCATGAc                                             >  7:148470/1‑148 (MQ=37)
                                                                                 cGTGCCCGCGAAACAGGCGCCGCGGTGCATGGCACCGGCTTTGCCGGTGTTCGCGGGCACGCCCGCTCCCACAGACTGATTGCCATTAAGCGGGCCTATCGCTACAGTGACGGCCCGTTTTTTGTTGCCTCGCTAACGAATCATGACt                                            >  1:178838/1‑148 (MQ=25)
                                                                                   tGCCCGCGAAACAGGCGCCGCGGTGCAGGGCACCGGCTTGGCCGGTGTTCGCGGGCACGCCCGCTCCCACAGACTGATTGCCATTAAGCGGGCCTATCGCTACAGTGACGGCCCGTTTTTTGTTGCCTCGCTAACGAATCATGACTcc                                          <  3:37322/148‑1 (MQ=25)
                                                                                      ccGCGAAACAGGCGCCGCGGTGCATGGCACCGGCTTTGCCGGTGTTCGCGGGCACGCCCGCTCCCACAGACTGATTGCCATTAAGCGGGCCTATCGCTACAGTGACGGCCCGTTTTTTGCTTGCCTCGCTAACGAATCATGACTCCCt                                        <  1:192401/148‑1 (MQ=25)
                                                                                         ataaaCCGGCGCCGCGGTCCTTGCCACGGGCTTGGCGGTTGTGCGCGGCCACGCCCGCTCCCACAGACTGATTGCCATTAAGCGGGCCTATCGCTACAGTGACGGCCCGTTTTTTGTTGCCTCGCTAACGAATCATGACTCCCTTAGa                                    <  1:318712/146‑1 (MQ=255)
                                                                                                  cgccgcGGTGCATGGCACCGGCTTTGCCGGTGTTCGCGGGCACGCCCGCTCCCACAGACTGATTGCCATTAAGCGGGCCTATCGCTACAGTGACGGCCCGTTTTTTGTTGCCTCGCTAACGAATCATGACTCCCTTAGAACGCTATCa                           <  2:52533/148‑1 (MQ=37)
                                                                                                    ccgcGGTGCATGGCACCGGCTTTGCCGGTGTTCGCGGGCACGCCCGCTCCCACAGACTGATTGCCATTAAGCGGGCCTATCGCTACAGTGACGGCCCGATTTTTGTTGCCTCGCTAACGAATAATGACTCCCTTAGAACGCTATCAAg                         >  2:126630/1‑148 (MQ=37)
                                                                                                       gggTGCAGGGCACCGGCTTTGCCGTTGTTCGCGGGCACGCCCGCTCCCACAGACTGATTGCCATTAAGCGGGCCTATCGCTACAGTGACGGCCCGTTTTTTGTTGCCTCGCTAACGAATCATGTCTCCCTTAGAACGCTATCATGCAg                      <  8:49950/147‑1 (MQ=37)
                                                                                                              tggcaccGGCTTTGCCGGTGTTCGCGGGCACGCCCGCTCCCACAGACTGATTGCCATTAAGCGGGCCTATCGCTACAGTGACGGCCCGTTTTTTGTTGCCTCGCTAACGAATCATGACTCCCTTAGAACGCTATCAAGCAGATCTg                 >  5:224010/7‑146 (MQ=38)
                                                                                                                        cTTTGCCGGTGTTCGCGGGCACGCCCGCTCCCACAGACTGATTGCCATTAAGCGGGCCTATCGCTACAGTGACGGCCCGTTTTTTGTTGCCTCGCTAACGAATCATGACTCCCTTAGAACGCTATCAAGCAGATCTGAAACGTCCCGa      >  1:239821/1‑148 (MQ=39)
                                                                                                                         tttGCCGGTGTTCGCGGGCACGCCCGCTCCCACAGACTGATTGCCATTAAGCGGGCCTATCGCTACAGAGACGGCCCGTTTTTTGTTGCCTCGCTAACGAATCATGACTCCCTTAGAACGCTATCAAGCAGATCTGAAACGTcccta      >  5:174283/1‑145 (MQ=39)
                                                                                                                             ccGGTGTTCGCGGGCACGCCCGCTCCCACAGACTGATTGCCATTAAGCGGGCCTATCGCTACAGTGACGGCCCGATTTTTGTTGCCTCGCTAACGAATCATGACTGCCTTAGCACGCTATCAAGCAGATCTGAAAAGTCCCGActtc  >  6:151022/1‑147 (MQ=255)
                                                                                                                                                  |                                                                                                                                 
AGCCCGTGATGCGGCCATTCAGAGCCTGAACGAAGCCTTGGCACCACAGCAGGACTGAGGCTGCCAAACGGTGGGAGCGGGCGTGCCCGCGAAACAGGCGCCGCGGTGCATGGCCACGGCTTTGCCGGTGTTTCGGGCACGCCGCTCCCACAGACTGATTGCCATTAAGCGGGCCTATCGCTACAGTGACGGCCCGTTTTTTGTTGCCTCGCTAACGAATCATGACTCCCTTAGAACGCTATCAAGCAGATCTGAAACGTCCCGACTTC  >  NC_002947/1499366‑1499634

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: