Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	NC_000913	13,436	T→C	V425A (GTG→GCG)	dnaK →	chaperone protein DnaK

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	13,436	0	T	C	100.0%	170.8 / NA	49	V425A (GTG→GCG)	dnaK	chaperone protein DnaK
Reads supporting (aligned to +/- strand):  ref base T (0/0);  new base C (25/24);  total (25/24)

TGGTGTTCTGACTGGTGACGTAAAAGACGTACTGCTGCTGGACGTTACCCCGCTGTCTCTGGGTATCGAAACCATGGGCGGTGTGATGACGACGCTGATCGCGAAAAACACCACTATCCCGACCAAGCACAGCCAGGTGTTCTCTACCGCTGAAGACAACCAGTCTGCGGTAACCATCCATGTGCTGCAGGGTGAACGTAAACGTGCGGCTGATAACAAATCTCTGGGTCAGTTCAACCTAGATGGTATCAACCCGGCACCGCGCGGCATGCC  >  NC_000913/13299‑13571                                                                                                                                          |                                                                                                                                        tggtgTTCTGACTGGTGACGTAAAAGACGTACTGCTGCTGGACGTTACCCCGCTGTCTCTGGGTATCGAAACCATGGGCGGTGTGATGACGACGCTGATCGCGAAAAACACCACTATCCCGACCAAGCACAGCCAGgcgt                                                                                                                                       >  2:118022/1‑140 (MQ=255)    tgtTCTGACTGGTGACGTAAAAGACGTACTGCTGCTGGACGTTACCCCGCTGTCTCTGGGTATCGAAACCATGGGCGGTGTGATGACGACGCTGATCGCGAAAACCACGACTATCCCAACCATGCACAGCCAGGCGTtct                                                                                                                                    >  2:235896/1‑140 (MQ=255)    tgtTCTGACTGGTGACGTAAAAGACGTACTGCTGCTGGACGTTACCCCGCTGTCTCTGGGTATCGAAACCATGGGCGGTGTGATGACGACGCTGATCGCGAAAAACACCACTATCCCGACCAAGCACAGCCAGGCGTtct                                                                                                                                    >  2:420255/1‑140 (MQ=255)           aCTGGTGACGTAAAGTACGTGCAGCTGCTGGACGTTACCCCGCTGTCTCTGGGTATCGAAACCATGGGCGGTGTGATGACGACGCTGATCGCGAAAAACACCACTATCCCGACCAAGCACAGCCAGGCGTTCTCTACCGc                                                                                                                             <  1:286846/140‑1 (MQ=255)           aCTGGTGACGTAAAAGACGTACTGCTGCTGGACGTTACCCCGCTGTCTCTGGGTATCGAAACCATGGGCGGTGTGATGACGACGCTGATCGCGAAAAACATCACCATCACGACCAAGCACAGCCAGGCGTTCTCTACCGc                                                                                                                             >  2:286846/1‑140 (MQ=255)           aCTGGTGACGTAAAAGACGTACTGCTGCTGGACGTTACCCCGCTGTCTCTGGGTATCGAAACCATGGGCGGTGTGATGACGACGCTGATCGCGAAAAACACCACTATCCCGACCAAGCACAGCCAGGCGTTCTCTACCGc                                                                                                                             <  2:4028/140‑1 (MQ=255)           aCTGGTGACGTAAAAGACGTACTGCTGCTGGACGTTACCCCGCTGTCTCTGGGTATCGAAACCATGGGCGGTGTGATGACGACGCTGATCGCGAAAAACACCACTATCCCGACCAAGCACAGCCAGGCGTTCTCTACCGc                                                                                                                             <  2:435119/140‑1 (MQ=255)               gTGACGTAAAAGACGTACTGCTGCTGGACGTTACCCCGCTGTCTCTGGGTATCGAAACCATGGGCGGTGTGATGACGACGCTGATCGCGAAAAACACCACTATCCCGACCAAGCACAGCCAGGCGTTCTCTACCGCTGaa                                                                                                                         <  1:6213/140‑1 (MQ=255)               gTGACGTAAAAGACGTACTGCTGCTGGACGTTACCCCGCTGTCTCTGGGTATCGAAACCATGGGCGGTGTGATGACGACGCTGATCGCGAAAAACACCACTATCCCGACCAAGCACAGCCAGGCGTTCTCTACCGCTGaa                                                                                                                         <  1:555493/140‑1 (MQ=255)                tGACGTAAAAGACGTACTGCTGCTGGACGTTACCCCGCTGTCTCTGGGTATCGAAACCATGGGCGGTGTGATGACGACGCTGATCGCGAAAAACACCACTATCCCGACCAAGCACAGCCAGGCGTTCTCTACCGCTGAAg                                                                                                                        >  2:193804/1‑140 (MQ=255)                             gacaggctgctgGACGTCACCCCGCTGTCTCTGGGAATCTAAACCATGCGCGGTGTGATGACGACGCTGATCGCGAAAAACACCACTATCCCGACCAAGCACAGCCAGGCGTTCTCTACCGCTGAAGACAACCAGTCTGc                                                                                                           <  2:270713/135‑1 (MQ=255)                                   ctgctgGACGTTACCCCGCTGTCTCTGGGTATCGAAACCATGGGCGGTGTGATGACGACGCTGATCGCGAAAAACACCACTATCCCGACCAAGCACAGCCAGGCGTTCTCTACCGCTGAAGACAACCAGTCTGCGGTAAc                                                                                                     <  2:176514/140‑1 (MQ=255)                                     gctgGACGTTACCCCGCTGTCTCTGGGTATCGAAACCATGGGCGGTGTGATGACGACGCTGATCGCGAAAAACACCACTATCCCGACCAAGCACAGCCAGGCGTTCTCTACCGCTGAAGACAACCAGTCTGCGGTAAcca                                                                                                   <  1:118022/140‑1 (MQ=255)                                      ctgGACGTTACCCCGCTGTCTCTGGGTATCGAAACCATGGGCGGTGTGATGACGACGCTGATCGCGAAAAACACCACTATCCCGACCAAGCACAGCCAGGCGTTCTCTACCGCTGAAGACAACCAGTCTGCGGTAAccat                                                                                                  <  1:181035/140‑1 (MQ=255)                                        ggACGTTACCCCGCTCTCTCTGCGTCTCGCAACCATGCGCGGTGTGATGACGACGCTGATCGCGAAAAACACCACTATCCCGACCAAGCACAGCCAGGCGTTCTCTACCGCTGAAGACAACCAGTCTGCGGTAAccatcc                                                                                                <  1:224060/140‑1 (MQ=255)                                         gACGTTACCCCGCTGTCTCTGGGTATCGAAACCATGGGCGGTGTGATGACGACGCTGATCGCGAAAAACACCACTATCCCGACCAAGCACAGCCAGGCGTTCTCTACCGCTGAAGACAACCAGTCTGc                                                                                                           <  2:697085/128‑1 (MQ=255)                                         gACGTTACCCCGCTGTCTCTGGGTATCGAAACCATGGGCGGTGTGATGACGACGCTGATCGCGAAAAACACCACTATCCCGACCAAGCACAGCCAGGCGTTCTCTACCGCTGAAGACAACCAGTCTGc                                                                                                           >  1:697085/1‑128 (MQ=255)                                         gACGTTACCCCGCTGTCTCTGGGTATCGAAACCATGGGCGGTGTGATGACGACGCTGATCGCGAAAAACACCACTATCCCGACCAAGCACAGCCAGGCGTTCTCTACCGCTGAAGACAACCAGTCTGCGGTAAccatcca                                                                                               >  1:685575/1‑140 (MQ=255)                                          aCGTTGCCCCGCTGTCTCTGGGTATCGAAACCATGGGCGGTGTGATGACGACGCTGATCGCGAAAAACACCACTATCCCGACCAAGCACAGCCAGGCGTTCTCTACCGCTGAAGACAACCAGTCTGCGGTAAccatccat                                                                                              <  2:419738/140‑1 (MQ=255)                                             ttACCCCGCTGTCTCTGGGTATCGAAACCATGGGCGGTGTGATGACGACGCTGATCGCGAAAAACACCACTATCCCGACCAAGCACAGCCAGGCGTTCTCTACCGCTGAAAACATCCTGTCTGCGGTAAccatcactgtc                                                                                           >  1:270713/1‑134 (MQ=255)                                                ccccGCTGTCTCTGGGTATCGAAACCATGGGCGGTGTGATGACGACGCTGATCGCGAAAAACACCACTATCCCGACCAAGCACAGCCAGGCGTTCTCTACCGCTGAAGACAACCAGTCTGCGGTAACCATCCATGtgctg                                                                                        <  1:193804/140‑1 (MQ=255)                                                     cTGTCTCTGGGTATCGAAACCATGGGCGGTGTGATGACGACGCTGATCGCGAAAAACACTACTATCCCGACCAAGCACAGCCAGGCGTTCTTTACCGCTGCAGAGAACCAGTCTGCGGTAACCATCCATGTGCTGCAggg                                                                                   >  2:271012/1‑140 (MQ=255)                                                     cTGTCTCTGGGTATCGAAACCATGGGCGGTGTGATGACGACGCTGATCGCGAAAAACACCACTATCCCGACCAAGCACAGCCAGGCGTTCTCTACCGCTGAATACAACCATTCTTCGATAACCAGCCATGTGCTGCAggg                                                                                   >  1:221682/1‑140 (MQ=255)                                                        tctcTGGGTATCGAAACCATGGTCGGTGTTATGACGACGCTGATCGCGAAAAACACCACTATCCCGACCAAGCACAGCCAGGCGTTCTCTACCGCTGAAGACAACCAGTCTGCGGTAACCATCCATGTGCTGCAGGGTGa                                                                                <  1:281109/140‑1 (MQ=255)                                                        tctcTGGGAATCGAAACCATGGGCGGTGTGATGACGACGCTGATCGCGAAAAACACCACTATCCCGACCAAGCACAGCCAGGCGTTCTCTAACGCTGAAGACAACCAGCCTGCGGTAACCATGCATGTGCTGCAGGGTGa                                                                                >  2:288729/1‑140 (MQ=255)                                                           ctGGGTATCGAAACCATGGGCGGTGTGATGACGACGCTGATCGCGAAAAACACCACTATCCCGACCAAGCACAGCCAGGCGTTCTCTACCGCTGAAGACAACCAGTCTGCGGTAACCATCCATGTGCTGCAGGGTGAACg                                                                             >  2:141050/1‑140 (MQ=255)                                                             gggTATCGAAACCATGGGCGGTGTGATGACGACGCTGATCGCGAAAAACACCACTATCCCGACCAAGCACAGCCAGGCGTTCTCTACCGCTGAAGACAACCAGTCTGCGGTAACAATCCATGTGATGCAGGGTGAACGTa                                                                           >  1:275972/1‑140 (MQ=255)                                                                  tCGAAACCATGGGCGGTGTGATGACGACGCTGATCGCGAAAAACACCACTATCCCGACCAAGCACAGCCAGGCGTTCTCTACCGCTGAAGACAACCAGTCTGTGGTAACCTTCCGTGCGCTGCAGGGTGAACGTAAACGt                                                                      >  1:217846/1‑140 (MQ=255)                                                                  tCGAAACCATGGGCGGTGTGATGACGACGCTGATCGCGAAAAACACCACTATCCCGACCAAGCACAGCCAGGCGTTCTCTACCGCTGAAGACAACCAGTCTGCGGTAACCATCCATGTGCTGCAGGGTGAACGTAAACGt                                                                      >  1:758445/1‑140 (MQ=255)                                                                  tCGAAACCATGGGCGGTGTGATGACGACGCTGATCGCGAAAAACACCACTATCCCGACCAAGCACAGCCAGGCGTTCTCTACCGCTGAAGACAACCAGTCCGCGGTAACCATCCATGGGCTGCAGGGTGAACGTAAACGt                                                                      >  2:376419/1‑140 (MQ=255)                                                                          aTGGGCGGTGTGATGACGACGCTGATCGCGAAAAACACCACTATCCCGACCAAGCACAGCCAGGCGTTCTCTACCGCTGAAGACAACCAGTCTGCGGTAACCATCCATGTGCTGTAGGGTGAACGTAAACGTGCGGCTGa                                                              >  1:343914/1‑140 (MQ=255)                                                                          aTGGGCGGTGTGATGACGACGCCGATCGCTAAAAACAGCACTATCCCGACCAAGCACAGCCAGGCGTTCTCTACCGCTGAAGACAACCAGTCTGCGGTAACCATCCATGTGCTGCAGGGTGAACGTAAACGTGCGGCTGa                                                              <  1:235896/140‑1 (MQ=255)                                                                           tGGGCGGTGTGATGACGACGCTGATCGCGAAAAACACCACTATCCCGACCAAGCACAGCCAGGCGTTCTCTACCGCTGAAGACAACCAGTCAGCGGTAACTATCGATGAGCTGCAGGATGAACGTAAACGTGCGGCTGAt                                                             >  2:223101/1‑140 (MQ=255)                                                                            gggCGGTGTGATGACGACGCTGATCGCGAAAAACACCACTATCCCGACCAAGCACAGCCAGGCGTTCTCTACCGCTGAAGACAACCAGTCTGCGGTAACCGTCCGTGTACTGCAGGGGGAACTTAAACGTGCGGCTGATa                                                            >  2:236982/1‑140 (MQ=255)                                                                                ggTGTGATGACGACGCTGATCGCGAAAAACACCACTATCCCGACCAAGCACAGCCAGGCGTTCTCTACCGCTGAAGACAACCAGTCTGCGGTAACCATCCATGTGCTGCAGGGTGAACGTAAACGTg                                                                     >  2:739190/1‑127 (MQ=255)                                                                                ggTGTGATGACGACGCTGATCGCGAAAAACACCACTATCCCGACCAAGCACAGCCAGGCGTTCTCTACCGCTGAAGACAACCAGTCTGCGGTAACCATCCATGTGCTGCAGGGTGAACGTAAACGTg                                                                     <  1:739190/127‑1 (MQ=255)                                                                                ggTGTGATGACGACGCTGATCGCGAAAAACACCACTATCCCGACCAAGCACAGCCAGGCGTTCTCTACCGCTGAAGACAACCAGTCTGCGGTAACCATCCATGTGCTGCAGGGTGAACGTAAACGTGCGGCTGATAACaa                                                        <  2:685575/140‑1 (MQ=255)                                                                                  tgtgATGACGACGCTGATCGCGAAAAACACCACTATCCCGACCAAGCACAGCCAGGCGTTCTCTACCGCTGAAGACAACCAGTCTGCGGTAACCATCCATGTGCTGCAGGGTGAACGTAAACGTGCGGCTGATAACAAAt                                                      >  1:92540/1‑140 (MQ=255)                                                                                   gtgATGACGACGCTGATCGCGAAAAACACCACTATCCCGACCAAGCACAGCCAGGCGTTCTCTACCGCTGAAGACAACCAGTCTGCGGTAACCATCCATGTGCTGCAGGGTGAACGTAAACGTGCGGCTGATAACAAAtc                                                     >  1:531740/1‑140 (MQ=255)                                                                                     gatgaCGACGCTGATAGCGAAACACACCACTATCCCGTCCAAGCACAGCCAGGCGTTCTCTACCGCTGAAGACAACCAGTCTGCGGTAACCATCCATGTGCTGCAGGGTGAACGTAAACGTGCGGCTGATAACAAAtctc                                                   <  1:236982/140‑1 (MQ=255)                                                                                          cgacgCTGATCGCGAAAAACACCACTATCCCGACCAAGCACAGCCAGGCGTTCTCTACCGCTGAAGACAACCAGTCTGCGGTAACCATCCATGTGCTGCAGGGTGAACGTAAACGTGCGGCTGATAACAAATCTCTGGGt                                              <  2:569313/140‑1 (MQ=255)                                                                                               cTGATCGCGAAAAACACCACTATCCCGACCAAGCACAGCCAGGCGTTCTCTACCGCTGAAGACAACCAGTCTGCGGTAACCATCCATGTGCTGCAGGGTGAACGTAAACGTGCGGCTGATAACAAATCTCTGGGTCAGtt                                         <  2:758445/140‑1 (MQ=255)                                                                                                        gagacagccacTATCTCGACCACGCTCAGCCAGGCGTTCTCTACCGCTGAAGACAACCAGTCTGCGGTAACCATCCATGTGCTGCAGGGTGAACGTAAACGTGCGGCTGATAACAAATCTCTGGGTCAGTTCAACCTAGa                                <  1:309540/133‑1 (MQ=255)                                                                                                              accacTATCCCGACCAAGCACAGCCAGGCGTTCTCTACCGCTGAAGACAACCAGTCTGCGGTAACCATCCATGTGCTGCAGGGTGAACGTAAACGTGCGGCTGATAACAAATCTCTGGGTCAGTTCAACCTAGATGGTAt                          <  2:678464/140‑1 (MQ=255)                                                                                                              accacTATCCCGACCAAGCACAGCCAGGCGTTCTCTACCGCTGAAGACAACCAGTCTGCGGTAACCATCCATGTGCTGCAGGGTGAACGTAAACGTGCGGCTGATAACAAATCTCTGGGTCAGTTCAACCTAGATGGTAt                          <  1:141050/140‑1 (MQ=255)                                                                                                                        cGACCAAGCACAGCCAGGCGTTCTCTACCGCTGAAGACAACCAGTCTGCGGTAACCATCCATGTGCTGCAGGGTGAACGTAAACGTGCGGCTGATAACAAATCTCTGGGTCAGTTCAACCTAGATGGTATCAACCCGGCa                >  1:483065/1‑140 (MQ=255)                                                                                                                          aCCAAGCACAGCCAGGCGTTCTCTACCGCTGAAGACAACCAGTCTGCGGTAACCATCCATGTGCTGCAGGGTGAACGTAAACGTGCGGCTGATAACAAATGTCTGGGTAAGTTCAACGTAGATGGTATCAAccctgtctc              >  2:309540/1‑134 (MQ=255)                                                                                                                              aGCACAGCCAGGCGTTCTCTACCGCTGAAGACAACCAGTCTGCGGTAACCATCCATGTGCTGCAGGGTGAACGTAAACGTGCGGCTGATAAAAAATCTCTGGGTGAGTACAACCTAGGTGGTCTCAACCCGGCACcgcgc          >  2:255087/1‑140 (MQ=255)                                                                                                                               gCACAGCCAGGCGTTCTCTACCGCTGAAGACAACCAGTCTGCGGTAACCATCCATGTGCTGCAGGGTGAACGTAAACGTGCGGCTGATAACAAATCTCTGGGTCAGTTCAACCTAGATGGTATCAACCCGGCACcgcgcg         >  2:742029/1‑140 (MQ=255)                                                                                                                                 acaGCCAGGCGTTCTCTACCGCTGATGACCACCAGTCTGCGGTAACCATCCATGTGCTGCAGGGTGAACGTAAACGTGCGGCTGATAACAAATCTCTGGGTCAGTTCAACCTAGATGGTATCAACCCGGCACCGCGCGGc       <  1:288729/140‑1 (MQ=255)                                                                                                                                      cAGGCGTTCTCTACCGCTGAAGACAACCAGTCTGCGGTAACCATCCATGTGCTGCAGGGTGAACGTAAACGTGCGGCTGATAACAAATCTCTGGGTCAGTTCAACCTAGATGGTATCAACCCGGCACCGCGCGGCATGcc  <  2:531740/140‑1 (MQ=255)                                                                                                                                          |                                                                                                                                        TGGTGTTCTGACTGGTGACGTAAAAGACGTACTGCTGCTGGACGTTACCCCGCTGTCTCTGGGTATCGAAACCATGGGCGGTGTGATGACGACGCTGATCGCGAAAAACACCACTATCCCGACCAAGCACAGCCAGGTGTTCTCTACCGCTGAAGACAACCAGTCTGCGGTAACCATCCATGTGCTGCAGGGTGAACGTAAACGTGCGGCTGATAACAAATCTCTGGGTCAGTTCAACCTAGATGGTATCAACCCGGCACCGCGCGGCATGCC  >  NC_000913/13299‑13571

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