Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	NC_000913	310	(T)9→8	intergenic (+55/‑27)	thrL → / → thrA	thr operon leader peptide/fused aspartate kinase/homoserine dehydrogenase 1

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	310	0	T	.	100.0%	117.2 / NA	27	intergenic (+55/‑27)	thrL/thrA	thr operon leader peptide/fused aspartate kinase/homoserine dehydrogenase 1
Reads supporting (aligned to +/- strand):  ref base T (0/0);  new base . (12/15);  total (12/15)

GCAGCTTCTGAACTGGTTACCTGCCGTGAGTAAATTAAAATTTTATTGACTTAGGTCACTAAATACTTTAACCAATATAGGCATAGCGCACAGACAGATAAAAATTACAGAGTACACAACATCCATGAAACGCATTAGCACCACCATTACCACCACCATCACCATTACCACAGGTAACGGTGCGGGCTGACGCGTACAGGAAACACAGAAAAAAGCCCGCACCTGACAGTGCGGGCTTTTTTTTTCGACCAAAGGTAACGAGGTAACAACCATGCGAGTGTTGAAGTTCGGCGGTACATCAGTGGCAAATGCAGAACGTTTTCTGCGTGTTGCCGATATTCTGGAAAGCAATGCCAGGCAGGGGCAGGTGGCCACCGTCCTCTCTGCCCCCGCCAAAATCACCAACCACCTGGTGGCGATGATTGAAAAAACCATTAGCGGCCAGGATGCTTTACCCAATATCAGCGATGCCGAACGTATTTTTGCCGA  >  NC_000913/66‑554                                                                                                                                                                                                                                                     |                                                                                                                                                                                                                                                     gcagcTTCTTAACTGGTTACCTGCCGTGAGTAAATTAAAATTTTATTGACTTAGTTCACTAAATACTTTAACCAATATAGGCATAGCGCACAGACAGATAAAAATTACAGAGTACACAACATCCATGAAACGCATTAGCACCACCATTACCACCACCATCACCATTACCACAGGTAACGGTGCGGGCTGACGCGTACAGGAAACACAGAAAAAAGCCCGCACCTGACAGTGCGGGCTTTTTTTT‑CGACCaa                                                                                                                                                                                                                                               <  1:554643/251‑1 (MQ=255)       tCTGAACTGGTTACCTGCCGTGAGTAAATTAAAATTTTATTGACTTAGGTCACTAAATACTTTAACCAATATAGGCATAGCGCACAGACAGATAAAAATTACAGAGTACACAACATCCATGAAACGCATTAGCACCACCATTACCACCACCATCACCATTACCACAGGTAACGGTGCGGGCTGACGCGTACAGGAAACACAGAAAAAAGCCCGCACCTGACAGTGCGGGCTTTTTTTT‑CGACCAAAGGTa                                                                                                                                                                                                                                          <  1:538925/250‑1 (MQ=255)                        ccGTGAGTAAATTAAAATTTTATTGACTTAGGTCACTAAATACTTTAACCAATATAGGCATAGCGCACAGACAGATAAAAATTACAGAGTACACAACATCCATGAAACGCATTAGCACCACCATTACCACCACCATCACCATTACCACAGGTAACGGTGCGGGCTGACGCGTACAGGAAACACAGAAAAAAGCCCGCACCTGACAGTGCGGGCTTTTTTTT‑CGACCAAAGGTAACGAGGTAACAACCATGc                                                                                                                                                                                                                        >  1:165473/1‑251 (MQ=255)                             aGTAAATTAAAATTTTATTGACTTAGGTCACTAAATACTTTAACCAATATAGGCATAGCGCACAGACAGATAAAAATTACAGAGTACACAACATCCATGAAACGCATTAGCACCACCATTACCACCACCATCACCAGTACCACAGGTAACGGTGCGGGCTGACGCGTTCAGGAAACACAGAAAAAAGCCCGCACCTGACAGTGCGGGCTTGTTTTG‑CGACCAAAGGTAACGAGGTAACAACCATGCGAgtg                                                                                                                                                                                                                   <  2:394193/251‑1 (MQ=255)                                            tATTGACTTAGGTCACTAAATACTTTAACCAATATAGGCATAGCGCACAGACAGATAAAAATTACAGAGTACACAACATCCATGAAACGCATTAGCACCACCATTACCACCACCATCACCATTACCACAGGTAACGGTGCGGGCTGACGCGTACAGGAAACACAGAAAAAAGCCCGCACCTGACAGTGCGGGCTTTTTTTT‑CGACCAAAGGTAACGAGGTAACAACCATGCGAGTGTTGAAGTTcggcggg                                                                                                                                                                                                    >  2:610272/1‑250 (MQ=255)                                                  cTTAGGACACTAAATACTTTAACCAATATAGGCATAGAGCACAGACAGATAAACATTACAGAGTACACAACATCCATGAAACGCATTAGCACCACCATTACCACCACCATCACCATTACCACAGGTAACGGTGCGGGCTGACGCGTACAGGAAACACAGAAAAAAGCCCGCACCTGACAGTGCGGGCTTTTTTTT‑CGACCAAAGGTAACGAGGTAACAACCATGCGAGTATTGAAGTTAGGCGGTACATcc                                                                                                                                                                                              >  2:112946/1‑250 (MQ=255)                                                                          aaTATAGGCATAGCGCACAGACAGATAAAAATTACAGAGTACACAACATCCATGAAACGCATTAGCACCACCATTACCACCACCATCACCATTACCACAGGTAACGGTGCGGGCTGACGCGTACAGGAAACACAGAAAAAAGCCCGCACCTGACAGTGCGGGCTTTTTTTT‑CGACCAAAGGTAACGAGGTAACAACCATGCGAGTGTTGAAGTTCGGCGGTACATCAGGGGCAAATGCAGAACGTTTTCTg                                                                                                                                                                      >  1:553025/1‑251 (MQ=255)                                                                            tataGGCATAGCGCACAGACAGATAAAAATTACAGAGTACACAACATCCATGAAACGCATTAGCACCACCATTACCACCACCATCACCATTACCACAGGTAACGGTGCGGGCTGACGCGTACAGGAAACACAGAAAAAAGCCCGCACCTGACAGTGCGGGCTTTTTTTT‑CGACCAAAGGTAACGAGGTAACAACCATGCGAGTGTTGAAGTTCGGCGGTACATCAGTGGCAAATGCAGAACGTTTTCTGCg                                                                                                                                                                    >  1:474774/1‑251 (MQ=255)                                                                               aGGCATAGCGCACAGACAGATAAAAATTACAGAGTACACAACATCCATGACACGCATTAGCACCACCATTACCACCACCATCACCATTACCACAGGTAACGGTGCGGGCTGACGCGTACAGGAAACACAGAAAAAAGCCCGCACCTGACAGTGCGGGCTTTTTTTT‑CGACCAAAGGTAACGAGGTAACAACCATGCGAGTGTTGAAGTTCGGCGGTACATCAGTGGCAAATGCAGAACGTTTTCTGCgggt                                                                                                                                                                 >  2:110506/1‑248 (MQ=255)                                                                                            agacagaTAAAAATTACAGAGTACACAACATCCATGAAACGCATTAGCACCACCATTACCACCACCATCACCATTACCACAGGTAACGGTGCGGGCTGACGCGTACAGGAAACACAGAAAAAAGCCCGCACCTGACAGTGCGGGCTTTTTTTT‑CGACCAAAGGTAACGAGGTAACAACCATGCGAGTGTTGAAGTTCGGCGGTACATCAGTGGCAAATGCAGAACGTTTTCTGCGGGTTGCCGATATTCTg                                                                                                                                                    <  1:469334/251‑1 (MQ=255)                                                                                                 gaTAAAAATTACAGAGTACACAACATCCATGAAACGCATTAGCACCACCATTACCACCACCATCACCATTACCACAGGTAACGGTGCGGGCTGACGCGTACAGGAAAGACAGAAAAAAGCCCGCACCTGACAGTGCGGGCTTTTTTTT‑CGACCAAAGGTAACGAGGTAACAACCATGCGAGTGTTGAAGTTCGGCGGTACATCAGTGGCAAATGCAGAACGTTTTCGGCGGGTTGCCGATATTCTGGAAAg                                                                                                                                               <  2:207376/251‑1 (MQ=255)                                                                                                       aaTTACAGAGTACACAACATCCATGAAACGCATTAGCACCACCATTACCACCACCATCACCATTACCACAGGTAACGGTGCGGGCTGACGCGTACAGGAAACACAGAAAAAAGCCCGCACCTGACAGTGCGGGCTTTTTTTT‑CGACCAAAGGTAACGAGGTAACAACCATGCGAGTGTTGAAGTTCGGCGGTACATCAGTGGCAAATGCAGAACGTTTTCTGCGGGTTGCCGATATTCTGGAAAGCAATGc                                                                                                                                         <  1:262852/251‑1 (MQ=255)                                                                                                             agagTCCACAACATCCATGAACCGCATTAGCACCACCATTACCACCACCATCACCATTCCCACAGGTAACGGTGCGGGCTGACGCGTACAGGAAACACAGAAAAAAGCCCGCACCTGACAGTGCGGGCTTTTTTTT‑CGACCAAAGGTAACGAGGTAACAACCATGCGAGTGTTGAAGTTCGGCGGTACATCAGTGGCAAATGCAGAACGTTTTCTGCGGGTTGCCGATATTCTGGAAAGCAATGCcaggca                                                                                                                                   <  2:446625/251‑1 (MQ=255)                                                                                                                                     cATTAGCACCACCCTTACCACCACCATCACCATTACCCCAGGTCACGGTGCGGGCTGACGCGTACCGGACACACAGAAAAAAGCCCGCAACTGACAGTGCGGGCTTTTTTTT‑CGACCAAAGGTAACGAGGGAACAACAATGCGAGTGTTGAAGTTCGGCGGTACATCAGTGGCAAATGCAGAACGTTTTCTGCGGGTTGCCGATATTCTGGAAAGCAATGCCAGGGAGGGGCCGGGGGCCACCGTCccctc                                                                                                           >  2:378655/1‑251 (MQ=255)                                                                                                                                     cATTAGCACCACCATTACCACCACCATCACCATTACCACAGGTAACGGTGCGGGCTGACGCGTACAGGAAACACAGAAAAAAGCCCGCACCTGACAGTGCGGGCTTTTTTTT‑CGACCAAAGGTAACGAGGTAACAACCATGCGAGTGTTGAAGTTCGGCGGTACATCAGTGGCAAATGCAGAACGTTTTCTGCGGGTTGCCGATATTCTGGAAAGCAATGCCAGGCAGGGGCAGGTGGCCACCGTCctctc                                                                                                           <  1:373958/251‑1 (MQ=255)                                                                                                                                           caccacCATTCCCACCACCATCACCATTACAACAGGTAACGGTGCGGGCTGACGCGTACAGGAAACACAGAAAAAAGCCCGCACCTGACAGTGCGGGCTTTTTTTT‑CGACCAAAGGTAACGAGGTAACAACCATGCGAGTGTTGAAGTTCGGCGGTACATCAGTGGCAAATGCAGAACGTTTTCTGCGGGTTGCCGATATTCTGGAAAGCAATGCCAGGCAGGGGCAGGTGGCCACCGTCCTCTCTGcccc                                                                                                     >  1:192194/1‑251 (MQ=255)                                                                                                                                            cccaccaTTACCTCCACCATCACCAGTACCACAAGTAACGGTGCGGGCTGACGCGTACAGGAAACTCAGAAAAAAGCCCGCACCTGACAGTGCGGGCTTTTTTTT‑CGACCAAAGGTAACGAGGTAACAACCAGGCGAGTGTTGAAGTTCGGCGGTACATCAGTGGCAAATGCAGAACGTTTTCTGCGGGTTGCCGATAGTCTTGAAAGCAATGCCAGGCAGGGGCAGGTGGCCACCGTCCTCTCTGccccc                                                                                                    <  1:278939/250‑1 (MQ=255)                                                                                                                                                      ccaccaccaTCACCATTACCACAGGTAACGGTGCGGGCTGACGCGTACAGGAAACACAGAAAAAAGCCCGCACCTGACAGTGCGGGCTTTTTTTT‑CGACCAAAGGTAACGAGGTAACAACCATGCGAGTGTTGAAGTTCGGCGGTACATCAGTGGCAAATGCAGAACGTTTTCTGCGGGTTGCCGATATTCTGGAAAGCAATGCCAGGCAGGGGCAGGTGGCCACCGTCCTCTCTGCCCCCGCCAAAATc                                                                                           >  1:595157/1‑250 (MQ=255)                                                                                                                                                              atcaccatTACCACAGGTAACGGTGCGGGCTGACGCGTACAGGAAACACAGAAAAAAGCCCGCACCTGACAGTGCGGGCTTTTTTTT‑CGACCAAAGGTAACGAGGTAACAACCATGCGAGTGTTGAAGTTCGGCGGTACATCAGTGGCAAATGCAGAACGTTTTCTGCGGGTTGCCGATATTCTGGAAAGCAATGCCAGGCAGGGGCAGGTGGCCACCGTCCTCTCTGCCCCCGCCAAAATCACCAaccac                                                                                  >  1:487806/1‑251 (MQ=255)                                                                                                                                                              atcaccatTACCACAGGTAACGGGGCGGGCTAACGCGGACAGGAACCACAGAAAAAAGCCCGCGCCTGACAGTGCGGGCTTTTTTTT‑CGACCAAAGGTAACGAGGTAACAACCATGCGAGTGTTGAAGTTCGGCGGTACATCAGTGGCAAATGCAGAACGTTTTCTGCGGGTTGCCGATATTCTGGAAAGCAATGCCAGGCAGGGGCAGGTGGCCACCTTCCTCTCTGCCCCCGCCAAAATCACCAaccac                                                                                  <  1:378655/251‑1 (MQ=255)                                                                                                                                                                caccctTACCACAGGAAACGGTGCGGGCTGACGCGTCCAGGAAACCCAGAAAAAAGCCCGCACCTGACAGTGCGTGCTTTTTTTA‑CGACCAAATGTATCGATGTAACAACCATGCGAGTGTTGAAGTTCGGCGGTACATCAGTGGCAACTTCAGAACGTTTTCTGCGGGTTGCCGATATTCTGTACAGCAATCCCAGTCAGGGGCAGGTGTCAACGGTCCTCTCTGCCCCCGCCAAAATCgccatccatct                                                                                <  2:141822/251‑4 (MQ=255)                                                                                                                                                                       aCCACAGGTAACGGTGCGGGCTGACGCGTACAGGAAACACAGAAAAAAGCCCGCACCTGACAGTGCGGGCTTTTTTTT‑CGACCAAAGGTAACGAGGTAACAACCATGCGAGTGTTGAAGTTCGGCGGTACATCAGTGGCAAATGCAGAACGTTTTCTGCGGGTTGCCGATATTCTGGAAAGCAATGCCAGGCAGGGGCAGGTGGCCACCGTCCTCTCTGCCCCCGCCAAAATCACCAACCACCTGGTGGc                                                                          <  2:165473/250‑1 (MQ=255)                                                                                                                                                                                        gggCTGACGCGTACAGGAAACACAGAAAAAAGCCCGCACCTGACAGTGCGGGCTTTTTTTT‑CGACCAAAGGTAACGAGGTAACAACCATGCGAGTGTTGAAGTTCGGCGGTACATCAGTGGCAAATGCAGAACGTTTTCTGCGGGTTGCCGATATTCTGGAAAGCAATGCCAGGCAGGGGCAGGTGGCCACCGTCCTCTCTGCCCCCGCCAAAATCACCAACCACCTGGTGGCGATGATTGAAAAAACCAt                                                        >  1:109114/1‑251 (MQ=255)                                                                                                                                                                                                gcgTACAGGAAACACAGAAAAAAGCCCGCACCTGACAGGGCGGGCTTTTTTTT‑CGACCAAAGGTAACGAGGTAACAACCATGCGAGTGTTGAAGTTCGGCGGTACATCAGTGGCAAATGCAGAACGTTTTCTGCGGGTTGCCGATATGCTGGAAAGCAATGCCAGGCAGGGGCAGGTGGCCACCGTCCTCTCTGCCCCCGCCAAAATCACCAACCACCTGGTGGCGATGATTGAAAAAACCATTAGCGGcc                                                <  1:610272/251‑1 (MQ=255)                                                                                                                                                                                                    aCAGGAAAAACAGAAAAAAGCCCGCACCTGACAGTGCGGGCTTTTTTTT‑CGACCAAAGGTAACGAGGTACCACCCATGCGAGTGTTGAAGTTCGGCGGTACATCAGTGGCAAATGCAGAACGTTTTCTGCGGGTTGCCGATATTCTGGAAAGCAATGCCAGGCAGGGGCAGGTGGCCACCGTCCTCTCTGCCCCCGCCAAAATCACCAACCACCTGGTGGCGATGATTGAAAAAACCATTAGCGGCCAgg                                             <  1:157336/250‑1 (MQ=255)                                                                                                                                                                                                         aaaCACAGAAAAAAGCCCGCCCCTGACAGTGCGGGCTTTTTTTT‑CGACCAAAGGTAACGAGGTAACAACCATGCGAGTGTTGAAGTTCGGCGGTACATCAGTGGCAAATGCAGAACGTTTTCTGCGGGTTGCCGATATTCTGGAAAGCAATGCCAGGCAGGGGCAGGTGGCCACCGTCCTCTCTGCCCCCGCCAAAATCACCAACCACCTGGTGGCGATGATTGAAAAAACCATTAGCGGCCAGGATGCtt                                       <  2:108034/251‑1 (MQ=255)                                                                                                                                                                                                                gAAAAAAGCCCGCACCTGACAGTGCGGGCTTTTTTTT‑CGACCAAAGGTAACGAGGTAACAACCATGCGAGTGTTGAAGTTCGGCGGAACATCAGTGGCAAATGCAGAACGTTTTCTGCGGGTTGCCGATATTCTGGAAAGCAATGCCAGGCAGGGGCAGGTGGCCACCGTCCTATCTGCCCCCGCCAAAATCACCAACCACCTGGTAGCGATGATTGGAAAAACCATTCGCGGCCAGGATGCTTTTCCCa                                 >  2:166635/1‑250 (MQ=255)                                                                                                                                                                                                                                                ttttttCGACCAAAGGTAAAGAGGTAACAACCATGCGAGTGTTGTAGTTCGGCGGTACATCAGTGGCAAATGCAGAACGTTTTCTGCGGGTTGCCGATATTCTGGAAAGCAATGCCAGGCAGGGGCAGGTGGCCACCGTCCTCTCTGCCCCCGCCAAAATCACCAACCACCTGGTGGCGATGATTGTAAAAACCATTAGCGGCCAGGATGCTTTACCCAATATCAGCGATGCAGAACTAATTTTTGCCGa  >  2:143448/1‑250 (MQ=255)                                                                                                                                                                                                                                                     |                                                                                                                                                                                                                                                     GCAGCTTCTGAACTGGTTACCTGCCGTGAGTAAATTAAAATTTTATTGACTTAGGTCACTAAATACTTTAACCAATATAGGCATAGCGCACAGACAGATAAAAATTACAGAGTACACAACATCCATGAAACGCATTAGCACCACCATTACCACCACCATCACCATTACCACAGGTAACGGTGCGGGCTGACGCGTACAGGAAACACAGAAAAAAGCCCGCACCTGACAGTGCGGGCTTTTTTTTTCGACCAAAGGTAACGAGGTAACAACCATGCGAGTGTTGAAGTTCGGCGGTACATCAGTGGCAAATGCAGAACGTTTTCTGCGTGTTGCCGATATTCTGGAAAGCAATGCCAGGCAGGGGCAGGTGGCCACCGTCCTCTCTGCCCCCGCCAAAATCACCAACCACCTGGTGGCGATGATTGAAAAAACCATTAGCGGCCAGGATGCTTTACCCAATATCAGCGATGCCGAACGTATTTTTGCCGA  >  NC_000913/66‑554

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