Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	NC_000913	734,601	C→T	pseudogene (382/1434 nt)	ybfO →	conserved protein, rhs‑like YbfO

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	734,601	0	C	T	100.0%	50.1 / NA	24	pseudogene (382/1434 nt)	ybfO	conserved protein, rhs‑like YbfO
Reads supporting (aligned to +/- strand):  ref base C (0/0);  new base T (13/11);  total (13/11)

GGTGTGGCGACGTGAACGGGACCTGACGGGCTGGATGTCGCTGTCACGGAAACCGCAAGTGACCTGGTACGGCTGGGACGGCGACCGCCTGACCACGATACAGAACGACAGAACCCGCATCCAGACGATTTATCAGCCGGGGAGCTTCACGCCACTCATCAGGGTTGAAACCGCCACCGGTGAGCTGGCGAAAACGCAGCGCCGCAGCCTGGCGGATACCCTTCAGCAGTCCGGCGGCGAAGACGGTGGCAGTGTGGTGTTCCCGCCGGTG  >  NC_000913/734465‑734735                                                                                                                                         |                                                                                                                                       ggTGTGGCGACGTGAACGGGACCTGACGGGCTGGATGTCGCTGTCACGGAAACCGCAAGTGACCTGGTACGGCTGGGACGGCGACCGCCTGACCACGATACAGAACGACAGAACCCGCATCCAGACGATTTATCAGTCgg                                                                                                                                     >  1:684529/1‑140 (MQ=37)  gtgtGGCGACGTGAACGGGACCTGACGGGCTGGATGTCGCTGTCACGGAAACCGCAAGTGACCTGGTACGGCTGGGACGGCAACCGCCTGACCACGATACAGAACGACAGAACCCGCATCCAGACGATTTATCAGTCggg                                                                                                                                    <  1:682545/140‑1 (MQ=19)   tgtgGCGACGTGAACGGGACCTGACGGGCTGGATGTCGCTGTCACGGAAACCGCAAGTGACCTGGTACGGCTGGGACGGCGACCGCCTGACCACGATACAGAACGACAGAACCCGCATCCAGACGATTTATCAGTCgggg                                                                                                                                   <  1:713929/140‑1 (MQ=18)       gCGACGTGAACGGGACCTGACGGGCTGGATGTCGCTGTCACGGAAACCGCAAGTGACCTGGTACGGCTGGGACGGCGACCGCCTGACCACGATACAGAACGACAGAACCCGCATCCAGACGATTTATCAGTCGGGGAGCt                                                                                                                               >  1:73524/1‑140 (MQ=17)          aCGTGAACGGGACCTGACGGGCTGGATGTCGCTGTCACGGAAACCGCAAGTGACCTGGTACGGCTGGGACGGCGACCGCCTGACCACGATACAGAACGACAGAACCCGCATCCAGACGATTTATCAGTCGGGGAGCTTCa                                                                                                                            <  2:878096/140‑1 (MQ=11)          aCGTGAACGGGACCTGACGGGCTGGATGTCGCTGTCACGGAAACCGCAAGTGACCTGGTACGGCTGGGACGGCGACCGCCTGACCACGATACAGAACGACAGAACCCGCATCCAGACGATTTATCAGTCGGGGAGCTTCa                                                                                                                            <  1:460672/140‑1 (MQ=11)                                                                                                      aGAACGACAGAACCCGCATCCAGACGATTTATCAGTCGGGGAGCTTCACGCCACTCATCAGGGTTGAAACCGCCACCGGTGAGCTGGCGAAAACGCAGCGCCGCAGCCTGGCGGATACCCTTCAGCAGTCCGGCGGCGaa                                <  2:641232/140‑1 (MQ=18)                                                                                                      aGAACGACAGAACCCGCATCCAGACGATTTATCAGTCGGGGAGCTTCACGCCACTCATCAGGGTTGAAACCGCCACCGGTGAGCTGGCGAAAACGCAGCGCCGCAGCCTGGCGGATACCCTTCAGCAGTCCGGCGGCGaa                                >  2:626559/1‑140 (MQ=17)                                                                                                         aCGACAGAACCCGCATCCAGACGATTTATCAGTCGGGGAGCTTCACGCCACTCATCAGGGTTGAAACCGCCACCGGTGAGCTGGCGAAAACGCAGCGCCGCAGCCTGGCGGATACCCTTCAGCAGTCCGGCGGCGAAGAc                             >  2:890287/1‑140 (MQ=17)                                                                                                           gACAGAACCCGCATCCAGACGATTTATCAGTCGGGGAGCTTCACGCCACTCATCAGGGTTGAAACCGCCACCGGTGAGCTGGCGAAAACGCAGCGCCGCAGCCTGGCGGATACCCTTCAGCAGTCCGGCGGCGAAGACgg                           <  1:431362/140‑1 (MQ=17)                                                                                                           gACAGAACCCGCATCCAGACGATTTATCAGTCGGGGAGCTTCACGCCACTCATCAGGGTTGAAACCGCCACCGGTGAGCTGGCGAAAACGCAGCGCCGCAGCCTGGCGGATACCCTTCAGCAGTCCGGCGGCGAAGACgg                           <  2:154648/140‑1 (MQ=17)                                                                                                            aCAGAACCCGCATCCAGACGATTTATCAGTCGGGGAGCTTCACGCCACTCATCAGGGTTGAAACCGCCACCGGTGAGCTGGCGAAAACGCAGCGCCGCAGCCTGGCGGATACCCTTCAGCAGTCCGGCGGCGAAGACGGt                          >  1:303811/1‑140 (MQ=17)                                                                                                             cAGAACCCGCATCCAGACGATTTATCAGTCGGGGAGCTTCACGCCACTCATCAGGGTTGAAACCGCCACCGGTGAGCTGGCGAAAACGCAGCGCCGCAGCCTGGCGGATACCCTTCAGCAGTCCGGCGGCGAAGa                              <  2:152452/135‑1 (MQ=18)                                                                                                             cAGAACCCGCATCCAGACGATTTATCAGTCGGGGAGCTTCACGCCACTCATCAGGGTTGAAACCGCCACCGGTGAGCTGGCGAAAACGCAGCGCCGCAGCCTGGCGGATACCCTTCAGCAGTCCGGCGGCGAAGa                              >  1:152452/1‑135 (MQ=17)                                                                                                             cAGAACCCGCATCCAGACGATTTATCAGTCGGGGAGCTTCACGCCACTCATCAGGGTTGAAACCGCCACCGGTGAGCTGGCGAAAACGCAGCGCCGCAGCCTGGCGGATACCCTTCAGCAGTCCGGCGGCGAAGACGGTg                         >  2:308284/1‑140 (MQ=17)                                                                                                             cAGAACCCGCATCCAGACGATTTATCAGTCGGGGAGCTTCACGCCACTCATCAGGGTTGAAACCGCCACCGGTGAGCTGGCGAAAACGCAGCGCCGCAGCCTGGCGGATACCCTTCAGCAGTCCGGCGGCGAAGACGGTg                         >  2:388582/1‑140 (MQ=17)                                                                                                             cAGAACCCGCATCCAGACGATTTATCAGTCGGGGAGCTTCACGCCACTCATCAGGGTTGAAACCGCCACCGGTGAGCTGGCGAAAACGCAGCGCCGCAGCCTGGCGGATACCCTTCAGCAGTCCGGCGGCGAAGACGGTg                         >  2:479069/1‑140 (MQ=17)                                                                                                             cAGAACCCGCATCCAGACGATTTATCAGTCGGGGAGCTTCACGCCACTCATCAGGGTTGAAACCGCCACCGGTGAGCTGGCGAAAACGCAGCGCCGCAGCCTGGCGGATACCCTTCAGCAGTCCGGCGGCGAAGACGGTg                         >  2:543406/1‑140 (MQ=17)                                                                                                             cAGAACCCGCATCCAGACGATTTATCAGTCGGGGAGCTTCACGCCACTCATCAGGGTTGAAACCGCCACCGGTGAGCTGGCGAAAACGCAGCGCCGCAGCCTGGCGGATACCCTTCAGCAGTCCGGCGGCGAAGACGGTg                         >  1:530968/1‑140 (MQ=17)                                                                                                             cAGAACCCGCATCCAGACGATTTATCAGTCGGGGAGCTTCACGCCACTCATCAGGGTTGAAACCGCCACCGGTGAGCTGGCGAAAACGCAGCGCCGCAGCCTGGCGGATACCCTTCAGCAGTCCGGCGGCGAAGACGGTg                         >  1:422809/1‑140 (MQ=17)                                                                                                              aGAACCCGCATCCAGACGATTTATCAGTCGGGGAGCTTCACGCCACTCATCAGGGTTGAAACCGCCACCGGTGAGCTGGCGAAAACGCAGCGCCGCAGCCTGGCGGATACCCTTCAGCAGTCCGGCGGCGAAGACGGTgg                        <  2:497797/140‑1 (MQ=18)                                                                                                              aGAACCCGCATCCAGACGATTTATCAGTCGGGGAGCTTCACGCCACTCATCAGGGTTGAAACCGCCACCGGTGAGCTGGCGAAAACGCAGCGCCGCAGCCTGGCGGATACCCTTCAGCAGTCCGGCGGCGAAGACGGTgg                        <  2:867948/140‑1 (MQ=18)                                                                                                                aaCCCGCATCCAGACGATTTATCAGTCGGGGAGCTTCACGCCACTCATCAGGGTTGAAACCGCCACCGGTGAGCTGGCGAAAACGCAGCGCCGCAGCCTGGCGGATACCCTTCAGCAGTCCGGCGGCGAAGACGGTGGCa                      >  1:723253/1‑140 (MQ=17)                                                                                                                  cccGCATCCAGACGATTTATCAGTCGGGGAGCTTCACGCCACTCATCAGGGTTGAAACCGCCACCGGTGAGCTGGCGAAAACGCAGCGCCGCAGCCTGGCGGATACCCTTCAGCAGTCCGGCGGCGAAGACGGTGGCAgt                    <  2:422809/140‑1 (MQ=18)                                                                                                                                aTTTATCAGTCGGGGAGCTTCACGCCACTCATCAGGGTTGAAACCGCCACCGGTGAGCTGGCGAAAACGCAGCGCCGCAGCCTGGCGGATACCCTTCAGCAGTCCGGCGGCGAAGACGGTGGCAGTGTGGTGTTCccgcc      <  1:82915/140‑1 (MQ=17)                                                                                                                                    aTCAGTCGGGGAGCTTCACGCCACTCATCAGGGTTGAAACCGCCACCGGTGAGCTGGCGAAAACGCAGCGCCGCAGCCTGGCGGATACCCTTCAGCAGTCCGGCGGCGAAGACGGTGGCAGTGTGGTGTTCCCGCCGGTg  <  1:543406/140‑1 (MQ=17)                                                                                                                                         |                                                                                                                                       GGTGTGGCGACGTGAACGGGACCTGACGGGCTGGATGTCGCTGTCACGGAAACCGCAAGTGACCTGGTACGGCTGGGACGGCGACCGCCTGACCACGATACAGAACGACAGAACCCGCATCCAGACGATTTATCAGCCGGGGAGCTTCACGCCACTCATCAGGGTTGAAACCGCCACCGGTGAGCTGGCGAAAACGCAGCGCCGCAGCCTGGCGGATACCCTTCAGCAGTCCGGCGGCGAAGACGGTGGCAGTGTGGTGTTCCCGCCGGTG  >  NC_000913/734465‑734735

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