Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	NC_000913	2,012,001	C→T	pseudogene (40/678 nt)	yedN ←	pseudogene, IpaH/YopM family

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	2,012,001	0	C	T	100.0%	129.1 / NA	41	pseudogene (40/678 nt)	yedN	pseudogene, IpaH/YopM family
Reads supporting (aligned to +/- strand):  ref base C (0/0);  new base T (19/22);  total (19/22)

ACATTTTTTTTGAAGTCATTCGATGCGCCTTCTTTTACTGAATTAGCCCAGTCGCTGAAACCTTTGCAAAGTCTGGAATCTGTAATTTTCGTAAAGTTGTCTAGAAACTCTGCAATATTAGTCGATATAAACTCTAATTTATGCCGAATGGGTGGGTTACCAGTTCCCTGATTTTTTCTATCATGATAAAACCTATCTTATTATAAATCGAGAAATCTTGAAACCAGAGACTATAATCACCCCATTCGAAATATTCTGAGTGCA  >  NC_000913/2011868‑2012131                                                                                                                                      |                                                                                                                                   acaTTTTTTTTGAAGTCATTCGATGCGCCTTCTTTTACTGAATTAGCCCAGTCGCTGAAACCTTTGCAAAGTCTGGAATCTGTAATTTTCGTAAAGTTGTCTAGAAACTCTGCAATATTAGTCGATATAAACTTTAAtt                                                                                                                               >  2:123215/1‑139 (MQ=255)   aTTTTTTTTGAAGTCATTCGATGCGCCTTCTTTTACTGAATTAGCCCAGTCGCTGAAACCTTTGCAAAGTCTGGAATCTGTAATTTTCGTAAAGTTGTCTAGAAACTCTGCAATATTAGTCGATATAAACTTTAATTTa                                                                                                                             >  1:90244/1‑139 (MQ=255)   aTTTTTTTTGAAGTCATTCGATGCGCCTTCTTTTACTGAATTAGCCCAGTCGCTGAAACCTTTGCAAAGTCTGGAATCTGTAATTTTCGTAAAGTTGTCTAGAAACTCTGCAATATTAGTCGATATAAACTTTAATTTa                                                                                                                             >  1:610535/1‑139 (MQ=255)          ttGAAGTCATTCGATGCGCCTTCTTTTACTGGATTAGCCCAGTCGCTGAAACCTTTGCAAAGTCTGGAATCTGTAATTTTCGTAAAGTTGTCTAGAAACTCTGCAATATTAGTCGATATAAACTTTAATTTATGCCGaa                                                                                                                      <  2:980664/139‑1 (MQ=255)              aGTCATTCGATGCGCCTTCTTTTACTGAATTAGCCCAGTCGCTGAAACCTTTGCAAAGTCTGGAATCTGTAATTTTCGTAAAGTTGTCTAGAAACTCTGCAATATTAGTCGATATAAACTTTAATTTATGCCGAAtggg                                                                                                                  <  1:1410949/139‑1 (MQ=255)               gTCATTCGATGCGCCTTCTTTTACTGAATTAGCCCAGTCGCTGAAACCTTTGCAAAGTCTGGAATCTGTAATTTTCGTAAAGTTGTCTAGAAACTCTGCAATATTAGTCGATATAAACTTTAATTTATGCCGAAtgggt                                                                                                                 <  1:674405/139‑1 (MQ=255)                             cttcttTTACTGAATTAGCCCAGTCGCTGAAACCTTTGCAAAGTCTGGAATCTGTAATTTTCGTAAAGTTGTCTAGAAACTCTGCAATATTAGTCGATATAAACTTTAATTTATGCCGAAtgggtg                                                                                                                >  2:1020882/1‑126 (MQ=255)                             cttcttTTACTGAATTAGCCCAGTCGCTGAAACCTTTGCAAAGTCTGGAATCTGTAATTTTCGTAAAGTTGTCTAGAAACTCTGCAATATTAGTCGATATAAACTTTAATTTATGCCGAAtgggtg                                                                                                                <  1:1020882/126‑1 (MQ=255)                                           ttAGCCCAGTCGCTGAAACCTTTGCAAAGTCTGGAATCTGTAATTTTCGTAAAGTTGTCTAGAAACTCTGCAATATTAGTCGATATAAACTTTAATTTATGCCGAATGGGTGGGTTACCAGTTCCCTGATTTTTTCTAt                                                                                     >  2:880260/1‑139 (MQ=255)                                            tAGCCCAGTCGCTGAAACCTTTGCAAAGTCTGGAATCTGTAATTTTCGTAAAGTTGTCTAGAAACTCTGCAATATTAGTCGATATAAACTTTAATTTATGCCGAATGGGTGGGTTACCAGTTCCCTGATTTTTTCTATc                                                                                    <  2:90244/139‑1 (MQ=255)                                               cccAGTCGCTGAAACCTTTGCAAAGTCTGGAATCTGTAATTTTCGTAAAGTTGTCTAGAAACTCTGCAATATTAGTCGATATAAACTTTAATTTATGCCGAATGGGTGGGTTACCAGTTCCCTGATTTTTTCTATCATg                                                                                 >  1:876729/1‑139 (MQ=255)                                                ccAGTCGCTGAAACCTTTGCAAAGTCTGGAATCTGTAATTTTCGTAAAGTTGTCTAGAAACTCTGCAATATTAGTCGATATAAACTTTAATTTATGCCGAATGGGTGGGTTACCAGTTCCCTGATTTTTTCTATCATGa                                                                                <  1:880260/139‑1 (MQ=255)                                                 cAGTCGCTGAAACCTTTGCAAATTCTGGAATCTGTAATTTTCGTAAAGTTGTCTAGAAACTCTGCAATATTAGTCGATATAAACTTTAATTTATGCCGAATGGGTGGGTTACCAGTTCCCTGATTTTTTCTATCATGAt                                                                               <  1:461906/139‑1 (MQ=255)                                                   gTCGCTGAAACCTTTGCAAAGTCTGGAATCTGTAATTTTCGTAAAGTTGTCTAGAAACTCTGCAATATTAGTCGATATAAACTTTAATTTATGCCGAAtggg                                                                                                                  <  2:602983/102‑1 (MQ=255)                                                   gTCGCTGAAACCTTTGCAAAGTCTGGAATCTGTAATTTTCGTAAAGTTGTCTAGAAACTCTGCAATATTAGTCGATATAAACTTTAATTTATGCCGAAtggg                                                                                                                  >  1:602983/1‑102 (MQ=255)                                                      gCTGAAACCTTTGCAAAGTCTGGAATCTGTAATTTTCGTAAAGTTGTCTAGAAACTCTGCAATATTAGTCGATATAAACTTTAATTTATGCCGAATGGGTGGGTTACCAGTTCCCTGATTTTTTCTATCATGATAAAAc                                                                          <  2:610535/139‑1 (MQ=255)                                                      gCTGAAACCTTTGCAAAGTCTGGAATCTGTAATTTTCGTAAAGTTGTCTAGAAACTCTGCAATATTAGTCGATATAAACTTTAATTTATGCCGAATGGGTGGGTTACCAGTTCCCTGATTTTTTCTATCATGATAAAAc                                                                          <  1:737622/139‑1 (MQ=255)                                                                ttGCAAAGTCTGGAATCTGTAATTTTCGTAAAGTTGTCTAGAAACTCTGCAATATTAGTCGATATAAACTTTAATTTATGCCGAATGGGTGGGTTACCAGTTCCCTGATTTTTTCTATCATGa                                                                                >  2:1352670/1‑123 (MQ=255)                                                                ttGCAAAGTCTGGAATCTGTAATTTTCGTAAAGTTGTCTAGAAACTCTGCAATATTAGTCGATATAAACTTTAATTTATGCCGAATGGGTGGGTTACCAGTTCCCTGATTTTTTCTATCATGa                                                                                <  1:1352670/123‑1 (MQ=255)                                                                ttGCAAAGTCTGGAATCTGTAATTTTCGTAAAGTTGTCTAGAAACTCTGCAATATTAGTCGATATAAACTTTAATTTATGCCGAATGGGTGGGTTACCAGTTCCCTGATTTTTTCTATCATGAt                                                                               <  1:10774/124‑1 (MQ=255)                                                                ttGCAAAGTCTGGAATCTGTAATTTTCGTAAAGTTGTCTAGAAACTCTGCAATATTAGTCGATATAAACTTTAATTTATGCCGAATGGGTGGGTTACCAGTTCCCTGATTTTTTCTATCATGAt                                                                               >  2:10774/1‑124 (MQ=255)                                                                 tGCAAAGTCTGGAATCTGTAATTTTCGTAAAGTTGTCTAGAAACTCTGCAATATTAGTCGATATAAACTTTAATTTATGCCGAATGGGTGGGTTACCAGTTCCCTGATTTTTTCTATCATGAt                                                                               >  1:1060090/1‑123 (MQ=255)                                                                 tGCAAAGTCTGGAATCTGTAATTTTCGTAAAGTTGTCTAGAAACTCTGCAATATTAGTCGATATAAACTTTAATTTATGCCGAATGGGTGGGTTACCAGTTCCCTGATTTTTTCTATCATGAt                                                                               <  2:1060090/123‑1 (MQ=255)                                                                                  gTAATTTTCGTAAAGTTGTCTAGAAACTCTGCAATATTAGTCGATATAAACTTTAATTTATGCCGAATGGGTGGGTTAc                                                                                                          >  1:1349607/1‑79 (MQ=255)                                                                                  gTAATTTTCGTAAAGTTGTCTAGAAACTCTGCAATATTAGTCGATATAAACTTTAATTTATGCCGAATGGGTGGGTTAc                                                                                                          <  2:1349607/79‑1 (MQ=255)                                                                                  gTAATTTTCGTAAAGTTGTCTAGAAACTCTGCAATATTAGTCGATATAAACTTTAATTTATGCCGAATGGGTGGGTTACCAGTTCCCTGATTTTTTCTATCATGAt                                                                               <  1:508386/106‑1 (MQ=255)                                                                                  gTAATTTTCGTAAAGTTGTCTAGAAACTCTGCAATATTAGTCGATATAAACTTTAATTTATGCCGAATGGGTGGGTTACCAGTTCCCTGATTTTTTCTATCATGAt                                                                               >  2:508386/1‑106 (MQ=255)                                                                                          cGTAAAGTTGTCTAGAAACTCTGCAATATTAGTCGATATAAACTTTAATTTATGCCGAATGGGTGGGTTACCAGTTCCCTGATTTTTTCTATc                                                                                    >  1:462025/1‑93 (MQ=255)                                                                                          cGTAAAGTTGTCTAGAAACTCTGCAATATTAGTCGATATAAACTTTAATTTATGCCGAATGGGTGGGTTACCAGTTCCCTGATTTTTTCTATc                                                                                    <  2:462025/93‑1 (MQ=255)                                                                                                   gTCTAGAAACTCTGCAATATTAGTCGATATAAACTTTAATTTATGCCGAATGGGTGGGTTACCAGTTCCCTGATTTTTTCTATCATGATa                                                                              <  1:195756/90‑1 (MQ=255)                                                                                                   gTCTAGAAACTCTGCAATATTAGTCGATATAAACTTTAATTTATGCCGAATGGGTGGGTTACCAGTTCCCTGATTTTTTCTATCATGATa                                                                              >  2:195756/1‑90 (MQ=255)                                                                                                              ctGCAATATTAGTCGATATAAACTTTAATTTATGCCGAATGGGTGGGTTACCAGTTCCCTGATTTTTTCTATCATGATAAAAcc                                                                         <  2:1448682/84‑1 (MQ=255)                                                                                                              ctGCAATATTAGTCGATATAAACTTTAATTTATGCCGAATGGGTGGGTTACCAGTTCCCTGATTTTTTCTATCATGATAAAAcc                                                                         >  1:1448682/1‑84 (MQ=255)                                                                                                               tGCAATATTAGTCGATATAAACTTTAATTTATGCCGAATGGGTGGGTTAc                                                                                                          <  2:447172/50‑1 (MQ=255)                                                                                                               tGCAATATTAGTCGATATAAACTTTAATTTATGCCGAATGGGTGGGTTAc                                                                                                          >  1:447172/1‑50 (MQ=255)                                                                                                                  aaTATTAGTCGATATAAACTTTAATTTATGCCGAATGGGTGGGTTACCAGTTCCCTGATTTTTTCTATCATGATAAAACCTATCTTATTATAAATCGAGAAATCTTGAAACCAGAGACTATAATCACCCCATTCGAAat              <  2:934006/139‑1 (MQ=255)                                                                                                                          tCGATATAAACTTTAATTTATGCCGAATGGGTGGGTTACCAGTTCCCTGATTTTTTCTATCATGATAAAACCTATCTTATTATAAATCGAGAAATCTTGAAACCAGAGACTATAATCACCCCATTCGAAATATTCTGAg      <  1:123215/139‑1 (MQ=255)                                                                                                                              tataAACTTTAATTTATGCCGAATGGGTGGGTTACCAGTTCCCTGATTTTTTCTATCATGATAAAACCTATCTTATTATAAATCGAGAAATCTTGAAACCAGAGACTATAATCACCCCATTCGAAATATTCTGAGTGCa  >  1:54658/1‑139 (MQ=255)                                                                                                                              tataAACTTTAATTTATGCCGAATGGGTGGGTTACCAGTTCCCTGATTTTTTCTATCATGATAAAACCTATCTTATTATAAATCGAGAAATCTTGAAACCAGAGACTATAATCACCCCATTCGAAATATTCTGAGTGCa  >  1:168306/1‑139 (MQ=255)                                                                                                                              tataAACTTTAATTTATGCCGAATGGGAGGGTTACCAGTTCCCTGATTTTTTCTATCATGATAAAACCTATCTTATTATAAATCGAGAAATCTTGAAACCag                                       <  2:182488/102‑1 (MQ=255)                                                                                                                              tataAACTTTAATTTATGCCGAATGGGAGGGTTACCAGTTCCCTGATTTTTTCTATCATGATAAAACCTATCTTATTATAAATCGAGAAATCTTGAAACCag                                       >  1:182488/1‑102 (MQ=255)                                                                                                                                      |                                                                                                                                   ACATTTTTTTTGAAGTCATTCGATGCGCCTTCTTTTACTGAATTAGCCCAGTCGCTGAAACCTTTGCAAAGTCTGGAATCTGTAATTTTCGTAAAGTTGTCTAGAAACTCTGCAATATTAGTCGATATAAACTCTAATTTATGCCGAATGGGTGGGTTACCAGTTCCCTGATTTTTTCTATCATGATAAAACCTATCTTATTATAAATCGAGAAATCTTGAAACCAGAGACTATAATCACCCCATTCGAAATATTCTGAGTGCA  >  NC_000913/2011868‑2012131

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