Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	NC_000913	697,151	A→G	intergenic (‑18/+362)	metT ← / ← asnB	tRNA‑Met/asparagine synthetase B

Read alignment evidence...
	seq id	position		change	freq	score	reads	annotation	genes	product
*	NC_000913	697,151	0	A→G	100.0%	126.3	35	intergenic (‑18/+362)	metT/asnB	tRNA‑Met/asparagine synthetase B
Reads supporting (aligned to +/- strand):  new base (19/16):  ref base (0/0):  total (19/16)

TCTGGTGAATGGTGCGGGAGGCGAGACTTGAACTCGCACACCTTGCGGCGCCAGAACCTAAATCTGGTGCGTCTACCAATTTCGCCACTCCCGCAAAAAAAGATGGTGGCTACGACGGGATTCGAACCTGTGACCCCATCATTATGAGTGATGTGCTCTAACCAACTGAGCTACGTAGCCATCTTTTTTTTCGCGATACCTTATCGGCGTTGCGGGGCGCATTATGCGTATAGAGCCTTGCAGCGTCAACCTCTTTTTCAAGGAAAATTGCTCGAAAGTGACTGTTTGGTTAGGTTGCGAACAGCGTGGCGCTATATTCGTCAATTATTGTTTACTTTGTGTTTTTTCCCACCCTACAGCCATTCTTTTGTCATACAGGATGAAATTCGGAATTTAACAATA  >  NC_000913/696954‑697355                                                                                                                                                                                                      |                                                                                                                                                                                                             tCTGGTGAATGGTGCGGGAGGCGAGACTTGAACTCGCACACCTTGCGGCGCCAGAACCTAAATCTGGTGCGTCTACCAATTTCGCCACTCCCGCAAAAAAAGATGGTGGCTACGACGGGATTCGAACCTGTGACCCCATCATTATGAGTGATGTGCTCTAACCAACTGAGCTACGTAGCCATCTTTTTTTTCGCGATGCCTTATCGGCGt                                                                                                                                                                                                  <  2:510855/210‑1 (MQ=255)   ggctGAATTGTGCCGGAAGCGAGAGCGGAACTCGCACCCCTTGCGGGGCCCGAACCTAAATCATGTGCGTGTAGCAATGGCGCCACTCCCGCAAAAAAAGATGGTGGCTACGACGGGATTCGAACCTGTGACCCCATCATTATGAGTGATGTGCTCTAACCAACTGAGCTACGTAGCCATCTTTTTTTTCGCGATGCCTTATCGGCGTTg                                                                                                                                                                                                <  2:45730/207‑1 (MQ=25)             tGCGGGAGGCGAGACTTGAACTCGCACACCTTGCGGCGCCAGAACCTAAATCTGGTGCGTCTACCAATTTCGCCACTCCCGCAAAAAAAGATGGTGGCTACGACGGGATTCGAACCTGTGACCCCATCATTATGAGTGATGTGCTCTAACCAACTGAGCTACGTAGCCATCTTTTTTTTCGCGATGCCTTATCGGCGTTGCGGGGCGCAt                                                                                                                                                                                      >  1:533495/1‑210 (MQ=255)               cGGGAGGCGAGACTTGAACTCGCACACCTTGCGGCGCCAGAACCTAAATCTGGTGCGTCTACCAATTTCGCCACTCCCGCAAAAAAAGATGGTGGCTACGACGGGATTCGAACCTGTGACCCCATCATTATGAGTGATGTGCTCTAACCAACTGAGCTACGTAGCCATCTTTTTTTTCGCGATGCCTTATCGGCGTTGCGGGGCGCATTa                                                                                                                                                                                    >  2:362125/1‑210 (MQ=255)                 ggAGGCGAGACTTGAACTCGCACACCTTGCGGCGCCAGAACCTAAATCTGGTGCGTCTACCAATTTCGCCACTCCCGCAAAAAAAGATGGTGGCTACGACGGGATTCGAACCTGTGACCCCATCATTATGAGTGATGTGCTCTAACCAACTGAGCTACGTAGCCATCTTTTTTTTCGCGATGCCTTATCGGCGTTGCGGGGCGCATTATg                                                                                                                                                                                  <  1:397593/210‑1 (MQ=255)                     gCGAGACTTGAACTCGCACACCTTGCGGCGCCAGAACCTAAATCTGGTGCGTCTACCAATTTCGCCACTCCCGCAAAAAAAGATGGTGGCTACGACGGGATTCGAACCTGTGACCCCATCATTATGAGTGATGTGCTCTAACCAACTGAGCTACGTAGCCATCTTTTTTTTCGCGATGCCTTATCGGCGTTGCGGGGCGCATTATGCGta                                                                                                                                                                              <  2:136583/210‑1 (MQ=255)                      cGAGACTTGAACTCGCACACCTTGCGGCGCCAGAACCTAAATCTGGTGCGTCTACCAATTTCGCCACTCCCGCAAAAAAAGATGGTGGCTACGACGGGATTCGAACCTGTGACCCCATCATTATGAGTGATGTGCTCTAACCAACTGAGCTACGTAGCCATCTTTTTTTTCGCGATGCCTTATCGGCGTTGCGGGGCGCATTATGCGtat                                                                                                                                                                             >  2:444809/1‑210 (MQ=255)                      cGAGACTTGAACTCGCACACCTTGCGGCGCCAGAACCTAAATCTGGTGCGTCTACCAATTTCGCCACTCCCGCAAAAAAAGATGGTGGCTACGACGGGATTCGAACCTGTGACCCCATCATTATGAGTGATGTGCTCTAACCAACTGAGCTACGTAGCCATCTTTTTTTTCGCGATGCCTTATCGGCGTTGCGGGGCGCATTATGCGtat                                                                                                                                                                             >  1:285417/1‑210 (MQ=255)                         gaCTTGAACTCGCACACCTTGCGGCGCCAGAACCTAAATCTGGTGCGTCTACCAATTTCGCCACTCCCGCAAAAAAAGATGGTGGCTACGACGGGATTCGAACCTGTGACCCCATCATTATGAGTGATGTGCTCTAACCAACTGAGCTACGTAGCCATCTTTTTTTTCGCGATGCCTTATCGGCGTTGCGGGGCGCATTATGCGTATaga                                                                                                                                                                          >  1:155781/1‑210 (MQ=255)                          aCTTGAACTCGCACACCTTGCGGCGCCAGAACCTAAATCTGGTGCGTCTACCAATTTCGCCACTCCCGCAAAAAAAGATGGTGGCTACGACGGGATTCGAACCTGTGACCCCATCATTATGAGTGATGTGCTCTAACCAACTGAGCTACGTAGCCATCTTTTTTTTCGCGATGCCTTATCGGCGTTGCGGGGCGCATTATGCGTATagag                                                                                                                                                                         <  1:444827/210‑1 (MQ=255)                          aCTTGAACTCGCAAACCTTGCGGCGCCAGAACCTAAATCTGGTGCGTCTACCAAATTCGCCACCCCCGCAAAAAAAGATGGTTGCTACGACGGGATTCGAACCTGTGACCCCATCATTATGAGTGATGTGCTCTAACCAACTGAGCTACGTAGCCATCTTTTTTTTCGCGATGCCTTATCGGCGTTGCGGGGCGCATTATGCGTATagag                                                                                                                                                                         <  1:267302/210‑1 (MQ=255)                           cTTGAACTCGCACACCTTGCGGCGCCAGAACCTAAATCTGGTGCGTCTACCAATTTCGCCACTCCCGCAAAAAAAGATGGTGGCTACGACGGGATTCGAACCTATGACCCCATCATTATGAGTGATGTGCTCTAACCAACTGAGCTACGTAGCCATCTTTTTTTTCGCGATGCCTTACCGGCGTTGCGTGGCGCATTATCCGTATTGAGc                                                                                                                                                                        >  1:14564/1‑210 (MQ=255)                                 cTCGCACACCTTGCGGCGCCAGAACCTAAATCTGGTGCGTCTACCAATTTCGCCACTCCCGCAAAAAAAGATGGTGGCTACGACGGGATTCGAACCTGTGACCCCATCATTATGAGTGATGTGCTCTAACCAACTGAGCTACGTAGCCATCTTTTTTTTCGCGATGCCTTATCGGCGTTGCGGGGCGCATTATGCGTATAGAGCCTTGCa                                                                                                                                                                  >  1:503659/1‑210 (MQ=255)                                         ccTTGCGGCGCCAGAACCTAAATCTGGTGCGTCTACCAATTTCGCCACTCCCGCAAAAAAAGATGGTGGCTACGACGGGATTCGAACCTGTGACCCCATCATTATGAGTGATGTGCTCTAACCAACTGAGCTACGTAGCCATCTTTTTTTTCGCGATGCCTTATCGGCGTTGCGGGGCGCATTATGCGTATAGAGCCTTGCAGCGTCAAc                                                                                                                                                          >  1:348229/1‑210 (MQ=255)                                                gcggCAGAACCTAAATCTGGTGGGTCTACCAATTTCGCCACTCCCGCAAAAAAAGATGGTGGCTACGACGGGATTCGAACCTGTGACCCCATCATTATGAGTGATGTGCTCTAACCAACTGAGCTACGTAGCCATCTTTTTTTTCGCGATGCCTTATCGGCGTTGCGGGGCGCATTATGCGTATAGAGCCTTGCAGCGTCAACCTCtttt                                                                                                                                                   <  1:333032/210‑1 (MQ=255)                                                     aGAACCTAAATCTGGTGCGTCTACCAATTTCGCCACTCCCGCAAAAAAAGATGGTGGCTACGACGGGATTCGAACCTGTGACCCCATCATTATGAGTGATGTGCTCTAACCAACTGAGCTACGTAGCCATCTTTTTTTTCGCGATGCCTTATCGGCGTTGCGGGGCGCATTATGCGTATAGAGCCTTGCAGCGTCAACCTCTTTTTCAAg                                                                                                                                              <  1:415745/210‑1 (MQ=255)                                                                    tGCGTCTACCAATTTCGCCACTCCCGCAAAAAAAGATGGTGGCTACGACGGGATTCGAACCTGTGACCCCATCATTATGAGTGATGTGCTCTAACCAACTGAGCTACGTAGCCATCTTTTTTTTCGCGATGCCTTATCGGCGTTGCGGGGCGCATTATGCGTATAGAGCCTTGCAGCGTCAACCTCTTTTTCAAGGAAAATTGCTCGaaa                                                                                                                               >  1:18862/1‑210 (MQ=255)                                                                          tACCAATTTCGCCACTCCCGCAAAAAAAGATGGTGGCTACGACGGGATTCGAACCTGTGACCCCATCATTATGAGTGATGTGCTCTAACCAACTGAGCTACGTAGCCATCTTTTTTTTCGCGATGCCTTATCGGCGTTGCGGGGCGCATTATGCGTATAGAGCCTTGCAGCGTCAACCTCTTTTTCAAGGAAAATTGCTCGAAAGTGACt                                                                                                                         >  2:445690/1‑210 (MQ=255)                                                                                tttCGCCACTCCCGAAAAAAAAGATGGTGGCTACGACGGGATTCGAACCTGTGACCCCATCATTATGAGTGATGTGCTCTAACCAACTGAGCTACGTAGCCATCTTTTTTTTCGCGATGCCTTATCGGCGTTGCGGGGCGCATTATGCGTATAGAGCCTTGCAGCGTCAACCTCTTTTTCAAGGAAAATTGCTCGAAAGTGACTGTTTgg                                                                                                                   <  1:314253/210‑1 (MQ=255)                                                                                   cGCCACTCCCGCAAAAAAAGATGGTGGCTACGACGGGATTCGAACCTGTGACCCCATCATTATGAGTGATGTGCCCTAACCAACTGAGCTACGTAGACATCTTTTTTTTCGCGATGCCTTATCGGCGTTGCGGGGCGCATTATGCGTATAGAGCCTTGCACCGTCAACCTCTTTTTCAAGCACAATTGCGCGAAAGTGACTGTTTGGttc                                                                                                                >  2:362826/1‑209 (MQ=38)                                                                                           ccGCAAAAAAAGATGGTGGCTACGACGGGATTCGAACCTGTGACCCCATCATTATGAGTGATGTGCTCTAACCAACTGAGCTACGTAGCCATCTTTTTTTTCGCGATGCCTTATCGGCGTTGCGGGGCGCATTATGCGTATAGAGCCTTGCAGCGTCAACCTCTTTTTCAAGGAAAATGGCTCGAAATTGACTGTTGGGTTGGGTTCCGa                                                                                                        >  1:17216/1‑210 (MQ=255)                                                                                                      gATGGTGGCTACGACGGGATTCGAACCTGTGACCCCATCATTATGAGTGATGTGCTCTAACCAACTGAGCTACGTAGCCATCTTTTTTTTCGCGATGCCTTATCGGCGTTGCGGGGCGCATTATGCGTATAGAGCCTTGCAGCGTCAACCTCTTTTTCAAGGAAAATTGCTCGAAAGTGACTGTTTGGTTAGGTTGCGAACAGCGTGgcg                                                                                             <  1:278628/210‑1 (MQ=255)                                                                                                       aTGGTGGCTACGACGGGATTCGAACCTGTGACCCCATCATTATGAGTGATGTGCTCTAACCAACTGAGCTACGTAGCCATCTTTTTTTTCGCGATGCCTTATCGGCGTTGCGGGGCGCATTATGCGTATAGAGCCTTGCAGCGTCAACCTCTTTTTCAAGGAAAATTGCTCGAAATTGACTGTTTGGTTAGGTTGCGAACAGCGTGgcgc                                                                                            >  1:287742/1‑210 (MQ=255)                                                                                                               tACGACGGGATTCGAACCTGTGACCCCATCATTATGAGTGATGTGCTCTAACCAACTGAGCTACGTAGCCATCTTTTTTTTCGCGATGCCTTATCGGCGTTGCGGGGCGCATTATGCGTATAGAGCCTTGCAGCGTCAACCTCTTTTTCAAGGAAAATTGCTCGAAAGTGACTGTTTGGTTAGGTTGCGAACAGCGTGGCGCTATATTCg                                                                                    >  1:49843/1‑210 (MQ=255)                                                                                                                 cgacgGGATTCGAACCTGTGACCCCATCATTATGAGTGATGTGCTCTAACCAACTGAGCTACGTAGCCATCTTTTTTTTCGCGATGCCTTATCGGCGTTGCGGGGCGCATTATGCGTATAGAGCCTTGCAGCGTCAACCTCTTTTTCAAGGAAAATTGCTCGAAAGTGACTGTTTGGTTAGGTTGCGAACAGCGTGGCGCTATATTCGTc                                                                                  <  1:445708/210‑1 (MQ=255)                                                                                                                       gATTCGAACCTGTGACCCCATCATTATGAGTGATGTGCTCTAACCAACTGAGCTACGTAGCCATCTTTTTTTTCGCGATGCCTTATCGGCGTTGCGGGGCGCATTATGCGTATAGAGCCTTGCAGCGTCAACCTCTTTTTCAAGGAAAATTGCTCGAAAGTGACTGTTTGGTTAGGTTGCGAACAGCGTGGCGCTATATTCGTCAattat                                                                            >  2:397006/1‑210 (MQ=255)                                                                                                                                     aCCCCATCATTATGAGTGATGTGCTCTAACCAACTGAGCTACGTAGCCATCTTTTTTTTCGCGATGCCTTATCGGCGTTGCGGGGCGCATTATGCGTATAGAGCCTTGCAGCGTCAACCTCTTTTTCAAGGAAAATTGCTCGAAAGTGACTGTTTGGTTAGGTTGCGAACAGCGTGGCGCTATATTCGTAAATTATTGTTTACTTTgtgt                                                              >  2:350734/1‑210 (MQ=255)                                                                                                                                      ccccATCATTATGAGTGATGTGCTCTAACCAACTGAGCTACGTAGCCATCTTTTTTTTCGCGATGCCTTATCGGCGTTGCGGGGCGCATTATGCGTATAGAGCCTTGCAGCGTCAACCTCTTTTTCAAGGAAAATTGCTCGAAAGTGACTGTTTGGTTAGGTTGCGAACAGCGTGGCGCTATATTCGTCAATTATTGTTTACTTTGTGtt                                                             >  2:372332/1‑210 (MQ=255)                                                                                                                                           tcatTATGAGTGATGTGCCCTAACCAACTTAGCTACGTAGCCATCTTTTTTTTCGCGATGCCTTATCGGCGTTGCGGGGCGCATTATGCGTATAGAGCCTTGCAGCGTCAACCTCTTTTTCAAGGAAAATTGCTCGAAAGTGACTGTTTGGTTAGGTTGCGAACAGCGTGGCGCTATATTCGTCAATTATTGTTTACTTTGTGTTTTTTc                                                        <  1:349995/210‑1 (MQ=255)                                                                                                                                                aTGAGTGATGTGCTCTAACCAACTGAGCTACGTAGCCATCTTTTTTTTCGCGATGCCTTATCGGCGTTGCGGGGCGCATTATGCGTATAGAGCCTTGCAGCGTCAACCTCTTTTTCAAGGAAAATTGCTCGAAAGTGACTGTTTGGTTAGGTTGCGAACAGCGTGGCGCTATATTCGTCAATTATTGTTTACTTTGTGTTTTTTCCCAcc                                                   <  2:364296/210‑1 (MQ=255)                                                                                                                                                                   cAACTGAGCTACGTAGCCATCTTTTGTTTCGCGATGCCTTATCGGCGTTGCGGGGCGCATTATGCGTATAGAGACTTGCAGCGTCAACCTCTTTTTCAAGGAAAATTGCTCGAAAGTGACTGTTTGGTTAGGTTGCGAACAGCGTGGCGCTATATTCGTCAATTATTGTTTACTTTGTGTTTTTTCCCACCCTACAGCCATTCTTTTGTc                                <  2:392620/210‑1 (MQ=255)                                                                                                                                                                    aaCTGAGCTACGTAGCCATCTTTTTTTTCGCGATGCCTTATCGGCGTTGCGGGGCGCATTATGCGTATAGAGCCTTGCAGCGTCAACCTCTTTTTCAAGGAAAATTGCTCGAAAGTGACTGTTTGGTTAGGTTGCGAACAGCGTGGCGCTATATTCGTCAATTATTGTTTACTTTGTGTTTTTTCCCACCCTACAGCCATTCTTTTGTCa                               <  2:77810/210‑1 (MQ=255)                                                                                                                                                                        gAGCTACGTAGCCATCTTTTTTTTCGCGATGCCTTATCGGCGTTGCGGGGCGCATTATGCGTATAGAGCCTTGCAGCGTCAACCTCTTTTTCAAGGAAAATTGCTCGAAAGTGACTGTTTGGTTAGGTTGCGAACAGCGTGGCGCTATATTCGTCAATTATTGTTTACTTTGTGTTTTTTCCCACCCTACAGCCATTCTTTTGTCATACa                           <  1:293095/210‑1 (MQ=255)                                                                                                                                                                                       cTTTTTTTTCGCGATGCCTTATCGGCGTTGCGGGGCGCATTATGCGTATAGAGCCTTGCAGCGTCAACCTCTTTTTCAAGGAAAATTGCTCGAAAGTGACTGTTTGGTTAGGTTGCGAACAGCGTGGCGCTATATTCCTCAATTACTGTTTACTTTGTGTTTTTTCCCACCCTACAGCCCTTCTTTCGTCATACCGTATGAAATTCGGaa            >  2:288132/1‑210 (MQ=255)                                                                                                                                                                                                 gcgATGCCTTATCGGCGTTGCGGGGCGCATTATGCGTATAGAGTCTTGAAGCGTCAACCTCTTTTTAAAGGAAAATTGCTCGAAAGTGACTGCTTGGTTAGGTTGCGAACAGCGTGGCGCTATATTCGTCAATTATTGTTTACTTTGTGTTTTTTCCCATCCTACAGCCATTCTTTTGTCATACAGGATGAAATTCGGAATTTAACAATa  >  1:102758/1‑210 (MQ=255)                                                                                                                                                                                                      |                                                                                                                                                                                                             TCTGGTGAATGGTGCGGGAGGCGAGACTTGAACTCGCACACCTTGCGGCGCCAGAACCTAAATCTGGTGCGTCTACCAATTTCGCCACTCCCGCAAAAAAAGATGGTGGCTACGACGGGATTCGAACCTGTGACCCCATCATTATGAGTGATGTGCTCTAACCAACTGAGCTACGTAGCCATCTTTTTTTTCGCGATACCTTATCGGCGTTGCGGGGCGCATTATGCGTATAGAGCCTTGCAGCGTCAACCTCTTTTTCAAGGAAAATTGCTCGAAAGTGACTGTTTGGTTAGGTTGCGAACAGCGTGGCGCTATATTCGTCAATTATTGTTTACTTTGTGTTTTTTCCCACCCTACAGCCATTCTTTTGTCATACAGGATGAAATTCGGAATTTAACAATA  >  NC_000913/696954‑697355

Alignment Legend
Aligned base mismatch/match (shaded by quality score): ATCG/ATCG < 3 ≤ ATCG/ATCG < 10 ≤ ATCG/ATCG < 24 ≤ ATCG/ATCG < 37 ≤ ATCG/ATCG < 38 ≤ ATCG/ATCG
Unaligned base: atcg    Masked matching base: atcg    Alignment gap: ‑    Deleted base: ‑
Reads not counted as support for junction
read_name Not counted due to insufficient overlap past the breakpoint.
read_name Not counted due to not crossing MOB target site duplication.