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

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	1	778,237	0	G	A	100.0%	122.8 / NA	37	E31K (GAA→AAA)	unknown	lipoprotein
Reads supporting (aligned to +/- strand):  ref base G (0/0);  new base A (25/12);  total (25/12)

ATTACCGACTAAAGCGCAACTGGAAAAAGCCATGTTCAGTGGTGAGCTTTCCGCCCAATCTATTTAATCCACAAAAACTGTGAATAACCAGGAGGATAATATGAAAAAATTTTTAGTTATGAATACGTTAATCCCTGTTTGCTTTTTAAGCGCGT‑GCAGCTCGTCAACGCCGGAACCGCCACAAATCGATGAAGACGCCCCGGCAATTTATTTAAATAATCAGGTGTACCGACAAACACCAACTGCAACAGTGCCAAAAAACGCCGGCGATCATAACGGGCAACCTTGGGTATATCAATATGTGAACTTAGATAGAAACGATTATGTGAATGAAACAGAAAAAGTTCGTTTTTTCTACTTTGCTCATCATGCAGATAGTATTGAAATCTACGGTCAGCCAGCGCGTACCGAAGCCTATAAA  >  1/778047‑778467                                                                                                                                                                                                |                                                                                                                                                                                                                                       attaccgaaaaaagaacaatgggAAAAAGTTATGTTCAGTGGCGAGCTGTCAGTCCAATAAATTTAATCCAAAAAAACTGTGAAAAACCAGGAGGATAATATGAAAAAAATTTTACTTATGAATATGTTAATCCCTGTTTGCTTTTTAAGCGCGT‑GCAGCTCGTCACCACCAGAACCGCCGCAAATTGATAAAGACGCCCCAGTAATCCATTTAAATAATCAGGTGTACCGACAAACACCAACTGCAACAg                                                                                                                                                                            <  1:896182/230‑1 (MQ=255)              gaacaattGGAAAAAGCTATGTTCAGTGGCGAGCTGTCAGTCCAATCAATTTAATCCACAAAAACTGTGAATAACCAGGAGGATAATATGAAAAAAATTTTACTTATGAATATGTTAATACCTGTTTGCTTTTTAAGCGCGT‑GCAGCTCGTCACCACCAGAACCGCCGCAAATTGATAAAGACGCCCCAGTAATCCATTTAAATAATCAGGTGTACCGACAAACACCAACTGCAACAGTGCCAAAAAACg                                                                                                                                                                >  1:342336/8‑250 (MQ=255)                   attGGAAAAAGCTATGTTCAGTGGCGAGCTGTCAGTCCAATCAATTTAATCCACAAAAACTGTGAATAACCAGGAGGATAATATGAAAAAAATTTTACTTATGAATATGTTAATACCTGTTTGCTTTTTAAGCGCGT‑GCAGCTCGTCACCACCAGAACCGCCGCAAATTGATAAAGACGccccagtaa                                                                                                                                                                                                                         >  1:882833/3‑183 (MQ=255)                           aaGCTATGTTCAGTGGCGAGCTGTCAGTCCAATCAATTTAATCCACAAAAACTGTGAATAACCAGGAGGATAATATGAAAAAAATTTTACTTATGAATATGTTAATACCTGTTTGCTTTTTAAGCGCGT‑GCAGCTCGTCACCACCAGAACCGCCGCAAATTGATAAAGACGccccagtaatcc                                                                                                                                                                                                                      >  1:8297/1‑175 (MQ=255)                              ctaTGTCCAGTGGCGAGCTGTCAGTCCAAACAATTTAACCCAAAAAAACGGTGAATAACAAGGAGAAAAATAAAAAAAAAATTTTACTTATAAATATTTTAATCCCTGTTTGTTTTTAAAGCGCGT‑GCAGCTCGTCCCCACCAGAACCGCCGCAAATTGATAAAGACGCCCCAGTAATCCATTTAAATAATCAGGTGTACCGACAAACACCAACTGCAACAGTGCCAAAAAACGCCGGCGATCATAATggg                                                                                                                                               <  1:880328/249‑1 (MQ=255)                                  gTTCAGTGGCGAGCTGTCAGTCCAATCAATTTAATCCACAAAAACTGTGAATAACCAGGAGGATAATATGAAAAAAATTTTACTTATGAATATGTTAATACCTGTTTGCTTTTTAAGCGCGT‑GCAGCTCGTCACCACCAGAACCGCCGCAAATTGATAAAGACGCCCCAGTAATCCATTTAAATAATCAGGTGTACCGACAAACACCAACTGCAACAGTGCCAAAAAACGCCGGCGATCATAATGGGCaa                                                                                                                                            >  1:145306/1‑250 (MQ=255)                                  gTTCAGTGGCGAGCTGTCAGTCCAATCAATTTAATCCACAAAAACTGTGAATAACCAGGAGGATAATATGAAAAAAATTCTACTTATGAATATGTTAATACCTGTTTGCTTTTTAAGCGCGT‑GCAGCTCGTCACCACCAGAACCGCCGCAAATTGATAAAGACGCCCCAGTAATCCATTTAAATAATCAGGTGTACCGACAAACACCAACTGCAACAGTGCCAAAAAACGCCGGCGATCATAATGGGCaa                                                                                                                                            >  1:638072/1‑250 (MQ=255)                                               ctgtcagtccaataaattaaaTCAACAAAAACTGTGAAAAACCAGGAGAAAAATAAAAAAAAAATTTTACTTATGAATATGTTAATCCCTGTTTCCTTTTAAGGCG‑GTGGCAGCTCGACACCACCAAAACCGCCGCAAATTGATAAAGACGCCCCAGTAATCCATTTAAATAATCAGGTGTACCGACAAACACCAACTGCAACAGTGCCAAAAAACGCCGGCGATCATAATGGGCAACCTTGGGTATATc                                                                                                                               <  1:408089/231‑1 (MQ=255)                                                        caatcaatTTAATCCACAAAAACTGTGAATAACCAGGAGGATAATATGAAAAAAATTTTACTTATGAATATGTTAATACCTGTTTGCTTTTTAAGCGCGT‑GCAGCTCGTCACCACCAGAACCGCCGCAAATTGATAAAGACGCCCCAGTAATCCATTTAAATAATCAGGTGTACCGACAAACACCAACTGCAACAGTGCCAAAAAACGCCGGCGATCATAATGGGCAACCTTGGGTATATCAATAtg                                                                                                                         >  1:393289/1‑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)                                                                                            gaggaTAATATAAAAAAAATTTTACTTATGAATATGTTAATCCCTGTTTGCTTTTTAAGCGCGT‑GCAGCTCGTCACCACCAGAACCGCCGCAAATTGATAAAGACGCCCCAGTAATCCATTTAAATAATCAGGTGTACCGACAAacac                                                                                                                                                                                        <  1:431333/148‑1 (MQ=255)                                                                                              ggaTAATATAAAAAAAATTTTACTTATAAATATGTTAATCCCTGTTTGCTTTTTAAGCGCGT‑GCAGCTCGTCCCCACCAGAACCGCCGCAAATTGATAAAGACGCCCCAGTAATCCATTTAAATAATCAGGTGTACCGACAAACACCAACTGCAACAGTGCCAAAAAACGCCGGCGATCATaatg                                                                                                                                                 <  1:461020/185‑3 (MQ=255)                                                                                                                 ttacttATGAATATGTTAATACCTGTTTGCTTTTTAAGCGCGT‑GCAGCTCGTCACCACCAGAACCGCCGCAAATTGATAAAGACGCCCCAGTAATCCATTTAAATAATCAGGTGTACCGACAAACACCAACTGCAACAGTGCCAAAAAACGCCGGCGATCATAATGGGCAACCTTGGGTATATCAATATGTGAACCTAGGTAGAAGCGATTATGTGAATGAAACAGAAAAAGTTCGTTTTTTCTACTTTg                                                             >  1:94806/5‑250 (MQ=255)                                                                                                                    cttATGAATATGTTAATACCTGTTTGCTTTTTAAGCGCGT‑GCAGCTCGTCACCACCAGAACCGCCGCAAATTGATAAAGACGCCCCAGTAATCCATTTAAATAATCAGGTGTACCGACAAACACCAACTGCAACAGTGCCAAAAAACGCCGGCGATCATAATGGGCAACCTTGGGTATATCaa                                                                                                                             >  1:224595/2‑183 (MQ=255)                                                                                                                    cttATGAATATGTTAATACCTGTTTGCTTTTTAAGCGCGT‑GCAGCTCGTCACCACCAGAACCGCCGCAAATTGATAAAGACGCCCCAGTAATCCATTTAAATAATCAGGTGTACCGACAAACACCAACTGCAACAGTGCCAAAAAACGCCGGCGATCATAATGGGCAACCTTGGGTATATCAATATGTGAACCTAGGTAGAAGCGATTATGTGAATGAAACAGAAAAAGTTCGTTTTTTCTACTTTGcgca                                                         >  1:892947/2‑248 (MQ=255)                                                                                                                        tGAATATGTTAATACCTGTTTGCTTTTTAAGCGCGT‑GCAGCTCGTCACCACCAGAACCGCCGCAAATTGATAAAGACGCCCCAGTAATCCATTTAAATAATCAGGTGTACCGACAAACACCAACTGCAACAGTGCCAAAAAACGCCGGCGATCATAATGGGCAACCTTGGGTATATCAATATGTGAACCTAGGTAGAAGCGATTATGTGAATGAAACAGAAAAAGTTCGTTTTTTCTACtt                                                               >  1:499998/1‑241 (MQ=255)                                                                                                                            tatgTTAATACCTGTTTGTTTTTTAAGCGCGT‑GCAGCTCGTCACAACCAGAACCGCCGCAAATTGATAAAGACGCCCCAGTAATCCATTTAAATAATCAGGTGTACCGACAAACACCAACTGCAACAGTGCCAAAAAACGCCGGCGATCATAATGGGCAACCTTGGGTATATCAATATGTGAACCTAGGTAGAAGCGATTATGTGAATGAAACAGAAAAAGTTCGTTTTTTCTACTTTGCGCATCATGCAg                                                 <  1:451406/248‑1 (MQ=255)                                                                                                                               gTTAATCCCTTTTTGTTTTTTAAGCGCGT‑GCAGCTCGCCCCCACCAAAACCGCCGCAAATTGATAAAGACGCCCCAGTAATCCATTTAAATAATCAGGTGTACCGACAAACACCAACTGCAACAGTGCCAAAAAACGCCGGCGATCATAATGGGCAACCTTGGGTATATCAATAtg                                                                                                                         <  1:705712/176‑1 (MQ=255)                                                                                                                               gTTAATACCTGTTTGCTTTTTAAGCGCGT‑GCAGCTCGTCACCACCAGAACCGCCGCAAATTGATAAAGACGccccagtaatcc                                                                                                                                                                                                                      >  1:640671/1‑75 (MQ=255)                                                                                                                               gTTAATACCTGTTTGCTTTTTAAGCGCGT‑GCAGCTCGTCACCACCAGAACCGCCGCAAATTGATAAAGACGCCCCAGTAATCCATTTAAATAATCAGGTGTACCGACAAACACCAACTGCAACAGTGCCAAAAAACGCCGGCGATCATaatg                                                                                                                                                 >  1:316859/1‑150 (MQ=255)                                                                                                                               gTTAATACCTGTTTGCTTTTTAAGCGCGT‑GCAGCTCGTCACCACCAGAACCGCCGCAAATTGATAAAGACGCCCCAGTAATCCATTTAAATAATCAGGTGTACCGACAAACACCAACTGCAACAGTGCCAAAAAACGCCGGCGATCATAATGGGCAACCTTGGGt                                                                                                                                    >  1:603154/1‑165 (MQ=255)                                                                                                                               gTTAATACCTGTTTGCTTTTTAAGCGCGT‑GCAGCTCGTCACCACCAGAACCGCCGCAAATTGATAAAGACGCCCCAGTAATCCATTTAAATAATCAGGTGTACCGACAAACACCAACTGCAACAGTGCCAAAAAACGCCGGCGATCATAATGGGCAACCTTGGGt                                                                                                                                    >  1:670508/1‑165 (MQ=255)                                                                                                                               gTTAATACCTGTTTGCTTTTTAAGCGCGT‑GCAGCTCGTCACCACCAGAACCGCCGCAAATTGATAAAGACGCCCCAGTAATCCATTTAAATAATCAGGTGTACCGACAAACACCAACTGCAACAGTGCCAAAAAACGCCGGCGATCATAATGGGCAACCTTGGGTATATCAATAtg                                                                                                                         >  1:588047/1‑176 (MQ=255)                                                                                                                               gTTAATACCTGTTTGCTTTTTAAGCGCGT‑GCAGCTCGTCACCACCAGAACCGCCGCAAATTGATAAAGACGCCCCAGTAATCCATTTAAATAATCAGGTGTACCGACAAACACCAACTGCAACAGTGCCAAAAAACGCCGGCGATCATAATGGGCAACCTTGGGTATATCAATATGTGAACCTAGGTAGAAGCGATTATGTGAATGAAACAGAAAAAGTTCGTTTTTTCTACTTTGCGCATCATGCAGATa                                              >  1:265026/1‑251 (MQ=255)                                                                                                                               gTTAATACCTGTTTGCTTTTTAAGCGCGT‑GCAGCTCGTCACCACCAGAACCGCCGCAAATTGATAAAGACGCCCCAGTAATCCATTTAAATAATCAGGTGTACCGACAAACACCAACTGCAACAGTGCCAAAAAACGCCGGCGATCATAATGGGCAACCTTGGGTATATCAATATGTGAACCTAGGTAGAAGCGATTATGTGAATGAAACAGAAAAAGTTCGTTTTTTCTACTTTGCGCATCATGCAGATa                                              >  1:492869/1‑251 (MQ=255)                                                                                                                               aggtgtataagagacattttaAAGCGCGT‑GCAGCTCGTCCCCACCAGAACCGCCGCAAATTGATAAAGACGCCCCAGTAATCCATTTAAATAATCAGGTGTACCGACAAACACCAACTGCAACAGTGCCAAAAAACGCCGGCGATCATAATGGGCAACCTTGGGTATATCAATATGTGAACCTAGGTAGAAGCGATTATGTGAATGAAACAGAAAAAGTTCGTTTTTTCTACTTTGCGCATCATGCAGAt                                               <  1:417446/234‑1 (MQ=255)                                                                                                                                 tAATACCTGTTTGCTTTTTAAGCGCGT‑GCAGCTCGTCACCACCAGAACCGCCGCAAATTGATAAAGACGCCCCAGTAATCCATTTAAATAATCAGGTGTACCGACAAACACCAACTGCAACAGTGCCAAAAAACGCCGGCGATCATAATGGGCAACCTTGGGTATATCAATATGTGAACCTAGGTAGAAGCGATTAtg                                                                                                 >  1:86519/1‑198 (MQ=255)                                                                                                                                             gCTTTTTAAGCGCGT‑GCAGCTCGTCACCACCAGAACCGCCGCAAATTGATAAAGACGCCCCAGTAATCCATTTAAATAATCAGGTGTACCGACAAACACCAACTGCAACAGTGCCAAAAAACGCCGGCGATCATAATGGGCAACCTTGGGt                                                                                                                                    >  1:834182/1‑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)                                                                                                                                                                                                |                                                                                                                                                                                                                                       ATTACCGACTAAAGCGCAACTGGAAAAAGCCATGTTCAGTGGTGAGCTTTCCGCCCAATCTATTTAATCCACAAAAACTGTGAATAACCAGGAGGATAATATGAAAAAATTTTTAGTTATGAATACGTTAATCCCTGTTTGCTTTTTAAGCGCGT‑GCAGCTCGTCAACGCCGGAACCGCCACAAATCGATGAAGACGCCCCGGCAATTTATTTAAATAATCAGGTGTACCGACAAACACCAACTGCAACAGTGCCAAAAAACGCCGGCGATCATAACGGGCAACCTTGGGTATATCAATATGTGAACTTAGATAGAAACGATTATGTGAATGAAACAGAAAAAGTTCGTTTTTTCTACTTTGCTCATCATGCAGATAGTATTGAAATCTACGGTCAGCCAGCGCGTACCGAAGCCTATAAA  >  1/778047‑778467

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