Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	NC_000913	315,061	G→T	pseudogene (705/872 nt)	eaeH →	pseudogene, attaching and effacing protein homology,factor, Not classified

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	315,061	0	G	T	100.0%	182.7 / NA	52	pseudogene (705/872 nt)	eaeH	pseudogene, attaching and effacing protein homology,factor, Not classified
Reads supporting (aligned to +/- strand):  ref base G (0/0);  new base T (27/25);  total (27/25)

TCCCGTAGTCATACCCGCATTGGTGTTGGTGCGGAATACTGGCGCGATTATCTGAAACTGAGCGCCAATGGTTATATTCGGGCTTCTGGCTGGAAAAAATCGCCGGATATTGAGGATTATCAGGAACGCCCGGCGAATGGCTGGGATATTCGTGCTGAGGGCTATTTACCCGCCTGGCCGCAGCTTGGCGCAAGCCTGATGTATGAACAGTATTATGGCGATGAAGTCGGGCTGTTTGGTAAAGATAAGCGCCAGAAAGACCCGCATG  >  NC_000913/314927‑315194                                                                                                                                       |                                                                                                                                      tCCCGTAGTCATACCCGCATTGGTGTTGGTGCGGAATACTGGCGCGATTATCTGAAACTGAGCGCCAATGGTTATATTCGGGCTTCTGGCTGGAAAAAATCGCCGGATATTGAGGATTATCAGGAACGCCCGGCTAATg                                                                                                                                   >  1:188002/1‑139 (MQ=255)     gTAGTCATACCCGCATTGGTGTTGGTGCGGAATACTGGCGCGATTATCTGAAACTGAGCGCCAATGGTTATATTCGGGCTTCTGGCTGGAAAAAATCGCCGGATATTGAGGATTATCAGGAACGCCCGGCTAATGGCTg                                                                                                                               >  2:374492/1‑139 (MQ=255)      tAGTCATACCCGCATTGGTGTTGGTGCGGAATACTGGCGCGATTATCTGAAACTGAGCGCCAATGGTTATATTCGGGCTTCTGGCTGGAAAAAATCGCCGGATATTGAGGATTATCAGGAACGCCCGGCTAATGGCTgg                                                                                                                              <  2:904595/139‑1 (MQ=255)      tAGTCATACCCGCATTGGTGTTGGTGCGGAATACTGGCGCGATTATCTGAAACTGAGCGCCAATGGTTATATTCGGGCTTCTGGCTGGAAAAAATCGCCGGATATTGAGGATTATCAGGAACGCCCGGCTAATGGCTgg                                                                                                                              <  2:188002/139‑1 (MQ=255)      tAGTCATACCCGCATTGGTGTTGGTGCGGAATACTGGCGCGATTATCTGAAACTGAGCGCCAATGGTTATATTCGGGCTTCTGGCTGGAAAAAATCGCCGGATATTGAGGATTATCAGGAACGCCCGGCTAATGGCTgg                                                                                                                              >  1:942165/1‑139 (MQ=255)      tAGTCATACCCGCATTGGTGTTGGTGCGGAATACTGGCGCGATTATCTGAAACTGAGCGCCAATGGTTATATTCGGGCTTCTGGCTGGAAAAAATCGCCGGATATTGAGGATTATCAGGAACGCCCGGCTAATGGCTgg                                                                                                                              <  1:803262/139‑1 (MQ=255)      tAGTCATACCCGCATTGGTGTTGGTGCGGAATACTGGCGCGATTATCTGAAACTGAGCGCCAATGGTTATATTCGGGCTTCTGGCTGGAAAAAATCGCCGGATATTGAGGATTATCAGGAACGCCCGGCTAATGGCTgg                                                                                                                              >  1:671250/1‑139 (MQ=255)             aCCCGCATTGGTGTTGGTGCGGAATACTGGCGCGATTATCTGAAACTGAGCGCCAATGGTTATATTCGGGCTTCTGGCTGGAAAAAATCGCCGGATATTGAGGATTATCAGGAACGCCCGGCTAATGGCTGGGATATTc                                                                                                                       >  1:309255/1‑139 (MQ=255)                   aTTGGTGTTGGTGCGGAATACTGGCGCGATTATCTGAAACTGAGCGCCAATGGTTATATTCGGGCTTCTGGCTGGAAAAAATCGCCGGATATTGAGGATTATCAGGAACGCCCGGCTAATGGCTGGGATATTCGTGCTg                                                                                                                 >  2:776283/1‑139 (MQ=255)                            ggtgCGGAATACTGGCGCGATTATCTGAAACTGAGCGCCAATGGTTATATTCGGGCTTCTGGCTGGAAAAAATCGCCGGATATTGACGATTATCAGGAACGCCCGGCTAATGGCTGGGATATTCGTGCTGAGGGCTAtt                                                                                                        >  2:537711/1‑139 (MQ=255)                                               aTTATCTGAAACTGAGCGCCAATGGTTATATTCGGGCTTCTGGCTGGAAAAAATCGCCGGATATTGAGGATTATCAGGAACGCCCGGCTAATGGCTGGGATATTCGTGCTGAGGGCTATTTACCCGCCTGGCCGCAGCt                                                                                     <  1:374492/139‑1 (MQ=255)                                                         aCTGAGCGCCAATGGTTATATTCGGGCTTCTGGCTGGAAAAAATCGCCGGATATTGAGGATTATCAGGAACGCCCGGCTAATGGCTGGGATATTCGTGCTGAGGGCTATTTACCCGCCTGGCCGCAGCTTGGCGCAAGc                                                                           <  1:184011/139‑1 (MQ=255)                                                         aCTGAGCGCCAATGGTTATATTCGGGCTTCTGGCTGGAAAAAATCGCCGGATATTGAGGATTATCAGGAACGCCCGGCTAATGGCTGGGATATTCGTGCTGAGGGCTATTTACCCGCCTGGCCGCAGCTTGGCGCAAGc                                                                           >  2:184011/1‑139 (MQ=255)                                                         aCTGAGCGCCAAAGGTTATATTCGGGCTTCTGGCTGGAAAAAATCGCCGGATATTGAGGATTATCAGGAACGCCCGGCTAATGGCTGGGATATTCGTGCTGAGGGCTATTTACCCGCCTGGCCGCAGCTTGGCGCAAGc                                                                           <  2:639964/139‑1 (MQ=255)                                                             aGCGCCAATGGTTATATTCGGGCTTCTGGCTGGAAAAAATCGCCGGATATTGAGGATTATCAGGAACGCCCGGCTAATGGCTGGGATATTCGTGCTGAg                                                                                                               >  2:285726/1‑99 (MQ=255)                                                             aGCGCCAATGGTTATATTCGGGCTTCTGGCTGGAAAAAATCGCCGGATATTGAGGATTATCAGGAACGCCCGGCTAATGGCTGGGATATTCGTGCTGAg                                                                                                               <  1:285726/99‑1 (MQ=255)                                                             aGCGCCAATGGTTATATTCGGGCTTCTGGCTGGAAAAAATCGCCGGATATTGAGGATTATCAGGAACGCCCGGCTAATGGCTGGGATATTCGTGCTGAGGGCTATTTACCCGCCTGGCCGCAGCTTGGCGCAAGCCTGa                                                                       >  1:987136/1‑139 (MQ=255)                                                             aGCGCCAATGGTTATATTCGGGCTTCTGGCTGGAAAAAATCGCCGGATATTGAGGATTATCAGGAACGCCCGGCTAATGGCTGGGATATTCGTGCTGAGGGCTATTTACCCGCCTGGCCGCAGCTTGGCGCAAGCCTGa                                                                       >  2:86256/1‑139 (MQ=255)                                                             aGCGCCAATGGTTATATTCGGGCTTCTGGCTGGAAAAAATCGCCGGATATTGAGGATTATCAGGAACGCCCGGCTAATGGCTGGGATATTCGTGCTGAGGGCTATTTACCCGCCTGGCCGCAGCTTGGCGCAAGCCTGa                                                                       >  1:725949/1‑139 (MQ=255)                                                             aGCGCCAATGGTTATATTCGGGCTTCTGGCTGGAAAAAATCGCCGGATATTGAGGATTATCAGGAACGCCCGGCTAATGGCTGGGATATTCGTGCTGAGGGCTATTTACCCGCCTGGCCGCAGCTTGGCGCAAGCCTGa                                                                       >  1:493165/1‑139 (MQ=255)                                                              gcgcCAATGGTTATATTCGGGCTTCTGGCTGGAAAAAATCGCCGGATATTGAGGATTATCAGGAACGCCCGGCTAATGGCTGGGATATTCGTGCTGAGGGCTATTTACCCGCCTGGCCGCAGCTTGGCGCAAGCCTGAt                                                                      <  1:611503/139‑1 (MQ=255)                                                                gcCAATGGTTATATTCGGGCTTCTGGCTGGAAAAAATCGCCGGATATTGAGGATTATCAGGAACGCCCGGCTAATGGCTGGGATATTCGTGCTGAGGGCTATTTACCCGCCTGGCCGCAGCTTGGCGCAAGCCTGATGt                                                                    <  1:776283/139‑1 (MQ=255)                                                                gcCAATGGTTATATTCGGGCTTCTGGCTGGAAAAAATCGCCGGATATTGAGGATTATCAGGAACGCCCGGCTAATGGCTGGGATATTCGTGCTGAGGGCTATTTACCCGCCTGGCCGCAGCTTGGCGCAAGCCTGATGt                                                                    <  2:556647/139‑1 (MQ=255)                                                                 ccAATGGTTATATTCGGGCTTCTGGCTGGAAAAAATCGCCGGATATTGAGGATTATCAGGAACGCCCGGCTAATGGCTGGGATATTCGTGCTGAGGGCTATTTACCCGCCTGGCCGCAGCTTGGCGCAAGCCTGATGTa                                                                   <  2:493165/139‑1 (MQ=255)                                                                      ggTTATATTCGGGCTTCTGGCTGGAAAAAATCGCCGGATATTGAGGATTATCAGGAACGCCCGGCTAATGGCTGGGATATTCGTGCTGAGGGCTATTTACCCGCCTGGCCGCAGCTTGGCGCAAGCCTGATGTATGAAc                                                              <  2:603789/139‑1 (MQ=255)                                                                            atTCGGGCTTCTGGCTGGAAAAAATCGCCGGATATTGAGGATTATCAGGAACGCCCGGCTAATGACTGGGa                                                                                                                            <  2:354906/71‑1 (MQ=255)                                                                            atTCGGGCTTCTGGCTGGAAAAAATCGCCGGATATTGAGGATTATCAGGAACGCCCGGCTAATGACTGGGa                                                                                                                            >  1:354906/1‑71 (MQ=255)                                                                                 ggCTTCTGGCTGGAAAAAATCGCCGGATATTGAGGATTATCAGGAACGCCCGGCTAATGGCTGGGATATTCGTGCTGAGGGCTATTTACCCGCCTGGCCGCAGCTTGGCGCAAGCCTGATGTATGAACAGTATTATGGc                                                   >  2:255064/1‑139 (MQ=255)                                                                                       tggctggAAAAAATCGCCGGATATTGAGGATTATCAGGAACGCCCGGCTAATGGCTGGGATATTCGTGCTGAGGGCTATTTACCCGCCTGGCCGCAGCTTGGCGCAAGCCTGATGt                                                                    <  1:316989/116‑1 (MQ=255)                                                                                       tggctggAAAAAATCGCCGGATATTGAGGATTATCAGGAACGCCCGGCTAATGGCTGGGATATTCGTGCTGAGGGCTATTTACCCGCCTGGCCGCAGCTTGGCGCAAGCCTGATGt                                                                    >  2:316989/1‑116 (MQ=255)                                                                                       tggctggAAAAAATCGCCGGATATTGAGGATTATCAGGAACGCCCGGCTAATGGCTGGGATATTCGTGCTGAGGGCTATTTACCCGCCTGGCCGCAGCTTGGCGCAAGCCTGATGTATGAACAGTATTATGGCGATGaa                                             <  2:670074/139‑1 (MQ=255)                                                                                             gAAAAAATCGCCGGATATTGAGGATTATCAGGAACGCCCGGCTAATGGCTGGGATATTCGTGCTGAGGGCTATTTACCCGCCTGGCCGCAGCTTGGCGCAAGCCTGATGTATGAACAGTATTATGGCGATGAAGTCggg                                       <  2:284798/139‑1 (MQ=255)                                                                                             gAAAAAATCGCCGGATATTGAGGATTATCAGGAACGCCCGGCTAATGGCTGGGATATTCGTGCTGAGGGCTATTTACCCGCCTGGCCGCAGCTTGGCGCAAGCCTGATGTATGAACAGTATTATGGCGATGAAGTCggg                                       <  2:671250/139‑1 (MQ=255)                                                                                             gAAAAAATCGCCGGATATTGAGGATTATCAGGAACGCCCGGCTAATGGCTGGGATATTCGTGCTGAGGGCTATTTACCCGCCTGGCCGCAGCTTGGCGCAAGCCTGATGTATGAACAGTATTATGGCGATGAAGTCggg                                       <  1:86256/139‑1 (MQ=255)                                                                                                                     ttATCAGGAACGCCCGGCTAATGGCTGGGATATTCGTGCTGAGGGCTATTTACCCGCCTGGCCGCAGCTTGGCGCAAGCCTGATGTATGAACAGTATTATGGCGATGAAGTCGGGCTGTTTGGTAAAGATAAGCGCCAg               >  2:364540/1‑139 (MQ=255)                                                                                                                      tATCAGGAACGCCCGGCTAATGGCTGGGATATTCGTGCTGAGGGCTATTTACCCGCCTGGCCGCAGCTTGGCGCAAGCCt                                                                         <  2:279018/80‑1 (MQ=255)                                                                                                                      tATCAGGAACGCCCGGCTAATGGCTGGGATATTCGTGCTGAGGGCTATTTACCCGCCTGGCCGCAGCTTGGCGCAAGCCt                                                                         >  1:279018/1‑80 (MQ=255)                                                                                                                      tATCAGGAACGCCCGGCTAATGGCTGGGATATTCGTGCTGAGGGCTATTTACCCGCCTGGCCGCAGCTTGGCGCAAGCCTGATGTATGAACAGTATTATGGCGATGAAGTCGGGCTg                                    >  2:720284/1‑117 (MQ=255)                                                                                                                      tATCAGGAACGCCCGGCTAATGGCTGGGATATTCGTGCTGAGGGCTATTTACCCGCCTGGCCGCAGCTTGGCGCAAGCCTGATGTATGAACAGTATTATGGCGATGAAGTCGGGCTg                                    <  1:720284/117‑1 (MQ=255)                                                                                                                      tATCAGGAACGCCCGGCTAATGGCTGGGATATTCGTGCTGAGGGCTATTTACCCGCCTGGCCGCAGCTTGGCGCAAGCCTGATGTATGAACAGTATTATGGCGATGAAGTCGGGCTGtt                                  >  1:222023/1‑119 (MQ=255)                                                                                                                      tATCAGGAACGCCCGGCTAATGGCTGGGATATTCGTGCTGAGGGCTATTTACCCGCCTGGCCGCAGCTTGGCGCAAGCCTGATGTATGAACAGTATTATGGCGATGAAGTCGGGCTGtt                                  <  2:222023/119‑1 (MQ=255)                                                                                                                      tATCAGGAACGCCCGGCTAATGGCTGGGATATTCGTGCTGAGGGCTATTTACCCGCCTGGCCGCAGCTTGGCGCAAGCCTGATGTATGAACAGTATTATGGCGATGAAGTCGGGCTGTTTGGTAAAGATAAGCGCCAGa              >  1:1111182/1‑139 (MQ=255)                                                                                                                       aTCAGGAACGCCCGGCTAATGGCTGGGATATTCGTGCTGAGGGCTATTTACCCGCCTGGCCGCAGCTTGGCGCAAGCCTGATGTATGAACAGTATTATGGCGATGAAGTCGGGCTGTTTGGTAAAGATAAGCGCCAGaa             >  2:375494/1‑139 (MQ=255)                                                                                                                       aTCAGGAACGCCCGGCTAATGGCTGGGATATTCGTGCTGAGGGCTATTTACCCGCCTGGCCGCAGCTTGGCGCAAGCCTGATGTATGAACAGTATTATGGCGATGAAGTCGGGCTGTTTGGTAAAGATAAGCGCCAGaa             >  1:334581/1‑139 (MQ=255)                                                                                                                        tCAGGAACGCCCGGCTAATGGCTGGGATATTCGTGCTGAGGGCTATTTACCCGCCTGGCCGCAGCTTGGCGCa                                                                              <  1:782579/73‑1 (MQ=255)                                                                                                                        tCAGGAACGCCCGGCTAATGGCTGGGATATTCGTGCTGAGGGCTATTTACCCGCCTGGCCGCAGCTTGGCGCa                                                                              >  2:782579/1‑73 (MQ=255)                                                                                                                              aCGCCCGGCTAATGGCTGGGATATTCGTGCTGAGGGCTATTTACCCGCCTGGCCGCAGCTTGGCGCAAGCCTGATGTATGAACAGTATTATGGCGATGAAGTCGGGCTGTTTGGTAAAGATAAGCGCCAGAAAGACCCg      <  2:1081357/139‑1 (MQ=255)                                                                                                                              aCGCCCGGCTAATGGCTGGGATATTCGTGCTGAGGGCTATTTACCCGCCTGGCCGCAGCTTGGCGCAAGCCTGATGTATGAACAGTATTATGGCGATGAAGTCGGGCTGTTTGGTAAAGATAAGCGCCAGAAAGACCCg      <  1:385893/139‑1 (MQ=255)                                                                                                                              aCGCCCGGCTAATGGCTGGGATATTCGTGCTGAGGGCTATTTACCCGCCTGGCCGCAGCTTGGCGCAAGCCTGATGTATGAACAGTATTATGGCGATGAAGTCGGGCTGTTTGGTAAAGATAAGCGCCAGAAAGACCCg      >  2:906811/1‑139 (MQ=255)                                                                                                                                gCCCGGCTAATGGCTGGGATATTCGTGCTGAGGGCTATTTACCCGCCTGGCCGCAGCTTGGCGCAAGCCTGATGTATg                                                                 >  2:274372/1‑78 (MQ=255)                                                                                                                                gCCCGGCTAATGGCTGGGATATTCGTGCTGAGGGCTATTTACCCGCCTGGCCGCAGCTTGGCGCAAGCCTGATGTATg                                                                 <  1:274372/78‑1 (MQ=255)                                                                                                                                  ccGGCTAATGGCTGGGATATTCGTGCTGAGGGCTATTTACCCGCCTGGCCGCAGCTTGGCGCAAGCCTGATGTATGAACAGTATTATGGCGATGAAGTCGGGCTGTTTGGTAAAGATAAGCGCCAGAAAGACCCGCATg  >  1:1021057/1‑139 (MQ=255)                                                                                                                                       |                                                                                                                                      TCCCGTAGTCATACCCGCATTGGTGTTGGTGCGGAATACTGGCGCGATTATCTGAAACTGAGCGCCAATGGTTATATTCGGGCTTCTGGCTGGAAAAAATCGCCGGATATTGAGGATTATCAGGAACGCCCGGCGAATGGCTGGGATATTCGTGCTGAGGGCTATTTACCCGCCTGGCCGCAGCTTGGCGCAAGCCTGATGTATGAACAGTATTATGGCGATGAAGTCGGGCTGTTTGGTAAAGATAAGCGCCAGAAAGACCCGCATG  >  NC_000913/314927‑315194

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