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%	136.0 / NA	43	pseudogene (40/678 nt)	yedN	pseudogene, IpaH/YopM family
Reads supporting (aligned to +/- strand):  ref base C (0/0);  new base T (18/25);  total (18/25)

TTTTTTTTGAAGTCATTCGATGCGCCTTCTTTTACTGAATTAGCCCAGTCGCTGAAACCTTTGCAAAGTCTGGAATCTGTAATTTTCGTAAAGTTGTCTAGAAACTCTGCAATATTAGTCGATATAAACTCTAATTTATGCCGAATGGGTGGGTTACCAGTTCCCTGATTTTTTCTATCATGATAAAACCTATCTTATTATAAATCGAGAAATCTTGAAACCAGAGACTATAATCACCCCATTCGAAATATTCT  >  NC_000913/2011871‑2012124                                                                                                                                   |                                                                                                                            ttttttttGAAGTCATTCGATGCGCCTTCTTTTACTGAATTAGCCCAGTCGCTGAAACCTTTGCAAAGTCTGGAATCTGTAATTTTCGTAAAGTTGTCTAGAAACTCTGCAATATTAGTCGATATAAACTTTAATTTAt                                                                                                                     >  1:896645/1‑139 (MQ=255)           aGTCATTCGATGCGCCTTCTTTTACTGAATTAGCCCAGTCGCTGAAACCTTTGCAAAGTCTGGAATCTGTAATTTTCGTAAAGTTGTCTAGAAACTCTGCAATATTAGTCGATATAAACTTTAATTTATGCCGAAtggg                                                                                                           <  1:1525058/139‑1 (MQ=255)           aGTCATTCGATGCGCCTTCTTTTACTGAATTAGCCCAGTCGCTGAAACCTTTGCAAAGTCTGGAATCTGTAATTTTCGTAAAGTTGTCTAGAAACTCTGCAATATTAGTCGATATAAACTTTAATTTATGCCGAAtggg                                                                                                           <  2:574940/139‑1 (MQ=255)           aGTCATTCGATGCGCCTTCTTTTACTGAATTAGCCCAGTCGCTGAAACCTTTGCAAAGTCTGGAATCTGTAATTTTCGTAAAGTTGTCTAGAAACTCTGCAATATTAGTCGATATAAACTTTAATTTATGCCGAAtggg                                                                                                           <  2:1469968/139‑1 (MQ=255)           aGTCATTCGATGCGCCTTCTTTTACTGAATTAGCCCAGTCGCTGAAACCTTTGCAAAGTCTGGAATCTGTAATTTTCGTAAAGTTGTCTAGAAACTCTGCAATATTAGTCGATATAAACTTTAATTTATGCCGAAtggg                                                                                                           <  1:495420/139‑1 (MQ=255)                       cgcCTTCTTTTACTGAATTAGCCCAGTCGCTGAAACCTTTGCAAAGTCTGGAATCTGTAATTTTCGTAAAGTTGTCTAGAAACTCTGCAATATTAGTCGATATAAACTTTAATTTATGCCGAAtg                                                                                                             <  2:835204/125‑1 (MQ=255)                       cgcCTTCTTTTACTGAATTAGCCCAGTCGCTGAAACCTTTGCAAAGTCTGGAATCTGTAATTTTCGTAAAGTTGTCTAGAAACTCTGCAATATTAGTCGATATAAACTTTAATTTATGCCGAAtg                                                                                                             >  1:835204/1‑125 (MQ=255)                                      aaTTAGCCCAGTCGCTGAAACCTTTGCAAAGTCTGGAATCTGTAATTTTCGTAAAGTTGTCTAGAAACTCTGCAATATTAGTCGATATAATCTTTAATTTATGCCGAAtggg                                                                                                           <  2:1017343/112‑1 (MQ=255)                                      aaTTAGCCCAGTCGCTGAAACCTTTGCAAAGTCTGGAATCTGTAATTTTCGTAAAGTTGTCTAGAAACTCTGCAATATTAGTCGATATAAACTTTAATTTATGCCGAAtggg                                                                                                           >  1:1017343/1‑112 (MQ=255)                                      aaTTAGCCCAGTCGCTGAAACCTTTGCAAAGTCTGGAATCTGTAATTTTCGTAAAGTTGTCTAGAAACTCTGCAATATTAGTCGATATAAACTTTAATTTATGCCGAATGGGTGGGTTACCAGTTCCCTGATTTTTTCt                                                                                >  2:923397/1‑139 (MQ=255)                                         tAGCCCAGTCGCTGAAACCTTTGCAAAGTCTGGAATCTGTAATTTTCGTAAAGTTGTCTAGAAACTCTGCAATATTAGTCGATATAAACTTTAATTTATGCCGAATGGGTGGGTTACCAGTTCCCTGATTTTTTCTATc                                                                             <  1:420397/139‑1 (MQ=255)                                              cAGTCGCTGAAACCTTTGCAAAGTCTGGAATCTGTAATTTTCGTAAAGTTGTCTAGAAACTCTGCAATATTAGTCGATATAAACTTTAATTTATGCCGAATGGGTGGGTTACCAGTTCCCTGATTTTTTCTATCATGAt                                                                        <  1:923397/139‑1 (MQ=255)                                                gTCGCTGAAACCTTTGCAAAGTCTGGAATCTGTAATTTTCGTAAAGTTGTCTAGAAACTCTGCAATATTAGTCGATATAAACTTTAATTTATGCCGAAtggg                                                                                                           >  2:9111/1‑102 (MQ=255)                                                gTCGCTGAAACCTTTGCAAAGTCTGGAATCTGTAATTTTCGTAAAGTTGTCTAGAAACTCTGCAATATTAGTCGATATAAACTTTAATTTATGCCGAAtggg                                                                                                           <  1:9111/102‑1 (MQ=255)                                                gTCGCTGAAACCTTTGCAAAGTCTGGAATCTGTAATTTTCGTAAAGTTGTCTAGAAACTCTGCAATATTAGTCGATATAAACTTTAATTTATGCCGAAt                                                                                                              <  2:74892/99‑1 (MQ=255)                                                gTCGCTGAAACCTTTGCAAAGTCTGGAATCTGTAATTTTCGTAAAGTTGTCTAGAAACTCTGCAATATTAGTCGATATAAACTTTAATTTATGCCGAAt                                                                                                              >  1:74892/1‑99 (MQ=255)                                                gTCGCTGAAACCTTTGCAAAGTCTGGAATCTGTAATTTTCGTAAAGTTGTCTAGAAACTCTGCAATATTAGTCGATATAAACTTTAATTTATGCCGAATGGGTGGGTTACCAGTTCCCTGATTTTTTCTATc                                                                             >  1:361057/1‑132 (MQ=255)                                                gTCGCTGAAACCTTTGCAAAGTCTGGAATCTGTAATTTTCGTAAAGTTGTCTAGAAACTCTGCAATATTAGTCGATATAAACTTTAATTTATGCCGAATGGGTGGGTTACCAGTTCCCTGATTTTTTCTATc                                                                             <  2:361057/132‑1 (MQ=255)                                                              tgtaaaGTCTGGAATCTGTAATTTTCGTAAAGTTGTCTAGAAACTCTGCAATATTAGTCGATATAAACTTTAATTTATGCCGAATgggtgggt                                                                                                      <  1:297427/90‑1 (MQ=255)                                                              tgtaaaGTCTGGAATCTGTAATTTTCGTAAAGTTGTCTAGAAACTCTGCAATATTAGTCGATATAAACTTTAATTTATGCCGAATgggtgggt                                                                                                      >  2:297427/4‑93 (MQ=255)                                                                     tCTGGAATCTGTAATTTTCGTAAAGTTGTCTAGAAACTCTGCAATATTAGTCGATATAAACTTTAATTTATGCCGAAtggg                                                                                                           <  1:300201/81‑1 (MQ=255)                                                                     tCTGGAATCTGTAATTTTCGTAAAGTTGTCTAGAAACTCTGCAATATTAGTCGATATAAACTTTAATTTATGCCGAAtggg                                                                                                           >  2:300201/1‑81 (MQ=255)                                                                     tCTGGAATCTGTAATTTTCGTAAAGTTGTCTAGAAACTCTGCAATATTAGTCGATATAAACTTTAATTTATGCCGAATGGGTGGGTTACCAGTTCCCTGATTTTTTCTATCATGATAAAACCTATCTTATTATAAATCg                                                 >  2:1650833/1‑139 (MQ=255)                                                                         gAATCTGTAATTTTCGTAAAGTTGTCTAGAAACTCTGCAATATTAGTCGATATAAACTTTAATTTATGCCGAATGGGTGGGTTACCAGTTCCCTGATTTTTTCTATCATGATAAAACCTATCTTATTATAAATCGAGaa                                             <  1:1655018/139‑1 (MQ=255)                                                                               gTAATTTTCGTAAAGTTGTCTAGAAACTCTGCAATATTAGTCGATATAAACTTTAATTTATGCCGAAtggg                                                                                                           >  1:883329/1‑71 (MQ=255)                                                                               gTAATTTTCGTAAAGTTGTCTAGAAACTCTGCAATATTAGTCGATATAAACTTTAATTTATGCCGAAtggg                                                                                                           <  2:883329/71‑1 (MQ=255)                                                                               gTAATTTTCGTAAAGTTGTCTAGAAACTCTGCAATATTAGTCGATATAAACTTTAATTTATGCCGAATGGGTGGGTTACCAGTTCCCTGATTTTTTCTATCATGa                                                                         >  1:214129/1‑105 (MQ=255)                                                                               gTAATTTTCGTAAAGTTGTCTAGAAACTCTGCAATATTAGTCGATATAAACTTTAATTTATGCCGAATGGGTGGGTTACCAGTTCCCTGATTTTTTCTATCATGa                                                                         <  2:214129/105‑1 (MQ=255)                                                                               gTAATTTTCGTAAAGTTGTCTAGAAACTCTGCAATATTAGTCGATATAAACTTTAATTTATGCCGAATGGGTGGGTTACCAGTTCCCTGATTTTTTCTATCATGATAAAACCTa                                                                >  1:903122/1‑114 (MQ=255)                                                                               gTAATTTTCGTAAAGTTGTCTAGAAACTCTGCAATATTAGTCGATATAAACTTTAATTTATGCCGAATGGGTGGGTTACCAGTTCCCTGATTTTTTCTATCATGATAAAACCTa                                                                <  2:903122/114‑1 (MQ=255)                                                                               gTAATTTTCGTAAAGTTGTCTAGAAACTCTGCAATATTAGTCGATATAAACTTTAATTTATGCCGAATGGGTGGGTTACCAGTTCCCTGATTTTTTCTATCATGATAAAACCTATCTTATTATAAATCGAGAAATCTTg                                       >  1:914600/1‑139 (MQ=255)                                                                               gTAATTTTCGTAAAGTTGTCTAGAAACTCTGCAATATTAGTCGATATAAACTTTAATTTATGCCGAATGGGTGGGTTACCAGTTCCCTGATTTTTTCTATCATGATAAAACCTATCTTATTATAAATCGAGAAATCTTg                                       <  1:293182/139‑1 (MQ=255)                                                                                     ttCGTAAAGTTGTCTAGAAACTCTGCAATATTAGTCGATATAAACTTTAATTTATGCCGAATGGGTGGGTTACCAGTTCCCTGATTTTTTCTATCATGATAAAACCTATCTTATTATAAATCGAGAAATCTTGAAACCa                                 <  2:896645/139‑1 (MQ=255)                                                                                         tAAAGTTGTCTAGAAACTCTGCAATATTAGTCGATATAAACTTTAATTTATGCCGAATGGGTGGGTTACCAGTTCCCTGATTTTTTCTATCATGATAAAACCTATCTTATTATAAATCGAGAAATCTTGAAACCAgaga                             <  1:1650833/139‑1 (MQ=255)                                                                                                gTCTAGAAACTCTGCAATATTAGTCGATATAAACTTTAATTTATGCCGAAtgggtgg                                                                                                        <  1:1012656/57‑1 (MQ=255)                                                                                                gTCTAGAAACTCTGCAATATTAGTCGATATAAACTTTAATTTATGCCGAAtgggtgg                                                                                                        >  2:1012656/1‑57 (MQ=255)                                                                                                   tAGAAACTCTGCAATATTAGTCGATATAAACTTTAATTTATGCCGAATGGGTGGGTTACCAGTTCCCTGATTTTTTCTATCATGATAAAACCTATCTTATTATaaa                                                    >  2:746090/1‑106 (MQ=255)                                                                                                   tAGAAACTCTGCAATATTAGTCGATATAAACTTTAATTTATGCCGAATGGGTGGGTTACCAGTTCCCTGATTTTTTCTATCATGATAAAACCTATCTTATTATaaa                                                    <  1:746090/106‑1 (MQ=255)                                                                                                           ctGCAATATTAGTCGATATAAACTTTAATTTATGCCGaa                                                                                                               >  2:549935/1‑39 (MQ=255)                                                                                                           ctGCAATATTAGTCGATATAAACTTTAATTTATGCCGaa                                                                                                               <  1:549935/39‑1 (MQ=255)                                                                                                            tGCAATATTAGTCGATATAAACTTTAATTTATGCCGAATGGGTGGGTTACCAGTTCCCTGATTTTTTCTATCATGATAAAACCTATCTTATTATAAATCGAGAAATCTTGAAACCag                                <  1:1546487/117‑1 (MQ=255)                                                                                                            tGCAATATTAGTCGATATAAACTTTAATTTATGCCGAATGGGTGGGTTACCAGTTCCCTGATTTTTTCTATCATGATAAAACCTATCTTATTATAAATCGAGAAATCTTGAAACCag                                >  2:1546487/1‑117 (MQ=255)                                                                                                                    tAGTCGATATAAACTTTAATTTATGCCGAATGGGTGGGTTACCAGTTCCCTGATTTTTTCTATCATGATAAAACCTATCTTATTATAAATCGAGAAATCTTGAAACCAGAGACTATAATCACCCCATTCGAAATATTCt  <  2:914600/139‑1 (MQ=255)                                                                                                                                   |                                                                                                                            TTTTTTTTGAAGTCATTCGATGCGCCTTCTTTTACTGAATTAGCCCAGTCGCTGAAACCTTTGCAAAGTCTGGAATCTGTAATTTTCGTAAAGTTGTCTAGAAACTCTGCAATATTAGTCGATATAAACTCTAATTTATGCCGAATGGGTGGGTTACCAGTTCCCTGATTTTTTCTATCATGATAAAACCTATCTTATTATAAATCGAGAAATCTTGAAACCAGAGACTATAATCACCCCATTCGAAATATTCT  >  NC_000913/2011871‑2012124

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