Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	1	778,215	G→A	T23P (ACG→CCA)	unknown →	lipoprotein

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	1	778,215	0	G	A	100.0%	102.5 / NA	36	T23P (ACG→CCA)	unknown	lipoprotein
Reads supporting (aligned to +/- strand):  ref base G (0/0);  new base A (26/10);  total (26/10)

GTATTACAACGACCCTTATTTCATCGATAAATTGAAAGCGGTTTCCCCACAACTCCGCGCATTAGGCAAGAAATTACCGACTAAAGCGCAACTGGAAAAAGCCATGTTCAGTGGTGAGCTTTCCGCCCAATCTATTTAATCCACAAAAACTGTGAATAACCAGGAGGATAATATGAAAAAATTTTTAGTTATGAATACGTTAATCCCTGTTTGCTTTTTAAGCGCGT‑GCAGCTCGTCAACGCCGGAACCGCCACAAATCGATGAAGACGCCCCGGCAATTTATTTAAATAATCAGGTGTACCGACAAACACCAACTGCAACAGTGCCAAAAAACGCCGGCGATCATAACGGGCAACCTTGGGTATATCAATATGTGAACTTAGATAGAAACGATTATGTGAATGAAACAGAAAAAGTTCGTTTTTTCTACTTTGCTCATCATGCAGATAGTATTGAAATCTACGGT  >  1/777975‑778440                                                                                                                                                                                                                                                  |                                                                                                                                                                                                                                  gtATTACAACGATCCTTATTTTATAGACAAACTCAAAGCGGTATCTCCGCAATTGCGAGCATTAGGTAAAAAATTACCGACAAAAGAACAATTGGAAAAAGCTATGTTCAGTGGCGAGCTGTCAGTCCAATCAATTTAATCCACAAAAACTGTGAATAACCAGGAGGATAATATGAAAAAAATTTTACTTATGAATATGTTAATACCTGTTTGCTTTTTAAGCGCGT‑GCAGCTCGTCACCACCAGAACCg                                                                                                                                                                                                                          >  1:659897/1‑250 (MQ=255)  tATTACAACGATCCTTATTTTATAGACAAACTCAAAGCGGTATCTCCGCAATTGCGAGCATTAGGTAAAAAATTACCGACAAAAGAACAATTGGAAAAAGCTATGTTCAGTGGCGAGCTGTCAGTCCAATCAATTTAATCCACAAAAACTGTGAATAACCAGGAGGATAATATGAAAAAAATTTTACTTATGAATATGTTAATACCTGTTTGCTTTTTAAGCGCGT‑GCAGCTCGTCACCACCAGAACCGc                                                                                                                                                                                                                         <  1:817973/250‑1 (MQ=255)             tccTTATTTTATAGACAAACTCAAAGCGGTATCTCCGCAATTGCGAGCATTAGGTAAAAAATTACCGACAAAAGAACAATTGGAAAAAGCTATGTTCAGTGGCGAGCTGTCAGTCCAATCAATTTAATCCACAAAAACTGTGAATAACCAGGAGGATAATATGAAAAAAATTTTACTTATGAATATGTTAATACCTGTTTGCTTTTTAAGCGCGT‑GCAGCTCGTCACCACCAGAACCGCCGCAAATTGAta                                                                                                                                                                                                             >  1:274096/2‑250 (MQ=255)                                                                         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:492869/1‑251 (MQ=255)                                                                                                                                                                                                       gTTAATACCTGTTTGCTTTTTAAGCGCGT‑GCAGCTCGTCACCACCAGAACCGCCGCAAATTGATAAAGACGCCCCAGTAATCCATTTAAATAATCAGGTGTACCGACAAACACCAACTGCAACAGTGCCAAAAAACGCCGGCGATCATAATGGGCAACCTTGGGTATATCAATATGTGAACCTAGGTAGAAGCGATTATGTGAATGAAACAGAAAAAGTTCGTTTTTTCTACTTTGCGCATCATGCAGATa                   >  1:265026/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)                                                                                                                                                                                                                                                  |                                                                                                                                                                                                                                  GTATTACAACGACCCTTATTTCATCGATAAATTGAAAGCGGTTTCCCCACAACTCCGCGCATTAGGCAAGAAATTACCGACTAAAGCGCAACTGGAAAAAGCCATGTTCAGTGGTGAGCTTTCCGCCCAATCTATTTAATCCACAAAAACTGTGAATAACCAGGAGGATAATATGAAAAAATTTTTAGTTATGAATACGTTAATCCCTGTTTGCTTTTTAAGCGCGT‑GCAGCTCGTCAACGCCGGAACCGCCACAAATCGATGAAGACGCCCCGGCAATTTATTTAAATAATCAGGTGTACCGACAAACACCAACTGCAACAGTGCCAAAAAACGCCGGCGATCATAACGGGCAACCTTGGGTATATCAATATGTGAACTTAGATAGAAACGATTATGTGAATGAAACAGAAAAAGTTCGTTTTTTCTACTTTGCTCATCATGCAGATAGTATTGAAATCTACGGT  >  1/777975‑778440

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