Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	NC_000913	4,261,497	T→G	intergenic (+191/‑172)	yjbM → / → dusA	uncharacterized protein/tRNA‑dihydrouridine synthase A

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	4,261,497	0	T	G	100.0%	109.9 / NA	31	intergenic (+191/‑172)	yjbM/dusA	uncharacterized protein/tRNA‑dihydrouridine synthase A
Reads supporting (aligned to +/- strand):  ref base T (0/0);  new base G (10/21);  total (10/21)

ACCCAATATTTATTGGGTTCTTTTTTCTCTATCAATACTATTAGCAGGGAGATATATCACCAGAGTTTAATGTGTGATTTTTTATTTATCGTCGAACCTGGATTGTTTATCATTGGCCTTAACAAAGTTAACGGCTAATAAGATTATTTCCATCACTTCGTGAGAGCTTCATGCCTTGAGAGGATCTCAATTTTCTTTTGCAATGAGACAGGCGCTTCCTGTTGTTATGGTATAGTACCCCGCTATTGAGCCTCCTGAATAGTGATGCTGAATAACATAACCCCATGATATATC  >  NC_000913/4261350‑4261643                                                                                                                                                    |                                                                                                                                                   aCCCAATATTTATTGGGTTCTTTTTTCTCTATCAATACTATTAGCAGGGAGATATATCACCAGAGTTTAATGTGTGATTTTTTATTTATCGTCGAACCTGGATTGTTTATCATTGGCCTTAACAAAGTTAACGGCTAATAAGATTATGTc                                                                                                                                                  <  2:174390/150‑1 (MQ=255)      atatTTATTGGGTTCTTTTTTCTCTATCAATACTATTAGCAGGGAGATATATCACCAGAGTTTAATGTGTGATTTTTTATTTATCGTCGAACCTGGATTGTTTATCATTGGCCTTAACAAAGTTAACGGCTAATAAGATTATGTCCATCAc                                                                                                                                            <  2:45964/151‑1 (MQ=255)      atatTTATTGGGTTCTTTTTTCTCTATCAATACTATTAGCAGGGAGATATATCACCAGAGTTTAATGTGTGATTTTTTATTTATCGTCGAACCTGGATTGTTTATCATTGGCCTTAACAAAGTTAACGGCTAATAAGATTATGTCCATCAc                                                                                                                                            <  1:27494/151‑1 (MQ=255)                  ttCTTTTTTCTCTATCAATACTATTAGCAGGGAGATATATCACCAGAGTTTAATGTGTGATTTTTTATTTATCGTCGAACCTGGATTGTTTATCATTGGCCTTAACAAAGTTAACGGCTAATAAGATTATGTCCATCACTTCGTGAGAGc                                                                                                                                 <  2:47646/150‑1 (MQ=255)                             ctATCAATACTATTAGCAGGGAGATATATCACCAGAGTTTAATGTGTGATTTTTTATTTATCGTCGAACCTGGATTGTTTATCATTGGCCTTAACAAAGTTAACGGCTAATAAGATTATGTCCATCACTTCGTGAGAGCTTCATGCCTTga                                                                                                                     >  1:52008/1‑151 (MQ=255)                                tCAATACTATTAGCAGGGAGATATATCACCAGAGTTTAATGTGTGATTTTTTATTTATCGTCGAACCTGGATTGTTTATCATTGGCCTTAACAAAGTTAACGGCTAATAAGATTATGTCCATCACTTCGTGAGAGCTTCATGCCTTGagag                                                                                                                  <  2:201375/151‑1 (MQ=255)                                 cAATACTATTAGCAGGGAGATATATCACCAGAGTTTAATGTGTGATTTTTTATTTATCGTCGAACCTGGATTGTTTATCATTGGCCTTAACAAAGTTAACGGCTAATAAGATTATGTCCATCACTTCGTGAGAGCTTCATGCCTTGAGAgg                                                                                                                 <  2:166817/151‑1 (MQ=255)                                   aTACTATTAGCAGGGAGATATATCACCAGAGTTTAATGTGTGATTTTTTATTTATCGTCGAACCTGGATTGTTTATCATTGGCCTTAACAAAGTTAACGGCTAATAAGATTATGTCCATCACTTCGTGAGAGCTTCATGCCTTGAGAGGAt                                                                                                               >  2:103269/1‑151 (MQ=255)                                          tAGCAGGGAGATATATCACCAGAGTTTAATGTGTGATTTTTTATTTATCGTCGAACCTGGATTGTTTATCATTGGCCTTAACAAAGTTAACGGCTAATAAGATTATGTCCATCACTTCGTGAGAGCTTCATGCCTTGAGAGGATCTCAAtt                                                                                                        >  1:316275/1‑151 (MQ=255)                                                ggAGATATATCACCAGAGTTTAATGTGTGATTTTTTATTTATCGTCGAACCTGGATTGTTTATCATTGGCCTTAACAAAGTTAACGGCTAATAAGATTATGTCCATCACTTCGTGAGAGCTTCATGCCTTGAGAGGATCTCAATTTTCttt                                                                                                  >  2:210905/1‑151 (MQ=255)                                                 gagaTATATCACCAGAGTTTAATGTGTGATTTTTTATTTATCGTCGAACCTGGATTGTTTATCATTGGCCTTAACAAAGTTAACGGCTAATAAGATTATGTCCATCACTTCGTGAGAGCTTCATGCCTTGAGAGGATCTCAATTTTCtttt                                                                                                 >  1:76844/1‑151 (MQ=255)                                                                           tgATTTTTTATTTATCGTCGAACCTGGATTGTTTATCATTGGCCTTAACAAAGTTAACGGCTAATAAGATTATGTCCATCACTTCGTGAGAGCTTCATGCCTTGAGAGGATCTCAATTTTCTTTTGCAATGAGACAGGCGCTTCCtgttgt                                                                       <  1:226390/151‑1 (MQ=255)                                                                                      ttatCGTCGAACCTGGATTGTTTATCATTGGCCTTAACAAAGTTAACGGCTAATAAGATTATGTCCATCACTTCGTGAGAGCTTCATGCCTTGAGAGGATCTCAATTTTCTTTTGCAATGAGACAGGCGCTTCCTGTTGTTATGGTATAGt                                                            <  1:231812/151‑1 (MQ=255)                                                                                        atCGTCGAACCTGGATTGTTTATCATTGGCCTTAACAAAGTTAACGGCTAATAAGATTATGTCCATCACTTCGTGAGAGCTTCATGCCTTGAGAGGATCTCAATTTTCTTTTGCAATGAGACAGGCGCTTCCTGTTGTTATGGTATAGTAc                                                          <  2:91631/151‑1 (MQ=255)                                                                                        atCGTCGAACCTGGATTGTTTATCATTGGCCTTAACAAAGTTAACGGCTAATAAGATTATGTCCATCACTTCGTGAGAGCTTCATGCCTTGAGAGGATCTCAATTTTCTTTTGCAATGAGACAGGCGCTTCCTGTTGTTATGGTATAGTAc                                                          <  2:25705/151‑1 (MQ=255)                                                                                                      aTTGTTTATCATTGGCCTTAACAAAGTTAACGGCTAATAAGATTATGTCCATCACTTCGTGAGAGCTTCATGCCTTGAGAGGATCTCAATTTTCTTTTGCAATGAGACAGGCGCTTCCTGTTGTTATGGTATAGTACCCCGCTATTGAGc                                             <  2:132757/150‑1 (MQ=255)                                                                                                         gTTTATCATTGGCCTTAACAAAGTTAACGGCTAATAAGATTATGTCCATCACTTCGTGAGAGCTTCATGCCTTGAGAGGATCTCAATTTTCTTTTGCAATGAGACAGGCGCTTCCTGTTGTTATGGTATAGTACCCCGCTATTGAGcctc                                          >  1:181936/1‑150 (MQ=255)                                                                                                         gTTTATCATTGGCCTTAACAAAGTTAACGGCTAATAAGATTATGTCCATCACTTCGTGAGAGCTTCATGCCTTGAGAGGATCTCAATTTTCTTTTGCAATGAGACAGGCGCTTCCTGTTGTTATGGTATAGTACCCCGCTATTGAGcctc                                          >  2:306920/1‑150 (MQ=255)                                                                                                              tCATTGGCCTTAACAAAGTTAACGGCTAATAAGATTATGTCCATCACTTCGTGAGAGCTTCATGCCTTGAGAGGATCTCAATTTTCTTTTGCAATGAGACAGGCGCTTCCTGTTGTTATGGTATAGTACCCCGCTATTGAGCCTCCTGAAt                                    <  1:150973/151‑1 (MQ=255)                                                                                                              tCATTGGCCTTAACAAAGTTAACGGCTAATAAGATTATGTCCATCACTTCGTGAGAGCTTCATGCCTTGAGAGGATCTCAATTTTCTTTTGCAATGAGACAGGCGCTTCCTGTTGTTATGGTATAGTACCCCGCTATTGAGCCTCCTGAAt                                    <  2:107422/151‑1 (MQ=255)                                                                                                                aTTGGCCTTAACAAAGTTAACGGCTAATAAGATTATGTCCATCACTTCGTGAGAGCTTCATGCCTTGAGAGGATCTCAATTTTCTTTTGCAATGAGACAGGCGCTTCCTGTTGTTATGGTATAGTACCCCGCTATTGAGCCTCCTGAATAg                                  <  1:29355/151‑1 (MQ=255)                                                                                                                aTTGGCCTTAACAAAGTTAACGGCTAATAAGATTATGTCCATCACTTCGTGAGAGCTTCATGCCTTGAGAGGATCTCAATTTTCTTTTGCAATGAGACAGGCGCTTCCTGTTGTTATGGTATAGTACCCCGCTATTGAGCCTCCTGAATAg                                  <  1:268051/151‑1 (MQ=255)                                                                                                                   ggCCTTAACAAAGTTAACGGCTAATAAGATTATGTCCATCACTTCGTGAGAGCTTCATGCCTTGAGAGGATCTCAATTTTCTTTTGCAATGAGACAGGCGCTTCCTGTTGTTATGGTATAGTACCCCGCTATTGAGCCTCCTGAATAGTGa                               >  2:55708/1‑151 (MQ=255)                                                                                                                      cTTAACAAAGTTAACGGCTAATAAGATTATGTCCATCACTTCGTGAGAGCTTCATGCCTTGAGAGGATCTCAATTTTCTTTTGCAATGAGACAGGCGCTTCCTGTTGTTATGGTATAGTACCCCGCTATTGAGCCTCCTGAATAGTGATGc                            <  2:316273/151‑1 (MQ=255)                                                                                                                                ttAACGGCTAATAAGATTATGTCCATCACTTCGTGAGAGCTTCATGCCTTGAGAGGATCTCAATTTTCTTTTGCAATGAGACAGGCGCTTCCTGTTGTTATGGTATAGTACCCCGCTATTGAGCCTCCTGAATAGTGATGCTGAataacat                  <  1:257571/151‑1 (MQ=255)                                                                                                                                ttAACGGCTAATAAGATTATGTCCATCACTTCGTGAGAGCTTCATGCCTTGAGAGGATCTCAATTTTCTTTTGCAATGAGACAGGCGCTTCCTGTTGTTATGGTATAGTACCCCGCTATTGAGCCTCCTGAATAGTGATGCTGAataacat                  <  1:204270/151‑1 (MQ=255)                                                                                                                                ttAACGGCTAATAAGATTATGTCCATCACTTCGTGAGAGCTTCATGCCTTGAGAGGATCTCAATTTTCTTTTGCAATGAGACAGGCGCTTCCTGTTGTTATGGTATAGTACCCCGCTATTGAGCCTCCTGAATAGTGATGCTGAataac                    <  2:19177/149‑1 (MQ=255)                                                                                                                                      gCTAATAAGATTATGTCCATCACTTCGTGAGAGCTTCATGCCTTGAGAGGATCTCAATTTTCTTTTGCAATGAGACAGGCGCTTCCTGTTGTTATGGTATAGTACCCCGCTATTGAGCCTCCTGAATAGTGATGCTGAATAACATAAccc             <  2:285245/150‑1 (MQ=255)                                                                                                                                           taaGATTATGTCCATCACTTCGTGAGAGCTTCATGCCTTGAGAGGATCTCAATTTTCTTTTGCAATGAGACAGGCGCTTCCTGTTGTTATGGTATAGTACCCCGCTATTGAGCCTCCTGAATAGTGATGCTGAATAACATAACCCCATGat       <  1:265198/151‑1 (MQ=255)                                                                                                                                             aGATTATGTCCATCACTTCGTGAGAGCTTCATGCCTTGAGAGGATCTCAATTTTCTTTTGCAATGAGACAGGCGCTTCCTGTTGTTATGGTATAGTACCCCGCTATTGAGCCTCCTGAATAGTGATGCTGAATAACATAACCCCATGatat     >  2:136286/1‑151 (MQ=255)                                                                                                                                              gATTATGTCCATCACTTCGTGAGAGCTTCATGCCTTGAGAGGATCTCAATTTTCTTTTGCAATGAGACAGGCGCTTCCTGTTGTTATGGCATAGTACCCCGCTATTGAGCCTCCTGAATAGTGATGCTGAATAACATAACCCCATGatata    >  2:140815/1‑151 (MQ=255)                                                                                                                                               attatgTCCATCACTTCGTGAGAGCTTCATGCCTTGAGAGGATCTCAATTTTCTTTTGCAATGAGACAGGCGCTTCCTGTTGTTATGGTATAGTACCCCGCTATTGAGCCTCCTGAATAGTGATGCTGAATAACATAACCCCATGatat     <  1:92198/149‑1 (MQ=255)                                                                                                                                                ttatgTCCATCACTTCGTGAGAGCTTCATGCCTTGAGAGGATCTCAATTTTCTTTTGCAATGAGACAGGCGCTTCCTGTTGTTATGGTATAGTACCCCGCTATTGAGCCTCCTGAATAGTGATGCTGAATAACATAACCCCATGATATATc  <  2:306273/151‑1 (MQ=255)                                                                                                                                                    |                                                                                                                                                   ACCCAATATTTATTGGGTTCTTTTTTCTCTATCAATACTATTAGCAGGGAGATATATCACCAGAGTTTAATGTGTGATTTTTTATTTATCGTCGAACCTGGATTGTTTATCATTGGCCTTAACAAAGTTAACGGCTAATAAGATTATTTCCATCACTTCGTGAGAGCTTCATGCCTTGAGAGGATCTCAATTTTCTTTTGCAATGAGACAGGCGCTTCCTGTTGTTATGGTATAGTACCCCGCTATTGAGCCTCCTGAATAGTGATGCTGAATAACATAACCCCATGATATATC  >  NC_000913/4261350‑4261643

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