Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	NC_000913	4,267,869	C→T	Y252Y (TAC→TAT)	tyrB →	tyrosine aminotransferase

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	4,267,869	0	C	T	100.0%	63.3 / NA	19	Y252Y (TAC→TAT)	tyrB	tyrosine aminotransferase
Reads supporting (aligned to +/- strand):  ref base C (0/0);  new base T (8/11);  total (8/11)

GATATTGCCTATCAAGGATTTGGTGCCGGTATGGAAGAGGATGCCTACGCTATTCGCGCCATTGCCAGCGCTGGATTACCCGCTCTGGTGAGCAATTCGTTCTCGAAAATTTTCTCCCTTTACGGCGAGCGCGTCGGCGGACTTTCTGTTATGTGTGAAGATGCCGAAGCCGCTGGCCGCGTACTGGGGCAATTGAAAGCAACAGTTCGCCGCAACTACTCCAGCCCGCCGAATTTTGGTGCGCAGGTGGTG  >  NC_000913/4267747‑4267998                                                                                                                           |                                                                                                                                  gATATTGCCTATCAAGGATTTGGTGCCGGTATGGAAGAGGATGCCTACGCTATTCGCGCCATTGCCAGCGCTGGATTACCCGCTCTGGTGAGCAATTCGTTCTCGAAAATTTTCTCCCTTTATGGCGAGCGCGTcggcgg                                                                                                                  >  2:581402/1‑140 (MQ=255) gATATTGCCTATCAAGGATTTGGTGCCGGTATGGAAGAGGATGCCTACGCTATTCGCGCCATTGCCAGCGCTGGATTACCCGCTCTGGTGAGCAATTCGTTCTCGAAAATTTTCTCACTTTATGGCGAGCGCGTcggcgg                                                                                                                  >  1:391184/1‑140 (MQ=255) gATATTGCCTATCAAGGATTTGGTGCCGGTATGGAAGAGGATGCCTACGCTATTCGCGCCATTGCCAGCGCTGGATTACCCGCTCTGGTGAGCAATTCGTTCTCGAAAATTTTCTCACTTTATGGCGAGCGCGTcggcgg                                                                                                                  <  2:391184/140‑1 (MQ=255)     ttGCCTATCAAGGATTTGGTGCCGGTATGGAAGAGGATGCCTACGCTATTCGCGCCATTGTCAGCGCTGGATTACCCGCTCTGGTGAGCAATTCGTTCTCGAAAATTTTCTCCCTTTATGGCGAGCGCGTCGGCGGACtt                                                                                                              <  2:514539/140‑1 (MQ=255)     ttGCCTATCAAGGATTTGGTGCCGGTATGGAAGAGGATGCCTACGCTATTCGCGCCATTGCCAGCGCTGGATTACCCGCTCTGGTGAGCAATTCGTTCTCGAAAATTTTCTCCCTTTATGGCGAGCGCGTCGGCGGACtt                                                                                                              <  1:131694/140‑1 (MQ=255)         cTATCAAGGATTTGGTGCCGGTATGGAAGAGGATGCCTACGCTATTCGCGCCATTGCCAGCGCTGGATTACCCGCTCTGGTGAGCAATTCGTTCTCGAAAATTTTCTCCCTTTATGGCGAGCGCGTCGGCGGACTTTCTg                                                                                                          <  1:160785/140‑1 (MQ=255)           aTCAAGGATTTGGTGCCGGTATGGAAGAGGATGCCTACGCTATTCGCGCCATTGCCAGCGCTGGATTACCCGCTCTGGTGAGCAATTCGTTCTCGAAAATTTTCTCCCTTTATGGCGAGCGCGTCGGCGGACTTTCTGtt                                                                                                        <  1:206458/140‑1 (MQ=255)             cAAGGATTTGGTGCCGGTATGGAAGAGGATGCCTACGCTATTCGCGCCATTGCCAGCGCTGGATTACCCGCTCTGGTGAGCAATTCGTTCTCGAAAATTTTCTCCCTTTATGGCGAGCGCGTCGGCGGACTTTCTGTTAt                                                                                                      >  1:143921/1‑140 (MQ=255)             cAAGGATTTGGTGCCGGTATGGAAGAGGATGCCTACGCTATTCGCGCCATTGCCAGCGCTGGATTACCCGCTCTGGTGAGCAATTCGTTCTCGAAAATTTTCTCCCTTTATGGCGAGCGCATCGGCGGACTTTCTGTTAt                                                                                                      >  1:56677/1‑140 (MQ=255)                    ttGGTGCCGGTATGGAAGAGGATGCCTACGCTATTCGCGCCATTGCCAGCGCTGGATTACCCGCTCTGGTGAGCAATTCGTTCTCGAAAATTTTCTCCCTTTATGGCGAGCGCGTCGGCGGACTTTCTGTTATGTGTGaa                                                                                               <  1:165877/140‑1 (MQ=255)                                  gAAGAGGATGCCTACGCTATTCGCGCCATTGCCAGCGCTGGATTACCCGCTCTGGTGAGCAATTCGTTCTCGAAAATTTTCTCCCTTTATGGCGAGCGCGTCGGCGGACTTTCTGTTATGTGTGAAGATGCCGAAGCCGc                                                                                 <  2:651360/140‑1 (MQ=255)                                  gAAGAGGATGCCTACGCTATTCGCGCCATTGCCAGCGCTGGATTACCCGCTCTGGTGAGCAATTCGTTCTCGAAAATTTTCTCCCTTTATGGCGAGCGCGTCGGCGGACTTTCTGTTATGTGTGAAGATGCCGAAGCCGc                                                                                 <  2:143921/140‑1 (MQ=255)                                       ggATGCCTACGCTATTCGCGCCATTGCCAGCGCTGGATTACCCGCTCTGGTGAGCAATTCGTTCTCGAAAATTTTCTCCCTTTATGGCGAGCGCGTCGGCGGACTTTCTGTTATGTGTGAAGATGCCGAAGCCGCTGGcc                                                                            <  1:581402/140‑1 (MQ=255)                                          tGCCTACGCTATTCGCGCCATTGCCAGCGCTGGATTACCCGCTCTGGTGAGCAATTCGTTCTCGAAAATTTTCTCCCTTTATGGCGAGCGCGTCGGCGGACTTTCTGTTATGTGTGAAGATGCCGAAGCCGCTGGCcgcg                                                                         >  2:499881/1‑140 (MQ=255)                                                      tCGCGCCATTGCCAGCGCTGGATTACCCGCTCTGGTGAGCAATTCGTTCTCGAAAATTTTCTCCCTTTATGGCGAGCGCGTCGGCGGACTTTCTGTTATGTGTGAAGATGCCGAAGCCGCTGGCCGCGTACTGGGGCAAt                                                             >  2:186429/1‑140 (MQ=255)                                                                       cTGGATTACCCGCTCTGGTGAGCAATTCGTTCTCGAAAATTTTCTCCCTTTATGGCGAGCGCATCGGCGGACTTTCTGTTATGTGTGAAGATGCCGAAGCCGCTGGCCGCGTACTGGGGCAATTGAAAGCAACAGTTcgc                                            <  2:56677/140‑1 (MQ=255)                                                                                  gCTCTGGTGAGCAATTCGTTCTCGAAAATTTTCTCCCTTTATGGCGAGCGCGTCGGCGGACTTTCTGTTATGTGTGAAGATGCCGAAGCCGCTGGCCGCGTACTGGGGCAATTGAAAGCAACAGTTCGCCGCAACTACTc                                 >  1:338863/1‑140 (MQ=255)                                                                                      tGGTGAGCAATTCGTTCTCGAAAATTTTCTCCCTTTATGGCGAGCGCGTCGGCGGACTTTCTGTTATGTGTGAAGATGCCGAAGCCGCTGGCCGCGTACTGGGGCAATTGAAAGCAACAGTTCGCCGCAACTACTCCAGc                             <  2:338863/140‑1 (MQ=255)                                                                                                                 tctcCCTTTATGGCGAGCGCGTCGGCGGACTTTCTGTTATGTGTGAAGATGCCGAAGCCGCTGGCCGCGTACTGGGGCAATTGAAAGCAACAGTTCGCCGCAACTACTCCAGCCCGCCGAATTTTGGTGCGCAGgtggtg  >  2:389118/1‑140 (MQ=255)                                                                                                                           |                                                                                                                                  GATATTGCCTATCAAGGATTTGGTGCCGGTATGGAAGAGGATGCCTACGCTATTCGCGCCATTGCCAGCGCTGGATTACCCGCTCTGGTGAGCAATTCGTTCTCGAAAATTTTCTCCCTTTACGGCGAGCGCGTCGGCGGACTTTCTGTTATGTGTGAAGATGCCGAAGCCGCTGGCCGCGTACTGGGGCAATTGAAAGCAACAGTTCGCCGCAACTACTCCAGCCCGCCGAATTTTGGTGCGCAGGTGGTG  >  NC_000913/4267747‑4267998

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