Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	1	885,622	G→A	V8I (GTT→ATT)	unknown →	hypothetical protein

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	1	885,622	0	G	A	100.0%	79.9 / NA	35	V8I (GTT→ATT)	unknown	hypothetical protein
Reads supporting (aligned to +/- strand):  ref base G (0/0);  new base A (18/17);  total (18/17)

CTTTGTCGATCCGTATCTGAGCGAGAAGTAAGAAAAATTCGGAATGTGTAAAACTATTTCGGAAAATCAGGGTCAAACGGGCATTCTTGCTTCAGTCGCCATTTATATTAAAGCCACATATAGTTACACGTAATTCAGACTTTTAAGGTTGCTGTAGCGAGGCAATGATTCATCTATTTTCAACCTGAAGGAATATCATGTTAAAAAAATACTTACTTGTTTTACCGATTTCTATTTTAACCGCCTGCGCGGGTTATTATCCTTATCCGCAATTTGAACCCAGTTATCC‑GTCCTTTTCCGTGCAGCATAATGACGATGAAGAACATGCCGTTTTCAGCATCCGCCGCACGCATGAACCGAAATATCAGGATCCAAGATACAAAATTTCCACGGAAAACATGGAAAAATTCATTTTAGCCAAAAATAACCTGGAATATTGTTTGTTACCGGAATTGGGCAAGCAA  >  1/885404‑885867                                                                                                                                                                                                                           |                                                                                                                                                                                                                                                       cttTGTCGATCCGTATCTGAGCGAGAAGTAAGAAAAAGTCGGAAAGTGTAAAACTATTTCGAAAAAACAGGGTCAAACGGGCTTTCTTGCTTCAGTTGCCAT‑‑‑TA‑‑AAAGCCATATATGGTTATACGTAATTCAGGCTTTTAAGGTTGCTGCAGCGGGACAATAATTAATCTATTTCCACCCTGAAGAAATATCATGTAAAAAAAATCCTTACTTATTTTACCGATTTCTATTTTAACCGCCTGTGCCGGtt                                                                                                                                                                                                                    <  1:189979/250‑1 (MQ=255)       cGATCCGTATCTGAGCGAGAAGTAAAAAAAAGTCGAAAAGTGTAAAACTATTTCGAAAAAACAGGGTCAACCGGGCATTCTTGCTTCAGTTGCCAT‑‑‑TA‑‑AAAGCCATATATGTTTATACGTAATTCAGGCTTTTAAGGTTGCTGCAGCGGGACAATGATTAATCTTTTTTCAACCTGAAGGAATATCATGTAAAAAAAATACTTATTTATTTTACCGATTTCTATTTTAACCGCCTGTGCCGGTTAttatcc                                                                                                                                                                                                             <  1:386685/251‑1 (MQ=255)        gATCCGTATCTGAGCGAGAAGTAAGAAAAAGGCGGAAAGTGTAAAACTATTTCGAAAAAACAGGGTCAAACGGGCTTTCTTGCTTCAGTTGCCAT‑‑‑TA‑‑AAAGCCATATATGGTTATACGTAATTCAGGCTTTTAAGGTTGCTGCAGCGGGACAATGATTAATCTATTTTCAACCTAAAGGAATATCATGTTAAAAAAATACTTACTTATTTTACCGATTTCTATTTTAACCGCCTGTGCCGGTTAttatcc                                                                                                                                                                                                             <  1:43973/250‑1 (MQ=255)            atccgtatctgagcgagaagtaagaaaaagtcggaaagtgtaaaactatttcggaaaaacagggtcaaacgggcattcttgcttcagttgccattA‑‑AAAGCCATATATGGTTATACGTAATTCAGGCTTTTAAGGTTGCTGCAGCGGGACAATGATTAATCTATTTTCACCCTGAAGGAATATCATGTTAAAAAAATACTTACTTATTTTACCGATTTCTATTTTAACCGCCTGTGCCGGTTAttatc                                                                                                                                                                                                              >  1:444596‑M2/96‑248 (MQ=255)            atccgtatctgagcgagaagtaagaaaaagtcggaaagtgtaaaactatttcggaaaaacagggtcaaacgggcattcttgcttcagttgccattA‑‑AAAGCCATATATGGTTATACGTAATTCAGGCTTTTAAGGTTGCTGCAGCGGGACAATGATTAATCTATTTTCAACCTGAAGGAATATCATGTTAAAAAAATACTTACTTATTTTCCCGATTTCTATTTTAACCGCCTGTGCCGGtt                                                                                                                                                                                                                    <  1:580545‑M2/147‑1 (MQ=255)             gTATCTGAGCGAGAAGTAAGAAAAAGTCGGAGAGTGTAAAACTATTTCGGAAAAACAGGGTCAAACGGGCATTCTTGCTTCAGTTGCCAT‑‑‑TA‑‑AAAGCCATATATGGTTATACGTAATTCAGGCTTTTAAGGTTGCTGCAGCGGGACAATGATTAATCTATTTTCAACCTGAAGGAATATCATGTTAAAAAAATACTTACTTATTTTACCg                                                                                                                                                                                                                                                >  1:42256/1‑210 (MQ=255)                 tatctgagcgagaagtaagaaaaagtcggaaagtgtaaaactatttcggaaaaacagggtcaaacgggcattcttgcttcagttgccattA‑‑AAAGCCATATATGGTTATACGTAATTCAGGCTTTTAAGGTTGCTGCAGCGGGACAATGATTAATCTATTTTCAACCTGAAGGAATATCATGTTAAAAAAATACTTACTTATTTTACCGATTTCTATTTTAACCGCCTGTGCCGGTTATTATCCTTATCCg                                                                                                                                                                                                      >  1:711569‑M2/91‑251 (MQ=255)                          ggagatgtgtataagaaaaattgtAAAACTATTTCGGAAAAACAGGGTCAAACGGGCATTCTTGCTTCAGTTGCCA‑‑TA‑A‑‑AAAGCCATATATGGTTATACGTAATTCAGGCTTTTAAGGTTGCTGCAGCGGAAAAATGATTAATCTTTTTTCAACCTGAAGGAATATCATTTTAAAAAAATACTTACTTATTTTACCGATTTCTATTTTAACCGCCTGTGCCGGTTATTATCCTTATCCGCAATTTGAAcc                                                                                                                                                                                           <  1:522820/229‑1 (MQ=255)                                  aaaaGTCGGAAAGTGAAAAACTATTTCGGAAAACCAGGGTAAAACGGGCATTCTTGCTTCAGTTGCCAT‑‑‑TA‑‑AAAGCCATATATGTTTATACGTAATTCAGGCTTTTAAGTTTGCTGCAGCGGGAAAATGATTAATCTATTTTCAACCTGAAGGAATATCATGTTAAAAAAATACTTTCTTATTTTACCGATTTCTATTTTAACCGCCTGTGCCGGTTATTATCCTTATCCGCAATTTGAACCTGGTTATcc                                                                                                                                                                                  <  1:223172/251‑1 (MQ=255)                                                 tAAAACTATATCGAAAAAACAGGGCCAACGGGGCATTCTTGCTTCAGTTGCCAT‑‑‑TA‑‑AAAGCCATATATGTTTATACGTAATCCAGGTTTTTAGGGTGGCGCCAGCGGAACAATGATTATTCTATTTTCACCCTAAAGGAATATCATTTTAAAAAAATACTTGTTTATTTTACCGTTTTGTATTTAAACCGCCTGTGCCGGTTATTATCCTTATCCGCAATTTGAACCTGGTTATCC‑ggtcccttccgtgc                                                                                                                                                                   <  1:508630/250‑14 (MQ=255)                                                 gtgtaaaactatttcggaaaaacagggtcaaacgggcattcttgcttcagttgccattA‑‑AAAGCCATATATGGTTATACGTAATTCAGGCTTTTAAGGTTGCTGCAGCGGGACAATGATTAATCTATTTTCAACCTGAAGGAATATCATGTTAAAAAAATACTTACTTATTTTACCGATTTCTATTTTAAc                                                                                                                                                                                                                                  <  1:771184‑M2/133‑1 (MQ=255)                                                 gtgtaaaactatttcggaaaaacagggtcaaacgggcattcttgcttcagttgccattA‑‑AAAGCCATATATGGTTATACGTAATTCAGGCTTTTAAGGTTGCTGCAGCGGGACAATGATTAATCTATTTTCAACCTGAAGGAATATCATGTTAAAAAAATACTTACTTATTTTACCGATTTCTATTTTAACCGCCTGTGCCGGTTATTATCCTTATCCGCAATTTGAACCTGGTTatcc‑ggtcc                                                                                                                                                                            <  1:396556‑M2/186‑5 (MQ=255)                                                          ttCGGAAAAACAGGTTAAAACGGGAATTCTTGCTTCATTTGCCAT‑‑‑TA‑‑AAAGCCATATATGTTTATACGTAATTCAGGTTTTTAAGGTTGCTGCAGCGGAACAATGATTAATCTATTTTCAACCTAAAGGAATATCATGTAAAAAAAATACTTACTTATTTTACCGATTTCTATTTTAACCGCCTGTGCCGGTTATTATCCTTATCCGCAATTTGAACCTGGTTATCCGGTCC‑CTTCCGTGCAGCATAATGa                                                                                                                                                         <  1:810714/251‑1 (MQ=255)                                                          ttCGGAAAAACAGGGTCAAACGGGCATTCTTGCTTCAGTTGCCAT‑‑‑TA‑‑AAAGCCATATATGGTTATACGTAATTCAGGTTTTTAAGGTTGCTGCAGCGGGACAATGATTAATCTATTTTCACCCTGAAGGAATATCATGTTAAAAAAATACTTACTTATTTTCCCGATTTCTATTTTAACCGCCTGTGCCGGTTATTATCCTTATCCGCAATTTGAACCTGGTTATCCGGTCC‑CTTCCGTGCAGCATAATGa                                                                                                                                                         <  1:710923/251‑1 (MQ=255)                                                              gaaaaacggggTCAAACGGGCATTCTTGCTTCAGTTGCCAT‑‑‑TA‑‑AAAGCCATATATGGTTATACGTAATTCAGGCTTTTAAGGTTGCTGCAGCGGAACAATGATTAATCTATTTTCACCCTGAAGGAATATCATGTTAAAAAAATACTTACTTATTTTACCGATTTCTTTTTTAACCGCCTGTGCCGGTTATTATCCTTATCCGCAATTTGAACCTGGTTATCCGGTCC‑CTTCCGTGCAGCATAATGa                                                                                                                                                         <  1:143318/239‑1 (MQ=255)                                                                      gggTCAAACGGGCATTCTTGCTTCAGTTGCCAT‑‑‑TA‑‑AAAGCCATATATGGTTATACGTAATTCAGGCTTTTAAGGTTGCTGCAGCGGGACAATGATTAATCTATTTTCACCCTGAAGGAATATCATGTTAAAAAAATACTTACTTATTTTACCGATTTCTATTTTAACCGCCTGTGCCGGTTAttatcc                                                                                                                                                                                                             <  1:330856/188‑1 (MQ=255)                                                                        gTCAAACGGGCATTCTTGCTTCAGTTGCCAT‑‑‑TA‑‑AAAGCCATATATGGTTATACGTAATTCAGGCTTTTAAGGTTGCTGCAGCGGGACAATGATTAATCTATTTTCAACCTGAAGGAATATCATGTTAAAAAAATACTTACTTATTTTACCGATTTCTATTTTAACCGCCTGTGCCGGTTATTATCCTTATCCGCAATTTGAACCTGGTTATCCGGTCC‑CTTCCGTGCAGCATAATGATGATGAAGAAAATg                                                                                                                                           >  1:327642/1‑251 (MQ=255)                                                                               gggCATTCTTGCTTCAGTTGCCAT‑‑‑TA‑‑AAAGCCATATATGGTTATACGTAATTCAGGCTTTTAAGGTTGCTGCAGCGGGACAATGATTAATCTATTTTCAACCTGAAGGAATATCATGTTAAAAAAATACTTACTTATTTTACCGATTTCTATTTTAACCGCCTGTGCCGGTTATTATCCTTATCCGCAATTTGAACCTGGTTATCCGGTCC‑CTTCCGTGCAGCATAATGATGATGAAGAAAATGCCGtttt                                                                                                                                    >  1:153639/1‑251 (MQ=255)                                                                                          tcttgcttcattttccattaaaaGCCATATATGGTTATACGTAATTCAGGCTTTTAAGGTTGCTGCAGCGGAACAATGATTAATCTATTTTCAACCTGAAGGAATATCATGTTAAAAAAATACTTACTTATTTTACCGATTGTTATTTTAGCCGCCTGTGCCGGTTATTATCCTTATCCGCAATTTGAACCTGGTTATCCGGTCC‑CTTCCGTGCAGCATaa                                                                                                                                                            <  1:258853/201‑1 (MQ=255)                                                                                                                 gCCATATATGGTTATACGTAATTCAGGCTTTTAAGGTTGCTGCAGCGGGACAATGATTAATCTATTTTCAACCTGAAGGAATATCATGTTAAAAAAATACTTACTTATTTTACCGATTTCTATTTTAACCGCCTGTGCCGGTTATTATCCTTATCCGCAATTTGAACCTGGTTATCCGGTCC‑CTTCCGTGCAGCATAATGa                                                                                                                                                         >  1:593899/1‑201 (MQ=255)                                                                                                                    atatatGGTTATACGTAATTCAGGCTTTTAAGGTTGCTGCAGCGGGAAAATGATTAATCTTTTTTCAACCTGAAGGAATATCATGTTAAAAAAATACTTACTTATTTTCCCGATTTCTATTTTAACCGCCTGTGCCGGTTATTATCCTTATCCGCAATTTGAACCTGGTTATCCGGTCC‑CTTCCGTGCAGCATAATGatg                                                                                                                                                       <  1:762781/198‑3 (MQ=255)                                                                                                                     tatatGGTTATACGTAATTCAGGCTTTTAAGGTTGCTGCAGCGGGACAATGATTAATCTATTTTCAACCTGAAGGAATATCATGTTAAAAAAATACTTACTTATTTTACCGATTTCTATTTTAACCGCCTGTGCCGGTTATTATCCTTATCCGCAATTTGAACCTGGTTATCCGGTCC‑CTTCCGTGCAGCATAATGATGATGAAGAAAATGCCGTTTTCAGCACCCGCCGCACGCATGAACCGAAATATCa                                                                                                   >  1:135765/2‑251 (MQ=255)                                                                                                                      atatGGTTATACGTAATTCAGGCTTTTAAGGTTGCTGCAGCGGGACAATGATTAATCTATTTTCAACCTGAAGGAATATCATGTTAAAAAAATACTTACTTATTTTACCGATTTCTATTTTAACCGCCTGTGCCGGTTATTATCCTTATCCGCAATTTGAACCTGGTTATCCGGTCC‑CTTCCGTGCAGCATAATGa                                                                                                                                                         >  1:819786/1‑196 (MQ=255)                                                                                                                      atatGGTTATACGTAATTCAGGCTTTTAAGGTTGCTGCAGCGGGACAATGATTAATCTATTTTCAACCTGAAGGAATATCATGTTAAAAAAATACTTACTTATTTTACCGATTTCTATTTTAACCGCCTGTGCCGGTTATTATCCTTATCCGCAATTTGAACCTGGTTATCCGGTCC‑CTTCCGTGCAGCATAATGATGATGAAGAAAATGCCGTTTTCAGCACCCGCCGCACGCATGAACCGAAATATCAg                                                                                                  >  1:693307/1‑251 (MQ=255)                                                                                                                      atatGGTTATACGTAATTCAGGCTTTTAAGGTTGCTGCAGCGGGACAATGATTAATCTATTTTCAACCTGAAGGAATATCATGTTAAAAAAATACTTACTTATTTTACCGATTTCTATTTTAACCGCCTGTGCCGGTTATTATCCTTATCCGCAATTTGAACCTGGTTATCCGGTCC‑CTTCCGTGCAGCATAATGATGATGAAGAAAATGCCGTTTTCAGCACCCGCCGCACGCATGAACCGAAATATCAg                                                                                                  >  1:231242/1‑251 (MQ=255)                                                                                                                           gTTATACGTAATTCAGGCTTTTAAGGTTGCTGCAGCGGGACAATGATTAATCTATTTTCAACCTGAAGGAATATCATGTTAAAAAAATACTTACTTATTTTACCGATTTCTATTTTAACCGCCTGTGCCGGTTATTATCCTTATCCGCAATTTGAACCTGGTTATCCGGTCC‑CTTCCGTGCAGCATAATGATGATGAAGAAAATGCCGTTTTCAGCACCCGCCGCACGCATGAACCGAAATATCAGGATc                                                                                              >  1:357354/1‑250 (MQ=255)                                                                                                                                acttAATTCAGGTTTTTAAGGTTGCTGCAGCGGGACAATGATTAATCTTTTTTCAACCTAAAGAAATATCATGTTAAAAAAATACTTACTTATTTTACCGTTTTCTTTTTTAACCGCCTGTGCCGGTTATTATCCTTATCCGCAATTTGAACCTGGTTATCCGGTCC‑CTTCCGTGCAGCATAATGa                                                                                                                                                         <  1:840948/183‑1 (MQ=255)                                                                                                                                acGTAATTCAGGCTTTTAAGGTTGCTGCAGCGGGACAATGATTAATCTATTTTCAACCTGAAGGAATATCATGTTAAAAAAATACTTACTTATTTTACCGATTTCTATTTTAACCGCCTGTGCCGGTTATTATCCTTATCCGCAATTTGAACCTGGTTATCCGGTCC‑CTTCCGTGCAGCATAATGATGATGAAGAAAATGCCGTTTTCAGCACCCGCCGCACGCATGAACCGAAATATCAGGATCCGAGAt                                                                                        >  1:215132/1‑251 (MQ=255)                                                                                                                                                        ctgcagcgggacAATGATTAATCTTTTTTCATCCTGAGGAAATATCATTTTAAAAAAATACTTTCTTATTTTCCCGATTTCTATTTTATCCGCCTGTGCCGGTTATTATCCTTATCCGCAATTTGAACCTGGTTATCCGGTCC‑CTTCCGTGCAGCATAAt                                                                                                                                                           <  1:253216/149‑1 (MQ=255)                                                                                                                                                             gcgggacAATGATTAATCTATTTTCAACCTGAAGGAATATCATGTTAAAAAAATACTTACTTATTTTACCGATTTCTATTTTAACCGCCTGTGCCGGTTATTATCCTTATCCGCAATTTGAACCTGGTTATCCGGTCC‑CTTCCGTGCAGCATAATGATGATGAAGAAAATGCCGTTTTCAGCACCCGCCGCACGCATGAACCGAAATATCAGGATCCGAGATACAAAATTTCCACAGAAGACATGGaaaaa                                                           >  1:772262/7‑251 (MQ=255)                                                                                                                                                                ggacAATGATTAATCTATTTTCAACCTGAAGGAATATCATGTTAAAAAAATACTTACTTATTTTACCGATTTCTATTTTAACCGCCTGTGCCGGTTATTATCCTTATCCGCAATTTGAACCTGGTTATCCGGTCC‑CTTCCGTGCAGCATAATGATGATGAAGAAAATGCCGTTTTCAGCACCCGCCGCACGCATGAACCGAAATATCAGGATCCGAGATACAAAATTTCCACAGAAGACATGGAAAAAtt                                                         >  1:23965/4‑250 (MQ=255)                                                                                                                                                                                aTTTTCAACCTGAAGGAATATCATGTTAAAAAAATACTTACTTATTTTACCGATTTCTATTTTAACCGCCTGTGCCGGTTATTATCCTTATCCGCAATTTGAACCTGGTTATCCGGTCC‑CTTCCGTGCAGCAt                                                                                                                                                              >  1:409835/1‑133 (MQ=255)                                                                                                                                                                                              ggAATATCATGTTAAAAAAATACTTACTTATTTTACCGATTTCTATTTTAACCGCCTGTGCCGGTTATTATCCTTATCCGCAATTTGAACCTGGTTATCCGGTCC‑CTTCCGTGCAGCATAATGATGATGAAGAAAATGCCGTTTTCAGCACCCGCCGCACGCATGAACCGAAATATCAGGATCCGAGATACAAAATTTCCACAGAAGACATGGAAAAATTCATTTTAGCCAAAAATAACTTGGAATAttgt                          >  1:417865/1‑251 (MQ=255)                                                                                                                                                                                                                   acttacttATTTTACCGATTTCTATTTTAACCGCCTGTGCCGGTTATTATCCTTATCCGCAATTTGAACCTGGTTATCCGGTCC‑CTTCCGTGCAGCATAATGATGATGAAGAAAATGCCGTTTTc                                                                                                                                   >  1:673166/1‑125 (MQ=255)                                                                                                                                                                                                                    cttacttaTTTTACCGATTTCTATTTTAACCGCCTGTGCCGGTTATTATCCTTATCCGCAATTTGAACCTGGTTATCCGGTCC‑CTTCCGTGCAGCATAATGATGATGAAGAAAATGCCGTTTTCAGCACCCGCCGCACGCATGAACCGAAATATCAGGATCCGAGATACAAAATTTCCACAGAAGACATGGAAAAATTCATTTTAGCCAAAAATAACTTGGAATATTGTTTGTTGCCGGAATTGGgcaag     >  1:586646/1‑250 (MQ=255)                                                                                                                                                                                                                      tacttATTTTACCGATTTCTATTTTAACCGCCTGTGCCGGTTATTATCCTTATCCGCAATTTGAACCTGGTTATCCGGTCC‑CTTCCGTGCAGCATAATGATGATGAAGAAAATGCCGTTTTCAGCACCCGCCGCACGCATGAACCGAAATATCAGGATCCGAGATACAAAATTTCCACAGAAGACATGGAAAAATTCATTTTAGCCAAAAATAACTTGGAATATTGTTTGTTGCCGGAATTGGgcaagcaa  >  1:664378/1‑251 (MQ=255)                                                                                                                                                                                                                           |                                                                                                                                                                                                                                                       CTTTGTCGATCCGTATCTGAGCGAGAAGTAAGAAAAATTCGGAATGTGTAAAACTATTTCGGAAAATCAGGGTCAAACGGGCATTCTTGCTTCAGTCGCCATTTATATTAAAGCCACATATAGTTACACGTAATTCAGACTTTTAAGGTTGCTGTAGCGAGGCAATGATTCATCTATTTTCAACCTGAAGGAATATCATGTTAAAAAAATACTTACTTGTTTTACCGATTTCTATTTTAACCGCCTGCGCGGGTTATTATCCTTATCCGCAATTTGAACCCAGTTATCC‑GTCCTTTTCCGTGCAGCATAATGACGATGAAGAACATGCCGTTTTCAGCATCCGCCGCACGCATGAACCGAAATATCAGGATCCAAGATACAAAATTTCCACGGAAAACATGGAAAAATTCATTTTAGCCAAAAATAACCTGGAATATTGTTTGTTACCGGAATTGGGCAAGCAA  >  1/885404‑885867

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