Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	1	1,264,384	A→G	H211R (CAC→CGC)	unknown →	Mobile element protein

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	1	1,264,384	0	A	G	100.0%	84.5 / NA	29	H211R (CAC→CGC)	unknown	Mobile element protein
Reads supporting (aligned to +/- strand):  ref base A (0/0);  new base G (16/13);  total (16/13)

ATACAAAGTATGAGTCGCCGAGGGAATTGCTATGATAATGCGGTGATTGAAAGCTTTTTTGCGATATTAAAATCTGAGTGTTTTTACTCACACACTTATACTTCGATTGCCGAATTACAGGCGGAAATTGAAGAATATTTGGTGTATTACAACCAAGAACGAATTAAACTTGATTTAAAAGGATTAAGCCCGGTGCAATACCGAGCTCAATATTTAAGTTAACTAACCTGTCCAGCTTTTTGGGGGCAGATCACTTTTTGCAGCAAGACAAAAGTGCGGTAGGTTCTTACCGTGTTTTTTACTTCAACGGTCTAAACTTTACCTGACTGTATTGCAG  >  1/1264293‑1264629                                                                                            |                                                                                                                                                                                                                                                      aTACAAAGTATGAGTCGCCGAGGGAATTGCTATGATAATGCGGTGATTGAAAGCTTTTTTGCGATATTAAAATCTGAGTGTTTTTACTCACGCACTTATACTTCGATTGCCGAATTACAGGCGGAAATTGAAGAATATTTGGTGTATTACAACCAAGAACGAATTAAACTTGATTTAAAAGGATTAAGCCCGGTGCAATACCGAGCTCAATATTTAAGTTAACTAACCTGTCCAGCTTTTTGGGGGCAGAt                                                                                        >  1:579945/1‑251 (MQ=14)   aCAAAGTATGAGTCGCCGAGGGAATTGCTATGATAATGCGGTGATTGAAAGCTTTTTTGCGATATTAAAATCTGAGTGTTTTTACTCACGCACTTATACTTCGATTGCCGAATTACAGGCGGAAATTGAAGAATATTTGGTGTATTACAACCAAGAACGAATTAAACTTGATTTAAAAGGATTAAGCCCGGTGCAATACCGAGCTCAATATTTAAGTTAACTAACCTGTCCAGCTTTTTgg                                                                                                >  1:725956/1‑241 (MQ=14)      aaGTATGAGTCGCCGAGGGAATTGCTATGATAATGCGGTGATTGAAAGCTTTTTTGCGATATTAAAATCTGAGTGTTTTTACTCACGCACTTATACTTCGATTGCCGAATTACAGGCGGAAATTGAAGAATATTTGGTGTATTACAACCAAGAACGAATTAAACTTGATTTAAAAGGATTAAGCCCGGTGCAATACCGAGCTCAATATTTAAGTTAACTAACCTGTCCAGCTTTTTgg                                                                                                <  1:299756/238‑1 (MQ=14)      aaGTATGAGTCGCCGAGGGAATTGCTATGATAATGCGGTGATTGAAAGCTTTTTTGCGATATTAAAATCTGAGTGTTTTTACTCACGCACTTATACTTCGATTGCCGAATTACAGGCGGAAATTGAAGAATATTTGGTGTATTACAACCAAGAACGAATTAAACTTGATTTAAAAGGATTAAGCCCGGTGCAATACCGAGCTCAATATTTAAGTTAACTAACCTGTCCAGCTTTTTGGGGGCAGATCACtt                                                                                   >  1:140879/1‑251 (MQ=14)      aaGTATGAGTCGCCGAGGGAATTGCTATGATAATGCGGTGATTGAAAGCTTTTTTGCGATATTAAAATCTGAGTGTTTTTACTCACGCACTTATACTTCGATTGCCGAATTACAGGCGGAAATTGAAGAATATTTGGTGTATTACAACCAAGAACGAATTAAACTTGATTTAAAAGGATTAAGCCCGGTGCAATACCGAGCTCAATATTTAAGTTAACTAACCTGTCCAGCTTTTTGGGGGCAGATCACtt                                                                                   >  1:753827/1‑251 (MQ=14)      aaGTATGAGTCGCCGAGGGAATTGCTATGATAATGCGGTGATTGAAAGCTTTTTTGCGATATTAAAATCTGAGTGTTTTTACTCACGCACTTATACTTCGATTGCCGAATTACAGGCGGAAATTGAAGAATATTTGGTGTATTACAACCAAGAACGAATTAAACTTGATTTAAAAGGATTAAGCCCGGTGCAATACCGAGCTCAATATTTAAGTTAACTAACCTGTCCAGCTTTTTGGGGGCAGATCACtt                                                                                   >  1:558978/1‑251 (MQ=14)      aaGTATGAGTCGCCGAGGGAATTGCTATGATAATGCGGTGATTGAAAGCTTTTTTGCGATATTAAAATCTGAGTGTTTTTACTCACGCACTTATACTTCGATTGCCGAATTACAGGCGGAAATTGAAGAATATTTGGTGTATTACAACCAAGAACGAATTAAACTTGATTTAAAAGGATTAAGCCCGGTGCAATACCGAGCTCAATATTTAAGTTAACTAACCTGTCCAGCTTTTTGGGGGCAGATCACtt                                                                                   >  1:535027/1‑251 (MQ=14)      aaGTATGAGTCGCCGAGGGAATTGCCATGATAATGCGGTGATTGTACGCTTTTTTCCGATATTAAAATCTGAGTGTTTTTACTCACGCACTTATACTTCGATTGCCGAATTACAGGCGCAAATTGAAGAATATTTGGTGTATTACAACCAAGAACGAATTAAACTTGATTTAAAAGGATTAAGCCCGGTGCAATACCGACCCCAATATTTAAGTCAACTAACCTGTCCAGCTTTTTGGGGGCAGATCACtt                                                                                   >  1:383649/1‑251 (MQ=11)                            tGCTATGATAATGCGGTGATTGAAAGCTTTTTTGCGATATTAAAATCTGAGTGTTTTTACTCACGCACTTATACTTCGATTGCCGAATTACAGGCGGAAATTGAAGAATATTTGGTGTATTACAACCAAGAACGAATTAAACTTGATTTAAAAGGATTAAGCCCGGTGCAATACCGAGCTCAATATTTAAGTTAACTAACCTGTCCAGCTTTTTggg                                                                                               <  1:282018/217‑1 (MQ=14)                            tGCTATGATAATGCGGTGATTGAAAGCTTTTTTGCGATATTAAAATCTGAGTGTTTTTACTCACGCACTTATACTTCGATTGCCGAATTACAGGCGGAAATTGAAGAATATTTGGTGTATTACAACCAAGAACGAATTAAACTTGATTTAAAAGGATTAAGCCCGGTGCAATACCGAGCTCAATATTTAAGTTAACTAACCTGTCCAGCTTTTTGGGGGCAGATCACTTTTTGCAGCAAGACAAAAGTGCg                                                             >  1:600070/1‑251 (MQ=255)                                aTGATAATGCGGTGATTGAAAGCTTTTTTGCGATATTAAAATCTGAGTGTTTTTACTCACGCACTTATACTTCGATTGCCGAATTACAGGCGGAAATTGAAGAATATTTGGTGTATTACAACCAAGAACGAATTAAACTTGATTTAAAAGGATTAAGCCCGGTGCAATACCGAGCTCAATATTTAAGTTAACTAACCTGTCCAGCTTTTTGGGGGCAGATCACTTTTTGCAGCAAGACAAAAGTGCGGTAg                                                         <  1:583007/251‑1 (MQ=255)                                   aTAATGCGGTGATTGAAAGCTTTTTTGCGATATTAAAATCTGAGTGTTTTTACTCACGCACTTATACTTCGATTGCCGAATTACAGGCGGAAATTGAAGAATATTTGGTGTATTACAACCAAGAACGAATTAAACTTGATTTAAAAGGATTAAGCCCGGTGCAATACCGGGCTCAATATTTAAGTTAACTAACCTGTCCAGCTTTTTGGGGGCAGATCACTTTTTGCAGCAAGACAAAAGTGCGGTAGGtt                                                      >  1:434663/1‑251 (MQ=255)                                   aTAATGCGGTGATTGAAAGCTTTTTTGCGATATTAAAATCTGAGTGTTTTTACTCACGCACTTATACTTCGATTGCCGAATTACAGGCGGAAATTGAAGAATATTTGGTGTATTACAACCAAGAACGAATTAAACTTGATTTAAAAGGATTAAGCCCGGTGCAATACCGAGCTCAATATTTAAGTTAACTAACCTGTCCAGCTTTTTGGGGGCAGATCACTTTTTGCAGCAAGACAAAAGTGCGGTAGGt                                                       <  1:628003/250‑1 (MQ=255)                                         cGGTGATTGAAAGCTTTTTTGCGATATTAAAATCTGAGTGTTTTTACTCACGCACTTATACTTCGATTGCCGAATTACAGGCGGAAATTGAAGAATATTTGGTGTATTACAACCAAGAACGAATTAAACTTGATTTAAAAGGATTAAGCCCGGTGCAATACCGAGCTCAATATTTAAGTTAACTAACCTGTCCAGCTTTTTGGGGGCAGATCACTTTTTGCAGCAAGACAAAAGTGCGGTAGGTTCTTAcc                                                <  1:772890/251‑1 (MQ=255)                                         cGGTGATTGAAAGCTTTTTTGCGATATTAAAATCTGAGTGTTTTTACTCACGCACTTATACTTCGATTGCCGAATTACAGGCGGAAATTGAAGAATATTTGGTGTATTACAACCAAGAACGAATTAAACTTGATTTAAAAGGATTAAGCCCGGTGCAATACCGAGCTCAATATTTAAGTTAACTAACCTGTCCAGCTTTTTGGGGGCAGATCACTTTTTGCAGCAAGACAAAAGTGCGGTAGGTTCTTAc                                                 <  1:232475/250‑1 (MQ=255)                                                     gCTTTTTTGCGATATTAAAATCTGAGTGTTTTTACTCACGCACTTATACTTCGATTGCCGAATTACAGGCGGAAATTGAAGAATATTTGGTGTATTACAACCAAGAACGAATTAAACTTGATTTAAAAGGATTAAGCCCGGTGCAATACCGAGCTCAATATTTAAGTTAACTAACCTGTCCAGCTTTTTGGGGGCAGATCACTTTTTGCAGCAAGACAAAAGTGCGGTAGGTTCTTAc                                                 >  1:132021/1‑238 (MQ=255)                                                            tGCGATATTAAAATCTGAGTGTTTTTACTCACGCACTTATACTTCGATTGCCGAATTACAGGCGGAAATTGAAGAATATTTGGTGTATTACAACCAAGAACGAATTAAACTTGATTTAAAAGGATTAAGCCCGGTGCAATACCGAGCTCAATATTTAAGTTAACTAACCTGTCCAGCTTTTTGGGGGCAGATCACTTTTTGCAGCAAGACAAAAGTGCGGTAGGTTCTTACCGTGTTTTTTACTTCAACgg                             <  1:773148/251‑1 (MQ=255)                                                                atatTAAAATCTGAGTGTTTTTACTCACGCACTTATACTTCGATTGCCGAATTACAGGCGGAAATTGAAGAATATTTGGTGTATTACAACCAAGAACGAATTAAACTTGATTTAAAAGGATTAAGCCCGGTGCAATACCGAGCTCAATATTTAAGTTAACTAACCTGTCCAGCTTTTTGGGGGCAGATCACTTTTTGCAGCAAGACAAAAGTGCGGTAGGTTCTTACCGTGTTTTTTACTTCAACGGTCTa                         >  1:87863/1‑251 (MQ=255)                                                                  atTAAAATCTGAGTGTTTTTACTCACGCACTTATACTTCGATTGCCGAATTACAGGCGGAAATTGAAGAATATTTGGTGTATTACAACCAAGAACGAATTAAACTTGATTTAAAAGGATTAAGCCCGGTGCAATACCGAGCTCAATATTTAAGTTAACTAACCTGTCCAGCTTTTTGGGGGCAGATCACTTTTTGCAGCAAGACAAAAg                                                                 >  1:362697/1‑209 (MQ=255)                                                                  atTAAAATCTGAGTGTTTTTACTCACGCACTTATACTTCGATTGCCGAATTACAGGCGGAAATTGAAGAATATTTGGTGTATTACAACCAAGAACGAATTAAACTTGATTTAAAAGGATTAAGCCCGGTGCAATACCGAGCTCAATATTTAAGTTAACTAACCTGTCCAGCTTTTTGGGGGCAGATCACTTTTTGCAGCAAGACAAAAGTg                                                               <  1:638822/211‑1 (MQ=255)                                                                  atTAAAATCTGAGTGTTTTTACTCACGCACTTATACTTCGATTGCCGAATTACAGGCGGAAATTGAAGAATATTTGGTGTATTACAACCAAGAACGAATTAAACTTGATTTAAAAGGATTAAGCCCGGTGCAATACCGAGCTCAATATTTAAGTTAACTAACCTGTCCAGCTTTTTGGGGGCAGATCACTTTTTGCAGCAAGACAAAAGTGCGGTAGGTTCTTACCg                                               >  1:658604/1‑227 (MQ=255)                                                                  atTAAAATCTGAGTGTTTTTACTCACGCACTTATACTTCGATTGCCGAATTACAGGCGGAAATTGAAGAATATTTGGTGTATTACAACCAAGAACGAATTAAACTTGATTTAAAAGGATTAAGCCCGGTGCAATACCGAGCTCAATATTTAAGTTAACTAACCTGTCCAGCTTTTTGGGGGCAGATCACTTTTTGCAGCAAGACAAAAGTGCGGTAGGTTCTTACCGTGTTTTTTACTTCAACGGTCTaaa                       >  1:600346/1‑251 (MQ=255)                                                                  atTAAAATCTGAGTGTTTTTACTCACGCACTTATACTTCGATTGCCGAATTACAGGCGGAAATTGAAGAATATTTGGTGTATTACAACCAAGAACGAATTAAACTTGATTTAAAAGGATTAAGCCCGGTGCAATACCGAGCTCAATATTTAAGTTAACTAACCTGTCCAGCTTTTTGGGGGCAGATCACTTTTTGCAGCAAGACAAAAGTGCGGTAGGTTCTTACCGTGTTTTTTACTTCAACGGTCTaaa                       >  1:664666/1‑251 (MQ=255)                                                                    tAAAATCTGAGTGTTTTTACTCACGCACTTATACTTCGATTGCCGAATTACAGGCGGAAATTGAAGAATATTTGGTGTATTACAACCAAGAACGAATTAAACTTGATTTAAAAGGATTAAGCCCGGTGCAATACCGAGCTCAATATTTAAGTTAACTAACCTGTCCAGCTTTTTGGGGGCAGATCACTTTTTGCAGCAAGACAAAAGTGCGGTAGGTTCTTACCGTGTTTTTTACTTCAACGGTCTAAACt                     <  1:155646/251‑1 (MQ=255)                                                                         tCTGAGTGTTTTTACTCACGCACTTATACTTCGATTGCCGAATTACAGGCGGAAATTGAAGAATATTTGGTGTATTACAACCAAGAACGAATTAAACTTGATTTAAAAGGATTAAGCCCGGTGCAATACCGAGCTCAATATTTAAGTTAACTAACCTGTCCAGCTTTTTGGGGGCAGATCACTTTTTGCAGCAAGACAAAAGTGCGGTAGGTTCTTACCGTGTTTTTTACTTCAACGGTCTAAACTTTAcc                <  1:880722/251‑1 (MQ=255)                                                                              gtgtTTTTACTCACGCACTTATACTTCGATTGCCGAATTACAGGCGGAAATTGAAGAATATTTGGTGTATTACAACCAAGAACGAATTAAACTTGATTTAAAAGGATTAAGCCCGGTGCAATACCGAGCTCAATATTTAAGTTAACTAACCTGTCCAGCTTTTTGGGGGCAGATCACTTTTTGCAGCAAGACAAAAGTGCGGTAGGTTCTTACCGTGTTTTTTACTTCAACGGTCTAAACTTTACCTGa             >  1:558036/1‑249 (MQ=255)                                                                                     tACTCACGCACTTATACTTCGATTGCCGAATTACAGGCGGAAATTGAAGAATATTTGGTGTATTACAACCAAGAACGAATTAAACTTGATTTAAAAGGATTAAGCCCGGTGCAATACCGAGCTCAATATTTAAGTTAACTAACCTGTCCAGCTTTTTGGGGGCAGATCACTTTTTGCAGCAAGACAAAAGTGCGGTAGGTTCTTACCGTGTTTTTTACTTCAACGGTCTAAACTTTACCTGACTGTATTGc    <  1:19603/251‑1 (MQ=255)                                                                                     tACTCACGCACTTATACTTCGATTGCCGAATTACAGGCGGAAATTGAAGAATATTTGGTGTATTACAACCAAGAACGAATTAAACTTGATTTAAAAGGATTAAGCCCGGTGCAATACCGAGCTCAATATTTAAGTTAACTAACCTGTCCAGCTTTTTGGGGGCAGATCACTTTTTGCAGCAAGACAAAAGTGCGGTAGGTTCTTACCGTGTTTTTTACTTCAACGGTCTAAACTTTACCTGACTGTATTGc    <  1:436228/251‑1 (MQ=255)                                                                                       cTCACGCACTTATACTTCGATTGCCGAATTACAGGCGGAAATTGAAGAATATTTGGTGTATTACAACCAAGAACGAATTAAACTTGATTTAAAAGGATTAAGCCCGGTGCAATACCGAGCTCAATATTTAAGTTAACTAACCTGTCCAGCTTTTTGGGGGCAGATCACTTTTTGCAGCAAGACAAAAGTGCGGTAGGTTCTTACCGTGTTTTTTACTTCAACGGTCTAAACTTTACCTGACTGTATTGCAg  <  1:532551/251‑1 (MQ=255)                                                                                            |                                                                                                                                                                                                                                                      ATACAAAGTATGAGTCGCCGAGGGAATTGCTATGATAATGCGGTGATTGAAAGCTTTTTTGCGATATTAAAATCTGAGTGTTTTTACTCACACACTTATACTTCGATTGCCGAATTACAGGCGGAAATTGAAGAATATTTGGTGTATTACAACCAAGAACGAATTAAACTTGATTTAAAAGGATTAAGCCCGGTGCAATACCGAGCTCAATATTTAAGTTAACTAACCTGTCCAGCTTTTTGGGGGCAGATCACTTTTTGCAGCAAGACAAAAGTGCGGTAGGTTCTTACCGTGTTTTTTACTTCAACGGTCTAAACTTTACCTGACTGTATTGCAG  >  1/1264293‑1264629

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