Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	NC_000913	2,041,467	A→G	K31K (AAA→AAG)	zinT →	metal‑binding protein ZinT

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	2,041,467	0	A	G	100.0%	122.5 / NA	34	K31K (AAA→AAG)	zinT	metal‑binding protein ZinT
Reads supporting (aligned to +/- strand):  ref base A (0/0);  new base G (17/17);  total (17/17)

CCTCAAATCACCAAAATAATTTATTTTATTTCAATAAGATAGATAAGTAGAACTGAGAAAGCCATGCTCTCGTTTCCTAAGAGTTGTTGCATTTTGCTATATGTTACAATATAACATTACACATCATATACATTAACTCTGGAGGAAACTGTTTTGGCGATTCGTCTTTACAAACTGGCTGTTGCTTTAGGTGTCTTTATTGTTAGCGCTCCTGCCTTTTCGCATGGTCATCACTCACACGGCAAACCCTTAACAGAGGTCGAACAAAAAGCTGCTAATGGTGTTTTTGATGATGCCAATGTACAAAACCGAACGCTCAGTGACTGGGATGGAGTCTGGCAATCCGTTTATCCTTTACTGCAAAGTGGCAAACTTGACCCCGTCTTTCAGAAGAAAGCGGATGCAGATAAAACTAAAACATTTGCTGAAATTAAAGATTATTATCACAAAGGTTATGCAACAGATATCGAGATGATTGGCATTGAGGAC  >  NC_000913/2041222‑2041710                                                                                                                                                                                                                                                      |                                                                                                                                                                                                                                                    cgcggtcttgctgggtacgggactgggtggagaaacactgtcgaagattgaatggctggcgctcggcgtaattgtcttcgcggtggtactggtcacgttgggaaaatatctcttcccggcaaaacccgtagttgcgccagttattcaggacgcatcaCGATTCGTCTTTACAAACTGGCTGTTGCTTTAGGTGTCTTTATTGTTAGCGCTCCTGCCTTTTCGCATGGTCATCACTCACACGGCAAGCCCtt                                                                                                                                                                                                                                                <  1:163444‑M2/94‑1 (MQ=255)                     gactgggtggagaaacactgtcgaagattgaatggctggcgctcggcgtaattgtcttcgcggtggtactggtcacgttgggaaaatatctcttcccggcaaaacccgtagttgcgccagttattcaggacgcatcaCGATTCGTCTTTACAAACTGGCTGTTGCTTTAGGTGTCTTTATTGTTAGCGCTCCTGCCTTTTCGCATGGTCATCACTCACACGGCAAGCCCTTAACAGAGGTCGAACAAAAAg                                                                                                                                                                                                                            <  1:313222‑M2/114‑1 (MQ=255)                                 aaacactgtcgaagattgaatggctggcgctcggcgtaattgtcttcgcggtggtactggtcacgttgggaaaatatctcttcccggcaaaacccgtagttgcgccagttattcaggacgcatcaCGATTCGTCTTTACAAACTGGCTGTTGCTTTAGGTGTCTTTATTGTTAGCGCTCCTGCCTTTTCGCATGGTCATCACTCACACGGCAAGCCCTTAACAGAGGTCGAACAAAAAGCTGCTAATGgtg                                                                                                                                                                                                                >  1:42445‑M2/126‑251 (MQ=255)                                   acactgtcgaagattgaatggctggcgctcggcgtaattgtcttcgcggtggtactggtcacgttgggaaaatatctcttcccggcaaaacccgtagttgcgccagttattcaggacgcatcaCGATTCGTCTTTACAAACTGGCTGTTGCTTTAGGTGTCTTTATTGTTAGCGCTCCTGCCTTTTCGCATGGTCATCACTCACACGGCAAGCCCTTAACAGAGGTCGAACAAAAAGCTGCTAATGGTGtt                                                                                                                                                                                                              >  1:39084‑M2/124‑251 (MQ=255)                                                                cggcgtaattgtcttcgcggtggtactggtcacgttgggaaaatatctcttcccggcaaaacccgtagttgcgccagttattcaggacgcatcaCGATTCGTCTTTACAAACTGGCTGTTGCTTTAGGTGTCTTTATTGTTAGCGCTCCTGCCTTTTCGCATGGTCATCACTCACACGGCAAGCCCTTAACAGAGGTCGAACAAAAAGCTGCTAATGGTGTTTTTGATGATGCCAATGTACAAAACCGAAc                                                                                                                                                                                 <  2:365089‑M2/157‑1 (MQ=255)                                                                        ttgtcttcgcggtggtactggtcacgttgggaaaatatctcttcccggcaaaacccgtagttgcgccagttattcaggacgcatcaCGATTCGTCTTTACAAACTGGCTGTTGCTTTAGGTGTCTTTATTGTTAGCGCTCCTGCCTTTTCGCATGGTCATCACTCACACGGCAAGCCCTTAACAGAGGTCGAACAAAAAGCTGCTAATGGTGTTTTTGATGATGCCAATGTACAAAACCGAACGCTCAGTg                                                                                                                                                                         >  1:65533‑M2/87‑251 (MQ=255)                                                                                        actggtcacgttgggaaaatatctcttcccggcaaaacccgtagttgcgccagttattcaggacgcatcaCGATTCGTCTTTACAAACTGGCTGTTGCTTTAGGTGTCTTTATTGTTAGCGCTCCTGCCTTTTCGCATGGTCATCACTCACACGGCAAGCCCTTAACAGAGGTCGAACAAAAAGCTGCTAATGGTGTTTTTGATGATGCCAATGTACAAAACCGAACGCTCAGTGACTGGGATGGAGTCTg                                                                                                                                                         >  2:339118‑M2/71‑251 (MQ=255)                                                                                        actggtcacgttgggaaaatatctcttcccggcaaaacccgtagttgcgccagttattcaggacgcatcaCGATTCGTCTTTACAAACTGGCTGTTGCTTTAGGAGTCTTTATTGTTAGCGCTCCTGCCTTTTCGCATGGTCATCACTCACACGGCAAGCCCTTAACAGAGGTCGAACAAAAAGCTGCTCATGGTGTTTTTGATGATGCCAATGTCCAAAACCGAACGCTCAGTGACTGGGATGGAGTCTg                                                                                                                                                         >  2:142925‑M2/71‑251 (MQ=255)                                                                                         ctggtcacgttgggaaaatatctcttcccggcaaaacccgtagttgcgccagttattcaggacgcatcaCGATTCGTCTTTACAAACTGGCTGTTGCTTTAGGTGTCTTTATTGTTAGCGCTCCTGCCTTTTCGCATGGTCATCACTCACACGGCAAGCCCTTAACAGAGGTCGAACAAAAAGCTGCTAATGGTGTTTTTGATGATGCCAATGTACAAAACCGAACGCTCAGTGACTGGGATGGAGTCTgg                                                                                                                                                        >  1:371731‑M2/70‑251 (MQ=255)                                                                                                                      ggcaaaacccgtagttgcgccagttattcaggacgcatcaCGATTCGTCTTTACAAACTGGCTGTTGCTTTAGGTGTCTTTATTGTTAGCGCTCCTGCCTTTTCGCATGGTCATCACTCACACGGCAAGCCCTTAACAGAGGTCGAACAAAAAGCTGCTAATGGTGTTTTTGATGATGCCAATGTACAAAACCGAACGCTCAGTGACTGGGATGGAGTCTGGCAATCCGTTTATCCTTTACTGCAAAGTgg                                                                                                                           <  1:288925‑M2/211‑1 (MQ=255)                                                                                                                               cgtagttgcgccagttattcaggacgcatcaCGATTCGTCTTTACAAACTGGCTGTTGCTTTAGGTGTCTTTATTGTTAGCGCTCCTGCCTTTTCGCATGGTCATCACTCACACGGCAAGCCCTTAACAGAGGTCGAACAAAAAGCTGCTAATGGTGTTTTTGATGATGCCAATGTACAAAACCGAACGCTCAGTGACTGGGATGGAGTCTGGCAATCCGTTTATCCTTTACTGCAAAGTGGCAAACTTGa                                                                                                                  <  2:25114‑M2/220‑1 (MQ=255)                                                                                                                                  agttgcgccagttattcaggacgcatcaCGATTCGTCTTTACAAACTGGCTGTTGCTTTAGGTGTCTTTATTGTTAGCGCTCCTGCCTTTTCGCATGGTCATCACTCACACGGCAAGCCCTTAACAGAGGTCGAACAAAAAGCTGCTAATGGTGTTTTTGCTGATGCCAATGTACAAAACCGAACGCTCAGTGACTGGGATGGAGTCTGGCAATCCGTTTATCCTTTACTGCAAAGTGGCACACttgaacc                                                                                                               >  2:88308‑M2/29‑248 (MQ=255)                                                                                                                                   gttgcgccagttattcaggacgcatcaCGATTCGTCTTTACAAACTGGCTGTCGCTTTAGGTGTCTTTATTGTTAGCGCTCCTGCCTTTTCGCATGGTCATCACTCACACGGCAAGCCCTTAACAGAGGTCGAACAAAAAGCTGCTAATGGTGTTTTTGATGATGCCAATGTACAAAACCGAACGCTCAGTGACTGGGATGGAGTCTGGCAATCCGTTTATCCTTTACTGCAAAGTGGCAAACTTGAcccc                                                                                                              >  1:286131‑M2/28‑251 (MQ=255)                                                                                                                                       cgccagttattcaggacgcatcaCGATTCGTCTTTACAAACTGGCTGTTGCTTTAGGTGTCTTTATTGTTAGCGCTCCTGCCTTTTCGCATGGTCATCACTCACACGGCAAGCCCTTAACAGAGGTCGAACAAAAAGCTGCTAATGGTGTTTTTGATGATGCCAATGTACAAAACCGAACGCTCAGTGACTGGGATGGAGTCTGGCAATCCGTTTATCCTTTACTGCAAAGTGGCAAACGTGACCCCGTCt                                                                                                          <  1:88308‑M2/228‑1 (MQ=255)                                                                                                                                                  caggacgcatcaCGATTCGTCTTTACAAACTAGCTGTTGCTTTAGGTGTCTTTATTGTTAACGCTCCTGCCTTTTCGCATGGTCCTCACTCACACGGCAAGCCCTTAACAGAGGTCGAACAAAAAGCTGCTACTGGTGTTTTTGATGATGCCAATGTACAAAACCGAACGCTCAGTGACTGGGATGGAGTCTGGCAATCCGTTTATCCTTTACTGCAAAGTGGCAAACTTGACCCCGTCTTTCAgaagaa                                                                                                >  1:381976‑M2/13‑250 (MQ=255)                                                                                                                                                          atcaCGATTCGTCTTTACAAACTGGCTGTTGCTTTAGGTGTCTTTATTGTTAGCGCTCCTGCCTTTTCGCATGGTCATCACTCACACGGCAAGCCCTTAACAGAGGTCGAACAAAAAGCTGCTAATGGTGTTTTTGATGATGCCAATGTACAAAACCGAACGCTCAGTGACTGGGATGGAGTCTGGCAATCCGTTTATCCTTTACTGCAAAGTGGCAAACTTGACCCCGTCTTTCAGAAGAAAGCGGATGc                                                                                       >  2:186729‑M2/5‑251 (MQ=255)                                                                                                                                                            caCGATTCGTCTTTACAAACTGGCTGTTGCTTTAGGTGTCTTTATTGTTAGCGCTCCTGCCTTTTCGCATGGTCATCACTCACACGGCAAGCCCTTAACAGAGGTCGAACAAAAAGCTGCTAATGGTGTTTTTGATGATGCCAATGTACAAAACCGAACGCTCAGTGACTGGGATGGAGTCTGGCAATCCGTTTATCCTTTACTGCAAAGTGGCAAACTTGACCCCGTCTTTCAGAAGAAAGCGGATGCAg                                                                                     <  1:342052‑M2/249‑1 (MQ=255)                                                                                                                                                              cGATTCGTCTTTACAAACTGGCTGTTGCTTTAGGTGTCTTTATTGTTAGCGCTCCTGCCTTTTCGCATGGTCATCACTCACACGGCAAGCCCTTAACAGAGGTCGAACAAAAAGCTGCTAATGGTGTTTTTGATGATGCCAATGTACAAAACCGAACGCTCAGTGACTGGGATGGAGTCTGGCAATCCGTTTATCCTTTACTGCAAAGTGGCAAGCTTGACCCCGTCTTTCAGAAGAAAGCGGATGCAGAt                                                                                   <  1:372220/251‑1 (MQ=255)                                                                                                                                                              cGATGCGTCTTTACAACCTGGCTGTTGCTTTAGGTGTCTGTAGTGTTAGCGCTCCTGCCTTTTCGCAGGGTCATCACTCACACGGCAAGCCCTTAACAGAGGGCGAACAAAAAGCTGCTAATGGTGTTTTTGATGATGCCAATGTACAAAACCGAACGCTCAGTGACTGGGATGGAGTCTGGCAATCCGTTTATCCTTTACTTCAAAGTGGCAAACGTGACCCCGTCGTTCAGAAGAAAGCGTAGGCAGAt                                                                                   <  2:336060/251‑1 (MQ=255)                                                                                                                                                                         tACAAACTGGCTGTTGCTTTAGGTGTCTTTATTGTTAGCGCTCCTGCCTTTTCGCATGGTCATCACTCACACGGCAAGCCCTTAACAGAGGTCGAACAAAAAGCTGCTAATGGTGTTTTTGATGATGCCAATGTACAAAACCGAACGCTCAGTGACTGGGATGGAGTCTGGCAATCCGTTTATCCTTTACTGCAAAGTGGCAAACTTGACCCCGTCTTTCAGAAGAAAGCGGATGCAGAtaaaactaaaa                                                                         <  1:279997/250‑1 (MQ=255)                                                                                                                                                                           cAAACTGGCTGTTGCTTTAGGTGTCTTTATTGTTAGCGCTCCTGCCTTTTCGCATGGTCATCACTCACACGGCAAGCCCTTAACAGAGGTCGAACAAAAAGCTGCTAATGGTGTTTTTGATGATGCCAATGTACAAAACCGAACGCTCAGTGACTGGGATGGAGTCTGGCAATCCGTTTATCCTTTACTGCAAAGTGGCAAACTTGACCCCGTCTTTCAGAAGAAAGCGGATGCAGATAAAACTAAAACAt                                                                      >  1:91658/1‑251 (MQ=255)                                                                                                                                                                              aCTGGCTGTTGCTTTAGGTGTCTTTATTGTTAGCGCTCCTGCCTTTTCGCATGGTCATCACTCACACGGCAAGCCCTTAACAGAGGTCGAACAAAAAGCTGCTAATGGTGTTTTTGATGATGCCAATGTACAAAACCGAACGCTCAGTGACTGGGATGGAGTCTGGCAATCCGTTTATCCTTTACTGCAAAGTGGCAAACTTGACCCCGTCTTTCAGAAGAAAGCGGATGCAGATAAAACTAAAACATTTg                                                                   <  1:287384/251‑1 (MQ=255)                                                                                                                                                                                    tGTTGCTTTAGGTGTCTTTATTGTTAGCGCTCCTGCCTTTTCGCATGGTCATCACTCACACGGCAAGCCCTTAACAGAGGTCGAACAAAAAGCTGCTAATGGTGTTTTTGATGATGCCAATGTACAAAACCGAACGCTCAGTGACTGGGATGGAGTCTGGCAATCCGTTTATCCTTTACTGCAAAGTGGCAAACTTGACCCCGTCTTTCAGAAGAAAGCGGATGCAGATAAAACTAAAACATTTGCTGaa                                                              >  1:285811/1‑250 (MQ=255)                                                                                                                                                                                     gTTGCTTTAGGTGTCTTTATTGTTAGCGCTCCTGCCTTTTCGCATGGTCATCACTCACACGGCAAGCCCTTAACAGAGGTCGAACAAAAAGCTGCTAATGGTGTTTTTGATGATGCCAATGTACAAAACCGAACGCTCAGTGACTGGGATGGAGTCTGGCAATCCGTTTATCCTTTACTGCAAAGTGGCAAACTTGACCCCGTCTTTCAGAAGAAAGCGGATGCAGATAAAACTAAAACATTTGCTGaa                                                              <  2:285190/249‑1 (MQ=255)                                                                                                                                                                                         cTTTAGGTGTCTTTATTGTTAGCGCTCCTGCCTTTTCGCATGGTCATCACTCACACGGCAAGCCCTTAACAGAGGTCGAACAAAAAGCTGCTAATGGTGTTTTTGATGATGCCAATGTACAAAACCGAACGCTCAGTGACTGGGATGGAGTCTGGCAATCCGTTTATCCTTTACTGCAAAGTGGCAAACTTGACCCCGTCTTTCAGAAGAAAGCGGATGCAGATAAAACTAAAACATTTGCTGAAATTaa                                                         >  1:416809/1‑250 (MQ=255)                                                                                                                                                                                          tttAGGTGTCTTTATTGTTAGCGCTCCTGCCTTTTCGCATGGTCATCACTCACACGGCAAGCCCTTAACAGAGGTCGAACAAAAAGCTGCTAATGGTGTTTTTGATGATGCCAATGTACAAAACCGAACGCTCAGTGACTGGGATGGAGTCTGGCAATCCGTTTATCCTTTACTGCAAAGTGGCAAACTTGACCCCGTCTTTCAGAAGAAAGCGGATGCAGATAAAACTAAAACATTTGCTGAAATTaa                                                         >  1:30769/1‑249 (MQ=255)                                                                                                                                                                                                 gtCTTTATTGTTAGCGCTCCTGCCGTTTCGCATGGTCATCACTCACACGGCAAGCCCTTAACAGAGGTCGAACAAAAAGCTGCTAATGGTGTTTTTGATGATGCCAATGTACAAAACCGAACGCTCAGTGACTGGGATGGAGTCTGGCAATCCGTTTATCCTTTACTGCAAAGTGGCAAACTTGACCCCGTCTTTCAGAAGAAAGCGGATGCCGATAAAACTAAAACCTTTGCTGAAATTAAAGattatt                                                 <  2:121619/250‑1 (MQ=255)                                                                                                                                                                                                  tCTTTATTGTTAGCGCTCCTGCCTTTTCGCATGGTCATCACTCACACGGCAAGCCCTTAACAGAGGTCGAACAAAAAGCTGCTAATGGTGTTTTTGATGATGCCAATGTACAAAACCGAACGCTCAGTGACTGGGATGGAGTCTGGCAATCCGTTTATCCTTTACTGCAAAGTGGCAAACTTGACCCCGTCTTTCAGAAGAAAGCGGATGCAGATAAAACTAAAACATTTGCTGAAATTAAAGATtattat                                               >  1:428502/1‑251 (MQ=255)                                                                                                                                                                                                       aTTGTTAGCGCTCCTGCCTTTTCGCATGGTCATCACTCACACGGCAAGCCCTTAACAGAGGTCGAACAAAAAGCTGCTAATGGTGTTTTTGATGATGCCAATGTACAAAACCGAACGCTCAGTGACTGGGATGGAGTCTGGCAATCCGTTTATCCTTTACTGCAAAGTGGCAAACTTGACCCCGTCTTTCAGAAGAAAGCGGATGCAGATAAAACTAAAACATTTGCTGAAATTAAAGATTATTATcaca                                           <  2:371731/250‑1 (MQ=255)                                                                                                                                                                                                            tAGCGCTCCTGCCTTTTCGCATGGTCATCACTCACACGGCAAGCCCTTAACAGAGGTCGAACAAAAAGCTGCTAATGGTGTTTTTGATGATGCCAATGTACAAAACCGAACGCTCAGTGACTGGGATGGAGTCTGGCAATCCGTTTATCCTTTACTGCAAAGTGGCAAACTTGACCCCGTCTTTCAGAAGAAAGCGGATGCAGATAAAACTAAAACATTTGCTGAAATTAAAGATTATTATCACAAAGGtt                                     >  1:18902/1‑251 (MQ=255)                                                                                                                                                                                                            tAGCGCTCCTGCCTTTTCGCATGGTCATCACTCACACGGCAAGCCCTTAACAGAGGTCGAACAAAAAGCTGCTAATGGTGTTTTTGATGATGCCAATGTACAAAACCGAACGCTCAGTGACTGGGATGGAGTCTGGCAATCCGTTTATCCTTTACTGCAAAGTGGCAAACTTGACCCCGTCTTTCAGAAGAAAGCGGATGCAGATAAAACTAAAACATTTGCTGAAATTAAAGATTATTATCACAAAGGt                                      <  1:100814/250‑1 (MQ=255)                                                                                                                                                                                                                                      agcactcacACGGCAAGCCCTTAACAGAGGTCGAACAAAAAGCTGCTAATGGTGTTTTTGATGATGCCAATGTACAAAACCGAACGCTCAGTGACTGGGATGGAGTCTGGCAATCCGTTTATCCTTTACTGCAAAGTGGCAAACTTGACCCCGTCTTTCAGAAGAAAGCGGATGCAGATAAAACTAAAACATTTGCTGAAATTAAAGATTATTATCACAAAGGTTATGCAACAGATATCGAGATGATTggt           <  2:65533/249‑2 (MQ=255)                                                                                                                                                                                                                                              cacGGCAAGCACTTAACAGACGTCGTACCAAAATCCGCTACTGGTGTTTTTGATGATGCCAATGTACAAAACCGAACGCTCAGTGACTGGGATGGAGTCTGGCAATCCGTTTATCCTTTACTGCAAAGTGGAAAACTTGCCCCCGTCTTTCAGAAGAAAGCGGATGCAGATAAAACTAAAACATTTGATGAAATTAAAGATTATTATCACAAAGGTTATGAAACAGATATCGAGATGATTGGCATTGAgg    >  2:148191/1‑250 (MQ=255)                                                                                                                                                                                                                                               acGGCAAGCCCTTAACAGAGGTCGAACAAAAAGCTGCTAATGGTGTTTTTGATGATGCCAATGTACAAAACCGAACGCTCAGTGACTGGGATGGAGTCTGGCAATCCGTTTATCCTTTACTGCAAAGTGGCAAACTTGACCCCGTCTTTCAGAAGAAAGCGGATGCAGATAAAACTAAAACATTTGCTGAAATTAAAGATTATTATCACAAAGGTTATGCAACAGATATCGAGATGATTGGCATTGAGGAc  <  2:299044/251‑1 (MQ=255)                                                                                                                                                                                                                                                      |                                                                                                                                                                                                                                                    CCTCAAATCACCAAAATAATTTATTTTATTTCAATAAGATAGATAAGTAGAACTGAGAAAGCCATGCTCTCGTTTCCTAAGAGTTGTTGCATTTTGCTATATGTTACAATATAACATTACACATCATATACATTAACTCTGGAGGAAACTGTTTTGGCGATTCGTCTTTACAAACTGGCTGTTGCTTTAGGTGTCTTTATTGTTAGCGCTCCTGCCTTTTCGCATGGTCATCACTCACACGGCAAACCCTTAACAGAGGTCGAACAAAAAGCTGCTAATGGTGTTTTTGATGATGCCAATGTACAAAACCGAACGCTCAGTGACTGGGATGGAGTCTGGCAATCCGTTTATCCTTTACTGCAAAGTGGCAAACTTGACCCCGTCTTTCAGAAGAAAGCGGATGCAGATAAAACTAAAACATTTGCTGAAATTAAAGATTATTATCACAAAGGTTATGCAACAGATATCGAGATGATTGGCATTGAGGAC  >  NC_000913/2041222‑2041710

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