Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 5,268,952 T→A 7.8% intergenic (+20/+16) PP_4643 → / ← radA xanthine/uracil permease family protein/DNA repair enzyme

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0029475,268,9520TA 7.8% 149.9 / 3.1 51intergenic (+20/+16)PP_4643/radAxanthine/uracil permease family protein/DNA repair enzyme
Reads supporting (aligned to +/- strand):  ref base T (21/26);  new base A (2/2);  total (23/28)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 7.63e-01

CTGCGTGATGGGGACTACTTCAAGGATGGCAAGCTGTTCGATGCCGAGGGTAACGAAGTGCCGATGCTGGAACTGGACGAGCACGGCAATGAAACGGTCAGGCGCACGGCCGTTGCCGAGCACTGACGGGTTGCAGGAGCGAGCTTGCTCGCTCCTGCACGCTATTCGAACAATGCATCCAGCGCCTGTTCCAGCCGGGTAACGGCAATCACCTGCAACCCCGCAGGCGCTTCTTTCGGCGCATTGCCCTTGGGCACGATGGCACGCTTGAAGCCGTGCTTGGCCGCTT  >  NC_002947/5268807‑5269095
                                                                                                                                                 |                                                                                                                                               
cTGCGTGATGGGGACTACTTCAAGGATGGCAAGCTGTTCGATGCCGAGGGTAACGAAGTGCCGATGCTGGAACTGGACGAGCACGGCAATGAAACGGTCAGGCGCACGGCCGTTGCCGAGCACTGACGGGTTGCAGGAGCGAGCTTgct                                                                                                                                              <  7:148845/149‑1 (MQ=255)
     tGATGGGGACTACTTCAAGGATGGCAAGCTGTTCGATGCCGAGGGTAACGAAGTGCCGATGCTGGAACTGGACGAGCACGGCAATGAAACGGTCAGGCGCACGGCCGTTGCCGAGCACTGACGGGTTGCAGGAGCGAGCTTgctcgctc                                                                                                                                         <  8:195246/149‑1 (MQ=255)
          gggACTACTTCATGGATGGCAAGGGGTTGGATGCCGAGGGTAACTAAGTGCCGATTCTGGAACTGCAGGTGCACGGCAATCTAACGGGCAGGCGCACGGCCGTTGCCGAGCACTGACGGGTTTCAGGAGCGAGCTTGCTCGCTCCTGc                                                                                                                                     <  2:158578/148‑1 (MQ=255)
               tactTCAAGGATGGCAAGCTGTTCGATGCCGAGGGTAACGAAGTGCCGATGCTGGAACTGGACGAGCACGGCAATGAAACGGTCAGGCGCACGGCCGTTGCCGAGCACTGACGGGTTGCAGGAGCGAGCTTGCTCGCTCCt                                                                                                                                       <  5:249051/141‑1 (MQ=255)
               tactTCAAGGATGGCAAGCTGTTCGATGCCGAGGGTAACGAAGTGCCGATGCTGGAACTGGACGAGCACGGCAATGAAACGGTCAGGCGCACGGCCGTTGCCGAGCACTGACGGGTTGCAGGAGCGAGCTTGCTCGCTCCt                                                                                                                                       >  6:249051/1‑141 (MQ=255)
                 ctacAAGGATGGCAAGCTGTTCGATGCCGAGGGTAACGAAGTGCCGATGCTGGAACTGGACGAGCACGGCAATGAAACGGTCAGGCGCACGGCCGTTGCCGAGCACTGACGGGTTGCAGGAGCGAGCTTGCTCGCTCCTGCACGCTAt                                                                                                                              >  6:238307/4‑148 (MQ=255)
                  ttCAAGGATGGCAAGCTGTTCGATGCCGAGGGTAACGAAGTGCCGATGCTGGAACTGGACGAGCACGGCAATGAAACGGTCAGGCGCACGGCCGTTGCCGAGCACTGACGGGTTGCAGGAGCGAGCTTGCTCGCTCCTGCACGCTAt                                                                                                                              >  6:114690/1‑147 (MQ=255)
                  ttCAAGGATGGCAAGCTGTTCGATGCCGAGGGTAACGAAGTGCCGATGCTGGAACTGGACGAGCACGGCAATGAAACGGTCAGGCGCACGGCCGTTGCCGAGCACTGACGGGTTGCAGGAGCGAGCTTGCTCGCTCCTGCACGCTAt                                                                                                                              <  5:114690/147‑1 (MQ=255)
                            gCAAGCTGTTCGATGCCGAGGGTAACGAAGTGCCGATGCTGGAACTGGACCAGCACGGCAATGAAACGGTCAGGCGCACGGCCGTTGCCGAGCACTGACGGGTTGCAGGAGCGAGCTTGCTCGCTCCTGCACGCTATTCGAACAATGc                                                                                                                   <  4:10909/148‑1 (MQ=255)
                                gCTGTTCGATGCCGAGGGTAACGAAGTGCCGATGCTGGAACTGGACGAGCACGGCAATGAAACGGTCAGGCGCACGGCCGTTGCCGAGCACTGACGGGTTGCAGGAGCGAGCTTGCTCGCTCCTGCACGCTATTCGAACAAt                                                                                                                     >  7:302601/1‑142 (MQ=255)
                                gCTGTTCGATGCCGAGGGTAACGAAGTGCCGATGCTGGAACTGGACGAGCACGGCAATGAAACGGTCAGGCGCACGGCCGTTGCCGAGCACTGACGGGTTGCAGGAGCGAGCTTGCTCGCTCCTGCACGCTATTCGAACAAt                                                                                                                     <  8:302601/142‑1 (MQ=255)
                                  tGTTCGATGCCGAGGGTAACGAAGTGCCGATGCTGGAACTGGACGAGCACGGCAATGAAACGGTCAGGCGCACGGCCGTTGCCGAGCACTGACGGGTTGCAGGAGCGAGCTTGCTCGCTCCTGCACGCTATTCGAACAATGCATCCAgc                                                                                                            >  6:125992/1‑149 (MQ=255)
                                    ttCGATGCCGAGGGTAACGAAGTGCCGATGCTGGAACTGGACGAGCACGGCAATGAAACGGTCAGGCGCACGGCCGTTGCCGAGCACTGACGGGTTGCAGGAGCGAGCAAGCTCGCTcc                                                                                                                                        <  1:316217/119‑1 (MQ=255)
                                    ttCGATGCCGAGGGTAACGAAGTGCCGATGCTGGAACTGGACGAGCACGGCAATGAAACGGTCAGGCGCACGGCCGTTGCCGAGCACTGACGGGTTGCAGGAGCGAGCAAGCTCGCTcc                                                                                                                                        >  2:316217/1‑119 (MQ=255)
                                           ccGAGGGTAACGAAGTGCCGATGCTGGAACTGGACGAGCACGGCAATGAAACGGTCAGGCGCACGGCCGTTGCCGAGCACTGACGGGTTGCAGGAGCGAGCTTGCTCGCTCCTGCACGCTATTCGAACAATGCATCCAGCGCCTGTTc                                                                                                    >  2:133619/1‑148 (MQ=255)
                                                 gTAACGAAGTGCCGATGCTGGAACTGGACGAGCACGGCAATGAAACGGTCAGGCGCACGGCCGTTGCCGAGCACTGACGGGTTGCAGGAGCGAGCTTGCTCGCTCCTGCACGCTATTCGAACAATGCATCCAGCGCCTGTTCCAGCCgg                                                                                             <  5:238307/149‑1 (MQ=255)
                                                    aCGAAGTGCCGATGCTGGAACTGGACGAGCACGGCAATGAAACGGTCAGGCGCACGGCCGTTGCCGAGCACTGACGGGTTGCAGGAGCGAGCAAGCTCGCTCCt                                                                                                                                       <  5:280543/104‑1 (MQ=255)
                                                    aCGAAGTGCCGATGCTGGAACTGGACGAGCACGGCAATGAAACGGTCAGGCGCACGGCCGTTGCCGAGCACTGACGGGTTGCAGGAGCGAGCAAGCTCGCTCCt                                                                                                                                       >  6:280543/1‑104 (MQ=255)
                                                       aaGTGCCGATGCTGGAACTGGACGAGCACGGCAATGAAACGGTCAGGCGCACGGCCGTTGCCGAGCACTGACGGGTTGCAGGAGCGAGCTTGCTCGCTCCTGCACGCTATTCGAACAATGCATCCAGCGCCTGTTCCAGCCGGGTAACg                                                                                       <  3:246786/149‑1 (MQ=255)
                                                        aGTGCCGATGCTGGAACTGGACGAGCACGGCAATGAAACGGTCAGGCGCACGGCCGTTGCCGAGCACTGACGGGTTGCAGGAGCGAGCTTGCTCGCTCCTGCACGCTATTCGAACAATGCATCCAGCGCCTGTTCCAGCCGGGTAACgg                                                                                      >  8:117706/1‑149 (MQ=255)
                                                            ccGATGCTGGAACTGGACGAGCACGGCAATGAAACGGTCAGGCGCACGGCCGTTGCCGAGCACTGACGGGTTGCAGGAGCGAGCTTGCTCGCTCCTGCACGCTATTCGAACAATGCATCCAGCGCCTGTTCCAGCCGGGTAACGGCAAt                                                                                  >  1:256375/1‑149 (MQ=255)
                                                              gATGCTGGAACTGGACGAGCACGGCAATGAAACGGTCAGGCGCACGGCCGTTGCCGAGCACTGACGGGTTGCAGGAGCGAGCTTGCTCGCTCCTGCACGCTATTCGAACAATGCATCCAGCGCCTGTTCCAGCCGGGTAACGGCAATCa                                                                                >  2:11964/1‑149 (MQ=255)
                                                               aTGCTGGAACTGGACGAGCACGGCAATGAAACGGTCAGGCGCACGGCCGTTGCCGAGCACTGACGGGTTGCAGGAGCGAGCTTGCTCGCTCCTGCACGCTATTCGAACAATGCATCCAGCGCCTGTTCCAGCCGGGTAACGGCAATCAc                                                                               <  7:117706/149‑1 (MQ=255)
                                                                tGCTGGAACTGGACGAGCACGGCAATGAAACGGTCAGGCGCACGGCCGTTGCCGAGCACTGACGGGTTGCAGGAGCGAGCTTGCTCGCTCCTGCACGCTATTCGAACAATGCATCCAgcgc                                                                                                          <  6:229925/121‑1 (MQ=255)
                                                                tGCTGGAACTGGACGAGCACGGCAATGAAACGGTCAGGCGCACGGCCGTTGCCGAGCACTGACGGGTTGCAGGAGCGAGCTTGCTCGCTCCTGCACGCTATTCGAACAATGCATCCAgcgc                                                                                                          >  5:229925/1‑121 (MQ=255)
                                                                    ggAACTGGACGAGCACGGCAATGAAACGGTCAGGCGCACGGCCGTTGCCGAGCACTGACGGGTTGCAGGAGCGAGCTTGCTCGCTCCTGCACGCTATTCGAACAAt                                                                                                                     <  4:178925/106‑1 (MQ=255)
                                                                    ggAACTGGACGAGCACGGCAATGAAACGGTCAGGCGCACGGCCGTTGCCGAGCACTGACGGGTTGCAGGAGCGAGCTTGCTCGCTCCTGCACGCTATTCGAACAAt                                                                                                                     >  3:178925/1‑106 (MQ=255)
                                                                       aCTGGACGAGCACGGCAATGAAACGGTCAGGCGCACGGCCGTTGCCGAGCACTGACGGGTTGCAGGAGCGAGCTTGCTCGCTCCTGCACGCTATTCGAACAATGCATCCAGCGCCTGTTCCAGCCGGGTAACGGCAATCACCTGCAAcc                                                                       >  6:16049/1‑149 (MQ=255)
                                                                        cTGGACGAGCACGGCAATGAAACGGTCAGGCGCACGGCCGTTGCCGAGCACTGACGGGTTGCAGGAGCGAGCTTGCTCGCTCCTGCACGCTATTCGAACAAt                                                                                                                     <  4:129466/102‑1 (MQ=255)
                                                                         tGGACGAGCACGGCAATGAAACGGTCAGGCGCACGGCCGTTGCCGAGCACTGACGGGTTGCAGGAGCGAGCTTGCTCGCTCCTGCACGCTATTCGAACAATGCATCCAGCGCCTGTTCCAGCCGGGTAACGGCAATCACCTGCAAcccc                                                                     >  8:83875/1‑149 (MQ=255)
                                                                           gACGAGCACGGCAATGAAACGGTCAGGCGCACGGCCGTTGCCGAGCACTGACGGGTTGCAGGAGCGAGCTTGCTCGCTCCTGCACGCTATTCGAACAATGCATCCAGCGCCTGTTCCAGCCCGGTAACGGCAATCACCTGCAACCCCg                                                                    >  3:84940/1‑148 (MQ=255)
                                                                                    ggCAATGAAACGGTCAGGCGCACGGCCGTTGCCGAGCACTGACGGGTTGCAGGAGCGAGCTTGCTCGCTCCTGCACGCTATTCGAACAATGCATCCAGCGCCTGTTCCAGCCGGGTAACGGCAATCACCTGCAACCCCGCAGGCGcttc                                                          <  4:84940/149‑1 (MQ=255)
                                                                                      cAATGAAACGGTCAGGCGCACGGCCGTTGCCGAGCACTGACGGGTTGCAGGAGCGAGCTTGCTCGCTCCTGCACGCTATTCGAACAATGCATCCAGCGCCTGTTCCAGCCGGGTAACGGCAATCACCTGCAACCCCGCAGGCGcttctt                                                        >  8:58211/1‑149 (MQ=255)
                                                                                      cAATGAAACGGTCAGGCGCACGGCCGTTGCCGAGCACTGACGGGTTGCAGGAGCGAGCTTGCTCGCTCCTGCACGCTATTCGAACAATGCATCCAGCGCCTGTTCCAGCCGGGTAACGGCAATCACCTGCAACCCCGCAGGCGcttctt                                                        >  2:209773/1‑149 (MQ=255)
                                                                                         tGAAACGGTCAGGCGCACGGCCGTTGCCGAGCACTGACGGGTTGCAGGAGCGAGCTTGCTCGCTCCTGCACGCTATTCGAACAATGCATCCAGCGCCTGTTCCAGCCGGGTAACGGCAATCACCTGCAACCCCGCAGGCGCTTCTTTCg                                                     <  1:11964/149‑1 (MQ=255)
                                                                                           aaaCGGTCAGGCGCACGGCCGTTGCCGAGCACTGACGGGTTGCAGGAGCGAGCTTGCTCGCTCCTGCACGCTATTCGAACAATGCATCCAGCGCCTGTTCCAGCCGGGTAACGGCAATCACCTGCAACCCCGCAGGCGCTTCTTTCGgc                                                   <  5:16049/149‑1 (MQ=255)
                                                                                           aaaCGGTCAGGCGCACGGCCGTTGCCGAGCACTGACGGGTTGCAGGAGCGAGCTTGCTCGCTCCTGCACGCTATTCGAACAATGCATCCAGCGCCTGTTCCAGCCGGGTAACGGCAATCACCTGCAACCCCGCAGGCGCTTCTTTCGgc                                                   >  1:37851/1‑149 (MQ=255)
                                                                                           aaaCGGTCAGGCGCACGGCCGTTGCCGAGCACTGACGGGTTGCAGGAGCGAGCTTGCTCGCTCCTGCACGCTATTCGAACAATGCATCCAGCGCCTGTTCCAGCCGGGTAACGGCAATCACCTGCAACCCCGCAGGCGCTTCTTTCGgc                                                   <  7:83875/149‑1 (MQ=255)
                                                                                                      cgcACGGCCGTTGCCGAGCACTGACGGGTTGCAGGAGCGAGCTTGCTCGCTCCTGCACGCTATTCGAACAATGCATCCAGCGCCTGTTCCAGCCGGGTAACGGCAATCACCTGCAACCCCGCAGGCGCTTCTTTCGGCGCATTGCCCtt                                        <  2:256375/149‑1 (MQ=255)
                                                                                                       gcACGGCCGTTGCCGAGCACTGACGGGTTGCAGGAGCGAGCTTGCTCGCTCCTGCACGCTATTCGAACAATGCATCCAGCGCCTGTTCCAGCCGGGTAACGGCAATCACCTGCAACCCCGCAGGCGCTTCTTTCGGCGCATTGCCCTTg                                       <  1:133619/149‑1 (MQ=255)
                                                                                                        cACGGCCGTTGCCGAGCACTGACGGGTTGCAGGAGCGAGCTTGCTCGCTCCTGCACGCTATTCGAACAATGCATCCAGCGCCTGTTCCAGCCGGGTAACGGCAATCACCTGCAACCCCGCAGGCGCTTCTTTCGGCGCATTGCCCTTgg                                      <  8:202572/149‑1 (MQ=255)
                                                                                                        cACGGCCGTTGCCGAGCACTGACGGGTTGCAGGAGCGAGCTTGCTCGCTCCTGCACGCTATTCGAACAATGCATCCAGCGCCTGTTCCAGCCGGGTAACGGCAATCACCTGCAACCCCGCAGGCGCTTCTTTCGGCGCATTGCCCTTgg                                      <  5:125992/149‑1 (MQ=255)
                                                                                                         aCGGCCGTTGCCGAGCACTGACGGGTTGCAGGAGCGAGCTTGCTCGCTCCTGCACGCTATTCGAACAATGCATCCAGCGCCTGTTCCAGCCGGGTAACGGCAATCACCTGCAACCCCGCAGGCGCTTCTTTCGGCGCATTGCCCTTggg                                     >  3:290606/1‑149 (MQ=255)
                                                                                                         aCGGCCGTTGCCGAGCACTGACGGGTTGCAGGAGCGAGCTTGCTCGCTCCTGCACGCTATTCGAACAATGCATCCAGCGCCTGTTCCAGCCGGGTAACGGCAATCACCTGCAACCCCGCAGGCGCTTCTTTCGGCGCATTGCCCTTggg                                     >  1:143250/1‑149 (MQ=255)
                                                                                                               gTTGCCGAGCACTGACGGGTTGCAGGAGCGAGCTTGCTCGCTCCTGCACGCTATTCGAACAATGCATCCAGCGCCTGTTCCAGCCGGGTAACGGCAATCACCTGCAACCCCGCAGGCGCTTCTTTCGGCGCATTGCCCTTGGGCACGAt                               <  2:143250/149‑1 (MQ=255)
                                                                                                               gTTGCCGAGCACTGACGGGTTGCAGGAGCGAGCTTGCTCGCTCCTGCACGCTATTCGAACAATGCATCCAGCGCCTGTTCCAGCCGGGTAACGGCAATCACCTGCAACCCCGCAGGCGCTTCTTTCGGCGCATTGCCCTTGGGCACGAt                               <  4:290606/149‑1 (MQ=255)
                                                                                                                   ccGAGCACTGACGGGTTGCAGGAGCGAGCTTGCTCGCTCCTGCACGCTATTCGAACAATGCATCCAGCGCCTGTTCCAGCCGGGTAACGGCAATCACCTGCAACCCCGCAGGCGCTTCTTTCGGCGCATTGCCCTTGGGCACGATGGc                            <  2:37851/148‑1 (MQ=255)
                                                                                                                          cTGACGGGTTGCAGGAGCGAGCTTGCTCGCTCCTGCACGCTATTCGAACAATGCATCCAGCGCCTGTTCCAGCCGGGTAACGGCAATCACCTGCAACCCCGCAGGCGCTTCTTTCGGCGCATTGCCCTTGGGCACGATg                              >  3:298927/1‑139 (MQ=255)
                                                                                                                          cTGACGGGTTGCAGGAGCGAGCTTGCTCGCTCCTGCACGCTATTCGAACAATGCATCCAGCGCCTGTTCCAGCCGGGTAACGGCAATCACCTGCAACCCCGCAGGCGCTTCTTTCGGCGCATTGCCCTTGGGCACGATg                              <  4:298927/139‑1 (MQ=255)
                                                                                                                                       ggAGCGAGCTTGCTCGCTCCTGCACGCTATTCGAACAATGCATCCAGCGCCTGTTCCAGCCGGGTAACGGCAATCACCTGCAACCCCGCAGGCGCTTCTTTCGGCGCATTGCCCTTGGGCACGATGGCACGCTTGAAGCCGTGCTTGGc       <  1:301712/149‑1 (MQ=255)
                                                                                                                                             agcTTGCTCGCTCCTGCACGCTATTCGAACAATGCATCCAGCGCCTGTTCCAGCCGGGTAACGGCAATCACCTGCAACCCCGCAGGCGCTTCTTTCGGCGCATTGCCCTTGGGCACGATGGCACGCTTGAAGCCGTGCTTGGCCGctt  >  7:310911/1‑148 (MQ=255)
                                                                                                                                                 |                                                                                                                                               
CTGCGTGATGGGGACTACTTCAAGGATGGCAAGCTGTTCGATGCCGAGGGTAACGAAGTGCCGATGCTGGAACTGGACGAGCACGGCAATGAAACGGTCAGGCGCACGGCCGTTGCCGAGCACTGACGGGTTGCAGGAGCGAGCTTGCTCGCTCCTGCACGCTATTCGAACAATGCATCCAGCGCCTGTTCCAGCCGGGTAACGGCAATCACCTGCAACCCCGCAGGCGCTTCTTTCGGCGCATTGCCCTTGGGCACGATGGCACGCTTGAAGCCGTGCTTGGCCGCTT  >  NC_002947/5268807‑5269095

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