Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	1	302,175	C→T	P181S (CCT→TCT)	unknown →	Lysine decarboxylase family

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	1	302,175	0	C	T	100.0%	83.1 / NA	32	P181S (CCT→TCT)	unknown	Lysine decarboxylase family
Reads supporting (aligned to +/- strand):  ref base C (0/0);  new base T (18/14);  total (18/14)

TGTATTGCTTTGCCGGGTGGGCCGGGCACGCTGGAAGAAATCAGCGAAGTCTATTCTTGGGCGCGAATCGGTAAAAATCCTCATCCTTGTATTCTATTTAACCAAGATGGTTTTTATGAACCGCTGAGAAGTCAATATCAAAGCATGGTCACTGCCGGTTTCTTGACGCAAAGTCATTTTGATAAGTTATTATTTAGTACAAATATTTTTGAAATAGAACGTTTTATTCTTGATTATCAAGTGCCTGAGATTCGTACTTATGAAGATGTGTAAAAAGTGCGGTCAAAAATAAAGTTGTTTTTCTATGAGAACTAACTCGGTGGATTTTTCGTACATTACATGATTTTTCTTATTACCCCAATCTAAGTTGTTACTTTTTTCGGCATCTTATTTTTTATTTTAAAATAAGTTTGCTTCTATGGATTTTTTTGTGTAAGATGGCGCGCATTATTATTTTTCACTAAAAAATAAGG  >  1/301932‑302404                                                                                                                                                                                                                                                    |                                                                                                                                                                                                                                      tGTATTGCTTTGCCGGGTGGGCCGGGCACGCTGGAAGAAATCAGCGAAGTCTATTCTTGGGCGCGAATCGGTAAAAATCCTCATCCTTGTATTCTATTTAACCAAGATGGTTTTTATGAACCGCTGAGAAGTCAATATCAAAGCATGGTCACTGCCGGTTTCTTGACGCAAAGTCATTTTGATAAGTTATTATTTAGTACAAATATTTTTGAAATAGAACGTTTTATTCTTGATTATCAAGTGTCTGAGAt                                                                                                                                                                                                                                <  1:698375/251‑1 (MQ=255)             ccGGGTGGGCCGGGCACGCTGGAAGAAATCAGCGAAGTCTATTCTTGGGCGCGAATCGGTAAAAATCCTCATCCTTGTATTCTATTTAACCAAGATGGTTTTTATGAACCGCTGAGAAGTCAATATCAAAGCATGGTCACTGCCGGTTTCTTGACGCAAAGTCATTTTGATAAGTTATTATTTAGTACAAATATTTTTGAAATAGAACGTTTTATTCTTGATTATCAAGTGTCTGAGATTCGTACTTATg                                                                                                                                                                                                                     <  1:237947/250‑1 (MQ=255)                             cGCTGGAAGAAATCAGCGAAGTCTATTCTTGGGCGCGAATCGGTAAAAATCCTCATCCTTGTATTCTATTTAACCAAGATGGTTTTTATGAACCGCTGAGAAGTCAATATCAAAGCATGGTCACTGCCGGTTTCTTGACGCAAAGTCATTTTGATAAGTTATTATTTAGTACAAATATTTTTGAAATAGAACGTTTTATTCTTGATTATCAAGTGTCTGAGATTCGTACTTATGAAGATGTGTAAAAAGTg                                                                                                                                                                                                    <  1:383232/251‑1 (MQ=255)                                 ggAAGAAATCAGCGAAGTCTATTCTTGGGCGCGAATCGGTAAAAATCCTCATCCTTGTATTCTATTTAACCAAGATGGTTTTTATGAACCGCTGAGAAGTCAATATCAAAGCATGGTCACTGCCGGTTTCTTGACGCAAAGTCATTTTGATAAGTTATTATTTAGTACAAATATTTTTGAAATAGAACGTTTTATTCTTGATTATCAAGTGTCTGAGATTCGTACTTATGAAGATGTGTAAAAAGTGCGGt                                                                                                                                                                                                >  1:816968/1‑251 (MQ=255)                                                   cTATTCTTGGGCGCGAATCGGTAAAAATCCTCATCCTTGTATTCTATTTAACCAAGATGGTTTTTATGAACCGCTGAGAAGTCAATATCAAAGCATGGTCACTGCCGGTTTCTTGACGCAAAGTCATTTTGATAAGTTATTATTTAGTACAAATATTTTTGAAATAGAACGTTTTATTCTTGATTATCAAGTGTCTGAGATTCGTACTTATGAAGATGTGTAAAAAg                                                                                                                                                                                                      >  1:698775/1‑227 (MQ=255)                                                            ggCGCGAATCGGTAAAAATCCTCATCCTTGTATTCTATTTAACCAAGATGGTTTTTATGAACCGCTGAGAAGTCAATATCAAAGCATGGTCACTGCCGGTTTCTTGACGCAAAGTCATTTTGATAAGTTATTATTTAGTACAAATATTTTTGAAATAGAACGTTTTATTCTTGATTATCAAGTGTCTGAGATTCGTACTTATGAAGATGTGTAAAAAGTGCGGTCAAAAATAAAGTTGTTTTTCTATgag                                                                                                                                                                      <  1:801480/250‑1 (MQ=255)                                                                                    tCCTTGTATTCTATTTAACCAAGATGGTTTTTATGAACCGCTGAGAAGTCAATATCAAAGCATGGTCACTGCCGGTTTCTTGACGCAAAGTCATTTTGATAAGTTATTATTTAGTACAAATATTTTTGAAATAGAACGTTTTATTCTTGATTATCAAGTGTCTGAGAt                                                                                                                                                                                                                                >  1:703526/1‑168 (MQ=255)                                                                                    tCCTTGTATTCTATTTAACCAAGATGGTTTTTATGAACCGCTGAGAAGTCAATATCAAAGCATGGTCACTGCCGGTTTCTTGACGCAAAGTCATTTTGATAAGTTATTATTTAGTACAAATATTTTTGAAATAGAACGTTTTATTCTTGATTATCAAGTGTCTGAGATTCGTACTTATGAAGATGTGTaaaa                                                                                                                                                                                                        >  1:392840/1‑192 (MQ=255)                                                                                    tCCTTGTATTCTATTTAACCAAGATGGTTTTTATGAACCGCTGAGAAGTCAATATCAAAGCATGGTCACTGCCGGTTTCTTGACGCAAAGTCATTTTGATAAGTTATTATTTAGTACAAATATTTTTGAAATAGAACGTTTTATTCTTGATTATCAAGTGTCTGAGATTCGTACTTATGAAGATGTGTAAAAAGTg                                                                                                                                                                                                    <  1:898640/196‑1 (MQ=255)                                                                                              cTATTTAACCAAGATGGTTTTTATGAACCGCTGAGAAGTCAATATCAAAGCATGGTCACTGCCGGTTTCTTGACGCAAAGTCATTTTGATAAGTTATTATTTAGTACAAATATTTTTGAAATAGAACGTTTTATTCTTGATTATCAAGTGTCTGAGATTCGTACTTATGAAGATGTGTAAAAAGTGCGGTCAAAAATAAAGTTGTTTTTCTATGAGAACTAACTCGGTGGATTTTTCGTACATTACATGAt                                                                                                                                   <  1:816601/251‑1 (MQ=255)                                                                                                aTTTAACCAAGATGGTTTTTATGAACCGCTGAGAAGTCAATATCAAAGCATGGTCACTGCCGGTTTCTTGACGCAAAGTCATTTTGATAAGTTATTATTTAGTACAAATATTTTTGAAATAGAACGTTTTATTCTTGATTATCAAGTGTCTGAGATTCGTACTTATGAAGATGTGTAAAAAGTGCGGTCAAAAATAAAGTTGTTTTTCTATGAGAACTAACTCGGTGGATTTTTCGTACATTACATGAttt                                                                                                                                 <  1:369495/251‑1 (MQ=255)                                                                                                         aGATGGTTTTTATGAACCGCTGAGAAGTCAATATCAAAGCATGGTCACTGCCGGTTTCTTGACGCAAAGTCATTTTGATAAGTTATTATTTAGTACAAATATTTTTGAAATAGAACGTTTTATTCTTGATTATCAAGTGTCTGAGATTCGTACTTATGAAGATGTGTAAAAAGTGCGGTCAAAAATAAAGTTGTTTTTCTATGAGAACTAACTCGGTGGATTTTTCGTACATTACATGATTTTTCttatta                                                                                                                        >  1:785063/1‑251 (MQ=255)                                                                                                                         ccGCTGAGAAGTCAATATCAAAGCATGGTCACTGCCGGTTTCTTGACGCAAAGTCATTTTGATAAGTTATTATTTAGTACAAATATTTTTGAAATAGAACGTTTTATTCTTGATTATCAAGTGTCTGAGATTCGTACTTATGAAGATGTGTAAAAAGTGCGGTCAAAAATAAAGTTGTTTTTCTATGAGaac                                                                                                                                                                   <  1:708723/192‑1 (MQ=255)                                                                                                                         ccGCTGAGAAGTCAATATCAAAGCATGGTCACTGCCGGTTTCTTGACGCAAAGTCATTTTGATAAGTTATTATTTAGTACAAATATTTTTGAAATAGAACGTTTTATTCTTGATTATCAAGTGTCTGAGATTCGTACTTATGAAGATGTGTAAAAAGTGCGGTCAAAAATAAAGTTGTTTTTCTATGAGAACTAACTCGGTGGATTTTTCGtacattacat                                                                                                                                      >  1:384778/1‑221 (MQ=255)                                                                                                                         ccGCTGAGAAGTCAATATCAAAGCATGGTCACTGCCGGTTTCTTGACGCAAAGTCATTTTGATAAGTTATTATTTAGTACAAATATTTTTGAAATAGAACGTTTTATTCTTGATTATCAAGTGTCTGAGATTCGTACTTATGAAGATGTGTAAAAAGTGCGGTCAAAAATAAAGTTGTTTTTCTATGAGAACTAACTCGGTGGATTTTTCGTACATTACATGATTTTTCTTATTACCCCAATCTAAgttgt                                                                                                        >  1:336994/1‑251 (MQ=255)                                                                                                                                   gTCAATATCAAAGCATGGTCACTGCCGGTTTCTTGACGCAAAGTCATTTTGATAAGTTATTATTTAGTACAAATATTTTTGAAATAGAACGTTTTATTCTTGATTATCAAGTGTCTGAGATTCGTACTTATGAAGATGTGTAAAAAGTGCGGTCAAAAATAAAGTTGTTTTTCTATGAGAACTAACTCGGTGGATTTTTCGtacattac                                                                                                                                        >  1:434534/1‑209 (MQ=255)                                                                                                                                                          gCCGGTTTCTTGACGCAAAGTCATTTTGATAAGTTATTATTTAGTACAAATATTTTTGAAATAGAACGTTTTATTCTTGATTATCAAGTGTCTGAGATTCGTAc                                                                                                                                                                                                                          >  1:385305/1‑104 (MQ=255)                                                                                                                                                                                          gTTATTATTTAGTACAAATATTTTTGAAATAGAACGTTTTATTCTTGATTATCAAGTGTCTGAGATTCGTACTTATGAAGATGTGTAAAAAGTGCGGTCAAAAATAAAGTTGTTTTTCTATGAGAACTAACTCGGTGGATTTTTCGTACATTACATGATTTTTCTTATTACCCCAATCTAAgt                                                                                                           >  1:789273/1‑183 (MQ=255)                                                                                                                                                                                          gTTATTATTTAGTACAAATATTTTTGAAATAGAACGTTTTATTCTTGATTATCAAGTGTCTGAGATTCGTACTTATGAAGATGTGTAAAAAGTGCGGTCAAAAATAAAGTTGTTTTTCTATGAGAACTAACTCGGTGGATTTTTCGTACATTACATGATTTTTCTTATTACCCCAATCTAAGTTGTTACTTTTTTCGGCATCTTATTTTTTATTTTAAAATAAGTTTGCTTCTATGGATTTTTTTGTGTaa                                       >  1:688774/1‑251 (MQ=255)                                                                                                                                                                                           ttattaTTTAGTACAAATATTTTTGAAATAGAACGTTTTATTCTTGATTATCAAGTGTCTGAGATTCGTACTTATGAAGATGTGTAAAAAGTGCGGTCAAAAATAAAGTTGTTTTTCTATGAGAACTAACTCGGTGGATTTTTCGTACATTACATGATTTTTCTTATTACCCCAATCTAAGTTGTTACTTTTTTCGGCATCttattttt                                                                                >  1:757287/1‑209 (MQ=255)                                                                                                                                                                                             attattTAGTACAAATATTTTTGAAATAGAACGTTTTATTCTTGATTATCAAGTGTCTGAGATTCGTACTTATGAAGATGTGTAAAAAGTGCGGTCAAAAATAAAGTTGTTTTTCTATGAGAACTAACTCGGTGGATTTTTCGTACATTACATGATTTTTCTTATTACCCCAATCTAAGTTGTTACTTTTTTCGGCATCTTATTTTTTATTTTAAAATAAGTTTGCTTCTATGGATTTTTTTGTGTAAGAt                                    <  1:460817/251‑1 (MQ=255)                                                                                                                                                                                               tattTAGTACAAATATTTTTGAAATAGAACGTTTTATTCTTGATTATCAAGTGTCTGAGATTCGTACTTATGAAGATGTGTAAAAAGTGCGGTCAAAAATAAAGTTGTTTTTCTATGAGAACTAACTCGGTGGATTTTTCGTACATTACATGATTTTTCTTATTACCCCAATCTAAGTTGTTACTTTTTTCGGCATCTTATTTTTTATTTTAAAATAAGTTTGCTTCTATGGATTTTTTTGTGTAAGATgg                                  <  1:556311/251‑1 (MQ=255)                                                                                                                                                                                                attTAGTACAAATATTTTTGAAATAGAACGTTTTATTCTTGATTATCAAGTGTCTGAGATTCGTACTTATGAAGATGTGTAAAAAGTGCGGTCAAAAATAAAGTTGTTTTTCTa                                                                                                                                                                          >  1:116392/1‑114 (MQ=255)                                                                                                                                                                                                attTAGTACAAATATTTTTGAAATAGAACGTTTTATTCTTGATTATCAAGTGTCTGAGATTCGTACTTATGAAGATGTGTAAAAAGTGCGGTCAAAAATAAAGTTGTTTTTCTATgag                                                                                                                                                                      <  1:737880/118‑1 (MQ=255)                                                                                                                                                                                                attTAGTACAAATATTTTTGAAATAGAACGTTTTATTCTTGATTATCAAGTGTCTGAGATTCGTACTTATGAAGATGTGTAAAAAGTGCGGTCAAAAATAAAGTTGTTTTTCTATGAGaactaact                                                                                                                                                              <  1:580471/126‑1 (MQ=255)                                                                                                                                                                                                                        aGAACGTTTTATTCTTGATTATCAAGTGTCTGAGATTCGTACTTATGAAGATGTGTAAAAAGTGCGGTCAAAAATAAAGTTGTTTTTCTATGAGAACTAACTCGGTGGATTTTTCGTACATTACATGATTTTTCTTATTACCCCAATCTAAGTTGTTACTTTTTTCGGc                                                                                           <  1:749545/169‑1 (MQ=255)                                                                                                                                                                                                                            cGTTTTATTCTTGATTATCAAGTGTCTGAGATTCGTACTTATGAAGATGTGTAAAAAGTGCGGTCAAAAATAAAGTTGTTTTTCTATGAGAACTAACTCGGTGGATTTTTCGTACATTACATGATTTTTCtt                                                                                                                            >  1:26352/1‑132 (MQ=255)                                                                                                                                                                                                                             gTTTTATTCTTGATTATCAAGTGTCTGAGATTCGTACTTATGAAGATGTGTAAAAAGTGCGGTCAAAAATAAAGTTGTTTTTCTATGAGAACTAACTc                                                                                                                                                             >  1:465674/1‑98 (MQ=255)                                                                                                                                                                                                                                  aTTCTTGATTATCAAGTGTCTGAGATTCGTACTTATGAAGATGTGTAAAAAGTGCGGTCAAAAATAAAGTTGTTTTTCTATGAGAACTAACTCGGTGGATTTTTCGTACATTACATGATTTTTCTTATTACCCCAATCTAAgtt                                                                                                          >  1:463878/1‑144 (MQ=255)                                                                                                                                                                                                                                       tGATTATCAAGTGTCTGAGATTCGTACTTATGAAGATGTGTAAAAAGTGCGGTCAAAAATAAAGTTGTTTTTCTATGAGAACTAACTCGGTGGATTTTTCGTACATTACATGATTTTTCTTATTACCCCAATCTAAGTTGTTACTTTTTTCGGCATCTTATTTTTTATTTTAAAATAAGTTTGCTTCTATGGATTTTTTTGTGTAAGATGGCGCGCattattat                     >  1:713253/1‑224 (MQ=255)                                                                                                                                                                                                                                         aTTATCAAGTGTCTGAGATTCGTACTTATGAAGATGTGTAAAAAGTGCGGTCAAAAATAAAGTTGTTTTTCTATGAGAACTAACTCGGTGGATTTTTCGTACATTACATGATTTTTCTTATTAc                                                                                                                       >  1:621053/1‑124 (MQ=255)                                                                                                                                                                                                                                         aTTATCAAGTGTCTGAGATTCGTACTTATGAAGATGTGTAAAAAGTGCGGTCAAAAATAAAGTTGTTTTTCTATGAGAACTAACTCGGTGGATTTTTCGTACATTACATGATTTTTCTTATTACCCCAATCTAAGTTGTTACTTTTTTCGGCATCTTATTTTTTATTTTAAAATAAGTTTGCTTCTATGGATTTTTTTGTGTAAGATGGCGCGCATTATTATTTTTCACTAAAAAATAAgg  <  1:861486/241‑1 (MQ=255)                                                                                                                                                                                                                                                    |                                                                                                                                                                                                                                      TGTATTGCTTTGCCGGGTGGGCCGGGCACGCTGGAAGAAATCAGCGAAGTCTATTCTTGGGCGCGAATCGGTAAAAATCCTCATCCTTGTATTCTATTTAACCAAGATGGTTTTTATGAACCGCTGAGAAGTCAATATCAAAGCATGGTCACTGCCGGTTTCTTGACGCAAAGTCATTTTGATAAGTTATTATTTAGTACAAATATTTTTGAAATAGAACGTTTTATTCTTGATTATCAAGTGCCTGAGATTCGTACTTATGAAGATGTGTAAAAAGTGCGGTCAAAAATAAAGTTGTTTTTCTATGAGAACTAACTCGGTGGATTTTTCGTACATTACATGATTTTTCTTATTACCCCAATCTAAGTTGTTACTTTTTTCGGCATCTTATTTTTTATTTTAAAATAAGTTTGCTTCTATGGATTTTTTTGTGTAAGATGGCGCGCATTATTATTTTTCACTAAAAAATAAGG  >  1/301932‑302404

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