Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 3,534,805 T→G 10.4% intergenic (+49/‑106) atoB → / → PP_3124 3‑oxoacid CoA‑transferase subunit B/short chain fatty acid transporter

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0029473,534,8050TG10.4% 142.1 / 8.1 57intergenic (+49/‑106)atoB/PP_31243‑oxoacid CoA‑transferase subunit B/short chain fatty acid transporter
Reads supporting (aligned to +/- strand):  ref base T (21/30);  new base G (4/2);  total (25/32)
Fisher's exact test for biased strand distribution p-value = 3.88e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 7.93e-01

TTCTGCGCGAGACCGCGCCTGGGGTCAGCGTCGAGGAAATCATCGAGAAAACCGCCGGCAAACTGATCGTGCCGGATGATGTGAAGGAAATGACCTTCTAAGGCATGACGTGATCTCACCCTGTAGGCGCGGGTTTACCCGCGAATGCGTCGGTACATTGACCGCCGCATTCGCGGGTAAGCCTGTTCCTACAGGTTTTTTGCGTTGCCTGGCAATACCGCTTCATCCCATAAAACCAATAAAAGGAATTCAACCGTGGCCGCTGAAATTCAAGACAGCCGCTCCGCGCGCTTT  >  NC_002947/3534656‑3534949
                                                                                                                                                     |                                                                                                                                                 
ttCTGCGCGAGACCGCGCCTGGGGTCAGCGTCGAGGAAATCATCGAGAAAACCGCCGGCAAACTGATCGTGCCGGATGATGTGAAGGAAATGACCTTCTAAGGCATGACGTGATCTCACCCTGTAGGCGCGGGTTTACCCGCGAATGCGt                                                                                                                                                   <  2:1826760/150‑1 (MQ=255)
 tCTGCGCGAGACCGCGCCTGGGGTCAGCGTCGAGGAAATCATCGAGAAAACCGCCGGCAAACTGATCGTGCCGGATGATGTGAAGGAAATGACCTTCTAAGGCATGACGTGATCTCACCCTGTAGGCGCGGGTTTACCCGCGAATGCGTc                                                                                                                                                  <  1:682601/150‑1 (MQ=255)
  cTGCGCGAGACCGCGCCTGGGGTCAGCGTCGAGGAAATCATCGAGAAAACCGCCGGCAAACTGATCGTGCCGGATGATGTGAAGGAAATGACCTTCTAAGGCATGACGTGATCTCACCCTGTAGGCGCGGGTTTACCCGCGAATGCGTCg                                                                                                                                                 <  1:482146/150‑1 (MQ=255)
      gcgAGACCGCGCCTGGGGTCAGCGTCGAGGAAATCATCGAGAAAACCGCCGGCAAACTGATCGTGCCGGATGATGTGAAGGAAATGACCTTCTAAGGCATGACGTGATCTCACCCTGTAGGCGCGGGTTTACCCGCGAATGCGTCGGTAc                                                                                                                                             <  2:1060895/150‑1 (MQ=255)
      gcgAGACCGCGCCTGGGGTCAGCGTCGAGGAAATCATCGAGAAAACCGCCGGCAAACTGATCGTGCCGGATGATGTGAAGGAAATGACCTTCTAAGGCATGACGGGATCTCACCCTGAAGGCGCGGGTTTACCCGCGAGTGCGTCGGTAc                                                                                                                                             >  1:196636/1‑150 (MQ=255)
        gagaCCGCGCCTGGGGTCAGCGTCGAGGAAATCATCGAGAAAACCGCCGGCAAACTGATCGTGCCGGATGATGTGAAGGAAATGACCTTCTAAGGCATGACGTGATCTCACCCTGTAGGCGCGGGTTTACCCGCGAATGCGTCGGTACAt                                                                                                                                           <  2:1129562/150‑1 (MQ=255)
         agaCCGCGCCTGGGGTCAGCGTCGAGGAAATCATCGAGAAAACCGCCGGCAAACTGATCGTGCCGGATGATGTGAAGGAAATGACCTTCTAAGGCATGACGTGATCTCACCCTGTAGGCGCGGGTTTACCCGCGAATGCGTCGGTACAtt                                                                                                                                          <  2:402317/150‑1 (MQ=255)
           aCCGCGCCTGGGGTCAGCGTCGAGGAAATCATCGAGAAAACCGCGGGCAAACTGATCGTGCCGGATGATGTGAAGGAAATGACCTTCTAAGGCATGACGTGATCTCACCCTGTAGGCGCGGGTTTACCCGCGAATGCGTCGGTACATTGa                                                                                                                                       >  1:631255/1‑150 (MQ=255)
           aCCGCGCCTGGGGTCAGCGTCGAGGAAATCATCGAGAAAACCGCCGGCAAACTGATCGTGCCGGATGATGTGAAGGAAATGACCTTCTAAGGCATGACGTGATCTCACCCTGTAGGCGCGGGTTTACCCGCGAATGCGTCGGTACATTGa                                                                                                                                       <  2:1839467/150‑1 (MQ=255)
                  cTGGGGTCAGCGTCGAGGAAATCATCGAGAAAACCGCCGGCAAACTGATCGTGCCGGATGATGTGAAGGAAATGACCTTCTAAGGCATGACGTGATCTCACCCTGTAGGCGCGGGTTTACCCGCGAATGCGTCGGTACAtt                                                                                                                                          <  2:495085/141‑1 (MQ=255)
                  cTGGGGTCAGCGTCGAGGAAATCATCGAGAAAACCGCCGGCAAACTGATCGTGCCGGATGATGTGAAGGAAATGACCTTCTAAGGCATGACGTGATCTCACCCTGTAGGCGCGGGTTTACCCGCGAATGCGTCGGTACAtt                                                                                                                                          >  1:495085/1‑141 (MQ=255)
                   tGGGGTCAGCGTCGAGGAAATCATCGAGAAAACCGCCGGCAAACTGATCGTGCCGGATGATGTGAAGGAAATGACCTTCTAAGGCATGACGTGATCTCACCCTGTAGGCGCGGGTTTACCCGCGAATGCGTCGGTACATTGACCGCCGCa                                                                                                                               <  1:1900151/150‑1 (MQ=255)
                             gTCGAGGAAATCATCGAGAAAACCGCCGGCAAACTGATCGTGCCGGATGATGTGAAGGAAATGACCTTCTAAGGCATGACGTGATCTCACCCTGTAGGCGCGGGTTTACCCGCGAATGCGTCGGTACATTGACCGCCGCATTCGCgg                                                                                                                        <  2:1604114/147‑1 (MQ=255)
                             gTCGAGGAAATCATCGAGAAAACCGCCGGCAAACTGATCGTGCCGGATGATGTGAAGGAAATGACCTTCTAAGGCATGACGTGATCTCACCCTGTAGGCGCGGGTTTACCCGCGAATGCGTCGGTACATTGACCGCCGCATTCGCgg                                                                                                                        >  1:1604114/1‑147 (MQ=255)
                                          tcGAGAAAACCGCCGGCAAACTGATCGTGCCGGATGATGTGAAGGAAATGACCTTCTAAGGCATGACGTGATCTCACCCTGTAGGCGCGGGTTTACCCGCGAATGCGTCGGTACATTGACCGCCGCATTCGCgg                                                                                                                        <  1:1909019/134‑1 (MQ=255)
                                          tcGAGAAAACCGCCGGCAAACTGATCGTGCCGGATGATGTGAAGGAAATGACCTTCTAAGGCATGACGTGATCTCACCCTGTAGGCGCGGGTTTACCCGCGAATGCGTCGGTACATTGACCGCCGCATTCGCgg                                                                                                                        >  2:1909019/1‑134 (MQ=255)
                                          tcGAGAAAACCGCCGGCAAACTGATCGTGCCGGATGATGTGAAGGAAATGACCTTCTAAGGCATGACGTGATCTCACCCTGTAGGCGCGGGTTTACCCGCGAATGCGTCGGTACATTGACCGCCGCATTCGCGGGTAAGCCTGTTCCTAc                                                                                                        >  2:203416/1‑150 (MQ=255)
                                                           aaaCTGATCGTGCCGGATGATGTGAAGGAAATGACCTTCTAAGGCATGACGTGATCTCACCCTGTAGGCGCGGGTTTACCCGCGAATGCGt                                                                                                                                                   <  1:1664765/91‑1 (MQ=255)
                                                           aaaCTGATCGTGCCGGATGATGTGAAGGAAATGACCTTCTAAGGCATGACGTGATCTCACCCTGTAGGCGCGGGTTTACCCGCGAATGCGt                                                                                                                                                   >  2:1664765/1‑91 (MQ=255)
                                                            aaCTGATCGTGGCGGATGATGTGAAGGAAATGACCTGCTAAGGCATGACGTGATCTCACCCTGTAGGCGCGGGTTTACCCGCGAATGCGTCGGTACATTGACCGCCGCATTCGCGGGTAAGCCTGTTCCTACAGGTTTTTTGCGTTGCCt                                                                                      <  2:519669/150‑1 (MQ=255)
                                                            aaCTGATCGTGCCGGATGATGTGAAGGAAATGACCTTCTAAGGCATGACGTGATCTCACCCTGTAGGCGCGGGTTTACCCGCGAATGCGTCGGTTCATTGACCGCCGCATTCGCGGGTAAGCCTGTTCCTACAGGTTTTTTGCGTTGCCt                                                                                      >  1:1755982/1‑150 (MQ=255)
                                                            aaCTGATCGTGCCGGATGATGTGAAGGAAATGACCTTCTAAGGCATGACGTGATCTCACCCTGTAGGCGCGGGTTTACCCGCGAATGCGTCGGTACATTGACCGCCGCATTCGCGGGTAAGCCTGTTCCTACAGGTTTTTTGCGTTGCCt                                                                                      <  1:203416/150‑1 (MQ=255)
                                                             aCTGATCTTGCCGGATAATGTGAAAGAAATGACCTTCTAAGGCATGACGTGATCCCACCCTGTAGGCGCGGGTTTACCCGCGAATGCGTCGGTACCTTGACCTCCGCATTCGCGGGTAAGCCTGTTCCTACAGGTTTTTTGCGTTGCCTg                                                                                     <  1:163566/150‑1 (MQ=38)
                                                             aCTGATCGTGCCGGATGATGTGAAGGAAATGACCTTCTAAGGCATGACGTGATCTCACCCTGTAGGCGCGGGTTTACCCGCGAATGCGTCGGTACATTGACCGCCGCATTCGCGGGTAAGCCTGTTCCTACAGGTTTTTTGCGTTGCCTg                                                                                     <  1:1798053/150‑1 (MQ=255)
                                                             aCTGATCGTGCCGGATGATGTGAAGGAAATGACCTTCTAAGGCATGACGTGATCTCACCCTGTAGGCGCGGGTTTACCCGCGAATGCGTCGGTACATTGACCGCCGCATTCGCGGGTAAGCCTGTTCCTACAGGTTTTTTGCGTTGCCTg                                                                                     >  2:1272057/1‑150 (MQ=255)
                                                               tGATCGTGCCGGATGATGTGAAGGAAATGACCTTCTAAGGCATGACGTGATCTCACCCTGTAGGCGCGGGTTTACCCGCGAATGCGTCGGTACATTGACCGCCGCATTCGCGGGTAAGCCTGTTCCTACAGGTTTTTTGCGTTGCCTGGc                                                                                   <  1:1210533/150‑1 (MQ=255)
                                                                            tgatgTGAAGGAAATGACCTTCTAAGGCATGACGTGATCTCACCCTGTAGGCGCGGGTTTACCCGCGAATGCGTCGGTACATTGACCGCCGCATTCGCGGGTAAGCCTGTTCCTACAGGTTTTTTGCGTTGCCTGGCAATACCGCTTCAt                                                                      <  2:564402/150‑1 (MQ=255)
                                                                            tgatgTGAAGGAAATGACCTTCTAAGGCATGACGTGATCCCGACCTGTGGGCGCGGGTTTACCCGCGAATGAATAATTTCATTGACCGCCGCATTCGCGGGTAAGCCTGTTCCTACAGGTTTTTTGCGTTGCCTGGCAATACCGCTTCAt                                                                      <  1:126546/150‑1 (MQ=255)
                                                                               tgtgAAGGAAATGACCTTCTAAGGCATGACGTGATCTCACCCTGTAGGCGCGGGTTTACCCGCGAATGCGTCgg                                                                                                                                                >  1:1433941/1‑74 (MQ=255)
                                                                               tgtgAAGGAAATGACCTTCTAAGGCATGACGTGATCTCACCCTGTAGGCGCGGGTTTACCCGCGAATGCGTCgg                                                                                                                                                <  2:1433941/74‑1 (MQ=255)
                                                                                gtgAAGGAAATGACCTTCTAAGGCATGACGTGATCTCACCCTGTAGGCGCGGGTTTACCCGCGAATGCGTCGGTACATTGACCGCCGCATTCGCGGGTAAGCCTGTTCCTACAGGTTTTTTGCGTTGCCTGGCAATACCGCTTCATCCCa                                                                  >  1:1181810/1‑150 (MQ=255)
                                                                                 tgAAGGAAATGACCTTCTAAGGCATGACGTGATCTCACCCTGTAGGCGCGGGTTTACCCGAGAATGCGTCGGTACATTGACCGCCGCATTCGCgg                                                                                                                        >  1:526238/1‑95 (MQ=255)
                                                                                 tgAAGGAAATGACCTTCTAAGGCATGACGTGATCTCACCCTGTAGGCGCGGGTTTACCCGAGAATGCGTCGGTACATTGACCGCCGCATTCGCgg                                                                                                                        <  2:526238/95‑1 (MQ=255)
                                                                                 tgAAGGAAATGACCTTCTAAGGCATGACGTGATCTCACCCTGTAGGCGCGGGCTTACCCGCGAATGCGGCGGTCAATGTACCGACGCATTCGCgg                                                                                                                        <  1:898216/95‑1 (MQ=18)
                                                                                 tgAAGGAAATGACCTTCTAAGGCATGACGTGATCTCACCCTGTAGGCGCGGGCTTACCCGCGAATGCGGCGGTCAATGTACCGACGCATTCGCgg                                                                                                                        >  2:898216/1‑95 (MQ=18)
                                                                                 tgAAGGAAATGACCTCCTAAGGCATGACGTGACCTCACCCTGTAGGCGCGGGTTTACCCGCGAATGCGTCGCTATATTTACCGCCGCATTCGCGGGTAAGCCTGTTCCTACAGGTTTTTTGCGTTGCCTGGCAATACCGCTTCATCCCAt                                                                 <  1:126145/150‑1 (MQ=39)
                                                                                     ggAAATGACCTTCTAAGGCATGACGTGATCTCACCCTGTAGGCGCGGGTTTACCCGCGAATGCGTCGGTACATTGACCGCCGCATTCGCGGGTAAGCCTGTTCCTACAGGtttttt                                                                                               <  2:1354772/116‑1 (MQ=255)
                                                                                     ggAAATGACCTTCTAAGGCATGACGTGATCTCACCCTGTAGGCGCGGGTTTACCCGCGAATGCGTCGGTACATTGACCGCCGCATTCGCGGGTAAGCCTGTTCCTACAGGtttttt                                                                                               >  1:1354772/1‑116 (MQ=255)
                                                                                         aTGACCTTCTAAGGCATGACGTGATCTCACCCTGTAGGCGCGGGTTTCTCTGCGAATGCGTCGGTACATTGACCGCCGCATTCGCGGGTAAGCCTGTTCCTACAGGTTTTTTGCGTTGCCTGGCAATACCGCTTCATCCCATAAAACCaa                                                         <  1:416605/150‑1 (MQ=255)
                                                                                              cTTCTAAGGCATGACGTGATCTCACCCTGTAGGCGCGGGTTTACCCGCGAATGCGTCGGTACATTGACCGCCGCATTCGCGGGTAAGCCTGTTCCTACAGGTTTTTTGCGTTGCCTGGCAATACCGCTTCATCCCATAAAACCaaagatc                                                    >  1:5359/1‑145 (MQ=255)
                                                                                              cTTCTAAGGCATGACGTGATCTCACCCTGTAGGCGCGGGTTTACCCGCGAATGCGTCGGTACATTGACCGCCGCATTCGCGGGTAAGCCTGTTCCTACAGGTTTTTTGCGTTGCCTGGCAATACCGCTTCATCCCATAAAACCAATaaaa                                                    <  1:89070/150‑1 (MQ=255)
                                                                                                   aaGGCATGACGTGATCTCACCCTGTAGGCGCGGGTTTACCCGCGAATGCGTCGGTACATTGACCGCCGCATTCGCGGGTAAGCCTGTTCCTACAGGTTTTTTGCGTTGCCTGGCAATACCGCTTCATCCCATAAAACCAATAAAAGGAAt                                               <  1:920217/150‑1 (MQ=255)
                                                                                                      gCATGACGTGATCTCACCCTGTAGGCGCGGGTTTACCCGCGAATGCGTCGGTACATTGACCGCCGCATTCGCGGGTAAGCCTGTTCCTACAGGTTTTTTGCGTTGCCTGGCAATACCGCTTCATCCCATAAAACCAATAAAAGGAATTCa                                            >  1:1711655/1‑150 (MQ=255)
                                                                                                           aCGTGATCTCACCCTGTAGGCGCGGGTTTACCCGCGAATGCGTCGGTACATTGACCGCCGCATTCGCGGGTAAGCCTGTTCCTACAGGTTTTTTGCGTTGCCTGGCAATACCGCTTCATCCCATAAAACCaa                                                         <  1:1431450/132‑1 (MQ=255)
                                                                                                           aCGTGATCTCACCCTGTAGGCGCGGGTTTACCCGCGAATGCGTCGGTACATTGACCGCCGCATTCGCGGGTAAGCCTGTTCCTACAGGTTTTTTGCGTTGCCTGGCAATACCGCTTCATCCCATAAAACCaa                                                         >  2:1431450/1‑132 (MQ=255)
                                                                                                                   tcACCCTGTAGGCGCGGGTTTACCCGCGAATGCGTCGGTACATTGACCGCCGCATTCGCGGGTAAGCCTGTTCCTACAGGTTTTTTGCGTTGCCTGGCAATACCGCTTCATCCCATAAAACCAATAAAAGGAATTCAACCGTGGCCGCTg                               >  1:910601/1‑150 (MQ=255)
                                                                                                                     aCCCTGTAGGCGCGGGTTTACCCGCGAATGCGTCGGTACATTGACCGCCGCATTCGCGGGTAAGCCTGTTCCTACAGGTTTTTTGCGTTGCCTGGCAATACCGCTTCATCCCATAAAACCAATAAAAGGAATTCAACCGTGGCCGCTGaa                             >  2:1430211/1‑150 (MQ=255)
                                                                                                                       ccTGTAGGCGCGGGTTTACCCGCGAATGCGTCGGTACATTGACCGCCGCATTCGCGGGTAAGCCTGTTCCTACAGGTTTTTTGCGTTGCCTGGCAATACCGCTTCATCCCATAAAACCAATAAAAGGAATTCAACCGTGGCCGCTGAAAt                           >  2:1814931/1‑150 (MQ=255)
                                                                                                                        cTGTAGGCGCGGGTTTACCCGCGAATGCGTCGGTACa                                                                                                                                            >  2:447072/1‑37 (MQ=255)
                                                                                                                        cTGTAGGCGCGGGTTTACCCGCGAATGCGTCGGTACa                                                                                                                                            <  1:447072/37‑1 (MQ=255)
                                                                                                                         tGTAGGCGCGGGTTTACCCGCGAATGCGTCGGTACATTGACCGCCGCATTCGCGGGTAAGCCTGTTCCTACAGGTTTTTTGCGTTGCCTGGCAATACCGCtt                                                                         <  2:939493/102‑1 (MQ=255)
                                                                                                                         tGTAGGCGCGGGTTTACCCGCGAATGCGTCGGTACATTGACCGCCGCATTCGCGGGTAAGCCTGTTCCTACAGGTTTTTTGCGTTGCCTGGCAATACCGCtt                                                                         >  1:939493/1‑102 (MQ=255)
                                                                                                                                          ccGCGAATGCGGCGGTCAATGTACCGACGCATTCGCGGGTAAGCCTGTTCCTACAGGTTTTTTGCGTTGCCTGGCAATACCGCTTCAt                                                                      <  1:1245205/88‑1 (MQ=18)
                                                                                                                                          ccGCGAATGCGGCGGTCAATGTACCGACGCATTCGCGGGTAAGCCTGTTCCTACAGGTTTTTTGCGTTGCCTGGCAATACCGCTTCAt                                                                      >  2:1245205/1‑88 (MQ=18)
                                                                                                                                          ccGCGAATGCGGCGGTCAATGTACCGACGCATTCGCGGGTAAGCCTGTTCCTACAGGTTTTTTGCGTTGCCTGGCAATACCGCTTCATCCCATAAAACCAATAAAAGGAATTCAACCGTGGCCGCTGAAATTCAAGACAGCCGCTCcgcg        >  2:1266920/1‑150 (MQ=38)
                                                                                                                                          ccGCGAATGCGGCGGTCAATGTACCGACGCATTCGCGGGTAAGCCTGTTCCTACAGGTTTTTTGCGTTGCCTGGCAATACCGCTTCATCCCATAAAACCAATAAAAGGAATTCAACCGTGGCCGCTGAAATTCAAGACAGCCGCTCcgcg        >  2:1257761/1‑150 (MQ=38)
                                                                                                                                           cgcgAATGCGTCGGTACATTGACCGCCGCATTCGCGGGTAAGCCTGTTCCTACAGGTTTTTTGCGTTGCCTGGCAATACCGCTTCATCCCATAAAACCAATAAAAGGAATTCAACCGTGGCCGCTGAAATTCAAGACAGCCGCTCcgcgc       >  1:622971/1‑150 (MQ=255)
                                                                                                                                             cgAATGCGTCGGTACATTGACCGCCGCATTCGCGGGTAAGCCTGTTCCTACAGGTTTTTTGCGTTGCCTGGCAATACCGCTTCATCCCATAAAACCAATAAAAGGAATTCAACCGTGGCCGCTGAAATTCAAGACAGCCGCTCCgcgcgc     <  2:213577/150‑1 (MQ=255)
                                                                                                                                                aTGCGTCGGTACATTGACCGCCGCATTCGCGGGTAAGCCTGTTCCTACAGGTTTTTTGCGTTGCCTGGCAATACCGCTTCATCCCATAAAACCAATAAAAGGAATTCAACCGTGGCCGCTGAAATTCAAGACAGCCGCTCCGCGCGCttt  <  1:1814931/150‑1 (MQ=255)
                                                                                                                                                aTGCGTCAGTAAATTGACCTCAGTTTTCGCGGGTAAGCCTGTTCCTACAGGTTTTTTGCGTTGCCTTGATACGCCGCTTCATCCCATAAAACCAATAAAAGGAATTCAACCGTGGCCGCTGAAATTCAAGACAGCCGCTCCGCGCGCttt  <  1:594358/150‑1 (MQ=255)
                                                                                                                                                     |                                                                                                                                                 
TTCTGCGCGAGACCGCGCCTGGGGTCAGCGTCGAGGAAATCATCGAGAAAACCGCCGGCAAACTGATCGTGCCGGATGATGTGAAGGAAATGACCTTCTAAGGCATGACGTGATCTCACCCTGTAGGCGCGGGTTTACCCGCGAATGCGTCGGTACATTGACCGCCGCATTCGCGGGTAAGCCTGTTCCTACAGGTTTTTTGCGTTGCCTGGCAATACCGCTTCATCCCATAAAACCAATAAAAGGAATTCAACCGTGGCCGCTGAAATTCAAGACAGCCGCTCCGCGCGCTTT  >  NC_002947/3534656‑3534949

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: