Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	NC_000913	3,651,423	C→T	pseudogene (133/746 nt)	yhiS →	putative uncharacterized protein YhiS

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	3,651,423	0	C	T	100.0%	100.3 / NA	30	pseudogene (133/746 nt)	yhiS	putative uncharacterized protein YhiS
Reads supporting (aligned to +/- strand):  ref base C (0/0);  new base T (18/12);  total (18/12)

CGCATTACATTTATTGTGTGTAGAAGAAAATCCATTCCTGTTCTAATTATCCAATTTAAACATCTTAGCATAAAACAAATGATGAATAAGGAATTTTATGTCAATTGACTTTACCCCAGGTATAATAAATACATATCACGGCGATATTTATAACTGCACAACAAATACCGATAATGCTAAAACGCCAGACACACCGAAGTGGCCTTGCGATAATTGGGAAGAACAACAACCCATTAATTCCACATTCTCTGGAGAAGGATATATCTCTGATCAATATGATTTAGCGCAACACCAACTTCAACAAATTAATGCGTGCCACACGAACACTACATATACAAATGCAGACTACTCAAAAGTTGTGGCTCAACTTGTAAGTCTTATTACTAACATTGAAACGATAAGTTCGACACAACTTACGCAACAGACACAATCAATACTTAATCAGATTAATAACATTAG  >  NC_000913/3651194‑3651652                                                                                                                                                                                                                                      |                                                                                                                                                                                                                                      cGCATTACATTTATTGTGTGTAGAAGAAAATCCATTCCTGTTCTAATTATCCAATTTAAACATCTTAGCATAAAACAAATGATGAATAAGGAATTTTATGTCAATTGACTTTACCCCAGGTATAATAAATACATATCACGGCGATATTTATAACTGCACAACAAATACCGATAATGCTAAAACGCCAGACACACCGAAGTGGCCTTGCGATAATTGGGAAGAACAACAATCCATTAATTCCACATTCTCTg                                                                                                                                                                                                                  >  2:61993/1‑251 (MQ=255)  gCATTACATTTATTGTGTGTAGAAGAAAATCCATTCCTGTTCTAATTATCCAATTTAAACATCTTAGCATAAAACAAATGATGAATAAGGAATTTTATGTCAATTGACTTTACCCCAGGTATAATAAATACATATCACGGCGATATTTATAACTGCACAACAAATACCGATAATGCTAAAACGCCAGACACACCGAAGTGGCCTTGCGATAATTGGGAAGAACAACAATCCATTAATTCCACATTCTCTgg                                                                                                                                                                                                                 <  1:331037/251‑1 (MQ=255)         aTTTATTGTGTGTAGAAGAAAATCCATTCCTGTTCTAATTATCCAATTTAAACATCTTAGCATAAAACAAATGATGAATAAGGAATTTTATGTCAATTGACTTTACCCCAGGTATAATAAATACATATCACGGCGATATTTATAACTGCACAACAAATACCGATAATGCTAAAACGCCAGACACACCGAAGTGGCCTTGCGATAATTGGGAAGAACAACAATCCATTAATTCCACATTCTCTGGAGAAGGa                                                                                                                                                                                                          <  1:232097/251‑1 (MQ=255)            tATTGTGTGTAGAAGAAAATCCATTCCTGTTCTAATTAACCAAGTTAAACATCTTAGCATAAAACAAATGATGAATAAGGAATTTTATGTCAATTGACTTTACCCCAGGTATAATAAATACATATCACAGCGATATTTATAACTGCACAACAAATGCCGATAATGCTAAAACGCCAGACACACCGAAGTGGCCTTGCGATAAGAAGGAAGAACAACAATCCATTAATTCCACATTCTCTGGAGAAGGata                                                                                                                                                                                                        <  2:160248/250‑1 (MQ=255)              ttGTGTGTAGAAGAAAATCCATTCCTGTTCTAATTATCCAATTTAAACATCTTAGCATAAAACAAATGATGAATAAGGAATTTTATGTCAATTGACTTTACCCCAGGTATAATAAATACATATCACGGCGATATTTATAACTGCACAACAAATACCGATAATGCTAAAACGCCAGACACACCGAAGTGGCCTTGCGATAATTGGGAAGAACAACAATCCATTAATTCCACATTCTCTGGAGAAGGatatat                                                                                                                                                                                                     >  2:163085/1‑251 (MQ=255)                 tgtgtAGAAGAAAATCCATTCCTGTTCTAATTATCCAATTTAAACATCTTAGCATAAAACAAATGATGAATAAGGAATTTTATGTCAATTGACTTTACCCCAGGTATAATAAATACATATCACGGCGATATTTATAACTGCACAACAAATACCGATAATGCTAAAACGCCAGACACACCGAAGTGGCCTTGCGATAATTGGGAAGAACAACAATCCATTAATTCCACATTCTCTGGAGAAGGATATAtctc                                                                                                                                                                                                  >  1:347408/1‑251 (MQ=255)                                   ttCCTGTTCTAATTATCCAATTTAAACATCTTAGCATAAAACAAATGATGAATAAGGAATTTTATGTCAATTGACTTTACCCCAGGTATAATAAATACATATCACGGCGATATTTATAACTGCACAACAAATACCGATAATGCTAAAACGCCAGACACACCGAAGTGGCCTTGCGATAATTGGGAAGAACAACAATCCATTAATTCCACATTCTCTGGAGAAGGATATATCTCTGATCAATATGATTTAgc                                                                                                                                                                                >  1:15941/1‑251 (MQ=255)                                     ccTGTTCTAATTATCCAATTTAAACATCTTAGCATAAAACAAATGATGAATAAGGAATTTTATGTCAATTGACTTTACCCCAGGTATAATAAATACATATCACGGCGATATTTATAACTGCACAACAAATACCGATAATGCTAAAACGCCAGACACACCGAAGTGGCCTTGCGATAATTGGGAAGAACAACAATCCATTAATTCCACATTCTCTGGAGAAGGATATATCTCTGATCAATATGATTTAgcgc                                                                                                                                                                              >  1:81641/1‑251 (MQ=255)                                                     aaTTTAAACATCTTAGCATAAAACAAATGATGAATAAGGAATTTTATGTCAATTGACTTTACCCCAGGTATAATAAATACATATCACGGCGATATTTATAACTGCACAACAAATACCGATAATGCTAAAACGCCAGACACACCGAAGTGGCCTTGCGATAATTGGGAAGAACAACAATCCATTAATTCCACATTCTCTGGAGAAGGATATATCTCTGATCAATATGATTTAGCGCAACACCAACTTcaac                                                                                                                                                               >  2:129400/1‑250 (MQ=255)                                                                cTTAGCATAAAACAAATGATGAATAAGGAATTTTATGTCAATTGACTTTACCCCAGGTATAATAAATACATATCACGGCGATATTTATAACTGCACAACAAATACCGATAATGCTAACACGCCAGACACACCGAAGTGGCCTTGCGATACTTGGGAAGAACAACAATCCATTAATTCCACATTCTCTGGAGACGGATATATCTCTGATCAATATGATTTAGCGCAACACCAACTTCAACAAATTAATGCGt                                                                                                                                                   >  1:215173/1‑251 (MQ=255)                                                                  tAGCATAAAACAAATGATGAATAAGGAATTTTATGTCAATTGACTTTACCCCAGGTATAATAAATACATATCACGGCGATATTTATAACTGCACAACAAATACCGATAATGCTAAAACGCCAGACACACCGAAGTGGCCTTGCGATAATTGGGAAGAACAACAATCCATTAATTCCACATTCTCTGGAGAAGGATATATCTCTGATCAATATGATTTAGCGCAACACCAACTTCAACAAATTAATGc                                                                                                                                                     <  1:4036/247‑1 (MQ=255)                                                                   aGCATAAAACAAATGATGAATAAGGAATTTTATGTCAATTGACTTTACCCCAGGTATAATAAATACATATCACGGCGATATTTATAACTGCACAACAAATACCGATAATGCTAAAACGCCAGACACACCGAAGTGGCCTTGCGATAATTGGGAAGAACAACAATCCATTAATTCCACATTCTCTGGAGAAGGATATATCTCTGATCAATATGATTTAGCGCAACACCAACTTCAACAAATTAAGGCGTGcc                                                                                                                                                <  2:15941/251‑1 (MQ=255)                                                                      aTAAAACAAATGATGAATAAGGAATTTTATGTCAATTGACTTTACCCCAGGTATAATAAATACATATCACGGCGATATTTATAACTGCACAACAAATACCGATAATGCTAAAACGCCAGACACACCGAAGTGGCCTTGCGATAATTGGGAAGAACAACAATCCATTAATTCCACATTCTCTGGAGAAGGATATATCTCTGATCAATATGATTTAGCGCAACACCAACTTCAACAAATTAATGCGTGCcaca                                                                                                                                             <  1:229736/251‑1 (MQ=255)                                                                        aaaacaccagatgaATCAGGATTTTCAAGTCAATTGACTTTACCCCAGGTATAATAAATACATATCACGGCGCTATTTATAACTGCATAACAACTACCGATAAAGCTAAAACGCCAGACACACCGAAGTGGACATGCGATAATTGGGAAGAACAACAATCCATTAATTCCACATTCTCTGGAGAAGGATATATCACTGATCAATATGATTTAGCGCAACACAAACTTCAACAAATTAATACGTGCcacccg                                                                                                                                           >  2:395974/10‑248 (MQ=255)                                                                        aaaaCAAATGATGAATAAGGAATTTTATGTCAATTGACTTTACCCCAGGTATAATAAATACATATCACGGCGATATTTATAACTGCACAACAAATACCGATAATGCTAAAACGCCAGACACACCGAAGTGGCCTTGCGATAATTGGGAAGAACAACAATCCATTAATTCCACATTCTCTGGAGAAGGATATATCTCTGATCAATATGATTTAGCGCAACACCAACTTCAACAAATTAATGCGTGCCacac                                                                                                                                            <  1:144843/250‑1 (MQ=255)                                                                            cAAATGAAGAATAAGGAATATTATGTCAATTGACTTTACCCCAGGAATAATAAATAGATATCACGGAGATATTTATAACTGCACAACAAATACAGATAATGCTACAACGCCAGAAACACCGAAGTGGACATGCGATAATTGCGGAGAACAACAATCCATCAATTCCACATTCTCTGGTGAAGGATATATCTCTGATCAATATGATTTAGCGCACCACCAACTTCAACAAATTAATGCGTGCCACACGAac                                                                                                                                        >  2:226747/1‑250 (MQ=255)                                                                                       tAAGGAATTTTATGTCAATTGACTTTACCCCAGGTATAATAAATACATATCACGGCGATATTTATAACTGCACAACAAATACCGATAATGCTAAAACGCCAGACACACCGAAGTGGCCTTGCGATAATTGGGAAGAACAACAATCCATTAATTCCACATTCTCTGGAGAAGGATATATCTCTGATCAATATGATTTAGCGCAACACAAACTTCAACAAATTAATGCGTGCCACACGAACACTACATATAc                                                                                                                             >  2:144331/1‑250 (MQ=255)                                                                                         aggcatcttATGTCAATTGGCTTTACCCCAGGTATAATAACTACATATAACGGCGATATTTATAACTGCACAACAAATACCGATAATGCTCACTCGCCAGACACACCGAAGTGGCCTTGCGATAATTGGGAAGAACAACCATCCATTACTTCCACATTCTCTGGCGAAGGATATATCTCTGATCAATATGATTTAGCGCAACACTAACTTCAACAAATTAATGCTTGCCACACGAACACTACATATACaa                                                                                                                           >  2:137886/8‑250 (MQ=255)                                                                                                                    ccAGGTATAATAAATACATATCACGGCGATATTTATAACTGCACAACAAATACCGATAATGCTAAAACGCCAGACACACCGAAGTGGCCTTGCGATAATTGGGAAGAACAACAATCCATTAATTCCACATTCTCTGGAGAAGGATATATCTCTGATCAATATGATTTAGCGCAACACCAACTTCAACAAATTAATGCGTGCCACACGAACACTACATATACAAATGCAGACTACTCAAAAGTTGTGGCTc                                                                                                >  1:68512/1‑250 (MQ=255)                                                                                                                     cAGGTATAATAAATACATATCACGGCGATATTTATAACTGCACAACAAATACCGATAATGCTAAAACGCCAGACACACCGAAGTGGCCTTGCGATAATTGGGAAGAACAACAATCCATTAATTCCACATTCTCTGGAGAAGGATATATCTCTGATCAATATGATTTAGCGCAACACCAACTTCAACAAATTAATGCGTGCCACACGAACACTACATATACAAATGCAGACTACTCAAAAGTTGTGGCTCa                                                                                               >  1:141477/1‑250 (MQ=255)                                                                                                                       ggTATACTAAATACCTATCACGGCGATATTTATAACTGCACAACACATACAGATAATGCTAAAACGCCAGACACACCGAAGCGGCCTTGCGATAATTGGGAAGAACAACAATCCATTAATTCCACATTCTCTGGCGAAGGATATATCTCTGATCAATATGATTTAGCGCAACACCAACTTCAACAACTTAATGCGTGCCACACGAACACTACATATACAAATGCAGACTACTCAAAAGCTGTGGCTCAACt                                                                                            >  1:315136/1‑251 (MQ=255)                                                                                                                        gTATAATAAATACATATCACGGCGATATTTATAACTGCACAACAAATACCGATAATGCTAAAACGCCAGACACACCGAAGTGGCCTTGCGATAATTGGGAAGAACAACAATCCATTAATTCCACATTCTCTGGAGAAGGATATATCTCTGATCAATATGATTTAGCGCAACACCAACTTCAACAAATTAATGCGTGCCACACGAACACTACATATACAAATGCAGACTACTCAAAAGTTGTGGCTCAACtt                                                                                           <  2:143166/251‑1 (MQ=255)                                                                                                                                           cGGCGATATTTATAACTGCACAACAAATACCGATAATGCTAAAACGCCAGACACACCGAAGTGGCCTTGCGATAATTGGGAAGAACAACAATCCATTAATTCCACATTCTCTGGAGAAGGATATATCTCTGATCAATATGATTTAGCGCAACACCAACTTCAACAAATTAATGCGTGCCACACGAACACTACATATACAAATGCAGACTACTCAAAAGTTGTGGCTCAACTTGTAAGTCTTATTACTAAc                                                                         >  1:212232/1‑250 (MQ=255)                                                                                                                                                                aacaaATACCGATAATGCTAAAACGCCAGACACACCGAAGTGGCCTTGCGATAATTGGGAAGAACAACAATCCATTAATTCCACATTCTCTGGAGAAGGATATATCTCTGATCAATATGATTTAGCGCAACACCAACTTCAACAAATTAATGCGTGCCACACGAACACTACATATACAAATGCAGACTACTCAAAAGTTGTGGCTCAACTTGTAAGTCTTATTACTAACATTGAAACGATAAGTTCGacac                                                   >  2:279053/1‑251 (MQ=255)                                                                                                                                                                   aaaTACCGATAATGCTAAAACGCCAGACACACCGAAGTGGCCTTGCGATAATTGGGAAGAACAACAATCCATTAATTCCACATTCTCTGGAGAAGGATATATCTCTGATCAATATGATTTAGCGCAACACCAACTTCAACAAATTAATGCGTGCCACACGAACACTACATATACAAATGCAGACTACTCAAAAGTTGTGGCTCAACTTGTAAGTCTTATTACTAACATTGAAACGATAAGTTCGACACAAc                                                >  2:202155/1‑251 (MQ=255)                                                                                                                                                                             aaTGCTAAAACGCCAGACACACCGAAGTGGCCTTGCGATAATTGGGAAGAACAACAATCCATTAATTCCACATTCTCTGGAGAAGGATATATCTCTGATCAATATGATTTAGCGCAACACCAACTTCAACAAATTAATGCGTGCCACACGAACACTACATATACAAATGCAGACTACTCAAAAGTTGTGGCTCAACTTGTAAGTCTTATTACTAACAt                                                                       <  1:64019/218‑1 (MQ=255)                                                                                                                                                                             aaTGCTAAAACGCCAGACACACCGAAGTGGCCTTGCGATAATTGGGAAGAACAACAATCCATTAATTCCACATTCTCTGGAGAAGGATATATCTCTGATCAATATGATTTAGCGCAACACCAACTTCAACAAATTAATGCGTGCCACACGAACACTACATATACAAATGCAGACTACTCAAAAGTTGTGGCTCAACTTGTAAGTCTTATTACTAACAt                                                                       >  2:64019/1‑218 (MQ=255)                                                                                                                                                                                                cacCGAAGTGGCCTTGCGATAATTGGGAAGAACAACAATCCATTAATTCCACATTCTCTGGAGAAGGATATATCTCTGATCAATATGATTTAGCGCAACACCAACTTCAACAAATTAATGCGTGCCACACGAACACTACATATACAAATGCAGACTACTCAAAAGTTGTGGCTCAACTTGTAAGTCTTATTACTAACATTGAAACGATAAGTTCGACACAACTTACGCAACAGACACAATCAATACTTAAt                   <  2:210860/251‑1 (MQ=255)                                                                                                                                                                                                                cGATAATTGGGAAGAACAACAATCCATTAATTCCACATTCTCTGGAGAAGGATATATCTCTGATCAATATGATTTAGCGCAACACCAACTTCAACAAATTAATGCGTGCCACACGAACACTACATATACAAATGCAGACTACTCAAAAGTTGTGGCTCAACTTGTAAGTCTTATTACTAACATTGAAACGATAAGTTCGACACAACTTACGCAACAGACACAATCAATACTTAATCAGATTAATAACATTa   <  2:191207/251‑1 (MQ=255)                                                                                                                                                                                                                 gATAATTGGGAAGAACAACAATCCATTAATTCCACATTCTCTGGAGAAGGATATATCTCTGATCAATATGATTTAGCGCAACACCAACTTCAACAAATTAATGCGTGCCACACGAACACTACATATACAAATGCAGACTACTCAAAAGTTGTGGCTCAACTTGTAAGTCTTATTACTAACATTGAAACGATAAGTTCGACACAACTTACGCAACAGACACAATCAATACTTAATCAGATTAATAACATTAg  <  1:260964/251‑1 (MQ=255)                                                                                                                                                                                                                                      |                                                                                                                                                                                                                                      CGCATTACATTTATTGTGTGTAGAAGAAAATCCATTCCTGTTCTAATTATCCAATTTAAACATCTTAGCATAAAACAAATGATGAATAAGGAATTTTATGTCAATTGACTTTACCCCAGGTATAATAAATACATATCACGGCGATATTTATAACTGCACAACAAATACCGATAATGCTAAAACGCCAGACACACCGAAGTGGCCTTGCGATAATTGGGAAGAACAACAACCCATTAATTCCACATTCTCTGGAGAAGGATATATCTCTGATCAATATGATTTAGCGCAACACCAACTTCAACAAATTAATGCGTGCCACACGAACACTACATATACAAATGCAGACTACTCAAAAGTTGTGGCTCAACTTGTAAGTCTTATTACTAACATTGAAACGATAAGTTCGACACAACTTACGCAACAGACACAATCAATACTTAATCAGATTAATAACATTAG  >  NC_000913/3651194‑3651652

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