Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	NC_000913	3,767,673	T→C	pseudogene (453/702 nt)	yibJ →	putative RHS domain‑containing protein YibJ

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	3,767,673	0	T	C	100.0%	127.9 / NA	37	pseudogene (453/702 nt)	yibJ	putative RHS domain‑containing protein YibJ
Reads supporting (aligned to +/- strand):  ref base T (0/0);  new base C (20/17);  total (20/17)

CTGTATTACAACCGCCACCGCTATTATGACCCGCTGCTGGGGCGGTATATCACTCAGGATCCGATTGGACTGAAGGGGGGATGGAATTTTTATCAGTATCCGTTGAATCCGGTCATAAATGTAGATCCGCAAGGTTTGGTTGATATAAATTTATACCCCGAAAGTGATCTTATCCATTCTGTAGCTGATGAGATTAATATCCCAGGCGTTTTCACAATCGGGGGGCATGGTACCCCCACATCTATTGAATCCGCAACGCGCAGTATCAT  >  NC_000913/3767539‑3767807                                                                                                                                       |                                                                                                                                       cTGTATTACAACCGCCACCGCTATTATGACCCGCTGCTGGGGCGGTATATCACTCAGGATCCGATTGGACTGAAGGGGGGATGGAATTTTTATCAGTATCCGTTGAATCCGGTCATAAATGTAGATCCGCAAGGCTTGGt                                                                                                                                   <  1:128904/140‑1 (MQ=255)         cAACCGCCACCGCTATTATGACCCGCTGCTGGGGCGGTATATCACTCAGGATCCGATTGGACTGAAGGGGGGATGGAATTTTTATCAGTATCCGTTGAATCCGGTCATAAATGTAGATCCGCAAGGCTTGGTTGATATaa                                                                                                                           >  1:92901/1‑140 (MQ=255)         cAACCGCCACCGCTATTATGACCCGCTGCTGGGGCGGTATATCACTCAGGATCCGATTGGACTGAAGGGGGGATGGAATTTTTATCAGTATCCGTTGAATCCGGTCATAAATGTAGATCCGCAAGGCTTGGTTGATATaa                                                                                                                           >  2:163078/1‑140 (MQ=255)              gCCACCGCTATTATGACCCGCTGCTGGGGCGGTATATCACTCAGGATCCGATTGGACTGCAGGGGGGATGGAATTTTTATCAGTATCCGTTGAATCCGGTCATAAATGTAGATCCGCAAGGCTTGGTTGATATAAATTta                                                                                                                      >  2:245708/1‑140 (MQ=255)              gCCACCGCTATTATGACCCGCTGCTGGGGCGGTATATCACTCAGGATCCGATTGGACTGAAGGGGGGATGGAATTTTTATCAGTATCCGTTGAATCCGGTCATAAATGTAGATCCGCAAGGCTTGGTTGATATAAATTta                                                                                                                      >  2:317203/1‑140 (MQ=255)              gCCACCGCTATTATGACCCGCTGCTGGGGCGGTATATCACTCAGGATCCGATTGGACTGAAGGGGGGATGGAATTTTTATCAGTATCCGTTGAATCCGGTCATAAATGTAGATCCGCAAGGCTTGGTTGATATAAATTta                                                                                                                      >  1:272237/1‑140 (MQ=255)              gCCACCGCTATTATGACCCGCTGCTGGGGCGGTATATCACTCAGGATCCGATTGGACTGAAGGGGGGATGGAATTTTTATCAGTATCCGTTGAATCCGGTCATAAATGTAGATCCGCAAGGCTTGGTTGATATAAATTta                                                                                                                      >  2:158477/1‑140 (MQ=255)              gCCACCGCTATTATGACCCGCTGCTGGGGCGGTATATCACTCAGGATCCGATTGGACTGAAGGGGGGATGGAATTTTTATCA‑‑‑‑‑‑‑‑‑GAATCCGGTCATAAATGTAGATCCGCAAGGCTTGGTTGATATAAATTTATACCCCGaa                                                                                                             >  2:124833/1‑140 (MQ=17)                      tattatGACCCGCTGCTGGGGCGGTATATCACTCAGGATCCGATTGGACTGAAGGGGGGATGGAATTTTTATCAGTATCCGTTGAATCCGGTCATAAATGTAGATCCGCAAGGCTTGGTTGATATAAATTTATACCCCGa                                                                                                              <  2:92901/140‑1 (MQ=255)                         tatGACCCGCTGCTGGGGCGGTATATCACTCAGGATCCGATTGGACTGAAGGGGGGATGGAATTTTTATCAGTATCCGTTGAATCCGGTCATAAATGTAGATCCGCAAGGCTTGGTTGATATAAATTTATACCCCGAAAg                                                                                                           <  1:317203/140‑1 (MQ=255)                             aCCCGCTGCTGGGGCGGTATATCACTCAGGATCCGATTGGACTGAAGGGGGGATGGAATTTTTATCAGTATCCGTTGAATCCGGTCATAAATGTAGATCCGCAAGGCTTGGTTGATATAAATTTATACCCCGAAAGTGAt                                                                                                       >  1:154893/1‑140 (MQ=255)                                cGCTGCTGGGGCGGTATATCACTCAGGATCCGATTGGACTGAAGGGGGGATGGAATTTTTATCAGTATCCGTTGAATCCGGTCATAAATGTAGATCCGCAAGGCTTGGTTGATATAAATTTATACCCCGAAAGTGATCtt                                                                                                    >  1:170717/1‑140 (MQ=255)                                cGCTGCTGGGGCGGTATATCACTCAGGATCCGATTGGACTGAAGGGGGGATGGAATTTTTATCAGTATCCGTTGAATCCGGTCATAAATGTAGATCCGCAAGGCTTGGTTGATATAAATTTATACCCCGAAAGTGATCtt                                                                                                    >  1:55678/1‑140 (MQ=255)                                       ggggCGGTATATCACTCAGGATCCGATTGGACTGAAGGGGGGATGGAATTTTTATCA‑‑‑‑‑‑‑‑‑GAATCCGGTCATAAATGTAGATCCGCAAGGCTTGGTTGATATAAATTTATACCCCGAAAGTGATCTTATCCATTCTGTAGCtg                                                                                    <  1:124833/140‑1 (MQ=21)                                         ggCGGTATATCACTCAGGATCCGATTGGACTGAAGGGGGGATGGAATTTTTATCAGTATCCGTTGAATCCGGTCATAAATGTAGATCCGCAAGGCTTGGTTGATATAAATTTATACCCCGAAAGTGATCTTATCCATTCt                                                                                           >  1:87157/1‑140 (MQ=255)                                                tatCACTCAGGATCCGATTGGACTGAAGGGGGGATGGAATTTTTATCAGTATCCGTTGAATCCGGTCATAAATGTAGATCCGCAAGGCTTGGTTGATATAAATTTATACCCCGAAAGTGATCTTATCCATTCTGTAGCtg                                                                                    <  1:299880/140‑1 (MQ=255)                                                        aGGATCCGATTGGACTGAAGGGGGGATGGAATTTTTATCAGTATCCGTTGAATCCGGTCATAAATGTAGATCCGCAAGGCTTGGTTGATATAAATTTATACCCCGAAAGTGATCTTATCCATTCTGTAGCTGATGAGAtt                                                                            >  2:134007/1‑140 (MQ=255)                                                         ggATCCGATTGGACTGAAGGGGGGATGGAATTTTTATCAGTATCCGTTGAATCCGGTCATAAATGTAGATCCGCAAGGCTTGGTTGATATAAATTTATACCCCGAAAGTGATCTTATCCATTCTGTAGCTGATGAGATTa                                                                           <  1:182319/140‑1 (MQ=255)                                                                         aaGGGGGGATGGAATTTTTATCAGTATCCGTTGAATCCGGTCATAAATGTAGATCCGCAAGGCTTGGTTGATATAAATTTATACCCCGAAAGTGATCTTATCCATTCTGTAGCTGATGAGATTAATATCCCAGGCGtttt                                                           >  1:202779/1‑140 (MQ=255)                                                                          aGGGGGGATGGAATTTTTATCAGTATCCGTTGAATCCGGTCATAAATGTAGATCCGCAAGGCTTGGTTGATATAAATTTATACCCCGAAAGTGATCTTATCCATTCTGTAGCTGATGAGATTAATATCCCAGGCGTTTTc                                                          >  1:57676/1‑140 (MQ=255)                                                                              gggATGGAATTTTTATCAGTATCCGTTGAATCCGGTCATAAATGTAGATCCGCAAGGCTTGGTTGATATAAATTTATACCCCGAAAGTGATCTTATCCATTCTGTAGCTGATGAGATTAATATCCCAGGCGTTTTCACaa                                                      >  2:8395/1‑140 (MQ=255)                                                                               ggATGGAATTTTTATCAGTATCCGTTGAATCCTGTCATAAATGTAGATCCGCAAGGCTTGGTTGATATAAATTTATACCCCGAAAGTGATCTTATCCATTCTGTAGCTGATGAGATTAATATCCCAGGCGTTTTCACAAt                                                     >  2:109449/1‑140 (MQ=255)                                                                               ggATGGAATTTTTATCAGTATCCGTTGAATCCGGTCATAAATGTAGATCCGCAAGGCTTGGTTGATATAAATTTATACCCCGAAAGTGATCTTATCCATTCTGTAGCTGATGAGATTAATATCCCAGGCGTTTTCACAAt                                                     <  1:245708/140‑1 (MQ=255)                                                                                gATGGAATTTTTATCAGTATCCGTTGAATCCGGTCATAAATGTAGATCCGCAAGGCTTGGTTGATATAAATTTATACCCCGAAAGTGATCTTATCCATTCTGTAGCTGATGAGATTAATATCCCAGGCGTTTTCACAATc                                                    >  1:36565/1‑140 (MQ=255)                                                                                    gAATTTTTATCAGTATCCGTTGAATCCGGTCATAAATGTTGATCCGCAAGGCTTGGTTGATATAAATTTATACCCCGAAAGTGATCTTATCCATTCTGTAGCTGATGAGATTAATATCCCAGGCGTTTTCACAATCgggg                                                <  2:100464/140‑1 (MQ=255)                                                                                      aTTTTTATCAGTATCCGTTGAATCCGGTCATAAATGTAGATCCGCAAGGCTTGGTTGATATAAATTTATACCCCGAAAGTGATCTTATCCATTCTGTAGCTGATGAGATTAATATCCCAGGCGTTTTCACAATCgggggg                                              <  1:280814/140‑1 (MQ=255)                                                                                           tATCAGTATCCGTTGAATCCGGTCATAAATGTAGATCCGCAAGGCTTGGTTGATATAAATTTATACCCCGAAAGTGATCTTATCCATTCTGTAGCTGATGAGATTAATATCCCAGGCGTTTTCACAATCGGGGGGCATgg                                         <  2:170717/140‑1 (MQ=255)                                                                                           tATCAGTATCCGTTGAATCCGGTCATAAATGTAGATCCGCAAGGCTTGGTCGATATAAATCTATACCCCGAAAGTGATCTTATCCATTCTGTAGCTGATGAGATTAATATCCCAGGCGTTTTCACAATCGGGGGGCATgg                                         <  2:318045/140‑1 (MQ=255)                                                                                                 tATCCGTTGAATCCGGTCATAAATGTAGATCCGCAAGGCTTGGTTGATATAAATTTATACCCCGAAAGTGATCTTATCCATTCTGTAGCTGATGAGATTAATATCCCAGGCGTTTTCACAATCGGGGGGCATGGTAcccc                                   <  2:154893/140‑1 (MQ=255)                                                                                                                    tAAATGTAGATCCGCAAGGCTTGGTTGATATAAATTTATACCCCGAAAGTGATCTTATCCATTCTGTAGCTGATGAGATTAATATCCCAGGCGTTTTCACAATCGGGGGGCATGGTACCCCCACATCTATTGAATCCGCa                <  1:163078/140‑1 (MQ=255)                                                                                                                        tGTAGATCCGCAAGGCTTGGTTGATATAAATTTATACCCCGAAAGTGATCTTATCCATTCTGTAGCTGATGAGATTAATATCCCAGGCGTTTTCACAATCGGGGGGCATGGTACCCCCACATCTATTGAATCCGCAAcgc            <  2:202779/140‑1 (MQ=255)                                                                                                                          tAGATCCGCAAGGCTTGGTTGATATAAATTTATACCCCGAAAGTGATCTTATCCATTCTGTAGCTGATGAGATTAATATCCCAGGCGTTTTCACAATCGGGGGGCATGGTACCCCCACATCTATTGAATCCGCAACgcgc          >  2:261431/1‑140 (MQ=255)                                                                                                                          tAGATCCGCAAGGCTTGGTTGATATAAATTTATACCCCGAAAGTGATCTTATCCATTCTGTAGCTGATGAGATTAATATCCCAGGCGTTTTCACAATCGGGGGGCATGGTACCCCCACATCTATTGAATCCGCAACgcgc          <  2:57676/140‑1 (MQ=255)                                                                                                                          tAGATCCGCAAGGCTTGGTTGATATAAATTTATACCCCGAAAGTGATCTTATCCATTCTGTAGCTGATGAGATTAATATCCCAGGCGTTTTCACAATCGGGGGGCATGGTACCCCCACATCTATTGAATCCGCAACgcgc          <  1:322336/140‑1 (MQ=255)                                                                                                                             aTCCGCAAGGCTTGGTTGATATAAATTTATACCCCGAAAGTGATCTTATCCATTCTGTAGCTGATGAGATTAATATCCCAGGCGTTTTCACAATCGGGGGGCATGGTACCCCCACATCTATTGAATCCGCaa               >  1:339873/1‑132 (MQ=255)                                                                                                                             aTCCGCAAGGCTTGGTTGATATAAATTTATACCCCGAAAGTGATCTTATCCATTCTGTAGCTGATGAGATTAATATCCCAGGCGTTTTCACAATCGGGGGGCATGGTACCCCCACATCTATTGAATCCGCaa               <  2:339873/132‑1 (MQ=255)                                                                                                                                  cAAGGCTTGGTTGATATAAATTTATACCCCGAAAGTGATCTTATCCATTCTGTAGCTGATGAGATTAATATCCCAGGCGTTTTCACAATCGGGGGGCATGGTACCCCCACATCTATTGAATCCGCAACGCGCAGTATCAt  >  2:134949/1‑140 (MQ=255)                                                                                                                                       |                                                                                                                                       CTGTATTACAACCGCCACCGCTATTATGACCCGCTGCTGGGGCGGTATATCACTCAGGATCCGATTGGACTGAAGGGGGGATGGAATTTTTATCAGTATCCGTTGAATCCGGTCATAAATGTAGATCCGCAAGGTTTGGTTGATATAAATTTATACCCCGAAAGTGATCTTATCCATTCTGTAGCTGATGAGATTAATATCCCAGGCGTTTTCACAATCGGGGGGCATGGTACCCCCACATCTATTGAATCCGCAACGCGCAGTATCAT  >  NC_000913/3767539‑3767807

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