Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	NC_000913	2,042,171	G→A	intergenic (+146/‑197)	zinT → / → yodB	metal‑binding protein ZinT/putative cytochrome b561 YodB

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	2,042,171	0	G	A	100.0%	87.0 / NA	29	intergenic (+146/‑197)	zinT/yodB	metal‑binding protein ZinT/putative cytochrome b561 YodB
Reads supporting (aligned to +/- strand):  ref base G (0/0);  new base A (14/15);  total (14/15)

CAGCAATCTTTATTAAATGAAATGGAAAACTGGCCAACGTATTATCCATATCAGTTGAGTAGCGAAGAAGTGGTCGAGGAAATGATGTCTCATTGAGTATTCTCATGATAACGCCTCGATGCCGCTTTAGTAAGTTATCATAACTGCCACTGGTCCATCCACAAACGCCACTGAACGCAAGCTAGCTACAGACACGCTCATCACTATGACGTGTCTGTATATTAATAAGCTAACCCGCATTGAGTTAACCAATAACGGATTCCATACACAATACGGCCAACAGAGAAAGATACCTGTGCTCACGCCATTGCTTATATTGGCTGGTTACAATATGCACTATCAATTTTTTAAATAACAAATTATTAAAACCACATCTGTATAAGGAATTTTTAAGGTTCGTGGGATAGCTTGACTGTGAAAATCACAGGAGCTACAAAAATGAACCGATTCTCAAAAACTCAA  >  NC_000913/2041930‑2042391                                                                                                                                                                                                                                                  |                                                                                                                                                                                                                             caGCAATCTTTATTAAATGAAATGGAAAACTGGCCAACGTATTATCCATATCAGTTGAGTTGCGAAGAAGGGGTCGAGGAAATGATGTCTCATTGAGTATTCCCATGATAACGCCTCGATGCCGTTTTAGTAAGTTATCATACCTGCCACTGGTCCATCCACAAACGCCACTGAACGCAAGCTAGCTACAGACACGCTCATCACTATGACGTGTCTGTATATTAATAAGCTAACCCGCATTAAGTTAACCa                                                                                                                                                                                                                     <  1:422249/251‑1 (MQ=255)     aaTCTTTATTAAATGAAATGGAAAACTGGCCAACGTATTATCCATATCAGTTGAGTAGCGAAGAAGTGGTAGAGGAAATGATGTCTCATTGAGTATTCTCATGATAACGCCTCGATGCCGCTTTAGTAAGTTATCATAACTGCCACTGGTCCATCCACAAACGCCACTGAACGCAAGCTAGCTACAGACACGCTCATCACTATGACGTGTCTGTATATTAATAAGCTAACCCGCATTAAGTTATCCAATaa                                                                                                                                                                                                                 <  1:3692/251‑1 (MQ=255)              taAATGAAATGGAAAACTGGCCAACGTATTATCCATATCAGTTGAGTAGCGAAGAAGTGGTCGAGGAAATGATGTCTCATTGAGTATTCTCATGATAACGCCTCGATGCCGCTTTAGTAAGTTATCATAACTGCCACTGGTCCATCCACAAACGCCACTGAACGCAAGCTAGCTACAGACACGCTCATCACTATGACGTGTCTGTATATTAATAAGCTAACCCGCATTAAGTTAACCAATAACGGATTcc                                                                                                                                                                                                         <  1:356302/250‑1 (MQ=255)                   gaaatgGAAAACTGGCCAACGTATTATCCATATCAGTTGAGTAGCGAAGAAGTGGTCGAGGAAATGATGTCTCATTGAGTATTCTCATGATAACGCCTCGATGCCGCTTTAGTAAGTTATCATAACTGCCACTGGTCCATCCACAAACGCCACTGAACGCAAGCTAGCTACAGACACGCTCATCACTATGACGTGTCTGTATATTAATAAGCTAACCCGCATTAAGTTAACCAATAACGGATTCCATacac                                                                                                                                                                                                   >  1:362540/1‑251 (MQ=255)                   aaaatgGAAAACTGGCCAACGTAGTAGCTATATAGGTTGATTAGCGAGGAAGTGGTCGAGGAAATGATGTCTCCGTGAGTATTCCCATGAAAACGCCTCGAGGCCGCTTTAGTAAGTTATCATAACTGCCACTGGTCCATCCACAAACGCCACTGAACGCAAGCTAGCTACAGACACGCTCATCACTATGACGTGTCCGTATATTAATAAGTTAACCCGCATTAAGTTAACCAATACCGGATTCCATacg                                                                                                                                                                                                    <  2:195755/249‑2 (MQ=255)                            aaCTGGCCAACGTATTATCCATATCAGTTGAGTAGCGAAGAAGTGGTCGAGGAAATGATGTCTCATTGAGTATTCTCATGATAACGCCTCGATGCCGCTTTAGTAAGTTATCATAACTGCCACTGGTCCATCCACAAACGCCACTGAACGCAAGCTAGCTACAGACACGCTCATCACTATGACGTGTCTGTATATTAATAAGCTAACCCGCATTAAGTTAACCAATAACGGATTCCATACACAATACGGcc                                                                                                                                                                                          >  2:232963/1‑251 (MQ=255)                             aCTGGCCAACGTATTATCCCTATCAGTTGAGTAGCGAAGAAGTGGTCGAGGAAATGATGTCTCATTGAGTATTCTCATGATAACGCCTCGATGCCGCTTTAGTAAGTTATCATAACTGCCACTGGTCCATCCACAAACGCCACTGAACGCAAGCTAGCTACAGACACGCTCATCACTATGACGTGTCTGTATATTAATAAGCTAACCCGCATTAAGTTAACCAATAACGGATTCCATACACAATACGGcc                                                                                                                                                                                          >  1:379041/1‑250 (MQ=255)                                       gcattatCCATATCAGTTGAGTAGCGAAGAAGTGGTCGAGGAAATGATGTCTCATTGAGTATTCTCATGATAACGCCTCGATGCCGCTTTAGTAAGTTATCATAACTGCCACTGGTCCATCCACAACCGCCACTGAACGAAAGCTAGCTACAGACACGCTCATCACTATGACGTGTCTGTATATTAATAAGCTAACCCGCGTTAAGTTAACCAATCACGGATTCCATAAACAATACGGCCAACAGAGaaa                                                                                                                                                                                >  1:222003/3‑250 (MQ=255)                                                     aGTTGAGTAGCGAAGAAGTGGTCGAGGAAATGATGTCTCATTGAGTATTCTCATGATAACGCCTCGATGCCGCTTTAGTAAGTTATCATAACTGCCACTGGTCCATCCACAAACGCCACTGAACGCAAGCTAGCTACAGACACGCTCATCACTATGACGTGTCTGTATATTAATAAGCTAACCCGCATTAAGTTAACCAATAACGGATTCCATACACAATACGGCCAACAGAGAAAGATACCTGTGCTCAc                                                                                                                                                                 >  1:470352/1‑251 (MQ=255)                                                                gaagaaGTGGTCGAGGAAATGATGTCTCATTGAGTATTCTCATGATAACGCCTCGATGCCGCTTTAGTAAGTTATCATAACTGCCACTGGTCCATCCACAAACGCCACTGAACGCAAGCTAGCTACAGACACGCTCATCACTATGACGTGTCTGTATATTAATAAGCTAACCCGCATTAAGTTAACCAATAACGGATTCCATACACAATACGGCCAACAGAGAAAGATACCTGTGCTCACGCCATTGCtt                                                                                                                                                       >  2:224144/1‑250 (MQ=255)                                                                        gggggggaaaaggttggccattttggTATTCACAGGATAACGCCTCGATGCCTCTTTAGAAAGTTACCATAACTGCCACTGGTCCATCCACAAACGCCACTGAACGCAAGCTAGCTACAGACACGCTCATCACTATGACGTGTCTGTATATTAATAAGCTAACCCGCATTAAGTTAACCAATAACGGATTCCATACTCAATACGGCCAACAGAGAAAGATACCTGTGCTCACGCCATTGCTTATATTGGc                                                                                                                                               <  2:34497/225‑1 (MQ=255)                                                                                  atgatgTCTCATTGAGTATTCTCATGATAACGCCTCGATGCCGCTTTAGTAAGTTATCATAACTGCCACTGGTCCATCCACAAACGCCACTGAACGCAAGCTAGCTACAGACACGCTCATCACTATGACGTGTCTGTATATTAATAAGCTAACCCGCATTAAGTTAACCAATAACGGATTCCATACACAATACGGCCAACAGAGAAAGATACCTGTGCTCACGCCATTGCTTATATTGGCTGGTTACAata                                                                                                                                    <  1:52076/251‑1 (MQ=255)                                                                                    tatgtctaattaggTATTCTCATGAAAACGCCCCGATGCCGCTTTAGTAAGTTACCATAACTGCCACTGGTCCATCCACAAACGCCACTGAACGCAAGCTAGCTACAGACACGCTCATCACTATGACGTGTCTGTATATTAATAAGCTAACCCGCATTAAGTTAACCAATAACGGATTCCATACACAATACGGCCAACAGAGAAAGATACCTGTGCTCACGCCATTGCTTATATTGGCTGGTTACAatat                                                                                                                                   <  2:362540/237‑1 (MQ=255)                                                                                     atgTCTCATTGAGTATTCTCATGATAACGCCTCGATGCCGCTTTAGTAAGTTATCATAACTGCCACTGGTCCATCCACAAACGCCACTGAACGCAAGCTAGCTACAGACACGCTCATCACTATGACGTGTCTGTATATTAATAAGCTAACCCGCATTAAGTTAACCAATAACGGATTCCATACACAATACGGCCAACAGAGAAAGATACCTGTGCTCACGCCATTGCTTATATTGGCTGGTTACAATATGc                                                                                                                                 <  1:430007/251‑1 (MQ=255)                                                                                          tcATTGAGTATTCTCATGATAACGCCTCGATGCCGCTTTAGTAAGTTATCATAACTGCCACTGGTCCATCCACAAACGCCACTGAACGCAAGCTAGCTACAGACACGCTCATCACTATGACGTGTCTGTATATTAATAAGCTAACCCGCATTAAGTTAACCAATAACGGATTCCATACACAATACGGCCAACAGATAAAGATACCTGTGCTCACGCCATTGCTTATATTGGCTGGTTACAATATGCACTAt                                                                                                                            >  1:142884/1‑251 (MQ=255)                                                                                              tGAGTATTCTCATGATAACGCCTCGATGCGGATTTAGTAAGATATCATAACTAACACATGTCCATCCACAACCGCCACTGCACGCAAGCTAGCTACAGACACGCTCATCACTATGACGTGTATGTATATTACTAAGCTAACCCGCATTAAGTTAACCAATAACGAATTCAATACACAATACGGACAACAAAGAAAGATACCTGTGCTCAAGCCATTGCATATATTGGCTGGATACAATATGTACTATAAAt                                                                                                                        >  1:316436/1‑251 (MQ=255)                                                                                                 gTATTCTCATGATAACGCCTCGATGCCGCTTTAGTAAGTTATCATAACTGCCACTGGTCCATCCACAAACGCCACTGAACGCAAGCTAGCTACAGACACGCTCATCACTATGACGTGTCTGTATATTAATAAGCTAACCCGCATTAAGTTAACCAATAACGGATTCCATACACAATACGGCCAACAGAGAAAGATACCTGTGCTCACGCCATTGCTTATATTGGCTGGTTACAATATGCACTATCAAtttt                                                                                                                     >  1:375933/1‑251 (MQ=255)                                                                                                                cGCCTCGATGCCGCTTTAGTAAGTTATCATAACTGCAACTGGTCCATCCACAAACGCCACTGAACGCAAGCTAGCTACAGACACGCTCATCACTATGACGTGTCTGTATATTAATAAGCTAACCCGCATTAAGTTAACCAATAACGGATTCCATACACAATACGGCCAACAGAGAAAGATACCTGTGCTCACGCCATTGCTTATATTGGCTGGTTACAATATGCACTATCAATTTTTTAAATAACAAatta                                                                                                      <  1:91390/251‑1 (MQ=255)                                                                                                                 gCCTCGATGCCGCTTTAGTAAGTTATCATAACTGCCACTGGTCCATCCACAAACGCCACTGAACGCAAGCTAGCTACAGACACGCTCATCACTATGACGTGTCTGTATATTAATAAGCTAACCCGCATTAAGTTAACCAATAACGGATTCCATACACAATACGGCCAACAGAGAAAGATACCTGTGCTCACGCCATTGCTTATATTGGCTGGTTACAATATGCACTATCAATTTTTTAAATAACAAattat                                                                                                     >  2:131464/1‑251 (MQ=255)                                                                                                                      gATGCCGCTTTAGTAAGTTATCAGAACTGCCACTGGTCCATCCACAAACGCCACTGAACGCAAGCTAGCTACAGACACGCTCATCACTATGACGTGTCTGTATATTAATAAGCTAACCCGCATTAAGTTAACCAATAACTGATTCCATACACAATACGGCCAACAGAGAAAGATACCTGTGCTCACGCCATTGCTTATATTGGCTGGTTACAATATGCACTATCAATTTTTTAAATAACAAATTATTaaaa                                                                                                <  1:230417/251‑1 (MQ=255)                                                                                                                                              aaCTGCCACTGGTCCAGCCACAAACGCCACTGAACGCAAGCTAGCTACAGACACGCTCATCACTATGACGTGTCTGTATATTAATAAGCTAACCCGCATTAAGTTAACCAATAACGGATTCCATACACAATACGGCCAACAGAGAAAGAAACATGTGCTCACGCCATTGCTTATATTGGCTGGTTACAATATGCACTATCAATTTTTTAAATAACAAATTATTAAAACCACATCTGTATAAGGAATTTTTa                                                                        <  1:102937/251‑1 (MQ=255)                                                                                                                                                cTGCCACTGGTCCATCCACAAACGCCACTGAACGCAAGCTAGCTACAGACACGCTCATCACTATGACGTGTCTGTATATTAATAAGCTAACCCGCCTTAAGTTAACCCATACCGGATTCCATACACAATACGGCCAACTGAGAAAGATACCTGTGCTCACGCCATTGCTTATATTGGCTGGTTACAATATGCACTATCAATTTTTTAAATAACAAATTATTAAAACCACATCTGTATCAGGAAttttaa                                                                        >  1:222219/1‑247 (MQ=255)                                                                                                                                                   ccACGGGTCCATCCACAAACGCTACTGAACGCAAGCTAGCTACATACACGCTCAGCACTATGACGTGACTGTATATTAATAAGCTAACCCGCATTAAGTTAACCAATAACGGATTCCATACACAATACTGCCAACAGAGAAAGATACCTTTTCTCACGCCATGGCTTATATTGTCTGGTTACAATATGCACTATCAATTTTTTAAATAACAAATTATTAAAACCACATCTGTATAAGGAAGTTGTAAGGtt                                                                   <  2:222003/251‑1 (MQ=255)                                                                                                                                                        ggTCCATCCACAAACGCCCCTTAACGCAAGCTAGCTACAGACACGCTCATCCCTATGACGTGTCTGTATAATAATAAGCTAACCCGCATTAAGTTAACCAATAACGGATTCCATACACAATACGGCCAACAGAGAAAGATACCTGTGCTTACGCCATTGCTTATATTGGCTGGTTACAATATGCACTATCAATTTTTTAAATAACAAATTATTAAAACCACATCTGTATAAGGAATTTTTAAGGTTCGTgg                                                              <  2:410118/251‑1 (MQ=255)                                                                                                                                                         gTCCATCCACAAACGCCACTGAACGCAAGCTAGCTACAGACACGCTCATCACTATGACGTGTCTGTATATTAATAAGCTAACCCGCATTAAGTTAACCAATAACGGATTCCATACACAATACGGCCAACAGAGAAAGATACCTGTGCTCACGCCATTGCTTATATTGGCTGGTTACAATATGCACTATCAATTTTTTAAATAACAAATTATTAAAACCACATCTGTATAAGGAATTTTTAATGTTCGTggg                                                             >  1:51787/1‑251 (MQ=255)                                                                                                                                                                             gAACGCAAGCTAGCTACAGACACGCTCATCACTATGACGTGTCTGTATATTAATAAGCTAACCCGCATTAAGTTAACCAATAACGGATTCCATACACAATACGGCCAACAGAGAAAGATACCTGTGCTCACGCCATTGCTTATATTGGCTGGTTACAATATGCACTATCAATTTTTTAAATAACAAATTATTAAAACCACATCTGTATAAGGAATTTTTAAGGTTCGTGGGATAGCTTGACTGTGAAAATc                                         <  2:118642/251‑1 (MQ=255)                                                                                                                                                                                                        atcaCTATGACGTGTCTGAATATTACTAAGCTAGCCCGCATTAAGTTAACCAATAAAGGATTCCATACACAATACGGCCAACAGAGAAAGATACCTGTGCTCACGCAATTGCTTATATTGGCTGGTTACAATATGCACTATCAATTTTTTAAATAACAAATTATTAAAACCACATCTGTATAAGGAATTTTTAAGGTTCGTGGGATAGCTTGACTGTGAAAATCCAAGGAGCTACAAAAATGAACCGATtc              >  1:176362/1‑251 (MQ=255)                                                                                                                                                                                                           aCTATGACGTGTCTGTATATTAATAAGCTAACCCGCATTAAGTTAACCAATAACGGATTCCATACACAATACGGCCAACAGAGAAAGATACCTGTGCTCACGCCATTGCTTATATTGGCTGGTTACAATATGCACTATCAATTTTTTAAATAACAAATTATTAAAACCACATCTGTATAAGGAATTTTTAAGGTTCGTGGGATAGCTTGACTGTGAAAATCACAGGAGATACAAAAATGAACCGATtctc            >  2:212914/1‑250 (MQ=255)                                                                                                                                                                                                                     gtCTGTATATTAATAAGCTAACCCGCATTAAGTTAACCAATAACGGATTCCATACACAATACGGCCAACAGATAAAGATACCTGTGCTCACGCCATTGCTTATATTGGCTGGTTACAATATGCACTATCAATTTTTTAAATAACAAATTATTAAAACCACATCTGTATAAGGAATTTTTAAGGTTCGTGGGATAGCTTGACTGTGAAAATCACAGGAGCTACAAAAATGAACCGATTCTCAAAAACTCaa  <  2:142884/250‑1 (MQ=255)                                                                                                                                                                                                                                                  |                                                                                                                                                                                                                             CAGCAATCTTTATTAAATGAAATGGAAAACTGGCCAACGTATTATCCATATCAGTTGAGTAGCGAAGAAGTGGTCGAGGAAATGATGTCTCATTGAGTATTCTCATGATAACGCCTCGATGCCGCTTTAGTAAGTTATCATAACTGCCACTGGTCCATCCACAAACGCCACTGAACGCAAGCTAGCTACAGACACGCTCATCACTATGACGTGTCTGTATATTAATAAGCTAACCCGCATTGAGTTAACCAATAACGGATTCCATACACAATACGGCCAACAGAGAAAGATACCTGTGCTCACGCCATTGCTTATATTGGCTGGTTACAATATGCACTATCAATTTTTTAAATAACAAATTATTAAAACCACATCTGTATAAGGAATTTTTAAGGTTCGTGGGATAGCTTGACTGTGAAAATCACAGGAGCTACAAAAATGAACCGATTCTCAAAAACTCAA  >  NC_000913/2041930‑2042391

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