Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 1,921,550 T→A 24.7% 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,5500TA24.7% 68.6 / 20.8 47intergenic (+25/+36)PP_1720/PP_1721alcohol dehydrogenase/phosphoserine phosphatase
Reads supporting (aligned to +/- strand):  ref base T (20/15);  new base A (5/6);  total (25/22)
Fisher's exact test for biased strand distribution p-value = 7.30e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.64e-01

GCGGTTGTCGCCGCAGCTGGTCGACACCTATCCGGTGGAGTTT‑GCCCAGGCGGCGTATGCCGAGCTGGAGACCAACCAGGTGTCGGGCAAGCTGGTGATGGTGATCGACCCTAGCCTGGTGTAAGCGGGGCTGGCCCCTTCGCGGGTTTACCCGCGAAGGGGCCGGCACAGCCCATAAAAACCTTCAGCTCCAGCTCAGAGTCGGCCACCCATTGATCTCGGCATGCTCGCGCAACACTGCGTCCGGGTTCACCACATGCGGGTGATCCACTCTCAACAACAACGGCAAGTCATTC  >  NC_002947/1921402‑1921697
                                                                                                                                                     |                                                                                                                                                   
gCGGTTGTCGCCGCAGCTGGTCGACACCTATCCGGTGGAGTTT‑GCCCAGGCGGCGTATGCCGAGCTGGAGACCAACCAGGTGTCGGGCAAGCTGGTGATGGTGATCGACCCTAGCCTGGTGTAAGCGGGGCTGGCCCCTTCGCGGGttt                                                                                                                                                     >  1:335987/1‑149 (MQ=255)
       tCGCCGCAGCTGGTCGACACCTATCCGGTGGAGTTT‑GCCCAGGCGGCGTATGCCGAGCTGGAGACCAACCAGGTGTCGGGCAAGCTGGTGATGGTGATCGACCCTAGCCTGGTGTAAGCGGGGCTGGCCCCTTCGCGGGTTTACCcgc                                                                                                                                               <  8:114651/148‑1 (MQ=255)
        cgccgcAGCTGGTCGACACCTATCCGGTGGAGTTTGGCCCAGGCGGCGTATGCCGAGCTGGAGACCAACCAGGTGTCGGGCAAGCTGGTGATGGTGATCGACCCTAGCCTGGTGTAAGCGGGGCTGGCCCCTTCGCGGGTTTACCcgc                                                                                                                                               >  1:240426/1‑148 (MQ=255)
         gccgcAGCTGGTCGACACCTATCCGGTGGAGTTT‑GCCCAGGCGGCGTATGCCGAGCTGGAGACCAACCAGGTGTCGGGCAAGCTGGTGATGGTGATCGACCCTAGCCTGGTGTAAGCGGGGCTGGCCCCTTCGCGGGTTTACCcgcg                                                                                                                                              >  1:95943/1‑147 (MQ=255)
         gccgcAGCTGGTCGACACCTATCCGGTGGAGTTT‑GCCCAGGCGGCGTATGCCGAGCTGGAGACCAACCAGGTGTCGGGCAAGCTGGTGATGGTGATCGACCCTAGCCTGGTGTAAGCGGGGCTGGCCCCTTCGCGGGTTTACCCGCGaa                                                                                                                                            <  8:249672/149‑1 (MQ=255)
         gccgcAGCTGGGCGACACCTATCCGGTGGAGTTT‑GCCCAGGCGGCGTATGCCGAGCTGGAGACCAACCAGGTGTCGGGCAAGCTGGTGATGGTGATCGACCCTAGCCTGGTGTAAGCGGGGCTGGCCCCTTCGCGGGTTTACCCGCGaa                                                                                                                                            <  7:239743/149‑1 (MQ=255)
           cgcAGCTGGTCGACACCTATCCGGTGGAGTTT‑GCCCAGGCGGCGTATGCCGAGCTGGAGACCAACCAGGTGTCGGGCAAGCTGGTGATGGTGATCGACCCTAGCCTGGTGTAAGCGGGNCTGGCCCCTTCGCGGGTAAACCCGCGAAgg                                                                                                                                          <  1:33075/149‑1 (MQ=255)
           cgcAGCTGGTCGACACCTATCCGGTGGAGTTT‑GCCCAGGCGGCGTATGCCGAGCTGGAGACCAACCAGGTGTCGGGCAAGCTGGTGATGGTGATCGACCCTAGCCTGGTGTAAGCGGGGCTGGCCCCTTCGCGGGGTTACCCGCGAAgg                                                                                                                                          >  5:243620/1‑149 (MQ=255)
            gcAGCTGGTCGACACCTATCCGGTGGAGTTT‑GCCCAGGCGGCGTATGCCGAGCTGGAGACCAACCAGGTGTCGGGCAAGCTGGTGATGGTGATCGACCCTAGCCTGGTGTAAGCGGGGCCGGCCCCTTCGCGGGTAAACCCGCGAAggg                                                                                                                                         <  1:218979/149‑1 (MQ=255)
                  ggTCGACACCTATCCGGTGGAGTTT‑GCCCAGGCGGCGTATGCCGAGCCGGAGACCAACCAGGTGTCGGGCAAGCTGGTGATGGTGATCGACCCTAGCCTGGTGTAAGCGGGGCTGGCCCCTTCGCGGGTTTACCCGCGAAGGGGCCgg                                                                                                                                    >  5:165824/1‑148 (MQ=255)
                                 ggTGGAGTTT‑GCCCAGGCGGCGTATGCCGAGCTGGAGACCAACCAGGTGTCGGGCAAGCTGGTGATGGTGATCGACCCTAGCCTGGTGTAAGCGGGGCTGGCCCCTTCGCGGGTTTACCCGCGAAGGGGCCGGCCCAGCCCATAAAAAc                                                                                                                    >  6:221366/1‑149 (MQ=255)
                                  ggggATTTT‑GCCCTGGCGGCGTATGCCGTGCTGTAGACCAACCAGGTGTCGGGCAAGCTGGTGATGGTGATCGACCGTTGCCTGGTGTAAGCGGGGCTGGCCCCTTGGCGGGTTTACCGGCGAAGGGGCCGGCACAGCCCATAAAAAcc                                                                                                                   <  2:116067/147‑1 (MQ=255)
                                      aGTTT‑GCCCTGGCGGCGTATTCCGAGCTGGTGACCTAGCTGGTGTCGGGGAAGTGGGTGATGTTGATCGACGATTGCGTGGTGTTTGGGGGGCTGGACGGTTCGCGGGTTTACCCGCGATGGTGCCGGCACAGCCCATAAAAACCTTc                                                                                                                <  2:240426/148‑1 (MQ=255)
                                      aGTTT‑GCCCAGGCGGCGTATGCCGAGCTGGAGACCAACCAGGTGTCGGGCAAGCTGGTGATGGTGATCGACCCTAGCCTGGTGTAAGCGGGGCTGGCCCCTTCGCGGGTTTACCCGCGAAGGGGCCGGCACAGCCCATAAAAACCtt                                                                                                                 >  8:202481/1‑147 (MQ=255)
                                            gCCCAGGCGGCGTATGCCGAGCTGGAGACCAACCAGGTGTCGGGCAAGCTGGTGATGGTGATCGACCCTAGCCTGGTGTAAGCGGGGCTGGCCCCTTCGCGGGTTTACCCGCGAAGGGGCCGGCACAGCCCATAAAAACCTTcagctc                                                                                                           >  7:131916/1‑148 (MQ=255)
                                                           gCCGAGCTGGAGACCAACCAGGTGTCGGGCAAGCTGGTGATGGTGATCGACCCTAGCCTGGTGTAAGCGGGGCTGGCCCCTTCGCGGGTTTACCCGCGAAGGGGCCGGCACAGCCCATAAAAACCTTCAGCTCCAGCTCAGAGTCGGcc                                                                                           >  4:49448/1‑149 (MQ=255)
                                                                 cTGGAGACCAACCAGGTGTCGGGCAAGCTGGTGATGGTGATCGACCCTAGCCTGGTGTAAGCGGGGCTGGCCCCTTCGCGGGTTTACCCGCGAAGGGGCCGGCACAGCCCATAAAAACCTTCAGCTCCAGCTCAGAGTCGGCCACCCAt                                                                                     <  7:202481/149‑1 (MQ=255)
                                                                    tttacgaacctGGTGGGGGGCACGGTGGTGATGTTTATCGACCTTACCCTGGTTTAAGCGGGTCTGGCCGCTTCGCGTTTTTACCCGCTAAGGGGCCGGCACAGCCCATAAAAACCTTCAGCTCCAGCTCAGAGTCGGCCACCCATTg                                                                                   <  1:309149/142‑1 (MQ=255)
                                                                      gaCCAACCAGGTGTCGGGCAAGCTGGTGATGGTGATCGACCCTAGCCTGGTGTAAGCGGGGCTGGCCCCTTCGCGGGTTTACCCGCGAAGGGGCCGGCACAGCCCATAAAAACCTTCAGCTACAGCTCAGAGTCGGCCCCCCAATGAt                                                                                 >  4:1874/1‑148 (MQ=255)
                                                                       accaaccaGGTGTCGGGCAAGCTGGAGAAGTTGATCGACCCTACCCTGTTCTAAGCGGGCCTGGCCCCTTCGCGCCTTTACCCACGAACGGGCCCACACAGCCCATAAAAACCTTCACCACCACCAAATATTCGCACACCCAttaaac                                                                                >  2:309149/1‑144 (MQ=255)
                                                                         caaccaGGTGTCGGGCAAGCTGGTGATGGTGATCGACCCTAGCCTGGTGTAAGCGGGGCTGGCCCCTTCGCGGGTTTACCCGCGAAGGGGCCGGCACAGCCCATAAAAACCTTCAGCTCCAGCTCAGAGTCGGCCACCCATTGATCTCg                                                                             >  3:48052/1‑149 (MQ=255)
                                                                         caaccaGGTGTCGGGCAAGCTGGTGATGGTGATCGACCCTAGCCTGGTGTAAGCGGGGCTGGCCCCTTCGCGGGTTTACCCGCGAAGGGGCCGGCACAGCCCATAAAAACCTTCAGCTCCAGCTCAGAGTCGGCCACCCATTGATCTCg                                                                             >  5:289669/1‑149 (MQ=255)
                                                                                     ggccAAGGTGGTGTTGGTGATCCCCCCTAGCCGGGTGTACCCGGGGCTGGCCCCTCCGCGGGTTCACCCGCGAAGGGGCCGGCACAGCCCATAAAAACCTTCAGCTCCAGCTCAGAGTCGGCCACCCATTGATCTCGGCATGCTCgcgc                                                                 <  8:131916/146‑1 (MQ=255)
                                                                                                      tgatCGACCCTAGCCTGGTGTAAGCGGGGCTGGCCCCTTCGCGGGTTTACCCGCGAAGGGGCCGGCACAGCCCATAAAAACCTTCAGCTCCAGCTCAGAGTAGGCCACCCATTGATCTCGGCATGCTCGCGCAACACTGCGTCCGGGtt                                                >  1:329068/1‑149 (MQ=255)
                                                                                                        atCGACCCTAGCCTGGTGTAAGCGGGGCTGGCCCCTTCGCGGGTTTACCCGCGAAGGGGCCGGCACAGCCCATAAAAACCTTCAGCTCCAGCTCAGAGTCGGCCACCCATTGATCTCGGCATGCTCGCGCAACACTGCGTCCGGGTTca                                              <  4:48052/149‑1 (MQ=255)
                                                                                                         tCGACCCTAGCCGGGTGTACGCGGGGCTGGCCCCTTCGCGGGTTTACCCGCGACGGGGCCGGCACAGCCCATAAAAACCTTCAGCTCCAGCTCAGAGTCGGCCACCCATTGATCTCGGCATGCTCGCGCAACACTGCGTCCGGGTTcac                                             <  3:1874/149‑1 (MQ=255)
                                                                                                              ccTAGCCTGGTGTAAGCGGGGCTGGCCCCTTCGCGGGTTTACCCGCGAAGGGGCCGGCACAGCCCATAAAAACCTTCAGCTCCAGCTCAGAGTCGGCCACCCATTGATCTCGGCATGCTCGCGCAACACTGCGTCCGGGTTCACCACAt                                        >  8:66284/1‑149 (MQ=255)
                                                                                                                     tGGTGTAAGCGGGGCTGGCCCCTTCGCGGGTTTACCCGCGAAGGGGCCGGCACAGCCCATAAAAACCTTCAGCTCCAGCTCAGAGTCGGCCACCCATTGATCTCGGCATGCTCGCGCAACACTGCGTCCGGGTTCACCACATGCGGGTg                                 >  5:251670/1‑149 (MQ=255)
                                                                                                                          tAAGCGGGGCTGGCCCCTTCGCGGGTTTACCCGCGAAGGGGCCGGCACAGCCCATAAAAACCTTCAGCTCCAGCTCAGAGTCGGCCACCCATTGATCTCGGCATGCTCGCGCAACACTGCGTCCGGGTTCACCACATGCGGGTGATCCa                            >  5:220358/1‑149 (MQ=255)
                                                                                                                              cGGGGCTGGCCCCTGCGCGGGTTACCCCGCGAAGGGGCCGGCACAGCCCATAAAAACCTTCAGCTCCAGCTCAGAGTCGGCCACCCATTGATCTCGGCATGCTCGCGCAACACTGCGTCCGGGTTCACCACATGCGGGTGATCCActct                        <  1:157101/149‑1 (MQ=255)
                                                                                                                               ggggCTGGCCCCTTCGCGGGTTTACCCGCGAAGGGGCCGGCACAGCCCATAAAAACCTTCAGCTCCAGCTCAGAGTCGGCCACCCATTGATCTCGGCATGCTCGCGCAACACTGCGTCCGGGTTCACCACATGCGGGTGATCCACtctc                       <  5:221366/149‑1 (MQ=255)
                                                                                                                               ggggCTGGCCCCTTCGCGGGTTTACCCGCGAAGGGGCCGGCACAGCCCATAAAAACCTTCAGCTCCAGCTCAGAGTCGGCCACCCATTGATCTCGGCATGCTCGCGCAACACTGCGTCCGGGTTCACCACATGCGGGTGATCCACtctc                       <  3:49448/149‑1 (MQ=255)
                                                                                                                               ggggCTGCCCCCTTCGCGGGTTTACCCGCTAAGGGGCCGGCACAGCCCATAAAAACCTTCAGCTCCAGCTCAGAGTCGGCCACCCATTGATCTCGGCATGCTCGCGCAACACTGCGTCTGGGTTCACCACATGCGGGTGATCCACtctc                       <  8:225589/149‑1 (MQ=255)
                                                                                                                                 ggCTGGCCCCTTCGCGGGTTTACCCGCGAAGGGGCCGGCACAGCCCATAAAAACCTTCAGCTCCAGCTCAGAGTCGGCCCCCCATTGATGTCGGCATGCGCGCGCCACCCTGGGTGCGGGTTCGCCACATGCGGGGGATCCGCtctcga                     >  6:366781/1‑147 (MQ=255)
                                                                                                                                  gCTGGCCCCTTCGCGGGTTTACCCGCGAAGGGGCCGGCACAGCCCATAAAAACCTTCAGCTCCAGCTCAGAGTCGGCCACCCATTGATCTCGGCATGCTCGCGCAACACTGCGTCCGGGTTCACCACATGCGGGTGATCCACTCTcaac                    <  6:165824/149‑1 (MQ=255)
                                                                                                                                        ccgttCGCGGGTTTACCCGCGATGGGGCCGGCACAGCCCATAAAAACCTTCAGCTCCAGCTCAGAGTCGGCCACCCATTGATCTCGGCATGCTCGCGCAACACTGCGTCCGGGTTCACCACATGCGGGTGATCCACTCTCaacaacaac              <  1:306932/146‑1 (MQ=255)
                                                                                                                                        cccTTCGCGGGTTTACCCGCGAAGGGGCCGGCACAGCCCATAAAAACCTTCAGCTCCAGCTCAGAGTCGGCCACCCATTGATCTCGGCATGCTCGCGCAACACTGCGTCCGGGTTCACCACATGCGGGTGATCCACTCTCaacaacaac              >  7:176710/1‑149 (MQ=255)
                                                                                                                                         ccTTCGCGGGTAAACCCGCGAAGGGGCCGGCACAGCCCATAAAAACCTTCAGCTCCAGCTCAGAGTCGGCCACCCATTGATCTCGGCATGCTCGCGCAACACTGCGTCCGGGTTCACCACATGCGGGTGATCCACTCTcaacaa                  >  3:20086/1‑144 (MQ=255)
                                                                                                                                         ccTTCGCGGGAAAACCCGCGAAGGGGCCGGCACAGCCCATAAAAACCTTCAGCTCCAGCTCAGAGTCGGCCACCCATTGATCTCGGCATGCTCGCGCAACACTGCGTCCGGGTTCACCACATGCGGGTGATCCACTCTcaacaa                  <  4:20086/144‑1 (MQ=255)
                                                                                                                                          cTTCGCGGGTAAACCCGCGAAGGGGCCGGCACAGCCCATAAAAACCTTCAGCTCCAGCTCAGAGTCGGCCACCCATTGATCTCGGCATGCTCGCGCAACACTGCGTCCGGGTTCACCACATGCGGGTGATCCACTCTCAACAACAACgg            >  6:189499/1‑149 (MQ=255)
                                                                                                                                          cTTCGCGGGTAAACCCGCGAAGGGGCCAGCACAGCCCATAAAAACCTTCAGCTCCAGCTCAGAGTCGGCCACCCATTGAtctc                                                                              >  4:99321/1‑83 (MQ=38)
                                                                                                                                          cTTCGCGGGTAAACCCGCGAAGGGGCCAGCACAGCCCATAAAAACCTTCAGCTCCAGCTCAGAGTCGGCCACCCATTGAtctc                                                                              <  3:99321/83‑1 (MQ=38)
                                                                                                                                           ttCGCGGGTTTACCCGCGAAGGGGCCGGCACAGCCCATAAAAACCTTCAGCTCCAGCTCAGAGTCGGCCACCCATTGATCTCGGCATGCTCGCGCAACACTGCGTCCGGGTTCACCACATGCGGGTGATCCACTCCCAACAACAACgg            >  7:69469/1‑148 (MQ=255)
                                                                                                                                           ttCGCGGGTAAACCCGCGAAGGGGCCGGCACAGCCCATAAAAACCTTCAGCTCCAGCTCAGAGTCGGCCACCCATTGATCTCGGCATGCTCGCGCAACACTGCGTCCGGGTTCACCACATGCGGGTGATCCACTCTCaacaacaac              <  5:295589/146‑1 (MQ=255)
                                                                                                                                           ttCGCGGGTAAACCCGCGAAGGGGCCGGCACAGCCCATAAAAACCTTCAGCTCCAGCTCAGAGTCGGCCACCCATTGATCTCGGCATGCTCGCGCAACACTGCGTCCGGGTTCACCACATGCGGGTGATCCACTCTCaacaacaac              >  6:295589/1‑146 (MQ=255)
                                                                                                                                              gcgGGTTTACCCGCGAAGGGGCCGGCACAGCCCATAAAAACCTTCAGCTCCAGCTCAGAGTCGGCCACCCATTGATCTCGGCATGCTCGCGCAACACTGCGTCCGGGTTCACCACATGCGGGTGATCCACTCTCAACAACAACGGCAAg        >  6:107670/1‑149 (MQ=255)
                                                                                                                                              gcgGGTAAACCCGCGAAGGGGCCGGCACAGCCCATAAAAACCTTCAGCTCCAGCTCAGAGTCGGCCACCCATTGATCTCGGCATGCTCGCGCAACACTGCGTCCGGGTTCACCACATGCGGGTGATCCACTCTCAACAACAACGGCAAg        >  6:344448/1‑149 (MQ=255)
                                                                                                                                               cgGGTTTACCCGCGAAGGGGCCGGCACAGCCCATAAAAACCTTCAGCTCCAGCTCAGAGTCGGCCACCCATTGATCTCGGCATGCTCGCGCAACACTGCGTCCGGGTTCACCACATGCGGGTGATCCACTCTCAACAACAACGGCAAGt       <  2:95943/149‑1 (MQ=255)
                                                                                                                                                     tACCCGCGAAGGGGCCGGCACAGCCCATAAAAACCTTCAGCTCCAGCTCAGAGTCGGCCACCCATTGATCTCGGCATGCTCGCGCAACACTGCGTCCGGGTTCACCACATGCGGGTGATCCACTCTCAACAACAACGGCAAGTCATTc  <  6:220358/148‑1 (MQ=255)
                                                                                                                                                     |                                                                                                                                                   
GCGGTTGTCGCCGCAGCTGGTCGACACCTATCCGGTGGAGTTT‑GCCCAGGCGGCGTATGCCGAGCTGGAGACCAACCAGGTGTCGGGCAAGCTGGTGATGGTGATCGACCCTAGCCTGGTGTAAGCGGGGCTGGCCCCTTCGCGGGTTTACCCGCGAAGGGGCCGGCACAGCCCATAAAAACCTTCAGCTCCAGCTCAGAGTCGGCCACCCATTGATCTCGGCATGCTCGCGCAACACTGCGTCCGGGTTCACCACATGCGGGTGATCCACTCTCAACAACAACGGCAAGTCATTC  >  NC_002947/1921402‑1921697

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