Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	NC_000913	443,051	C→T	*197* (TAG→TAA)	yajL ←	oxidative‑stress‑resistance chaperone

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	443,051	0	C	T	100.0%	135.5 / NA	43	*197* (TAG→TAA)	yajL	oxidative‑stress‑resistance chaperone
Reads supporting (aligned to +/- strand):  ref base C (0/0);  new base T (22/21);  total (22/21)

GAGCGCGGGGTGATGAGCCGTCTGCAGGCGCTCTATCTGCGCGAGCAGGGCTTTAACAATGTGAAGGTATATCGCCCGTAATTTGTGGTTTTTACGTCGCATCTGGTCAGATGCGACGTTTGCCGCATCCGACACTACTCGTAATAATTATAAATCCCTGCCGCCATCACCAGTTGTGATGCCACTTCATGGGCTTTTTCACGCCCAACCAACAGGTCGATAATTTTCAGACCAAAGTCGATAGCTGTACCCGGCCCCTGGCTGGTC  >  NC_000913/442917‑443183                                                                                                                                       |                                                                                                                                     gAGCGCGGGGTGATGAGCCGTCTGCAGGCGCTCTATCTGCGCGAGCAGGGCTTTAACAATGTGAAGGTATATCGCCCGTAATTTGTGGTTTTTACGTCGCATCTGGTCAGATGCGACGTTTGCCGCATCCGACattact                                                                                                                                  <  2:25971/139‑1 (MQ=255)       ggggTGATGAGCCGTCTGCAGGCGCTCTATCTGCGCGAGCAGGGCTTTAACAATGTGAAGGTATATCGCCCGTAATTTGTGGTTTTTACGTCGCATCTGGTCAGATGCGACGTTTGCCGCATCCGACATTACTCGtaat                                                                                                                            <  1:1041047/139‑1 (MQ=255)                               ctctATCTGCGCGAGCAGGGCTTTAACAATGTGAAGGTATATCGCCCGTAATTTGTGGTTTTTACGTCGCATCTGGTCAGATGCGACGTTTGCCGCATCCGACATTACTCGTAATAATTATAAATCCCTGCCGCCATCa                                                                                                    >  2:1575283/1‑139 (MQ=255)                               ctctATCTGCGCGAGCAGGGCTTTAACAATGTGAAGGTATATCGCCCGTAATTTGTGGTTTTTACGTCGCATCTGGTCAGATGCGACGTTTGCCGCATCCGACATTACTCGTAATAATTATAAATCCCTGCCGCCATCa                                                                                                    >  2:1606362/1‑139 (MQ=255)                                      tGCGCGAGCAGGGCTTTAACAATGTGAAGGTATATCGCCCGTAATTTGTGGTTTTTACGTCGCATCTGGTCAGATGCGACGTTTGCCGCATCCGACATTACTCGTAATAATTATAAATCCCTGCCGCCATCACCAGTtg                                                                                             >  2:70898/1‑139 (MQ=255)                                              cAGGGCTTTAACAATGTGAAGGTATATCGCCCGTAATTTGTGGTTTTTACGTCGCATCTGGTCAGATGCGACGTTTGCCGCATCCGACATTACTCGtaat                                                                                                                            <  2:129269/100‑1 (MQ=255)                                              cAGGGCTTTAACAATGTGAAGGTATATCGCCCGTAATTTGTGGTTTTTACGTCGCATCTGGTCAGATGCGACGTTTGCCGCATCCGACATTACTCGtaat                                                                                                                            >  1:129269/1‑100 (MQ=255)                                                           aTGTGAAGGTATATCGCCCGTAATTTGTGGTTTTTACGTCGCATCTGGTCAGATGCGACGTTTGCCGCATCCGACATTACTCGtaat                                                                                                                            >  2:1173766/1‑87 (MQ=255)                                                           aTGTGAAGGTATATCGCCCGTAATTTGTGGTTTTTACGTCGCATCTGGTCAGATGCGACGTTTGCCGCATCCGACATTACTCGtaat                                                                                                                            <  1:1173766/87‑1 (MQ=255)                                                           aTGTGAAGGTATATCGCCCGTAATTTGTGGTTTTTACGTCGCATCTGGTCAGATGCGACGTTTGCCGCATCCGACATTACTCGTAATAATTATAAATCCCTGCCGCCATCACCAGTtgt                                                                                            <  1:816618/119‑1 (MQ=255)                                                           aTGTGAAGGTATATCGCCCGTAATTTGTGGTTTTTACGTCGCATCTGGTCAGATGCGACGTTTGCCGCATCCGACATTACTCGTAATAATTATAAATCCCTGCCGCCATCACCAGTtgt                                                                                            >  2:816618/1‑119 (MQ=255)                                                             gtgAAGGTATATCGCCCGTAATTTGTGGTTTTTACGTCGCATCTGGTCAGATGCGACGTTTGCCGCATCCGACATTACTCGTAATAATTATAAATCCCTGCCGCCATCACCAGTTGTGATGCCACTTCATGGGCttttt                                                                      >  1:967092/1‑139 (MQ=255)                                                             gtgAAGGTATATCGCCCGTAATTTGTGGTTTTTACGTCGCATCTGGTCAGATGCGACGTTTGCCGCATCCGACATTACTCGTAATAATTATAAATCCCTGCCGCCATCACCAGTTGTGATGCCACTTCATGGGCttttt                                                                      <  2:967092/139‑1 (MQ=255)                                                             gtgAAGGTATATCGCCCGTAATTTGTGGTTTTTACGTCGCATCTGGTCAGATGCGACGTTTGCCGCATCCGACATTACTCGTAATAATTATAAATCCCTGCCGCCATCACCAGTTGTGATGCCACTTCATGGGCttttt                                                                      <  1:70898/139‑1 (MQ=255)                                                                       atCGCCCGTAATTTGTGGTTTTTACGTCGCATCTGGTCAGATGCGACGTTTGCCGCATCCGACATTACTCGTAATAATTATAAATCCCTGCCGCCATCACCAGTTGTGATGCCACTTCATGGGCTTTTTCAc                                                                   <  2:113374/132‑1 (MQ=255)                                                                       atCGCCCGTAATTTGTGGTTTTTACGTCGCATCTGGTCAGATGCGACGTTTGCCGCATCCGACATTACTCGTAATAATTATAAATCCCTGCCGCCATCACCAGTTGTGATGCCACTTCATGGGCTTTTTCAc                                                                   >  1:113374/1‑132 (MQ=255)                                                                       atCGCCCGTAATTTGTGGTTTTTACGTCGCATCTGGTCAGATGCGACGTTTGCCGCATCCGACATTACTCGTAATAATTATAAATCCCTGCCGCCATCACCAGTTGTGATGCCACTTCATGGGCTTTTTCACGc                                                                 <  2:1570433/134‑1 (MQ=255)                                                                       atCGCCCGTAATTTGTGGTTTTTACGTCGCATCTGGTCAGATGCGACGTTTGCCGCATCCGACATTACTCGTAATAATTATAAATCCCTGCCGCCATCACCAGTTGTGATGCCACTTCATGGGCTTTTTCACGc                                                                 >  1:1570433/1‑134 (MQ=255)                                                                        tCGCCCGTAATTTGTGGTTTTTACGTCGCATCTGGTCAGATGCGACGTTTGCCGCATCCGACATTACTCGTAATAATTATAAATCCCTGCCGCCATCAc                                                                                                   >  2:215662/1‑99 (MQ=255)                                                                        tCGCCCGTAATTTGTGGTTTTTACGTCGCATCTGGTCAGATGCGACGTTTGCCGCATCCGACATTACTCGTAATAATTATAAATCCCTGCCGCCATCAc                                                                                                   <  1:215662/99‑1 (MQ=255)                                                                        tCGCCCGTAATTTGTGGTTTTTACGTCGCATCTGGTCAGATGCGACGTTTGCCGCATCCGACATTACTCGTAATAATTATAAATCCCTGCCGCCATCACCAGTTGTGATGCCACTTCATGGGCTTTTTCACGCccaacc                                                           >  1:779205/1‑139 (MQ=255)                                                                          gCCCGTAATTTGTGGTTTTTACGTCGCATCTGGTCAGATGCGACGTTTGCCGCATCCGACattact                                                                                                                                  <  1:1119418/66‑1 (MQ=255)                                                                          gCCCGTAATTTGTGGTTTTTACGTCGCATCTGGTCAGATGCGACGTTTGCCGCATCCGACattact                                                                                                                                  >  2:1119418/1‑66 (MQ=255)                                                                          gCCCGTAATTTGTGGTTTTTACGTCGCATCTGGTCAGATGCGACGTTTGCCGCATCCGACATTACTCGTAATAATTATAAATCCCTGCCGCCATCACCAGTTGTGATGCCACTTCATGGGCTTTTTCACGCccaaccaa                                                         <  1:1833229/139‑1 (MQ=255)                                                                           cccGTAATTTGTGGTTTTTACGTCGCATCTGGTCAGATGCGACGTTTGCCGCATCCGACATTACTCGtaata                                                                                                                           >  1:1130121/1‑72 (MQ=255)                                                                           cccGTAATTTGTGGTTTTTACGTCGCATCTGGTCAGATGCGACGTTTGCCGCATCCGACATTACTCGtaata                                                                                                                           <  2:1130121/72‑1 (MQ=255)                                                                             cGTAATTTGTGGTTTTTACGTCGCATCTGGTCAGATGCGACGTTTGCCGCATCCGACattact                                                                                                                                  >  1:222895/1‑63 (MQ=255)                                                                             cGTAATTTGTGGTTTTTACGTCGCATCTGGTCAGATGCGACGTTTGCCGCATCCGACattact                                                                                                                                  <  2:222895/63‑1 (MQ=255)                                                                             cGTAATTTGTGGTTTTTACGTCGCATCTGGTCAGATGCGACGTTTGCCGCATCCGACATTACTCGtaataa                                                                                                                          >  2:982927/1‑71 (MQ=255)                                                                             cGTAATTTGTGGTTTTTACGTCGCATCTGGTCAGATGCGACGTTTGCCGCATCCGACATTACTCGtaataa                                                                                                                          <  1:982927/71‑1 (MQ=255)                                                                              gTAATTTGTGGTTTTTACGTCGCATCTGGTCAGATGCGACGTTTGCCGCATCCGACATTACTCGTAATAATTATAAATCCCTgccgcc                                                                                                        >  1:247020/1‑88 (MQ=255)                                                                              gTAATTTGTGGTTTTTACGTCGCATCTGGTCAGATGCGACGTTTGCCGCATCCGACATTACTCGTAATAATTATAAATCCCTgccgcc                                                                                                        <  2:247020/88‑1 (MQ=255)                                                                              gTAATTTGTGGTTTTTACGTCGCATCTGGTCAGATGCGACGTTTGCCGCATCCGACATTACTCGTAATAATTATAAATCCCTGCCGCCATCACCAGTTGTGATGcc                                                                                      >  1:458355/1‑106 (MQ=255)                                                                              gTAATTTGTGGTTTTTACGTCGCATCTGGTCAGATGCGACGTTTGCCGCATCCGACATTACTCGTAATAATTATAAATCCCTGCCGCCATCACCAGTTGTGATGcc                                                                                      <  2:458355/106‑1 (MQ=255)                                                                               tAATTTGTGGTTTTTACGTCGCATCTGGTCAGATGCGACGTTTGCCGCATCCGACATTACTCGTAATAATTATAAATCCCTGCCGCCATCACCAGTTGTGATGCCACTTCATGGGCTTTTTCACGCCCAACCAACAGGt                                                    <  1:1575283/139‑1 (MQ=255)                                                                                  tttGTGGTTTTTACGTCGCATCTGGTCAGATGCGACGTTTGCCGCATCCGACATTACTCGTAATAATTATAAATCCCTg                                                                                                             <  2:212563/79‑1 (MQ=255)                                                                                  tttGTGGTTTTTACGTCGCATCTGGTCAGATGCGACGTTTGCCGCATCCGACATTACTCGTAATAATTATAAATCCCTg                                                                                                             >  1:212563/1‑79 (MQ=255)                                                                                      tgGTTTTTACGTCGCATCTGGTCAGATGCGACGTTTGCCGCATCCGACATTACTCGTAATAATTATAAATcc                                                                                                                >  2:1708951/1‑72 (MQ=255)                                                                                      tgGTTTTTACGTCGCATCTGGTCAGATGCGACGTTTGCCGCATCCGACATTACTCGTAATAATTATAAATcc                                                                                                                <  1:1708951/72‑1 (MQ=255)                                                                                      tgGTTTTTACGTCGCATCTGGTCAGATGCGACGTTTGCCGCATCCGACATTACTCGTAATAATTATAAATCCCTGCCGCCATc                                                                                                     <  2:1421330/83‑1 (MQ=255)                                                                                      tgGTTTTTACGTCGCATCTGGTCAGATGCGACGTTTGCCGCATCCGACATTACTCGTAATAATTATAAATCCCTGCCGCCATc                                                                                                     >  1:1421330/1‑83 (MQ=255)                                                                                                   gCATCTGGTCAGATGCGACGTTTGCCGCATCCGACATTACTCGTAATAATTATAAATCCCTGCCGCCATCACCAGTTGTGATGCCACTTCATGGGCTTTTTCACGCCCAACCAACAGGTCGATAAtttt                                          <  1:755402/129‑1 (MQ=255)                                                                                                   gCATCTGGTCAGATGCGACGTTTGCCGCATCCGACATTACTCGTAATAATTATAAATCCCTGCCGCCATCACCAGTTGTGATGCCACTTCATGGGCTTTTTCACGCCCAACCAACAGGTCGATAAtttt                                          >  2:755402/1‑129 (MQ=255)                                                                                                        tGGTCAGATGCGACGTGTGCCGCATCCGACATTACTCGTAATAATTATAAATCCCTGCCGCCATCACCAGTTGTGATGCCACTTCATGGGCTTTTTCACGCCCAACCAACAGGTCGATAATTTTCAGACCAAAGTCGAt                           <  1:1034138/139‑1 (MQ=255)                                                                                                               aTGCGACGTTTGCCGCATCCGACATTACTCGTAATAATTATAAATCCCTgccg                                                                                                          >  2:40964/1‑53 (MQ=255)                                                                                                               aTGCGACGTTTGCCGCATCCGACATTACTCGTAATAATTATAAATCCCTgccg                                                                                                          <  1:40964/53‑1 (MQ=255)                                                                                                                         tGCCGCATCCGACATTACTCGTAATAATTATAAATCCCTGCCGCCATCACCAGTTGTGATGCCACTTCATGGGCTTTTTCACGCCCAACCAACAGGTCGATAATTTTCAGACCAAAGTCGATAGCTGTACCCGGCCCct          >  1:968461/1‑139 (MQ=255)                                                                                                                                 ccGACATTACTCGTAATAATTATAAATCCCTGCCGCCATCACCAGTTGTGATGCCACTTCATGGGCTTTTTCACGCCAAACCAACAGGCCGAAAATTTTCAGACCAACTTCTATACCTATCCCCCGCCCctgccgggcc  >  1:185670/1‑132 (MQ=255)                                                                                                                                       |                                                                                                                                     GAGCGCGGGGTGATGAGCCGTCTGCAGGCGCTCTATCTGCGCGAGCAGGGCTTTAACAATGTGAAGGTATATCGCCCGTAATTTGTGGTTTTTACGTCGCATCTGGTCAGATGCGACGTTTGCCGCATCCGACACTACTCGTAATAATTATAAATCCCTGCCGCCATCACCAGTTGTGATGCCACTTCATGGGCTTTTTCACGCCCAACCAACAGGTCGATAATTTTCAGACCAAAGTCGATAGCTGTACCCGGCCCCTGGCTGGTC  >  NC_000913/442917‑443183

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