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.C97.3% 103.2 / ‑2.9 37intergenic (+83/‑81)trpS/zapEtryptophan‑‑tRNA ligase/nucleoside triphosphate hydrolase domain‑containing protein
Reads supporting (aligned to +/- strand):  ref base . (1/0);  new base C (15/21);  total (16/21)
Fisher's exact test for biased strand distribution p-value = 4.32e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 7.13e-01

GCTGAAGCCCGTGATGCGGCCATTCAGAGCCTGAACGAAGCCTTGGCACCACAGCAGGACTGAGGCTGCCAAACGGTGGGAGCGGGCGTGCCCGCGAAACAGGCGCCGCGGTGCATGGCCACGGCTTTGCCGGTGTTTCGGGCACGCCGCTCCCACAGACTGATTGCCATTAAGCGGGCCTATCGCTACAGTGACGGCCCGTTTTTTGTTGCCTCGCTAACGAATCATGACTCCCTTAGAACGCTATCAAGCAGATCTGAAACGTCCCGACTTCTTCCATGAC  >  NC_002947/1499361‑1499643
                                                                                                                                                     |                                                                                                                                           
gCTGAAGCCCGTGATGCGGCCATTCAGAGCCTGAACGAAGCCTTGGCACCACAGCAGGACTGAGGCTGCCAAACGGTGGGAGCGGGCGTTCCCCCGCAAAAGGCGCCCCCGTGCATTGCACCCGCTATTCCGGTGTTACCGGGCAcccc                                                                                                                                           >  6:134239/1‑149 (MQ=255)
         cGTGATGCGGCCATTCAGAGCCTGAACGAAGCCTTGGCACCACAGCAGGACTGAGGCTGCCAAACGGTGGGAGCGGGCGTGCCCGCGCAAAAAGCGCCCCCGTGCATGGCACCGGCTTTGCCGGTGTTCGCGGGCCCGcccccccccc                                                                                                                                   >  4:336869/1‑146 (MQ=21)
               gCGGCCATTCAGAGCCTGAACGAAGCCTTGGCACCACAGCAGGACTGAGGCTGCCAAACGGTGGGAGCGGGCGTGCCCGCGAAACAGGCGCCGCGGTGCATGGCACCGGCTTTGCCGGTGTTCGCGGGCACGCCCGCTCCCACAGACt                                                                                                                              >  7:276278/1‑148 (MQ=21)
                 ggCCATTCAGAGCCTGAACGAAGCCTTGGCACCACAGCAGGACTGAGGCTGCCAAACGGTGGGAGCGGGCGTGCCCGCGAAACAGGCGCCGCGGTGCATGGCACCGGCTTTGCCGGTGTTCGCGGGCACGCCCGCTCCCACAGACTGAt                                                                                                                          >  5:267629/1‑149 (MQ=21)
                  gCCTTTCAGAGCCTGAACGAAGCGTTGGCACCAACGCAGGACTGAGGCTGCCAAACGGGGGGCGCGGGCGTGCCCGCTAAACAGGCGCCGCGGTGCTGGGCACGGGCTTTGCGGGTGTTCGCGGGCACGCCCGCTCCCACAGACTGAtt                                                                                                                         <  8:221099/149‑1 (MQ=21)
                  gCCATTCAGAGCCTGAACGAAGCCTTGGCACCACAGCAGGACTGAGGCTGCCAAACGGTGGGAGCGGGCGTGCCCGCGAAACAGGCGCCGCGGTGCATGGCACCGGCTTTGCCGGTGTTCGCGGGCACGCCCGCTCCCACAGACTGAtt                                                                                                                         <  8:276278/149‑1 (MQ=21)
                    cATTCAGAGCCTGAACGAAGCCTTGGCACCACAGCAGGACTGAGGCTGCCAAACGGTGGGAGCGGGCGTGCCCGCGAAACAGGCGCCGCGGTGCATGGCACCGGCTTTGCCGGTGTTCGCGGGCACGCCCGCTCCCACAGACTGATTGc                                                                                                                       <  3:84971/149‑1 (MQ=21)
                         agagCCTCAAGTCAGCCTTGCCACCACTCCCGTAGTGAGGCTCCAAAAGGGTGGCTGCGGGCGGGCCCGCTAAACTGGCGCCGGGGTGATTGCCCCCGGCTTTGCGGGTGTTCGGGGGCAGGCCCGCTCCCACAGACTGATTGCCAtt                                                                                                                   <  8:17180/148‑1 (MQ=255)
                             ccTGACCGAAGCCTTGGCACCACAGCAGGACTGAGGCTGCCAAACGGTGGGAGCGGGCGGGCCCGCGAAACAGGCGCCGCGGTGCATGGCACCGGCTTTGCCGGTGTTCGCGGGCACGCCCGCTCCCACAGACTGATTGCCa                                                                                                                     <  7:143567/142‑1 (MQ=21)
                             ccTGAACGAAGCCTTGGCACCACAGCAGGACTGAGGCTGCCAAACGGTGGGAGCGGGCGTGCCCGCGAAACAGGCGCCGCGGTGCATGGCACCGGCTTTGCCGGTGTTCGCGGGCACGCCCGCTCCCACAGACTGATTGCCAtt                                                                                                                   >  6:265611/1‑144 (MQ=21)
                             ccTGAACGAAGCCTTGGCACCACAGCAGGACTGAGGCTGCCAAACGGTGGGAGCGGGCGTGCCCGCGAAACAGGCGCCGCGGTGCATGGCACCGGCTTTGCCGGTGTTCGCGGGCACGCCCGCTCCCACAGACTGATTGCCAtt                                                                                                                   <  5:265611/144‑1 (MQ=21)
                             ccTGAACGAAGCCTTGGCACCACAGCAGGACTGAGGCTGCCAAACGGTGGGAGCGGGCGTGCCCGCGAAACAGGCGCCGCGGTGCATGGCACCGGCTTTGCCGGTGTTCGCGGGCACGCCCGCTCCCACAGACTGATTGCCATTAAGCg                                                                                                             <  5:207128/149‑1 (MQ=21)
                               tGAACGAAGCCTTGGCACCACAGCAGGACTGAGGCTGCCAAACGGTGGGAGCGGGCCTGGCCCCGAAACACGGGCCGCGGTGCAATGCACCCGCGTTGGCGGTTTTAGCGGGCACGCCCGCACCCACACCCAGAATTGCCATAAGCggg                                                                                                           >  6:322620/1‑149 (MQ=255)
                                                      cagGACTGAGGCTGCCAAACGGTGGGAGCGGGCGTGCCCGCGAAACAGGCGCCGCGGTGCATGGCACCGGCTTTGCCGGTGTTAGCGGGCACGCCCGCTCCCACAGACTGATTGCCATTAAGCGGGCCTATCGCTACAGTGACGGCCCg                                                                                    >  2:267452/1‑149 (MQ=25)
                                                         gACTGAGGCTGCCAAACGGTGCGAGCGGGCGTGCCCGCGAAACAAGCGCCCCGGTGCATGGCACCGGCTTTGCCCCTGTTGGCGGGCCCCCCCGCGCCCAAAGACAGAATTCCATTATACCGCCCAATAGCTTCACTGACGGCCCGgtt                                                                                 >  8:203971/1‑149 (MQ=9)
                                                              tggCTGCCATACGTGGGTAGGGGGCGTGCCCGCGAAACATGCGCGGCGGTGCTTGGCACGGGCTTTGCCGGTGTTCGCGGGCACGCCCGCTCCCACAGACTGATTGCCATTAAGCGGGCCTATCGCTACAGTGACGGCCCGtttttt                                                                              <  7:203971/146‑1 (MQ=21)
                                                                gCTGCCAAACGGTGGCTGCGGGCGTGCCCGCGAAACAGGCGCGGCGGTGCATGGCACCGGCTTTGCCGGTTTTCGCGGTCACGCCCGCTCCCACAGACTGATTGCCATTAAGCGGGCCTATCGCTACAGTGACGGCCCGTTTTttgttg                                                                          <  6:279518/149‑1 (MQ=21)
                                                                         cGGTGGGAGCGGGCGTGCCCGCGAAACAGGCGCCGCGGTGCATGGCACCGGCTTTGCCGGTGTTCGCGGGCACGCCCGCTCCCACAGACTGATTGCCATTAAGCGGGCCTATCGCTACAGTGACGGCCCGTTTTTTGTTGCCTCGCTaa                                                                 <  3:336869/149‑1 (MQ=21)
                                                                                aGCGGGCGTGCCCGCGAAACAGGCGCCGCGGTGCATGGCACCGGCTTTGCCGGTGTTCGCGGGCACGCCCGCTCCCACAGACTGATTGCCATTAAGCGGGCCTATCGCTACAGTGACGGCCCGTTTTTTGTTGCCTCGCTAACGAATCa                                                          <  6:267629/149‑1 (MQ=25)
                                                                                 gggggCGTGCCGGGGAAACTGGCGCGGCGGTGCTTGGCAGCGGCTTTGCGGGTGTTCGCGGTCACGCCCGCTCCCACAGACTGATTGCCATTAATGCGGGCCTATCGCTACAGTGACGGCCCGTTTTTTGTTGCCTCGCTAACGAATc                                                           <  1:97227/146‑1 (MQ=255)
                                                                                     gCGTGCCCGCGAAACAGGCGCCGCGGTGCATGGCACCGGCTTTGCCGGTGTTCGCGGGCACGCCCGCTCCCACAGACTGATTGCCATTAAGCGGGCCTATCGCTACAGTGACGGCCCGTTTTTTGTTGCCTCGCTAACGAATCATGACt                                                     >  8:161259/1‑149 (MQ=25)
                                                                                          cccGCGAAACAGGCGCCGCGGTGCATGGCACCGGCTTTGCCGGTGTTCGCGGGCACGCCCGCTCCCACAGACTGATTGCCATTAAGCGGGCCTATCGCTACAGTGACGGCCCGTTTTTTGTTGCCTCGCTAACGAATCATGACTCCCt                                                 >  1:181507/1‑148 (MQ=37)
                                                                                           ccGCGAAACAGGCGCCGCGGTGCATGGCACCGGCTTTGCCGGTGTTCGCGGGCACGCCCGCTCCCACAGACTGATTGCCATTAAGCGGGCCTATCGCTACAGTGACGGCCCGTTTTTTGTTGCCTCGCTAACGAATCATGACTCCCt                                                 >  2:321829/1‑147 (MQ=37)
                                                                                                       cgccgcGGTGCATGGCACCGGCTTTGCCGGTGTTCGCGGGCACGCCCGCTCCCACAGACTGATTGCCATTAAGCGGGCCTATCGCTACAGTGACGGCCCGTTTTTTGTTGCCTCGCTAACGAATCATGACTCCCTTAGAACGCTATCaa                                   <  1:321829/149‑1 (MQ=37)
                                                                                                        gccgcGGTGCATGGCACCGGCTTTGCCGGTGTTCGCGGGCACGCCCGCTCCCACAGACTGATTGCCATTAAGCGGGCCTATCGCTACAGTGACGGCCCGTTTTTTGTTGCCTCGCTAACGAATCATGACTCCCTTAGAACGCTATCAAg                                  <  1:267452/149‑1 (MQ=37)
                                                                                                                        ccGGCTTTGCCGGTGTTCGCGGGCACGCCCGCTCCCACAGACTGATTGCCATTAAGCGGGCCTATCGCTACAGTGACGGCCCGTTTTTTGTTGCCTCGCTAACGAATCATGACTCCCTTAGAACGCTATCAAGCAGATCTGAAACGTcc                  >  1:146356/2‑149 (MQ=39)
                                                                                                                           ggCTTTGCCGGTGTTCGCGGGCACGCCCGCTCCCACAGACTGATTGCCATTAAGCGGGcc                                                                                                         >  1:345355/1‑60 (MQ=14)
                                                                                                                           ggCTTTGCCGGTGTTCGCGGGCACGCCCGCTCCCACAGACTGATTGCCATTAAGCGGGcc                                                                                                         <  2:345355/60‑1 (MQ=14)
                                                                                                                            gCTTTGCCGGTGTTCGCGGGCACGCCCGCTCCCACAGACTGATTGCCATTAAGCGGGCCTATCGCTACAGTGACGGCCCGTTTTTTGTTGCCTCGCTAACGAATCATGACTCCCTTAGAACGCTATCAAGc                                 >  4:35079/1‑131 (MQ=38)
                                                                                                                            gCTTTGCCGGTGTTCGCGGGCACGCCCGCTCCCACAGACTGATTGCCATTAAGCGGGCCTATCGCTACAGTGACGGCCCGTTTTTTGTTGCCTCGCTAACGAATCATGACTCCCTTAGAACGCTATCAAGc                                 <  3:35079/131‑1 (MQ=38)
                                                                                                                            gCTTTGCCGGTGTTCGCGGGCACGCCCGCTCCCACAGACTGATTGCCATTAAGCGGGCCTATCGCTACAGTGACGGCCCGTTTTTTGTTGCCTCGCTAACGAATCATGACTCCCTTAGAACGCTATCAAGCAGATCTGAAACGTCCCGa               >  7:192728/1‑149 (MQ=39)
                                                                                                                              tttGCCGGTGTTCGCGGGCACGCCCGCTCCCACAGACTGATTGCCATTAAGCGGGCCTATCGCTACAGTGACGGCCCGTTTTTTGTTGCCTCGCTAACGAATCa                                                          <  1:110545/104‑1 (MQ=37)
                                                                                                                              tttGCCGGTGTTCGCGGGCACGCCCGCTCCCACAGACTGATTGCCATTAAGCGGGCCTATCGCTACAGTGACGGCCCGTTTTTTGTTGCCTCGCTAACGAATCa                                                          >  2:110545/1‑104 (MQ=37)
                                                                                                                                  ccGGTGTTCGCGGGCACGCCCGCTCCCACAGACTGATTGCCATTAAGCGGGCCTATCGCTACAGTGACGGCCCGTTTTTTGTTGCCTCGCTAACGAATCATGACTCCCTTAGAACGCTATCAAGCAGATCTGAAACGTCCCGActtct          <  5:134239/148‑1 (MQ=255)
                                                                                                                                          cgcGGGCACGCCCGCTCCCACAGACTGATTGCCATTAAGCGGGCCTATCGCTACAGTGACGGCCCGTTTTTTGTTGCCTCGCTAACGAATCATGACTCCCTTAGAACGCTATCAAGCAGATCTGAAACGTCCCGACTTCTTCCATGAc  <  7:161259/146‑1 (MQ=255)
                                                                                                                                           gcGGGCACGCCCGCTCCCACAGACTGATTGCCATTAAGCGGGCCTATCGCTACAGTGACGGCCCGTTTTTTGTTGCCTCGCTAACGAATCATGACTCCCTTAGAACGCTATCAAGCAGATCTGAAACGTcc                  >  6:351441/2‑131 (MQ=255)
                                                                                                                                           gcGGGCACGCCCGCTCCCACAGACTGATTGCCATTAAGCGGGCCTATCGCTACAGTGACGGCCCGTTTTTTGTTGCCTCGCTAACGAATCATGACTCCCTTAGAACGCTATCAAGCAGATCTGAAACGTcc                  <  5:351441/130‑1 (MQ=255)
                                                                                                                                           gcGGGCACGCCCGCTCCCACAGACTGATTGCCATTAAGCGGGCCTATCGCTACAGTGACGGCCCGTTTTTTGTTGCCTCGCTAACGAATCATGACTCCCTTAGAACGCTATCAAGCAGATCTGAAACGTCACGACTTCTTCCATg    <  5:373746/144‑1 (MQ=255)
                                                                                                                                           gcGGGCACGCCCGCTCCCACAGACTGATTGCCATTAAGCGGGCCTATCGCTACAGTGACGGCCCGTTTTTTGTTGCCTCGCTAACGAATCATGACTCCCATAGAACGCTATCAAGCAGATCTGAAACGTCCCGACTTCTTCCATg    >  6:373746/2‑145 (MQ=255)
                                                                                                                                                     |                                                                                                                                           
GCTGAAGCCCGTGATGCGGCCATTCAGAGCCTGAACGAAGCCTTGGCACCACAGCAGGACTGAGGCTGCCAAACGGTGGGAGCGGGCGTGCCCGCGAAACAGGCGCCGCGGTGCATGGCCACGGCTTTGCCGGTGTTTCGGGCACGCCGCTCCCACAGACTGATTGCCATTAAGCGGGCCTATCGCTACAGTGACGGCCCGTTTTTTGTTGCCTCGCTAACGAATCATGACTCCCTTAGAACGCTATCAAGCAGATCTGAAACGTCCCGACTTCTTCCATGAC  >  NC_002947/1499361‑1499643

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: