Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	1	91,724	T→C	K229K (AAA→AAG)	unknown ←	hypothetical protein

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	1	91,724	0	T	C	100.0%	83.6 / NA	29	K229K (AAA→AAG)	unknown	hypothetical protein
Reads supporting (aligned to +/- strand):  ref base T (0/0);  new base C (17/12);  total (17/12)

AATGTTATATTATATAATATAACCGCGCAACCGCTATTTTGTAAGTTTTGTGTAAATATTTAAAAGCCCCTTAAACATCGTTTAAAGGGCTTTTGAAGTTAAAATTTATAAACCAGATTATCCTCATACTCGCCATGATAGTTGAGGTTTGTCTTAATCGCTATCCGCTTAGCATTGTCAAACGCTTGCCAGTTATCTACTTTATCCGTTTGTATATACTCAATAAAGCGCACAAACTCGCTTTTGGTTGACGAGAAGAAGATATACGGCGGGCGGGTGATGTTGACTAAGCGCAGGAAGTCGA  >  1/91482‑91785                                                                                                                                                                                                                                                   |                                                              gtctaatattatattatataatataactgcgcaaccggtattttgtaagttttgtgtaaatatttaaaagccccttaaacatcgtttaaagggcttttgaagttaaaatttataaaccagattatcctcatactcgccatgatagttgaggtttgtcttaatcgctatccgcttagcattAACGCTTGCCAGTTATCTACTTTATCCGTTTGTATATACTCAATAAAGCGCACAAACTCGCTCTTAGTTGa                                                       >  1:248287‑M2/181‑251 (MQ=255)      atattatattatataatataactgcgcaaccggtattttgtaagttttgtgtaaatatttaaaagccccttaaacatcgtttaaagggcttttgaagttaaaatttataaaccagattatcctcatactcgccatgatagttgaggtttgtcttaatcgctatccgcttagcattAACGCTTGCCAGTTATCTACTTTATCCGTTTGTATATACTCAATAAAGCGCACAAACTCGCTCTTAGTTGACgaga                                                  >  1:73114‑M2/176‑251 (MQ=255)      atattatattatataatataactgcgcaaccggtattttgtaagttttgtgtaaatatttaaaagccccttaaacatcgtttaaagggcttttgaagttaaaatttataaaccagattatcctcatactcgccatgatagttgaggtttgtcttaatcgctatccgcttagcattAACGCTTGCCAGTTATCTACTTTATCCGTTTGTATATACTCAATAAAGCGCACAAACTCGCTCTTAGTTGACgaga                                                  >  1:504031‑M2/176‑251 (MQ=255)      atattatattatataatataactgcgcaaccggtattttgtaagttttgtgtaaatatttaaaagccccttaaacatcgtttaaagggcttttgaagttaaaatttataaaccagattatcctcatactcgccatgatagttgaggtttgtcttaatcgctatccgcttagcattAACGCTTGCCAGTTATCTACTTTATCCGTTTGTATATACTCAATAAAGCGCACAAACTCGCTCTTAGTTGACgaga                                                  >  1:117604‑M2/176‑251 (MQ=255)                atataatataactgcgcaaccggtattttgtaagttttgtgtaaatatttaaaagccccttaaacatcgtttaaagggcttttgaagttaaaatttataaaccagattatcctcatactcgccatgatagttgaggtttgtcttaatcgctatccgcttagcattAACGCTTGCCAGTTATCTACTTTATCCGTTTGTATATACTCAATAAAGCGCACAAACTCGCTCTTAGTTGacgagaagaatat                                           >  1:278575‑M2/166‑245 (MQ=255)                atataatataactgcgcaaccggtattttgtaagttttgtgtaaatatttaaaagccccttaaacatcgtttaaagggcttttgaagttaaaatttataaaccagattatcctcatactcgccatgatagttgaggtttgtcttaatcgctatccgcttagcattAACGCTTGCCAGTTATCTACTTTATCCGTTTGTATATACTCAATAAAGCGCACAAACTCGCTCTTAGTTGACGAGAAGAAtatata                                        >  1:822762‑M2/166‑251 (MQ=255)                atataatataactgcgcaaccggtattttgtaagttttgtgtaaatatttaaaagccccttaaacatcgtttaaagggcttttgaagttaaaatttataaaccagattatcctcatactcgccatgatagttgaggtttgtcttaatcgctatccgcttagcattAACGCTTGCCAGTTATCTACTTTATCCGTTTGTATATACTCAATAAAGCGCACAAACTCGCTCTTAGTTGACGAGAAGAAtatata                                        >  1:808279‑M2/166‑251 (MQ=255)                atataatataactgcgcaaccggtattttgtaagttttgtgtaaatatttaaaagccccttaaacatcgtttaaagggcttttgaagttaaaatttataaaccagattatcctcatactcgccatgatagttgaggtttgtcttaatcgctatccgcttagcattAACGCTTGCCAGTTATCTACTTTATCCGTTTGTATATACTCAATAAAGCGCACAAACTCGCTCTTAGTTGACGAGAAGAAtatata                                        >  1:767160‑M2/166‑251 (MQ=255)                 tataatataactgcgcaaccggtattttgtaagttttgtgtaaatatttaaaagccccttaaacatcgtttaaagggcttttgaagttaaaatttataaaccagattatcctcatactcgccatgatagttgaggtttgtcttaatcgctatccgcttagcattAACGCTTGCCAGTTATCTACTTTATCCGTTTGTATATACTCAATAAAGCGCACAAACTCGCTCTTAGTTGACGAGAAGAATATATAc                                       <  1:140861‑M2/87‑1 (MQ=255)                   taatataactgcgcaaccggtattttgtaagttttgtgtaaatatttaaaagccccttaaacatcgtttaaagggcttttgaagttaaaatttataaaccagattatcctcatactcgccatgatagttgaggtttgtcttaatcgctatccgcttagcattAACGCTTGCCAGTTATCTACTTTATCCGTTTGTATATACTCAATAAAGCGCACAAACTCGCTCTTAGTTGACGAGAAGAAtatata                                        <  1:628600‑M2/86‑1 (MQ=255)                   taatataactgcgcaaccggtattttgtaagttttgtgtaaatatttaaaagccccttaaacatcgtttaaagggcttttgaagttaaaatttataaaccagattatcctcatactcgccatgatagttgaggtttgtcttaatcgctatccgcttagcattAACGCTTGCCAGTTATCTACTTTATCCGTTTGTATATACTCAATAAAGCGCACAAACTCGCTCTTAGTTGACGAGAAGAAtatata                                        >  1:41569‑M2/163‑248 (MQ=255)                   taatataactgcgcaaccggtattttgtaagttttgtgtaaatatttaaaagccccttaaacatcgtttaaagggcttttgaagttaaaatttataaaccagattatcctcatactcgccatgatagttgaggtttgtcttaatcgctatccgcttagcattAACGCTTGCCAGTTATCTACTTTATCCGTTTGTATATACTCAATAAAGCGCACAAACTCGCTCTTAGTTGACGAGAAGAATATATAcgg                                     >  1:895333‑M2/163‑251 (MQ=255)                     atataactgcgcaaccggtattttgtaagttttgtgtaaatatttaaaagccccttaaacatcgtttaaagggcttttgaagttaaaatttataaaccagattatcctcatactcgccatgatagttgaggtttgtcttaatcgctatccgcttagcattAACGCTTGCCAGTTATCTACTTTATCCGTTTGTATATACTCAATAAAGCGCACAAACTCGCTCTTAGTTGACGAGAAGAAtatata                                        <  1:553053‑M2/86‑1 (MQ=255)                     atataactgcgcaaccggtattttgtaagttttgtgtaaatatttaaaagccccttaaacatcgtttaaagggcttttgaagttaaaatttataaaccagattatcctcatactcgccatgatagttgaggtttgtcttaatcgctatccgcttagcattAACGCTTGCCAGTTATCTACTTTATCCGTTTGTATATACTCAATAAAGCGCACAAACTCGCTCTTAGTTGACGAGAAGAATATATAcggcg                                   >  1:809348‑M2/161‑251 (MQ=255)                        taactgcgcaaccggtattttgtaagttttgtgtaaatatttaaaagccccttaaacatcgtttaaagggcttttgaagttaaaatttataaaccagattatcctcatactcgccatgatagttgaggtttgtcttaatcgctatccgcttagcattAACGCTTGCCAGTTATCTACTTTATCCGTTTGTATATACTCAATAAAGCGCACAAACTCGCTCTTAGTTGACGAGAAGAATATATAc                                       >  1:507963‑M2/158‑244 (MQ=255)                           ctgcgcaaccggtattttgtaagttttgtgtaaatatttaaaagccccttaaacatcgtttaaagggcttttgaagttaaaatttataaaccagattatcctcatactcgccatgatagttgaggtttgtcttaatcgctatccgcttagcattAACGCTTGCCAGTTATCTACTTTATCCGTTTGTATATACTCAATAAAGCGCACAAACTCGCTCTTAGTTGACGAGAAGAATATATAc                                       <  1:516350‑M2/87‑1 (MQ=255)                                 aaccggtattttgtaagttttgtgtaaatatttaaaagccccttaaacatcgtttaaagggcttttgaagttaaaatttataaaccagattatcctcatactcgccatgatagttgaggtttgtcttaatcgctatccgcttagcattAACGCTTGCCAGTTATCTACTTTATCCGTTTGTATATACTCAATAAAGCGCACAAACTCGCTCTTAGTTGACGAGAAGAATATATACGGCGGGCGGGTGATGt                       >  1:763292‑M2/149‑251 (MQ=255)                                                        gtaaatatttaaaagccccttaaacatcgtttaaagggcttttgaagttaaaatttataaaccagattatcctcatactcgccatgatagttgaggtttgtcttaatcgctatccgcttagcattAACGCTTGCCAGTTATCTACTTTATCCGTTTGTATATACTCAATAAAGCGCACAAACTCGCTCTTAGTTGACGAGAAGAAtatata                                        >  1:6909‑M2/126‑211 (MQ=255)                                                                       ccccttaaacatcgtttaaagggcttttgaagttaaaatttataaaccagattatcctcatactcgccatgatagttgaggtttgtcttaatcgctatccgcttagcattAACGCTTGCCAGTTATCTACTTTATCCGTTTGTATATACTCAATAAAGCGCACAAACTCGCTCTTAGTTGACGAGAAGAATATATACGGCGGGcgggtgatgttgactagcaattctttgacta  >  1:730979‑M2/111‑219 (MQ=255)                                                                        cccttaaacatcgtttaaagggcttttgaagttaaaatttataaaccagattatcctcatactcgccatgatagttgaggtttgtcttaatcgctatccgcttagcattAACGCTTGCCAGTTATCTACTTTATCCGTTTGTATATACTCAATAAAGCGCACAAACTCGCTCTTAGTTGACGAGAAGAATATATAc                                       <  1:261590‑M2/87‑1 (MQ=255)                                                                        cccttaaacatcgtttaaagggcttttgaagttaaaatttataaaccagattatcctcatactcgccatgatagttgaggtttgtcttaatcgctatccgcttagcattAACGCTTGCCAGTTATCTACTTTATCCGTTTGTATATACTCAATAAAGCGCACAAACTCGCTCTTAGTTGACGAGAAGAATATATACGGCGGGCGGGTGATGTTGACt                  <  1:372950‑M2/108‑1 (MQ=255)                                                                                                                       agattatcctcatactcgccatgatagttgaggtttgtcttaatcgctatccgcttagcattAACGCTTGCCAGTTATCTACTTTATCCGTTTGTATATACTCAATAAAGCGCACAAACTCGCTCTTAGTTGACGAGAAGAATATATAc                                       <  1:881290‑M2/87‑1 (MQ=255)                                                                                                                             tcctcatactcgccatgatagttgaggtttgtcttaatcgctatccgcttagcattAACGCTTGCCAGTTATCTACTTTATCCGTTTGTATATACTCAATAAAGCGCACAAACTCGCTCTTAGTTGACGAGAAGAAtatata                                        >  1:146535‑M2/57‑142 (MQ=255)                                                                                                                             tcctcatactcgccatgatagttgaggtttgtcttaatcgctatccgcttagcattAACGCTTGCCAGTTATCTACTTTATCCGTTTGTATATACTCAATAAAGCGCACAAACTCGCTCTTAGTTGACGAGAAGAATATATAcg                                      >  1:172254‑M2/57‑144 (MQ=255)                                                                                                                                         ccatgatagttgaggtttgtcttaatcgctatccgcttagcattAACGCTTGCCAGTTATCTACTTTATCCGTTTGTATATACTCAATAAAGCGCACAAACTCGCTCTTAGTTGACGAGAAGAATATATAc                                       <  1:290215‑M2/87‑1 (MQ=255)                                                                                                                                               tagttgaggtttgtcttaatcgctatccgcttagcattAACGCTTGCCAGTTATCTACTTTATCCGTTTGTATATACTCAATAAAGCGCACAAACTCGCTCTTAGTTGACGAGAAGAATATATAc                                       <  1:896805‑M2/87‑1 (MQ=255)                                                                                                                                                                                           tGCCAGTTATCTACTTTATCCGTTTGTATATACTCAATAAAGCGCACAAACTCGCTCTTAGTTGACgag                                                   <  1:612023/69‑1 (MQ=255)                                                                                                                                                                                                gTTATCTACTTTATCCGTTTGTATATACTCAATAAAGCGCACAAACTCGCTCTTAGTTGACGAGAAGAATATATAc                                       <  1:555716/76‑1 (MQ=255)                                                                                                                                                                                                gTTATCTACTTTATCCGTTTGTATATACTCAATAAAGCGCACAAACTCGCTCTTAGTTGACGAGAAGAATATATAc                                       <  1:320648/76‑1 (MQ=255)                                                                                                                                                                                                                                                   |                                                              AATGTTATATTATATAATATAACCGCGCAACCGCTATTTTGTAAGTTTTGTGTAAATATTTAAAAGCCCCTTAAACATCGTTTAAAGGGCTTTTGAAGTTAAAATTTATAAACCAGATTATCCTCATACTCGCCATGATAGTTGAGGTTTGTCTTAATCGCTATCCGCTTAGCATTGTCAAACGCTTGCCAGTTATCTACTTTATCCGTTTGTATATACTCAATAAAGCGCACAAACTCGCTTTTGGTTGACGAGAAGAAGATATACGGCGGGCGGGTGATGTTGACTAAGCGCAGGAAGTCGA  >  1/91482‑91785

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