Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 1,236,976 C→T 9.6% intergenic (+15/+23) ybaK → / ← PP_1078 Cys‑tRNAPro deacylase/ABC transporter ATP‑binding protein

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0029471,236,9760CT 9.6% 144.1 / 4.8 52intergenic (+15/+23)ybaK/PP_1078Cys‑tRNAPro deacylase/ABC transporter ATP‑binding protein
Reads supporting (aligned to +/- strand):  ref base C (24/23);  new base T (2/3);  total (26/26)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 6.10e-01

TGACTCGGCGCAGAATTTTGCAACTATCCATGTCAGCGCCGGCAGGCGTGGCCTGGAAGTGGAGTTGGCCGCAGCGGTGCTGGCGGAGTATACGCAGGGCACGTTTGCCGGGATTGGTCGAGGCTGAGAACAGGGGGCCGCCAAGCGGCCCTCTGGTTCTGAAGTCAAGTCGCCGAACTGTATCGACGCACGCCATTCTCCTGGCGTGGCAACTGTGCCGCAACGCTACCCGGTACCGGGAATAGCACGAGGTGATCGGCCTTCACTGCAATCTCCACATTC  >  NC_002947/1236835‑1237116
                                                                                                                                             |                                                                                                                                            
tgaCTCGGCGCAGAATTTTGCAACTATCCATGTCAGCGCCGGCAGGCGTGGCCTGGAAGTGGAGTTGGCCGCAGCGGTGCTGGCAGAGTACACGCAGGGCACGTTTGCCGGGATTGGTCGAGGCTGAGAACAGGGGGCCGCCAAGCgg                                                                                                                                        <  1:99531/148‑1 (MQ=255)
      ggCGCAGAATTTTGCAACTATCCATGTCAGCGCCGGCAGGCGTGGCCTGGAAGTGGAGTTGGCCGCAGCGGTGCTGGCGGAGTATACGCAGGGCACGTTTGCCGGGATTGGTCGAGGCTGAGAACAGGGGGCCGCCAAGCGGCCCTCTg                                                                                                                                 <  5:316826/149‑1 (MQ=255)
      ggCGCAGAATTTTGCAACTATCCATGTCAGCGCCGGCAGGCGTGGCCTGGAAGTGGAGTTGGCCGCAGCGGTGCTGGCGGAGTATACGCAGGGCACGTTTGCCGGGATTGGTCGAGGCTGAGAACAGAGGGCCGCTTGGCGGCCccct                                                                                                                                  <  6:113534/148‑1 (MQ=255)
          cAGAATTTTGCAACTATCCATGTCAGCGCCGGCAGGCGTGGCCTGGAAGTGGAGTTGGCCGCAGCGGTGCTGGCGGAGTATACGCAGGGCACGTTTGCCGGGATTGGTCGAGGCTGAGAACAGGGGGCCGCCAAGCGGCCCTCTGGTTc                                                                                                                             >  7:120011/1‑149 (MQ=255)
          cAGAATTTTGCAACTATCCATGTCAGCGCCGGCAGGCGTGGCCTGGAAGTGGAGTTGGCCGCAGCGGTGCTGGCGGAGTATACGCAGGGCACGTTTGCCGGGATTGGTCGAGGCTGAGAACAGAGGGCCGCTTGGCGGCCCCCTgttct                                                                                                                             <  5:159598/149‑5 (MQ=255)
           aGAATTTTGCAAATATCCATGTCAGCGCCGGCAGGCGTGGCCTGGAAGTGGAGGTGGCCGCAGCGGTGCTGGCGGAGTATACGCAGGGCACGTTTGCCGGGATTGGTCGAGGCTGAGAACAGGGGGCCGCCAAGCGGCCCTCTGGTTCt                                                                                                                            <  1:36234/149‑1 (MQ=255)
              aTTTTGCAACTATCCATGTCAGCGCCGGCAGGCGTGGCCTGGAAGTGGAGTTGGCCGCAGCGGTGCTGGCGGAGTATACGCAGGGCACGTTTGCCGGGATTGGTCGAGGCTGAGAACAGGGGGCCGCCAAGCGGCCCTCTGGTTCTGaa                                                                                                                         <  5:80559/149‑1 (MQ=255)
              aTTTTGCAACTATCCATGTCAGCGCCGGCAGGCGTGGCCTGGAAGTGGAGTTGGCCGCAGCGGTGCTGGCGGAGTATACGCAGGGCACGTTTGCCGGGATTGGTCGAGGCTGAGAACAGGGGGCCGCCAAGCGGCCCTCTGGTTCTGaa                                                                                                                         <  5:85316/149‑1 (MQ=255)
                    cAACTATCCATGTCAGCGCCGGCAGGCGTGGCCTGGAAGTGGAGTTGGCCGCAGCGGTGCTGGCGGAGTATACGCAGGGCACGTTTGCCGGGATTGGTCGAGGCTGAGAACAGGGGGCCGCCAAGCGGCCCTCTGGTTCTGaagtcaa                                                                                                                    <  7:15638/148‑1 (MQ=255)
                     aaCTATCCATGTCAGCGCCGGCAGGCGTGGCCTGGAAGTGGAGTTGGCCGCAGCGGTGCTGGCGGAGTATACGCAGGGCACGTTTGCCGGGATTGGTCGAGGCTGAGAACAGGGGGCCGCCAAGCGGCCCTCTGGTTCTGaagtcaag                                                                                                                   >  4:167435/1‑148 (MQ=255)
                     aaCTATCCATGTCAGCGCCGGCAGGCGTGGCCTGGAAGTGGAGTTGGCCGCAGCGGTGCTGGCGGAGTATACGCAGGGCACGTTTGCCGGGATTGGTCGAGGCTGAGAACAGAGGGCCGCTTGGCGGCCCCCTgttc                                                                                                                              <  5:264056/137‑4 (MQ=255)
                     aaCTATCCATGTCAGCGCCGGCAGGCGTGGCCTGGAAGTGGAGTTGGCCGCAGCGGTGCTGGCGGAGTATACGCAGGGCACGTTTGCCGGGATTGGTCGAGGCTGAGAACAGAGGGCCGCTTGGCGGCCCCCTgttc                                                                                                                              >  6:264056/1‑134 (MQ=255)
                                    cgcCGGCAGGCGTGGCCTGGAAGTGGAGTTGGCCGCAGCGGTGCTGGCGGAGTATACGCAGGGCACGTTTGCCGGGATTGGTCGAGGCTGAGAACAGGGGGCCGCCAAGc                                                                                                                                          <  6:37108/110‑1 (MQ=255)
                                    cgcCGGCAGGCGTGGCCTGGAAGTGGAGTTGGCCGCAGCGGTGCTGGCGGAGTATACGCAGGGCACGTTTGCCGGGATTGGTCGAGGCTGAGAACAGGGGGCCGCCAAGc                                                                                                                                          >  5:37108/1‑110 (MQ=255)
                                                tGGCCTGGAAGTGGAGTTGGCCGCAGCGGTGCTGGCGGAGTATACGCAGGGCACGTTTGCCGGGATTGGTCGAGGCTGAGAACAGGGGGCCGCCAAGCGGCCCTCTGGTTCTGAAGTCAAGTCGCCGAACTGTATCGACGCACGCCAtt                                                                                       >  4:27204/1‑149 (MQ=255)
                                                  gCCTGGAAGTGGAGTTGGCCGCAGCGGTGCTGGCGGAGTATACGCAGGGCACGTTTGCCGGGATTGGTCGAGGCTGAGAACAGGGGGCCGCCAAGCGGCCCTCTGGTTCTGAAGTCAAGTCGCCGAACTGTATCGACGCACGCCATtct                                                                                     >  3:63320/1‑149 (MQ=255)
                                                    cTGGAACTGGAGTTGGCCGCAGCGGTGCTGGAGGAGTATACGCACGGCACGTTTGCCGGGATTGGTCGAGGCTGAGAACAAGGAGCCGCCAAGCGGCCCTCTGGTTCTGAAGTCAAGTCGCCGAACTGTATCGACGCACGCCATTCTcc                                                                                   >  7:206505/1‑149 (MQ=255)
                                                       gAAGTGGAGTTGGCCGCAGCGGTGCTGGCGGAGTATACGCAGGGCACGTTTGCCGGGATTGGTCGAGGCTGAGAACAGGGGGCCGCCAAGCGGCCCTCTGGTTCTGAAGTCAAGTCGCCGAACTGTATCGACGCACGCCATTCTCCTg                                                                                 <  8:120011/148‑1 (MQ=255)
                                                                  ggCCGCAGCGGTGCTGGCGGAGTATACGCAGGGCACGTTTGCCGGGATTGGTCGAGGCTGAGAACAGGGGGCCGCCAAGCGGCCCTCTGGTTCTGAAGTCAAGTCGCCGAACTGTATCGACGCACGCCATTCTCCTGGCGTGGCAACtg                                                                     >  1:322869/1‑149 (MQ=255)
                                                                  ggCCGCAGCGGTGCTGGCGGAGTATACGCAGGGCACGTTTGCCGGGATTGGTCGAGGCTGAGAACAGGGGGCCGCCAAGCGGCCCTCTGGTTCTGAAGTCAAGTCGCCGAACTGTATCGACGCACGCCATTCTCCTGGCGTGGCAACtg                                                                     >  3:150378/1‑149 (MQ=255)
                                                                            gTGCTGGCGGAGTATACGCAGGGCACGTTTGCCGGGATTGGTCGAGGCTGAGAACAGGGGGCCGCCAAGCGGCCCTCTGGTTCTGAAGTCAAGTCGCCGAACTGTATCGACGCACGCCATTCTCCTGGCGTGGCAACTGTGCCGCAACg                                                           >  4:249662/1‑149 (MQ=255)
                                                                                tGGCGGAGTATACGCAGGGCACGTTTGCCGGGATTGGTCGAGGCTGAGAACAGGGGGCCGCCAAGCGGCCCTCTGGTTCTGAAGTCAAGTCGCCGAACTGTATCGACGCACGCCATTCTCCTGGCGTGGCAACTGTGCCGCAACGCTAc                                                       >  4:132894/1‑149 (MQ=255)
                                                                                 ggCGGAGTATACGCAGGGCACGTTTGCCGGGATTGGTCGAGGCTGAGAACAGGGGGCCGCCAAGCGGCCCTCTGGTTCTGAAGTCAAGTCGCCGAACTGTATCGACGCACGCCATTCTCCTGGCGTGGCAACTGTGCCGCAACGCTAcc                                                      <  8:130622/149‑1 (MQ=255)
                                                                                  gCGGAGTATACGCAGGGCACGTTTGCCGGGATTGGTCGAGGCTGAGAACAGGGGGCCGCCAAGCGGCCCTCTGGTTCTGAAGTCAAGTCGCCGAACTGTATCGACGCACGCCATTCTCCTGGCGTGGCAACTGTGCCGCAACGCTAccc                                                     >  3:61140/1‑149 (MQ=255)
                                                                                  gCGGAGTATACGCAGGGCACGTTTGCCGGGATTGGTCGAGGCTGAGAACAGGGGGCCGCCAAGCGGCCCTCTGGTTCTGAAGTCAAGTCGCCGAACTGTATCGACGCACGCCATTCTCCTGGCGTGGCAACTGTGCCGCAACGCTAccc                                                     >  3:31722/1‑149 (MQ=255)
                                                                                  gCGGAGTATACGCAGGGCACGTTTGCCGGGATTGGTCGAGGCTGAGAACAGGGGGCCGCCAAGCGGCCCTCTGGTTCTGAAGTCAAGTCGCCGAACTGTATCGACGCACGCCATTCTCCTGGCGTGGCAACTGTGCCGCAACGCTAcc                                                      <  3:167435/148‑1 (MQ=255)
                                                                                        tataCGCAGGGCACGTTTGCCGGGATTGGTCGAGGCTGAGAACAGGGGGCCGCCAAGCGGCCCTCTGGTTCTGAAGTCAAGTCGCCGAACTGTATCGACGCACGCCATTCTCCTGGCGTGGCAACTGTGCCGCAACGCTACCCGGTAcc                                               >  1:299450/1‑149 (MQ=255)
                                                                                          taCGCAGGGCACGTTTGCCGGGATTGGTCGAGGCTGAGAACAGGGGGCCGCCAAGCGGCCCTCTGGTTCTGAAGTCAAGTCGCCGAACTGTATCGACGCACGCCATTCTCCTGGCGTGGCAACTGTGCCGCAACGCTACCCGGTACCgg                                             <  2:299450/149‑1 (MQ=255)
                                                                                              cAGGGCACGTTTGCCGGGATTGGTCGAGGCTGAGAACAGGGGGCCGCCAAGCGGCCCTCTGGTTCTGAAGTCAAGTCGCCGCACTGTATCGACGCACGCCATTCTCCTGGCGTGGCAACTGTGCCGCAACGCTACCCGGTACCGGGaa                                          >  3:99889/1‑148 (MQ=255)
                                                                                                  gCACGTTTGCCGGGTTTGGTCGAGGCTGAGAACAGGGGGCCGCCAAGCGGCCCTCTGGTTCTGAAGTCAAGTCGCCGAACTGTATCGACGCACGCCATTCTCCTGGCGTGGCAACTGTGCCGCAACGCTACCCGGTACCGGGAATAGCa                                     <  8:206505/149‑1 (MQ=255)
                                                                                                  gCACGTTTGCCGGGATTGGTCGAGGCTGAGAACAGGGGGCCGCCAAGCGGCCCTCTGGTTCTGAAGTCAAGTCGCCGAACTGTATCGACGCACGCCATTCTCCTGGCGTGGCAACTGTGCCGCAACGCTACCCGGTACCGGGAATAGCa                                     >  3:94800/1‑149 (MQ=255)
                                                                                                  gCACGTTTGCCGGGATTGGTCGAGGCTGAGAACAGGGGGCCGCCAAGCGGCCCTCTGGTTCTGAAGTCAAGTCGCCGAACTGTATCGACGCACGCCATTCTCCAGGCGTGGCAACTGTGCCGCAACGCTACCCGGTACCGGGAATAGCa                                     <  6:319582/149‑1 (MQ=255)
                                                                                                    aCGTTTGCCGGGATTGGTCGAGGCTGAGAACAGGGGGCCGCCAAGCGGCCCTCTGGTTCTGAAGTCAAGTCGCCGAACTGTATCGACGCACGCCATTCTCCTGGCGTGGCAACTGTGCCGCAACGCTACCCGGTACCggg                                            >  7:176189/1‑140 (MQ=255)
                                                                                                    aCGTTTGCCGGGATTGGTCGAGGCTGAGAACAGGGGGCCGCCAAGCGGCCCTCTGGTTCTGAAGTCAAGTCGCCGAACTGTATCGACGCACGCCATTCTCCTGGCGTGGCAACTGTGCCGCAACGCTACCCGGTACCggg                                            <  8:176189/140‑1 (MQ=255)
                                                                                                     cGTTTGCCGGGATTGGTCGAGGCTGAGAACAGGGGGCCGCCAAGCGGCCCTCTGGTTCTGAAGTCAAGTCGCCGAACTGTATCGACGCACGCCATTCTCCTGGCGTGGCAACTGTGCCGCAACGCTACCCGGTACCGGGAATAGCACGa                                  <  4:63320/149‑1 (MQ=255)
                                                                                                      gTTTGCCGGGATTGGTCGAGGCTGAGAACAGGGGGCCGCCAAGCGGCCCTCTGGTTCTGAAGTCAAGTCGCCGAACTGTATCGACGCACGCCATTCTCCTGGCGTGGCAACTGTGCCGCAACGCTACCCGGTACCGGGAATAGCACGAg                                 <  3:27204/149‑1 (MQ=255)
                                                                                                       tttGCCGGGATTGGTCGAGGCTGAGAACAGGGGGCCGCCAAGCGGCCCTCTGGTTCTGAAGTCAAGTCGCCGAACTGTATCGACGCACGCCATTCTCCTGGCGTGGCAACTGTGCCGCAACGCTACCCGGTACCGGGAATAGCACGAgg                                <  4:61140/149‑1 (MQ=255)
                                                                                                          gCCGGGATTGGTCGAGGCTGAGAACAGGGGGCCGCCAAGCGGCCCTCTGGTTCTGAAGTCAAGTCGCCGAACTGTATCGACGCACGCCATTCTCCTGGCGTGGCAACTGTGCCGCAACGCTACCCGGTACCgg                                             <  2:232509/133‑1 (MQ=255)
                                                                                                          gCCGGGATTGGTCGAGGCTGAGAACAGGGGGCCGCCAAGCGGCCCTCTGGTTCTGAAGTCAAGTCGCCGAACTGTATCGACGCACGCCATTCTCCTGGCGTGGCAACTGTGCCGCAACGCTACCCGGTACCgg                                             >  1:232509/1‑133 (MQ=255)
                                                                                                            cGGGATTGGTCGAGGCTGAGAACAGGGGGCCGCCAAGCGGCCCTCTGGTTCTGAAGTCAAGTCGCCGAACTGTATCGACGCACGCCATTCTCCTGGCGTGGCAACTGTGCCGCAACGCTACCCGGTACCAGGAATAGCAAGAGGTGATc                           >  8:248732/1‑149 (MQ=255)
                                                                                                                   ggTCGAGGCTGAGAACAGGGGGCCGCCAAGCGGCCCTCTGGTTCTGAAGTCAAGTCGCCGAACTGTATCGACGCACGCCATTCTCCTGGCGTGGCAACTGTGCCGCAACGCTACCCGGTACCGGGAATAGCACGAGGTGATCGGCCtt                     <  4:150378/148‑1 (MQ=255)
                                                                                                                   ggTCGAGGCTGAGAACAGGGGGCCGCCAAGCGGCCCTCTGGTTCTGAAGTCAAGTCGCCGAACTGTATCGACGCACGCCATTCTCCTGGCGTGGCAACTGTGCCGCAACGCTACCCGGTACCGGGAATAGCACGAGGTGATCGGCCtt                     <  2:322869/148‑1 (MQ=255)
                                                                                                                     tCGAGGCTGAGAACAGGGGGCCGCCAAGCGGCCCTCTGGTTCTGAAGTCAAGTCGCCGAACTGTATCGACGCACGCCATTCTCCTGGCGTGGCAACTGTGCCGCAACGCTACCCGGTACCg                                              <  4:79267/121‑1 (MQ=255)
                                                                                                                     tCGAGGCTGAGAACAGGGGGCCGCCAAGCGGCCCTCTGGTTCTGAAGTCAAGTCGCCGAACTGTATCGACGCACGCCATTCTCCTGGCGTGGCAACTGTGCCGCAACGCTACCCGGTACCg                                              >  3:79267/1‑121 (MQ=255)
                                                                                                                         ggCTGAGAACAGGGGGCCGCCAAGCGGCCCTCTGGTTCTGAAGTCAAGTCGCCGAACTGTATCGACGCACGCCATTCTCCTGGCGTGGCAACTGTGCCGCaa                                                             >  4:118980/1‑102 (MQ=255)
                                                                                                                         ggCTGAGAACAGGGGGCCGCCAAGCGGCCCTCTGGTTCTGAAGTCAAGTCGCCGAACTGTATCGACGCACGCCATTCTCCTGGCGTGGCAACTGTGCCGCaa                                                             <  3:118980/102‑1 (MQ=255)
                                                                                                                            tGAGAACAGGGGGCCGCCAAGCGGCCCTCTGGTTCTGAAGTCAAGTCGCCGAACTGTATCGACGCACGCCATTCTCCTGGCGTGGCAACTGTGCCGCAACGCTAccc                                                     <  4:64984/107‑1 (MQ=255)
                                                                                                                            tGAGAACAGGGGGCCGCCAAGCGGCCCTCTGGTTCTGAAGTCAAGTCGCCGAACTGTATCGACGCACGCCATTCTCCTGGCGTGGCAACTGTGCCGCAACGCTAccc                                                     >  3:64984/1‑107 (MQ=255)
                                                                                                                                aaCAGGGGGCCGCCAAGCGGCCCTCTGGTTCTGAAGTCAAGTCGCCGAACTGTATCGACGCACGCCATTCTCCTGGCGTGGCAACTGTGCCGCAACGCTACCCGGTACCGGGAATAGCACGAGGTGATCGGCCTTCACTGCAATCTCca       >  1:243643/1‑149 (MQ=255)
                                                                                                                                  cAGAGGGCCGCTTGGCGGCCCTCTGGTTCTGAAGTCAAGTCGCCGAACTGTATCGACGCACGCCATTCTCCTGGCGTGGCAACTGTGCCGCAACGCTACCCGGTACCGGGAATAGCACGAGGTGATCGGCCTTCACTGCAATCTCcaca     >  1:310955/1‑149 (MQ=255)
                                                                                                                                   aGGGGGCCGCCAAGCGGCCCTCTGGTTCTGAAGTCAAGTCGCCGAACTGTATCGACGCACGCCATTCTCCTGGCGTGGCAAATGTGCCGCAACGCTACCCGGTACCGGGAATAGCACGAGGTGATCGGCCTTCACTGCAATCTCCACAt    >  6:178573/1‑149 (MQ=255)
                                                                                                                                     ggggCCGCCAAGCGGCCCTCTGGTTCTGAAGTCAAGTCGCCGAACTGTATCGACGCACGCCATTCTCCTGGCGTGGCAACTGTGACGCAACGCTACCCGGTACCGGGAATAGCACGAGGTGATCGGCCTTCACTGCAATCTCCACATtc  >  8:291249/1‑149 (MQ=255)
                                                                                                                                             |                                                                                                                                            
TGACTCGGCGCAGAATTTTGCAACTATCCATGTCAGCGCCGGCAGGCGTGGCCTGGAAGTGGAGTTGGCCGCAGCGGTGCTGGCGGAGTATACGCAGGGCACGTTTGCCGGGATTGGTCGAGGCTGAGAACAGGGGGCCGCCAAGCGGCCCTCTGGTTCTGAAGTCAAGTCGCCGAACTGTATCGACGCACGCCATTCTCCTGGCGTGGCAACTGTGCCGCAACGCTACCCGGTACCGGGAATAGCACGAGGTGATCGGCCTTCACTGCAATCTCCACATTC  >  NC_002947/1236835‑1237116

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