Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 524,178 C→G 20.4% intergenic (+45/‑771) 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,1780CG20.4% 83.9 / 3.7 42intergenic (+45/‑771)tyrS/PP_16SCtyrosine‑‑tRNA ligase/16S ribosomal RNA
Reads supporting (aligned to +/- strand):  ref base C (4/29);  new base G (7/2);  total (11/31)
Fisher's exact test for biased strand distribution p-value = 3.56e-04
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

CAGTGGTCGACAAGGACTTCATGTTTGTCCTTGGTGCTACCCATGTGTGCCAAGCGGGCAAGAAGTCGTTTGGCCGGGTTACCCTCAAGGCTGAGTGATTGCCTGGATGGTGTAGCTATATAGAAGCGGGGAGGCCGGTAGGCCTCCCCGCTTTGTTTTTGGGGTTTGGTGGTTGCTGTTGTAGGGGGAGGGGGCTGTTTTGAGAGGGTGTTGGCGGTGAGATGTGTCCGCCACAACCCTTGCAGCGCCTGTGAGATCGAGCGCCGCCCGCGTCCGGCTTGAGACATTCAC  >  NC_002947/524036‑524326
                                                                                                                                               |                                                                                                                                                     
cAGTGGTCGACAAGGACTTCATGTTTGTCCTTGGTGCTACCCATGTGTGCCAAGCGGGCAAGAAGTCGTTTGGCCGGGTTACCCTCAAGGCTGAGTGATTGCCTGGATGGTGTAGCTATATCGAAGCGGGGAGGCCGGTAGGCCTccc                                                                                                                                                  <  1:210950/148‑1 (MQ=255)
  gTGGTCGACAAGGACTTCATGTTTGTCCTTGGTGCTACCCATGTGTGCCAAGCGGGCAAGAAGTCGTTTGGCCGGGTTACCCTCAAGGCTGAGTGATTGCCTGGATGGTGTAGCTATATAGAAGCGGGGGGGCCGCGAGGccagccccc                                                                                                                                               >  1:359902/1‑142 (MQ=255)
  gTGGTCGACAAGGACTTCATGTTTGTCCTTGGTGCTACCCATGTGTGCCAAGCGGGCAAGAAGTCGTTTGGCCGGGTTACCCTCAAGGCTGAGTGATTGCCTGGATGGTGTAGCTATATAGAAGCGGGGAGGCCGGTAGGCCTCCCCGc                                                                                                                                               <  4:251947/149‑1 (MQ=255)
            aGGACTTCATGTTTGTCCTTGGTGCTACCCATGTGTGCCAAGCGGGCAAGAAGTCGTTTGGCCGGGTTACCCTCAAGGCTGAGTGATTGCCTGGATGGTGTAGCTATATAGAAGCGGGGAGGCCGGTAGGCCTCCCCGCTTTGTTTTTg                                                                                                                                    <  1:11841/149‑1 (MQ=255)
                 ttCATGTTTGTCCTTGGTGCTACCCATGTGTGCCAAGCGGGCAAGAAGTCGTTTGGCCGGGTTACCCTCAAGGCTGAGTGATTGCCTGGATGGTGTAGCTATATAGAAGCGGGGGGGCCGGATGGGCTGCCCCCATTGGTTTTTGGGGt                                                                                                                                >  4:327230/1‑149 (MQ=255)
                     tGTTTGTCCTTGGTGCTACCCATGTGTGCCAAGCGGGCAAGAAGTCGGTTGGCCGGGTTACCCTCAAGGCTGAGTGATTGCCTGGATGGTGTAGCTATATAGAAGCGGGGAGGCCGGTAGGCCTCCCCGCTTTGTTTTTGGGGTTtggt                                                                                                                           <  4:86338/149‑1 (MQ=255)
                      gTTTGTCCTTGGTGCTACCCATGTGTGCCAAGCGGGCAAGAAGTCGTTTGGCCGGGTTACCCTCAAGGCTGAGTGATTGCCTGGATGGTGTAGCTATATAGAAGCGGGGGGGGGGGTTGGCCGGCCCGCTGTGTTTTTGGGGTTtgggg                                                                                                                          >  6:88461/1‑147 (MQ=255)
                      gTTTGTCCTTGGTGCTACCCATGTGTGCCAAGCGGGCAAGAAGTCGTTTGGCCGGGTTACCCTCAAGGCTGAGTGATTGCCTGGATGGTGTAGCTATATAGAAGCGGGGGGGGCGGGGGGGCGGCCCGCGATGGTTTTTGGGGTTtggt                                                                                                                           >  1:394612/1‑149 (MQ=255)
                          ttctttGGGGATACCCTTGTGTGCAAACCTGGAAATACGTCTTTGGGCGGGGTTCCCCTCATGGGTGAGTTATGGCCTGTTGGGTTTAGATATATACACGGGGGGCGGCCGTTAGGCCCCCCCGCTTTGTTTTTGGGGTTTGGTGGtt                                                                                                                       <  6:56258/144‑1 (MQ=255)
                          gTCCTTGGTGCTACCCATGTGTGCCAAGCGGGCAAGAAGTCGTTTGGCCGGGTTACCCTCAAGGCTGAGTGATTGCCTGGATGGTGTTGCTATATAGAAGCGGGGAGGCCGGTAGGCCTCCCCGCTTTGTTTTTGGGGTTTGGTGGtt                                                                                                                       <  3:327230/148‑1 (MQ=255)
                           tCCTTGGTGCTACCCATGTGTGCCAAGCGGGCAAGAAGTCGTTTGGCCGGGTTACCCTCAAGGCTGAGTGATTGCCTGGATGGTGTAGCTATATAGAAGCGGGGAGGCCGGTAGGCCTCCCCGCTTTGTTTTTGGGGTTTGGTGGTTg                                                                                                                      <  5:312138/148‑1 (MQ=255)
                             cTTGGTGATACCCATGTGTGCCAAGCGGGCAATAAGTCGTTTGGCGGGGTTACCCTCAAGGCTGAGTGATTGCCTGGATGGTGTAGCTATTTAGAAGGGGGGAGGCCGGTAGGCCTCCCCGCTTTGTTTTTGGGGTTTGGTGGTTGCtg                                                                                                                   <  4:333984/149‑1 (MQ=39)
                              ttGGTGCTACCCATGTGTGCCAAGCGGGCAAGAAGTCGTTTGGCCGGGTTACCCTCAAGGCTGAGTGATTGCCTGGATGGTGTAGCTATATAGAAGCGGGGGGGGGGGGAGGGCCCGCCCCCGTGTTTTTTggggttgggggttgggtg                                                                                                                  >  2:145794/1‑145 (MQ=255)
                                   gCTACCCATGTGTGCCAAGCGGGCAAGAAGTCGTTTGGCCGGGTTACCCTCAAGGCTGAGTGATTGCCTGGATGGTGTAGCTATATAGAAGCGGGGCGGGCGGGGGGGCAGCCCGCGGTGTGTTTGGGGGTTGGTGGTTGGGGTTGTag                                                                                                             >  4:288035/1‑149 (MQ=255)
                                     tACCCATGTGTGCCAAGCGGGCAAGAAGTCGTTTGGCCGGGTTACCCTCAAGGCTGAGTGATTGCCTGGATGGTGTAGCTATATAGAAGCGGGGAGGCCGGTAGGCCTCCCCGCTTTGTTTTTGGGGTTTGGTGGTTGCTGTTGTaggg                                                                                                           <  1:169076/149‑1 (MQ=255)
                                        ccATGTGTGCCAAGCGGGCAAGAAGTCGTTTGGCCGGGTTACCCTCAAGGCTGAGTGATTGCCTGGATGGTGTAGCTATATAGAAGCGGGGCAGCGGGGAGGCCTCCCCCCGTTTTTTTTGTGGGATGTGGGTGTCTGTTGTagtggga                                                                                                        >  4:243653/1‑149 (MQ=255)
                                          aTGTGTGCCAAGCGGGCAAGAAGTCGTTTGGCCGGGTTACCCACAAGGCTGAGTGATTGCCTGGATGGTGTAGCTATATAGAAGCGGGGGGGCCGCTAGGGCGGACCCCGTTTTTTTTTGGGGTTGGTGGTTGGGGTTGTtggggggcg                                                                                                      >  1:272473/1‑146 (MQ=255)
                                             tgtgCCAAGCGGGCAAGAAGTCGTTTGGCCGGGTTACCCTCAAGGCTGAGTGATTGCCTGGATGGTGTAGCTATATAGAAGCGGGGAGGCCGGTAGGCCTCCCCGCTTTGTTTTTGGGGTTTGGTGGTTGCTGTTGTagggggagggg                                                                                                    <  5:139397/148‑1 (MQ=255)
                                               tgCCAAGCGGGCAAGAAGTCGTTTGGCCGGGTTACCCTCAAGGCTGAGTGATTGCCTGGATGGTGTAGCTATATAGAAGCGGGGGGGCCGGGAGGGGGGCCCCCCGTTTTTTTTTGGTTTGGTGGGTGGTGTTGTtgtggggggggggg                                                                                                 >  5:5656/1‑147 (MQ=255)
                                                       gggCAAGAAGTCGTTTGGTCGGGTTACCCTCAAGGCTGAGTGATTGCCTGGATGGTGTAGCTATATAGAAGCGGGGGAGCCCGTAGTGCTGCCCGCGTTGTTTTTTGGGGTTTGTGGTTGGGGTTGTATTGGGGGGGTGGGGtttt                                                                                            >  1:266741/1‑146 (MQ=255)
                                                        ggCAATAAGTCGTTTGGCCGGGTTACCCTCACGGCTGAGTGATTGCCTGGATGGGGTAGCTATATAGAAGCGGGGAGGCCGGTAGGCCTCCCCGCTTTGTTTTTGGGGTTTGGTGGTTGCTGTTGTAGGGGGAGGGGGCTGTTTTgaga                                                                                        <  1:145794/149‑1 (MQ=255)
                                                        ggCAAGAATTGGTTTGGTCGGGTTACCCTCAAGGCTGAGTGATTGCGTGGATGGTGTAGCTATATAGAAGCGGGGAGGCCGGTAGGCCTCCCCGCTTTGTTTTTGGGGTTTGGTGGTTGCTGTTGTAGGGGGAGGGGGCTGTTTTgaga                                                                                        <  2:266741/149‑1 (MQ=255)
                                                          cAAGCAGTCGTTTGGCCGGTTGACCGTCAAGGAGGCGTGCTTGCCTGGATGGGGGAGCTATATAGAAGCGGGGCGGACGGTAGGCCTCCCCGCTTTGTTTTTGGGGTTTGGTGGTTGCTGTTGTAGGGGGAGGGGGCTGTTTTGAGAgg                                                                                      <  5:117096/149‑1 (MQ=37)
                                                          cAAGAAGGCGTTTGGCCGGGTTACCCTCAAGGCTGAGTGTTGGCCTGGATGGTGTAGATATATAGACGCGGGGAGGCCTGTAGGCCTCCCCGCTTTGTTTTTGGGGTTTGGTGGTTGCTGTTGTAGGGGGAGGGGGCTGTTTTGAGAgg                                                                                      <  3:30018/149‑1 (MQ=38)
                                                          aaagaatTCGTTAGCCCGGGGTACCCACCAGGATGTCTGATTGCCTGGATGGTGTAGGTATATAGAAGCGGGGAGGCCGGTAGGCCTCCCCGCTTTGTTTTTGGGGTTTGGTGGTTGCTGTTGTAGGGGGAGGGGGCTGTTTTGAGAgg                                                                                      <  3:243653/148‑1 (MQ=37)
                                                                 tCGTTTGGCCGGGTTACCCTCAAGGCTTAGTGATTGCCTGGATGGTGTAGCTATATAGAAGCGGGGAGGCCGGTAGGCCTCCCCGCTTTGTTTTTGGGGTTTGGTGGTTGCTGTTGTAGGGGGAGGGGGCTGTTTTGAGAGGGTGTTgg                                                                               <  3:64547/149‑1 (MQ=255)
                                                                          ggggTTCCCCGCAAGGCTTCGTGATTGCTGGGTGGGTGCAGCTAAATACACCCGGGGCGTCGGCTAGGCCGCCGGGCTTTGTTTTTGGGGTTTGGTGGTTGCTGTTGTAGGGGGAGGGGGCTGTTTTGAGAGTGTGTTGGCGGTGAGAt                                                                      <  8:239312/148‑1 (MQ=255)
                                                                          ctggTTAGCCTCATGGGTGAGTTTTTGCCTGTTTGGGGTAGCTATTTCGTAGCGGGGAGGCGTTTTGGGCTGCGCGCTTTGTTTTTGGGGTTTGGTGGTTGCTGTTGTAGGGGGTGGGGGCTGTTTTGAGAGGGTGTTGGCGGTGAGAt                                                                      <  3:158365/147‑1 (MQ=19)
                                                                          cgtgTTACCCTCCAGGCTGTGTGATTGCCGGGTTGGTGTAGTTATTTAGAAGCGGGTAGGCGGGTAGGCCTCCCCGCTTTGTTTTTGGGGTTTGGTGGTTGCTGTTGTAGGGGGAGGGGGCTGTTTTGAGAGGGTGTTGGCGGTGAGAt                                                                      <  2:359902/146‑1 (MQ=37)
                                                                          cGGGTTACCCTCAAGGCTGAGTGATTGCCTGGATGGTGTAGCTATATAGAAGCGGGGAGGCCGGTAGGCCTCCCCGCTTTGTTTTTGGGGTTTGGTGGTTGCTGTTGTAGGGGGAGGGGGCTGTTTTGAGAGGGTGTTGGCGGTGAGAt                                                                      <  5:88461/149‑1 (MQ=255)
                                                                            ggTTACCCTCAACGCTGAGTGATTGCCTGGCTGGTGTAGCTATATATAAGCGGGGAGGCCGGTAGGCCTCCCCGCTTTGTTTTTGGGGTTTGGTGGTTGCTGTTGTAGGGGGAGGGGGCTGTTTTGAGAGGGTGTTGGCGGTGAGAtgt                                                                    <  7:257772/149‑1 (MQ=255)
                                                                                        ggCTGAGTGATTGCCTGGATGGTGTAGCTATTTAGAAGCGGGGATGGCGGGAGGGCGACCCGCGTTGGTTTTTGGGTTTGTGGGGTGGTGTTGTAAGGGGAGCGGGGTGTTTTGACGTGCTGTTGGGTGTGCTATGTGTCGGCCCCAAc                                                        >  8:280555/1‑149 (MQ=255)
                                                                                             tgtgtTTGCCTGTAGTGTTTAGCTATTTAGAAGCTGGGAGGCCGGTAGGCCTCCCCGCTTTGTTTTTGGGGTTTGGTGGTTGCTGTTGTAGGGGGAGGGGGCTGTTTTGAGAGGGTGTTGGCGGTGAGATGTGTCCGCCACAACCCTTg                                                   <  8:320787/148‑1 (MQ=21)
                                                                                                  ttGCCTGGATGGTGTAGCTATATAGAAGCGGGGAGGCCGGGAGGCCTCCCCGGGTTGTTTTTGGGGTTTGGTGGTTGGTGGTGTAGGGGGAGGGGGCTGGTTTGAGAGGGTGTTGGCGGTGAGATTTGTGCGCCCCAACCCATGCAacg                                              >  5:378695/1‑149 (MQ=255)
                                                                                                             gggtAGCTATATAGAAGCGGGGAGGCCGGTAGGCCTCCCCGCTTTGTTTTTGGGGTTTGGTGGTTGCTGTTGTAGGGGGAGGGGGCTGTTTTGAGAGGGTGTTGGCGGt                                                                           <  6:390080/107‑1 (MQ=255)
                                                                                                              tgtAGCTATATAGAAGCGGGGTGGCCGGTAGGCCTCCCCGCTTTGTTTTTGGGGTTTGGTGGTTGCTGTTGTAGGGGGAGGGGGCTGTTTTGAGAGGGTGTTGGCGGTGAGATGTGTCCGCCACATCCCTTGCAGCGCCTGTGAGATc                                   <  2:105203/148‑1 (MQ=255)
                                                                                                              tgtAGCTATATAGAAGCGGGGAGGCCGGGAGGCCTCCCCGCTTTGTTTTTGGGGTTTGGTGGTTGCTGTTGTAGGGGGAGGGGGCTGTTTTGAGAGGGTGTTGGCGGTGAGATGTGTCCGCCACAACCCTTGCAGCGCCTGTGAGATc                                   <  7:280555/148‑1 (MQ=255)
                                                                                                                          gcaGCGGGGAGGCCGGTAGGCCTCCCCGCTTTGTTTTAGGGGTTTGGTGGTTGCTGTTTTAGGGGGAGGGGGCTGTTTTGAGAGGGTGTTGGCGGTGATATGTGTCCGCCACAACCCTTGCAGCGCCTGTGAGATCGAGCGccgcccgc                      <  2:272473/147‑1 (MQ=255)
                                                                                                                               ggggCGGCCGGTAGGCCTCCCCGCTTTGTTTTTGGGGTTTGGTGGTGGCTGTTGTAGGGGGAGGGGGCTGTTTTGAGAGGGTGTTGGCGGTGAGATGTGTCCGCCACAACCCTTGCAGCGCCTGTGAGATCGAGCGCCGCCCGCGTCCg                 <  5:210166/149‑1 (MQ=255)
                                                                                                                                   tggCCGGTAGGCCCCCCCGCTTTGTTTTTGGGGTGTGGTGGTTGCTGTTCTAGGGGGAGGGGGCTGTTTTGAGAGGGTGTTGGCGGTGAGATGTGTCCGCCACAACCCTTGCAGCGCCTGTGAGATCGAGCGCCGCCCGCGTCCGGCtt             <  7:138798/148‑1 (MQ=255)
                                                                                                                                    ggCCGGTAGGCCTCCCCGCTTTGTTTTTGGGGTTTGGTGGTTGCTGTTGTAGGGGGAGGGGGCTGTTTTGAGAGGGGGTTGGCGGTGAGATGTGTCCGCCACAACCCTTGCAGCGCCTGTGAGATCGAGCGCCGCCCGCGTCCGGCTTg            <  6:378695/149‑1 (MQ=255)
                                                                                                                                     gcgggTCGGCCGCCCCGCTTTGTTTTTGGGGGTTGGGGGGTGCTGTTGTAGGGGGAGGGGGCTGTTTTGAGAGGGTGTTGGCGGGGAGATGTGTCCGCCACAACCCTTGCAGCGCCTGTGAGATCGAGCGCCGCCCGCGTCCGGCTTga           <  5:231638/146‑1 (MQ=21)
                                                                                                                                      gcctggaGGGCTCGCCGCTTTGTTTTTGGCGTTTGGGGGTTGCTGTTGTAGGGGGAGGGGGCTGTTTTGAGAGGGTGTTGGCGGTGAGATGTGTCCTCCACAACCCTTGCAGCGCCTGTGAGATCGAGCGCCGCCCGCGTCCGGCTTga           <  5:236908/143‑1 (MQ=21)
                                                                                                                                               ccTCCCCGCTTTGTTTTTGGGGTTTGGTGGTTGCTGTTGTAGGGGGAGGGGGCTGTTTTGAGAGGGTGTTGGCGGTGAGATGTGTCCGCCACAACCCTTGCAGCGCCTGTGAGATCGAGCGCCGCCCGCGTCCGGCTTGAGACATTCAc  >  2:28419/1‑149 (MQ=255)
                                                                                                                                               |                                                                                                                                                     
CAGTGGTCGACAAGGACTTCATGTTTGTCCTTGGTGCTACCCATGTGTGCCAAGCGGGCAAGAAGTCGTTTGGCCGGGTTACCCTCAAGGCTGAGTGATTGCCTGGATGGTGTAGCTATATAGAAGCGGGGAGGCCGGTAGGCCTCCCCGCTTTGTTTTTGGGGTTTGGTGGTTGCTGTTGTAGGGGGAGGGGGCTGTTTTGAGAGGGTGTTGGCGGTGAGATGTGTCCGCCACAACCCTTGCAGCGCCTGTGAGATCGAGCGCCGCCCGCGTCCGGCTTGAGACATTCAC  >  NC_002947/524036‑524326

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