Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 524,336 A→G 8.0% intergenic (+203/‑613) 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,3360AG 8.0% 135.9 / 2.5 50intergenic (+203/‑613)tyrS/PP_16SCtyrosine‑‑tRNA ligase/16S ribosomal RNA
Reads supporting (aligned to +/- strand):  ref base A (24/22);  new base G (2/2);  total (26/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 = 8.43e-01

TTGGTGGTTGCTGTTGTAGGGGGAGGGGGCTGTTTTGAGAGGGTGTTGGCGGTGAGATGTGTCCGCCACAACCCTTGCAGCGCCTGTGAGATCGAGCGCCGCCCGCGTCCGGCTTGAGACATTCACCGCCTATGAAACCGAGCGCCGCGCGGGCGGCGCTCGGTCTCATAGGCGCTGAAACTGTCGTGGCGAACACTTGGCAGCCTCGATACGGTCTTTTCCCCTTATATTTCACTTTCTTGAAATTAAAGGTTGACGCCCTCTCGAATCCCCTTATAATG  >  NC_002947/524201‑524481
                                                                                                                                       |                                                                                                                                                 
ttGGTGGTTGCTGTGGTAGGGGGAGGGGGCTGTTTTGAGAGGGTGTTGGCGGTGAGATGTGTCCGCCACAACCCTTGCAGCGCCTGTGAGATCGAGCGCCGCCCGCGTCCGGCTTGAGACATTCACCGCCTATGAAACCGAGcgccgcc                                                                                                                                      >  7:274465/1‑148 (MQ=255)
   gtggtTGCTGTTGTAGGGGGAGGGGGCTGTTTTGAGAGGGTGTTGGCGGTGAGATGTGTCCGCCACAACCCTTGCAGCGCCTGTGAGATCGAGCGCCGCCCGCGTCCGGCTTGAGACATTCACCGCCTATGAGACCGAGCGCcgcccgc                                                                                                                                   <  1:182868/149‑2 (MQ=255)
   gtggtTGCTGTTGTAGGGGGAGGGGGCTGTTTTGAGAGGGTGTTGGCGGTGAGATGTGTCCGCCACAACCCTTGCAGCGCCTGTGAGATCGAGCGCCGCCCGCGTCCGGCTTGAGACATTCACCGCCTATGAAACCGAGCGCCGCgcgg                                                                                                                                   >  2:108517/1‑149 (MQ=255)
           tgttgtAGGGGGAGGGGGCTGTTTTGAGAGGGTGTTGGCGGTGAGATGTGTCCGCCACAACCCTTGCAGCGCCTGTGAGATCGAGCGCCGCCCGCGTCCGGCTTGAGACATTCACCGCCTATGAAACCGAGCGCCGCGCGGGCGGCGCt                                                                                                                           >  7:11031/1‑149 (MQ=255)
            gttgtCGGGGGAGGGGGCTGTTTTGAGAGGGTGTTGGCGGTGAGATGTGTCCGCCACAACCCTTGCAGCGCCTGTGAGATCGAGCGCCGCCCGCGTCCGGCTTGAGACATTCACCGCCTATGAAACCGAGCGCCGCGCGGGCGGCGCTc                                                                                                                          <  7:261624/149‑1 (MQ=255)
             ttgtAGGGGGAGGGGGCTGTTTTGAGAGGGTGTTGGCGGTGAGATGTGTCCGCCACAACCCTTGCAGCGCCTGTGAGATCGAGCGCCGCCCGCGTCCGGCTTGAGACATTCACCGCCTATGAAACCGAGCGCCGCGCGGGCGGCGCTCg                                                                                                                         <  6:150455/149‑1 (MQ=255)
             ttgtAGGGGGAGGGGGCTGTTTTGAGAGGGTGTTGGCGGTGAGATGTGTCCGCCACAACCCTTGCAGCGCCTGTGAGATCGAGCGCCGCCCGCGTCCGGCTTGAGACATTCACCGCCTATGAAACCGAGCGCCGCGCGGGCGGCGCTCg                                                                                                                         <  4:151547/149‑1 (MQ=255)
               gtAGGGGGAGGGGGCTGTTTTGAGAGGGTGTTGGCGGTGAGATGTGTCCGCCACAACCCTTGCAGCGCCTGTGAGATCGAGCGCCGCCCGCGTCCGGCTTGAGACATTCACCGCCTATGAAACCGAGCGCCGCGCGGGCGGCGCTCgg                                                                                                                        <  2:291965/148‑1 (MQ=255)
                  gggggagggggCTGTTTTGAGAGGGTGTTGGCGGTGAGATGTGTCCGCCACAACCCTTGCAGCGCCTGTGAGATCGAGCGCCGCCCGCGTCCGGCTTGAGACATTCACCGCCTATGAAACCGAGCGCCGCGCGGGCGGCGCTCGGtct                                                                                                                     >  1:57018/1‑148 (MQ=255)
                   gggtagggggCTGGTTTGTGCTGGTGTGGGCGGTGAGATGTGTCCGCCACATCCCTTGCAGCGCCTGTGAGATCGAGCGCCGCGCGCGTCCGGCTTGAGACATTCCCCGCCTATGAAACCGAGCGCCGCGCGGGCGGCGCTCGGTCTCa                                                                                                                   <  7:134360/149‑1 (MQ=37)
                     ggagggggCTGTTTTGAGAGGGTGTTGGCGGTGAGATGTGTCCGCCACAACCCTTGCAGCGCCTGTGAGATCGAGCGCCGCCCGCGTCCGGCTTGAGACATTCACCGCCTATGAGACCGAGCGCCGCCCGCGCGGCGCTCGGTTTCATa                                                                                                                 <  6:114946/149‑1 (MQ=255)
                       agggggCTGTTTTGAGAGGGTGTTGGCGGTGAGATGTGTCCGCCACAACCCTTGCAGCGCCTGTGAGATCGAGCGCCGCCCGCGTCCGGCTTGAGACATTCACCGCCTATGAAACCGAGCGCCGCGCGGGCGGCGCTCGGTCTCATAgg                                                                                                               <  7:25339/149‑1 (MQ=255)
                       agggggCTGTTTTGAGAGGGTGTTGGCGGTGAGATGTGTCCGCCACAACCCTTGCAGCGCCTGTGAGATCGAGCGCCGCCCGCGTCCGGCTTGAGACATTCACCGCCTATGAAACCGAGCGCCGCGCGGGCGGCGCTCGGTCTCATAgg                                                                                                               >  8:248146/1‑149 (MQ=255)
                            gCTGTTTTGAGAGGGTGTTGGCGGTGAGATGTGTCCGCCACAACCCTTGCAGCGCCTGTGAGATCGAGCGCCGCCCGCGTCCGGCTTGAGACATTCACCGCCTATGAAACCGAGCGCCGCGCGGGCGGCGCTCGGTCTCATAGGCGCt                                                                                                           >  3:60735/1‑148 (MQ=255)
                              tGTTTTGAGAGGGTGTTGGCGGTGAGATGTGTCCGCCACAACCCTTGCAGCGCCTGTGAGATCGAGCGCCGCCCGCGTCCGGCTTGAGACATTCACCGCCTATGAAACCGAGCGCCGCGCGGGCGGCGCTCGGTCTCATAGGCGCTGa                                                                                                         >  4:21420/1‑148 (MQ=255)
                                       agGGTGTTGGCGGTGAGATGTGTCCGCCACAACCCTTGCAGCGCCTGTGAGATCGAGCGCCGCCCGCGTCCGGCTTGAGACATTCACCGCCTATGAAACCGAGCGCCGCGCGGGCGGCGCTCGGTCTCATAGGCGCTGAAACTGTCGTg                                                                                               <  2:57018/149‑1 (MQ=255)
                                        gggTGTTGGCGGTGAGATGTGTCCGCCACAACCCTTGCAGCGCCTGTGAGATCGAGCGCCGCCCGCGTCCGGCTTGAGACATTCACCGCCTATGAAACCGAGCGCCGCGCGGGCGGCGCTCGGTCTCATAGGCGCTGAAACTGTCGTgg                                                                                              <  5:632/149‑1 (MQ=255)
                                         ggTGTTGGCGGTGAGATGTGTCCGCCACAACCCTTGCAGCGCCTGTGAGATCGAGCGCCGCCCGCGTCCGGCTTGAGACATTCACCGCCTATGAAACCGAgcg                                                                                                                                           >  2:286030/1‑103 (MQ=255)
                                         ggTGTTGGCGGTGAGATGTGTCCGCCACAACCCTTGCAGCGCCTGTGAGATCGAGCGCCGCCCGCGTCCGGCTTGAGACATTCACCGCCTATGAAACCGAgcg                                                                                                                                           <  1:286030/103‑1 (MQ=255)
                                             ttGGCGGTGAGATGTGTCCGCCACAACCCTTGCAGCGCCTGTGAGATCGAGCGCCGCCCGCGTCCGGCTTGAGACATTCACCGCCTATGAAACCGAGCGCCGCGCGGGCGGCGCTCGGTCTCATAGGCGCTGAAACTGTCGTGGCGaa                                                                                          >  6:255953/1‑148 (MQ=255)
                                                  ggTGAGATGTGTCCGCCACAACCCTTGCAGCGCCTGTGAGATCGAGCGCCGCCCGCGTCCGGCTTGAGACATTCACCGCCTATGAAACCGAGCGCCGCGCGGGCGGCGCTCGGTCTCATAGGCGCTGAAACTGTCGTGGCGAACACTTg                                                                                    <  1:167748/149‑1 (MQ=255)
                                                   gTGAGATGTGTCCGCCACAACCCTTGCAGCGCCTGTGAGATCGAGCGCCGCCCGCGTCCGGCTTGAGACATTCACCGCCTATGAAACCGAGCGCCGCGCGGGCGGCGCTCGGTCTCATAGGCGCTGAAACTGTCGTGGCGAACACTTgg                                                                                   <  1:108517/149‑1 (MQ=255)
                                                   gTGAGATGTGTCCGCCACAACCCTTGCAGCGCCTGTGAGATCGAGCGCCGCCCGCGTCCGGCTTGAGACATTCACCGCCTATGAAACCGAGCGCCGCGCGGGCGGCGCTCGGTCTCATAGGCGCTGAAACTGTCGTGGCAAACACTTgg                                                                                   <  4:60735/149‑1 (MQ=255)
                                                        aTGTGTCCGCCACAACCCTTGCAGCGCCTGTGAGATCGAGCGCCGCCCGCGTCCGGCTTGAGACATTCACCGCCTATGAAACCGAGCGCCGCGCGGGCGGCGCGCGGTCTCATCGGCGCGGAAACTGTCGTGGCGAACAGTTGgccgcc                                                                              >  8:302418/1‑149 (MQ=37)
                                                           tgtCCGCCACAACCCTTGCAGCGCCTGTGAGATCGAGCGCCGCCCGCGTCCGGCTTGAGACATTCACCGCCTATGAAACCGAGCGCCGCGCGGGCGGCGCTCGGTCTCATAGGCGCTGAAACTGTCGTGGCGAACACTTGGCAGCCTc                                                                            <  8:155995/148‑1 (MQ=255)
                                                                   acaACCCTTGCAGCGCCTGTGAGATCGAGCGCCGCCCGCGTCCGGCTTGAGACATTCACCGCCTATGAAACCGAGCGCCGCGCGGGCGGCGCTCGGTCTCATAGGCGCTGAAACTGTCGTGGCGAACACTTGGCAGCCTCGATACGGTc                                                                   >  8:126076/1‑149 (MQ=255)
                                                                          ttGCAGCGCCTGTGAGATCGAGCGCCGCCCGCGTCCGGCTTGAGACATTCACCGCCTATGAAACCGAGCGCCGCGCGGGCGGCGCTCGGTCTCATAGGCGCTGAAACTGTCGTGGCGAACACTTGGCAGCCTCGATACGGTCTTTTccc                                                            >  7:294403/1‑149 (MQ=255)
                                                                            gccgcgcCTGTTATATCGAGCGCCGCCCGCGTCCGGCTTGAGACATTCACCGCCTATGAACCCGCGCGCCGCGCGGGCGGCGCTCGGTCTCATAGGCGCTGAAACTGTCGTGGCGAACACTTGGCAGCCTCGATACGGTCTTTTccc                                                            <  8:294403/144‑1 (MQ=25)
                                                                            gCAGCGCCTGTGAGATCGAGCGCCGCCCGCGTCCGGCTTGAGACATTCACCGCCTATGAAACCGAGCGCCGCGCGGGCGGCGCTCGGTCTCATAGGCGCTGAAACTGTCGTGGCGAACACTTGGCAGCCTCGATACGGTCTTTTcccc                                                           <  7:126076/148‑1 (MQ=255)
                                                                            gCAGCGCCTGTGAGATCGAGCGCCGCCCGCGTCCGGCTTGAGACATTCACCGCCTATGAAACCGAGCGCCGCGCGGGCGGCGCTCGGTCTCATAGGCGCTGAAACTGTCGTGGCGAACACTTGGCAGCCTCGATACGGTCTTTTCCCCt                                                          >  3:102597/1‑149 (MQ=255)
                                                                               gcgcCTGTGAGATCGAGCGCCGCCCGCGTCCGGCTTGAGACATTCACCGCCTATGAAACCGAGCGCCGCGCGGGCGGCGCTCGGTCTCATAGGCGCTGAAACTGTCGTGGCGAACACTTGGCAGCCTCGATACGGTCTTTTCCCCTta                                                        >  2:268356/1‑148 (MQ=255)
                                                                                   cTGTGAGATCGAGCGCCGCCCGCGTCCGGCTTGAGACATTCACCGCCTATGAAACCGAGCGCCGCGCGGGCGGCGCTCGGTCTCATAGGCGCTGAAACTGTCGTGGCGAACACTTGGCAGCCTCGATACGGTCTTTTCCCCTTATAtt                                                    >  2:90355/1‑148 (MQ=255)
                                                                                     gtgAGATCGAGCGCCGCCCGCGTCCGGCTTGAGACATTCACCGCCTATGAAACCGCGCGCCGCGCGGGCGGCGCTCGGTCTCATAGGCGCTGAAACTGTCGTGGCGAACACTTGGCAGCCTCGATACGGTCTTTTCCCCTTATATTTCa                                                 <  6:236872/149‑1 (MQ=255)
                                                                                      tgAGATCGAGCGCCGCCCGCGTCCGGCTTGAGACATTCACCGCCTATGAAACCGAGCGCCGCGCGGGCGGCGCTCGGTCTCATAGGCGCTGAAACTGTCGTGGCGAACACTTGGCAGCCTCGATACGGTCTTTTCCCCTTATATTTCAc                                                >  6:192709/1‑149 (MQ=255)
                                                                                      tgAGATCGAGCGCCGCCCGCGTCCGGCTTGAGACATTCACCGCCTATGAAACCGAGCGCCGCGCGGGCGGCGCTCGGTCTCATAGGCGCTGAAACTGCCGTGGCGAACACTTGGCAGCCTCGATACGGTCTTTTCCCCTTATATTTc                                                  <  1:76790/147‑1 (MQ=255)
                                                                                        agaTCGAGCGCAGCCCGCGTCCGGCTTGAGACATTCACCGCCTATGAAACCGAGCGCCGCGCGGGCGGCGCTCGGTCTCATAGGCGCTGAAACTGTCGTGGCGAACACTTGGCAGCCTCGATACGGTCTTTTCCCCTTATATTTCACtt                                              >  8:238009/1‑149 (MQ=255)
                                                                                         gaTCGAGCGCCGCCCGCGTCCGGCTTGAGACATTCACCGCCTATGAAACCGAGCGCCGCGCGGGCGGCGCTCGGTCTCATAGGCGCTGAAACTCTCGTGGCGAACACTTGGCAGCCTCGCTACGGTCTTTTCCCCTTATATTTCACttt                                             >  5:51/1‑149 (MQ=255)
                                                                                              aGCGCCGCCCGCGTCCGGCTTGAGACATTCACCGCCTATGAAACCGAGCGCCGCGCGGGCGGCGCTCGGTCTCATAGGCGCTGAAACTGTCGTGGCGAACACTTGGCAGCCTCGATACGGTCTTTTCCCCTTATATTTCACTTTCTTGa                                        <  8:11031/149‑1 (MQ=255)
                                                                                                    gcacgcgTCCGGCTTGAGACATTCACCGCATATGAAACCGATCGCCGCGCGGGCGGCGCTCGGTATCATAGGCGCTGAAAATGTCGTGGCGAACACTTGGCAGCCTGGATACGGTCTTTTCCCCTTATATTTCACTTTCTTGAAATTaa                                  <  7:302418/146‑1 (MQ=18)
                                                                                                     cccgcGTCCGGCTTGAGACATTCACCGCCTATGAAACCGAGCGCCGCGCGGGCGGCGCTCGGTCTCATAGGCGCTGAAACTGTCGTGGCGAACACTTGGCAGCCTCGATACGGTCTTTTCCCCTTATATTTCACTTTCTTGAAATTaaa                                 <  7:248146/149‑1 (MQ=255)
                                                                                                       cgcgTCCGGCTTGAGACATTCACCGCCTATGAAACCGAGCGCCGCGCGGGCGGCGCTCGGTCTCATAGGCGCTGAAACTGTCGTGGCGAACACTTGGCAGCCTCGATACGGTCTTTTCCCCTTATATTTCACTTTCTTGAAATTAAAgg                               <  7:238009/149‑1 (MQ=255)
                                                                                                                        aTTCACCGCCTATGAAACCGAGCGCCGCGCGGGCGGCGCTCGGTCTCATAGGCGCTGAAACTGTCGTGGCGAACACTTGGCAGCCTCGATACGGGCTGTTCCCCTTATATTTCACTTTCTTGAAATTAAAGGATGACGCCCTCTCGAAt              >  6:149190/1‑149 (MQ=255)
                                                                                                                          tCACCGCCTATGAAACCGAGCGCCGCGCGGGCGGCGCTCGGTCTCATAGGCGCTGAAACTGTCGTGGCGAACACTTGGCAGCCTCGATACGGTCTTTTCCCCTTATATTTCACTTTCTTGAAATTAAAGGTTGACGCGCTCTCGAATcc            >  4:270827/1‑149 (MQ=255)
                                                                                                                             ccGCCTATGAAACCGAGCGCCGCGCGGGCGGCGCTCGGTCTCATAGGCGCTGAAACTGTCGTGGCGAACACTTGGCAGCCTCGATACGGTCTTTTCCCCTTATATTTCACTTTCTTGAAATTAAAGGTTGACGCCCTCTCGAATCCCCt         >  8:31245/1‑149 (MQ=255)
                                                                                                                              cGCCTATGAGACCGAGCGCCGCCCGCGCGGCGCTCGGTTTCATAGGCGGTGAAACTGTCGTGGAGAACACTTGGCAGCCTCGATACGGTCTTTTCCCCTTATATTTCACTTTCTTGAAATTAAAGGTTGACGCCCTCTCGAATCCACt         >  5:126118/1‑148 (MQ=11)
                                                                                                                                ccTATGAAACCGAGCGCCGCGCGGGCGGCGCTCGGTCTCATAGGCGCTGAAACTGTCGTGGCGAACACTTGGCAGCCTCGATACGGTCTTTTcccc                                                           >  2:170953/1‑96 (MQ=255)
                                                                                                                                ccTATGAAACCGAGCGCCGCGCGGGCGGCGCTCGGTCTCATAGGCGCTGAAACTGTCGTGGCGAACACTTGGCAGCCTCGATACGGTCTTTTcccc                                                           <  1:170953/96‑1 (MQ=255)
                                                                                                                                ccTATGAAACCGAGCGCCGCGCGGGCGGCGCTCGGTCTCATAGGCGCTGAAACTGTCGTGGCGAACACTTGGCAGCCTCGATACGGTCTTTTCCCCTTATATTTCACTTTCTTGAAATTAAAGGTTGACGCCCTCTCGAATCCCCTta       >  6:291783/1‑148 (MQ=255)
                                                                                                                                 cTATGAGACCGAGCGCCGCCCGCGCGGCGCTCGGTTTCATAGGCGCTGAAACTGTCGTGGCGAACACTTGGCAGCCTCGATACGGTCTTTTCACCTTATATTTCACTTTCTTGAAATTAAAGGTTGACGCCCTCTCGAATCCCCTtat      >  2:161526/1‑148 (MQ=11)
                                                                                                                                    tGAAACCGAGCGCCGCCCGGGCGGCGCTCGGTCTCATAGGCGCTGAAACAGACGTGGCGAACACATGGCAGCCTCGATACGGTCTTTTCCCCTTATAATTCACTATCTTCAAATTAAAGGTTGACGCCATCACGAATCCACTAATAATg  >  8:195813/1‑149 (MQ=14)
                                                                                                                                       |                                                                                                                                                 
TTGGTGGTTGCTGTTGTAGGGGGAGGGGGCTGTTTTGAGAGGGTGTTGGCGGTGAGATGTGTCCGCCACAACCCTTGCAGCGCCTGTGAGATCGAGCGCCGCCCGCGTCCGGCTTGAGACATTCACCGCCTATGAAACCGAGCGCCGCGCGGGCGGCGCTCGGTCTCATAGGCGCTGAAACTGTCGTGGCGAACACTTGGCAGCCTCGATACGGTCTTTTCCCCTTATATTTCACTTTCTTGAAATTAAAGGTTGACGCCCTCTCGAATCCCCTTATAATG  >  NC_002947/524201‑524481

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