Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 3,146,657 G→A 6.9% intergenic (+50/+16) PP_2761 → / ← sfnA ribose ABC transporter permease/acyl‑CoA dehydrogenase/oxidase family protein

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0029473,146,6570GA 6.9% 168.6 / 2.8 58intergenic (+50/+16)PP_2761/sfnAribose ABC transporter permease/acyl‑CoA dehydrogenase/oxidase family protein
Reads supporting (aligned to +/- strand):  ref base G (27/27);  new base A (2/2);  total (29/29)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.38e-01

GGCTGGAATTTCTTCTGGCAGCAGGTGGCCACCGGCGTTCTGATTCTCTCGGTTCTGCTGTTCAGTTTTGCCAGTCGGCGACTCGGCGCTGGATAGTTCTGGCGTTCCGGGCGGAACGCCAGTGTCGTGACGAGGCCTGGGGTATGCCCCAGGCCATGCTTTCAGGACTGCTTGGCACCGCCACCGATCTGCCATACATAGGGCGGCTCGCTGCCATTGACGTGCCAGTCGCCGATGATGCGTTCCTTGTAGATCACCGGGTTGTGCGAGGCGGCAGTACGGGCGTTGCGCC  >  NC_002947/3146512‑3146803
                                                                                                                                                 |                                                                                                                                                  
ggCTGGAATTTCTTCTGGCAGCAGGTGGCCACCGGCGTTCTGATTCTCTCGGTTCTGCTGTTCAGTTTTGCCAGTCGGCGACTCGGCGCTGGATAGTTCTGGCGTTCCGGTCGGAACGCCAGTGTCGTGACGAGGTCTGGGGTATGctcc                                                                                                                                                >  2:1275396/1‑150 (MQ=255)
     gAATTTCTTCTGGCAGCAGGTGGCCACCGGCGTTCTGATTCTCTCGGTTCTGCTGTTCAGTTTTGCCAGTCGGCGACTCGGCGCTGGATAGTTCTGGCGTTCCGGGCGGAACGCCAGTGTCGTGACGAGGCCTGGGGTATGCCCCAGGcc                                                                                                                                           >  1:517697/1‑150 (MQ=255)
     gAATTTCTTCTGGCAGCAGGTGGCCACCGGCGTTCTGATTCTCTCGGTTCTGCTGTTCAGTTTTGCCAGTCGGCGACTCGGCGCTGGATAGTTCTGGCGTTCCGGGCGGAACGCCAGTGTCGTGACGAGGCCTGGGGTATGCCCCAGGcc                                                                                                                                           >  1:517609/1‑150 (MQ=255)
         ttcttcTGGCAGCAGGTGGCCACCGGCGTTCTGATTCTCTCGGTTCTGCTGTTCAGTTTTGCCAGTCGGCGACTCGGCGCTGGATAGTTCTGGCGTTCCGGGCGGAACGCCAGTGTCGTGACGAGGCCTGGGGTATGCCCCAGGCCATGc                                                                                                                                       <  2:648425/150‑1 (MQ=255)
           cttctGGCAGCAGGTGGCCACCGGCGTTCTGATTCTCTCGGTTCTGCTGTTCAGTTTTGCCAGTCGGCGACTCGGCGCTGGATAGTTCTGGCGTTCCGGGCGGAACGCCAGTGTCGTGACGAGGCCTGGGGTATGCCCCAGGCCATGCtt                                                                                                                                     <  2:517609/150‑1 (MQ=255)
           cttctGGCAGCAGGTGGCCACCGGCGTTCTGATTCTCTCGGTTCTGCTGTTCAGTTTTGCCAGTCGGCGACTCGGCGCTGGATAGTTCTGGCGTTCCGGGCGGAACGCCAGTGTCGTGACGAGGCCTGGGGTATGCCCCAGGCCATGCtt                                                                                                                                     <  2:517697/150‑1 (MQ=255)
            ttctGGCAGCAGGTGGCCACCGGCGTTCTGATTCTCTCGGTTCTGCTGTTCAGTTTTGCCAGTCGGCGACTCGGCGCTGGATAGTTCTGGCGTTCCGGGCGGAACGCCAGTGTCGTGACGAGGCCTGGGGTATGCCCCAGGCCATGCttt                                                                                                                                    >  2:1250216/1‑150 (MQ=255)
               tGGCAGCAGGTGGCCACCGGCGTTCTGACTCTCTCGGGTCTGCTGTTCAGTTTTGCCAGTCGGCGACTCGGCGCTGGAGAGTTCCGGCGTTCCGGGCGGAACGCCAGTGTCGTGACGAGGCCTGGGGTATGCCCCAGGCCATGCTTTCAg                                                                                                                                 <  2:184204/150‑1 (MQ=255)
                   agcagGTGGCTACCGGCGTTCTGATTCTCTCGGTTCTGCTGTTCAGTTTTGCCAGTCGGCGACTCGGCGCTGGATAGTTCTGGCGTTCCGGGCGGAACGCCAGTGTCGTGACGAGGCCTGGGGTATGCCCCAGGCCATGCTTTCAGGACt                                                                                                                             >  2:660661/1‑150 (MQ=255)
                   agcagGTGGCCACCGGCGTTCTGATTCTCTCGGTTCTGCTGTTCAGTTTTGCCAGTCGGCGACTCGGCGCTGGATAGTTCTGGCGTTCCGGGCGGAACGCCAGTGTCGTGACGAGGCCTGGGGTATGCCCCAGGCCATGCTTTCAGGACt                                                                                                                             <  2:889674/150‑1 (MQ=255)
                      agGTGGCCACCGGCGTTCTGATTCTCTCGGTTCTGCTGTTCAGTTTTGCCAGTCGGCGACTCGGCGCTGGATAGTTCTGGCGTTCCGGGCGGAACGCCAGTGTCGTGACGAGGCCTGGGGTATGCCCCAGGCCATGCTTTCAGGACTGCt                                                                                                                          <  1:1069383/150‑1 (MQ=255)
                        gTGGCCACCGGCGTTCTGATTCTCTCGGTTCTGCTGTTCAGTTTTGCCAGTCGGCGACTCGGCGCTGGATAGTTCTGGCGTTCCGGGCGGAACGCCAGTGTCGTGACGAGGCCTGGGGTATGCCCCAGGCCATGCTTTCAGGACTGCTTg                                                                                                                        >  1:992630/1‑150 (MQ=255)
                           gCCACCGGCGTTCTGATTCTCTCGGTTCTGCTGTTCAGTTTTGCCAGTCGGCGACTCGGCGCTGGATAGTTCTGGCGTTCCGGGCGGAACGCCAGTGTCGTGACGAGGCCTGGGGTATGCCCCAGGCCATGCTTTCAGGACTg                                                                                                                            <  2:68479/143‑1 (MQ=255)
                           gCCACCGGCGTTCTGATTCTCTCGGTTCTGCTGTTCAGTTTTGCCAGTCGGCGACTCGGCGCTGGATAGTTCTGGCGTTCCGGGCGGAACGCCAGTGTCGTGACGAGGCCTGGGGTATGCCCCAGGCCATGCTTTCAGGACTg                                                                                                                            >  1:68479/1‑143 (MQ=255)
                             cACCGGCGTTCTGATTCTCTCGGTTCTGCTGTTCAGTTTTGCCAGTCGGCGACTCGGCGCTGGATAGTTCTGGCGTTCCGGGCGGAACGCCAGTGTCGTGACGAGGCCTGGGGTATGCCCCAGGCCATGCTTTCAGGACTGCTTGGCAcc                                                                                                                   >  2:996416/1‑150 (MQ=255)
                                  gCGTTCTGATTCTCTCGGTTCTGCTGTTCAGTTTTGCCAGTCGGCGACTCGGCGCTGGATAGTTCTGGCGTTCCGGGCGGAACGCCAGTGTCGTGACGAGGCCTGGGGTATGCCCCAGGCCATGCTTTCAGGACTGCTTGGCACCGCCAc                                                                                                              >  2:1062157/1‑150 (MQ=255)
                                         gATTCTCTCGGTTCTGCTGTTCAGTTTTGCCAGTCGGCGACTCGGCGCTGGATAGTTCTGGCGTTCCGGGCGGAACGCCAGTGTCGTGACGAGGCCTGGGGTATGCCCCAGGCCATGCTTTCAGGACTGCTTGGCACCGCCACCGATCTg                                                                                                       >  2:60412/1‑150 (MQ=255)
                                              tctcGGTTCTGCTGTTCAGTTTTGCCAGTCGGCGACTCGGCGCTGGATAGTTCTGGCGTTCCGGGCGGAACGCCAGTGTCGTGACGAGGCCTGGGGTATGCCCCAGGCCATGCTTTCAGGACTGCTTGGCACCGCCACCGATCTGCcata                                                                                                  >  2:770687/1‑150 (MQ=255)
                                                  ggTTCTGCTGTTCAGTTTTGCCAGTCGGCGACTCGGCGCTGGATAGTTCTGGCGTTCCGGTCGGAACGCCAGTGTCGTGACGAGGCCTGGGGTATGCCCCAGGCCATGCTTTCAGGACTGCTTGGCACCGCCACCGATCTGCcatacata                                                                                              >  2:113855/1‑150 (MQ=255)
                                                      ctgctgTTCAGTTTTGCCAGTCGGCGACTCGGCGCTGGATAGTTCTGGCGTTCCGGGCGGAACGCCAGTGTCGTGACGAGGCCTGGGGTATGCCCCAGGCCATGCTTTCAGGACTGCTTGGCACCGCCACCGATCTGCCATACATAGggc                                                                                          <  2:992630/150‑1 (MQ=255)
                                                       tgctgTTCAGTTTTGCCAGTCGGCGACTCGGCGCTGGATAGTTCTGGCGTTCCGGGCGGAACGCCAGTGTCGTGACGAGGCCTGGGGTATGCCCCAGGCCATGCTTTCAGGACTGCTTGGCACCGCCACCGATCt                                                                                                        >  1:220056/1‑135 (MQ=255)
                                                       tgctgTTCAGTTTTGCCAGTCGGCGACTCGGCGCTGGATAGTTCTGGCGTTCCGGGCGGAACGCCAGTGTCGTGACGAGGCCTGGGGTATGCCCCAGGCCATGCTTTCAGGACTGCTTGGCACCGCCACCGATCt                                                                                                        <  2:220056/135‑1 (MQ=255)
                                                       tgctgTTCAGTTTTGCCAGTCGGCGACTCGGCGCTGGATAGTTCTGGCGTTCCGGGCGGAACGCCAGTGTCGTGACGAGGCCTGGGGTATGCCCCAGGCCATGCTTTCAGGACTGCTTGGCACCGCCACCGATCTGCcatacata                                                                                              >  2:753192/1‑145 (MQ=255)
                                                       tgctgTTCAGTTTTGCCAGTCGGCGACTCGGCGCTGGATAGTTCTGGCGTTCCGGGCGGAACGCCAGTGTCGTGACGAGGCCTGGGGTATGCCCCAGGCCATGCTTTCAGGACTGCTTGGCACCGCCACCGATCTGCcatacata                                                                                              <  1:753192/145‑1 (MQ=255)
                                                        gctgTTCAGTTTTGCCAGTCGGCGACTCGGCGCTGGATAGTTCTGGCGTTCCGGGCGGAACGCCAGTGTCGTGACGAGGCCTGGGGTATGCCCCAGGCCATGCTTTCAGGa                                                                                                                               <  2:1035802/111‑1 (MQ=255)
                                                        gctgTTCAGTTTTGCCAGTCGGCGACTCGGCGCTGGATAGTTCTGGCGTTCCGGGCGGAACGCCAGTGTCGTGACGAGGCCTGGGGTATGCCCCAGGCCATGCTTTCAGGa                                                                                                                               >  1:1035802/1‑111 (MQ=255)
                                                           gTTCAGTTTTGCCAGTCGGCGACTCGGCGCTGGATAGTTCTGGCGTTCCGGGCGGAACGCCAGTGTCGTGACGAGGCCTGGGGTATGCCCCAGGCCATGCTTTCAGGACTGCTTGGCACCGCCACCGATCTGCcat                                                                                                   <  1:932612/136‑1 (MQ=255)
                                                           gTTCAGTTTTGCCAGTCGGCGACTCGGCGCTGGATAGTTCTGGCGTTCCGGGCGGAACGCCAGTGTCGTGACGAGGCCTGGGGTATGCCCCAGGCCATGCTTTCAGGACTGCTTGGCACCGCCACCGATCTGCcat                                                                                                   >  2:932612/1‑136 (MQ=255)
                                                             tCAGTTTTGCCAGTCGGCGACTCGGCGCTGGATAGTTCTGGCGTTCCGGGCGGAACGCCAGTGTCGTGACGAGGCCTGGGGCATGCCCCAGGCCATGCTTTCAGGACTGCTTGGCACCGCCACCGATCTGCCATACATAGGGCGGCTCGc                                                                                   >  1:1106816/1‑150 (MQ=255)
                                                                gTTTTGCCAGTCGGCGACTCGGCGCTGGATAGTTCTGGCGTTCCGGGCGGAACGCCAGTGTCGTGACGAGGCCTGGGGTATGCCCCAGGCCATGCTTTCAGGACTGCTTGGCACCGCCACCGATCTGCCATACATAGGGCGGCTCGCTGc                                                                                >  2:1242010/1‑150 (MQ=255)
                                                                 ttttGCCAGTCGGCGACTCGGCGCTGGATAGTTCTGGCGTTCCGGGCGGAACGCCAGTGTCGTGACGAGGCCTGGGGTATGCCCCAGGCCATGCTTTCAGGACTGCTTGGCACCGCCACCGATCTGCCATACATAGGGCGGCTCGCTGcc                                                                               <  1:1242010/150‑1 (MQ=255)
                                                                 ttttGCCAGTCGGCGACTCGGCGCTGGATAGTTCTGGCGTTCCGGGCGGAACGCCAGTGTCGTGACGAGGCCTGGGGCATGCCCCAGGCCATGCTTTCAGGACTGCTTGGCACCGCCACCGATCTGCCATACATAGGGCGGCTCGCTGcc                                                                               <  2:1106816/150‑1 (MQ=255)
                                                                             gCGACTCGGCGCTGGATAGTTCTGGCGTTCCGGGCGGAACGCCAGTGTCGTGACGAGGCCTGGGGTATGCCCCAGGCCATGCTTTCAGGACTGCTTGGCACCGCCACCGATCTGCCATACATAGGGCGGCTCGCTGCCATTGACGTGCCa                                                                   <  1:113855/150‑1 (MQ=255)
                                                                             gCGACTCGGCGCTGGATAGTTCTGGCGTTCCGGGCGGAACGCCAGTGTCGTGACGAGGCCTGGGGCATACCCCAGGcctc                                                                                                                                         >  2:1005566/1‑78 (MQ=39)
                                                                             gCGACTCGGCGCTGGATAGTTCTGGCGTTCCGGGCGGAACGCCAGTGTCGTGACGAGGCCTGGGGCATACCCCAGGcctc                                                                                                                                         <  1:1005566/80‑3 (MQ=39)
                                                                                   cGGCGCTGGATAGTTCTGGCGTTCCGGGCGGAACGCCAGTGTCGTGACGAGGCCTGGGGTATGCCCCAGGCCATGCTTTCAGGACTGCTTGGCACCGCCACCGATCTGCCATACATAGGGCGGCTCGCTGCCATTGACGTGCCAGTCGcc                                                             <  1:996416/150‑1 (MQ=255)
                                                                                            aTAGTTCTGGCGTTCCGGGCGGAACGCCAGTGTCGTGACGAGGCCTGGGGTATGCCCCAGGCCATGCTTTCAGGACTGCTTGGCACCGCCACCGATCTGCCATACATAGGGCGGCTCGCTGCCATTGACGTGCCAGTCGCCGATGATGCg                                                    >  2:147264/1‑150 (MQ=255)
                                                                                              gacacTGGCGTTCCGCCCGGAACGCCAGTGTCGTGACGAGGCCTGGGGTATGCCCCAGGCCATGCTTTCAGGACTGCTTGGCACCGCCACCGATCTGCCATACATAGGGCGGCTCGCTGCCATTGACGTGCCAGTCGCCGATGATGCGtt                                                  >  1:1192458/5‑150 (MQ=255)
                                                                                               acacTGGCGTTCCGCCCGGAACGCCAGTGTCGTGACGAGGCCTGGGGTATGCCCCAGGCCATGCTTTCAGGACTGCTTGGCACCGCCACCGATCTGCCATACATAGGGCGGCTCGCTGCCATTGACGTGCCAGTCGCCGATGATGCGTTc                                                 >  1:770723/4‑150 (MQ=255)
                                                                                                   tGGCGTTCCGTGCGGAAAGCCAGTGTCGTCACGAGGCCCGGGGGATGCCCCAGGCCATGCTTTCAGGACTGCTTGGCACCGCCACCGATCTGCCATACATAGGGCGGCTCGCTGCCATTGACGTGCCAGTCGCCGATGATGCGTTCCTTg                                             <  1:401280/150‑1 (MQ=255)
                                                                                                   tGGCGTTCCGGGCGGAACGCCAGTGTCGTGACGAGGCCTGGGGTATGCCCCAGGCCATGCTTTCAGGACTGCTTGGCACCGCCACCGATCTGCCATACATAGGGCGGCTCGCTGCCATTGACGTGCCAGTCGCCGATGATGCGTTCCTTg                                             <  1:147264/150‑1 (MQ=255)
                                                                                                    ggCGTTCCGGGCGGAACGCCAGTGTCGTGACGAGGCCTGGGGTATGCCCCAGGCCATGCTTTCAGGACTGCTTGGCACCGCCACCGATCTGCCATACATAGGGCGGCTCGCTGCCATTGACGTGCCAGTCGCCGATGATGCGTTCCTTGt                                            <  1:60412/150‑1 (MQ=255)
                                                                                                    ggCGTTCCGCCCGGAACGCCAGTGTCGTGACGAGGCCTGGGGTATGCCCCAGGCCATGCTTTCGGGACTGCTTGGCACCGc                                                                                                                 >  1:760181/1‑81 (MQ=38)
                                                                                                    ggCGTTCCGCCCGGAACGCCAGTGTCGTGACGAGGCCTGGGGTATGCCCCAGGCCATGCTTTCGGGACTGCTTGGCACCGc                                                                                                                 <  2:760181/81‑1 (MQ=38)
                                                                                                     gCGTTCCGGGCGGAACGCCAGTGTCGTGCCGAGGCCTGGGGTATGCCCCAGGCCATGCTTTCAGGACTGCTTGGCACCGCCACCGATCTGCCATACATAGGGCGGCTCGCTGCCATTGACGTGCCAGTCGCCGATGATGCGTTCCTTGTa                                           <  1:166249/150‑1 (MQ=255)
                                                                                                     gCGTTCCGCCCGGAACGCCAGTGTCGTGACGAGGCCTGGGGTATGCCCCAGGCCATGCTTTCAGGACTGCTTGGCACCGCCACCGATCTGCCATACATAGGGCGGCTCGCTGCCATTGACGTGCCAGTCGCCGATGATGCGTTCCTTGTa                                           <  2:770723/150‑1 (MQ=255)
                                                                                                       gTTCCGGGCGGAACGCCAGTGTCGTGACGAGGCCTGGGGTATGCCCCAGGCCa                                                                                                                                          >  2:469615/1‑53 (MQ=255)
                                                                                                       gTTCCGGGCGGAACGCCAGTGTCGTGACGAGGCCTGGGGTATGCCCCAGGCCa                                                                                                                                          <  1:469615/53‑1 (MQ=255)
                                                                                                              gCGGAACGCCAGTGTCGTGACGAGGCCTGGGGTATGCCCCAGGCCATGCTTTCAGGACTGCTTGGCACCGCCACCGATCTGCCATACATAGGGCGGCTCGCTGCCATTGACGTGCCAGTCGCCGATGATGCGTTCCTTGTAGATCACCgg                                  <  1:660661/150‑1 (MQ=255)
                                                                                                                     gCCAGTGTCGTGACGAGGCCTGGGGTATGCCCCAGGCCATGCTTTCTGGAGTAATCCCCAGCGTCACCTATCTGACATACATAGGGCGGCTCGCTGCCATTGACGTGCCAGTCGCCGATGATGCGTTCCTTGTAGATCACCGGGTTGTGc                           <  1:1250216/150‑1 (MQ=255)
                                                                                                                     gCCAGTGTCGTGACGAGGCCTGGGGTATGCCCCAGGCCATGCTTTCAGGACTGCTTGGCACCGCCACCGATCTGCCATACATAGGGCGGCTCGCTGCCATTGACGTGCCAGTCGCCGATGATGCGTTCCTTGTAGATCACCGGGTTGTGc                           <  1:770687/150‑1 (MQ=255)
                                                                                                                     gCCAGTGTCGTGACGAGGCCTGGGGCATACCCCAGGc                                                                                                                                            >  2:732806/1‑37 (MQ=18)
                                                                                                                     gCCAGTGTCGTGACGAGGCCTGGGGCATACCCCAGGc                                                                                                                                            <  1:732806/37‑1 (MQ=18)
                                                                                                                            tCGTGACGAGGCCTGGGGTATGCCCCAGGCCATGCTTTCAGGACTGCTTGGCACCGCCACCGATCTGCCATACATAGGGCGGCTCGCTGCCATTGACGTGCCAGTCGCCGATGATGCGTTCCTTGt                                            >  1:65750/1‑126 (MQ=255)
                                                                                                                            tCGTGACGAGGCCTGGGGTATGCCCCAGGCCATGCTTTCAGGACTGCTTGGCACCGCCACCGATCTGCCATACATAGGGCGGCTCGCTGCCATTGACGTGCCAGTCGCCGATGATGCGTTCCTTGt                                            <  2:65750/126‑1 (MQ=255)
                                                                                                                                   gAGGCCTGGGGTATGCCCCAGGCCATGCTTTCAGGACTGCTTGGCACCGCCACCGATCTGCCATACATAGGGCGGCTCGCTGCCATTGACGTGCCAGTCGCCGATGATGCGTTCCTTGTAGATCACCGGGTTGTGCGAGGCGGCAGTACg             >  2:610721/1‑150 (MQ=255)
                                                                                                                                        cTGGGGTATGCCCCAGGCCATGCTTTCAGGACTGCTTGGCACCGCCACCGATCTGCCATACATAGGGCGGCTCGCTGCCATTGACGTGCCAGTCGCCGATGATGCGTTCCTTGTAGATCACCGGGTTGTGCGAGGCGGCAGTACGGGCGt        >  1:507615/1‑150 (MQ=255)
                                                                                                                                              tATGCCCCAGGCCATGCTTTCAGGACTGCTTGGCACCGCCACCGATCTGCCATACATAGGGCGGCTCGCTGCCATTGACGTGCCAGTCGCCGATGATGCGTTCCTTGTAGATCACCGGGTTGTGCGAGGCGGCAGTACGGGCGTTGCGcc  >  2:1061995/1‑150 (MQ=255)
                                                                                                                                                 |                                                                                                                                                  
GGCTGGAATTTCTTCTGGCAGCAGGTGGCCACCGGCGTTCTGATTCTCTCGGTTCTGCTGTTCAGTTTTGCCAGTCGGCGACTCGGCGCTGGATAGTTCTGGCGTTCCGGGCGGAACGCCAGTGTCGTGACGAGGCCTGGGGTATGCCCCAGGCCATGCTTTCAGGACTGCTTGGCACCGCCACCGATCTGCCATACATAGGGCGGCTCGCTGCCATTGACGTGCCAGTCGCCGATGATGCGTTCCTTGTAGATCACCGGGTTGTGCGAGGCGGCAGTACGGGCGTTGCGCC  >  NC_002947/3146512‑3146803

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