Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 1,921,550 T→A 11.9% intergenic (+25/+36) PP_1720 → / ← PP_1721 alcohol dehydrogenase/phosphoserine phosphatase

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0029471,921,5500TA11.9% 115.4 / 6.1 50intergenic (+25/+36)PP_1720/PP_1721alcohol dehydrogenase/phosphoserine phosphatase
Reads supporting (aligned to +/- strand):  ref base T (23/21);  new base A (2/4);  total (25/25)
Fisher's exact test for biased strand distribution p-value = 6.67e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 8.56e-01

GGTTGTCGCCGCAGCTGGTCGACACCTATCCGGTGGAGTTTGCCCAGGCGGCGTATGCCGAGCTGGAGACCAACCAGGTGTCGGGCAAGCTGGTGATGGTGATCGACCCTAGCCTGGTGTAAGCGGGGCTGGCCCCTTCGCGGGTTTACCCGCGAAGGGGCCGGCACAGCCCATAAAAACCTTCAGCTCCAGCTCAGAGTCGGCCACCCATTGATCTCGGCATGCTCGCGCAACACTGCGTCCGGGTTCACCACATGCGGGTGATCCACTCTCAACAACAACGGCAAGT  >  NC_002947/1921404‑1921692
                                                                                                                                                  |                                                                                                                                              
ggTTGTCGCCGCAGCTGGTCGACACCTATCCGGTGGAGTTTGCCCAGGCGGCGTATGCCGAGCTGGAGACCAACCAGGTGTCGGGCAAGCTGGTGATGGTGATCGACCCTAGCCTGGTGTAAGCGGGGCTGGCCCCTTCGCGGGTTTAc                                                                                                                                              >  8:216435/1‑149 (MQ=255)
     tCGCCGCAGCTGGTCGACACCTATCCGGTGGAGTTTGCCCAGGCGGCGTATGCCGAGCTGGAGACCAACCAGGTGTCGGGCAAGCTGGTGATGGTGATCGACCCTAGCCTGGTGTAAGCGGGGCTGGCCCCTTCGCGGGTTTACCcgc                                                                                                                                          <  2:204946/148‑1 (MQ=255)
     tCGCCGCAGCTGGTCGACACCTATCCGGTGGAGTTTGCCCAGGCGGCGTATGCCGAGCTGGAGACCAACCAGGTGTCGGGCAAGCTGGTGATGGTGATCGACCCTAGCCTGGTGTAAGCGGGGCTGGCCCCTTCACGGGTTTACCcgc                                                                                                                                          >  1:356967/1‑148 (MQ=255)
       gccgcAGCTGGTCGACACCTATCCGGTGGAGTTTGCCCAGGCGGCGTATGCCGAGCTGGAGACCAACCAGGTGTCGGGCAAGCTGGTGATGGTGATCGACCCTAGCCTGGTGTAAGCGGGGCTGGCCCCTTCGCGGGTAAACCCGCGaa                                                                                                                                       <  4:208478/149‑1 (MQ=255)
         cgcAGCTGGTCGACACCTATCCGGTGGAGTTTGCCCAGGCGGCGTATGCCGAGCTGGAGACCAACCAGGTGTCGGGCAAGCTGGTGATGGTGATCGACCCTAGCCTGGTGTAAGCGGGGCTGGCCCCTTCGCGGGTTTACCCGCGAAgg                                                                                                                                     >  1:80106/1‑149 (MQ=255)
          gcAGCTGGTCGACACCTATCCGGTGGAGTTTGCCCAGGCGGCGTATGCCGAGCTGGAGACCAACCAGGTGTCGGGCAAGCTGGTGATGGTGATCGACCCTAGCCTGGTGTAAGCGGGGCTGGCCCCTTCGCGGGTAAACCCGCGAAggg                                                                                                                                    <  3:61127/149‑1 (MQ=255)
          gcAGCTGGTCGACACCTATCCGGTGGAGTTTGCCCAGGCGGCGTACGCCGAGCTGGAGACCAACCAGGTGTCGGGCAAGCTGGTGATGGTGATCGACCCTAGCCTGGTGTAAGCGGGGCTGGCCCCTTCGCGGGTTTACCCGCGAAggg                                                                                                                                    >  8:117321/1‑149 (MQ=255)
             gCTGGTCGACACCTATCCGGTGGAGTTTGCCCAGGCGGCGTATGCCGAGCTGGAGACCAACCAGGTGTCGGGCAAGCTGGTGATGGTGATCGACCCTAGCCTGGTGTAAGCGGGGCTGGCCCCTTCGCGGGTTTACCCGCGAAGGGGcc                                                                                                                                 >  3:117484/1‑149 (MQ=255)
                  tCGACACCTATCCGGTGGAGTTTGCCCAGGCGGCGTATGCCGAGCTGGAGACCAACCAGGTGGCGGGCAAGCTGGTGATGGTGAGCGACCCTAGCCTGGTGTAAGCGGGGCTGGCCCCTTCGCGGGTTTACCCGCGAAGGGGCCGGcac                                                                                                                            <  7:277571/149‑1 (MQ=255)
                   cGACACCTATCCGGTGGAGTTTGCCCAGGCGGCGTATGCCGAGCTGGAGACCAACCAGGTGTCGGGCAAGCTGGTGATGGTGATCGACCCTAGCCTGGTGTAAGCGGGGCTGGCCCCTTCGCGGGTTTACCCGCGAAGGGGCCGGcaca                                                                                                                           <  4:279485/149‑1 (MQ=255)
                       acCTATCCGGTGGAGTTTGCCCAGGCGGCGTATGCCGAGCTGGAGACCAACCAGGTGGCGGGCAAGCTGGTGATGGTGATCGACCCTAGCCTGGTGTAAGCGGGGCTGGCCCCTTCGCGGGTTTACCCGCGAAGGGGCCGGCACAGccc                                                                                                                       <  2:131107/149‑1 (MQ=255)
                         cTATCCGGTGGAGTTTGCCCAGGCGGCGTATGCCGAGCTGGAGACCAACCAGGTGTCGGGCAAGCTGGTGATGGTGATCGACCCTAGCCTGGTGTAAGCGGGGCTGGCCCCTTCGCGGGTTTACCCGCGAAGGGGCCGGCACAGCCCa                                                                                                                      >  7:278206/1‑148 (MQ=255)
                         cTATCCGGTGGAGTTTGCCCAGGCGGCGTATGCCGAGCTGGAGACCAACCAGGTGTCGGGCAAGCTGGTGATGGTGATCGACCCTAGCCTGGTGTAAGCGGGGCTGGCCCCTTCGCGGGTTTACCCGCGAAGGGGCCGGCACAGCCCa                                                                                                                      >  3:31549/1‑148 (MQ=255)
                            tCCGGTGGAGTTTGCCCAGGCGGCGTATGCCGAGCTGGAGACCAACCAGGTGTCGGGCAAGCTGGTGATGGTGATCCACCCTAGCCTGGTGTAAGCGGGGCTGGCCCCTTCGCTGGTTTACCCGCGAAGGGGCCGGCACAGCCCATaa                                                                                                                   <  1:53388/148‑1 (MQ=255)
                                   gAGTTTGCCCAGGCGGCGTATGCCGAGCTGGAGACCAACCAGGTGTCGGGCAAGCTGGTGATTGTGATCGACCCTAGCCTGGTGTAAGCGGGGCTGGCCCCTTCGCGGGTTTACCCGCGAAGGGGCCGGCACAGCCCATAAAAACCTTc                                                                                                           <  5:313128/149‑1 (MQ=255)
                                          cccAGGTGGCGTATGCGTAGGTGTAGACAAACATGGTGTGGGGAAAGTTGGTGATGTTGATCGACCGTAGCTTGGTGTATGCGGGGCTGGCCCCTTCGCGGGTTTACCCGCGAAGGGGCCGGCACAGCCCATAAAAACCTTcagctcca                                                                                                    <  2:127976/149‑1 (MQ=255)
                                          cccAGGCGGCGTATGCCGAGCTGGAGACCAACCAGGTGTCGGGCAAGCTGGTGATGGTGATCGACCCTAGCCTGGTGTAAGCGGGGCTGGCCCCTTCGCGGGTTTACCCGCGAAGGGGCCGGCACAGCCCATAAAAACCTTcagctc                                                                                                      >  3:175596/1‑147 (MQ=255)
                                           ccAGGCGGCGTATGCCGAGCTGGAGACCAACCAGGTGTCGGGCAAGCTGGTGATGGTGATCGACCCTAGCCTGGTGTAAGCGGGGCTGGCCCCTTCGCGGGTTTACCCGCGAAGGGGCCGGCACAGCCCATAAAAACCTTcagctccag                                                                                                   >  3:218778/1‑149 (MQ=255)
                                                cggcgTATGCCGAGCGGGCGACAAACCTGGGGTGGGGCAAGCTGGTTATGGTGTGCGCCCTTTGCGGGGTTTACCGGGGGGGGGCCCGTGCGGGGGTTTACCCGGGACGGGGCCGGCACAGCCCATAAAAACCTTCAGCTCCAGCTCag                                                                                              <  4:177630/149‑1 (MQ=255)
                                                       tGCCGAGCTGGAGACCAACCAGGTGTCGGGCAAGCTGGTGATGGTGATCGACCCTAGCCTGGTGTAAGCGGGGCTGGCCCCTTCGCGGGTTTACCCGCGAAGGGGCCGGCACAGCCCATAAAAACCTTCAGCTCCAGCTCAGAGTCg                                                                                         >  6:220655/1‑147 (MQ=255)
                                                       tGCCGAGCTGGAGACCAACCAGGTGTCGGGCAAGCTGGTGATGGTGATCGACCCTAGCCTGGTGTAAGCGGGGCTGGCCCCTTCGCGGGTTTACCCGCGAAGGGGCCGGCACAGCCCATAAAAACCTTCAGCTCCAGCTCAGAGTCg                                                                                         <  5:220655/147‑1 (MQ=255)
                                                          cGAGCTGGAGACCAACCAGGTGTCGGGCAAGCTGGTGATGGTGATCGACCCTAGCCTGGTGTAAGCGGGGCTGGCCCCTTCGCGGGTTTACCCGCGAAGGGGCCGGCACAGCCCATAAAAACCTTCAGCTCCAGCTCAGAGTCGGCCAc                                                                                    >  7:223081/1‑149 (MQ=255)
                                                           gAGCTGGAGACCAACCAGGTGTCGGGCAAGCTGGTGATGGTGATCGACCCTAGCCTGGTGTAAGCGGGGCTGGCCCCTTCGCGGGTTTACCCGCGAAGGGGCCGGCACAGCCCATAAAAACCTTCAGCTCCAGCTCAGAGTCGGCCAcc                                                                                   <  7:117321/149‑1 (MQ=255)
                                                              cTGGAGACCAACCAGGTGTCGGGCAAGCTGGTGATGGTGATCGACCCTAGCCTGGTGTAAGCGGGGCTGGCCCCTTCGCGGGTTTACCCGCGAAGGGGCCGGCACAGCCCATAAAAACCTTCAGCTCCAGCTCAGAGTCGGCCACCCa                                                                                 <  4:117484/148‑1 (MQ=255)
                                                                        accaGGTGTCGGGCAAGCTGGTGATGGTGATCGACCCTAGCCTGGTGTAAGCGGGGCTGGCCCCTTCGCGGGTTTACCCGCGAAGGGGCCGGCACAGCCCATAAAAACCTTCAGCTCCAGCTCAGAGTCGGCCACCCATTGATCTCgg                                                                       >  5:283232/1‑148 (MQ=255)
                                                                        accaGGTGTCGGGCAAGCTGGTGATGGTGATCGACCCTAGCCTGGTGTAAGCGGGGCTGGCCCCTTCGCGGGTTTACCCGCGAAGGGGCCGGCACAGCCCATAAAAACCATCAGCTCCAGCTCAGAGTCGGCCACCCATTGATCTCgg                                                                       >  2:258138/1‑148 (MQ=255)
                                                                         ccaGGTGTGGGGCAAGGTGGTGATGGTGATCGAGCCTAGGCTTGTGTAAGGGGGTCGGGCCGGTTGGGGGGTTTACCCGCGAAGGGGCCGGCACAGACCATAATATCCTTCAGCTCCAGCTCAGAGTCGGCCACCCATTGATCTCgg                                                                       <  2:80106/147‑1 (MQ=255)
                                                                            ggTGTCGGGCAAGCTGGTGATGGTGATCGACCCTAGCCTGGTGTAAGCGGGGCTGGCCCCTTCGCGGGTTTACCCGCGAAGGGGCCGGCACAGCCCATAAAAACCTTCAGCTCCAGCTCAGAGTCGGCCACCCATTGATCTCGGCATGc                                                                  >  2:249764/1‑149 (MQ=255)
                                                                            ggTGTCGGGCAAGCTGGTGATGGTGATCGACCCTAGCCTGGTGTAAGCGGGGCTGGCCCCTTCGCGGGTTTACCCGCGAAGGGGCCGGCACAGCCCATAAAAACCTTCAGCTCCAGCTCAGAGTCGGCCACCCATTGATCTCGGCATGc                                                                  <  6:283232/149‑1 (MQ=255)
                                                                                     cAAGCTGGTGATGGTGATCGACCCTAGCCTGGTGTAAGCGGGGCTGGCCCCTTCGCGGGTTTACCCGCGAAGGGGCCGGCACAGCCCATAAAAACCTTCCGCTCCAGCTCAGAGTCGGCCACCCATTGATGTCGGCATGCTCGCGCaa                                                          >  4:379965/1‑148 (MQ=255)
                                                                                     cAAGCTGGTGATGGTGATCGACCCTAGCCTGGTGTAAGCGGGGCTGGCCCCTTCGCGGGTTTACCCGCGAAGGGGCCGGCACAGCCCATAAAAACCTTCAGCTCCAGCTCAGAGTCGGCCCCCCAATGATCTCGGCATGCTAGCCCaa                                                          >  4:206920/1‑148 (MQ=255)
                                                                                        gCTGGTGATGGTGATCGACCCTAGCCTGGTGTAAGCGGGGCTGGCCCCTTCGCGGGTTTACCCGCGAAGGGGCCGGCACAGCCCATAAAAACCTTCAGCTCCAGCTCAGAGTCGGCCACCCATTGATCTCGGCATGCTCGCGCAACACt                                                      <  8:223081/149‑1 (MQ=255)
                                                                                             tgatggtgatCGACCCTAGCCTGGTGTAAGCGGGGCTGGCCCCTTCACGGGTTTACCCGCGAAGGGGCCGGCACAGCCCATAAAAACCTTCAGCTCCAGCTCAGAGTCGGCCACCCATTGATCTCGGCATGCTCGCGCAACACTGCGTc                                                 <  2:356967/149‑1 (MQ=255)
                                                                                                  gtgatCGCCCCTAGCCGGGTGTAAGCGGGGCTGGCCCGTTCGGGGGTTAACCGGCGAAGGGGCCGGCACAGCCCATAAAAACCTTCAGCTCCAGCTCAGAGTCGGCCACCCATTGATCTCGGCATGCTCGCGCAACACTGCGTCCggg                                             <  3:206920/148‑1 (MQ=255)
                                                                                                      tCGACCCTAGCCTGGTGTAAGCGGGGCTGGCCCCTTCGCGGGTTTACCCGCGAAGGGGCCGGCACAGCCCATAAAAACCTTCAGCTCCAGCTCAGAGTCGGCCACCCATTGATCTCGGCATGCTCGCGCAACACTGCGTCCGGGTTcac                                        <  4:218778/149‑1 (MQ=255)
                                                                                                         aCCCTAGCCTGGTGTAAGCGGGGCTGGCCCCTTCGCGGGTTTACCCGCGAAGGGGCCGGCACAGCCCATAAAAACCTTCAGCTCCAGCTCAGAGTCGGCCACCCATTGATCTCGGCATGCTCGCGCAACACTGCGTCCGGGTTcaccac                                     <  6:166396/149‑1 (MQ=255)
                                                                                                         aCCCTAGCCTGGTGTAAGCGGGGCTGGCCCCTTCGCGGGTTTACCCGCGAAGGGGCCGGCACAGCCCATAAAAACCTTCAGCTCCAGCTCAGAGTCGGCCACCCATTGATCTCGGCATGCTCGCGCAACACTGCGTCCGGGTTcacca                                      >  5:166396/1‑148 (MQ=255)
                                                                                                            cTAGCCTGGTGTAAGCGGGGCTGGCCCCTTCGCGGGTTTACCCGCGAAGGGGCCGGCACAGCCCATAAAAACCTTCAGCTCCAGCTCAGAGTCGGCCACCCATTGATCTCGGCATGCTCGCGCAACACTGCGTCCGGGTTCACCACATg                                  >  3:341562/1‑149 (MQ=255)
                                                                                                               gCCTGGTGTAAGCGGGGCTGGCCCCTTCGCGGGTTTACCCGCGAAGGGGCCGGCACAGCCCATAAAAACCTTCAGCTCCAGCTCAGAGTCGGCCACCCATTGATCTCGGCATGCTCGCGCAACACTGCGTCCGGGTTCACCACATGCgg                               <  8:278206/149‑1 (MQ=255)
                                                                                                                ccTGGTGTAAGCGGGGCTGGCCCCTTCGCGGGTTTACCCGCGAAGGGGCCGGCACAGCCCATAAAAACCTTCAGCTCCAGCTCAGAGTCGGCCACCCATTGATCTCGGCATGCTCGCGCAACACTGCGTCCGGGTTCACCACATGCggg                              <  7:216435/149‑1 (MQ=255)
                                                                                                                   gtggtATGCGGGGCTGGCCCCTGCGCGGGGTTTCCCGGGATGGGGCCGGCACAGCCCATAAAAACCTTCAGCTCCAGCTCAGAGTCGGCCACCCATTGATCTCGGCATGCTCGCGCAACACTGCGTCCGGGTTCACCACATGCGGGTGa                           <  4:31549/146‑1 (MQ=255)
                                                                                                                     tgtAAGCGGGGCTGGCCCCTTCGCGGGTTTACCCGCGAAGGGGCCGGCACAGCCCATAAAAACCTTCAGCTCCAGCTCAGAGTCGGCCACCCATTGATCTCGGCATGCTCGCGCCACACTGCGTCCGGGTTCACCACATGCGGGTGATc                         >  6:145880/1‑149 (MQ=255)
                                                                                                                                     cccTTCGCGGGTTTACCCGCGAAGGGGCCGGCACAGCCCATAAAAACCTTCAGCTCCAGCGCAGAGTCGGCCACCCATTGATCTCGGCATGCTCGCGCAACACTGCGTCCGGGTTCACCACATGCGGGTGATCCACTCTCaacagcaac         >  8:249301/1‑149 (MQ=255)
                                                                                                                                      ccTTCGCGGGTTTACCCGCGAAGGGGCCGGCACAACCCATAAAAACCTTCAGCTCCAGCTCAGAGTCGGCCACCCATTGATCTCGGCATGCTCGCGCAACACTGCGTCCGGGTTcacc                                       >  7:320561/1‑118 (MQ=255)
                                                                                                                                      ccTTCGCGGGTTTACCCGCGAAGGGGCCGGCACAACCCATAAAAACCTTCAGCTCCAGCTCAGAGTCGGCCACCCATTGATCTCGGCATGCTCGCGCAACACTGCGTCCGGGTTcacc                                       <  8:320561/118‑1 (MQ=255)
                                                                                                                                      ccTTCGCGGGTAAGCCCGCGAAGGGGCCGGCACAGCCCATAAAAACCTTCAGCTCCAGCTCAGAGTCGGCCCCCCATTGATCTCGGCATGCTCGCGCAACACTGCGTCCGGGTTCACCACATGCGGGTGATCCACTCTCAACAACAACg        >  2:191490/1‑149 (MQ=255)
                                                                                                                                       ctccgcgGGTTACCCCGCGACGGGGCCGGCACAGCCCATAAAAACCTTCAGCTCCAGCTCAGAGTCGGCCACCCATTGATCTCGGCATGCTCGCGCAACACTGCGTCCGGGTTCACCACATGCGGGTGATCCACTCTCAACAACAACgg       <  4:175596/146‑1 (MQ=255)
                                                                                                                                        ttCGCGGGTAAACCCGCGAAGGGGCCGGCACAGCCCATAAAAACCTTCAGCTCCAGCTCAGAGTCGGCCACCCATTGATCTCGGCATGCTCGCGCAACACTGCGTCCGGGTTCACCACATGCGGGTGATCCCCTCTACACAACAACgg       >  4:26488/1‑148 (MQ=255)
                                                                                                                                            cgGGTTTACCCGCGAAGGGGCCGGCACAGCCCATAAAAACCTTCAGATCCAGCTCAGAGACGGACACCCATTGATCTCGGCATGATCCAGCAACACTGCGTCCGGGTTCAGCACATGCGGGAGCTCCACTCTCAACAACAAAGGCAAg   >  4:336358/1‑148 (MQ=255)
                                                                                                                                             gggTTTACCCGCGAAGGGGCCGGCACAGCCCATAAAAACCTTCAGCTCCAGCTCAGAGTCGGCCACCCATTGATATCGGCATGCTCGCGCAACACTGCGTCCGGGTTCACCACATGCTGGTGATCCACTCTCAACAACAACGGCAAGt  >  8:113367/1‑148 (MQ=255)
                                                                                                                                                  |                                                                                                                                              
GGTTGTCGCCGCAGCTGGTCGACACCTATCCGGTGGAGTTTGCCCAGGCGGCGTATGCCGAGCTGGAGACCAACCAGGTGTCGGGCAAGCTGGTGATGGTGATCGACCCTAGCCTGGTGTAAGCGGGGCTGGCCCCTTCGCGGGTTTACCCGCGAAGGGGCCGGCACAGCCCATAAAAACCTTCAGCTCCAGCTCAGAGTCGGCCACCCATTGATCTCGGCATGCTCGCGCAACACTGCGTCCGGGTTCACCACATGCGGGTGATCCACTCTCAACAACAACGGCAAGT  >  NC_002947/1921404‑1921692

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: