Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	NC_000913	4,044,390	(T)6→5	coding (192/1431 nt)	yihF →	DUF945 family protein

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	4,044,385	0	T	.	100.0%	81.1 / NA	19	coding (187/1431 nt)	yihF	DUF945 family protein
Reads supporting (aligned to +/- strand):  ref base T (0/0);  new base . (8/11);  total (8/11)

TACGCAAATTCAGCCTGGTGTCGAAAAATTTATTAAAGATTTTAACGATGCTAAAAAGAAAGGTGAACATGCCTACGATATGACGTTAAGTTATCAAAATTTTGACAAAGGCTTTTTTAATTCCCGTTTTCAAATGCAAATGACATTCGATAACGGTGCACCCGATCTCAATATCAAGCCAGGCCAGAAAGTTGTATTTGATGTGGATGTTGAGCACGGTCCGTTGCCCATCACAATGTTAAT  >  NC_000913/4044273‑4044515                                                                                                                 |                                                                                                                                   tACGCAAATTCAGCCTGGTGTCGAAAAATTTATTAAAGATTTTAACGATGCTAAAAAGAAAGGTGAACATGCCTACGATATGACGTTAAGTTATCAAAATTTTGACAAAGGC‑TTTTTAATTCCCGTTTTCAAATGCAAATGACATTCGa                                                                                               <  1:53472/149‑1 (MQ=255) tACGCAAATTCAGCCTGGTGTCGAAAAATTTATTAAAGATTTTAACGATGCTAAAAAGAAAGGTGAACATGCCTACGATATGACGTTAAGTTATCAAAATTTTGACAAAGGC‑TTTTTAATTCCCGTTTTCAAATGCAAATGACATTCGa                                                                                               <  2:77943/149‑1 (MQ=255)    gCAAATTCAGCCTGGTGTCGAAAAATTTATTAAAGATTTTAACGATGCTAAAAAGAAAGGTGAACATGCCTACGATATGACGTTAAGTTATCAAAATTTTGACAAAGGC‑TTTTTAATTCCCGTTTTCAAATGCAAATGACATTCGATaa                                                                                            <  1:85083/149‑1 (MQ=255)                       gAAAAATTTATTAAAGATTTTAACGATGCTAAAAAGAAAGGTGAACATGCCTACGATATGACGTTAAGTTATCAAAATTTTGACAAAGGC‑TTTTTAATTCCCGTTTTCAAATGCAAATGACATTCGATAACGGTGCACCCGATCTCaa                                                                          >  5:62093/1‑148 (MQ=255)                           aaTTTATTAAAGATTTTAACGATGCTAAAAAGAAAGGTGAACATGCCTACGATATGACGTTAAGTTATCAAAATTTTGACAAAGGC‑TTTTTAATTCCCGTTTTCAAATGCAAATGACATTCGATAACGGTGCACCCGATCTCAATATCa                                                                     >  7:109105/1‑149 (MQ=255)                                 ttaAAGATTTTAACGATGCTAAAAAGAAAGGTGAACATGCCTACGATATGACGTTAAGTTATCAAAATTTTGACAAAGGC‑TTTTTAATTCCCGTTTTCAAATGCAAATGACATTCGATAACGGTGCACCCGATCTCAATATCAAgccag                                                               >  4:129278/1‑149 (MQ=255)                                   aaaGATTTTAACGATGCTAAAAAGAAAGGTGAACATGCCTACGATATGACGTTAAGTTATCAAAATTTTGACAAAGGC‑TTTTTAATTCCCGTTTTCAAATGCAAATGACATTCGATAACGGTGCACCCGAtc                                                                              <  1:119922/132‑1 (MQ=255)                                   aaaGATTTTAACGATGCTAAAAAGAAAGGTGAACATGCCTACGATATGACGTTAAGTTATCAAAATTTTGACAAAGGC‑TTTTTAATTCCCGTTTTCAAATGCAAATGACATTCGATAACGGTGCACCCGAtc                                                                              >  2:119922/1‑132 (MQ=255)                                                      aaaaGAAAGGTGAACATGCCTACGATATGACGTTAAGTTATCAAAATTTTGACAAAGGC‑TTTTTAATTCCCGTTTTCAAATGCAAATGACATTCGATAACGGTGCACCCGATCTCAATATCAAGCCAGGCCAGAAAGTTGTATTTGa                                            <  8:109105/147‑1 (MQ=255)                                                               gTGAACATGCCTACGATATGACGTTAAGTTATCAAAATTTTGACAAAGGC‑TTTTTAATTCCCGTTTTCAAATGCAAATGACATTCGATAACGGTGCACCCGATCTCAATATCa                                                                     >  1:47598/1‑113 (MQ=255)                                                               gTGAACATGCCTACGATATGACGTTAAGTTATCAAAATTTTGACAAAGGC‑TTTTTAATTCCCGTTTTCAAATGCAAATGACATTCGATAACGGTGCACCCGATCTCAATATCa                                                                     <  2:47598/113‑1 (MQ=255)                                                                 gAACATGCCTACGATATGACGTTAAGTTATCAAAATTTTGACAAAGGC‑TTTTTAATTCCCGTTTTCAAATGCAAATGACATTCGATAACGGTGCACCCGATCTCAATATCAAGCCAGGCCAGAAAGTTGTATTTGATGTGGATGtt                                  >  7:97209/1‑146 (MQ=255)                                                                 gAACATGCCTACGATATGACGTTAAGTTATCAAAATTTTGACAAAGGC‑TTTTTAATTCCCGTTTTCAAATGCAAATGACATTCGATAACGGTGCACCCGATCTCAATATCAAGCCAGGCCAGAAAGTTGTATTTGATGTGGATGTTGa                                <  8:97209/148‑1 (MQ=255)                                                                  aaCATGCCTACGATATGACGTTAAGTTATCAAAATTTTGACAAAGGC‑TTTTTAATTCCCGTTTTCAAATGCAAATGACATTCGATAACGGTGCACCCGATCTCAATATCAAGCCAGGCCAGAAAGtt                                                    >  5:99895/1‑127 (MQ=255)                                                                  aaCATGCCTACGATATGACGTTAAGTTATCAAAATTTTGACAAAGGC‑TTTTTAATTCCCGTTTTCAAATGCAAATGACATTCGATAACGGTGCACCCGATCTCAATATCAAGCCAGGCCAGAAAGtt                                                    <  6:99895/127‑1 (MQ=255)                                                                   aCATGCCTACGATATGACGTTAAGTTATCAAAATTTTGACAAAGGC‑TTTTTAATTCCCGTTTTCAAATGCAAATGACATTCGATAACGGTGCACCCGATCTCAATATCAAGCCAGGCCAGAAAGTTGTATTTGATGTGGATGTTGAGCa                             <  3:129278/149‑1 (MQ=255)                                                                         cTACGATATGACGTTAAGTTATCAAAATTTTGACAAAGGC‑TTTTTAATTCCCGTTTTCAAATGCAAATGACATTCGATAACGGTGCACCCGATCTCAATATCAAGCCAGGCCAGAAAGTTGTATTTGATGTGGATGTTGAGCACGGTcc                       >  7:48985/1‑149 (MQ=255)                                                                              atatGACGTTAAGTTATCAAAATTTTGACAAAGGC‑TTTTTAATTCCCGTTTTCAAATGCAAATGACATTCGATAACGGTGCACCCGATCTCAATATCAAGCCAGGCCAGAAAGTTGTATTTGATGTGGATGTTGAGCACGGTCCGTTGc                  <  6:62093/149‑1 (MQ=255)                                                                                              tCAAAATTTTGACAAAGGC‑TTTTTAATTCCCGTTTTCAAATGCAAATGACATTCGATAACGGTGCACCCGATCTCAATATCAAGCCAGGCCAGAAAGTTGTATTTGATGTGGATGTTGAGCACGGTCCGTTGCCCATCACAATGTTAAt  <  8:48985/149‑1 (MQ=255)                                                                                                                 |                                                                                                                                   TACGCAAATTCAGCCTGGTGTCGAAAAATTTATTAAAGATTTTAACGATGCTAAAAAGAAAGGTGAACATGCCTACGATATGACGTTAAGTTATCAAAATTTTGACAAAGGCTTTTTTAATTCCCGTTTTCAAATGCAAATGACATTCGATAACGGTGCACCCGATCTCAATATCAAGCCAGGCCAGAAAGTTGTATTTGATGTGGATGTTGAGCACGGTCCGTTGCCCATCACAATGTTAAT  >  NC_000913/4044273‑4044515

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