Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	NC_000913	3,517,171	G→A	L46F (CTT→TTT)	damX ←	cell division protein DamX

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	3,517,171	0	G	A	100.0%	100.0 / NA	33	L46F (CTT→TTT)	damX	cell division protein DamX
Reads supporting (aligned to +/- strand):  ref base G (0/0);  new base A (13/20);  total (13/20)

GAGGGGGCTTTTAGCGCAGAACCGATACCGATGATCAACAGCAGTAGAACCAGAATGCCGACGCCCATCATCATATACTGACGAGAAGCGGGTTTACTGGCTGCTTTTTTGCGCTTACGCGGACGACGCTCTACGCGCTCTTCATCCACGGTTTCATCTTCGGATTCGTCAATTTCTTCTTCGATTTCCGGTTCCTCATTGCGCTCTTTTTGCGCACGAGTCGGACGGCGATCGTCAGTGTCATCAAGTTCAATATCATCAAAATTGATCTGCGGTTCGCCACGTTCAGTACGCTCAGAACGTTCAGAAGATTGACGAGAACGACCAGTACGACGATCGCTGGGATCGGGTTTCAGCTCGTCTTCTGGTTTGAATTCATCCATTTAACACCCCACTAAAAGGTTAATGCTTACCACGTTGCAATTAACCTGAAGCTAATAGACCGCTTGATAAGCGGCCTGACCTTTCTTGTTGTTACGCTGATTGACAA  >  NC_000913/3516924‑3517413                                                                                                                                                                                                                                                        |                                                                                                                                                                                                                                                   gAGGGGGCTTTTAGCGCAGAACCGATCCCGATGATCAACAGCAGTAGAACCAGAATGCCGACGCCCATCATCATATACGGACGAGAAGCGGGTTTACTGGCTGCTTTTTTGCGCTTACGCGGACGACGCTCTACGCGCTCTTCATCCACGGTTTCATCTTCGGATTCGTCAATTTCTTCTTCGATTTCCGGTTCCTCATTGCGCTCTTTTTGCGCACGAGTCGGACGGCGATCGTCAGTGTCATCAAATTc                                                                                                                                                                                                                                                 <  2:74125/251‑1 (MQ=255)    ggggCTTTTAGCGCAGAACCGATACCGATGATCAACAGCAGTAGAACCAGAATGCCGACGCCCATCATCATATACTGACGAGAAGCGGGTTTACTGGCTGCTTTTTTGCGCTTACGCGGACGACGCTCTACGCGCTCTTCATCCACGGTTTCATCTTCGGATTCGTCAATTTCTTCTTCGATTTCCGGTTCCTCATTGCGCTCTTTTTGCGCACGAGTCGGACGGCGATCGTCAGTGTCATCAAATTCAat                                                                                                                                                                                                                                              >  1:256678/1‑251 (MQ=255)          tttAGCGCAGAACCGATACCGATGATCAACAGCAGTAGAACCAGAATGCCGACGCCCATCATCATATACTGACGAGAAGCGGGTTTACTGGCTGCTTTTTTGCGCTTACGCGGACGACGCTCTACGCGCTCTTCATCCACGGTTTCATCTTCGGATTCGTCAATTTCTTCTTCGATTTCCGGTTCCTCATTGCGCTCTTTTTGCGCACGAGTCGGACGGCGATCGTCAGTGTCATCAAATTCAATatcat                                                                                                                                                                                                                                         <  1:180899/250‑1 (MQ=255)                     cccgataccgatGATCAACAGCAGTAGAACCAGAATGCCGACGCCCATCATCATATACTGACGAGAAGCGGGTTTACTGGCTGCTTTTTTGCGCTTACGCGGACGACGCTCTACGCGCTCTTCATCCACGGTTTCATCTTCGGATTCGTCAATTTCTTCTTCGATTTCCGGTTCCTCATTGCGCTCTTTTTGCGCACGAGTCGGACGGCGATCGTCAGTGTCATCAAATTCAATATCATCAAAATTGATCt                                                                                                                                                                                                                             <  1:14419/250‑1 (MQ=255)                        gataccgatGATCAACAGCAGTAGAACCAGAATGCCGACGCCCATCATCATATACTGACGAGAAGCGGGTTTACTGGCTGCTTTTTTGCGCTTACGCGGACGACGCTCTACGCGCTCTTCATCCACGGTTTCATCTTCGGATTCGTCAATTTCTTCTTCGATTTCCGGTTCCTCATTGCGCTCTTTTTGCGCACGAGTCGGACGGCGATCGTCAGTGTCATCAAATTCAATATCATCAAAATTGATCTGCg                                                                                                                                                                                                                          >  1:205227/1‑251 (MQ=255)                        gataccgatGATCAACAGCAGTAGAACCAGAATGCCGACGCCCATCAGCATATACTGACGAGAAGCGGGTTTACTGGCTGCTTTTTTGCGCTTACGCGGACGACGCTCTACGCGCTCTTCATCCACGGTTTCATCTTCGGATTCGTCAATTTCTTCTTCGATTTCCGGTTCCTCATTGCGCTCTTGTTGCGCACGAGTCGGACGGCGATCGTCAGTGTCATCAAATTCAATATCATCAAAAGTGATCTGCg                                                                                                                                                                                                                          <  2:250866/251‑1 (MQ=255)                           accgatGATCAACAGCAGTAGAACCAGAATGCCGACGCCCATCATCATATACTGACGAGAAGCGGGTTTACTGGCTGCTTTTTTGCGCTTACGCGGACGACGCTCTACGCGCTCTTCATCCACGGTTTCATCTTCGGATTCGTCAATTTCTTCTTCGATTTCCGGTTCCTCATTGCGCTCTTTTTGCGCACGAGTCGGACGGCGATCGTCAGTGTCATCAAATTCAATATCATCAAAATTGATCTGCGGtt                                                                                                                                                                                                                       >  2:101164/1‑251 (MQ=255)                           accgatGATCAACAGCAGTAGAACCAGAATGCCGACGCCCATCATCATATACTGACGAGAAGCGGGTTTACTGGCTGCTTTTTTGCGCTTACGCGGACGACGCTCTACGCGCTCTTCATCCACGGTTTCATCTTCGGATTCGTCAATTTCTTCTTCGATTTCCGGTTCCTCATTGCGCTCTTTTTGCGCACGAGTCGGACGGCGATCGTCAGTGTCATCAAATTCAATATCATCAAAATTGATCTGCGGt                                                                                                                                                                                                                        <  2:141923/250‑1 (MQ=255)                                    cAACAGCAGTAGAACCAGAATGCCGACGCCCATCATCATATACTGACGAGAAGCGGGTTTACTGGCTGCTTTTTTGCGCTTACGCGGACGACGCTCTACGCGCTCTTCATCCACGGTTTCATCTTCGGATTCGTCAATTTCTTCTTCGATTTCCGGTTCCTCATTGCGCTCTTTTTGCGCACGAGTCGGACGGCGATCGTCAGTGTCATCAAATTCAATATCATCAAAATTGATCTGCGGTTCGCCACGtt                                                                                                                                                                                                              <  1:36921/251‑1 (MQ=255)                                              aGAACCAGAATGCCGTCGCCCATCATCATATACTTACGAGAAGCGGGTTTACTGGCTGCTTTTTTGCGCTTACTCGGACGACGCTCTACGCGCTCTTCATCCACGGTTTCATCTTCGGAGTCGTCAATTTCTTCTTCGATTTCCGGTTCCTCATTGCGCGCGTTTTGCGCACGAGTCGGACGGCGATCGGCAGTGTCATCAAATTCAATGTCATCAAAATTGAGCTGGGGGTCTCCACGTTCAGTTCGCTc                                                                                                                                                                                                    <  2:94321/251‑1 (MQ=255)                                                   cAGAATGCCGACGCCCATCATCATATACTGACGAGAAGCGGGTTTACTGGCTGCTTTTTTGCGCTTACGCGGACGACGCTCTACGCGCTCTTCATCCACGGTTTCATCTTCGGATTCGTCAATTTCTTCTTCGATTTCCGGTTCCTCATTGCGCTCTTTTTGCGCACGAGTCGGACGGCGATCGTCAGTGTCATCAAATTCAATATCATCAAAATTTATCTGCGGTTCGCCACGTTCAGTACGCTCAGAAc                                                                                                                                                                                               <  1:101164/251‑1 (MQ=255)                                                                        catATACTGACGAGAAGCGGGTTTACTGGCTGCTTTTTTGCGCTTACGCGGACGACGCTCTACGCGCTCTTCATCCACGGTTTCATCTTCGGATTCGTCAATTTCTTCTTCGATTTCCGGTTCCTCATTGCGCTCTTTTTGCGCACGAGTCGGACGGCGATCGTCAGTGTCATCAAATTCAATATCATCAAAATTGATCTGCGGTTCGCCACGTTCAGTACGCTCAGAACGTTCAGACGATTGACGCGAAc                                                                                                                                                                          >  1:239947/1‑251 (MQ=255)                                                                         atatAATGACGAGAAGCGGGTTTACTGGCTGCTTTTTTGCGCTTACGCGGACGACGCTCTACGCGCTCTTCATCCACGGTTTCATCTTCGGATTCGTCAATTTCTTCTTCGATTTCCGGTTCCTCATTGCGCTCTTTTTGCGCACGAGTCGGACGGCGATCGTCAGTGTCATCAAATTCAATATCATCAAAATTGATCTGCGGTTCGCCACGTTCAGTACGCTCAGAACGTTCAGAAGATTGACGAGAACg                                                                                                                                                                         >  2:82456/1‑251 (MQ=255)                                                                                gACGAGAAGCGGGTTTACTGGCTGCTTTTTTGCGCTTACGCGGACGACGCTCTACGCTCTCTTCATCCACGGTTTCATCTTCGGATTCGTCAATTTCTTCTTCGATTTCCGGTTCCTCATTGCGCTCTTTTTGCGCACGAGTCGGACGGCGATCGTCAGTGTCATCAAATTCAATATCATCAAAATTGATCTGCGGTTCGCCACGTTCAGTACGCTCAGAACGTTCAGAAGATTGACGAGAACGACCAGt                                                                                                                                                                   <  1:152751/250‑1 (MQ=255)                                                                                                aCTGGCTGCTTTTTTGCGCTTACGCGGACGACGCTCTACGCGCTCTTCATCCACGGTTTCATCTTCGGATTCGTCAATTTCTTCTTCGATTTCCGGTTCCTCATTGCGCTCTTTTTGCGCACGAGTCGGACGGCGATCGTCAGTGTCATCAAATTCAATATCATCAAAATTGATCTGCGGTTCGCCACGTTCAGTACGCTCAGAACGTTCAGAAGATTGACGAGAACGACCAGTACGACGATCGCTGGGAt                                                                                                                                                  >  1:266968/1‑251 (MQ=255)                                                                                                   ggCTGATTTTTTGCGCTTACGCGGACGACGCTCTACGCGCTCTTCATCCACGGTTTCATCTTCGGATTCGTCAATTTCTTCTTCGATTTCCGGTTCCTCCTTGCGCTCTTTTTGCGCACGAGTCGGACGGCGATCGTCAGTGTCATCAAATTCAATATCATCAAAATTGATCTGCGGTTCGCCACGTTCAGTACGCTAAGAACGTTCAGCAGATTGCCGAGAACGACACGTACGACGATCGCTGGGATCgg                                                                                                                                               >  2:15410/1‑251 (MQ=255)                                                                                                       gctTTTTTGCGCTTACGCGGACGACGCTCTACGCGCTCTTCATCCACGGTTTCATCTTCGGATTCGTCAATTTCTTCTTCGATTTCCGGTTCCTCATTGCGCTCTTTTTGCGCACGAGTCGGACGGCGATCGTCAGTGTCATCAAATTCAATATCATCAAAATTGATCTGCGGTTCGCCACGTTCAGTACGCTCAGAACGTTCAGAAGATTGACGAGAACGACCAGTACGACGATCGCTGGGATCGGGtt                                                                                                                                            <  2:48147/250‑1 (MQ=255)                                                                                                                      cgcgGACGACGCTCTACGCGCTCTTCATCCACGGTTTCATCTTCTGATTCGTCAATTTCTTCTTCGATTTCCGGTTCCTCATTGCGCTCTTTTTGCGCACGAGTCGGACGGCGATCGTCAGGGTCATCAAATTCAATAGCATCAAAATTGATCTGCGGTTCGCCACGTTCAGTACGCTCAGAACGTTCAGAAGATTTACGAGAACGACCAGTACGACGATCGCTGGGATCGGGTTTCAGCGCGTCTTCTgg                                                                                                                            <  2:124626/251‑1 (MQ=255)                                                                                                                      cgcgGACGACGCTCTACGCGCTCTTCATCCACGGTTTCATCTTCGGATTCGTCAATTTCTTCTTCGATTTCCGGTTCCTCATTGCGCTCTTTTTGCGCACGAGTCGGACGGCGATCGTCAGTGTCATCAAATTCAATATCATCAAAAGTGATCTGCGGTTCGCCACGTTCAGTACGCTCAGAACGTTCAGAAGATTGACGAGAACGACCAGTGCGACGATCGCTGGGATCGGGTTTCAGCTCGTCTTCTgg                                                                                                                            <  2:53401/251‑1 (MQ=255)                                                                                                                             gacgCTCTACGCGCTCTTCATCCACGGTTTCATCTTCGGATTCGTCAATTTCTTCTTCGATTTCCGGTTCCTCATTGCGCTCTTTTTGCGCACGAGTCGGACGGCGATCGTCAGTGTCATCAAATTCAATATCATCAAAATTGATCTGCGGTTCGCCACGTTCAGTACGCTCAGAACGTTCAGAAGATTGACGAGAACGACCAGTACGACGATCGCTGGGATCGGGTTTCAGCTCGTCTTCTGGTTTGAAt                                                                                                                     >  2:60606/1‑251 (MQ=255)                                                                                                                                                        tttCATCTTCGGATTCGTCAATTTCTTCTTCGATTTCCGGTTCCTCATTGCGCTCTTTTTGCGCACGAGTCGGACGGCGATCGTCAGTGTCATCAAATTCAATATCATCAAAATTGATCTGCGGTTCGCCACGTTCAGTACGCTCAGAACGTTCAGAAGATTGACGAGAACGACCAGTACGACGATCGCTGGGATCGGGTTTCAGCTCGTCTTCTGGTTTGAATTCATCCATTTAACACCCCACTACAAgg                                                                                          >  2:32344/1‑251 (MQ=255)                                                                                                                                                           cATCTTCGGATTCGTCAATTTCTTCTTCGATTTCCGGTTCCTCATTGCGCTCTTTTTGCGCACGAGGCGGACGGCGATCGTCAGTGTCATCAAATTCAATATCATCAAAATTGATCTGCGGTTCGCCACGTTCAGTACGCTCAGAACGTTCAGAAGATTGACGAGAACGACCAGTACGACGATCGCTGGGATCGGGTTTCAGCTCGTCTTCTGGTTTGAATTCATCCATTTAACACCCCACTAAAAGGtt                                                                                        <  1:143103/250‑1 (MQ=255)                                                                                                                                                                     ttCGTCTATTTCCTCGTCGAGTTCCGGTTCCTCATTGCGCTCTATTTGCGCACGAGTCGGACGGCGATCGTCAGTGTCATCAAATTCAATATCATCAAAATTGAACTGCGGTTCTCCACGTTCAGTACGCTCATAACGGTCAGAAGATTGACGAGAACGACCAGTACGACATTCGCTTGGATCGCGCTTCAGCTCGTCTTCTGGTTTGAATTCATCCAGTTAACACCCCACTAAAAGGGTAATGCttgcc                                                                              <  2:43599/250‑4 (MQ=255)                                                                                                                                                                       cGTCAATTTCTTCTTCGATTTCCGGTTCCTCATTGCGCTCTTTTTGCGCACGAGTCGGACGGCGATCGTCAGTGTCATCAAATTCAATATCATCAAAATTGACCTGCGGTTCGCCACGTTCAGTACGCTCAGAACGTTCAGAAGATTGACGAGAACGACCAGTACGACGATCGCTGGGATCGGGTTTCAGCTCGTCTTCTGGTTTGAATTCATCCATTTAACACCCCACTAAAAGGTTAATGCTTACCACg                                                                           <  1:129615/251‑1 (MQ=255)                                                                                                                                                                         tCAATTTCTTCTTCGATTTCCGGTTCCTCATTGCGCTCTTTTTGCGCACGAGTCGGACGGCGATCGTCAGTGTCATCAAATTCAATATCATCAAAATTGATCTGCGGTTCGCCACGTTCAGTACGCTCAGAACGTTCAGAAGATTGACGAGAACGACCAGTACGACGATCGCTGGGATCGGGTTTCAGCTCGTCTTCTGGTTTGAATTCATCCATTTAACACCCCACTAAAAGGTTAATGCTTACCACGtt                                                                         >  1:266475/1‑251 (MQ=255)                                                                                                                                                                                        ttttCCGGTTCCTCATTGCGCTCTTTATGAGCACGAGTCGGACGGCGATCGTCAGTGTCATCCAATTCAATATCATCAAAATTGATCTGCGGTTCGCCACGTTCAGTACGCTCAGAACGTTCAGATGATTGACGAGAACGACCAGTACGACGATCGCTGGGATCGGGTTTCAGCTCGTCTTCTGGTTTGAATTCATCCATTTAACACCCCACTACAAGGTTAATGCTTACCACGTTGCAATTAATCTGa                                                            >  2:127400/2‑249 (MQ=255)                                                                                                                                                                                                  ccTCATTGCGCTCTTTTTGCGCACGAGTCGGACGGCGATCGTCAGTGTCATCAAATTCAATATCATCAAAATTGATCTGCGGTTCGCCACGTTCAGTACGCTCAGAACGTTCAGAAGATTGACGAGAACGACCAGTACGACGATCGCTGGGATCGGGTTTCAGCTCGTCTTCTGGTTTGAATTCATCCATTTAACACCCCACTAAAAGGTTAATGCTTACCACGTTGCAATTAACCTGAAGCTAATAGAcc                                                <  2:205227/251‑1 (MQ=255)                                                                                                                                                                                                     cATTGCGCTCTTTTTGCGCTCGAGTCGGACGGCGATCGTCAGTGTCATCAAATTCAATATCATCAAAATTGAGCTGCGGTTCGCCACGTTCAGTACGCTCAGAACGTTCAGAAGATTGACGAGAACGACCAGTACGACGATCGCTGGGATCGGGTTTCAGCTCGTCTTCTGGTTTGAATTCATCCATTTAACACCCCACTAAAAGGTTAATGCTTACCACGTTGCAATTAACCTGAAGCTAATAGACCGc                                              <  2:5103/250‑1 (MQ=255)                                                                                                                                                                                                      aTTGCGCTCTTTTTGCGCACGAGTCGGACGGCGATCGTCAGTGTCATCAAATTCAATATCATCAAAATTGATCTGCGGTTCGCCACGTTCAGTACGCTCAGAACGTTCAGAAGATTGACGAGAACGACCAGTACGACGATCGCTGGGATCGGGTTTCAGCTCGTCTTCTGGTTTGAATTCATCCATTTAACACCCCACTAAAAGGTTCATGCTTACCACGTTGCAATTAACCTGAAGCTAATAGACCGCtt                                            >  2:47561/1‑251 (MQ=255)                                                                                                                                                                                                                                            cAGTGTCATCAAATTCAATATCATCAAAATTGATCTGCGGTTCGCCACGTTCAGTACGCTCAGAACGTTCAGAAGATTGACGAGAACGACCAGTACGACGATCGCTGGTATCGGGTTTCAGCTCGTCTTCTGGTTTGAATTCATCCATTTAACACCCCACTAAAAGGTTAATGCTTACCTCGTCGCAATTAACCTGAAGCGAATAGACCGCTTGAGAAGCGTCCTGACCTTTCGTGTTGTTACGCTGATTg      <  2:239264/251‑1 (MQ=255)                                                                                                                                                                                                                                                gtCATCAAATTCAATATCATCAAAATTGATCTGCGGTTCGCCACGTTCAGTACGCTCAGAACGTTCAGAAGATTGACGAGAACGACCAGTACGACGATCGCTGGGATCGGGTTTCAGCTCGTCTTCTGGTTTGAATTCATCCa                                                                                                              >  2:266175/1‑143 (MQ=255)                                                                                                                                                                                                                                                gtCATCAAATTCAATATCATCAAAATTGATCTGCGGTTCGCCACGTTCAGTACGCTCAGAACGTTCAGAAGATTGACGAGAACGACCAGTACGACGATCGCTGGGATCGGGTTTCAGCTCGTCTTCTGGTTTGAATTCATCCa                                                                                                              <  1:266175/143‑1 (MQ=255)                                                                                                                                                                                                                                                gtCATCAAATTCAATATCATCAAAATTGATCTGCGGTTCGCCACGTTCAGTACGCTCAGAACGTTCAGAAGATTGACGAGAACGACCAGTACGACGATCGCTGGGATCGGGTTTCAGCTCGTCTTCTGGTTTGAATTCATCCATTTAACACCCCACTAAAAGGTTAATGCTTACCACGTTGCAATTAACCTGAAGCTAATAGACCGCTTGATAAGCGGCCTGACCTTTCTTGTTGTTACGCTGATTGACaa  <  1:66350/251‑1 (MQ=255)                                                                                                                                                                                                                                                        |                                                                                                                                                                                                                                                   GAGGGGGCTTTTAGCGCAGAACCGATACCGATGATCAACAGCAGTAGAACCAGAATGCCGACGCCCATCATCATATACTGACGAGAAGCGGGTTTACTGGCTGCTTTTTTGCGCTTACGCGGACGACGCTCTACGCGCTCTTCATCCACGGTTTCATCTTCGGATTCGTCAATTTCTTCTTCGATTTCCGGTTCCTCATTGCGCTCTTTTTGCGCACGAGTCGGACGGCGATCGTCAGTGTCATCAAGTTCAATATCATCAAAATTGATCTGCGGTTCGCCACGTTCAGTACGCTCAGAACGTTCAGAAGATTGACGAGAACGACCAGTACGACGATCGCTGGGATCGGGTTTCAGCTCGTCTTCTGGTTTGAATTCATCCATTTAACACCCCACTAAAAGGTTAATGCTTACCACGTTGCAATTAACCTGAAGCTAATAGACCGCTTGATAAGCGGCCTGACCTTTCTTGTTGTTACGCTGATTGACAA  >  NC_000913/3516924‑3517413

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