Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	NC_000913	13,865	C→T	A568V (GCG→GTG)	dnaK →	chaperone Hsp70, with co‑chaperone DnaJ

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	13,865	0	C	T	100.0%	186.2 / NA	48	A568V (GCG→GTG)	dnaK	chaperone Hsp70, with co‑chaperone DnaJ
Reads supporting (aligned to +/- strand):  ref base C (0/0);  new base T (25/23);  total (25/23)

AGCTAACGCCGAAGCTGACCGTAAGTTTGAAGAGCTGGTACAGACTCGCAACCAGGGCGACCATCTGCTGCACAGCACCCGTAAGCAGGTTGAAGAAGCAGGCGACAAACTGCCGGCTGACGACAAAACTGCTATCGAGTCTGCGCTGACTGCACTGGAAACTGCTCTGAAAGGTGAAGACAAAGCCGCTATCGAAGCGAAAATGCAGGAACTGGCACAGGTTTCCCAGAAACTGATGGAAATCGCCCAGCAGCAACATGCCCAGCAGCAGACTGCCGGTGCTGAT  >  NC_000913/13722‑14007                                                                                                                                                |                                                                                                                                               aGCTAACGCCGAAGCTGACCGTAAGTTTGAAGAGCTGGTACAGACTCGCAACCAGGGCGACCATCTGCTGCACAGCACCCGTAAGCAGGTTGAAGAAGCAGGCGACAAACTGCCGGCTGACGACAAAACTGCTATCGAGTCTGTGCTGACt                                                                                                                                         >  2:195974/1‑151 (MQ=255)     aaCGCCGAAGCTGACCGTAAGTTTGAAGAGCTGGTACAGACTCGCAACCAGGGCGACCATCTGCTGCACAGCACCCGTAAGCAGGTTGAAGAAGCAGGCGACAAACTGCCGGCTGACGACAAAACTGCTATCGAGTCTGTGCTGACTGCAc                                                                                                                                     >  2:3483/1‑151 (MQ=255)      aCGCCGAAGCTGACCGTAAGTTTGAAGAGCTGGTACAGACTCGCAACCAGGGCGACCATCTGCTGCACAGCACCCGTAAGCAGGTTGAAGAAGCAGGCGACAAACTGCCGGCTGACGACAAAACTGCTATCGAGTCTGTGCTGACTGCACt                                                                                                                                    <  1:231824/151‑1 (MQ=255)        gCCGAAGCTGACCGTAAGTTTGAAGAGCTGGTACAGACTCGCAACCAGGGCGACCATCTGCTGCACAGCACCCGTAAGCAGGTTGAAGAAGCAGGCGACAAACTGCCGGCTGACGACAAAACTGCTATCGAGTCTGTGCTGACTGCACTgg                                                                                                                                  >  1:676514/1‑151 (MQ=255)          cGAAGCTGACCGTAAGTTTGAAGAGCTGGTACAGACTCGCAACCAGGGCGACCATCTGCTGCACAGCACCCGTAAGCAGGTTGAAGAAGCAGGCGACAAACTGCCGGCTGACGACAAAACTGCTATCGAGTCTGTGCTGACTGCACTGGa                                                                                                                                 >  1:356762/1‑150 (MQ=255)             aGCTGACCGTAAGTTTGAAGAGCTGGTACAGACTCGCAACCAGGGCGACCATCTGCTGCACAGCACCCGTAAGCAGGTTGAAGAAGCAGGCGACAAACTGCCGGCTGACGACAAAACTGCTATCGAGTCTGTGCTGACTGCACTGGAAACt                                                                                                                             <  2:415161/151‑1 (MQ=255)               cTGACCGTAAGTTTGAAGAGCTGGTACAGACTCGCAACCAGGGCGACCATCTGCTGCACAGCACCCGTAAGCAGGTTGAAGAAGCAGGCGACAAACTGCCGGCTGACGACAAAACTGCTATCGAGTCTGTGCTGACTGCACTGGAAACTGc                                                                                                                           >  1:768610/1‑151 (MQ=255)                       aaGTTTGAAGAGCTGGTACAGACTCGCAACCAGGGCGACCATCTGCTGCACAGCACCCGTAAGCAGGTTGAAGAAGCAGGCGACAAACTGCCGGCTGACGACAAAACTGCTATCGAGTCTGTGCTGACTGCACTGGAAACTGCTCTGaa                                                                                                                     >  1:357410/1‑149 (MQ=255)                       aaGTTTGAAGAGCTGGTACAGACTCGCAACCAGGGCGACCATCTGCTGCACAGCACCCGTAAGCAGGTTGAAGAAGCAGGCGACAAACTGCCGGCTGACGACAAAACTGCTATCGAGTCTGTGCTGACTGCACTGGAAACTGCTCTGAAAg                                                                                                                   <  1:176711/151‑1 (MQ=255)                             gAAGAGCTGGTACAGACTCGCAACCAGGGCGACCATCTGCTGCACAGCACCCGTAAGCAGGTTGAAGAAGCAGGCGACAAACTGCCGGCTGACGACAAAACTGCTATCGAGTCTGTGCTGACTGCACTGGAAACTGCTCTGAAAGGTGa                                                                                                               >  2:22047/1‑149 (MQ=255)                              aaGAGCTGGTACAGACTCGCAACCAGGGCGACCATCTGCTGCACAGCACCCGTAAGCAGGTTGAAGAAGCAGGCGACAAACTGCCGGCTGACGACAAAACTGCTATCGAGTCTGTGCTGACTGCACTGGAAACTGCTCTGAAAGGTGAAg                                                                                                             <  2:308175/150‑1 (MQ=255)                               agagCTGGTACAGACTCGCAACCAGGGCGACCATCTGCTGCACAGCACCCGTAAGCAGGTTGAAGAAGCAGGCGACAAACTGCCGGCTGACGACAAAACTGCTATCGAGTCTGTGCTGACTGCACTGGAAACTGCTCTGAAAGGTGAAGAc                                                                                                           >  2:104094/1‑151 (MQ=255)                                        aCAGACTCGCAACCAGGGCGACCATCTGCTGCACAGCACCCGTAAGCAGGTTGAAGAAGCAGGCGACAAACTGCCGGCTGACGACAAAACTGCTATCGAGTCTGTGCTGACTGCACTGGAAACTGCTCTGAAAGGTGAAGACAAAGCCGCt                                                                                                  <  1:3483/151‑1 (MQ=255)                                          aGACTCGCAACCAGGGCGACCATCTGCTGCACAGCACCCGTAAGCAGGTTGAAGAAGCAGGCGACAAACTGCCGGCTGACGACAAAACTGCTATCGAGTCTGTGCTGACTGCACTGGAAACTGCTCTGAAAGGTGAAGACAAAGCCGCTAt                                                                                                <  2:325678/151‑1 (MQ=255)                                           gACTCGCAACCAGGGCGACCATCTGCTGCACAGCACCCGTAAGCAGGTTGAAGAAGCAGGCGACAAACTGCCGGCTGACGACAAAACTGCTATCGAGTCTGTGCTGACTGCACTGGAAACTGCTCTGAAAGGTGAAGACAAAGCCGCTAt                                                                                                <  2:213217/150‑1 (MQ=255)                                                  aaCCAGGGCGACCATCTGCTGCACAGCACCCGTAAGCAGGTTGAAGAAGCAGGCGACAAACTGCCGGCTGACGACAAAACTGCTATCGAGTCTGTGCTGACTGCACTGGAAACTGCTCTGAAAGGTGAAGACAAAGCCGCTATCGAAGCg                                                                                         >  1:343527/1‑150 (MQ=255)                                                  aaCCAGGGCGACCATCTGCTGCACAGCACCCGTAAGCAGGTTGAAGAAGCAGGCGACAAACTGCCGGCTGACGACAAAACTGCTATCGAGTCTGTGCTGACTGCACTGGAAACTGCTCTGAAAGGTGAAGACAAAGCCGCTATCGAAGCg                                                                                         >  1:159273/1‑150 (MQ=255)                                                       gggCGACCATCTGCTGCACAGCACCCGTAAGCAGGTTGAAGAAGCAGGCGACAAACTGCCGGCTGACGACAAAACTGCTATCGAGTCTGTGCTGACTGCACTGGAAACTGCTCTGAAAGGTGAAGACAAAGCCGCTATCGAAGCGAAAATg                                                                                   <  1:745422/151‑1 (MQ=255)                                                             ccATCTGCTGCACAGCACCCGTAAGCAGGTTGAAGAAGCAGGCGACAAACTGCCGGCTGACGACAAAACTGCTATCGAGTCTGTGCTGACTGCACTGGAAACTGCTCTGAAAGGTGAAGACAAAGCCGCTATCGAAGCGAAAATGCAGGaa                                                                             <  2:7409/151‑1 (MQ=255)                                                             ccATCTGCTGCACAGCACCCGTAAGCAGGTTGAAGAAGCAGGCGACAAACTGCCGGCTGACGACAAAACTGCTATCGAGTCTGTGCTGACTGCACTGGAAACTGCTCTGAAAGGTGAAGACAAAGCCGCTATCGAAGCGAAAATGCAGGaa                                                                             <  1:695411/151‑1 (MQ=255)                                                              cATCTGCTGCACAGCACCCGTAAGCAGGTTGAAGAAGCAGGCGACAAACTGCCGGCTGACGACAAAACTGCTATCGAGTCTGTGCTGACTGCACTGGAAACTGCTCTGAAAGGTGAAGACAAAGCCGCTATCGAAGCGAAAATGCAGGaa                                                                             <  1:196506/150‑1 (MQ=255)                                                                  tgctgcACAGCACCCGTAAGCAGGTTGAAGAAGCAGGCGACAAACTGCCGGCTGACGACAAAACTGCTATCGAGTCTGTGCTGACTGCACTGGAAACTGCTCTGACAGGTGAAGACAAAGCCGCTATCGAAGCGAAAATGCAGGAACTGGc                                                                        >  1:484968/1‑151 (MQ=255)                                                                   gctgcACAGCACCCGTAAGCAGGTTGAAGAAGCAGGCGACAAACTGCCGGCTGACGACAAAACTGCTATCGAGTCTGTGCTGACTGCACTGGAAACTGCTCTGAAAGGTGAAGACAAAGCCGCTATCGAAGCGAAAATGCAGGAACTgg                                                                         <  1:569407/149‑1 (MQ=255)                                                                     tgcACAGCACCCGTAAGCAGGTTGAAGAAGCAGGCGACAAACTGCCGGCTGACGACAAAACTGCTATCGAGTCTGTGCTGACTGCACTGGAAACTGCTCTGAAAGGTGAAGACAAAGCCGCTATCGAAGCGAAAATGCAGGAACTGGcac                                                                      >  1:625805/1‑150 (MQ=255)                                                                     tgcACAGCACCCGTAAGCAGGTTGAAGAAGCAGGCGACAAACTGCCGGCTGACGACAAAACTGCTATCGAGTCTGTGCTGACTGCACTGGAAACTGCTCTGAAAGGTGAAGACAAAGCCGCTATCGAAGCGAAAATGCAGGAACTGGcac                                                                      >  2:124554/1‑150 (MQ=255)                                                                     tgcACAGCACCCGTAAGCAGGTTGAAGAAGCAGGCGACAAACTGCCGGCTGACGACAAAACTGCTATCGAGTCTGTGCTGACTGCACTGGAAACTGCTCTGAAAGGTGAAGACAAAGCCGCTATCGAAGCGAAAATGCAGGAACTGGcac                                                                      >  1:575918/1‑150 (MQ=255)                                                                         caGCACCCGTAAGCAGGTTGAAGAAGCAGGCGACAAACTGCCGGCTGACGACAAAACTGCTATCGAGTCTGTGCTGACTGCACTGGAAACTGCTCTGAAAGGTGAAGACAAAGCCGCTATCGAAGCGAAAATGCAGGAACTGGCACAGGtt                                                                 >  2:314626/1‑151 (MQ=255)                                                                               ccGTAAGCAGGTTGAAGAAGCAGGCGACAAACTGCCGGCTGACGACAAAACTGCTATCGAGTCTGTGCTGACTGCACTGGAAACTGCTCTGAAAGGTGAAGACAAAGCCGCTATCGAAGCGAAAATGCAGGAACTGGCACAGGTTTccc                                                             >  2:679417/1‑149 (MQ=255)                                                                                   aaGCAGGTTGAAGAAGCAGGCGACAAACTGCCGGCTGACGACAAAACTGCTATCGAGTCTGTGCTGACTGCACTGGAAACTGCTCTGAAAGGTGAAGACAAAGCCGCTATCGAAGCGAAAATGCAGGAACTGGCACAGGTTTCCCAGaaa                                                        <  1:527676/150‑1 (MQ=255)                                                                                           tGAAGAAGCAGGCGACAAACTGCCGGCTGACGACAAAACTGCTATCGAGTCTGTGCTGACTGCACTGGAAACTGCTCTGAAAGGTGAAGACAAAGCCGCTATCGAAGCGAAAATGCAGGAACTGGCACAGGTTTCCCAGAAACTGATGGa                                                >  1:600701/1‑150 (MQ=255)                                                                                                 agCAGGCGACAAACTGCCGGCTGACGACAAAACTGCTATCGAGTCTGTGCTGACTGCACTGGAAACTGCTCTGAAAGGTGAAGACAAAGCCGCTATCGAAGCGAAAATGCAGGAACTGGCACAGGTTTCCCAGAAACTGATGGAAATCGcc                                         >  2:183089/1‑151 (MQ=255)                                                                                                   cAGGCGACAAACTGCCGGCTGACGACAAAACTGCTATCGAGTCTGTGCTGACTGCACTGGAAACTGCTCTGAAAGGTGAAGACAAAGCCGCTATCGAAGCGAAAATGCAGGAACTGGCACAGGTTTCCCAGAAACTGATGGAAATCGCCca                                       <  1:568412/151‑1 (MQ=255)                                                                                                    aGGCGACAAACTGCCGGCTGACGACAAAACTGCTATCGAGTCTGTGCTGACTGCACTGGAAACTGCTCTGAAAGGTGAAGACAAAGCCGCTATCGAAGCGAAAATGCAGGAACTGGCACAGGTTTCCCAGAAACTGATGGAAATCGccc                                        >  1:319212/1‑149 (MQ=255)                                                                                                     ggCGACAAACTGCCGGCTGACGACAAAACTGCTATCGAGTCTGTGCTGACTGCACTGGAAACTGCTCTGAAAGGTGAAGACAAAGCCGCTATCGAAGCGAAAATGCAGGAACTGGCACAGGTTTCCCAGAAACTGATGGAAATCGCCcagc                                     <  1:93528/151‑1 (MQ=255)                                                                                                           caaCTGCCGCCTGGCGACAAAACTTATATCGACTCTGTGCTGATTGCACTGCAATCACGTCTGAAAGGCGAAGACAAATCCGCTATCGTACCGCAAATGCGGGAACTGGCACCGGCTGCCCACAAACTGATGGAAATCGCCCAgcagca                                 >  1:724441/2‑149 (MQ=255)                                                                                                             aCTGCCGGCTGACGACAAAACTGCTATCGAGTCTGTGCTGACTGCACTGGAAACTGCTCTGAAAGGTGAAGACAAAGCCGCTATCGAAGCGAAAATGCAGGAACTGGCACAGGTTTCCCAGAAACTGATGGAAATCGCCCAGCAGCAACa                              >  1:355729/1‑150 (MQ=255)                                                                                                                gCCGGCTGACGACAAAACTGCTATCGAGTCTGTGCTGACTGCACTGGAAACTGCTCTGAAAGGTGAAGACAAAGCCGCTATCGAAGCGAAAATGCAGGAACTGGCACAGGTTTCCCAGAAACTGATGGAAATCGCCCAGCAGCAACATg                            <  1:721005/149‑1 (MQ=255)                                                                                                                    gCTGACGACAAAACTGCTATCGAGTCTGTGCTGACTGCACTGGAAACTGCTCTGAAAGGTGAAGACAAAGCCGCTATCGAAGCGAAAATGCAGGAACTGGCACAGGTTTCCCAGAAACTGATGGAAATCGCCCAGCAGCAACATGCCcagc                      <  2:319245/151‑1 (MQ=255)                                                                                                                      tGACGACAAAACTGCTATCGAGTCTGTGCTGACTGCACTGGAAACTGCTCTGAAAGGTGAAGACAAAGCCGCTATCGAAGCGAAAATGCAGGAACTGGCACAGGTTTCCCAGAAACTGATGGAAATCGCCCAGCAGCAACATGCCcagc                      >  2:718711/1‑149 (MQ=255)                                                                                                                       gacgacAAAACTGCTATCGAGTCTGTGCTGACTGCACTGGAAACTGCTCTGAAAGGTGAAGACAAAGCCGCTATCGAAGCGAAAATGCAGGAACTGGCACAGGTTTCCCAGAAACTGATGGAAATCGCCCAGCAGCAACATGCCcagcag                    <  1:747281/150‑1 (MQ=255)                                                                                                                          gacAAAACTGCTATCGAGTCTGTGCTGACTGCACTGGAAACTGCTCTGAAAGGTGAAGACAAAGCCGCTATCGAAGCGAAAATGCAGGAACTGGCACAGGTTTCCCAGAAACTGATGGAAATCGCCCAGCAGCAACATGCCcagc                      >  1:330229/1‑145 (MQ=255)                                                                                                                           acAAAACTGCTATCGAGTCTGTGCTGACTGCACTGGAAACTGCTCTGAAAGGTGAAGACAAAGCCGCTATCGAAGCGAAAATGCAGGAACTGGCACAGGTTTCCCAGAAACTGATGGAAATCGCCCAGCAGCAACATGCCCAGCAGCAGAc               <  2:350719/151‑1 (MQ=255)                                                                                                                            cAAAACTGCTATCGAGTCTGTGCTGACTGCACTGGAAACTGCTCTGAAAGGTGAAGACAAAGCCGCTATCGAAGCGAAAATGCAGGAACTGGCACAGGTTTCCCAGAAACTGATGGAAATCGCCCAGCAGCAACATGCCCAGCAGCAGACt              <  1:182678/151‑1 (MQ=255)                                                                                                                               aaCTGCTATCGAGTCTGTGCTGACTGCACTGGAAACTGCTCTGAAAGGTGAAGACAAAGCCGCTATCGAAGCGAAAATGCAGGAACTGGCACAGGTTTCCCAGAAACTGATGGAAATCGCCCAGCAGCAACATGCCCAGCAGCAGACTGcc           >  2:607927/1‑151 (MQ=255)                                                                                                                                      aTCGAGTCTGTGCTGACTGCACTGGAAACTGCTCTGAAAGGTGAAGACAAAGCCGCTATCGAAGCGAAAATGCAGGAACTGGCACAGGTTTCCCAGAAACTGATGGAAATCGCCCAGCAGCAACATGCCCAGCAGCAGACTGCCGGTGCTg    <  2:683609/151‑1 (MQ=255)                                                                                                                                      aTCGAGTCTGTGCTGACTGCACTGGAAACTGCTCTGAAAGGTGAAGACAAAGCCGCTATCGAAGCGAAAATGCAGGAACTGGCACAGGTTTCCCAGAAACTGATGGAAATCGCCCAGCAGCAACATGCCCAGCAGCAGACTGCCGGTGCTg    <  1:625688/151‑1 (MQ=255)                                                                                                                                        cGAGTCTGTGCTGACTGCACTGGAAACTGCTCTGAAAGGTGAAGACAAAGCCGCTATCGAAGCGAAAATGCAGGACCTGGCACAGGTTTCCCAGAAACTGATGGAAATCGCCCAGCAGCAACATGCCCAGCAGCAGACTGCTGGCGCTGAt  >  1:640985/1‑151 (MQ=255)                                                                                                                                        cGAGTCTGTGCTGACTGCACTGGAAACTGCTCTGAAAGGTGAAGACAAAGCCGCTATCGAAGCGAAAATGCAGGAACTGGCACAGGTTTCCCAGAAACTGATGGAAATCGCCCAGCAGCAACATGCCCAGCAGCAGACTGCCGGTGCTGAt  <  2:301122/151‑1 (MQ=255)                                                                                                                                                |                                                                                                                                               AGCTAACGCCGAAGCTGACCGTAAGTTTGAAGAGCTGGTACAGACTCGCAACCAGGGCGACCATCTGCTGCACAGCACCCGTAAGCAGGTTGAAGAAGCAGGCGACAAACTGCCGGCTGACGACAAAACTGCTATCGAGTCTGCGCTGACTGCACTGGAAACTGCTCTGAAAGGTGAAGACAAAGCCGCTATCGAAGCGAAAATGCAGGAACTGGCACAGGTTTCCCAGAAACTGATGGAAATCGCCCAGCAGCAACATGCCCAGCAGCAGACTGCCGGTGCTGAT  >  NC_000913/13722‑14007

Alignment Legend
Aligned base mismatch/match (shaded by quality score): ATCG/ATCG < 3 ≤ ATCG/ATCG < 16 ≤ ATCG/ATCG < 33 ≤ ATCG/ATCG < 38 ≤ ATCG/ATCG < 39 ≤ ATCG/ATCG
Unaligned base: atcg    Masked matching base: atcg    Alignment gap: ‑    Deleted base: ‑