Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	NC_000913	3,006,380	A→G	D40G (GAC→GGC)	ygeW →	putative carbamoyltransferase YgeW

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	3,006,380	0	A	G	100.0%	115.0 / NA	33	D40G (GAC→GGC)	ygeW	putative carbamoyltransferase YgeW
Reads supporting (aligned to +/- strand):  ref base A (0/0);  new base G (10/23);  total (10/23)

ACTACCGGGACTTTATGATGAAAACTGTTAATGAGCTGATTAAGGATATCAATTCGCTGACCTCTCACCTTCACGAGAAAGATTTTTTGTTAACGTGGGAACAGACGCCAGATGAACTGAAACAAGTACTGGACGTTGCCGCAGCATTAAAAGCACTGCGTGCTGAAAACATCTCAACCAAAGTCTTTAATAGTGGATTAGGTATTTCCGTATTCCGCGACAACTCCACCCGTACCCGCTTCTCTTATGCTTCCGCGCTTAACCTGCT  >  NC_000913/3006248‑3006515                                                                                                                                     |                                                                                                                                        aCTACCGGGACTTTATGATGAAAACTGTTAATGAGCTGATTAAGGATATCAATTCGCTGACCTCTCACCTTCACGAGAAAGATTTTTTGTTAACGTGGGAACAGACGCCAGATGAACTGAAACAAGTACTGGGCGTTGcc                                                                                                                                  <  2:147212/140‑1 (MQ=255)  cTACCGGGACTTTATGATGAAAACTGTTAATGAGCTGATTAAGGATATCAATTCGCTGACCTCTCACCTTCACGAGAAAGATTTTTTGTTAACGTGGGAACAGACGCCAGATGAACTGAAACAAGTACTGGGCGTTGCCg                                                                                                                                 <  2:79726/140‑1 (MQ=255)   tACCGGGACTTTATGATGAAAACTGTTAATGAGCTGATTAAGGATATCAATTCGCTGACCTCTCACCTTCACGAGAAAGATTTTTTGTTAACGTGGGAACAGACGCCAGATGAACTGAAACAAGTACTGGGCGTTGCCgc                                                                                                                                <  1:89684/140‑1 (MQ=255)    aCCGGGACTTTATGATGAAAACTGTTAATGAGCTGATTAAGGATATCAATTCGCTGACCTCTCACCTTCACGAGAAAGATTTTTTGTTAACGTGGGAACAGACGCCAGATGAACTGAAACAAGTACTGGGCGTTGCCgca                                                                                                                               <  1:121040/140‑1 (MQ=255)     ccGGGACTTTATGATGAAAACTGTTAATGAGCTGATTAAGGATATCAATTCGCTGACCTCTCACCTTCACGAGAAAGATTTTTTGTTAACGTGGGAACAGACGCCAGATGAACTGAAACAAGTACTGGGCGTTGCCgcag                                                                                                                              >  2:113461/1‑140 (MQ=255)      cGGGACTTTATGATGAAAACTGTTAATGAGCTGATTAAGGATATCAATTCGCTGACCTCTCACCTTCACGAGAAAGATTTTTTGTTAACGTGGGAACAGACGCCAGATGAACTGAAACAAGTACTGGGCGTTGCCgcagc                                                                                                                             <  2:106487/140‑1 (MQ=255)           cTTTATGATGAAAACTGTTAATGAGCTGATTAAGGATATCAATTCGCTGACCTCTCACCTTCACGAGAAAGATTTTTTGTTAACGTGGGAACAGACGCCAGATGAACTGAAACAAGTACTGGGCGTTGCCGCAGCATTaa                                                                                                                        <  1:191768/140‑1 (MQ=255)            tttATGATGAAAACTGTTAATGAGCTGATTAAGGATATCAATTCGCTGACCTCTCACCTTCACGAGAAAGATTTTTTGTTAACGTGGGAACAGACGCCAGATGAACTGAAACAAGTACTGGGCGTTGCCgca                                                                                                                               <  2:134334/132‑1 (MQ=255)            tttATGATGAAAACTGTTAATGAGCTGATTAAGGATATCAATTCGCTGACCTCTCACCTTCACGAGAAAGATTTTTTGTTAACGTGGGAACAGACGCCAGATGAACTGAAACAAGTACTGGGCGTTGCCgca                                                                                                                               >  1:134334/1‑132 (MQ=255)                                   gCTGATTAAGGATATCAATTCGCTGACCTCTCACCTTCACGAGAAAGATTTTTTGTTAACGTGGGAACAGACGCCAGATGAACTGAAACAAGTACTGGGCGTTGCCGCAGCATTAAAAGCACTGCGTGCTGAAAACAtct                                                                                                <  2:63830/140‑1 (MQ=255)                                   gCTGATTAAGGATATCAATTCGCTGACCTCTCACCTTCACGAGAAAGATTTTTTGTTAACGTGGGAACAGACGCCAGATGAACTGAAACAAGTACTGGGCGTTGCCGCAGCATTAAAAGCACTGCGTGCTGAAAACAtct                                                                                                <  1:741230/140‑1 (MQ=255)                                    cTGATTAAGGATATCAATTCGCTGACCTCTCACCTTCACGAGAAAGATTTTTTGTTAACGTGGGAACAGACGCCAGATGAACTGAAACAAAAACTGGGCGCTGCTGCAGCATTAAAACCACTACGTGCTGAAAACAtctc                                                                                               >  2:259012/1‑140 (MQ=255)                                         tAAGGATATCAATTCTCTGACCCCTCACCCTCACGAGAAAGATTTTTTGTTAACGTGGGAACAGACGCCAGATGAACTGAAACAAATACTGGGCGTTGCCGCAGCATTAAAAGCACTGCGTGCTGAAAACATCTCAACCa                                                                                          <  1:259012/140‑1 (MQ=255)                                         tAAGGATATCAATTCGCTGACCTCTCACCTTCACGAGAAAGATTTTTTGTTAACGTGGGAACAGACGCCAGATGAACTGAAACAAGTACTGGGCGTTGCCGCAGCATTAAAAGCACTGCGTGCTGAAAACATCTCAACCa                                                                                          <  2:486705/140‑1 (MQ=255)                                         tAAGGATATCAATTCGCTGACCTCTCACCTTCACGAGAAAGATTTTTTGTTAACGTGGGAACAGACGCCAGATGAACTGAAACAAGTACTGGGCGTTGCCGCAGCATTAAAAGCACTGCGTGCTGAAAACATCTCAACCa                                                                                          <  2:471688/140‑1 (MQ=255)                                           aGGATATCAATTCGCTGACCTCTCACCTTCACGAGAAAGATTTTTTGTTAACGTGGGAACAGACGCCAGATGAACTGAAACAAGTACTGGGCGTTGCCGCAGCATTAAAAGCACTGCGTGCTGAAAACATCTCAACCaaa                                                                                        <  2:840964/140‑1 (MQ=255)                                           aGGATATCAATTCGCTGACCTCTCACCTTCACGAGAAAGATTTTTTGTTAACGTGGGAACAGACGCCAGATGAACTGAAACAAGTACTGGGCGTTGCCGCAGCATTAAAAGCACTGCGTGCTGAAAACATCTCAACCaaa                                                                                        <  1:113461/140‑1 (MQ=255)                                                 tCAATTCGCTGACCTCTCACCTTCACGAGAAAGATTTTTTGTTAACGTGGGAACAGACGCCAGATGAACTGAAACAAGTACTGGGCGTTGCCGCATGATTAAAAGGACTGCGTGGTGAAAACTTCTCGACTAGCGTCttt                                                                                  >  1:535442/1‑140 (MQ=255)                                                 tCAATTCGCTGACCTCTCACCTTCACGAGAAAGATTTTTTGTTAACGTGGGAACAGACGCCAGATGAACTGAAACAAGTACTGGGCGTTGCCGCAGCATTAAAAGCACTGCGTGCTGAAAACATCTCAACCAAAGTCttt                                                                                  >  2:845929/1‑140 (MQ=255)                                                 tCAATTCGCTGACCTCTCACCTTCACGAGAAAGATTTTTTGTTAACGTGGGAACAGACGCCAGATGAACTGAAACAAGTACTGGGCGTTGCCGCAGCATTAAAAGCACTGCGTGCTGAAAACATCTCAACCAAAGTCttt                                                                                  >  2:95845/1‑140 (MQ=255)                                                 tCAATTCGCTGACCTCTCACCTTCACGAGAAAGATTTTTTGTTAACGTGGGAACAGACGCCAGATGAACTGAAACAAGTACTGGGCGTTGCCGCAGCATTAAAAGCACTGCGTGCTGAAAACATCTCAACCAAAGTCttt                                                                                  >  1:142437/1‑140 (MQ=255)                                                 tCAATTCGCTGACCTCTCACCTTCACGAGAAAGATTTTTTGTTAACGTGGGAACAGACGCCAGATGAACTGAAACAAGTACTGGGCGTTGCAGCAGCATTAAAAGCACCGCGTGCTGCAAACGTCTCAACCAAAGTCttt                                                                                  >  2:283004/1‑140 (MQ=255)                                                    aTTCGCTGACCTCTCACCTTCACGAGAAAGATTTTTTGTTAACGTGGGAACAGACGCCAGATGAACTGAAACAAGTACTGGGCGTTGCCGCAGCATTAAAAGCACTGCGTGCTGAAAACATCTCAACCAAAGTCTTTAAt                                                                               >  1:677539/1‑140 (MQ=255)                                                     ttCGCTGACCTCTCACCTTCACGAGAAAGATTTTTTGTTAACGTGGGAACAGACGCCAGATGAACTGAAACAAGTACTGGGCGTTGCCGCAGCATTAAAAGCACTGCGTGCTGAAAACATCTCAACCAAAGTCTTTAATa                                                                              <  1:95845/140‑1 (MQ=255)                                                     ttCGCTGACCTCTCACCTTCACGAGAAAGATTTTTTGTTAACGTGGGAACAGACGCCAGATGAACTGAAACAAGTACTGGGCGTTGCCGCAGCATTAAAAGCACTGCGTGCTGAAAACATCTCAACCAAAGTCTTTAATa                                                                              <  1:72553/140‑1 (MQ=255)                                                       cGCTGACCTCTCACCTTCACGAGAAAGATTTTTTGTTAACGTGGGAACAGACGCCAGATGAACTGAAACAAGTACTGGGCGTTGCCGCAGCATTAAAAGCACTGCGTGCTGAAAACATCTCAACCAAAGTCTTTAATAGt                                                                            <  1:845929/140‑1 (MQ=255)                                                                                  aTTTTTTGTTAACGTGGGAACAGACCCCACATGAACTTAAACAAGTACTGGGCGTTGCCGCAGCATTAAAAGCACTGCGTGCTGAAAACATCTCAACCAAAGTCTTTAATAGTGGATTAGGTATTTCCGTATTCCGCGAc                                                 <  1:336663/140‑1 (MQ=255)                                                                                               gTGGGAACAGACGCCTGATGAACTGCAACAAGTACTGCGCGTTGCCGCAGCATTAAAAGCACTGCGTGCTGAAAACATCTCAACCAAAGTCTTTAATAGTGGATTAGGTATTTCCGTATTCCGCGACAACTCCACCCGTa                                    <  1:283004/140‑1 (MQ=255)                                                                                                   gAACAGACGCCAGATGAACTGAAACAAGTACTGGGCGTTGCCGCAGCATTAAAAGCACTGCGTGCTGAAAACATCTCAACCAAAGTCTTTAATAGTGGATTAGGTATTTCCGTATTCCGCGACAACTCCACCCGTACCCg                                <  1:566609/140‑1 (MQ=255)                                                                                                                  gAACTGAAACAAGTACTGGGCGTTGCCGCAGCATTAAAAGCACTGCGTGCTGAAAACATCTCAACCAAAGTCTTTAATAGTGGATTAGGTATTTCCGTATTCCGCGACAACTCCACCCGTACCCGCTTCTCTTATGCTTc                 <  2:142437/140‑1 (MQ=255)                                                                                                                   aaCTGAAACAAGTACTGGGCGTTGCCGCAGCATTAAAAGCACTGCGTGCTGAAAACATCTCAACCAAAGTCTTTAATAGTGGATTAGGTATTTCCGTATTGCGCCAAAATTCCACCCGTACCGACTTCTCTGATGCTTcc                >  2:415610/1‑140 (MQ=255)                                                                                                                                 cTGGGCGTTGCCGCAGCATTAAGAGTACTACGTGCTGAAAACATCTCAACCAAAGTCTTTAATAGTGGATTAGGTATTTCCGTATTCCGCGACAACTCCACCCGTACCCGCTTCTCTTATGCTTCCGCGCTTAACCTGCt  <  1:415610/140‑1 (MQ=255)                                                                                                                                 cTGGGCGTTGCCGCAGCATTAAAAGCACTGCGTGCTGAAAACATCTCAACCAAAGTCTTTAATAGTGGATTAGGTATTTCCGTATTCCGCGACAACCCCCCCCGTCCCCGCTTCTCTTATGCTTCCGCGCTTAACCTGCt  <  2:628477/140‑1 (MQ=255)                                                                                                                                     |                                                                                                                                        ACTACCGGGACTTTATGATGAAAACTGTTAATGAGCTGATTAAGGATATCAATTCGCTGACCTCTCACCTTCACGAGAAAGATTTTTTGTTAACGTGGGAACAGACGCCAGATGAACTGAAACAAGTACTGGACGTTGCCGCAGCATTAAAAGCACTGCGTGCTGAAAACATCTCAACCAAAGTCTTTAATAGTGGATTAGGTATTTCCGTATTCCGCGACAACTCCACCCGTACCCGCTTCTCTTATGCTTCCGCGCTTAACCTGCT  >  NC_000913/3006248‑3006515

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