Predicted mutation
evidence seq id position mutation freq annotation gene description
RA ACB122 5,508,974 C→G 18.0% F46L (TTC→TTG)  dnaK_2 → Chaperone protein DnaK

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*ACB1225,508,9740CG18.0% 128.2 / 5.2 57F46L (TTC→TTG) dnaK_2Chaperone protein DnaK
Reads supporting (aligned to +/- strand):  ref base C (11/30);  new base G (8/2);  total (23/34)
Fisher's exact test for biased strand distribution p-value = 3.17e-03
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

TATCCCCGGCCCGCGCCCTGGGCATCGACTTTGGCACCTCCAACTCCACAGTCGGTTGG‑CATCGGCCTGGCGTCGAGTCGCTG‑ATTGCCCTGGAAGACGGCAAGATCACCTTGCCCTCGGTGGTGTTCTTCAATATCGAAGAGCGGCGCCCGGTGTATGGCCGCCTGGCGCTGCACGAGTACCTGGAAGGCTACGAAGGCCGCCTGATGCGCTCGCTGAAAAGCCTGCTGGGCTCGAAGCTGATCAAGCACGACACCAGC  >  ACB122/5508847‑5509106
                                                                                                                                 |                                                                                                                                    
tATCCCCGGCCCGCGCCCTGGGCATCGACTTTGGCACCTCCAACTCCACAGTCGGTTGG‑CATCGGCCTGGCGTCGAGTCGCGG‑ATTGGCCGGGAAGACGGCCAGATCACCTTGCCCGCGGGGGGGTTGTTGAat                                                                                                                                >  7:249432/1‑134 (MQ=255)
 aTCCCCGGCCCGCGCCCTGGGCATCGACTTTGGCACCTCCAACTCCACAGTCGGTTGG‑CATCGGCCTGGCGTCGAGTCGCTG‑ATTG‑CCTGGAAGACGGCAAGATCACCGTGGCCGCGGGGGTGTTGGTGAatatg                                                                                                                             >  7:272306/1‑134 (MQ=255)
     ccGGCCCGCGCCCTGGGCATCGACTTTGGCACCTCCAACTCCACAGTCGGTTGG‑CATCGGCCTGGCGTCGAGTCGCTG‑ATTGCCCTGGAAGACCGCAAGATCCCCGTGCCCGCGGTGGGGTGTTTGAATATc                                                                                                                             >  7:173919/1‑132 (MQ=255)
     ccGGCCCGCGCCCTGGGCATCGACTTTGGCACCTCCAACTCCACAGTCGGTTGG‑CATCGGCCTGGCGTCGAGTCGCGG‑GTTTGCCGGGGAGCCCGCCAGCTCACCCTGCCCGCGGGGGTGTTGGTGAATATAGa                                                                                                                           >  7:129435/1‑134 (MQ=255)
     ccGGCCCGCGCCCTGGGCATCGACTTTGGCACCTCCAACTCCACAGTCGGTTGG‑CATCGGCCTGGCGTCGAGTCGCGG‑GTTGCCCTGGCAGACGGCCAGATCCCCTTGCCCGCGGGGGGGGTCGTGAATATCGat                                                                                                                          >  4:205657/1‑134 (MQ=255)
      cggccggcGCCGTGGGCAGCGACTTTGGCACCTCCAACTCCACAGTCGGTTGG‑CATCGGCCTGGCGTCGAGTCGCTG‑ATTGCCCTGGAAGACGGCAAGATCACCTTGCCCTCGGTGGTGTTCTTCAATATCGAag                                                                                                                         <  2:107976/135‑1 (MQ=255)
      cGGCCCGCGCCCTGGGCATCGACTTTGGCACCTCCAACTCCACAGTCGGTTGG‑CATCGGCCTGGCGTCGAGTCGCTG‑CTTGGCCTGGGAGACGGCCAGGTGACCGTGGCCGGGGTGGGGGGGTTCAatatggg                                                                                                                           >  1:174070/1‑130 (MQ=255)
       ggCCCGCGCCCTGGGCATCGACTTTGGCACCTCCAACTCCACAGTCGGTTGG‑CATCGGCCTGGCGTCGAGTCGCTG‑ATTGCCCTGGAAGACGGCAAGATCACCTTGCCCTCGGTGGTGTTCTTCAATATCGAaga                                                                                                                        <  8:284907/135‑1 (MQ=255)
            gcgcCGTGGGCATCTACTTTGGCACCTCCAACTCCACAGGCGGTTGG‑CATCGGCCTGGCGTCGAGTGGCTG‑ATTGCCCTGGAAGGCGGCAAGATCACATTGCCCTCGGTGTTGTTGTTCAATATCGAATAgcggc                                                                                                                   <  2:283046/135‑1 (MQ=255)
            gcgcCCTGGGCATCGACTTTGGCACCTCCAACTCCACAGTCGGTTGG‑CATCGGCCTGGCGTCGAGTCGCTG‑ATTGCCCTGGAAGACGGCAAGATCACCTTGCCCTCGGTGGTGTTCTTCAATATCGATGAgcggc                                                                                                                   <  3:205657/135‑1 (MQ=255)
            gcgcCCTGGGCATCGACTTTGGCACCTCCAACTCCACAGTCGGTTGG‑CATCGGCCTGGCGTCGAGTCGCTG‑ATTGCCCTGGAAGACGGCAAGATCACCTTGCCCTCGGTGGTGTTCTTCAATATCGAAGAgcggc                                                                                                                   <  5:153973/135‑1 (MQ=255)
            gcgcCCTGGGCATCGACTTTGGCACCTCCAACTCCACAGTCGGTTGG‑CATCGGCCTGGCGTCGAGTCGCTG‑ATTGCCCTGGAAGACGGCAAGATCACCTTGCCCTCGGTGGTGTTCTTCAATATCGAAGAgcggc                                                                                                                   <  6:114155/135‑1 (MQ=255)
            gcgcCCTGGGCATCGACTTTGGCACCTCCAACTCCACAGTCGGTTGG‑CATCGGCCTGGCGTCGAGTCGCTG‑ATTGCCCTGGAAGACGGCAAGATCACCTTGCCCTCGGTGGTGTTCTTCAATATCGAAGAgcggc                                                                                                                   <  5:64885/135‑1 (MQ=255)
            gcgcCCTGGGCATCGACTTTGGCACCTCCAACTCCACAGTCGGTTGG‑CATCGGCCTGGCGTCGAGTCGCTG‑ATTGCCCTGGAAGACGGCAAGATCACCTTGCCCTCGGTGGTGTTCTTCAATATCGAAGAgcggc                                                                                                                   <  3:42684/135‑1 (MQ=255)
            gcgcCCTGGGCATCGACTTTGGCACCTCCAACTCCACAGTCGGTTGG‑CATCGGCCTGGCGTCGAGTCGCTG‑ATTGCCCAGGAAGACCGCAAGCTAACCCTTCCCGGGGTTGGGTTATTGAATATCGCAGAgccgc                                                                                                                   >  8:5704/1‑135 (MQ=255)
            gcgcCCTGCGCATCGACTTTGGCACCTCCAACTCCACAGGCGGTTGG‑CATCGGCCGGGCGTCGAGTCGCTG‑ATTGCCCTGGAAGACGGCAAGATCACCTTGCCCTCGGTGGTGTTCTTCAATATCGAAGAgcggc                                                                                                                   <  6:221035/135‑1 (MQ=255)
            gcgcCCGGGGCAGGGACTTTGGCACCTCCAACTCCACTGT‑GGTTGGCCAGCGGCGTGGGGTCGAGTCGGTG‑ATTCCCCTGGAAGACGGCAAGATCAGCTTGCCCTCGGTGGTTTTCTTCAATATCGAAGAgcggc                                                                                                                   <  2:137600/135‑1 (MQ=255)
               cccTGGGCATCGACTTTGGCACCTCCAACTCCACAGTCGGTTGG‑CATCGGCCTGGCGTCGAGTCGCTG‑ATTGCCCTGGAAGACGGCAAGATCACCTTGCCCTCGGTGGTGTTCTTCAATATCGAAGAGCGGCGcc                                                                                                                <  6:117496/135‑1 (MQ=255)
                ccagggCATCGACTTTGGCACCTCCCACTCCACAGTCGGTTGG‑CCTCGGCCTGGCGTCGAGTCGGGGCTTTGCCGTGCACG‑CCGCAAGATCACATTGCGCGCGGGGGTGTTGTTGAAGAACGAAGCGCGGCcccc                                                                                                               >  8:51055/4‑135 (MQ=255)
                ccTGGGCATCGACTTTGGCACCTCCAACTCCACAGTCGGTTGG‑CATCGGCCTGGCGTCGAGTCGCTG‑GTTGGCCTGGTAGACGGCAAGGAAACCGTTCCCACGGGGGGGTTGTTTAATTTTGAAGAGCGGCGgcc                                                                                                               >  5:292938/1‑135 (MQ=255)
                ccTGGGCATCGACTTTGGCACCTCCAACTCCACAGTCGGTTGG‑CATCGGCCTGGCGTCGAGTCGCTG‑ATTGCCCTGGAAGACGGCAAGATAACCTTGCCCTCGGTGGGGTTATTCAATATCGAAGAGCGGGGccc                                                                                                               >  8:153724/1‑135 (MQ=255)
                ccTGGGCATCGACTTTGGCACCTCCAACTCCACAGTCGGTTGG‑CATCGGCCTGGCGTCGAGTCGCTG‑ATTGCCCTGGAAGACCGCAAGATCACCTTGGCCGCGGTGGTGTTGGTGCATATAGAAGAGCGGCGccc                                                                                                               >  4:292539/1‑135 (MQ=255)
                      cATCGACTTTGGCACCTCCAACTCCACAGTCGGTTGG‑CATCGGCCTGGCGTCGAGTCGCTG‑ATTGCCCGGGAAGACGGCAAGATCACCATTCGCTCGGTGGTGTTCTTCAATATCGAAGAGCGGCGCCCGGGGTa                                                                                                         >  7:243312/1‑135 (MQ=255)
                           aCTTTGGCACCTCCAACTCCACAGTCGGTTGG‑CATCGGCCTGGCGTCGAGTCGCTG‑ATTGCCCTGGAAGACGGCAAGATCACCTTGCCCTCGGTGGTGTTCTTCAATATCGAAGAGCGGCGCCCGGTGTATGgc                                                                                                     <  7:153724/134‑1 (MQ=255)
                            cTTTGGCACCTCCAACTCCACAGTCGGTTGG‑CATCGGCCTGGCGTCGAGTCGCTG‑ATTGCCCTGGAAGACGGCAAGATCACCTTGCCCTCGGTGGTGttct                                                                                                                                     <  7:225267/101‑1 (MQ=255)
                            cTTTGGCACCTCCAACTCCACAGTCGGTTGG‑CATCGGCCTGGCGTCGAGTCGCTG‑ATTGCCCTGGAAGACGGCAAGATCACCTTGCCCTCGGTGGTGttct                                                                                                                                     >  8:225267/1‑101 (MQ=255)
                            cTTTGGCACCTCCAACTCCACAGTCGGTTGG‑CATCGGCCTGGCGTCGAGTCGCGG‑ATTGCCCTGGAAGACGGCAAGCTAACCGTGCCCTCGGTGGTGTTCTTAAATATCGAAGAGCGGCGCCCCGGGTATGgccg                                                                                                   >  1:9449/1‑135 (MQ=255)
                                      tCCAACTCCACAGGCGGTTGG‑CATCGGCCTGGCGTCGAGTCGCTG‑ATTGCCCTGGAAGACGGCAAGATCACCTTGCCCTCGGTGGTGTTCTTCAATATCGAAGAGCGGCGCCCGGTGTATGGCCGCCTGGCGCTg                                                                                         <  7:41028/135‑1 (MQ=255)
                                        caccTCCCCAGTCGGGTGG‑CAGCGGCCTGGCGTCGAGGCGCTG‑ATTGCCCTGGAAGACGGCAAGATCACCTTGCCCTCGGTGGGGTTCTTCAAAATCGAAGAGCGGCGCCCGGTGTATGGCCGCCTGGCGCTGCa                                                                                       <  6:6952/132‑1 (MQ=255)
                                          aCTCCACAGTCGGTTGG‑CATCGGCCTGGCGTCGAGTCGCTG‑CTTGCCCTGGAAGACGGCAAGATGACCGTGCCCTCGGTGGTGTTCGTCAATATAGAAGAGCGGCGCCCCGTGGTTTGCCGCCTGGGGCTGCa                                                                                       >  2:215052/1‑133 (MQ=255)
                                                 aGTCGGTTGG‑CATCGGCCTGGCGTCGAGTCGCTG‑ATTGCCCTGGAAGACGGCAAGATCACCTTGCCCTCGGTGGTGTTGTTCAATATCGAAGAGCGGCGCCCGGTGTATGGCCGCCTGGCGCTGCACGCGTAccg                                                                              >  4:214427/1‑134 (MQ=255)
                                                  gTCGGTTGG‑CATCGGCCTGGCGTCGAGTCGCGG‑CATGCCCGGGAACACGGCCAGAACACCATGCCCGAGGTGGGGGTCTTCAATATCGCAGAGCGGCGCCCGGGGTAATGCCGCCCGGCGCGGGCCGCGTAccag                                                                             >  6:23722/1‑133 (MQ=255)
                                                   tCGGTTGG‑CATCGGCCTGGCGTCGAGTCGCTG‑ATTTCCCTGGAAGACCGCAATAACACCTTGCCCGCGGTGGTGTTATTCAATATCGAAGAGCGGCCCCCGGTGTATGGCCGCCTGGCGCTGCACGAGTACCTg                                                                             >  8:223036/1‑134 (MQ=255)
                                                   tCGGTTGG‑CATCGGCCTGGCGTCGAGTCGCTG‑ATTGCCCTGGAAGACGGCAAGATCACCTTGCCCTCGGTGGTGTTCTTCAATATCGAAGAGCGGCGCCCGGTGTATGGCCGCCTGGCGCTGCACGAGTACCTgg                                                                            <  6:121337/135‑1 (MQ=255)
                                                   tCGGTTGG‑CATCGGCCTGGCGTCGAGTCGCTG‑ATTGCCCTGGAAGACGGCAAGATCACCTTGCCCTCGGTGGTGTTCTTCAATATCGAAGAGCGGCGCCCGGTGTATGGCCGCCGGGCGCTGCACGAGTACCTgg                                                                            <  2:79248/135‑1 (MQ=255)
                                                   tCGGTTGG‑CATCGGCCTGGCGTCGAGGCGCTG‑ATTGCCCTGGAAGACGGCAAGATCACCTTGCCCTCGGTGGTGTTCTTCAATATCGAAGAGCGGCGCCCGGTGTATGGCCGCCTGGCGCTGCACGAGTACCTgg                                                                            <  8:158128/135‑1 (MQ=255)
                                                       ttGG‑CAGCGGCCTGGCGTCGAGTCGCTG‑ATTGCCCTGGAAGACGGCAAGATCACCTTGCCCTCGGTGGTGTTCTTCAATATCGAAGAGCGGCGCCCGGTGTATGGCCGCCTGGCGCTGCACGAGTACCTGGAAgg                                                                        <  5:158649/135‑1 (MQ=255)
                                                        tGG‑CATCGGCCTGGCGTCGAGTCGCTG‑ATTGCCCTGGAAGACGGCAAGATCACCTTGCCCTCGGTGGTGTTCTTCAATATCGAAGAGCGGCGCCCGGTGTATGGCCGCCTGGCGCTGCACGAGTACCTGGAAGGc                                                                       <  5:220677/135‑1 (MQ=255)
                                                            cATCGGCCTGGCGTCGAGTCGCTG‑ATTGCCCTGGAAGACGGCAAGATCACCTTGCCCGCGGTGGTGTTCTTCAATATCGAAGAGCGGCGCCCGGTGTATGGCCGCCTGGCGCTGCACGCGTACCTGGAAGGCTAc                                                                    >  3:141355/1‑135 (MQ=255)
                                                             aTCGGCCTGGCGTCGAGTCGCTG‑ATTGCCCTGGAAGACGGCAAGATCACCTTGCCCTCGGTGGTGTTCTTCAATATCGAAGAGCGGCGCCCGGTGTATGGCCGCCTGGCGCTGCACGAGTACCTGGAAGGCTACg                                                                   <  8:88597/135‑1 (MQ=255)
                                                             aTCGGCCTGGCGTCGAGTCGCTG‑ATTGCCCTGGAAGACGGCAAGATCACCTTGCCCTCGGTGGTGTTCTTCAATATCGAAGAGCGGCGCCCGGTGTATGGCCGCCTGGCGCTGCACGAGTACCTGGAAGGCTACg                                                                   <  7:45555/135‑1 (MQ=255)
                                                             aTCGGCCTGGCGTCGAGGCGCTG‑ATTGCCCTGGAAGACGGCAAGATCACCTTGCCCTCGGTGGTGTTCTTCAATATCGAAGAGCGGCGCGCGGTGTATGGCCGCCTGGCGCTGCACGAGTACCTGGAAGGCTACg                                                                   <  8:20131/135‑1 (MQ=255)
                                                                                  tG‑ATTGCCCTGGAAGACGGCAAGATCACCTTGCCCTCGGTGGTGTTCTTCAATATCGAAGAGCGGCGCCCGGTGTATGGCCGCCTGGCGCTGCACGAGTACCTGGAAGGCTACGAAGGCCGCCTGATGcgctcgc                                              <  5:59920/135‑1 (MQ=255)
                                                                                   g‑ATTGCCCTGGAAGACGGCAAGATCACCTTGCCCTCGGTGGTGTTCTTCAATATCGAAGAGCGGCGCCCGGTGTATGGCCGCCTGGCGCTGCACGCGTACCTGGAAGGCTACGAAGGCCGCCTGATGcgctcgct                                             >  3:275155/1‑135 (MQ=255)
                                                                                               aaGACGGCAAGATCACCTTGCCCTCGGTGGTGTTCTTCAATATCGAAGAGCGGCGCCCGGTGTATGGCCGCCTGGCGCTGCACGAGTACCTGGAAGGCTACGAAGGCCGCCTGATGCGCTCGCTGAAAAGCct                                    <  1:251537/133‑1 (MQ=255)
                                                                                                       aaGATCACCTTGCCCTCGGTGGTGTTCTTCAATATCGAAGAGCGGCGCCCGGTGTATGGCCGCCTGGCGCTGCACGAGTACCTGGAAGGCTACGAAGGCCGCCTGATGCGCTCGCTGAAAAGCCTGCTGGGCTc                           >  7:49925/1‑134 (MQ=255)
                                                                                                       aaGATCACCTTGCCCTCGGTGGTGTTCTTCAATATCGAAGAGCGGCGCCCGGTGTATGGCCGCCTGGCGCTGCACGAGTACCTGGAAGGCTACGAAGGCCGCCTGATGCGCTCGCTGAAAAGCCTGCTGGGCTCg                          <  8:173919/135‑1 (MQ=255)
                                                                                                              ccTGGCCCGCGGGGGTGTTATTCAATATCGAAGAGCGGCGCCGGTTGTTTGGCCGCCTGGCGCTGCACGAGTACCTGGAAGGCTACGAAGGCCGCCTGATGCGCTCGCTGAAAAGCCTGCTGGGCTCGAAGCTGa                   <  8:129435/135‑1 (MQ=255)
                                                                                                                   cccTCGGTGGTGTTGTTCAATAGCGAAGAGCGGCGCCCGGTGTATGGCCGCCTGGCGCGGCACGAGTACCTGGAAGGCTACGAAAGCCGCGTGATGCGCTCGCTGAAAAGCCTGCTGGGCTCGAAGCTGAGCaa               <  6:284722/134‑1 (MQ=255)
                                                                                                                   cccTCGGTGGTGTTCTTCAATATCGAAGAGCGGCGCCCGGTGTATGGCCGCCTGGCGCTGCACGAGTACCTGGAAGGCTACGAAGGCCGCCTGATGCGCTCGCTGAAAAGCCTGCTGGGCTCGAAGCTGATCaa               <  6:292938/134‑1 (MQ=255)
                                                                                                                   cccTCGGTGGTGTTCTTCAATATCGAAGAGCGGCGCCCGGTGTATGGCCGCCTGGCGCTGCACGAGTACCTGGAAGGCTACGAAGGCCGCCTGATGCGCTCGCTGAAAAGCCTGCTGGGCTCGAAGCTGATCaa               <  8:49925/134‑1 (MQ=255)
                                                                                                                   cccTCGGGGGTGTTAGTCAATCTCGAAGAGCGGCGCCCGGTGTATGGCCGCCTGGCGCTGCACGAGTACCTGGACGGCTACGAAGGCCGCCTGATGCGCTCGCTGAAAAGCCTGCTGGGCTCGAAGCTGATCaa               <  5:23722/134‑1 (MQ=255)
                                                                                                                       cGGTGGTGTTCTTCAATATCGAAGAGCGGCGCCCGGTGTa                                                                                                         <  4:278915/40‑1 (MQ=255)
                                                                                                                       cGGTGGTGTTCTTCAATATCGAAGAGCGGCGCCCGGTGTa                                                                                                         >  3:278915/1‑40 (MQ=255)
                                                                                                                          tggtggTCTTCAATATCGCAGAGCGGCGCCCGGTGTATGGCCGCCTGGCGCTGCACGAGTACCTGGAAGGCTACGAAGGCCGCCTGATGCGCTCGCTGAAAAGCCTGCTGGGCTCGAAGCTGATCAAGCACGaca       <  8:243312/135‑1 (MQ=255)
                                                                                                                          tggtgTTCTTCAATATCGAAGAGCGGCGCCCGGTGTATGGCCGCCTGGCGCTGCACGAGTACCTGGAAGGCTACGAAGGCCGCCTGATGCGCTCGCTGAAAAGCCTGCTGGGCTCGAAGCTGATCAAGCACGcca       >  5:284245/1‑135 (MQ=255)
                                                                                                                          tggtgTTCTTCAATATCGAAGAGCGGCGCCCGGTGGATGGCCGCCTGGCGCTGCACGAGTACCTGGAAGGCTACGAAGGCCGCCTGATGCGCTCGCTGAAAAGCCTGCTGGGCTCGAAGCTGATCAAGCACGaca       <  4:277576/135‑1 (MQ=255)
                                                                                                                          gggtgTTCTTCAATCTCTCTGAGCGGCGCCCGGTGTATGGCCGCCTGGCGCTGCACGAGTACCTGGAAGGCTACTAAGGCCGCCTTATGCGCTCGCTGAAAAGCCTGCTGGGCTCGAAGCTGATCAAGCACGaca       <  8:124160/134‑1 (MQ=255)
                                                                                                                           ggtgTTCTTCAATATCGAAGAGCGGCGCCCGGTGTATGGCCGCCTGGCGCTGCACGAGTACCTGGAAGGCTACGAAGGCCGCCTGATGCGCTCGCTGAAAAGCCTGCTGGGCTCGAAGCTGATCAAGCACGaca       >  3:126373/1‑134 (MQ=255)
                                                                                                                               ttcttcAATATCGAAGAGCGGCGCCCGGTGTCTGGCCGCCGGGCGCTGCACGAGTACCTGGAAGGCTACGAAGGCCGCCTGATGCGCTCGCTGAAAAGCCTGCTGGGCTCGAAGCTGATCAAGCACGACACCAGc  <  4:222984/135‑1 (MQ=255)
                                                                                                                                 |                                                                                                                                    
TATCCCCGGCCCGCGCCCTGGGCATCGACTTTGGCACCTCCAACTCCACAGTCGGTTGG‑CATCGGCCTGGCGTCGAGTCGCTG‑ATTGCCCTGGAAGACGGCAAGATCACCTTGCCCTCGGTGGTGTTCTTCAATATCGAAGAGCGGCGCCCGGTGTATGGCCGCCTGGCGCTGCACGAGTACCTGGAAGGCTACGAAGGCCGCCTGATGCGCTCGCTGAAAAGCCTGCTGGGCTCGAAGCTGATCAAGCACGACACCAGC  >  ACB122/5508847‑5509106

Alignment Legend
Aligned base mismatch/match (shaded by quality score): ATCG/ATCG < 3 ≤ ATCG/ATCG < 14 ≤ ATCG/ATCG < 15 ≤ ATCG/ATCG < 27 ≤ ATCG/ATCG < 36 ≤ ATCG/ATCG
Unaligned base: atcg    Masked matching base: atcg    Alignment gap: ‑    Deleted base: ‑