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%	101.0 / NA	34	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 (19/15);  total (19/15)

AAGCCTGATGTATGAACAGTATTATGGCGATGAAGTCGGGCTGTTTGGTAAAGATAAGCGCCAGAAAGACCCGCATGCTATTTCTGCCGAGGTGACCTATACGCCAGTGCCTCTTACCCAGCAATAGTGGACACGCGGCTAAGTGAGTAAACTCTCAGTCAGAGGTGACTCACATGACAAAAACAGTATCAACCAGTAAAAAACCCCGTAAACAGCATTCGCCTGAATTTCGCAGTGAAGCCCTGAAGC  >  NC_000913/315118‑315366                                                                                                                               |                                                                                                                           aaGCCTGATGTATGAACAGTATTATGGCGATGAAGTCGGGCTGTTTGGTAAAGATAAGCGCCAGAAAGACCCGCATGCTATTTCTGCCGAGGTGACCTATACGCCAGTGCCTCTTACCCAGCAATAATGGACACGCGGc                                                                                                                >  2:243086/1‑139 (MQ=255)   gCCTGATGTATGAACAGTATTATGGCGATGAAGTCGGGCTGTTTGGTAAAGATAAGCGCCAGAAAGACCCGCATGCTATTTCTGCCGAGGTGACCTATACGCCAGTGCCTCTTACCCAGCAATAATGGACACGCGGCTa                                                                                                              <  2:1228809/139‑1 (MQ=255)     cTGATGTATGAACAGTATTATGGCGATGAAGTCGGGCTGTTTGGTAAAGATAAGCGCCAGAAAGACCCGCATGCTATTTCTGCCGAGGTGACCTATACGCCAGTGCCTCTTACCCAGCAATAATGGACACGCGGCTAAg                                                                                                            <  2:704594/139‑1 (MQ=255)           tATGAACAGTATTATGGCGATGAAGTCGGGCTGTTTGGTAAAGATAAGCGCCAGAAAGACCCGCATGCTATTTCTGCCGAGGTGACCTATACGCCAGTGCCTCTTACCCAGCAATAATGGACACGCGGCTAAGTGAGTa                                                                                                      <  1:537001/139‑1 (MQ=255)             tGAACAGTATTATGGCGATGAAGTCGGGCTGTTTGGTAAAGATAAGCGCCAGAAAGACCCGCATGCTATTTCTGCCGAGGTGACCTATACGCCAGTGCCTCTTACCCAGCAATAATGGACACGCGGCTAAGTg                                                                                                          >  1:107803/1‑133 (MQ=255)             tGAACAGTATTATGGCGATGAAGTCGGGCTGTTTGGTAAAGATAAGCGCCAGAAAGACCCGCATGCTATTTCTGCCGAGGTGACCTATACGCCAGTGCCTCTTACCCAGCAATAATGGACACGCGGCTAAGTg                                                                                                          <  2:107803/133‑1 (MQ=255)                  aGTATTATGGCGATGAAGTCGGGCTGTTTGGTAAAGATAAGCGCCAGAAAGACCCGCATGCTATTTCTGCCGAGGTGACCTATACGCCAGTGCCTCTTACCCAGCAATAATGGACACGCGGCTAAGTGAGTAAACtctc                                                                                               >  1:694287/1‑139 (MQ=255)                   gTATTATGGCGATGAAGTCGGGCTGTTTGGTAAAGATAAGCGCCAGAAAGACCCGCATGCTATTTCTGCCGAGGTGACCTATACGCCAGTGCCTCTTACCCAGCAATAATGGACACGCGGCTAAGTGAGTAAACTCtca                                                                                              >  2:181769/1‑139 (MQ=255)                    tattatGGCGATGAAGTCGGGCTGTTTGGTAAAGATAAGCGCCAGAAAGACCCGCATGCTATTTCTGCCGAGGTGACCTATACGCCAGTGCCTCTTACCCAGCAATAATGGACACGCGGCTAAGTGAGTAAACTCtcag                                                                                             >  1:533523/1‑139 (MQ=255)                         tGGCGATGAAGTCGGGCTGTTTGGTAAAGATAAGCGCCAGAAAGACCCGCATGCTATTTCTGCCGAGGTGAACTATACGCCAGTGCCTCTTACCCAGCAATAATGGACACGCGGCTAAGTGAGTAAACTCTCAGTCaga                                                                                        >  2:1074094/1‑139 (MQ=255)                            cGATGAAGTCGGGCTGTTTGGTAAAGATAAGCGCCAGAAAGACCCGCATGCTATTTCTGCCGAGGTGACCTATACGCCAGTGCCTCTTACCCAGCAATAATGGACACGCGGCTAAGTGAGTAAACTCTCAGTCaga                                                                                        <  1:1296402/136‑1 (MQ=255)                            cGATGAAGTCGGGCTGTTTGGTAAAGATAAGCGCCAGAAAGACCCGCATGCTATTTCTGCCGAGGTGACCTATACGCCAGTGCCTCTTACCCAGCAATAATGGACACGCGGCTAAGTGAGTAAACTCTCAGTCaga                                                                                        >  2:1296402/1‑136 (MQ=255)                              aTGAAGTCGGGCTGTTTGGTAAAGATAAGCGCCAGAAAGACCCGCATGCTATTTCTGCCGAGGTGACCTATACGCCAGTGCCTCTTACCCAGCAATAATGGACACGCGGCTAAGTGAGTAAACTCTCAGTCAGAGGTGa                                                                                   >  1:882245/1‑139 (MQ=255)                                             ttGGTAAAGATAAGCGCCAGAAAGACCCGCATGCTATTTCTGCCGAGGTGACCTATACGCCAGTGCCTCTTACCCAGCAATAATGGACACGCGGCTAAGTGAGTAAACTCTCAGTCAGAGGTGACTCACATGACaaaaa                                                                    <  1:243086/139‑1 (MQ=255)                                             ttGGTAAAGATAAGCGCCAGAAAGACCCGCATGCTATTTCTGCCGAGGTGACCTATACGCCAGTGCCTCTTACCCAGCAATAATGGACACGCGGCTAAGTGAGTAAACTCTCAGTCAGAGGTGACTCACATGACaaaaa                                                                    <  1:498359/139‑1 (MQ=255)                                                      aTAAGCGCCAGAAAGACCCGCATGCTATTTCTGCCGAGGTGACCTATACGCCAGTGCCTCTTACCCAGCAATAATGGACACGCGGCTAAGTGAGTAAACTCTCAGTCAGAGGTGACTCACATGACAAAAACAGTATCaa                                                           >  2:989913/1‑139 (MQ=255)                                                           cgcCAGAAAGACCCGCATGCTATTTCTGCCGAGGTGACCTATACGCCAGTGCCTCTTACCCAGCAATAATGGACACGCGGCTAAGTGAGTAAACTCTCAGTCAGAGGTGACTCACATGACAAAAACAGTATCAACCAGt                                                      >  2:655552/1‑139 (MQ=255)                                                               aGAAAGACCCGCATGCTATTTCTGCCGAGGTGACCTATACGCCAGTGCCTCTTACCCAGCAATAATGGACACGCGGCTAAGTGAGTAAActc                                                                                                 <  1:1125138/92‑1 (MQ=255)                                                               aGAAAGACCCGCATGCTATTTCTGCCGAGGTGACCTATACGCCAGTGCCTCTTACCCAGCAATAATGGACACGCGGCTAAGTGAGTAAActc                                                                                                 >  2:1125138/1‑92 (MQ=255)                                                                          cATGCTATTTCTGCCGAGGTGACCTATACGCCAGTGCCTCTTACCCAGCAATAATGGACACGCGGCTAAGTGAGTAAACTCTCAGTCAGAGGTGACTCACATGACAAAAACAGTATCAACCAGTAAAAAACCCCGTaaa                                       <  1:1048447/139‑1 (MQ=255)                                                                                 tttCTGCCGAGGTGACCTATACGCCAGTGCCTCTTACCCAGCAATAATGGACACGCGGCTAAGTGAGTAAACTCtcagtca                                                                                          <  2:248897/81‑1 (MQ=255)                                                                                 tttCTGCCGAGGTGACCTATACGCCAGTGCCTCTTACCCAGCAATAATGGACACGCGGCTAAGTGAGTAAACTCtcagtca                                                                                          >  1:248897/1‑81 (MQ=255)                                                                                           ggTGACCTATACGCCAGTGCCTCTTACCCAGCAATAATGGACACGCGGCTAAGTGAGTAAACTCTCAGTCAGAGGTGACTCACATGACaaaaa                                                                    >  1:644476/1‑93 (MQ=255)                                                                                           ggTGACCTATACGCCAGTGCCTCTTACCCAGCAATAATGGACACGCGGCTAAGTGAGTAAACTCTCAGTCAGAGGTGACTCACATGACaaaaa                                                                    <  2:644476/93‑1 (MQ=255)                                                                                                   ataCGCCAGTGCCTCTTACCCAGCAATAATGGACACGCGGCTAAGTGAGTAAACTCTCAGTCAGAGGTGACTCACATGACAAAAACAGTATCAACCAGTAAAAAACCCCGTAAACAGCATTCGCCTGAATTTCGCAGTg              <  1:655552/139‑1 (MQ=255)                                                                                                    taCGCCAGTGCCTCTTACCCAGCAATAATGGACACGCGGCTAAGTGAGTAAACTCTCAGTCAGAGGTGACTCACATGACAAAAACAGTATCAACCAGTAAAAAACCCCGTAAACAGCATTCGCCTGAATTTCGCAGTGa             >  1:651252/1‑139 (MQ=255)                                                                                                         cAGTGCCTCTTACCCAGCAATAATGGACACGCGGCTAAGTGAGTAAACTCTCAGTCAGAGGTGACTCACATGACAAAAACAGTATCAAc                                                          <  2:951937/89‑1 (MQ=255)                                                                                                         cAGTGCCTCTTACCCAGCAATAATGGACACGCGGCTAAGTGAGTAAACTCTCAGTCAGAGGTGACTCACATGACAAAAACAGTATCAAc                                                          >  1:951937/1‑89 (MQ=255)                                                                                                         cAGTGCCTCTTACCCAGCAATAATGGACACGCGGCTAAGTGAGTAAACTCTCAGTCAGAGGTGACTCACATGACAAAAACAGTATCAACCAGTAAAAAACCCCGTAAACAGCATTCGCCTGAATTTCGCAGTGAAGccc        >  1:787536/1‑139 (MQ=25)                                                                                                           gTGCCTCTTACCCAGCAATAATGGACACGCGGCTAAGTGAGTAAACTCTCAGTCAGAGGTGACTCACATGCCAAAAACAGTATCAACCAGTAAAAAACCCCGTAAACAGCATTCGCCTGAATTTCGCAGTg              <  2:449590/131‑1 (MQ=16)                                                                                                           gTGCCTCTTACCCAGCAATAATGGACACGCGGCTAAGTGAGTAAACTCTCAGTCAGAGGTGACTCACATGCCAAAAACAGTATCAACCAGTAAAAAACCCCGTAAACAGCATTCGCCTGAATTTCGCAGTg              >  1:449590/1‑131 (MQ=16)                                                                                                           gTGCCTCTTACCCAGCAATAATGGACACGCGGCTAAGTGAGTAAACTCTCAGTCAGAGGTGACTCACATGACAAAAACAGTATCAACCAGTAAAAAACCCCGTAAACAGCATTCGCCTGAATTTCGCAGTGAAGCCCTg      >  2:1042899/1‑139 (MQ=18)                                                                                                              ccTCTTACCCAGCAATAATGGACACGCGGCTAAGTGAGTAAACTCTCAGTCAGAGGTGACTCACATGACAAAAACAGTATCAACCAGTAAAAAACCCCGTAAACAGCATTCGCCTGAATTTCGCAGTGAAGCCCTGAAg   >  2:235747/1‑139 (MQ=14)                                                                                                               ctctTACCCAGCAATAATGGACACGCGGCTAAGTGAGTAAACTCTCAGTCAGAGGTGACTCACATGACAAAAACAGTATCAACCAGTAAAAAACCCCGTAAACAGCATTCGCCTGAATTTCGCAGTGAAGCCCTGAAGc  <  1:989913/139‑1 (MQ=11)                                                                                                                               |                                                                                                                           AAGCCTGATGTATGAACAGTATTATGGCGATGAAGTCGGGCTGTTTGGTAAAGATAAGCGCCAGAAAGACCCGCATGCTATTTCTGCCGAGGTGACCTATACGCCAGTGCCTCTTACCCAGCAATAGTGGACACGCGGCTAAGTGAGTAAACTCTCAGTCAGAGGTGACTCACATGACAAAAACAGTATCAACCAGTAAAAAACCCCGTAAACAGCATTCGCCTGAATTTCGCAGTGAAGCCCTGAAGC  >  NC_000913/315118‑315366

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: ‑