Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	NC_000913	4,585,372	G→A	Q464* (CAA→TAA)	hsdR ←	type I restriction enzyme EcoKI endonuclease component

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	4,585,372	0	G	A	100.0%	121.9 / NA	37	Q464* (CAA→TAA)	hsdR	type I restriction enzyme EcoKI endonuclease component
Reads supporting (aligned to +/- strand):  ref base G (0/0);  new base A (20/17);  total (20/17)

TTCGCCAAGAGAACGGCGGTCGACAAGGAAGAGAATGCGTTTAAAACGCTGGGACTGGATCAGGCGGAACATCATGGCGATTGCCGTACGGGTTTTACCGGTACCGGTCGCCATCGCCAGCAGGATCTCTTGTTGCCCCTTGACGATTGCCTTTTCAACCGCGCGGACGGCATCTTCCTGATAATAACGCAGGCCCAGCTCGCTCATGCCAGGGTTATCGGCAAACCACTGATTCTGTTTTTGCGGTTCGCTGCCGAGCATTTCCAGCAG  >  NC_000913/4585238‑4585507                                                                                                                                       |                                                                                                                                        ttCGCCAAGAGAACGGCGGTCGACAAGGAAGAGAATGCGTTTAAAACGCTGGGACTGGATCAGGCGGAACATCATGGCGATTGCCGTACGGGTTTTACCGGTACCGGTCGCCATCGCCAGCAGGATCTCTTGTTACCCCt                                                                                                                                    <  2:658349/140‑1 (MQ=255)     ccAAGAGAACGGCGGTCGACAAGGAAGAGAATGCGTTTTAAACGCTGGGACTGGATCAGGCGGAACATCATGGCGATTGCCGTACGGGTTTTACCGGTACCGGTCGCCATCGCCAGCAGGATCTCTTGTTACCCCTTGAc                                                                                                                                >  2:706034/1‑140 (MQ=255)         gagaACGGCGGTCGACAAGGAAGAGAATGCGTTTAAAACGCTGGGACTGGATCAGGCGGAACATCATGGCGATTGCCGTACGGGTTTTACCGGTACCGGTCGCCATCGCCAGCAGGATCTCTTGTTACCCCTTGACGAtt                                                                                                                            <  1:801524/140‑1 (MQ=255)           gaACGGCGGTCGACAAGGAAGAGAATGCGTTTAAAACGCTGGGACTGGATCAGGCGGAACATCATGGCGATTGCCGTACGGGTTTTACCGGTACCGGTCGCCATCGCCAGCAGGATCTCTTGTTACCCCTTGACGATTGc                                                                                                                          <  2:576656/140‑1 (MQ=255)                   gTCGACAAGGAAGAGAATGCGTTTAAAACGCTGGGACTGGATCAGGCGGAACATCATGGCGATTGCCGTACGGGTTTTACCGGTACCGGTCGCCATCGCCAGCAGGATCTGTTGGTAGCCCTTGCCGATTGCCTTTTCaa                                                                                                                  >  1:294087/1‑140 (MQ=255)                    tCGACAAGGAAGAGAATGCGTTTAAAACGCTGGGACTGGATCAGGCGGAACATCATGGCGATTGCCGTACGGGTTTTACCGGTACCGGTCGCCATCGCCAGCAGGATCTCTTGTTACCCCTTGACGAtt                                                                                                                            >  2:710455/1‑129 (MQ=255)                    tCGACAAGGAAGAGAATGCGTTTAAAACGCTGGGACTGGATCAGGCGGAACATCATGGCGATTGCCGTACGGGTTTTACCGGTACCGGTCGCCATCGCCAGCAGGATCTCTTGTTACCCCTTGACGAtt                                                                                                                            <  1:710455/129‑1 (MQ=255)                       aCAAGGAAGAGAATGCGTTTAAAACGCTGGGACTGGATCAGGCGGAACATCATGGCGATTGCCGTACGGGTTTTACCGGTACCGGTCGCCATCGCCAGCAGGATCTCTTGTTACCCCTTGACGATTGCCTTTTCAACcgc                                                                                                              <  2:687492/140‑1 (MQ=255)                        cAAGGAAGAGAATGCGTTTAAAACGCTGGGACTGGATCAGGCGGAACATCATGGCGATTGCCGTACGGGTTTTACCGGTACCGGTCGCCATCGCCAGCAGGATCTCTTGTTACCCCTTGACGATTGCCTTTTCAACcgcg                                                                                                             <  2:632755/140‑1 (MQ=255)                            gAAGAGAATGCGTTTAAAACGCTGGGACTGGATCAGGCGGAACATCATGGCGATTGCCGTACGGGTTTTACCGGTACCGGTCGCCATCGCCAGCAGGATCTCTTGTTACCCCTTGACGATTGCCTTTTCAACCGCGCGGa                                                                                                         <  2:112308/140‑1 (MQ=255)                                 gaATGCGTTTAAAACGCTGGGACTGGATCAGGCGGAACATCATGGCGATTGCCGTACGGGTTTTACCGGTACCGGTCGCCATCGCCAGCAGGATCTCTTGTTACCCCTTGACGATTGCCTTTTCAACCGCGCGGACGGCa                                                                                                    <  2:512804/140‑1 (MQ=255)                                  aaTGCGTTTAAAACGCTGGGACTGGATCAGGCGGAACATCATGGCGATTGCCGTACGGGTTTTACCGGTACCGGTCGCCATCGCCAGCAGGATCTCTTGTTACCCCTTGACGATTGCCTTTTCAACCGCGCGGATGGCAt                                                                                                   >  2:561518/1‑140 (MQ=255)                                  aaTGCGTTTAAAACGCTGGGACTGGATCAGGCGGAACATCATGGCGATTGCCGTACGGGTTTTACCGGTACCGGTCGCCATCGCCAGCAGGATCTCTTGTTACCCCTTGACGATTGCCTTTTCAACCGCGCGGACGGCAt                                                                                                   >  2:656887/1‑140 (MQ=255)                                      cGTTTAAAACGCTGGGACTGGATCAGGCGGAACATCATGGCGATTGCCGTACGGGTTTTACCGGTACCGGTCGCCATCGCCAGCAGGATCTCTTGTTACCCCTTGACGATTGCCTTTTCAACCGCGCGGACGGc                                                                                                     <  2:838849/134‑1 (MQ=255)                                      cGTTTAAAACGCTGGGACTGGATCAGGCGGAACATCATGGCGATTGCCGTACGGGTTTTACCGGTACCGGTCGCCATCGCCAGCAGGATCTCTTGTTACCCCTTGACGATTGCCTTTTCAACCGCGCGGACGGc                                                                                                     >  1:838849/1‑134 (MQ=255)                                       gTTTAAAACGCTGGGACTGGATCAGGCGGAACATCATGGCGATTGCCGTACGGGTTTTACCGGTACCGGTCGCCATCGCCAGCAGGATCTCTTGTTACCCCTTGACGATTGCCTTTTCAACCGCGCGGACGGCATCTTcc                                                                                              >  2:643700/1‑140 (MQ=255)                                          tAAAACGCTGGGACTGGATCAGGCGGAACATCATGGCGATTGCCGTACGGGTTTTACCGGTACCGGTCGCCATCGCCAGCAGGATCTCTTGTTACCCCTTGACGATTGCCTTTTCAACCGCGCGGACGGCATCtt                                                                                                >  2:790443/1‑135 (MQ=255)                                          tAAAACGCTGGGACTGGATCAGGCGGAACATCATGGCGATTGCCGTACGGGTTTTACCGGTACCGGTCGCCATCGCCAGCAGGATCTCTTGTTACCCCTTGACGATTGCCTTTTCAACCGCGCGGACGGCATCtt                                                                                                <  1:790443/135‑1 (MQ=255)                                                  tGGGACTGGATCAGGCGGAACATCATGGCGATTGCCGTACGGGTTTTACCGGTACCGGTCGCCATCGCCAGCAGGATCTCTTGTTACCCCTTGACGATTGCCTTTTCAACCGCGCGGACGGCATCTTCCTGATAATAACg                                                                                   <  1:81951/140‑1 (MQ=255)                                                        tGGATCAGGCGGAACATCATGGCGATTGCCGTACGGGTTTTACCGGTACCGGTCGCCATCGCCAGCAGGATCTCTTGTTACCCCTTGACGATTGCCTTTTCAACCGCGCGGACGGCATCTTCCTGATAATAACGCAGGcc                                                                             <  2:651133/140‑1 (MQ=255)                                                                gCGGAACATCATGGCGATTGCCGTACGGGTTTTACCGGTACCGGTCGCCATCGCCAGCAGGATCTCTTGTTACCCCTTGACGATTGCCTTTTCAACCGCGCGGACGGCATCTTCCTGATAATAACGCAGGCCCAgctcgc                                                                     >  1:90872/1‑140 (MQ=255)                                                                 cGGAACATCATGGCGATTGCCGTACGGGTTTTACCGGTACCGGTCGCCATCGCCAGCAGGATCTCTTGTTACCCCTTGACGATTGCCTTTTCAACCGCGCGGACGGCATCTTCCTGATAATAACGCAGGCCCAgctcgct                                                                    <  2:769521/140‑1 (MQ=255)                                                                 cGGAACATCATGGCGATTGCCGTACGGGTTTTACCGGTACCGGTCGCCATCGCCAGCAGGATCTCTTGTTACCCCTTGACGATTGCCTTTTCAACCGCGCGGACGGCATCTTCCTGATAATAACGCAGGCCCAgctcgct                                                                    <  1:656887/140‑1 (MQ=255)                                                                 cGGAACATCATGGCGATTGCCGTACGGGTTTTACCGGTACCGGTCGCCATCGCCAGCAGGATCTCTTGTTACCCCTTGACGATTGCCTTTTCAACCGCGCGGACGGCATCTTCCTGATAATAACGCAGGCCCAgctcgct                                                                    <  1:706034/140‑1 (MQ=255)                                                                    aaCATCATGGCGATTGCCGTACGGGTTTTACCGGTACCGGTCGCCATCGCCAGCAGGATCTCTTGTTACCCCTTGACGATTGCCTTTTCAACCGCGCGGACGGCATCTTCCTGATAATAACGCAGGCCCAGCTCGCTCAt                                                                 >  2:614656/1‑140 (MQ=255)                                                                          atGGCGATTGCCGTACGGGTTTTACCGGTACCGGTCGCCATCGCCAGCAGGATCTCTTGTTACCCCTTGACGATTGCCTTTTCAACCGCGCGGACGGCATCTTCCTGATAATAACGCAGGCCCAGCTCGCTCATGCCAgg                                                           >  2:544579/1‑140 (MQ=255)                                                                              cGATTGCCGTACGGGTTTTACCGGTACCGGTCGCCATCGCCAGCAGGATCTCTTGTTACCCCTTGACGATTGCCTTTTCAACCGCGCGGACCGCATCTTCCTGAAAATAACGCAGGCGCAGCCCGCTCATGCCAGGGTTa                                                       >  2:300369/1‑140 (MQ=255)                                                                                aTTGCCGTACGGGTTTTACCGGTACCGGTCGCCATCGCCAGCAGGATCTCTTGTTACCCCTTGACGATTGCCTTTTCAACCGCGCGGACGGCATCTTCCTGATAATAACGCAGGCCCAGCTCGCTCATGCCAGGGTTATc                                                     >  2:101617/1‑140 (MQ=255)                                                                                 ttGCCGTACGGGTTTTACCGGTACCGGTCGCCATCGCCAGCAGGATCTCTTGTTACCCCTTGACGATTGCCTTTTCAACCGCGCGGACGGCATCTTCCTGATAATAACGCAGGCCCAgctcgct                                                                    <  1:757194/124‑1 (MQ=255)                                                                                 ttGCCGTACGGGTTTTACCGGTACCGGTCGCCATCGCCAGCAGGATCTCTTGTTACCCCTTGACGATTGCCTTTTCAACCGCGCGGACGGCATCTTCCTGATAATAACGCAGGCCCAgctcgct                                                                    >  2:757194/1‑124 (MQ=255)                                                                                       tACGGGTTTTACCGGTACCGGTCGCCATCGCCAGCAGGATCTCTTGTTACCCCTTGACGATTGCCTTTTCAACCGCGCGGACGGCATCTTCCTGATAATAACGCAGGCCCAGCTCGCTCATGCCAGGGTTATCGGCAAAc                                              >  1:204062/1‑140 (MQ=255)                                                                                                           gtCGCCATCGCCAGCAGGATCTCTTGTTACCCCTTGACGATTGCCTTTTCAACCGCGCGGACGGCATCTTCCTGATAATAACGCAGGCCCAGCTCGCTCATGCCAGGGTTATCGGCAAACCACTGATTCTGTTTTTGCgg                          >  1:707475/1‑140 (MQ=255)                                                                                                           gtCGCCATCGCCAGCAGGATCTCTTGTTACCCCTTGACGATTGCCTTTTCAACCGCGCGGACGGCATCTTCCTGATAATAACGCAGGCCCAGCTCGCTCATGCCAGGGTTATCGGCAAACCACTGATTCTGTTTTTGCgg                          >  2:777514/1‑140 (MQ=255)                                                                                                           gtCGCCATCGCCAGCAGGATCTCTTGTTACCCCTTGACGATTGCCTTTTCAACCGCGCGGACGGCATCTTCCTGATAATAACGCAGGCCCAGCTCGCTCATGCCAGGGTTATCGGCAAACCACTGATTCTGTTTTTGCgg                          >  1:837877/1‑140 (MQ=255)                                                                                                             cgccatcgccaGCAGGATCTCTTGTTACCCCTTGACGATTGCCTTTTCAACCGCGCGGACGGCATCTTCCTGATAATAACGCAGGCCCAGCTCGCTCATGCCAGGGTTATCGGCAAACCACTGATTCTGTTTTTGCGGtt                        <  2:10373/140‑1 (MQ=255)                                                                                                                catcgccaGCAGGATCTCTTGTTACCCCTTGACGATTGCCTTTTCAACCGCGCGGACGGCATCTTCCTGATAATAACGCAGGCCCAGCTCGCTCATGCCAGGTTTATCGGGAAAGCAGTGATTCAGTTTTTGCGGTTCGc                     >  1:235941/1‑140 (MQ=255)                                                                                                                catcgccaGCAGGATCTCTTGTTACCCCTTGACGATTGCCTTTTCAACCGCGCGGACGGCATCTTCCTGATAATAACGCAGGCCCAGCTCGCTCATGCCAGGGTTATCGGTAAACCAATGATTCCGTTTTTGCGGTTCGc                     >  1:354279/1‑140 (MQ=255)                                                                                                                                  ttgttaCCCCTTGACGATTGCCTTTTCAACCGCGCGGACGGCATCTTCCTGATAATAACGCAGGCCCAGCTCGCTCATGCCAGGGTTATCGGCAAACCACTGATTCTGTTTTTGCGGTTCGCTGCCAAGCATTTCcagca   <  1:777514/140‑1 (MQ=255)                                                                                                                                   tgttaCCCCTTGACGATTGCCTTTTCAACCGCGCGGACGGCATCTTCCTGATAATAACGCAGGCCCAGCTCGCTCATGCCAGGGTTATCGGCAAACCACTGATTCTGTTTTTGCGGTTCGCTGCCGAGCATTTCcagcag  <  1:614656/140‑1 (MQ=255)                                                                                                                                       |                                                                                                                                        TTCGCCAAGAGAACGGCGGTCGACAAGGAAGAGAATGCGTTTAAAACGCTGGGACTGGATCAGGCGGAACATCATGGCGATTGCCGTACGGGTTTTACCGGTACCGGTCGCCATCGCCAGCAGGATCTCTTGTTGCCCCTTGACGATTGCCTTTTCAACCGCGCGGACGGCATCTTCCTGATAATAACGCAGGCCCAGCTCGCTCATGCCAGGGTTATCGGCAAACCACTGATTCTGTTTTTGCGGTTCGCTGCCGAGCATTTCCAGCAG  >  NC_000913/4585238‑4585507

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