Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	NC_000913	1,220,971	G→A	pseudogene (1371/2648 nt)	ycgH →	putative transporter component, N‑terminal fragment

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	1,220,971	0	G	A	100.0%	128.8 / NA	37	pseudogene (1371/2648 nt)	ycgH	putative transporter component, N‑terminal fragment
Reads supporting (aligned to +/- strand):  ref base G (0/0);  new base A (20/17);  total (20/17)

ACAGTATTATCAATCTTTCGTATGATGATGGTCAAACATTTACCCAAGGAAAAACATTAACCGTAAAAGGTAATTATGTCGGTAATAATGGTCAGCTTAATATCCGCACCGTATTAGGTGATGATAAATCGGCTACGGACAGACTTATTGTTGAGGGTAATACTTCGGGTTCAACTACCGTCTATGTGAAAAATGCTGGAGGAAGCGGCGCGGCCACGCTAAACGGGATCGAACTCATAACTGTGAATGGCGATGAATC  >  NC_000913/1220841‑1221099                                                                                                                                   |                                                                                                                                 aCAGTATTATCAATCTTTCGTATGATGATGGTCAAACATTTACCCAAGGAAAATCATTAACCGTAAAAGGTAATTATGTCGGTAATAATGGTCAGCTTAATATCCGCACCGTATTAGGTGATGATAAATCAGc                                                                                                                                >  2:284966/1‑133 (MQ=255)  cAGTATTATCAATCTTTCGTATGATGATGGTCAAACATTTACCCAAGGAAAAACATTAACCGTAAAAGGTAATTATGTCGGTAATAATGGTCAGCTTAATATCCGCACCGTATTAGGTGATGATAAATCAGCt                                                                                                                               <  2:512997/133‑1 (MQ=255)       ttatCAATCTTTCGTATGATGATGGTCAAACATTTACCCAAGGAAAAACATTAACCGTAAAAGGTAATTATGTCGGTAATAATGGTCAGCTTAATATCCGCACCGTATTAGGTGATGATAAATCAGCTACGGa                                                                                                                          <  2:374356/133‑1 (MQ=255)               cTTTCGTATGATGATGGTCAAACATTTACCCAAGGAAAAACATTAACCGTAAAAGGTAATTATGTCGGTAATAATGGTCAGCTTAATATCCGCACCGTATTAGGTGATGATAAATCAGCTACGGACAGACTTa                                                                                                                  >  2:87989/1‑133 (MQ=255)               cTTTCGTATGATGATGGTCAAACATTTACCCAAGGAAAAACATTAACCGTAAAAGGTAATTATGTCGGTAATAATGGTCAGCTTAATATCCGCACCGTATTAGGTGATGATAAATCAGCTACGGACAGACTTa                                                                                                                  >  2:58997/1‑133 (MQ=255)               cTTTCGTATGATGATGGTCAAACATTTACCCAAGGAAAAACATTAACCGTAAAAGGTAATTATGTCGGTAATAATGGTCAGCTTAATATCCGCACCGTATTAGGTGATGATAAATCAGCTACGGACAGACTTa                                                                                                                  <  1:393197/133‑1 (MQ=255)               cTTTCGTATGATGATGGTCAAACATTTACCCAAGGAAAAACATTAACCGTAAAAGGTAATTATGTCGGTAATAATGGTCAGCTTAATATCCGCACCGTATTAGGTGATGATAAATCAGCTACGGACAGACTTa                                                                                                                  >  2:393197/1‑133 (MQ=255)               cTTTCGTATGATGATGGTCAAACATTTACCCAAGGAAAAACATTAACCGTAAAAGGTAATTATGTCGGTAATAATGGTCAGCTTAATATCCGCACCGTATTAGGTGATGATAAATCAGCTACGGACAGACTTa                                                                                                                  <  1:4760/133‑1 (MQ=255)                 ttCGTATGATGATGGTCAAACATTTACCCAAGGAAAAACATTAACCGTAAAAGGTAATTATGTCGGTAATAATGGTCAGCTTAATATCCGCACCGTATTAGGTGATGATAAATCAGCTACGGACAGACTTAtt                                                                                                                >  1:19852/1‑133 (MQ=255)                  tCGTATGATGATGGTCAAACATTTACCCAAGGAAAAACATTAACCGTAAAAGGTAATTATGTCGGTAATAATGGTCAGCTTAATATCCGCACCGTATTAGGTGATGATAAATCAGCTACGGACAGACTTAttg                                                                                                               >  1:300792/1‑133 (MQ=255)                    gTATGATGATTGTCAGACATTTACCCAATGAAGAACATTAAGCGTAAAAGGTAATTATGTCGGTAATAATGGTCAGCTTAATATCCGCACCGTATTAGGTGATGATAAATCAGCTACGGACAGACTTAttgtt                                                                                                             <  2:185419/133‑1 (MQ=255)                    gTATGATGATGGTCAAACATTTACCCAAGGAAAAACATTAACCGTAAAAGGTAATTATGTCGGTAATAATGGTCAGCTTAATATCCGCACCGTATTAGGTGATGATAAATCAGCTACGGACAGACTTAttgtt                                                                                                             <  2:383780/133‑1 (MQ=255)                    gTATGATGATGGTCAAACATTTACCCAAGGAAAAACATTAACCGTAAAAGGTAATTATGTCGGTAATAATGGTCAGCTTAATATCCGCACCGTATTAGGTGATGATAAATCAGCTACGGACAGACTTAttgtt                                                                                                             <  1:32779/133‑1 (MQ=255)                    gTATGATGATGGTCAAACATTTACCCAAGGAAAAACATTAACCGTAAAAGGTAATTATGTCGGTAATAATGGTCAGCTTAATATCCGCACCGTATTAGGTGATGATAAATCAGCTACGGACAGACTTAttgtt                                                                                                             <  1:67877/133‑1 (MQ=255)                      atgatgatgGTCAAACATTTACCCAAGGAAAAACATTAACCGTAAAAGGTAATTATGTCGGTAATAATGGTCAGCTTAATATCCGCACCGTATTAGGTGATGATAAATCAGCTACGGACAGACTTATTGTTGa                                                                                                           <  1:87989/133‑1 (MQ=255)                                      aTTTACCCAAGGAAAAACATTAACCGTAAAAGGTAATTATGTCGGTAATAATGGTCAGCTTAATATCCGCACCGTATTAGGTGATGATAAATCAGCTACGGACAGACTTATTGTTGAGGGTAATACTTCGGGt                                                                                           <  2:277112/133‑1 (MQ=255)                                        ttACCCAAGGGAAAAGATTAACCGTAAAAGGTAATTATGTCGGTAATAATGGTCAGCTTAATATCCGCACCGTATTAGGTGATGATAAATCAGCTACGGACAGACTTATTGTTGAGGGTAATACTTCGGGTTc                                                                                         <  2:165882/133‑1 (MQ=255)                                           cccAAGGAAAAACATAAACCGTAAAATGGAATTATGTCGGTAATAATGGTCAGCTTAATATCCGCACCGTATTAGGTGATGATAAATCAGCTACGGACAGACTTATTGTTGAGGGTAATACTTCGGGTTCAAc                                                                                      <  2:227883/133‑1 (MQ=255)                                               aGGAAAAACATTAACCGTAAAAGGTAATTATGTCGGTAATAATGGTCAGCTTAATATCCGCACCGTATTAGGTGATGATAAATCAGCTACGGACAGACTTATAGTTGAGGATAATACATCGGGTACAACTAcc                                                                                  >  1:238856/1‑133 (MQ=255)                                                 gAAAAACATTAACCGTAAAAGGTAATTATGTCGGTAATAATGGTCAGCTTAATATCCGCACCGTATTAGGTGATGATAAATCAGCTACGGACAGACTTATTGTTGAGGGTAATACTTCGGGTTCAACTACCGt                                                                                <  2:1744/133‑1 (MQ=255)                                                          tAACCGTAAAAGGTAATTATGTCGGTAATAATGGTCAGCTTAATATCCGCACCGTATTAGGTGATGATAAATCAGCTACGGACAGACTTATTGTTGAGGGTAATACTTCGGGTTCAACTACCGTCTATGTGaa                                                                       >  1:74508/1‑133 (MQ=255)                                                                 aaaaGGTAATTATGTCGGTAATAATGGTCAGCTTAATATCCGCACCGTATTAGGTGATGATAAATCAGCTACGGACAGACTTATTGTTGAGGGTAATACTTCGGGTTCAACTACCGTCTATGTGAAAAATGCt                                                                >  2:321039/1‑133 (MQ=255)                                                                               tCGGTAATAATGGTCAGCTTAATATCCGCACCGTATTAGGTGATGATAAATCAGCTACGGACAGACTTATTGTTGAGGGTAATACTTCGGGTTCAACTACCGTCTATGTGAAAAATGCTGGAGGAAgcggcgc                                                  >  1:422704/1‑133 (MQ=255)                                                                                      taatGGTCAGCTTAATATCCGCACCGTACTAGGTGATGATAAATCAGCTACGGACAGACTTATTGTTGAGGGTAATACTTCGGGTTCAACTACCGTCTATGTGAAAAATGCTGGAGGAAGCGGCGCGGCCACg                                           <  2:300792/133‑1 (MQ=255)                                                                                         tGGTCAGCTTAATATCCGCACCGTATTAGGTGATGATAAATCAGCTACGGACAGACTTATTGTTGAGGGTAATACTTCGGGTTCAACTACCGTCTATGTGAAAAATGCTGGAGGAAGCGGCGCGGCCACGCTa                                        >  2:356265/1‑133 (MQ=255)                                                                                           gTCAGCTTAATATCCGCACCGTATTAGGTGATGATAAATCAGCTACGGACAGACTTATTGTTGAGGGTAATACTTCGGGTTCAACTACCGTCTATGTGAAAAATGCTGGAGGAAGCGGCGCGGCCACGCTaaa                                      <  1:58997/133‑1 (MQ=255)                                                                                               gCTTAATATCCGCACCGTATTAGGTGATGATAAATCAGCTACGGACAGACTTATTGTTGAGGGTAATACTTCGGGTTCAACTACCGTCTATGTGAAAAATGCTGGAGGAAGCGGCGCGGCCACGCTAAACggg                                  >  1:464788/1‑133 (MQ=255)                                                                                                 ttAATATCCGCACCGTATTAGGTGATGATAAATCAGCTACGGACAGACTTATTGTTGAGGGTAATACTTCGGGTTCAACTACCGTCTATGTGAAAAATGCTGGAGGAAGCGGCGCGGCCACGCTAAACGGGAt                                >  2:48661/1‑133 (MQ=255)                                                                                                     tatCCGCACCGTATTAGGTGATGATAAATCAGCTACGGACAGACTTATTGTTGAGGGTAATACTTCGGGTTCAACTACCGTCTATGTGAAAAATGCTGGAGGTAGCGGCGGGGCCACTCTAAACTGGATCGaa                            >  1:166278/1‑133 (MQ=255)                                                                                                              cGTATTAGGTGATGATAAATCAGCTACGGACAGACTTATTGTTGAGGGTAATACTTCGGGTTCAACTACCGTCTATGTGAAAAATGCTGGAGGAAGCGGCGCGGCCACGCTAAACGGGATCGAACTCATAACt                   >  1:442296/1‑133 (MQ=255)                                                                                                                    aGGTGATGATAAATCAGCTACGGACAGACTTATTGTTGAGGGTAATACTTCGGGTTCAACTACCGTCTATGTGAAAAATGCTGGAGGAAGCGGCGCGGCCACGCTAAACGGGATCGAACTCATAACTGTGAAt             <  1:321039/133‑1 (MQ=255)                                                                                                                        gatgatAAGTCAGCTACGGACAGACTTATTGTTGAGGGTAATACTTCGGGTTCAACTACCGTCTATGTGAAAAATGCTGGAGGAAGCGGCGCGGCCACGCTAAACGGGATCGAACTCATAACTGTGAATGGCg         >  1:495792/1‑133 (MQ=255)                                                                                                                        gatgatAAATCAGCTACGGACAGACTTATTGTTGAGGGTAATACTTCGGGTTCAACTACCGTCTATGTGAAAAATGCTGGAGGAAGCGGCGCGGCCACGCTAAACGGGATCGAACTCATAACTGTGAATGGCg         >  2:95695/1‑133 (MQ=255)                                                                                                                           gatAAATCAGCTACGGACAGACTTATTGTTGAGGGTAATACTTCGGGTTCAACTACCGTCTATGTGAAAAATGCTGGAGGAAGCGGCGCGGCCACGCTAAACGGGATCGAACTCATAACTGTGAATGGCGATg      >  2:294558/1‑133 (MQ=255)                                                                                                                           gatAAATCAGCTACGGACAGACTTATTGTTGAGGGTAATACTTCGGGTTCAACTACCGTCTATGTGAAAAATGCTGGAGGAAGCGGCGCGGCCACGCTAAACGGGATCGAACTCATAACTGTGAATGGCGATg      >  2:21983/1‑133 (MQ=255)                                                                                                                            atAAATCAGCCACGGACAGACTTATTGTTGAGCGTAATACTGCGGGTTCAACTACCGTCTATGTGAAAAATGCTGGAGGAAGCGGCGCGGCCACGCTAAACGGGATCGAACTCATAACTGTGAATGGCGATGa     <  2:218198/133‑1 (MQ=255)                                                                                                                               aaTCAGCTACGGACAGACTTATTGTTGAGGGTAATACTTCGGGTTCAACTACCGTCTATGTGAAAAATGCTGGAGGAAGCGGCGCGGCCACGCTAAACGGGAACGAACTCATAACTGAGAATGGCGATGAAtc  >  1:218198/1‑133 (MQ=255)                                                                                                                                   |                                                                                                                                 ACAGTATTATCAATCTTTCGTATGATGATGGTCAAACATTTACCCAAGGAAAAACATTAACCGTAAAAGGTAATTATGTCGGTAATAATGGTCAGCTTAATATCCGCACCGTATTAGGTGATGATAAATCGGCTACGGACAGACTTATTGTTGAGGGTAATACTTCGGGTTCAACTACCGTCTATGTGAAAAATGCTGGAGGAAGCGGCGCGGCCACGCTAAACGGGATCGAACTCATAACTGTGAATGGCGATGAATC  >  NC_000913/1220841‑1221099

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