Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	NC_000913	315,244	G→A	intergenic (+16/‑47)	eaeH → / → insE1	pseudogene, attaching and effacing protein homology,factor, Not classified/IS3 transposase A

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	315,244	0	G	A	100.0%	151.8 / NA	48	intergenic (+16/‑47)	eaeH/insE1	pseudogene, attaching and effacing protein homology,factor, Not classified/IS3 transposase A
Reads supporting (aligned to +/- strand):  ref base G (0/0);  new base A (26/22);  total (26/22)

GCTTGGCGCAAGCCTGATGTATGAACAGTATTATGGCGATGAAGTCGGGCTGTTTGGTAAAGATAAGCGCCAGAAAGACCCGCATGCTATTTCTGCCGAGGTGACCTATACGCCAGTGCCTCTTACCCAGCAATAGTGGACACGCGGCTAAGTGAGTAAACTCTCAGTCAGAGGTGACTCACATGACAAAAACAGTATCAACCAGTAAAAAACCCCGTAAACAGCATTCGCCTGAATTTCGCAGTGAAGCCCTGA  >  NC_000913/315109‑315363                                                                                                                                        |                                                                                                                        gCTTGGCGCAAGCCTGATGTATGAACAGTATTATGGCGATGAAGTCGGGCTGTTTGGTAAAGATAAGCGCCAGAAAGACCCGCATGCTATTTCTGCCGAGGTGACCTATACGCCAGTGCCTCTTACCCAGCAATAATgg                                                                                                                      >  2:1262703/1‑139 (MQ=255)  cTTGGCGCAAGCCTGATGTATGAACAGTATTATGGCGATGAAGTCGGGCTGTTTGGTAAAGATAAGCGCCAGAAAGACCCGCATGCTATTTCTGCCGAGGTGACCTATACGCCAGTGCCTCTTACCCAGCAATAATGGa                                                                                                                     >  2:299266/1‑139 (MQ=255)     ggCGCAAGCCTGATGTATGAACAGTATTATGGCGATGAAGTCGGGCTGTTTGGTAAAGATAAGCGCCAGAAAGACCCGCATGCTATTTCTGCCGAGGTGACCTATACGCCAGTGCCTCTTACCCAGCAATAATGGacac                                                                                                                  <  2:1205939/139‑1 (MQ=255)              cTGATGTATGAACAGTATTATGGCGATGAAGTCGGGCTGTTTGGTAAAGATAAGCGCCAGAAAGACCCGCATGCTATTTCTGCCGAGGTGACCTATACGCCAGTGCCTCTTACCCAGCAATAATGGACACGCGGCTAAg                                                                                                         <  2:1039567/139‑1 (MQ=255)                 aTGTATGAACAGTATTATGGCGATGAAGTCGGGCTGTTTGGTAAAGATAAGCGCCAGAAAGACCCGCATGCTATTTCTGCCGAGGTGACCTATACGCCAGTGCCTCTTACCCAGCAATAATGGACACGCGGCTAAGTGa                                                                                                      >  1:1276577/1‑139 (MQ=255)                        aaCAGTATTATGGCGATGAAGTCGGGCTGTTTGGTAAAGATAAGCGCCAGAAAGACCCGCATGCTATTTCTGCCGAGGTGACCTATACGCCAGTGCCTCTTACCCAGCAATAATGGACACGCGGCTAAGTGAGTAAAct                                                                                               <  2:1276577/139‑1 (MQ=255)                             tattatGGCGATGAAGTCGGGCTGTTTGGTAAAGATAAGCGCCAGAAAGACCCGCATGCTATTTCTGCCGAGGTGACCTATACGCCAGTGCCTCTTACCCAGCAATAATGGACACGCGGCTAAGTGAGTa                                                                                                   >  1:1051957/1‑130 (MQ=255)                             tattatGGCGATGAAGTCGGGCTGTTTGGTAAAGATAAGCGCCAGAAAGACCCGCATGCTATTTCTGCCGAGGTGACCTATACGCCAGTGCCTCTTACCCAGCAATAATGGACACGCGGCTAAGTGAGTa                                                                                                   <  2:1051957/130‑1 (MQ=255)                             tattatGGCGATGAAGTCGGGCTGTTTGGTAAAGATAAGCGCCAGAAAGACCCGCATGCTATTTCTGCCGAGGTGACCTATACGCCAGTGCCTCTTACCCAGCAATAATGGACACGCGGCTAAGTGAGTAAACTCtcag                                                                                          >  1:292501/1‑139 (MQ=255)                                     cGATGAAGTCGGGCTGTTTGGTAAAGATAAGCGCCAGAAAGACCCGCATGCTATTTCTGCCGAGGTGACCTATACGCCAGTGCCTCTTACCCAGCAATAATGGACACGCGGCTAAGTGAGTAAACTCTCAGTCAGAGGt                                                                                  >  2:1375741/1‑139 (MQ=255)                                              cGGGCTGTTTGGTAAAGATAAGCGCCAGAAAGACCCGCATGCTATTTCTGCCGAGGTGACCTATACGCCAGTGCCTCTTACCCAGCAATAATGGACACGCGGCTAAGTGAGTAAACTCTCAGTCAGAGGTGACTCACAt                                                                         <  1:299266/139‑1 (MQ=255)                                               gggCTGTTTGGTAAAGATAAGCGCCAGAAAGACCCGCATGCTATTTCTGCCGAGGTGACCTATACGCCAGTGCCTCTTACCCAGCAATAATGGACACGCGGCTAAGTGAGTAAACTCTCAGTCAGAGGTGACTCACATg                                                                        >  1:251450/1‑139 (MQ=255)                                                             aGATAAGCGCCAGAAAGACCCGCATGCTATTTCTGCCGAGGTGACCTATACGCCAGTGCCTCTTACCCAGCAATAATGGACACGCGGCTAAGTGAGTAAACTCTCAGTCAGAGGTGACTCACATGACAAAAACAGTATc                                                          <  1:1262703/139‑1 (MQ=255)                                                               aTAAGCGCCAGAAAGACCCGCATGCTATTTCTGCCGAGGTGACCTATACGCCAGTGCCTCTTACCCAGCAATAATGGACACGCGGCTAAGTGAGTAAACTCTCAGTCAGAGGTGACTCACAt                                                                         >  2:788335/1‑122 (MQ=255)                                                               aTAAGCGCCAGAAAGACCCGCATGCTATTTCTGCCGAGGTGACCTATACGCCAGTGCCTCTTACCCAGCAATAATGGACACGCGGCTAAGTGAGTAAACTCTCAGTCAGAGGTGACTCACAt                                                                         <  1:788335/122‑1 (MQ=255)                                                               aTAAGCGCCAGAAAGACCCGCATGCTATTTCTGCCGAGGTGACCTATACGCCAGTGCCTCTTA‑CCAGCAATAATGGACACGCGGCTAAGTGAGTaaa                                                                                                 >  1:421379/1‑97 (MQ=255)                                                               aTAAGCGCCAGAAAGACCCGCATGCTATTTCTGCCGAGGTGACCTATACGCCAGTGCCTCTTA‑CCAGCAATAATGGACACGCGGCTAAGTGAGTaaa                                                                                                 <  2:421379/97‑1 (MQ=255)                                                               aTAAGCGCCAGAAAGACCCGCATGCTATTTCTGCCGAGGTGACCTATACACCAGTGCCTCTTACCCAGCAATAATGGACACGCGGCTAAGTGAGTAAACTCTCAGTCAGAGGTGACTCACATGACAAAAACAGTATCaa                                                        >  1:617472/1‑139 (MQ=255)                                                                    cgcCAGAAAGACCCGCATGCTATTTCTGCCGAGGTGACCTATACGCCAGTGCCTCTTACCCAGCAATAATGGACACGCGGCTAAGTGAGTAAActc                                                                                              <  2:1275021/96‑1 (MQ=255)                                                                    cgcCAGAAAGACCCGCATGCTATTTCTGCCGAGGTGACCTATACGCCAGTGCCTCTTACCCAGCAATAATGGACACGCGGCTAAGTGAGTAAActc                                                                                              >  1:1275021/1‑96 (MQ=255)                                                                      ccAGAAAGACCCGCATGCTATTTCTGCCGAGGTGACCTATACGCCAGTGCCTCTTACCCAGCAATAATGGACACGCGGCTAAGTGAGTAAACTCTCAGTCAGAGGTGACTCACATGACAAAAACAGTATCAACCAGTaa                                                 >  2:421156/1‑139 (MQ=255)                                                                      ccAGAAAGACCCGCATGCTATTTCTGCCGAGGTGACCTATACGCCAGTGCCTCTTACCCAGCAATAATGGACACGCGGCTAAGTGAGTAAACTCTCAGTCAGAGGTGACTCACATGACAAAAACAGTATCAACCAGTaa                                                 >  2:1456374/1‑139 (MQ=255)                                                                           aaGACCCGCATGCTATTTCTGCCGAGGTGACCTATACGCCAGTGCCTCTTACCCAGCAATAATGGACACGCGGCTAAGTGAGTAAACTCTCAGTCAGAGGTGACTCACATGACAAAAACAGTATCAACCAGTAAAAAAc                                            <  1:421156/139‑1 (MQ=255)                                                                           aaGACCCGCATGCTATTTCTGCCGAGGTGACCTATACGCCAGTGCCTCTTACCCAGCAATAATGGACACGCGGCTAAGTGAGTAAACTCTCAGTCAGAGGTGACTCACATGACAAAAACAGTATCAACCAGTAAAAAAc                                            <  1:1456374/139‑1 (MQ=255)                                                                           aaGACCCGCATGCTATTTCTGCCGAGGTGACCTATACGCCAGTGCCTCTTACCCAGCAATAATGGACACGCGGCTAAGTGAGTAAACTCTCAGTCAGAGGTGACTCACATGACAAAAACAGTATCAACCAGTAAAAAAc                                            <  2:292501/139‑1 (MQ=255)                                                                             gACCCGCATGCTATTTCTGCCGAGGTGACCTATACGCCAGTGCCTCTTACCCAGCAATAATGGACACGCGGCTAAGTGAGTAAACTCTCAGTCAGAGGTGACTCACATGAc                                                                      <  1:1639724/111‑1 (MQ=255)                                                                             gACCCGCATGCTATTTCTGCCGAGGTGACCTATACGCCAGTGCCTCTTACCCAGCAATAATGGACACGCGGCTAAGTGAGTAAACTCTCAGTCAGAGGTGACTCACATGAc                                                                      >  2:1639724/1‑111 (MQ=255)                                                                               cccGCATGCTATTTCTGCCGAGGTGACCTATACGCCAGTGCCTCTTACCCAGCAATAATGGACACGCGGCTAAGTGAGTAAActct                                                                                             <  1:519265/86‑1 (MQ=255)                                                                               cccGCATGCTATTTCTGCCGAGGTGACCTATACGCCAGTGCCTCTTACCCAGCAATAATGGACACGCGGCTAAGTGAGTAAActct                                                                                             >  2:519265/1‑86 (MQ=255)                                                                                   cATGCTATTTCTGCCGAGGTGACCTATACGCCAGTGCCTCTTACCCAGCAATAATGGACACGCGGCTAAGTGAGTAAACTCTCAGTCAGAGGTGACTCACATGACAAAAACAGTATCAACCAGTAAAAAACCCCGTaaa                                    >  2:816652/1‑139 (MQ=255)                                                                                    aTGCTATTTCTGCCGAGGTGACCTATACGCCAGTGCCTCTTACCCAGCAATAATGGACACGCGGCTAAGTGAGTAAACTCtcagt                                                                                         >  1:1361646/1‑85 (MQ=255)                                                                                    aTGCTATTTCTGCCGAGGTGACCTATACGCCAGTGCCTCTTACCCAGCAATAATGGACACGCGGCTAAGTGAGTAAACTCtcagt                                                                                         <  2:1361646/85‑1 (MQ=255)                                                                                    aTGCTATTTCTGCCGAGGTGACCTATACGCCAGTGCCTCTTACCCAGCAATAATGGACACGCGGCTAAGTGAGTAAACTCTCAGTCAGAGGTGACTCACATGACAAAAACAGTATCAACCAGTAAAAAACCCCGTAAAc                                   >  2:517679/1‑139 (MQ=255)                                                                                    aTGCTATTTCTGCCGAGGTGACCTATACGCCAGTGCCTCTTACCCAGCAATAATGGACACGCGGCTAAGTGAGTAAACTCTCAGTCAGAGGTGACACACATGACAAAAACAGTATCAACCAGTAAAAAACCCCGTAAAc                                   >  2:1126316/1‑139 (MQ=255)                                                                                                ccGAGGTGACCTATACGCCAGTGCCTCTTACCCAGCAATAATGGACAcgcg                                                                                                               >  2:212247/1‑51 (MQ=255)                                                                                                ccGAGGTGACCTATACGCCAGTGCCTCTTACCCAGCAATAATGGACAcgcg                                                                                                               <  1:212247/51‑1 (MQ=255)                                                                                                  gAGGTGACCTATACACCAGTGCCTCTTACCCAGCAATAATGGACACGCGGCTAAGTGAGTAAACTCTCAGTCAGAGGTGACTCACATGACAAAAACAGTATCAACCAGTAAAAAACCCCGTAAACAGCATTCGCCTGaa                     <  2:617472/139‑1 (MQ=255)                                                                                                                  cAGTGCCTCTTACCCAGCAATAATGGACACGCGGCTAAGTGAGTAAACTCtcagt                                                                                         >  2:208646/1‑55 (MQ=255)                                                                                                                  cAGTGCCTCTTACCCAGCAATAATGGACACGCGGCTAAGTGAGTAAACTCtcagt                                                                                         <  1:208646/55‑1 (MQ=255)                                                                                                                    gTGCCTCTTACCCAGCAATAATGGACACGCGGCTAAGTGAGTAAACTCTCAGTCAGAGGTGACTCACATGACaaaaa                                                                 >  2:1670557/1‑77 (MQ=255)                                                                                                                    gTGCCTCTTACCCAGCAATAATGGACACGCGGCTAAGTGAGTAAACTCTCAGTCAGAGGTGACTCACATGACaaaaa                                                                 <  1:1670557/77‑1 (MQ=255)                                                                                                                    gTGCCTCTTACCCAGCAATAATGGACACGCGGCTAAGTGAGTAAACTCTCAGTCAGAGGTGACTCACATGACaaaa                                                                  >  1:812032/1‑76 (MQ=255)                                                                                                                    gTGCCTCTTACCCAGCAATAATGGACACGCGGCTAAGTGAGTAAACTCTCAGTCAGAGGTGACTCACATGACaaaa                                                                  <  2:812032/76‑1 (MQ=255)                                                                                                                     tGCCTCTTACCCAGCAATAATGGACACGCGGCTAAGTGAGTAAACTCTCAGTCAGAGGTGACTCACATGACaaa                                                                   <  1:882071/74‑1 (MQ=25)                                                                                                                     tGCCTCTTACCCAGCAATAATGGACACGCGGCTAAGTGAGTAAACTCTCAGTCAGAGGTGACTCACATGACaaa                                                                   >  2:882071/1‑74 (MQ=25)                                                                                                                     tGCCTCTTACCCAGCAATAATGGACACGCGGCTAAGTGAGTAAACTCTCAGTCAGAGGTGACTCACATGACAAAAACAGTATCAACCAGTAAAAAACCCCGTAAACAGCATTCGCCTGAATTTCGCAGTGAAGCCCTGa  >  1:905102/1‑139 (MQ=18)                                                                                                                       ccTCTTACCCAGCAATAATGGACACGCGGCTAAGTGAGTAAACTCTCAGTCAGAGGTGACTCACATGACaaaaa                                                                 >  2:1727436/1‑74 (MQ=17)                                                                                                                       ccTCTTACCCAGCAATAATGGACACGCGGCTAAGTGAGTAAACTCTCAGTCAGAGGTGACTCACATGACaaaaa                                                                 <  1:1727436/74‑1 (MQ=17)                                                                                                                                        |                                                                                                                        GCTTGGCGCAAGCCTGATGTATGAACAGTATTATGGCGATGAAGTCGGGCTGTTTGGTAAAGATAAGCGCCAGAAAGACCCGCATGCTATTTCTGCCGAGGTGACCTATACGCCAGTGCCTCTTACCCAGCAATAGTGGACACGCGGCTAAGTGAGTAAACTCTCAGTCAGAGGTGACTCACATGACAAAAACAGTATCAACCAGTAAAAAACCCCGTAAACAGCATTCGCCTGAATTTCGCAGTGAAGCCCTGA  >  NC_000913/315109‑315363

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