Predicted mutation
evidence seq id position mutation freq annotation gene description
RA CP009273 14,485:1 +G 71.8% coding (318/1131 nt) dnaJ → chaperone Hsp40, DnaK co‑chaperone

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*CP00927314,4851.G71.8% 92.0 / 39.8 39coding (318/1131 nt)dnaJchaperone Hsp40, DnaK co‑chaperone
Reads supporting (aligned to +/- strand):  ref base . (6/5);  new base G (15/13);  total (21/18)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 5.49e-01

GCCATGAAATACCACCCGGACCGTAACCAGGGTGACAAAGAGGCCGAGGCGAAATTTAAAGAGATCAAGGAAGCTTATGAAGTTCTGACCGACTCGCAAAAACGTGCGGCATACGATCAGTATGGTCATGCTGCGTTTGAGCAAGGTGGCATGGGCGGCGGCGGTTTTGGCGGCGGCGCAGACTTCAGCGATATTTTTGGTGACGTTTTCGGCGATATTTTTGGCGGCGGACGTGGTCGTCAACGTGCGGCGCGCGGTGCTGATTTACGCTATAACATGGAGCTCACCCTCGAAGAAGCTGTACGTGGCGTGACCAAAGAGATCCGCATTCCGACTCTGGAAGAGTGTGACGTTTGCCACGGTAGCGGTGCAAAACCAGGTACACAGCCGCAGACTTGTCCGACCTGTCATGGTTCTGGTCAGGTGCAGATGCGCCAGGGATTCTTCGCTGTACAGCAGAC  >  CP009273/14252‑14712
                                                                                                                                                                                                                                          |                                                                                                                                                                                                                                   
gCCATGAAATACCACCCGGACCGTAACCAGGGTGACAAAGAGGCCGAGGCGAAATTTAAAGAGATCAAGGAAGCTTATGAAGTTCTGACCGACTCGCAAAAACGTGCGGCATACGATCAGTATGGTCATGCTGCGTTTGAGCAAGGTGGCATGGGCGGCGGCGGTTTTGGCGGCGGCGCAGACTTCAGCGATATTTTTGGTGCCGTTTTCGGCGATATTTTTGGCGGCGGACGTGGGTCGTCAACGTGCGGCGCGCGGTGCTGATTTACGCtat                                                                                                                                                                                              >  1:334835/1‑274 (MQ=255)
     gAAATACCACCCGGACCGTAACCAGGGTGACAAAGAGGCCGAGGCGAAATTTAAAGAGATCAAGGAAGCTTATGAAGTTCTGACCGACTCGCAAAAACGTGCGGCATACGATCAGTATGGTCATGCTGCGTTTGAGCAAGGTGGCATGGGCGGCGGCGGTTTTGGCGGCGGCGCAGACTTCAGCGATATTTTTGGTGACGTTTTCGGCGATATTTTTGGCGGCGGACGTGGGTCGTCAACGTGCGGCGCGGGGTGCTGATTTACGCTATAACATg                                                                                                                                                                                        >  1:7372/1‑275 (MQ=255)
                 ggACCGTAACCAGGGTGACAAAGAGGCCGAGGCGAAATTTAAAGAGATCAAGGAAGCTTATGAAGTTCTGACCGACTCGCAAAAACGTGCGGCATACGATCAGTATGGTCATGCTGCGTTTGAGCAAGGTGGCATGGGCGGCGGCGGTTTTGGCGGCGGCGCAGACTTCAGCGATATCTTTGGTGACGTTTTCGGCGATATTTTTGGCGGCGGACGTGGGTCGTCAACGTGCGGCTCGCGGTGCTGATTTACGCTATAACATGTAGCTCACCCTc                                                                                                                                                                            >  1:424411/1‑275 (MQ=255)
                   aCCGTAACCAGGGTGACAAAGAGGCCGAGGCGAAATTTAAAGAGATCAAGGAAGCTTATGAAGTTCTGACCGACTCGCAAAAACGTGCGGCATACGATCAGTATGGTCATGCTGCGTTTGAGCAAGGTGGCATGGGCGGCGGCGGTTTTGGCGGCGGCGCAGACTTCAGCGATATTTTTGGTGACGTTTTCGGCGATATTTTTGGCGGCGGACGTg                                                                                                                                                                                                                                    <  1:797992/216‑1 (MQ=255)
                                           ccGAGGCGAAATTTAAAGAGATCAAGGAAGCTTATGAAGTTCTGACCGACTCGCAAAAACGTGCGGCATACGATCAGTATGGTCATGCTGCGTTTGAGCAAGGTGGCATGGGCGGCGGCGGTTTTGGCGGCGGCGCAGACTTCAGCGATATTTTTGGTGACGTTTTCGGCGATATTTTTGGCGGCGGACGTGGGTCGTCAACGTgcggcg                                                                                                                                                                                                                   >  2:790303/1‑210 (MQ=255)
                                               ggCGAAATTTAAAGAGATCAAGGAAGCTTATGAAGTTCTGACCGACTCGCAAAAACGTGCGGCATACGATCAGTATGGTCATGCTGCGTTTGAGCAAGGTGGCATGGGCGGCGGCGGTTTTGGCGGCGGCGCAGACTTCAGCGATATTTTTGGTGACGTTTTCGGCGATATTTTTGGCGGCGGACGTGGGTCGTCAACGTGCGGcgcgcg                                                                                                                                                                                                               <  1:790309/210‑1 (MQ=255)
                                                          aaGAGAGCAAGGACGATTATGTAGTTCTGACCGACTCTCAAAAACGTGCGGCATACGATCAGTATGTTCATGCTGCGTTTGAGGAAGGTGGNATGGGCGGCGGCGGTTTTGGCGGCGGCGCAGACTTCAGCGATATTTTTGGTGCCGTTTTCGGCGATATTTTTGGCGGCGGACGTGGGTCGTCAACGTGCGGCGCGCGGTGCTGATTTACGCTATAACATGGAGCTCACCCTCGAAGAAGCTGTACGTGCAGTGACCATAGAGATCCGCATTcc                                                                                                                                   <  1:569325/275‑1 (MQ=255)
                                                                       aGCTTATGAAGTTCTGACCGACTCGCAAAAACGTGCGGCATACGATCAGTATGGTCATGCTGCGTTTGAGCAAGGTGGCATGGGCGGCGGCGGTTTTGGCGGCGGCGCAGACTTCAGCGATATTTTTGGTGACGTTTTCGGCGATATTTTTGGCGGCGGACGTGGGCCGTCAACGTGCGGCGCGCGGTGCTGATTTACGCTATGACATGGAGCTCACCCTCGAAGAAGCTGTACGTGGCGAGAACAAAGAGAACAGAATTCCGAAGCAGGAagag                                                                                                                      >  1:75535/1‑275 (MQ=255)
                                                                             tGAAGTTATGACCGAATCGCAAAAACGTGCGGCATACGATCAGTATGGTAATGCTGCGTTTGAGCAAGGTGGCATGGGCGGCGGCGGTTTTGGCGGCGGCGCAGACTTCAGCGATATTTTTGGTGACGTTTTCGGCGATATTTTTGGCGGCGGACGTGGGTCGTCAACGTGCGgcgcg                                                                                                                                                                                                                 <  1:416393/178‑1 (MQ=255)
                                                                                                             cATACGATCAGTATGGTCATGCTGCGTTTGAGCAAGGTGGCATGGGCGGCGGCGGTTTTGGCGGCGGCGCAGACTTCAGCGATATTTTTGGTGACGTTTTCGGCGATATTTTTGGCGGCGGACGTGGTCGTCAACGTgcgg                                                                                                                                                                                                                     >  2:257872/1‑141 (MQ=255)
                                                                                                                 cGATCAGTATGGTCATGCTGCGTTTGAGCAAGGTGGCATGGGCGGCGGCGGTTTTGGCGGCGGCGCAGACTTCAGCGATATTTTTGGTGACGTTTTCGGCGATATTTTTGGCGGCGGACGTGGTCGTCAACGTGCGgcgcg                                                                                                                                                                                                                 <  1:257875/141‑1 (MQ=255)
                                                                                                                               aggctgcGTNTTAGCAAGTTGGCATGGGGGGGGGCGGGTTTGGCGGCGGCGCTGGTTTCAGCGATATTTTTGGGGAGGTTTTCGGCGATATTTTTGGCTGCGGACGTGGGTCGTCAACGTGCGGCGCGCGGTGCTGATTTACGCTATGACATGGAGTTCACCCTCGAAGAAGCTGTACGTGGCGTGACCAAAGAGATCCGCATTCCGACTCTGGAAGAGTGTGACGTTTGCCACGGTAGCGGTGCAAAACCAGGTACACAGCCGCAGACTTGt                                                                <  2:75535/271‑1 (MQ=255)
                                                                                                                                 gctgcGTTTGAGCAAGGTGGCATGGGCGGCGGCGGTTTTGGCGGCGGCGCAGACTTCAGCGATATTTTTGGTGACGTTTTCGGCGATATTTTTGGCGGCGGACGTGGGTCGTCAACGTGCGGCGCGCGGTGCTGATTTACGCTATAACATGGAGCTCACCCTCGAAGAAGCTGTACGTGGc                                                                                                                                                          >  2:544517/1‑181 (MQ=255)
                                                                                                                                  ctgcGTTTGAGCAAGGTGGCATGGGCGGCGGCGGTTTTGGCGGCGGCGCAGACTTCAGCGATATTTTTGGTGACGTTTTCGGCGATATTTTTGGCGGCGGACGTGGGTCGTCAACGTGCGGCGCGCGGTGCTGATTTACGCTATAACATgg                                                                                                                                                                                       >  1:757901/1‑151 (MQ=255)
                                                                                                                                  ctgcGTTTGAGCAAGGTGGCATGGGCGGCGGCGGTTTTGGCGGCGGCGCAGACTTCAGCGATATTTTTGGTGACGTTTTCGGCGATATTTTTGGCGGCGGACGTGGTCGTCAACGTGCGgcgcgc                                                                                                                                                                                                                >  2:768616/1‑125 (MQ=255)
                                                                                                                                  ctgcGTTTGAGCAAGGTGGCATGGGCGACGGCGGTTTTGGCGGCGGCGCAGACTTCAGCGATATTTTTGGTGACGTTTTCGGCGATCTTTTTGGCGTCGGACGTGGTCGTCAACGTGCGGCGCGCGGTGCTGATTTACGCTATAACATGGAGCTCACACTCGACGCAGCTGTACGTGGCGTGACCAAAGAGATCCGCATTCCGACTCTGGCAGAGTGTGACGTTTGCCACGGTAGCGGTGCAAAACCAGGTCCACAGCCGCAGACTTGTCCGAcc                                                          >  1:502397/1‑275 (MQ=255)
                                                                                                                                     cGTTTGAGCACGGTGGCATGGGCGGCGGCGGTTTTGGCGGCGGCGCAGACTTCAGCGATATTTTTGGTGACGTTTTCGGCGATATTTTTGGCGGCGGACGTGGGTCGTCAACGTGCGGCGCGCGGTGCTGATTTACGCTATAACATTGAGCTCACCCTCGAAGAAGCTGTACGTGGCGTGa                                                                                                                                                      <  1:544520/181‑1 (MQ=255)
                                                                                                                                     cGTTTGAGCAAGGTGGCATGGGCGGCGGCGGTTTTGGCGGCGGCGCAGACTTCAGCGATATTTTTGGTGACGTTTTCGGCGATATTTTTGGCGGCGGACGTGGGTCGTCAACg                                                                                                                                                                                                                          >  2:732749/1‑113 (MQ=255)
                                                                                                                                      gTTTGAGCAAGGTGGCATGGGCGGCGGCGGTTTTGGCGGCGGCGCAGACTTCAGCGATATTTTTGGTGACGTTTTCGGCGATATTTTTGGCGGCGGACGTGGGTCGTCAACGTgcgg                                                                                                                                                                                                                     >  2:128195/1‑117 (MQ=255)
                                                                                                                                      gTTTGAGCAAGGTGGCATGGGCGGCGGCGGTTTTGGCGGCGGCGCAGACTTCAGCGATATTTTTGGTGACGTTTTCGGCGATATTTTTGGCGGCGGACGTGGGTCGTCAACGTGCGGCGCGCGGTGCTGATTTACGCTATAACATGGAGCt                                                                                                                                                                                   <  2:757895/151‑1 (MQ=255)
                                                                                                                                      gTTTGAGCAAGGTGGCATGGGCGGCGGCGGTTTTGGCGGCGGCGCAGACTTCAGCGATATTTTTGGTGACGTTTTCGGCGATATTTTTGGCGGCGGACGTggg                                                                                                                                                                                                                                  <  1:25103/103‑2 (MQ=255)
                                                                                                                                      gTTTGAGCAAGGTGGCATGGGCGGCGGCGGTTTTGGCGGCGGCGCAGACTTCAGCGATATTTTTGGTGACGTTTTCGACGATATTTTTGGCGGCGGACGTGGTCGTCAACGTGCGGCGCGCGGTg                                                                                                                                                                                                            <  1:768622/125‑1 (MQ=255)
                                                                                                                                       tttGAGCAAGGTGGCATGGGCGGCGGCGGTTTTGGCGGCGGCGCAGACTTCAGCGATATTTTTGGTGACGTTTTCGGCGATATTTTTGGCGGCGGACGTGGGTCGTCAACGTGCGGCGCGCGGTGc                                                                                                                                                                                                           >  2:17168/1‑126 (MQ=255)
                                                                                                                                       tttGAGCAAGGTGGCATGGGCGGCGGCGGTTTTGGCGGCGGCGCAGACTTCAGCGATATTTTTGGTGACGTTTTCGGCGATATTTTTGGCGGCGGACGTGGTCGTCAACGTGCGGCGCGCGGTGCTGATTTACGCTATAACATGGAGCTCACCCTCgaag                                                                                                                                                                        >  1:209596/1‑160 (MQ=255)
                                                                                                                                       tttGAGCAAGATGGCATGGGCGGCGGCGGTTTTGGCGGCGGCGCAGACTTCAGCGATATTTTTGGTGACGTTTTCGGCGATATTTTTGGCGGCGGACGTGGTCGTCAACGTGCGGCGCGCGGTGCTGa                                                                                                                                                                                                        >  1:812763/1‑128 (MQ=255)
                                                                                                                                         tGAGCAAGGTGGCATGGGCGGCGGCGGTTTTGGCGGCGGCGCAGACTTCAGCGATATTTTTGGTGACGTTTTCGGCGATATTTTTGGCGGCGGACGTGGGTCGTCAACGTgcg                                                                                                                                                                                                                      <  1:732755/113‑1 (MQ=255)
                                                                                                                                          gAGCAAGGTGGCATGGGCGGCGGCGGTTTTGGCGGCGGCGCAGACTTCAGCGATATTTTTGGTGACGTTTTCGGCGATATTTTTGGCGGCGGACGTGGGTCGTCAACGTGCGgcgcg                                                                                                                                                                                                                 <  1:128196/117‑1 (MQ=255)
                                                                                                                                           aGCAAGGTGGCATGGGCGGCGGCGGTTTTGGCGGTGGCGCAGACTTCTGCGATATTTTTGGTGACGTTTTCGGCGATATTTTTGGCGGCGGACGTGGTCGTCAACGTGCGGCGCGCGGTGCTGATTTa                                                                                                                                                                                                    <  2:812756/128‑1 (MQ=255)
                                                                                                                                           aGCAAGGTGGCATGGGCGGCGGCGGTTTTGGCGGCGGCGCAGACTTCAGCGATATTTTTGGTGACGTTTTCGGCGATATTTTTGGCGGCGGACGTGGGTCGTCAACGTGCGGCGCGCGGTGCTGAt                                                                                                                                                                                                       <  1:17168/126‑1 (MQ=255)
                                                                                                                                           aGCAAGGTGGCATGGGCGGCGGCGGTTTTGGCGGCGGCGCAGACTTCAGCGATATTTTTGGTGACGTTTTCGGCGATATTTTTGGCGGCGGACGTg                                                                                                                                                                                                                                    <  2:124150/96‑1 (MQ=255)
                                                                                                                                           aGCAAGGTGGCATGGGCGGCGGCGGTTTTGGCGGCGGCGCAGACTTCAGCGATATTTTTGGTGACGTTTTCGGCGATATTTTTGGCGGCGGACGTGGTCGTCAACGTGCGGCGCGCGGTGCTGATTTACGCTATAACATGGAGCTCACCCTCGAAGAAGc                                                                                                                                                                    <  2:209593/160‑1 (MQ=255)
                                                                                                                                                         ggcggcggcGGTTTTGGCGGCGGCGCAGACTTCAGCGATATTTTTGGTGACGTTTTCGGCGATATTTTTGGCGGCGGACGTGGGTCGTCAACGTGCGGCGCGCGGTGCTGAtt                                                                                                                                                                                                      >  2:186592/1‑113 (MQ=255)
                                                                                                                                                             gcggcggTTTTGGCGGCGGCGCAGACTTCAGCGATATTTTTGGTGACGTTTTCGGCGATATTTTTGGCGGCGGACGTGGGTCGTCAACGTGCGGCGCGCGGTGCTGATTTACg                                                                                                                                                                                                  <  1:186595/113‑1 (MQ=255)
                                                                                                                                                                gggtttttGGAGGCGGGGGACAATTCGGCGGTCTTTTTGGTGNCGTTTTCGGCGATAGTTGTGGCGGCGGACGTGGGTCGTCAACGTGCGGCGCGCGGTGCTGATTTACGCTATAACATGGAGCTCACCCTCGAAGAAGCTGTACGTGGCGTGACCAAAGAGATCCGCATTCCGACTCTGGAAGAGTGTGACGTTTGCCACGGTAGCGGTGCAAAACCAGGTACACAGCCGCAGACTTGTCCGACCTGTCATGGTTCTGGTCAGGTGCAGa                                 <  2:517620/267‑1 (MQ=255)
                                                                                                                                                                   gTTTTGGCGGCGGCGCAGACTTCAGCGATATTTTTGGTGACGTTTTCGGCGATATTTTTGGCGGCGGACGTGGGTCGTCAACGTGCGGCGCGCGGTGCTGATTTACGCTATAACATGGAGCTCACCCTCGAAGAAGCTGTACGTGGCGTGACCAAAGAGATCCGCATTCCGACTCTGGAAGAGTGTGACGTTTGCCACGGTAGCGGTGCAAAACCAGGTACACAGCCGCAGACTTGTCCGACCTGTCATGGTTATGGTCAGGTg                                     >  1:517623/1‑264 (MQ=255)
                                                                                                                                                                      ttGGCGGCGGCGCAGACTTCAGCGATATTTTTGGTGACGTTTTCGGCGATATTTTTGGCGGCGGACGTGGTCGTCAACGTGCGGCGCGCGGTGCTGATTTACGCTATAACATGGAGCTCACCCTCGAAGAAGCTGTACGt                                                                                                                                                             >  2:491038/1‑140 (MQ=255)
                                                                                                                                                                          cggcggcgCAGACTTCAGCGATATTTTTGGTGACGTTTTCGGCGATATTTTTGGCGGCGGACGTGGTCGTCAACGTGCGGCGCGCGGTGCTGATTTACGCTATAACATGGAGCTCACCCTCGAAGAAGCTGTACGTGGCg                                                                                                                                                         <  1:491041/140‑1 (MQ=255)
                                                                                                                                                                                            tgATATTTGTGGTGACGTGTTCGGCGATATTTTGGGCGGCGGACGTGGGTCGTCAACGTGCGGCGCGCGGTGCTGATTTACGCTATAACATGGAGCTCACCCTCGAAGAAGCTGTAAGTGGCGTGACCAAAGAGATCCGCATTGCGACTCTGGAAGAGTGTGACGTTTGCCACGGTAGCGGTGCAAAACCAGGTACACAGCCGCAGACCTGTCCGACCTGTCCTGGTTCTGGTCAGGTGCAGATGCGCCAGGGATTCTTCGCTGTACAGCAGAc  <  1:381328/273‑1 (MQ=255)
                                                                                                                                                                                            cGATATTTGTGGTGACGTTTTCGGCGATATTTTTGGCGGCGGACGTGGGTCGTCAACGTGCGGCGCGCGGTGCTGATTTACGCTATAACATGGAGCTCACCCTCGAAGAAGCTGTACGTGGCGTGACCAAAGAGATCCGCATTCCGACTCTGCAAGAGTGTGACGTTTGCCACg                                                                                                      >  1:725791/1‑174 (MQ=255)
                                                                                                                                                                                               tatTTTTGGTGACGTTTTCGGCGATATTTTTGGCGGCGGACGTGGGTCGTCAACGTGCGGCGCGCGGTGCTGATTTACGCTATAACATGGAGCTCa                                                                                                                                                                                 >  1:545477/1‑96 (MQ=255)
                                                                                                                                                                                                   tttGGTGACGTTTTCGGCGATATTTTTGGCGGCGGACGTGGGTCGTCAACGTGCGGCGCGCGGTGCTGATTTACGCTATAACATGGAGCTCACCCt                                                                                                                                                                             <  2:545474/96‑1 (MQ=255)
                                                                                                                                                                                                                                   cggACGTGGGTCGTCAACGTGCGGCGCGCGGTGCTGATTTACGCTATAACATGGAGCTCAccc                                                                                                                                                                              >  1:172576/1‑63 (MQ=255)
                                                                                                                                                                                                                                          |                                                                                                                                                                                                                                   
GCCATGAAATACCACCCGGACCGTAACCAGGGTGACAAAGAGGCCGAGGCGAAATTTAAAGAGATCAAGGAAGCTTATGAAGTTCTGACCGACTCGCAAAAACGTGCGGCATACGATCAGTATGGTCATGCTGCGTTTGAGCAAGGTGGCATGGGCGGCGGCGGTTTTGGCGGCGGCGCAGACTTCAGCGATATTTTTGGTGACGTTTTCGGCGATATTTTTGGCGGCGGACGTGGTCGTCAACGTGCGGCGCGCGGTGCTGATTTACGCTATAACATGGAGCTCACCCTCGAAGAAGCTGTACGTGGCGTGACCAAAGAGATCCGCATTCCGACTCTGGAAGAGTGTGACGTTTGCCACGGTAGCGGTGCAAAACCAGGTACACAGCCGCAGACTTGTCCGACCTGTCATGGTTCTGGTCAGGTGCAGATGCGCCAGGGATTCTTCGCTGTACAGCAGAC  >  CP009273/14252‑14712

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