Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 524,606 G→A 9.3% intergenic (+473/‑343) tyrS → / → PP_16SC tyrosine‑‑tRNA ligase/16S ribosomal RNA

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_002947524,6060GA 9.3% 105.0 / 2.9 44intergenic (+473/‑343)tyrS/PP_16SCtyrosine‑‑tRNA ligase/16S ribosomal RNA
Reads supporting (aligned to +/- strand):  ref base G (17/22);  new base A (2/2);  total (20/24)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 7.37e-01

CCCCACTTCCAGCGACATCGGAACGAAAAACTCCTTGAGATTCAATGAGTTAAGTAGTTAAGGCAAAAGCTGAAAGGGCTACGATCGAATGATCGAAAGCGGTTGAGATGATGGTTGACAGCGCTTCCAAGTGCTGTATGATTCGCCTCCCGCTACGAGAGATCGCAGCGAGCCAAGTGTTTGAAGCTAAACGAGTTTCTCGCAGAAAACTTCAAAATAAACGCTTGACAGGCCCTGAGGAAAGCGTAGAATGCGCGCCTCGGTTG  >  NC_002947/524484‑524749
                                                                                                                          |                                                                                                                                               
ccccACTTCCAGCGACATCGGAACGAAAAACTCCTTGAGATTCAATGAGTTAAGTAGTTAAGGCAAAAGCTGAAAGGGCTACGATCGAATGATCGAAAGCGGTTGAGATGATGGTTGACAGCGCTTCCAAGTGCTGTATGATTCGCCTc                                                                                                                       <  8:150209/149‑1 (MQ=32)
   cACTTCCAGCGACATCGGAACGAAAAACTCCTTGAGATTCAATGAGTTAAGTAGTTAAGGCAAAAGCTGAAAGGGCTACGATCGAATGATCGAAAGCGGTTGAGATGATGGTTGACAGCGCTTCCAAGTGCTGTATGATTCGCCTCCCg                                                                                                                    >  5:298111/1‑149 (MQ=32)
         cAGCGACATCGGAACGAAAAACTCCTTGAGATTCAATGAGTTAAGTAGTTAAGGCAAAAGCTGAAAGGGCTACGATCGAATGATCGAAAGCGGTTGAGATGATGGTTGACAGCGCTTCCAAGTGCTGTGTGATTCGCCTCCCGCTACga                                                                                                              <  4:279558/149‑1 (MQ=17)
         cAGCGACATCGGAACGAAAAACTCCTTGAGATTCAATGAGTTAAGTAGTTAAGGCAAAAGCTGAAAGGGCTACGATCGAATGATCGAAAGCGGTTGAGATGATGGTTGACAGCGCTTCCAAGTGCTGTATGATTCGCCTCCCGCTACga                                                                                                              <  4:51975/149‑1 (MQ=32)
            atctaTCGGAAGGAAAAACTCCTTGAGATTCTATGTGTTAAGTAGTTAAGGCAATAGCTGAAAGGTCTTCGATCGAATGATCGAATGCGGTTGAGATGATGGTTGACAGCGCTTCCAAGTGCTGTATGATTCGCCTCCCGCTACgagag                                                                                                           <  7:199322/145‑1 (MQ=11)
              aCATCGGAACGAAAAACTCCTTGAGATTCAATGAGTTAAGTAGTTAAGGCAAAAGCTGAAAGGGCTACGATCGAATGATCGAAAGCGGTTGAGATGATGGTTGACAGCGCTTCCAAGTGCTGTATGATTCGCCTCCCGCTACGAGAGAt                                                                                                         <  5:12407/149‑1 (MQ=32)
                aTCGGAACGAAAAACTCCTTGAGATTCAATGAGTTAAGTAGTTAAGGCAAAAGCTGAAAGGGCTACGATCGAATGATCGAAAGCGGTTGAGATGATGGTTGACAGCGCTTCCAAGTGCTGTATGATTCGCCTCCCGCTACGAGAGATCg                                                                                                       <  1:69164/149‑1 (MQ=32)
                 tCGGAACGAAAAACTCCTAGAGATTCAATGAGTTAAGTAGTTAAGGCAAAAGCTGAAAGGGCTACGATCGAATGATCGAAAGCGCTTGAGATGATGGATGACAGCCATTCCAAGTGCTGTATGATTAGCCTCCAGCTACgagac                                                                                                           >  8:199322/1‑143 (MQ=11)
                    gAACGAAAAACTCCTTGAGATTCAATGAGTTAAGTAGTTAAGGCAAAAGCTGAAAGGGCTACGATCGAATGATCGAAAGCGGTTGAGATGATGGTTGACAGCGCTTCCAAGTGCTGTATGATTCGCCTCCCGCTACGAGAGATCGCa                                                                                                     >  8:252740/1‑147 (MQ=32)
                    gAACGAAAAACTCCTTGAGATTCAATGAGTTAAGTAGTTAAGGCAAAAGCTGAAAGGGCTACGATCGAATGATCGAAAGCGGTTGAGATGATGGTTGACAGCGCTTCCAAGTGCTGTATGATTCGCCTCCCGCTACGAGAGATCGCAGc                                                                                                   >  2:94865/1‑149 (MQ=32)
                                    gagaTTCAATGAGTTAAGTAGTTAAGGCAAAAGCTGAAAGGGATACGATCGAATGATCGAAAGCGGTTGAGATGATGGTTGACAGCGCTTCCAAGTGCTGTATGATTCGCCTCCCGCTACGAGAGATCGCAGCGAGCCAAGTGTTTg                                                                                     >  5:23601/1‑147 (MQ=17)
                                     agaTTCAATGAGTTAAGTAGTTAATGCAAAAGCTGAAAGGGCTACGATCGAATGATCGAAAGCGGTTGAGATGATGGTTGACAGCGCTTCCAAGTGCTGTATGATTCGCCTCCCGCTACGAGAGATCGCAGCGAGCCAAGTGTTTGaa                                                                                   <  3:59493/148‑1 (MQ=17)
                                      gaTTCAATGAGTTAAGTAGTTAAGGCAAAAGCTGAAAGGGCTACGATCGAATGATCGAAAGCGGTTGAGATGATGGTTGACAGCGCTTCCAAGTGCTGTATGATTCGCCTCCCGCTACGAGAGATCGCAGCGAGCCaa                                                                                            <  6:4050/138‑1 (MQ=32)
                                      gaTTCAATGAGTTAAGTAGTTAAGGCAAAAGCTGAAAGGGCTACGATCGAATGATCGAAAGCGGTTGAGATGATGGTTGACAGCGCTTCCAAGTGCTGTATGATTCGCCTCCCGCTACGAGAGATCGCAGCGAGCCaa                                                                                            >  5:4050/1‑138 (MQ=32)
                                           aaTGAGTTAAGTAGTTAAGGCAAAAGCTGAAAGGGCTACGATCGAATGATCGAAAGCGGTTGAGATGATGGTTGACAGCACTTCCAAGTGCTGTATGATTCGCCTCCCGCTACGAGAGATCGCAGCGAGCCAAGTGTTTGAAGCTAAAc                                                                            >  1:89352/1‑149 (MQ=17)
                                               aGTTAAGTAGTTAAGGCAAAAGCTGAAAGGGCTACGATCGAATGATCGAAAGCGGTTGAGATGATGGTTGACAGCGCTTCCAAGTGCTGTATGATTCGCCTCCCGCTACGAGAGATCGCAGCGAGCCAAGTGTTTGAAGCTAAACGAg                                                                         <  2:261732/148‑1 (MQ=32)
                                                   aaGTAGTTAAGGCAAAAGCTGAAAGGGCTACGATCGAATGATCGAAAGCGGTTGAGATGATGGTTGACAGCGCTTCCAAGTGCTGTATGATTCGCCTCCCGCTACGAGAGATCGCAGCGAGCCAAGTGTTTGAAGCTAAACGAGTTtct                                                                    >  7:22525/1‑149 (MQ=32)
                                                        gtTAAGGCAAAAGCTGAAAGGGCTACGATCGAATGATCGAAAGCGGTTGAGATGATGGTTGACAGCACTTCCAAGTGCTGTATGATTCGCCTCCCGCTACGAGAGATCGCAGCGAGCCAAGTGTTTGAAGCTAAACGAGTTTCTCGCa                                                                <  2:89352/148‑1 (MQ=17)
                                                         ttAAGGCAAAAGCTGAAAGGGCTACGATCGAATGATCGAAAGCGGTTGAGATGATGGTTGACAGCGCTTCCAAGTGCTGTATGATTCGCCTCCCGCTACGAGAGATCGCAGCGAGCCAAGTGTTTGAAGCTAAACGAGTTTCTCGCAg                                                               >  8:22543/1‑148 (MQ=32)
                                                            acgCAAAAGCTGAACGGGCTACGATCGAATGATCGAAAGCGGTTGAGATGATGGTTGACAGCGCTTCCAAGTGCTGTATGATTCGCCTCCCGCTACTAGAGATCGCAGCGAGCCAAGTGTTTGAAGCTAAACGAGTTTCTCGCAGaaaa                                                           <  6:105659/147‑1 (MQ=12)
                                                                 aaaGCTGAAAGGGCTACGATCGAATGATCGAAAGCGGTTGAGATGATGGTTGACAGCGCTTCCAAGTGCTGTATGATTCGCCTCCCGCTACGAGAGATCGCAGCGAGCCAAGTGTTTGAAGCTAAACGAGTTTCTCGCAGAAAACTTCa                                                      <  6:298111/149‑1 (MQ=34)
                                                                             gCTACGATCGAATGATCGAAAGCGGTTGAGATGATGGTTGACAGCGCTTCCAAGTGCTGTATGATTCGCCTCCCGCTACGAGAGATCGCAGCGAGCCAAGTGTTTGAAGCTAAACGAGTTTCTCGCAGAAAACTTCAAAATAAACGCtt                                          <  4:180856/149‑1 (MQ=34)
                                                                                    tCGAATGATCGAAAGCGGTTGAGATGATGGTTGACAGCGCTTCCAAGTGCTGTATGATTCGCCTCACGCTACGAGAGATCGCAGCGAGCCAAGTGTTTGAAGCTAAACGAGTTTCTCTCAGAAAACTTCAAAATAAACGCTTGACAGGc                                   <  7:252740/149‑1 (MQ=21)
                                                                                       aaTGATCGAAAGCGGTTGAGATGATGGTTGACAGCGCTTCCAAGTGCTGTATGATTCGCCTCCCGCTACGAGAGATCGCAGCGAGCCAAGTGTTTGAAGCTAAACGAGTTTCTCGCAGAAAACTTCAAAATAAACGCTTGACAGGCCCt                                <  4:50960/149‑1 (MQ=37)
                                                                                           aTCGAAAGCGGTTGAGATGATGGTTGACAGCGCTTCCAAGTGCTGTATGATTCGCCTCCCGCTACGAGAGATCGCAGCGAGCCAAGTGTTTGAAGCTAAACGAGTTTCTCGCAGAAAACTTCAAAATAAACGCTTGACAGGCCCTGAg                             <  7:162097/148‑1 (MQ=37)
                                                                                            tCGAAAGCGGTTGAGATGATGGTTGACAGCGCTTCCAAGTGCTGTATGATTCGCCTCCCGCTACGAGAGATCGCAGCGAGCCAAGTGTTTGAAGCTAAACGAGTTTCTCGCAGAAAACTTCAAAATAAACGCTTGACAGGCCCTGAgg                            >  2:36818/1‑148 (MQ=37)
                                                                                               aaaGCGGTTGAGATGATGGTTGACAGCGCTTCCAAGTGCTGTATGATTCGCCTCCCGCTACGAGAGATCGCAGCGAGCCAAGTGTTTGAAGCTAAACGAGTTTCTCGCAGAAAACTTCAAAATAAACGCTTGACAGGCCCTGAGGAAAg                        >  2:48273/1‑149 (MQ=35)
                                                                                                  gCGGTTGAGATGATGGTTGACAGCGCTTCCAAGTGCTGTATGATTCGCCTCCCGCTACGAGAGATCGCAGCGAGCCAAGTGTTTGAAGCTAAACGAGTTTCTCGCAGAAAACTTCAAAATAAACGCTTGACAGGCCCTGAGGAAAGCGt                     <  4:165574/149‑1 (MQ=34)
                                                                                                     gTTGAGATGATGGTTGACAGCGCTTCCAAGTGCTGTATGATTCGCCTCCCGCTACGAGAGATCGCAGCGAGCCAAGTGTTTGAAGCTAAACGAGTTTCTCGCAGAAAACTTCAAAATAAACGCTTGACAGGCCCTGAGGAAAGCGTAGa                  <  7:22543/149‑1 (MQ=34)
                                                                                                     gTTGAGATGATGGTTGACAGAGCTTCCAAGTGCTGTATGATTCGCCTCCCGCTACGAGAGATCGCAGAGAGCCAAGTGTTTGAAGCTAAACGAGTTTATCGCAGAAAACTTAAAAATAAACTCTTGACAGGCCCTGAGTAAAGAGTAg                   >  8:226710/1‑148 (MQ=11)
                                                                                                         agaTGATGGTTGACAGCGCTTCCAAGTGCTGTATGATTCGCCTCCCGCTACGAGAGATCGCAGCGAGCCAAGTGTTTGAAGCTAAACGAGTTTCTCGCAGAAAACTTCAAAATAAACGCTTGACAGGCCCTGAGGAAAGCGTAGAATgc              >  5:66658/1‑149 (MQ=34)
                                                                                                            tgatgGTTGACAGCGCTTCCAAGTGCTGTATGATTCGCCTCCCGCTACGAGAGATGGCAGCGAGCCAAGTGTTTGAAGCTAAACGAGTTTCTCGCAGAAAACTTCAAAATAAACGCTTGACAGGCCCTGAGGAAAGCGTAGAATgcgc            <  5:257984/148‑1 (MQ=21)
                                                                                                             gatgGTTGACAGCGCTTCCAAGTGCTGTATGATTCGCCTCCCGCTACGAGAGATCGCAGCGAGCCAAGTGTTTGAAGCTAAACGAGTTTCTCGCAGAAAACTTCAAAATAAACGCTTGACAGGCCCTGAGGAAAGCGTa                    <  1:65797/139‑1 (MQ=34)
                                                                                                             gatgGTTGACAGCGCTTCCAAGTGCTGTATGATTCGCCTCCCGCTACGAGAGATCGCAGCGAGCCAAGTGTTTGAAGCTAAACGAGTTTCTCGCAGAAAACTTCAAAATAAACGCTTGACAGGCCCTGAGGAAAGCGTa                    >  2:65797/1‑139 (MQ=34)
                                                                                                                atcataCAGCACTTGGAAGCGCTGTATGATTCGCCTCCCGCTACGAGAGATCGCAGCGAGCCAAGTGTTTGAAGCTAAACGAGTTTCTCGCAGAAAACTTCAAAATAAACGCTTGACAGGCCCTGAGGAAAGCGTAGAATGCGCGCCTc       >  5:4296/6‑149 (MQ=11)
                                                                                                                   ataCAGCACTTGGAAGCGCTGTATGATTCGCCTCCCGCTACGAGAGATCGCAGCGAGCCAAGTGTTTGAAGCTAAACGAGTTTCTCGCAGAAAACTTCAAAATAAACGCTTGACAGGCCCTGAGGAAAGCGTAGAATGCGCGCCTCgg     <  6:4296/146‑1 (MQ=255)
                                                                                                                    gACAGCGCTTCCAAGTGCTGTATGATTCGCCTCCCGCTACGAGAGATCGCAGCGAGCCAAGTGTTTGAAGCTAAACGAGTTTCTCGCAGAAAACTTCAAAATAAACGCTTGACAGGCCCTGAGGAAAGCGTAGAATGCGCGCCTCGGtt   >  7:41009/1‑149 (MQ=34)
                                                                                                                     aCAGCGCTTCCAAGTGCTGTATGATTCGCCTCCCGCTACGAGAGATCGCAGCGAGCCAAGTGTTTGAAGCTAAACGAGTTTCTCGCAGAAAACTTCAAAATAAACGCTTGACAGGCCCTGAGGAAAGCGTAGaa                 >  6:162437/1‑134 (MQ=34)
                                                                                                                     aCAGCGCTTCCAAGTGCTGTATGATTCGCCTCCCGCTACGAGAGATCGCAGCGAGCCAAGTGTTTGAAGCTAAACGAGTTTCTCGCAGAAAACTTCAAAATAAACGCTTGACAGGCCCTGAGGAAAGCGTAGaa                 <  5:162437/134‑1 (MQ=34)
                                                                                                                     aCAGCGCTTCCAAGTGCTGTATGATTCGCCTCCCGCTACGAGAGATCGCAGCGAGCCAAGTGTTTGAAGCTAAACGAGTTTCTCGCAGAAAACTTCAAAATAAACGCTTGACAGGCCCTGAGGAAAGCGTAGAATGCGCGCCTCGGtt   >  6:218894/1‑148 (MQ=34)
                                                                                                                     aCAGCGCTTCCAAGTGCTGTATGATTCGCCTCCCGCTACGAGAGATCGCAGCGAGCCAAGTGTTTGAAGCTAAACGAGTTTCTCGCAGAAAACTTCAAAATAAACGCTTGACAGGCCCTGAGGAAAGCGTAGAATGCGCGCCTCGGtt   >  7:79841/1‑148 (MQ=34)
                                                                                                                     aCAGCGCTTCCAAGTGCTGTATGATTCGCCTCCCGCTACGAGAGATCGCAGCGAGCCAAGTGTTTGAAGCTAAACGAGTTTCTCGCAGAAAACTTCAAAATAAACGCTTGACAGGCCCTGAGGAAAGCGTAGAATGCGCGCCTCGGtt   <  5:218894/148‑1 (MQ=34)
                                                                                                                     aCAGCGCTTCCAAGTGCTGTATGATTCGCCTCCCGCTACGAGAGATCGCAGCGAGCCAAGTGTTTGAAGCTAAACGAGTTTCTCGCAGAAAACTTCAAAATAAACGCTTGACAGGCCCTGAGGAAAGCGTAGAATGCGCGCCTCGGTTg  <  6:66658/149‑1 (MQ=34)
                                                                                                                     aCAGCGCTTCCAAGTGCTGTATGATTCGCCTCCCGCTACGAGAGATCGCAGCGAGCCAAGTGTTTGAAGCTAAACGAGTTTCTCGCAGAAAACTTCAAAATAAAAGCTTGACAGGCCCTGAGGAAAGCGGAGAATGCGCGCCTCGGtt   >  3:50277/1‑148 (MQ=21)
                                                                                                                          |                                                                                                                                               
CCCCACTTCCAGCGACATCGGAACGAAAAACTCCTTGAGATTCAATGAGTTAAGTAGTTAAGGCAAAAGCTGAAAGGGCTACGATCGAATGATCGAAAGCGGTTGAGATGATGGTTGACAGCGCTTCCAAGTGCTGTATGATTCGCCTCCCGCTACGAGAGATCGCAGCGAGCCAAGTGTTTGAAGCTAAACGAGTTTCTCGCAGAAAACTTCAAAATAAACGCTTGACAGGCCCTGAGGAAAGCGTAGAATGCGCGCCTCGGTTG  >  NC_002947/524484‑524749

Alignment Legend
Aligned base mismatch/match (shaded by quality score): ATCG/ATCG < 3 ≤ ATCG/ATCG < 14 ≤ ATCG/ATCG < 21 ≤ ATCG/ATCG < 32 ≤ ATCG/ATCG < 36 ≤ ATCG/ATCG
Unaligned base: atcg    Masked matching base: atcg    Alignment gap:     Deleted base: