Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	1	1,454,670	G→A	E31K (GAA→AAA)	unknown →	lipoprotein

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	1	1,454,670	0	G	A	100.0%	91.0 / NA	28	E31K (GAA→AAA)	unknown	lipoprotein
Reads supporting (aligned to +/- strand):  ref base G (0/0);  new base A (18/10);  total (18/10)

ATTTTATTGACAAATTGAAAGCGGTTTCCCCGCAATTACGCGCACTAGGCAAGAAATTACCGACTAAAGCGCAACTGGAAAAAGCCATGTTCAGTGGTGAGCTCTCCGCCCAATCCATTTAAGTCAGTTCTATATCCACAAAAACTGTGAATAACCAGGAGGATAATATGAAAAAATTTTTACTTATGAATACGTTAATCTCTGCTTGTTTTTTAAGCGCGT‑GCAACTCGTCACCGCCGGAACCGCCACAAATTAATGAAGACGCGCCGGCAATCTATTTAAATAATCAGGCGTACCGACAAACACCAACTGCAACAGTGCCAAAAAACGCCGGCGATCATAATGGGCAACCTTGGGTATATCAATATGTGAACCTAGGTAGAAGCGATTATGTGAATGAAACAGAAAAAGTTCGTTTTTTCTACTTTGCGCATCATGCAGATAGTATTGAAATCTACGGTCAGCCGGCGCGTACCGAAGCCTATAAA  >  1/1454413‑1454900                                                                                                                                                                                                                                                                   |                                                                                                                                                                                                                                       tttttataGACAAACTCAAAGCGGTATCTCCGCAATTGCGAGCATTAGGTAAAAAATTACCGACAAAAGAACAATTGGAAAAAGCTATGTTCAGTGGCGAGCTGTCAGTCCAATCAATTTA‑‑‑‑‑‑‑‑‑‑‑‑ATCCACAAAAACTGTGAATAACCAGGAGGATAATATGAAAAAAATTTTACTTATGAATATGTTAATCCCTGTTTGTTTTTTAAGCGCGT‑GCAGCTCGTCACCACCAGAACCGCCGCAAATTGATAAAGAc                                                                                                                                                                                                                                   <  1:497178/250‑1 (MQ=255)                                                                                                                  ctgtcagtccaataaattaaatcaacaAAAACTGTGAAAAACCAGGAGAAAAATAAAAAAAAAATTTTACTTATGAATATGTTAATCCCTGTTTCCTTTTAAGGCG‑GTGGCAGCTCGACACCACCAAAACCGCCGCAAATTGATAAAGACGCCCCAGTAATCCATTTAAATAATCAGGTGTACCGACAAACACCAACTGCAACAGTGCCAAAAAACGCCGGCGATCATAATGGGCAACCTTGGGTATATc                                                                                                                               <  1:408089/226‑1 (MQ=255)                                                                                                                           caatcaatttaatCCACAAAAACTGTGAATAACCAGGAGGATAATATGAAAAAAATTTTACTTATGAATATGTTAATACCTGTTTGCTTTTTAAGCGCGT‑GCAGCTCGTCACCACCAGAACCGCCGCAAATTGATAAAGACGCCCCAGTAATCCATTTAAATAATCAGGTGTACCGACAAACACCAACTGCAACAGTGCCAAAAAACGCCGGCGATCATAATGGGCAACCTTGGGTATATCAATAtg                                                                                                                         >  1:393289/12‑247 (MQ=255)                                                                                                                                                aCTGTGAATAACCAGGAGGATAATATGAAAAAAATTTTACTTATGAATATGTTAATACCTGTTTGCTTTTTAAGCGCGT‑GCAGCTCGTCACCACCAGAACCGCCGCAAATTGATAAAGACGCCCCAGTAATCCATTTAAATAATCAGGTGTACCGACAAACACCAACTGCAACAGTGCCAAAAAACGCCGGCGATCATAATGGGCAACCTTGGGTATATCAATATGTGAACCTAGGTAGAAGCGATTAtgt                                                                                                >  1:58895/1‑251 (MQ=255)                                                                                                                                                 cTGTGAATAACCAGGAGGATAATATGAAAAAAATTTTACTTATGAATATGTTAATACCTGTTTGCTTTTTAAGCGCGT‑GCAGCTCGTCACCACCAGAACCGCCGCAAATTGATAAAGACGCCCCAGTAATCCATTTAAATAATCAGGTGTACCGACAAACACCAACTGCAACAGTGCCAAAAAACGCCGGCGATCATAATGGGCAACCTTGGGTATATCAATATGTGAACCTAGGTAGAAGCGATTAtgtg                                                                                               >  1:398930/1‑251 (MQ=255)                                                                                                                                                            cAGGAGGATAATATGAAAAAAATTTTACTTATGAATATGTTAATACCTGTTTGCTTTTTAAGCGCGT‑GCAGCTCGTCACCACCAGAACCGCCGCAAATTGATAAAGACGCCCCAGTAATCCATTTAAATAATCAGGTGTACCGACAAACACCAACTGCAACAGTGCCAAAAAACGCCGGCGATCATAATGGGCAACCTTGGGTATATCAATATGTGAACCTAGGTAGAAGCGATTATGTGAATGAAACAGa                                                                                    >  1:22401/1‑251 (MQ=255)                                                                                                                                                                 ggaTAATATAAAAAAAATTTTACTTATAAATATGTTAATCCCTGTTTGCTTTTTAAGCGCGT‑GCAGCTCGTCCCCACCAGAACCGCCGCAAATTGATAAAGACGCCCCAGTAATCCATTTAAATAATCAGGTGTACCGACAAACACCAACTGCAACAGTGCCAAAAAACGCCGGCGATCATAATg                                                                                                                                                 <  1:461020/185‑1 (MQ=255)                                                                                                                                                                                    ttACTTATGAATATGTTAATACCTGTTTGCTTTTTAAGCGCGT‑GCAGCTCGTCACCACCAGAACCGCCGCAAATTGATAAAGACGCCCCAGTAATCCATTTAAATAATCAGGTGTACCGACAAACACCAACTGCAACAGTGCCAAAAAACGCCGGCGATCATAATGGGCAACCTTGGGTATATCAATATGTGAACCTAGGTAGAAGCGATTATGTGAATGAAACAGAAAAAGTTCGTTTTTTCTACTTTg                                                             >  1:94806/1‑250 (MQ=255)                                                                                                                                                                                       cTTATGAATATGTTAATACCTGTTTGCTTTTTAAGCGCGT‑GCAGCTCGTCACCACCAGAACCGCCGCAAATTGATAAAGACGCCCCAGTAATCCATTTAAATAATCAGGTGTACCGACAAACACCAACTGCAACAGTGCCAAAAAACGCCGGCGATCATAATGGGCAACCTTGGGTATATCaa                                                                                                                             >  1:224595/1‑183 (MQ=255)                                                                                                                                                                                       cTTATGAATATGTTAATACCTGTTTGCTTTTTAAGCGCGT‑GCAGCTCGTCACCACCAGAACCGCCGCAAATTGATAAAGACGCCCCAGTAATCCATTTAAATAATCAGGTGTACCGACAAACACCAACTGCAACAGTGCCAAAAAACGCCGGCGATCATAATGGGCAACCTTGGGTATATCAATATGTGAACCTAGGTAGAAGCGATTATGTGAATGAAACAGAAAAAGTTCGTTTTTTCTACTTTGCGca                                                         >  1:892947/1‑251 (MQ=255)                                                                                                                                                                                           tGAATATGTTAATACCTGTTTGCTTTTTAAGCGCGT‑GCAGCTCGTCACCACCAGAACCGCCGCAAATTGATAAAGACGCCCCAGTAATCCATTTAAATAATCAGGTGTACCGACAAACACCAACTGCAACAGTGCCAAAAAACGCCGGCGATCATAATGGGCAACCTTGGGTATATCAATATGTGAACCTAGGTAGAAGCGATTATGTGAATGAAACAGAAAAAGTTCGTTTTTTCTACtt                                                               >  1:499998/1‑241 (MQ=255)                                                                                                                                                                                               tatgttaatacctgtttGTTTTTTAAGCGCGT‑GCAGCTCGTCACAACCAGAACCGCCGCAAATTGATAAAGACGCCCCAGTAATCCATTTAAATAATCAGGTGTACCGACAAACACCAACTGCAACAGTGCCAAAAAACGCCGGCGATCATAATGGGCAACCTTGGGTATATCAATATGTGAACCTAGGTAGAAGCGATTATGTGAATGAAACAGAAAAAGTTCGTTTTTTCTACTTTGCGCATCATGCAg                                                 <  1:451406/236‑1 (MQ=255)                                                                                                                                                                                                  gttaatccctttttGTTTTTTAAGCGCGT‑GCAGCTCGCCCCCACCAAAACCGCCGCAAATTGATAAAGACGCCCCAGTAATCCATTTAAATAATCAGGTGTACCGACAAACACCAACTGCAACAGTGCCAAAAAACGCCGGCGATCATAATGGGCAACCTTGGGTATATCAATAtg                                                                                                                         <  1:705712/164‑1 (MQ=255)                                                                                                                                                                                                  gttaatacctgtttgctttttAAGCGCGT‑GCAGCTCGTCACCACCAGAACCGCCGCAAATTGATAAAGACGCCCCAGTAATCCATTTAAATAATCAGGTGTACCGACAAACACCAACTGCAACAGTGCCAAAAAACGCCGGCGATCATAATGGGCAACCTTGGGt                                                                                                                                    >  1:670508/17‑165 (MQ=255)                                                                                                                                                                                                  gttaatacctgtttgctttttAAGCGCGT‑GCAGCTCGTCACCACCAGAACCGCCGCAAATTGATAAAGACGCCCCAGTAATCCATTTAAATAATCAGGTGTACCGACAAACACCAACTGCAACAGTGCCAAAAAACGCCGGCGATCATAATGGGCAACCTTGGGt                                                                                                                                    >  1:603154/17‑165 (MQ=255)                                                                                                                                                                                                  gttaatacctgtttgctttttAAGCGCGT‑GCAGCTCGTCACCACCAGAACCGCCGCAAATTGATAAAGACGCCCCAGTAATCCATTTAAATAATCAGGTGTACCGACAAACACCAACTGCAACAGTGCCAAAAAACGCCGGCGATCATAATGGGCAACCTTGGGTATATCAATAtg                                                                                                                         >  1:588047/17‑176 (MQ=255)                                                                                                                                                                                                  gttaatacctgtttgctttttAAGCGCGT‑GCAGCTCGTCACCACCAGAACCGCCGCAAATTGATAAAGACGCCCCAGTAATCCATTTAAATAATCAGGTGTACCGACAAACACCAACTGCAACAGTGCCAAAAAACGCCGGCGATCATAATGGGCAACCTTGGGTATATCAATATGTGAACCTAGGTAGAAGCGATTATGTGAATGAAACAGAAAAAGTTCGTTTTTTCTACTTTGCGCATCATGCAGATa                                              >  1:492869/17‑251 (MQ=255)                                                                                                                                                                                                  gttaatacctgtttgctttttAAGCGCGT‑GCAGCTCGTCACCACCAGAACCGCCGCAAATTGATAAAGACGCCCCAGTAATCCATTTAAATAATCAGGTGTACCGACAAACACCAACTGCAACAGTGCCAAAAAACGCCGGCGATCATAATGGGCAACCTTGGGTATATCAATATGTGAACCTAGGTAGAAGCGATTATGTGAATGAAACAGAAAAAGTTCGTTTTTTCTACTTTGCGCATCATGCAGATa                                              >  1:265026/17‑251 (MQ=255)                                                                                                                                                                                                  aggtgtataagagacattttaAAGCGCGT‑GCAGCTCGTCCCCACCAGAACCGCCGCAAATTGATAAAGACGCCCCAGTAATCCATTTAAATAATCAGGTGTACCGACAAACACCAACTGCAACAGTGCCAAAAAACGCCGGCGATCATAATGGGCAACCTTGGGTATATCAATATGTGAACCTAGGTAGAAGCGATTATGTGAATGAAACAGAAAAAGTTCGTTTTTTCTACTTTGCGCATCATGCAGAt                                               <  1:417446/234‑1 (MQ=255)                                                                                                                                                                                                    taatacctgtttgctttttAAGCGCGT‑GCAGCTCGTCACCACCAGAACCGCCGCAAATTGATAAAGACGCCCCAGTAATCCATTTAAATAATCAGGTGTACCGACAAACACCAACTGCAACAGTGCCAAAAAACGCCGGCGATCATAATGGGCAACCTTGGGTATATCAATATGTGAACCTAGGTAGAAGCGATTAtg                                                                                                 >  1:86519/15‑198 (MQ=255)                                                                                                                                                                                                                gctttttAAGCGCGT‑GCAGCTCGTCACCACCAGAACCGCCGCAAATTGATAAAGACGCCCCAGTAATCCATTTAAATAATCAGGTGTACCGACAAACACCAACTGCAACAGTGCCAAAAAACGCCGGCGATCATAATGGGCAACCTTGGGt                                                                                                                                    >  1:834182/3‑151 (MQ=255)                                                                                                                                                                                                                   ttttAAGCGCGT‑GCAGCTCGTCACCACCAGAACCGCCGCAAATTGATAAAGACGCCCCAGTAATCCATTTAAATAATCAGGTGTACCGACAAACACCAACTGCAACAGTGCCAAAAAACGCCGGCGATCATAATGGGCAACCTTGGGTATATCAATATGTGAACCTAGGTAGAAGCGATTATGTGAATGAAACAGAAAAAGTTCGTTTTTTCTACTTTGCGcatc                                                       >  1:878525/1‑225 (MQ=255)                                                                                                                                                                                                                   ttttAAGCGCGT‑GCAGCTCGTCACCACCAGAACCGCCGCAAATTGATAAAGACGCCCCAGTAATCCATTTAAATAATCAGGTGTACCGACAAACACCAACTGCAACAGTGCCAAAAAACGCCGGCGATCATAATGGGCAACCTTGGGTATATCAATATGTGAACCTAGGTAGAAGCGATTATGTGAATGAAACAGAAAAAGTTCGTTTTTTCTACTTTGCGCATCATGCAGATAGTATTGAAATCTACGGt                             >  1:80966/1‑251 (MQ=255)                                                                                                                                                                                                                   tttaaaGCGCGT‑GCAGCTCGTCACCACCAGAACCGCCGCAAATTGATAAAGACGCCCCAGTAATCCATTTAAATAATCAGGTGTACCGACAAACACCAACTGCAACAGTGCCAAAAAACGCCGGCGATCATAATGGGCAACCTTGGGTATATCAata                                                                                                                           <  1:96786/153‑1 (MQ=255)                                                                                                                                                                                                                                        tcaccaccaaaaCCGCCGCAAATTGATAAAGACGCCCCAGTAATCCATTTAAATAATCAGGTGTACCGACAAACACCAACTGCAACAGTGCCAAAAAACGCCGGCGATCATAATGGGCAACCTTGGGTATATCAATATGTGAACCTAGGTAGAAGCGATTATGTGAATGAAACAGAAAAAGTTCGTTTTTTCTACTTTGCGCATCATGCAGATAGTATTGAAATCTACGGTCAGCCGGCGCGTACCGAAg          <  1:562168/240‑1 (MQ=255)                                                                                                                                                                                                                                           ccaccagAACCGCCGCAAATTGATAAAGACGCCCCAGTAATCCATTTAAATAATCAGGTGTACCGACAAACACCAACTGCAACAGTGCCAAAAAACGCCGGCGATCATAATGGGCAACCTTGGGTATATCAATATGTGAACCTAGGTAGAAGCGATTATGTGAATGAAACAGAAAAAGTTCGTTTTTTCTACTTTGCGCATCATGCAGATa                                              >  1:465040/7‑211 (MQ=255)                                                                                                                                                                                                                                               caaaaCCGCCGCAAATTGATAAAGACGCCCCAGTAATCCATTTAAATAATCAGGTGTACCGACAAACACCAACTGCAACAGGGCCAAAAAACGCCGGCGATCATAATGGGCAACCTTGGGTATATCAATATGTGAACCTAGGTAGAAGCGATTATGTGAATGAAACAGAAAAAGTTCGTTTTTTCTACTTTGCGCATCATGCAGATAGTATTGAAATCTACGGTCAGCCGGCGCGTACCGAAGCCTATaaa  <  1:456126/248‑1 (MQ=255)                                                                                                                                                                                                                                                  aaCCGCCGCAAATTGATAAAGACGCCCCAGTAATCCATTTAAATAATCAGGTGTACCGACAAACACCAACTGCAACAGTGCCAAAAAACGCCGGCGATCATAATGGGCAACCTTGGGTATATCAATATGTGAACCTAGGTAGAAGCGATTATGTGAATGAAACAGAAAAAGTTCGTTTTTTCTACTTTGCGCATCATGCAGATAGTATTGAAATCTACGGTCAGCCGGCGCGTACCGAAGCCTATaaa  <  1:46102/248‑1 (MQ=255)                                                                                                                                                                                                                                                                   |                                                                                                                                                                                                                                       ATTTTATTGACAAATTGAAAGCGGTTTCCCCGCAATTACGCGCACTAGGCAAGAAATTACCGACTAAAGCGCAACTGGAAAAAGCCATGTTCAGTGGTGAGCTCTCCGCCCAATCCATTTAAGTCAGTTCTATATCCACAAAAACTGTGAATAACCAGGAGGATAATATGAAAAAATTTTTACTTATGAATACGTTAATCTCTGCTTGTTTTTTAAGCGCGT‑GCAACTCGTCACCGCCGGAACCGCCACAAATTAATGAAGACGCGCCGGCAATCTATTTAAATAATCAGGCGTACCGACAAACACCAACTGCAACAGTGCCAAAAAACGCCGGCGATCATAATGGGCAACCTTGGGTATATCAATATGTGAACCTAGGTAGAAGCGATTATGTGAATGAAACAGAAAAAGTTCGTTTTTTCTACTTTGCGCATCATGCAGATAGTATTGAAATCTACGGTCAGCCGGCGCGTACCGAAGCCTATAAA  >  1/1454413‑1454900

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