Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	NC_000913	3,411,689	G→T	E13* (GAA→TAA)	yhdJ →	DNA adenine methyltransferase

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	3,411,689	0	G	T	100.0%	84.8 / NA	32	E13* (GAA→TAA)	yhdJ	DNA adenine methyltransferase
Reads supporting (aligned to +/- strand):  ref base G (0/0);  new base T (15/17);  total (15/17)

GAACAGCCCCTGTTGGACATGGTGATGCAATACACCCGTGGTAACCAGACCCGTGCTGCGCTGATGATGGGCATCAACCGTGGTACGCTGCGTAAAAAATTGAAAAAATACGGCATGAACTAATTCAGGTTAGCTAAATGCTTGATTAAAAAGG‑CGCTACTCGGCATGGGGAAGCGCCTTTTTTATAGGTGTCACAAAGGGAGTGACCATGAGAACAGGATGTGAACCGACCCGGTTTGGTAATGAAGCTAAGACCATTATTCACGGTGATGCGCTTGCCGAACTTAAAAAGATCCCCGCCGAAAGTGTCGATCTGATCTTTGCCGACCCACCATATAACATCGGTAAAAATTTTGATGGTCTGATCGAAGCCTGGAAAGAAGATCTGTTTATCGACTGGCTGTTTGAAGTGATTGCAGAGTGCCACCGCGTTCTGAAAAAGCAGGGCAGCATGTACATCATGAACAGTACGGAAAACATGCCCTT  >  NC_000913/3411445‑3411930                                                                                                                                                                                                                                                      |                                                                                                                                                                                                                                                  gAACAGCCCCTGTTGGACATGGTGATGCAATACACCCGTGGTAACCAGACCCGTGCTGCGCTGATGATGGGCATCAACCGTGGTACGCTGCGTAAAAAATTGAAAAAATACGGCATGAACTAATTCAGGTTAGCTAAATGCTTGATTAAAAAGG‑CGCTACTCGGCATGGGGAAGCGCCTTTTTTATAGGTGTCACAAAGGGAGTGACCATGAGAACATGATGTGAACCGACCCTGTTTGGTAATTAAGCTa                                                                                                                                                                                                                                             <  1:192451/251‑1 (MQ=255)    cAGCCCCTGTTGGACATGGTGATGCAATACACCCGTGGTAACCAGACCCGTGCTGCGCTGATGATGGGCATCAACCGTGGTACGCTGCGTAAAAAATTGAAAAAATACGGCATGAACTAATTCAGGTTAGCTAAATGCTTGATTAAAAAGG‑CGCTACTCGGCATGGGGAAGCGCCTTTTTTATAGGTGTCACAAAGGGAGTGACCATGAGAACAGGATGTGAACCGACCCGGTTTGGTAATTAAGCTa                                                                                                                                                                                                                                             >  2:58519/1‑248 (MQ=255)            gTTGGACATGGTGATGCAATACACCCGTGGTAACCAGACCCGTGCTGCGCTGATGATGGGCATCAACCGTGGTACGCTGCGTAAAAAATTGAAAAAATACGGCATGAACTAATTCAGGTTAGCTAAATGCTTGATTAAAAAGG‑CGCTACTCGGCATGGGGAAGCGCCTTTTTTATAGGTGTCACAAAGGGAGTGACCATGAGAACAGGATGTGAACCGACCCGGTTTGGTAATTAAGCTAAGACCattatt                                                                                                                                                                                                                                  <  2:230085/251‑1 (MQ=255)                   aTGGTGATGCAATACACCCGTGGTAACCAGACCCGTGCTGCGCTGATGATGGGCATCAACCGTGGTACGCTGCGTAAAAAATTGAAAAAATACGGCATGAACTAATTCAGGTTAGCTAAATGCTTGATTAAAAAGG‑CGCTACTCGGCATGGGGAAGCGCCTTTTTTATAGGTGTCACAAAGGGAGTGACCATGAGAACAGGATGTGAACCGACCCGGTTTGGTAATTAAGCTAAGACCATTATTCACGGTg                                                                                                                                                                                                                           >  1:115596/1‑251 (MQ=255)                                         gTAACCAGACCCGTGCTGCGCTGATGATGGGCATCAACCGTGGTACGCTGCGTAAAAAATTGAAAAAATACGGCATGAACTAATTCAGGTTAGCTAAATGCTTGATTAAAAAGG‑CGCTACTCGGCATGGGGAAGCGCCTTTTTTATAGGTGTCACAAAGGGAGTGACCATGAGCACAGGATGTGAACCGACCCGGTTTGGTAATTAAGCTAAGACCATTATTCACGGTGATGCGCTTGCCGAACTTaaaa                                                                                                                                                                                                      >  2:173293/1‑250 (MQ=255)                                                                         aTCAAACGTGGTACGCTGCGTAAAAAATTGAAAAAATACGGCATGAACTAAGTCAGGTTAGCTAAATGCTTG‑TTAAAAAGGGCTGTACTCGGCAGGGGGAGGCGCCTTTTTTATAGGTGTCACAAAGGGAGTGACCGTGAGAACAGGTTGTGAACCGACCCGGTTGGGTAATTAAGCTAAGACCATGATACACGGGGATGCGCTGGCGGAACTAAAAAAGATCCCCGCCGAAAGTGTCGATCTGATCTTTg                                                                                                                                                                     <  2:209639/251‑1 (MQ=255)                                                                              ccGTGGTACGCTGCGTAAAAAATTGAAACAATACGGCATGCACTAATTCAGGTTAGCTAAATGCTTGATTAAAAAGG‑CGCTACTCGGCATGGGGAAGCGCCTTTTTTATAGGTGTCACAAAGGGAGTGACCATGAGAACAGGATGTGAACCGACCCGGTTTGGTAATTAAGCTAAGACCATTATTCACGGTGATGCGCTTGCCGAACTTAAAAAGATCCCCGCCGAAAGTGTCGATCTGCTCTTTGCCGAc                                                                                                                                                                >  2:210382/1‑251 (MQ=255)                                                                               cGTGGTACGCTGCGTAAAAAATTGAAAAAATACGGCATGAACTAATTCAGGTAAGCAAAATGCTTGATTAAAAAGG‑CGCTACTCGGCATGGGGAAGCGCCTTTTTTAAAGGTGTCACAAAGGGAGTGACCATGAGAACAGGATGTGAACCGACCCGGTTTGGTAATTAAGCTAAGACCATTATTCAGGGTGATGCGCTTGCCGAACTTAAAAAGATCCCCGCCGAAAGTGTCGATCTGATCTTTGCCGAcc                                                                                                                                                               <  2:84523/251‑1 (MQ=255)                                                                                     aCGCTGCGTAAAAAATTGAAAACATACGGAATGAACTAATTCCGGTTAGCTAAATGCTTGATTAAAAAGG‑CGCTACTCGGCATGGGGAAGCGCCTTTTTTATAGGTGTCACAAAGGGAGTGACCATGAGAACAGGATGTGAACCGACCCGGTTTGGTAATTAAGCTAAGACCATTATTCACGGTGATGCGCTTGCCGAACTTAAAAAGATCCCCGCCGAAAGAGTCGATCTGATCTTTGCCGACCCACCat                                                                                                                                                         >  2:231649/1‑251 (MQ=255)                                                                                          gCGTAAAAAATTGAAAAAATACGGCATGAACTAATTCAGGTTAGCTAAATGCTTGATTAAAAAGG‑CGCTACTCGGCATGGGGAAGCGCCTTTTTTATAGGTGTCACAAAGGGAGTGACCATGAGAACAGGATGTGAACCGACCCGGTTTGGTAATTAAGCTAAGACCATTATTCACGGTGATGCGCTTGCCGAACTTAAAAAGATCCCCGCCGAAAGTGTCGATCTGATCTTTGCCTACCCACCAtata                                                                                                                                                      <  1:181861/249‑1 (MQ=255)                                                                                          gCGTAAAAAATTGAAAAAATACGGCATGAACTAATTCAGGTTAGCTAAATGCTTGATTAAAAAGG‑CGCTACTCGGCATGGGGAAGCGCCTTTTTTATAGGTGTCACAAAGGGAGTGACCATGAGAACAGGATGTGAACCGACCCGGTTTGGTAATTAAGCTAAGACCATTATTCACGGTGATGCGCTTGCCGAACTTAAAAAGATCCCCGCCGAAAGTGTCGATCTGATCTTTGCCGACCCACCATATaa                                                                                                                                                     <  1:171676/250‑1 (MQ=255)                                                                                                          aaaTACGGCATGAACTAATTCAGGTTAGCTAAATGCTTGATTAAAAAGG‑CGCTACTCGGCATGGGGAAGCGCCTTTTTTATAGGTGTCACAAAGGGAGTGACCATGAGAACAGGATGTGAACCGACCCGGTTTGGTAATTAAGCTAAGACCATTATTCACGGTGATGCGCTTGCCGAACTTAAAAAGATCCCCGCCGAAAGTGTCGATCTGATCTTTGCCGACCCACCATATAACATCGGTAAAAATTTTg                                                                                                                                    <  1:147411/251‑1 (MQ=255)                                                                                                                ggCATGAAATAATTCAGGTTAGCTAAATGCTTGCTTTCAAGGG‑CGCTCCTCGGCATGGGGAAGAGCCTTTTTTATAGGTGTCACACAGGGAGTGACCATGAGAACAGGATGTGAACCGACCCGGTTTGGTAATTAAGCTAAGAACATTATTCACGGTGATGCGCTTGCCGAACTTAACAAGATCCCCGCCGAAAGAGTCGATCTGATCTTTGCCGACCCACCATATAACATCGGTAAAACTTTTGATgatc                                                                                                                              >  2:238452/1‑248 (MQ=255)                                                                                                                      aaCTAATTCAGGTTAGCTAAATGCTTGATTAAAAAGG‑CGCTACTCGGCATGGGGAAGCGCCTTTTTTATAGGTGTCACAAAGGGAGTGACCATGAGAACAGGATGTGAACCGACCCGGTTTGGTAATTAAGCTAAGACCATTATTCACGGTGATGCGCTTGCCGAACTTAAAAAGATCCCCGCCGAAAGTGTCGATCTGATCTTTGCCGACCCACCATATAACATCGGTAAAAATTTTGATGGTCTGATCg                                                                                                                        >  1:238266/1‑251 (MQ=255)                                                                                                                           aTTCAGGTTAGCTAAATGCTTGATTAAAAAGG‑CGCTACTCGGCATGGGGAAGCGCCTTTTTTATAGGTGTCACAAAGGGAGTGACCATGAGAACAGGATGTGAACCGACCCGGTTTGGTAATTAAGCTAAGACCATTATTCACGGTGATGCGCTTGCCGAACTTAAAAAGATCCCCGCCGAAAGTGTCGATCTGATCTTTGCCGACCCACCATATAACATCGGTAAAAATTTTGATGGTCTGATCGAAGcc                                                                                                                   >  2:21796/1‑251 (MQ=255)                                                                                                                                 gTTAGCTAAATGCTGGAGTAAAAAGG‑CGCGACTCGGCATGGGGAAGCGCCTTTTTTATAGGTTTCACAAAGGGAGTGACCATGAGAACAGGATTTGAACCGACCCGTTGTGTTAATTAAGCTAAGACCATGATTCGCGGGGATGCGCTTGCCGAGCTTAAAAAGATCCCCGCCGAAAGTGTCGATCTGATCTTTGCCGACCCACCATATAACATCGGTAAAAATTTTGATGGTCTGATCGAAGCCTGGaaa                                                                                                             <  1:143949/251‑1 (MQ=255)                                                                                                                                     gCTAAATGCTTGATTAAAAAGG‑CGCTACTCGGCATGGGTAAGCGCCTTTTTTATAGGTTTCACAAAGGGAGTGACCAGGAGAACAGGATGTGAACCGACCGGGTTTGGTAATTAAGCTAAGACCATTATTCACGGTGATGCGTTGGCCGAACTTAAAAAGATCCCCGCCGAAAGCGTCGATCTGATCTTTGCCGACCCACCATATAACATCGGTAAAAATTTTGATGGTCTGATCGAAGCCTGGAaagaag                                                                                                         <  2:115088/251‑1 (MQ=255)                                                                                                                                      cTAAATGCTTGATTAAAAAGG‑CGCTACTCGGCCTGGGGAAGCGCCTTTTTTATAGGTGTCACAAAGGGAGTGACCATGAGAACAGGATGTGAACCGACCCGGTTTGGTAATTAAGCTAAGACCATTATTCACGGTGATGCGCTTGCCGAACTTAAAAAGATCCCCGCCGAAAGTGTCGATCTGATCTTTGCCGACCCACCATATAACATCGGTAAAAATTTTGATGGTCTGATCGAAGCCTGGAaaga                                                                                                           >  1:56931/1‑248 (MQ=255)                                                                                                                                      cTAAATGCTTGATTAAAAAGG‑CGCTACTCGGCATGGGGAAGCGCCTTTTTTATAGGTGTCACAAAGGGAGTGACCATGAGAACAGGATGTGAACCGACCCGGTTTGGTAATTAAGCTAAGACCATTATTCACGGTGATGCGCTTGCCGAACTTAAAAAGATCCCCGCCGAACGTGTCGATCTGATCTTTGCCGACCCACCATATAACATCGGTACAAATTTTGATGGTCTGATCGCAGCCTGGAaaga                                                                                                           >  2:184422/1‑248 (MQ=255)                                                                                                                                       tAAATGCTTGATTAAAAAGG‑CGCTACTCGGCATGGGGAAGCGCCTTTTTTATAGGTGTCACAAAGGGAGTGACCATGAGAACAGGATGTGAACCGACCCGGTTTGGTAATTAAGCTAAGACCATTATTCACGGTGATGCGCTTGCCGAACTTAAAAAGATCCCCGCCGAAAGTGTCGATCTGATCTTTGCCGACCCACCATATAACATCGGTAAAAATTTTGATGGTCTGATCGAAGCTGGGAaaga                                                                                                           >  2:105735/1‑247 (MQ=255)                                                                                                                                              ttttttAAAAAGG‑CTCTACCCGGCATGGGGAAGCGCTTTTTTTATAGGTGTCACAAAGGGAGTACCCATGAGAACAGGTTGTGAACCGACCCGTTTTGGTAATTAACCTAAGCCCATTATTCACGGGGATGCGCTTGCCGAACTTAAAAAGATCCCCGCCGAAAGTGTCGATCTGATCTTTGCCGACCCACCATATAACATCGGTAAAAATTTTGATGGTCTGCTCGAAGCCTGGAAAGAAGATCTGttt                                                                                                 <  2:3189/246‑1 (MQ=255)                                                                                                                                                gATTAAAACGG‑CGCTACTCGGCATGGGGAAGCGCTTTTTTTATAGGTGTCCCAAAGGGAGTGACCATGAGAACAGGATGTGAACCGACCCGGTTTGGTAATTAAGCTAAGACCATTATTCACGGTGATGCGCTTGCCGAACTTAAAAAGATCCCCGCCGAAAGTGTCGATCTGATCTTTGCCGACCCACCATATAACATCGGTAAAAATTTTGATGGTCTGATCGAAGCCTGGAAAGAAGATCTGTTTATc                                                                                              <  1:231959/251‑1 (MQ=255)                                                                                                                                                                               gcgcCTTTTTTATAGGTGTCACAAAGGGAGTGACCATGAGAACAGGATGTGAACCGACCCGGTTTGGTAATTAAGCTAAGACCATTATTCACGGTGATGCGCTTGCCGAACTTAAAAAGATCCCCGCCGAAAGTGTCGATCTGATCTTTGCCGACCCACCATATAACATCGGTAAAAATTTTGATGGTCTGATCGAAGCCTGGAAAGAAGATCTGTTTATCGACTGGCTGTTTGAAGTGATTGCAGAGTGc                                                                >  1:197516/1‑251 (MQ=255)                                                                                                                                                                                           aaggtgtagcaAAGGGAGTAACCATGAGAACAGGATGTGAACCGACCCGTTTGGGTAATTAAGCTAAGCCCATTATTCACGGTGATGCGCTTTCCGAACTTAAAAAGTTCCCCGCCGAAAGTGTCGATCTGATCTTTGCCGACCCACCATATAACATCGGTAAAAATTTTGATGGTCTGATCGAAGCCTGGAAAGAAGATCTGTTTATCGACTGGCTGTTTGAAGTGATTGCAGAGTGCCACCGCGTTCTg                                                    <  2:209388/242‑1 (MQ=255)                                                                                                                                                                                             gggctCACAAAGGGAGTAACAATGAGAACAGGATGTGAACCGACCCGGTTTGGTAATTAAGCTAAAACCATTATTCACGGTGATGCGCTTGCCGAACTTAAAAAGATCCCCGCCGAAAGTGTCGATCTGATCTTTGCCGACCCACCATATAACATCGGTAAAAATTTTGATGGTCTGATCGAAGCCTGGAAAGAAGATCTGTTTATCGACGGGCTGTTTGAAGTGATTGCAGAGTGCCACCGCGTTCTGaa                                                  <  1:199236/247‑1 (MQ=255)                                                                                                                                                                                                  aaaaaggggAGTGACAATGAGACAAGAATGGTAACCGACCAGGTTTGGTAATTAACCTAAGACAATTATGTACGGTGATGCGCTGGCCGAACTTAAAAAGATCTCCGCCGAAAGTGTCGAGCTGATCTTTGCCGACCCACCATATAACATCGGTAAAAATTTTGATGGTCTGATCGAAGCCTGGAAAGAAGATCTGTGTATCGACTGGCTGTTTGAAGTGATTGCAGAGTTCCACCGCGTTCTGAAAAAGc                                             <  2:221854/245‑1 (MQ=255)                                                                                                                                                                                                   acaAAGGGAGGGACCATAAGAACAGGATGTGAACCAACCCGGTTTGGTAATTAAGCTAATACCATTATTCACGGTGATGCGCTTGCTGAACTTAAAAAGATCCCCGCCGAAAGTGTCGATCTGATCTGTGCCGACCCACCATATAACATCGGTAAAAATTTTGATGGTCTGATCGAAGCCTGGAAAGAATATCGGTTTATCGGGTGGCTGTGTTAAGTTATTGCAGAGTGCCACCGCGTTCTGAAAAAGCa                                            <  1:18681/251‑1 (MQ=255)                                                                                                                                                                                                                          aGGATGTGAACCGACCCGGTTTGGTAATTAAGCTAAGACCATTATTCACGGTGATGCGCTTGCCGAACTTAAAAAGATCCCCGCCGAAAGTGTCGATCTGATCTTTGCCGACCCACCATATAACATCGGTAAAAATTTTGATGGTCTGATCGAAGCCTGGAAAGAAGATCTGTTTATCGACTGGCTGTTTGAAGTGATTGCAGAGTGCCACCGCGTTCTGCAAAAGCAGGGCAGCATGTACATCATGAAcc                     >  2:87572/1‑250 (MQ=255)                                                                                                                                                                                                                           gtaggggAACTGACCCGGTTTGGTAATTAAGCAAAGACCATTATTCACGGTGATGCGCTTGCCGAACTTAAAAAGATCCCCGCCGAAAGTGTCGATCTGATCTTTGCCGACCCACCATATAACATCGGTAAAAATTTTGATGGTCTGATCGAAGCCTGGAAAGAAGATCTGTTTATCGACTGGCTGTTTGAAGTGATTGCAGAGTGCCACCGCGTTCTGAAAAAGCAGGGCAGCATGTACATCATGAACAg                    <  2:199413/245‑1 (MQ=255)                                                                                                                                                                                                                             atgggAACCGACCCGGTTTGGTAATTAAGCTAAGACCATGATGCACGGTGATGCGCTTGCCGAACTTAAAAAGATCCCCGCCGAAAGTGTCGATCTGATCTTTGCCGACCCACCATATAACATCGGTAAAAATTTTGATGGTCTGATCGAAGCCTGGAAAGAAGATCTGTTTATCGACTGGCTGTTTGAAGTGATTGCAGAGTGCCACCGCGTTCTGAAAAAGCAGGGCAGCATGTACATCATGAACAGTa                  <  2:54863/247‑1 (MQ=255)                                                                                                                                                                                                                                tgAACCGACCCGGTTTGGTAATTAAGCTAAGACCATTATTCACGGTGATGCGCTTGCCGAACTTAAAAAGATCCCCGCCGAAAGTGTCGATCTGATCTTTGCCGACCCACCATATAACATCGGTAAAAATTTTGATGGTCTGATCGAAGCCTGGAAAGAAGATCTGTTTATCGACTGGCTGTTTGAAGTGATTGCAGAGTGCCACCGCGTTCTGAAAAAGCAGGGCAGCATGTACATCATGAACAGTACgg               >  1:210521/1‑251 (MQ=255)                                                                                                                                                                                                                                             tttGGTAATTAAGCTAAGACCATTATTCACGGTGATGCGCTTGCCGAACTTAAAAAGATCCCCGCCGAAAGTGTCGATCTGATCTTTGCCGACCCACCATATAACATCGGTAAAAATTTTGATGGTCTGATCGAAGCCTGGAAAGAAGATCTGTTTATCGACTGGCTGTTTGAAGTGATTGCAGAGTGCCACCGCGTTCTGAAAAAGCAGGGCAGCATGTACATCATGAACAGTACGGAAAACATGCCCtt  >  1:63432/1‑251 (MQ=255)                                                                                                                                                                                                                                                      |                                                                                                                                                                                                                                                  GAACAGCCCCTGTTGGACATGGTGATGCAATACACCCGTGGTAACCAGACCCGTGCTGCGCTGATGATGGGCATCAACCGTGGTACGCTGCGTAAAAAATTGAAAAAATACGGCATGAACTAATTCAGGTTAGCTAAATGCTTGATTAAAAAGG‑CGCTACTCGGCATGGGGAAGCGCCTTTTTTATAGGTGTCACAAAGGGAGTGACCATGAGAACAGGATGTGAACCGACCCGGTTTGGTAATGAAGCTAAGACCATTATTCACGGTGATGCGCTTGCCGAACTTAAAAAGATCCCCGCCGAAAGTGTCGATCTGATCTTTGCCGACCCACCATATAACATCGGTAAAAATTTTGATGGTCTGATCGAAGCCTGGAAAGAAGATCTGTTTATCGACTGGCTGTTTGAAGTGATTGCAGAGTGCCACCGCGTTCTGAAAAAGCAGGGCAGCATGTACATCATGAACAGTACGGAAAACATGCCCTT  >  NC_000913/3411445‑3411930

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