Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 524,180 T→G 20.3% intergenic (+47/‑769) 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,1800TG20.3% 73.6 / 2.7 47intergenic (+47/‑769)tyrS/PP_16SCtyrosine‑‑tRNA ligase/16S ribosomal RNA
Reads supporting (aligned to +/- strand):  ref base T (8/25);  new base G (6/2);  total (20/27)
Fisher's exact test for biased strand distribution p-value = 1.20e-02
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

TGCAGTGGTCGACAAGGACTTCATGTTTGTCCTTGGTGCTACCCATGTGTGCCAAGCGGGCAAGAAGTCGTTTGGCCGGGTTACCCTCAAGGCTGAGTGATTGCCTGGATGGTGTAGCTATATAGAAGCGGGGAGGCCGGTAGGCCTCCCCGC‑‑TTTGTTTTTGGGGTTTGGTGGTTGCTGTTGTAGGGGGAGGGGGCTGTTTTGAGAGGGTGTTGGCGGTGAGATGTGTCCGCCACAACCCTTGCAGCGCCTGTGAGATCGAGCGCCGCCCGCGTCCGGCTTGAGACATTCACC  >  NC_002947/524034‑524327
                                                                                                                                                    |                                                                                                                                                         
tGCAGTGGTCGACAAGGACTTCATGTTTGTCCTTGGTGCTACCCATGTGTGCCAAGCGGGCAAGAAGTCGTTTGGCCGGGTTACCCTCAAGGCTGAGTGATTGCCTGGATGGTGTAGCTATATAGAAGCGGGGAGGCCGGTAGGCCTcc                                                                                                                                                         <  3:8523/149‑1 (MQ=255)
tGCAGTGGTCGACAAGGACTTCATGTTTGTCCTTGGTGCTACCCATGTGTGCCAAGCGGGCAAGAAGTCGTTTGGCCGGGTTACCCTCAAGGCTGAGTGATTGCCTGGATGGTGTAGCTATATAGAAGCGGGGAGGCCGGTAGGCCAcc                                                                                                                                                         >  2:5496/1‑149 (MQ=255)
              aGGACTTCATGTTTGTCCTTGGTGCTACCCATGTGTGCCAAGCGGGCAAGAAGTCGTTTGGCCGGGTTACCCTCAAGGCTGAGTGATTGCCTGGATGGTGTAGCTATATAGAAGCGGGGAGGCCGGTAGGCCGCCCCGC‑‑TTTGTTTTTg                                                                                                                                       <  3:228548/149‑1 (MQ=255)
                   ttCATGTTTGTCCTTGGTGCTACCCATGTGTGCCAAGCGGGCAAGAAGTCGTTTGGCCGGGTTACCCTCAAGGCTGAGTGATTGCCTGGATGGTGTAGCTATATAGAAGCGGGGGGGGCCGGTGGGCGCCCCCC‑‑CTTGTTTTTTGGGtt                                                                                                                                 >  1:129986/1‑149 (MQ=255)
                   ttCATGTTTGTCCTAGGTGCTACCCATGTGTGCCAACCGGGCAAGAAGTCGTTAGGCCGGGTTACCCTCAAGGCTGAGTGATTCCCTGGATGGTGTAGCTATATAGAAGCGGGGATGGCCGTAGGCCAACCCCCT‑‑TTTGTTTTTggctt                                                                                                                                   >  7:79262/1‑146 (MQ=255)
                     cATGTTTGTCCTTGGTGCTACCCATGTGTGCCAAGCGGGCAAGAAGTCGTTTGGCCGGGTTACCCTCAAGGCTGAGTGATTGCCTGGATGGTGTAGCTATATAGAAGCGGGGCGGCCGGTAGGCCTCCCCGC‑‑GTTGTGTTTTGGGGttg                                                                                                                                >  5:236872/1‑148 (MQ=255)
                     cATGTTTGTCCTTGGTGCTACCCATGTGTGCCAAGCGGGCAAGAAGTCGTTTGGCCGGGTTACCCTCAAGGCTGAGTGATTGCCTGGATGGTGTAGCTATATAGAAGCGGGGAGGCCGGTAGGCCTCCCCGC‑‑TTTGTTTTTGGGGtt                                                                                                                                 <  1:28504/147‑1 (MQ=255)
                     cATGTTTGTCCTTGGTGCTACCCATGTGTGCCAAGCGGGCAAGAAGTCGTTTGGCCGGGTTACCCTCAAGGCTGAGTGATTGCCTGGATGGTGTAGCTATATAGAAGCGGGGAGGCCGGTAAGCCAAACCCC‑‑ATTGTTTTTTGTGtttt                                                                                                                              >  6:103914/1‑148 (MQ=255)
                               cTTGGTGCTACCCATGTGTGCCAAGCGGGCAAGAAGTCGTTTGGCCGGGTTACCCTCAAGGCTGAGTGATTGCCTGGATGGTGTAGCTATATAGAAGCGGGGGGGGCGGGTGGGCGGCCCGC‑‑GGTGTTTttgggggtttgggggtggtg                                                                                                                     >  2:79300/1‑145 (MQ=255)
                               cTTGGTGCTACCCATGTGTGCCAAGCGGGCAAGAAGTCGTTTGGCCGGGTTACCCTCAAGGCTGAGTGATTGCCTGGATGGTGTAGCTATATAGAAGCGGGGGGGGCGGGGTGGCGGCCCGC‑‑GTTTTTTTTTTGGTTttgtgggttgtg                                                                                                                    >  1:291965/1‑143 (MQ=255)
                               cTTGGTGCTACCCATGTGTGCCAAGCGGGCAAGAAGTCGTTTGGCCGGGTTACCCTCAAGGCTGAGTGATTGCCTGGATGGTGTAGCTATATAGAAGCGGGGAGGCCGGTAGGCCTACCCGC‑‑GTTGTTTTTTGGGTTTGGTGGTTGGtg                                                                                                                    >  2:27155/1‑149 (MQ=255)
                                      cTACCCTTGTGTGCCAAGCGGGCAAGAAGTCGTTTGGCCGGGTTTCCCTCAAGGCTGAGTGATTGCCTGGATGGTGTAGCTATTTAGAAGCGGGGAGGCCGGTAGGCCTCCCCGC‑‑TTTGTTTTTGGGGTTTGGTGGTTGCTGTTGTagg                                                                                                             <  1:107129/149‑1 (MQ=255)
                                      cTACCCATGTGTGCCAAGCGGGCAAGAAGTCGTTTGGCCGGGTTACCCTCAAGGCTGAGTGATTGCCTGGATGGTGTAGCTATATAGAAGCGGGGAGGCCGGTTGGCCACCCCGC‑‑GTTGTTTTTTGGGTTTTGTGGTTTCagttgtttg                                                                                                             >  2:67914/1‑146 (MQ=255)
                                      cTACCCATGTGTGCCAAGCGGGCAAGAAGTCGTTTGGCCGGGTTACCCTCAAGGCTGAGTGATTGCCTGGATGGTGTAGCTATATAGAAGCGGGGAGGCCGGTAGGCCTCCCCGC‑‑TTTGTTTTTGGGGTTTGGTGGTTGCTGTTGTagg                                                                                                             <  3:91446/149‑1 (MQ=255)
                                         cccATGTGTGCCAAGCGGGCAAGAAGTCGTTTGGCCGGGTTACCCTCAAGGCTGAGTGATTGCCTGGATGGTGTAGCTATATAGAAGCGGGGAGGCCGGTAGGCCTCCCCGC‑‑TTTGTTTTTGGGGTTTGGTGGTTGCTGTTGTaggg                                                                                                            <  3:88891/147‑1 (MQ=255)
                                              gtgtgCCAAGCGGGCAAGAAGTCGTTTGGCCGGGTTACCCTCAAGGCTGAGTGATTGCCTGGATGGTGTAGCTATATAGAAGCGGGGGGGGCGGGTGGCCGCCCCGC‑‑GTTGTTTTTTGGGTTTGGTGGTTGGGGTTGTtgggggggggg                                                                                                     >  8:184796/1‑149 (MQ=255)
                                               tgtgCCAAGCGGGCAAGAAGTCGTTTGGCCGGGTTACCCTCAAGGCTGAGTGATTGCCTGGATGGTGTAGCTATATAGAAGCGGGGAGGCCGGTAGGCCTCCCCGC‑‑TTTGTTTTTGGGGTTTGGTGGTTGCTGTTGTagggggaggggg                                                                                                    <  1:5496/149‑1 (MQ=255)
                                                         gggCATGAAGTGGTTTGGCCGGGTTACCCTCAAGGCTGAGTGATTGCCTGGATTGTGTAGCTTTATAGAAGCGGGGAGGCCGGTAGGCCTCCCCGC‑‑TTTGTTTTTg                                                                                                                                       <  7:11892/106‑1 (MQ=38)
                                                         gggCAAGAAGTCGTTTGGCCGGGTTACCCTCAAGGCTGAGTGATTGCCTGGATGGTGTAGCTATATAGAAGCGGGGAGGCCGGTAGGCCTCCCCGC‑‑TTTGTTTTTgg                                                                                                                                      <  1:17439/107‑1 (MQ=255)
                                                          ggCAAGAAGTCGTTTGGCCGGGTTACCCTCAAGGCTGAGTGATTGCCTGGATGGTGTAGCTATATAGAAGCGGGGCGGCCGGTTGGCCGGCCCGG‑‑GTTTTTTTTGGGGTTTTGGGTGTTGTGTTGTTGGGGGAGGGGGCTGtttttt                                                                                           >  5:150455/1‑145 (MQ=255)
                                                          ggCAAGAAGTCGTTTGGCCGGGTTACCCTCAAGGCTGAGTGATTGCCTGGATGGTGTAGCTATATAGAAGCGGGGAGGCCGGTAGGCCTCCCCGC‑‑TTTGTTTTTGGGGTTTGGTGGTTGCTGTTGTAGGGGGAGGGGGCGGTTTTgaga                                                                                         <  2:69671/149‑1 (MQ=255)
                                                          ggCAAGAAGTCGTTTGGCCGGGCTACCCTCAAGGCTGAGGGATTGCCTGGATGGGGTAGCTATATAGAAGCGGGGAGGCCGGTAGGCCTCCCCGC‑‑TTTGTTTTTGGGGTTTGGTGGTTGCTGTTGTAGGGGGAGGGGGCTGTTTTgaga                                                                                         <  8:306025/149‑1 (MQ=255)
                                                           gCAAGAAGTCGTTTGGCCGGGTTACCCTCAAGGCTGAGTGATTGCCTGGATGGTGTAGCTATATAGAAGCGGGGGGGCGGGGTGTCCCCCCCCCCATTTTTTTTTGGGTTTGTGGTGGTGGTTTGTTTGGGGGGGGGGGGGTttttatg                                                                                           >  6:25859/1‑146 (MQ=255)
                                                              tgtagTCGTTTGGCGGGGTTACCGTCACGTCTTTGTGATTGCCTTGTTGGTGTAGCTATATAGAAGCGGGGAGGCCGGTAGGCCTCCCCGC‑‑TTTGTTTTTGGGGTTTGGTGGTTGCTGTTGTAGGGGGAGGGGGCTGTTTTGAGAGGgt                                                                                     <  5:18006/146‑1 (MQ=37)
                                                                   tctggtGGCCGGGTTAGCCTCACGGCTGTCTGATTGCCTGGATTGTGTCGCTATATAGAAGCGGGGAGGCCGCTTGGCCTCCCCGC‑‑TTTGTTTTTGGGGTTTGGTGGTTGCTGTTGTAGGGGGAGGGGGCTGTTTTGAGAGGGTGTTgg                                                                                <  8:235086/144‑1 (MQ=37)
                                                                        tGGCCGGGTTACCCTCAAGGCTGTGTGATTGCCTGGATGGTGTTGCTATTTAGTAGCGGGGCGGCCGGTAGGCCGCCCCGC‑‑TTTGTTTTTGGGGTTTGGTGGTTGCTGTTGTAGGGGGAGGGGGCTGTTTTGAGAGGGTGTTGGCGGt                                                                            <  5:25859/148‑1 (MQ=38)
                                                                         ggCCGGGTTACCCTCAAGGCTGAGTGATTGCCTGGATGGTGTAGCTATATAGAAGCGGGGAGGCCGGTAGGCCTCCCCGC‑‑TTTGTTTTTGGGGTTTGGTGGTTGCTGTTGTAGGGGGTGGGGGCTGTTTTGAGAGGGTGTTGGCGGTg                                                                           <  1:27155/148‑1 (MQ=255)
                                                                          gcggggTTACCCTCAAGGCGGTGTTTTTGGCTGGATGGTTTTGCTATTTTGATGCGGGGCGGCGGGTAGGCCTCCCCGC‑‑TTTGTTTTTGGGGTTTGGGGTTTGGAGTTGTAGGGGGAGGGGGCTGTTTTGAGAGGGTGTTGGCGGTgag                                                                         <  1:67914/146‑1 (MQ=25)
                                                                          gCCGGGTTACCCTCAAGGCTGAGTGATTGCCTGGATGGTGTAGCTATATAGAAGCGGGGAGGCCGGTAGGCCTCCCCGC‑‑TTTGTTTTTGGGGTTTGGTGGTTGCTGTTGTAGGGGGAGGGGGCTGTTTTGAGAGGGTGTTGGCGGTga                                                                          <  5:103914/148‑1 (MQ=255)
                                                                            cGGGTTACCCTCAAGGCTGAGTGATTGCCTGGATGGTGTAGCTATATAGAAGCGGGGAGGCCGGTAGGCCTCCCCGC‑‑TTTGTTTTTGGGGTTTGGTGGTTGCTGTTGTAGGGGGAGGGGGCTGTTTTGAGAGGGTGTTGGCGGTGAGAt                                                                       <  1:197052/149‑1 (MQ=255)
                                                                                        aaGGCTGAGTGATTGCCTGGATGGTGTAGCTATATAGAAGCGGGGAGGCCGGTTGGCCTCCCCGC‑‑TTTGTTTTTTGGGTTTTGTGGTTGCTGTTGTATGGGGAGGGGGCTGTTTTTAGAGGGTGTTTGGGGTGAGATGTGTACGCcac                                                            >  5:328381/1‑148 (MQ=255)
                                                                                         tggCTGAGTGATTGCCTGGAGGTGGTCGCTATATAGAAGCGGGGAGGCCGGTAGGCCTCCCCGC‑‑TTTGTTTTTGGGGTTTGGTGGGTGCTGTTGTAGGGGGAGGGGGCTGTTTTGAGAGGGTGTTGGCGGTGAGATGTGTCCGCCACaa                                                          <  5:64123/148‑1 (MQ=25)
                                                                                           gCTGAGTGATTGCCTGGATGGTGTAGCTATATAGAAGCGGGGGGGGCCGGAGGGCGCCCCCCG‑‑TTTTTTTTTGGGGTTTTGGTTTGGGGTTTTTGGGGGGGGGGGCGGGTTTTTGTGGGTGTTTGGGGGGGGATGTGTGCGGCACaaa                                                         >  6:64123/1‑147 (MQ=255)
                                                                                              gAGTGATTGCCTGGATGGTGTAGCTATATAGAAGCGGGGAGGCCGGTAGGCCTCCCCGC‑‑TTTGTTTTTGGGGTTTGGTGGTTGCt                                                                                                                     <  5:64759/85‑1 (MQ=255)
                                                                                                     tGCCTGGATGGGGTAGCTATATAGACGCGGGGCGGCCGGTAGGCCTCCCCGC‑‑TTTGTTTTTGGGGTTTGGGGGGTGCTGTTGTAGGGGGAGGGGGCTGTTTTGAGAGGGTGTTGGCGGTGAGATGTGTCCGCCACAACCCTTGCAgcgc                                              <  7:77396/149‑1 (MQ=25)
                                                                                                        cTGGTTGGTGTCGCTATATAGAAGCGGGGAGGCGGGTAGGCCTCCCCGC‑‑TTTGTTTTTGGGGTTTGGTGGTTGCTGTTGTAGGGGGAGGGGGCTGTTTTGAGAGGGTGTTGGCGGTGAGAt                                                                       <  1:131000/121‑1 (MQ=255)
                                                                                                           gATGGTGTAGCTATATAGAAGCGGGGAGGCCGGTAGGCCTCCCCGC‑‑TTTGTTTTTGGGGTTTGGTGGTTGCTGTTGTAGGGGGAGGGGGCTGTTTTGAGAGCGTGTTGGCGGTGAGATGTGTCCGCCACAACCCTTGCAGCGCCtgt                                          <  1:229552/147‑1 (MQ=255)
                                                                                                                             aaGCGGGGAGGCCGGTAGGCCGCCCCCC‑‑GTTGTTTTTTGGGTTTGGTGGTTGGTGTTGTAGGGGGAGGGGGCTGTTTTTAGAGGGTGTTGGCGGTGAGATGTTTCCGCCACAACCCATGCAGCGCCAGTGAGATAGAGCCCCGCCcccg                      >  2:167748/1‑149 (MQ=255)
                                                                                                                                cGGGGAGGCCGGTAGGCCTCCCCGC‑‑TTTGTTTTTGGGGTTTGGTGGTTGCTGTTGTAGGGGGAGGGGGCTGTTTTGAGAGGGTGTTGGCGGTGAGATGTGTCCGCCACAACCCTTGCAGCGCCTGTGAGATCGAGCGCCGCCCGCGTcc                   >  5:126637/1‑149 (MQ=255)
                                                                                                                                        gCCGGTAGGCCTCCCCGC‑‑TTTGTTTTTGGGGTTTGGTGGTTGCTGTTGTAGGGGGAGGGGGCTGTTTTGAGAGGGTGTTGGCGGTGAGATGTGTCCGCCACAACCCTTGCAGCGCCTGTGAGATCGAGCGCCGCCCGCGTCCGGCtt              >  8:25339/1‑147 (MQ=255)
                                                                                                                                          cGGTAGGCCTCCCCGC‑‑TTTGTTTTTGGGGTTTGGTGGTTGCTGTTGTAGGGGGAGGGGGCTGTTTTGAGAGGGTGTTGGCGGTGAGATGTGTCCGCCACAACCCTTGCAGCGCCTGTGAGATCGAGCGCCGCCCGCGTCCGGCTTgaga          <  7:184796/149‑1 (MQ=255)
                                                                                                                                          cGGTAGGCCTCCCCGC‑‑TTTGTTTTTGGGGTTTGGTGGTTGCTGTTGTAGGGGGAGGGGGCTGTTTTGAGAGGGTGTTGGCGGTGAGATGTGTCCGCCACAACCCTTGCAGCGCCTGTGAGATCGAGCGCCGCCCGCGTCCGGCTTgaga          <  6:328381/149‑1 (MQ=255)
                                                                                                                                          cGGTAGGCCTCCCCGC‑‑TTTGTTTTTGGGGTTTGGTGGGTGCTGGTGTAGGGGGAGGGGGCTGTTTTGAGAGGGTGTTGGCGGTGAGATGTGTCCGCCACAACCCTTGCAGCGCCTGTGAGATCGAGCGCCGCCCGCGTCCGGCTTgaga          <  1:79300/149‑1 (MQ=255)
                                                                                                                                           ggTAGGCCTCCCCGC‑‑TTTGTTTTTGGGGTTTGGTGGTTGCTGTTGTAGGGGGAGGGGGCTGTTTTGAGAGGGTGTTGGCGGTGAGATGTGTCCGCCACAACCCTTGCAGCGCCTGTGAGATCGAGCGCCGCCCGCGTCCGGCTTgaga          >  2:217272/1‑148 (MQ=255)
                                                                                                                                                gCCTCCCCGC‑‑TTTGTTTTTGGGGTTTGGTGGTTGCTGTTGTAGGGGGAGGGGGCTGTTTTGAGAGGGTGTTGGCGGTGAGATGTGTCCGCCACAACCCTTGCAGCGCCTGTGAGATCGAGCGCCGCCCGCGTCCGGCTTGAGACATTCa    >  6:74515/1‑149 (MQ=255)
                                                                                                                                                 ccTCCCCGC‑‑TTTGTTTTTGGGGTGTGGTGGTTGCTGTTGTAGGGGGTGGGGGCTGTTTTGTGAGGGTGTTGGCGGTTAGATGTGTCCGGCAGAACCCTTGCAGCGCCTGTGAGATCGAGCGCCGCCCGCGTCCGGCTTGAGACTTTCa  <  8:79262/149‑2 (MQ=25)
                                                                                                                                                  cTCCCCGC‑‑TTTGTTTTTGGGGTTTGGTGGTTGCTGTTGTAGGGGGAGGGGGCTGTTTTGAGAGGGTGTTGGCGGTGAGATGTGTCCGCCACAACCCTTGCAGCGCCTGTGAGATCGAGCGCCGCCCGCGTCCGGCTTGAGACATTCAcc  >  6:632/1‑149 (MQ=255)
                                                                                                                                                    |                                                                                                                                                         
TGCAGTGGTCGACAAGGACTTCATGTTTGTCCTTGGTGCTACCCATGTGTGCCAAGCGGGCAAGAAGTCGTTTGGCCGGGTTACCCTCAAGGCTGAGTGATTGCCTGGATGGTGTAGCTATATAGAAGCGGGGAGGCCGGTAGGCCTCCCCGC‑‑TTTGTTTTTGGGGTTTGGTGGTTGCTGTTGTAGGGGGAGGGGGCTGTTTTGAGAGGGTGTTGGCGGTGAGATGTGTCCGCCACAACCCTTGCAGCGCCTGTGAGATCGAGCGCCGCCCGCGTCCGGCTTGAGACATTCACC  >  NC_002947/524034‑524327

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: