Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	2	147,875	T→G	D51E (GAT→GAG)	unknown →	hypothetical protein

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	2	147,875	0	T	G	100.0%	94.0 / NA	33	D51E (GAT→GAG)	unknown	hypothetical protein
Reads supporting (aligned to +/- strand):  ref base T (0/0);  new base G (22/11);  total (22/11)

AATTAGAGTTGCAAATTAGAACTCAGTTACAACATGAGTGGGCAACCATAGTTGAAATATTGGGCTCTTTGATGAAAACATCATTTAAAGTTGGACAAGGTGATGAGGGAATAAGATTATTTTTGAAATTATGTAGTGCATTATTTTCTTATCAGGAAGATACTCCAATAATAGAAGAGTTTAATAATTTTTCTGTTGTGCAAATATGTGAAAAAATAAAAGATTTAGATTCAAAATTAGGTGTAATAGATAAACTAAAAAGCGTTAGTGCTATTGATGTAGAACATAACCCTAGCGCAGCATATTATTTGTTAATATTAGATTTAAAAAAACATAATACTGAAATTAAAACTTTTAATGAGAATGAGAAATTTAAAGCTGAGCAGGTATACTCCCATTATGAACAAGATTCTTTTAAG  >  2/147652‑148070                                                                                                                                                                                                                                |                                                                                                                                                                                                    aaTTAGAGTTGCAAATTAGAACTCAGTTACAACATGAGTGGGCAACCATAGTTGAAATATTGGGCTCTTTGATGAAAACATCATTTAAAGTTGGACAAGGTGATGAGGGAATAAGATTATTTTTGAAATTATGTAGTGCATTATTTTCTTATCAGGAAGATACTCCAATAATAGAAGAGTTTAATAATTTTTCTGTTGTGCAAATATGTGAAAAAATAAAAGAGTTAGATTCAAAATTAGGTGTAATAGAt                                                                                                                                                                          >  1:738530/1‑251 (MQ=255) aaTTAGAGTTGCAAATTAGAACTCAGTTACAACATGAGTGGGCAACCATAGTTGAAATATTGGGCTCTTTGATGAAAACATCATTTAAAGTTGGACAAGGTGATGAGGGAATAAGATTATTTTTGAAATTATGTAGTGCATTATTTTCTTATCAGGAAGATACTCCAATAATAGAAGAGTTTAATAATTTTTCTGTTGTGCAAATATGTGAAAAAATAAAAGAGTTAGATTCAAAATTAGGTGTAATAGAt                                                                                                                                                                          <  1:631028/251‑1 (MQ=255)                     aCTCAGTTACAACATGAGTGGGCAACCATAGTTGAAATATTGGGCTCTTTGATGAAAACATCATTTAAAGTTGGACAAGGTGATGAGGGAATAAGATTATTTTTGAAATTATGTAGTGCATTATTTTCTTATCAGGAAGATACTCCAATAATAGAAGAGTTTAATAATTTTTCTGTTGTGCAAATATGTGAAAAAATAAAAGAGTTAGATTCAAAATTAGGTGTAATAGATAAACTAAAAAGCg                                                                                                                                                             <  1:641728/244‑1 (MQ=255)                        cAGTTACAACATGAGTGGGCAACCATAGTTGAAATATTGGGCTCTTTGATGAAAACATCATTTAAAGTTGGACAAGGTGATGAGGGAATAAGATTATTTTTGAAATTATGTAGTGCATTATTTTCTTATCAGGAAGATACTCCAATAATAGAAGAGTTTAATAATTTTTCTGTTGTGCAAATATGTGAAAAAATAAAAGAGTTAGATTCAAAATTAGGTGTAATAGATAAACTAAAAAGCGTTAGTg                                                                                                                                                       <  1:530787/247‑1 (MQ=255)                          gTTACAACATGAGTGGGCAACCATAGTTGAAATATTGGGCTCTTTGATGAAAACATCATTTAAAGTTGGACAAGGTGATGAGGGAATAAGATTATTTTTGAAATTATGTAGTGCATTATTTTCTTATCAGGAAGATACTCCAATAATAGAAGAGTTTAATAATTTTTCTGTTGTGCAAATATGTGAAAAAATAAAAGAGTTAGATTCAAAATTAGGTGTAATAg                                                                                                                                                                            >  1:69353/1‑224 (MQ=255)                          gTTACAACATGAGTGGGCAACCATAGTTGAAATATTGGGCTCTTTGATGAAAACATCATTTAAAGTTGGACAAGGTGATGAGGGAATAAGATTATTTTTGAAATTATGTAGTGCATTATTTTCTTATCAGGAAGATACTCCAATAATAGAAGAGTTTAATAATTTTTCTGTTGTGCAAATATGTGAAAAAATAAAAGAGTTAGATTCAAAATTAGGTGTAATAGATAAACTAAAAAGCGTTAGTGCTATTg                                                                                                                                                 >  1:321221/1‑251 (MQ=255)                          gTTACAACATGAGTGGGCAACCATAGTTGAAATATTGGGCTCTTTGATGAAAACATCATTTAAAGTTGGACAAGGTGATGAGGGAATAAGATTATTTTTGAAATTATGTAGTGCATTATTTTCTTATCAGGAAGATACTCCAATAATAGAAGAGTTTAATAATTTTTCTGTTGTGCAAATATGTGAAAAAATAAAAGAGTTAGATTCAAAATTAGGTGTAATAGATAAACTAAAAAGCGTTAGTGCTATTg                                                                                                                                                 >  1:382253/1‑251 (MQ=255)                           ttACAACATGAGTGGGCAACCATAGTTGAAATATTGGGCTCTTTGATGAAAACATCATTTAAAGTTGGACAAGGTGATGAGGGAATAAGATTATTTTTGAAATTATGTAGTGCATTATTTTCTTATCAGGAAGATACTCCAATAATAGAAGAGTTTAATAATTTTTCTGTTGTGCAAATATGTGAAAAAATAAAAGAGTTAGATTCAAAATTAGGTGTAATAGATAAACTAAAAAGCGTTAGTGCTATTg                                                                                                                                                 <  1:65236/250‑1 (MQ=255)                              caacaTGAGTGGGCAACCATAGTTGAAATATTGGGCTCTTTGATGAAAACATCATTTAAAGTTGGACAAGGTGATGAGGGAATAAGATTATTTTTGAAATTATGTAGTGCATTATTTTCTTATCAGGAAGATACTCCAATAATAGAAGAGTTTAATAATTTTTCTGTTGTGCAAATATGTGAAAAAATAAAAGAGTTAGATTCAAAATTAGGTGTAATAGATAAACTAAAAAGCGTTAGTGCTATTGATGt                                                                                                                                             >  1:182696/1‑251 (MQ=255)                                acaTGAGTGGGCAACCATAGTTGAAATATTGGGCTCTTTGATGAAAACATCATTTAAAGTTGGACAAGGTGATGAGGGAATAAGATTATTTTTGAAATTATGTAGTGCATTATTTTCTTATCAGGAAGATACTCCAATAATAGAAGAGTTTAATAATTTTTCTGTTGTGCAAATATGTGAAAAAATAAAAGAGTTAGATTCAAAATTAGGTGTAATAGATAAACTAAAAAGCGTTAGTGCTATTGATGt                                                                                                                                             <  1:267522/249‑1 (MQ=255)                                   tGAGTGGGCAACCATAGTTGAAATATTGGGCTCTTTGATGAAAACATCATTTAAAGTTGGACAAGGTGATGAGGGAATAAGATTATTTTTGAAATTATGTAGTGCATTATTTTCTTATCAGGAAGATACTCCAATAATAGAAGAGTTTAATAATTTTTCTGTTGTGCAAATATGTGAAAAAATAAAAGAGTTAGATTCAAAATTAGGTGTAATAGATAAACTAAAAAGCGTTAGTGCTATTGATGTAGAAc                                                                                                                                        <  1:367770/251‑1 (MQ=255)                                        gggCAACCATAGTTGAAATATTGGGCTCTTTGATGAAAACATCATTTAAAGTTGGACAAGGTGATGAGGGAATAAGATTATTTTTGAAATTATGTAGTGCATTATTTTCTTATCAGGAAGATACTCCAATAATAGAAGAGTTTAATAATTTTTCTGTTGTGCAAATATGTGAAAAAATAAAAGAGTTAGATTCAAAATTAGGTGTAATAGATAAACTAAAAAGCGTTAGTGCTATTGATGTAGAACATaa                                                                                                                                    >  1:842812/1‑250 (MQ=255)                                              ccATAGTTGAAATATTGGGCTCTTTGATGAAAACATCATTTAAAGTTGGACAAGGTGATGAGGGAATAAGATTATTTTTGAAATTATGTAGTGCATTATTTTCTTATCAGGAAGATACTCCAATAATAGAAGAGTTTAATAATTTTTCTGTTGTGCAAATATGTGAAAAAATAAAAGAGTTAGATTCAAAATTAGGTGTAATAGATAAACTAAAAAGCGTTAGTGCTATTGATGTAGAACATAACTCTAg                                                                                                                              >  1:552175/1‑250 (MQ=255)                                              ccATAGTTGAAATATTGGGCTCTTTGATGAAAACATCATTTAAAGTTGGACAAGGTGATGAGGGAATAAGATTATTTTTGAAATTATGTAGTGCATTATTTTCTTATCAGGAAGATACTCCAATAATAGAAGAGTTTAATAATTTTTCTGTTGTGCAAATATGTGAAAAAATAAAAGAGTTAGATTCAAAATTAGGTGTAATAGATAAACTAAAAAGCGTTAGTGCTATTGATGTAGAACATAACTCTAg                                                                                                                              >  1:397539/1‑250 (MQ=255)                                                      gAAATATTGGGCTCTTTGATGAAAACATCATTTAAAGTTGGACAAGGTGATGAGGGAATAAGATTATTTTTGAAATTATGTAGTGCATTATTTTCTTATCAGGAAGATACTCCAATAATAGAAGAGTTTAATAATTTTTCTGTTGTGCAAATATGTGAAAAAATAAAAGAGTTAGATTCAAAATTAGGTGTAATAGATAAACTAAAAAGCGTTAGTGCTATTGATGTAGAACATAACTCTAgcgc                                                                                                                           >  1:790113/1‑245 (MQ=255)                                                         atatTGGGCTCTTTGATGAAAACATCATTTAAAGTTGGACAAGGTGATGAGGGAATAAGATTATTTTTGAAATTATGTAGTGCATTATTTTCTTATCAGGAAGATACTCCAATAATAGAAGAGTTTAATAATTTTTCTGTTGTGCAAATATGTGAAAAAATAAAAGAGTTAGATTCAAAATTAGGTGTAATAGATAAACTAAAAAGCGTTAGTGCTATTGATGTAGAACATAACTCTAGCGCAGCAtatta                                                                                                                  >  1:662483/1‑251 (MQ=255)                                                          tatTGGGCTCTTTGATGAAAACATCATTTAAAGTTGGACAAGGTGATGAGGGAATAAGATTATTTTTGAAATTATGTAGTGCATTATTTTCTTATCAGGAAGATACTCCAATAATAGAAGAGTTTAATAATTTTTCTGTTGTGCAAATATGTGAAAAAATAAAAGAGTTAGATTCAAAATTAGGTGTAATAGATAAACTAAAAAGCGTTAGTGCTATTGATGTAGAACATAACTCTAGCGCAGCAtattat                                                                                                                 <  1:417914/251‑1 (MQ=255)                                                                                                   ggTGATGAGGGAATAAGATTATTTTTGAAATTATGTAGTGCATTATTTTCTTATCAGGAAGATACTCCAATAATAGAAGAGTTTAATAATTTTTCTGTTGTGCAAATATGTGAAAAAATAAAAGAGTTAGATTCAAAATTAGGTGTAATAGATAAACTAAAAAGCGTTAGTGCTATTGATGTAGAACATAAc                                                                                                                                   <  1:635820/192‑1 (MQ=255)                                                                                                                       attTTTGAAATTATGTAGTGCATTATTTTCTTATCAGGAAGATACTCCAATAATAGAAGAGTTTAATAATTTTTCTGTTGTGCAAATATGTGAAAAAATAAAAGAGTTAGATTCAAAATTAGGTGTAATAGATAAACTAAAAAGCGTTAGTGCTATTg                                                                                                                                                 >  1:772825/1‑158 (MQ=255)                                                                                                                              aaaTTATGTAGTGCATTATTTTCTTATCAGGAAGATACTCCAATAATAGAAGAGTTTAATAATTTTTCTGTTGTGCAAATATGTGAAAAAATAAAAGAGTTAGATTCAAAATTAGGTGTAATAGATAAACTAAAAAGCGTTAGTGCTATTGATGTAGAACATAACTCTAgcgc                                                                                                                           >  1:856748/1‑173 (MQ=255)                                                                                                                              aaaTTATGTAGTGCATTATTTTCTTATCAGGAAGATACTCCAATAATAGAAGAGTTTAATAATTTTTCTGTTGTGCAAATATGTGAAAAAATAAAAGAGTTAGATTCAAAATTAGGTGTAATAGATAAACTAAAAAGCGTTAGTGCTATTGATGTAGAACATAACTCTAGCGCAGCAtattat                                                                                                                 <  1:807560/183‑1 (MQ=255)                                                                                                                              aaaTTATGTAGTGCATTATTTTCTTATCAGGAAGATACTCCAATAATAGAAGAGTTTAATAATTTTTCTGTTGTGCAAATATGTGAAAAAATAAAAGAGTTAGATTCAAAATTAGGTGTAATAGATAAACTAAAAAGCGTTAGTGCTATTGATGTAGAACATAACTCTAGCGCAGCATATTATTTGTTAATATTAGATTTAAAAAAACATAATACTGAAATTAAAACTTTTAATGAGAATGAGAAATTTaa                                             >  1:491089/1‑251 (MQ=255)                                                                                                                                     gTAGTGCATTATTTTCTTATCAGGAAGATACTCCAATAATAGAAGAGTTTAATAATTTTTCTGTTGTGCAAATATGTGAAAAAATAAAAGAGTTAGATTCAAAATTAGGTGTAATAGATAAACTAAAAAGCGTTAGTGCTATTGATGTAGAACATAACTCTAGCGCAGCAtattat                                                                                                                 <  1:278819/176‑1 (MQ=255)                                                                                                                                            attattTTCTTATCAGGAAGATACTCCAATAATAGAAGAGTTTAATAATTTTTCTGTTGTGCAAATATGTGAAAAAATAAAAGAGTTAGATTCAAAATTAGGTGTAATAGATAAACTAAAAAGCGTTAGTGCTATTGATGTAGAACATAACTCTAGCGCAGCAtattat                                                                                                                 >  1:512116/1‑169 (MQ=255)                                                                                                                                                                   cTCCAATAATAGAAGAGTTTAATAATTTTTCTGTTGTGCAAATATGTGAAAAAATAAAAGAGTTAGATTCAAAATTAGGTGTAATAGATAAACTAAAAAGCGTTAGTGCTATTGATGTAGAACATAACTCTAGCGCAGCAtattat                                                                                                                 >  1:159244/1‑146 (MQ=255)                                                                                                                                                                    tCCAATAATAGAAGAGTTTAATAATTTTTCTGTTGTGCAAATATGTGAAAAAATAAAAGAGTTAGATTCAAAATTAGGTGTAATAGATAAACTAAAAAGCGTTAGTGCTATTGATGTAGAACATAACTCTAGCGCAGCAtattat                                                                                                                 >  1:94131/1‑145 (MQ=255)                                                                                                                                                                     ccAATAATAGAAGAGTTTAATAATTTTTCTGTTGTGCAAATATGTGAAAAAATAAAAGAGTTAGATTCAAAATTAGGTGTAATAGAt                                                                                                                                                                          >  1:874222/1‑87 (MQ=255)                                                                                                                                                                         taataGAAGAGTTTAATAATTTTTCTGTTGTGCAAATATGTGAAAAAATAAAAGAGTTAGATTCAAAATTAGGTGTAATAGATAAACTAAAAAGCGTTAGTGCTATTGATGTAGAACATAActct                                                                                                                                >  1:565176/1‑122 (MQ=255)                                                                                                                                                                         taataGAAGAGTTTAATAATTTTTCTGTTGTGCAAATATGTGAAAAAATAAAAGAGTTAGATTCAAAATTAGGTGTAATAGATAAACTAAAAAGCGTTAGTGCTATTGATGTAGAACATAACTCTAGCGCAGCATATTATTTGTTAATATTAGATTTAAAAAAACATAATACTGAAATTAAAACTTTTAATGAGAATGAGAAATTTAAAGCTGAGCAGGTATACTCCCATTATGAACAAGATTCTTTTAAg  >  1:148895/1‑251 (MQ=255)                                                                                                                                                                           ataGAAGAGTTTAATAATTTTTCTGTTGTGCAAATATGTGAAAAAATAAAAGAGTTAGATTCAAAATTAGGTGTAATAGATAAACTAAAAAGCGTTAGTGCTATTGATGTAGaa                                                                                                                                         <  1:154492/114‑1 (MQ=255)                                                                                                                                                                                         taatTTTTCTGTTGTGCAAATATGTGAAAAAATAAAAGAGTTAGATTCAAAATTAGGTGTAATAGATAAACTAAAAAGCGTTAGTGCTATTGATGTAGAACATAACTCTAGCGCAGCATATTATTTGTTAATATTAGATTTAAAAAAACATAATACTGAAATTAAAACTTTTAATGAGAATGAGAAATTTAAAGCTGAGCAGGTATACTc                           >  1:583978/1‑210 (MQ=255)                                                                                                                                                                                                   gttgtGCAAATATGTGAAAAAATAAAAGAGTTAGATTCAAAATTAGGTGTAATAGATAAACTAAAAAGCGTTAGTGCTATTGATGTAGAACATAACTCTAGCGCAGCAtattat                                                                                                                 >  1:760546/1‑114 (MQ=255)                                                                                                                                                                                                   gttgtGCAAATATGTGAAAAAATAAAAGAGTTAGATTCAAAATTAGGTGTAATAGATAAACTAAAAAGCGTTAGTGCTATTGATGTAGAACATAACTCTAGCGCAGCATATTATTTGTTAATATTAGATTTAAAAAAACATAATACTGAAATTAAAACTTTTAATGAGAATGAGAAATTTAAAGCTGAGCAGGta                                >  1:58929/1‑195 (MQ=255)                                                                                                                                                                                                                                |                                                                                                                                                                                                    AATTAGAGTTGCAAATTAGAACTCAGTTACAACATGAGTGGGCAACCATAGTTGAAATATTGGGCTCTTTGATGAAAACATCATTTAAAGTTGGACAAGGTGATGAGGGAATAAGATTATTTTTGAAATTATGTAGTGCATTATTTTCTTATCAGGAAGATACTCCAATAATAGAAGAGTTTAATAATTTTTCTGTTGTGCAAATATGTGAAAAAATAAAAGATTTAGATTCAAAATTAGGTGTAATAGATAAACTAAAAAGCGTTAGTGCTATTGATGTAGAACATAACCCTAGCGCAGCATATTATTTGTTAATATTAGATTTAAAAAAACATAATACTGAAATTAAAACTTTTAATGAGAATGAGAAATTTAAAGCTGAGCAGGTATACTCCCATTATGAACAAGATTCTTTTAAG  >  2/147652‑148070

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