Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 524,365 C→T 8.6% intergenic (+232/‑584) 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,3650CT 8.6% 125.1 / 2.1 47intergenic (+232/‑584)tyrS/PP_16SCtyrosine‑‑tRNA ligase/16S ribosomal RNA
Reads supporting (aligned to +/- strand):  ref base C (23/18);  new base T (2/2);  total (25/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.07e-01

GGGGGAGGGGGCTGTTTTGAGAGGGTGTTGGCGGTGAGATGTGTCCGCCACAACCCTTGCAGCGCCTGTGAGATCGAGCGCCGCCCGCGTCCGGCTTGAGACATTCACCGCCTATGAAACCGAGCGCCGCGCGGGCGGCGCTCGGTCTCATAGGCGCTGAAACTGTCGTGGCGAACACTTGGCAGCCTCGATACGGTCTTTTCCCCTTATATTTCACTTTCTTGAAATTAAAGGTTGACGCCCTCTCGAATCCCCTTATAATGCGCCCCACTTCCAGCGACATCGGAACGAAA  >  NC_002947/524219‑524511
                                                                                                                                                  |                                                                                                                                                  
gggggagggggCTGTTTTGAGAGGGTGTTGGCGGTGAGATGTGTCCGCCACAACCCTTGCAGCGCCTGTGAGATCGAGCGCCGCCCGCGTCCGGCTTGAGACATTCACCGCCTATGAAACCGAGCGCCGCGCGGGCGGCGCTCGGtct                                                                                                                                                   >  1:57018/1‑148 (MQ=255)
 gggtagggggCTGGTTTGTGCTGGTGTGGGCGGTGAGATGTGTCCGCCACATCCCTTGCAGCGCCTGTGAGATCGAGCGCCGCGCGCGTCCGGCTTGAGACATTCCCCGCCTATGAAACCGAGCGCCGCGCGGGCGGCGCTCGGTCTCa                                                                                                                                                 <  7:134360/149‑1 (MQ=37)
   ggagggggCTGTTTTGAGAGGGTGTTGGCGGTGAGATGTGTCCGCCACAACCCTTGCAGCGCCTGTGAGATCGAGCGCCGCCCGCGTCCGGCTTGAGACATTCACCGCCTATGAGACCGAGCGCCGCCCGCGCGGCGCTCGGTTTCATa                                                                                                                                               <  6:114946/149‑1 (MQ=255)
     agggggCTGTTTTGAGAGGGTGTTGGCGGTGAGATGTGTCCGCCACAACCCTTGCAGCGCCTGTGAGATCGAGCGCCGCCCGCGTCCGGCTTGAGACATTCACCGCCTATGAAACCGAGCGCCGCGCGGGCGGCGCTCGGTCTCATAgg                                                                                                                                             >  8:248146/1‑149 (MQ=255)
     agggggCTGTTTTGAGAGGGTGTTGGCGGTGAGATGTGTCCGCCACAACCCTTGCAGCGCCTGTGAGATCGAGCGCCGCCCGCGTCCGGCTTGAGACATTCACCGCCTATGAAACCGAGCGCCGCGCGGGCGGCGCTCGGTCTCATAgg                                                                                                                                             <  7:25339/149‑1 (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)
                           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)
                                                                                                                        cGAGCGCCGCGCGGGCGGCGCTCGGTCTCATAGGCGCTGAAACTGTCGTGGCGAACACTTGGCAGCCTCGATACGGTCTTTTCCCCTTATATTTCACTTTCTTGAAATTAAAGGTTGACGCCCTCTCGAATCCCCTTATAATGCGcccc                          >  6:280471/1‑149 (MQ=255)
                                                                                                                        cGAGCGCCGCGCGGGCGGCGCTCGGTCTCATAGGCGCTGAAACTGTCGTGGCGAACACTTGGCAGCCTCGATACGGTCTTTTCCCCTTATATTTCACTTTCTTGAAATTAAAGGTTGACGCCCTCTCGAATCCCCTTATAATGCGcccc                          <  5:192709/149‑1 (MQ=255)
                                                                                                                          aGCGCCGCGCGGGCGGCGCTCGGTCTCATAGGCGCTGAAACTGTCGTGGCGAACACTTGGCAGCCTCGATACGGTCTTTTCCCCTTATATTTCACTTTCTTGAAATTAAAGGTTGACGCCCTCTCGAATCCCCTTATAATGCGCCCCAc                        <  1:90355/149‑1 (MQ=255)
                                                                                                                              ccgcGCGGGCGGCGCTCGGTCTCATAGGCGCTGAAACTGTCGTGGCGAACACTTGGCAGCCTCGATACGGTCTTTTCCCCTTATATTTCACTTTCTTGAAATTAAAGGTTGACGCCCTCTCGAATCCCCTTATAATGCGCCCCACTTcc                    >  4:278190/1‑149 (MQ=255)
                                                                                                                                 cccgcgcggcgCTCGGTTTCATAGGCGGTGAAACTGTCGTGGCGAACACTTGGCAGCCTCGATACGGTCTTTTCCCCTTATATTTCACTTTCTTGAAATTAAAGGTTGACGCCCTCTCGAATCCCCTTATAATGCGCCCCACTTCCAGc                 <  6:126118/144‑1 (MQ=255)
                                                                                                                                          cgcTCGGTCTCATAGGCGCTGAAACTGTCGTGGCGAACACTTGGCAGCCTCGATACGGTCTTTTCCCCTTATATTTCACTTTCTTGAAATTAAAGGTTGACGCCCTCTCGAATCCCCTTATAATGCGCCCCACTTCCAGCGACATCgg         >  1:10878/1‑148 (MQ=255)
                                                                                                                                           gcTCTGTATAATAGTCGATGAAACTGTCGTTGCGAACACTTGGCAGACTCGATACGGTCTTTTCCCCTTATATTTCACTTTCTTGAAATTAAAGGTTGACGCCCTCTCGAATCCCCTTATAATGCGCCCCACTTCCAGCGACATCGGaa       <  5:255953/149‑1 (MQ=14)
                                                                                                                                             tCGGTCTCATAGGCGCTGAAACTGTCGTGGCGAACACTTGGCAGCCTCGATACGGTCTTTTCCCCTTATATTTCACTTTCTTGAAATTAAAGGTTGACGCCCTCTCGAATCCCCTTATAATGCGCCCCACTTCCAGCGACATCGGAAc      >  3:49968/1‑148 (MQ=255)
                                                                                                                                                 tctcATAGGCGCTGAAACTGTCGTGGCGAACACTTGGCAGCCTCGATACGGTCTTTTCCCCTTATATTTCACTTTCTTGAAATTAAAGGTTGACGCCCTCTCGAATCCCCTTATAATGCGCCCCACTTCCAGCGACATCGGAACGaaa  >  6:149701/1‑148 (MQ=255)
                                                                                                                                                  |                                                                                                                                                  
GGGGGAGGGGGCTGTTTTGAGAGGGTGTTGGCGGTGAGATGTGTCCGCCACAACCCTTGCAGCGCCTGTGAGATCGAGCGCCGCCCGCGTCCGGCTTGAGACATTCACCGCCTATGAAACCGAGCGCCGCGCGGGCGGCGCTCGGTCTCATAGGCGCTGAAACTGTCGTGGCGAACACTTGGCAGCCTCGATACGGTCTTTTCCCCTTATATTTCACTTTCTTGAAATTAAAGGTTGACGCCCTCTCGAATCCCCTTATAATGCGCCCCACTTCCAGCGACATCGGAACGAAA  >  NC_002947/524219‑524511

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