Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	NC_000913	12,086	C→T	intergenic (‑300/‑77)	yaaI ← / → dnaK	DUF2541 domain‑containing protein YaaI/chaperone protein DnaK

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	12,086	0	C	T	100.0%	220.9 / NA	56	intergenic (‑300/‑77)	yaaI/dnaK	DUF2541 domain‑containing protein YaaI/chaperone protein DnaK
Reads supporting (aligned to +/- strand):  ref base C (0/0);  new base T (33/23);  total (33/23)

GTAGTCAACCGCAGTGAGTGAGTCTGCAAAAAAATGAAATTGGGCAGTTGAAACCAGACGTTTCGCCCCTATTACAGACTCACAACCACAT  >  NC_000913/12042‑12132                                             |                                               gtagtCAACCGCAGTGAGTGAGTCTGCAAAAAAATGAAATTGGGTAGtt                                            <  1:2581172/49‑1 (MQ=255) gtagtCAACCGCAGTGAGTGAGTCTGCAAAAAAATGAAATTGGGTAGTTg                                           <  1:1870553/50‑1 (MQ=255)  tagtCAACCGCAGTGAGTGAGTCTGCAAAAAAATGAAATTGGGTAGTTGa                                          <  1:4374100/50‑1 (MQ=255)    gtCAACCGCAGTGAGTGAGTCTGCAAAAAAATGAAATTGGGTAGTTGaaa                                        >  1:937933/1‑50 (MQ=255)     tCAACCGCAGTGAGTGAGTCTGCAAAAAAATGAAATTGGGTAGTTGAAAc                                       <  1:217733/50‑1 (MQ=255)     tCAACCGCAGTGAGTGAGTCTGCAAAAAAATGAAATTGGGTAGTTGAAAc                                       >  1:4220003/1‑50 (MQ=255)       aaCCGCAGTGAGTGAGTCTGCAAAAAAATGAAATTGGGTAGTTGAAAcc                                      <  1:2449145/49‑1 (MQ=255)       aaCCGCAGTGAGTGAGTCTGCAAAAAAATGAAATTGGGTAGTTGAAAcc                                      <  1:1390904/49‑1 (MQ=255)       aaCCGCAGTGAGTGAGTCTGCAAAAAAATGAAATTGGGTAGTTGAAAc                                       >  1:3992154/1‑48 (MQ=255)        aCCGCAGTGAGTGAGTCTGCAAAAAAATGAAATTGGGTAGTTGAAACCAg                                    <  1:3065787/50‑1 (MQ=255)        aCCGCAGTGAGTGAGTCTGCAAAAAAATGAAATTGGGTAGTTGAAACCAg                                    <  1:70561/50‑1 (MQ=255)          cGCAGTGAGTGAGTCTGCAAAAAAATGAAATTGGGTAGTTGAAACCAGAc                                  >  1:812377/1‑50 (MQ=255)           gCAGTGAGTGAGTCTGCAAAAAAATGAAATTGGGTAGTTGAAACCAGACg                                 >  1:326362/1‑50 (MQ=255)            cAGTGAGTGAGTCTGCAAAAAAATGAAATTGGGTAGTTGAAACCAGACGt                                >  1:1703736/1‑50 (MQ=255)              gtgagtgaGTCTGCAAAAAAATGAAATTGGGTAGTTGAAACCAGACGttt                              >  1:3752862/1‑50 (MQ=255)                gagtgagtCTGCAAAAAAATGAAATTGGGTAGTTGAAACCAGACGTTTCg                            >  1:688602/1‑50 (MQ=255)                  gtgagtCTGCAAAAAAATGAAATTGGGTAGTTGAAACCAGACGTTTCGcc                          >  1:703481/1‑50 (MQ=255)                   tgagtCTGCAAAAAAATGAAATTGGGTAGTTGAAACCAGACGTTTCGccc                         >  1:509108/1‑50 (MQ=255)                   tgagtCTGCAAAAAAATGAAATTGGGTAGTTGAAACCAGACGTTTCGccc                         >  1:2947766/1‑50 (MQ=255)                    gagtCTGCAAAAAAATGAAATTGGGTAGTTGAAACCAGACGTTTCGcccc                        >  1:477255/1‑50 (MQ=255)                     agtCTGCAAAAAAATGAAATTGGGTAGTTGAAACCAGACGTTTCGcccc                        <  1:3740700/49‑1 (MQ=255)                     agtCTGCAAAAAAATGAAATTGGGTAGTTGAAACCAGACGTTTCGCCCCt                       <  1:95578/50‑1 (MQ=255)                      gtCTGCAAAAAAATGAAATTGGGTAGTTGAAACCAGACGTTTCGCCCCTa                      >  1:1026343/1‑50 (MQ=255)                         tGCAAAAAAATGAAATTGGGTAGTTGAAACCAGACGTTTCGCCCCTATTa                   >  1:2159572/1‑50 (MQ=255)                          gCAAAAAAATGAAATTGGGTAGTTGAAACCAGACGTTTCGCCCCTATTa                   >  1:2620579/1‑49 (MQ=255)                          gCAAAAAAATGAAATTGGGTAGTTGAAACCAGACGTTTCGCCCCTATTa                   >  1:1494151/1‑49 (MQ=255)                           cAAAAAAATGAAATTGGGTAGTTGAAACCAGACGTTTCGCCCCTAtt                    <  1:2611092/47‑1 (MQ=255)                           cAAAAAAATGAAATTGGGTAGTTGAAACCAGACGTTTCGCCCCTATTa                   >  1:4178211/1‑48 (MQ=255)                           cAAAAAAATGAAATTGGGTAGTTGAAACCAGACGTTTCGCCCCTATTa                   <  1:3811427/48‑1 (MQ=255)                           cAAAAAAATGAAATTGGGTAGTTGAAACCAGACGTTTCGCCCCTATTa                   <  1:2166696/48‑1 (MQ=255)                           cAAAAAAATGAAATTGGGTAGTTGAAACCAGACGTTTCGCCCCTATTAc                  <  1:2067864/49‑1 (MQ=255)                           cAAAAAAATGAAATTGGGTAGTTGAAACCAGACGTTTCGCCCCTATTACa                 <  1:4322136/50‑1 (MQ=255)                            aaaaaaaTGAAATTGGGTAGTTGAAACCAGACGTTTCGCCCCTATTa                   >  1:630359/1‑47 (MQ=255)                             aaaaaaTGAAATTGGGTAGTTGAAACCAGACGTTTCGCCCCTATTACAGa               >  1:1140151/1‑50 (MQ=255)                              aaaaaTGAAATTGGGTAGTTGAAACCAGACGTTTCGCCCCTATTACAGAc              >  1:3470247/1‑50 (MQ=255)                               aaaaTGAAATTGGGTAGTTGAAACCAGACGTTTCGCCCCTATTACAGACt             >  1:2472104/1‑50 (MQ=255)                                aaaTGAAATTGGGTAGTTGAAACCAGACGTTTCGCCCCTATTACAGACt             >  1:3383341/1‑49 (MQ=255)                                aaaTGAAATTGGGTAGTTGAAACCAGACGTTTCGCCCCTATTACAGACt             >  1:334061/1‑49 (MQ=255)                                aaaTGAAATTGGGTAGTTGAAACCAGACGTTTCGCCCCTATTACAGACt             >  1:2466387/1‑49 (MQ=255)                                aaaTGAAATTGGGTAGTTGAAACCAGACGTTTCGCCCCTATTACAGACt             >  1:2047356/1‑49 (MQ=255)                                 aaTGAAATTGGGTAGTTGAAACCAGACGTTTCGCCCCTATTACAGACTca           <  1:1394602/50‑1 (MQ=255)                                 aaTGAAATTGGGTAGTTGAAACCAGACGTTTCGCCCCTATTACAGACTca           <  1:2712073/50‑1 (MQ=255)                                 aaTGAAATTGGGTAGTTGAAACCAGACGTTTCGCCCCTATTACAGACTca           >  1:2802515/1‑50 (MQ=255)                                  aTGAAATTGGGTAGTTGAAACCAGACGTTTCGCCCCTATTACAGACTcac          <  1:591818/50‑1 (MQ=255)                                  aTGAAATTGGGTAGTTGAAACCAGACGTTTCGCCCCTATTACAGACTca           >  1:4050207/1‑49 (MQ=255)                                   tGAAATTGGGTAGTTGAAACCAGACGTTTCGCCCCTATTACAGACTcac          >  1:2241609/1‑49 (MQ=255)                                     aaaTTGGGTAGTTGAAACCAGACGTTTCGCCCCTATTACAGACTCACaa        <  1:435028/49‑1 (MQ=255)                                     aaaTTGGGTAGTTGAAACCAGACGTTTCGCCCCTATTACAGACTCACaa        <  1:3148015/49‑1 (MQ=255)                                     aaaTTGGGTAGTTGAAACCAGACGTTTCGCCCCTATTACAGACTCACAAc       <  1:918295/50‑1 (MQ=255)                                        ttGGGTAGTTGAAACCAGACGTTTCGCCCCTATTACAGACTCACAACca     >  1:1052710/1‑49 (MQ=255)                                        ttGGGTAGTTGAAACCAGACGTTTCGCCCCTATTACAGACTCACAACca     >  1:383245/1‑49 (MQ=255)                                        ttGGGTAGTTGAAACCAGACGTTTCGCCCCTATTACAGACTCACAACca     >  1:1408738/1‑49 (MQ=255)                                         tGGGTAGTTGAAACCAGACGTTTCGCCCCTATTACAGACTCACAACcaca   >  1:1778157/1‑50 (MQ=255)                                          gggTAGTTGAAACCAGACGTTTCGCCCCTATTACAGACTCACAACcaca   <  1:2651511/49‑1 (MQ=255)                                          gggTAGTTGAAACCAGACGTTTCGCCCCTATTACAGACTCACAACcaca   <  1:1120297/49‑1 (MQ=255)                                          gggTAGTTGAAACCAGACGTTTCGCCCCTATTACAGACTCACAACCACat  >  1:987229/1‑50 (MQ=255)                                             |                                               GTAGTCAACCGCAGTGAGTGAGTCTGCAAAAAAATGAAATTGGGCAGTTGAAACCAGACGTTTCGCCCCTATTACAGACTCACAACCACAT  >  NC_000913/12042‑12132

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