Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	1	1,285,295	A→T	intergenic (+447/+54)	unknown → / ← unknown	hypothetical protein/hypothetical protein

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	1	1,285,295	0	A	T	100.0%	59.4 / NA	22	intergenic (+447/+54)	unknown/unknown	hypothetical protein/hypothetical protein
Reads supporting (aligned to +/- strand):  ref base A (0/0);  new base T (14/8);  total (14/8)

TTAGGTGTTAGAACCTTAGATATCTGATGGCACTGGAGTAAAAACCTGATATAGATCGTGTGTTGAACGATGTTGTACATTTGATGCTAAAGTGCCCTACTTGGGTATGCTAATTACGTAAGAGCGGTAATTATATAAAACGAATTTAAATAAACAGAAAATTATATGAAATGGGGAAGTATCCCACTCGATTATTTATCGGCGATATAATGCAATGATTTATAAATGAATAAAAATGATATTGTTTATTCTCTACCGTCATTTTTACCGTCACCAGAAGAAAATGCCATCAAGAACGTTCAGCAAGCCAAGTATCAATATCAGATTTTCTCCACGCTACACGGCCCATAGACGCTTTAACACTTTTGGGGAATTTACCTTGTTTAATCATACGCCAAATTGTTGTGTCGCTCAAGCTGGTGACTTCAATCACCTTTGCTTTACTTAGTAGCGCTGGTTGTTGTTCCATTATTCTTCTC  >  1/1285060‑1285538                                                                                                                                                                                                                                            |                                                                                                                                                                                                                                                    ttAGGTGTTAGAACCTTAGATATCTGATGGCACTGGAGTAAAAACCTGATATAGATCGTGTGTTGAACGATGTTGTACATTTGATGCTAAAGTGCCCTACTTGGGTATGCTAATTACGTAAGAGCGGTAATTATATAAAACGAATTTAAATAAACAGAAAATTATATGAAATGGGGAAGTATCCCACTCGATTATTTATCGGCGACATAATGCAATGATTTATCAATGAATAAAATTGATATTGTTTATtc                                                                                                                                                                                                                                      <  1:282296/251‑1 (MQ=255)                                                            gtgtTGAACGATGTTGTACATTTGATGCTAAAGTGCCCTACTTGGGTATGCTAATTACGTAAGAGCGGTAATTATATAAAACGAATTTAAATAAACAGAAAATTATATGAAATGGGGAAGTATCCCACTCGATTATTTATCGGCGACATAATGCAATGATTTATCAATGAATAAAATTGATATTGTTTATTCTATACCGTCATTTTTACCGTCTTTGGAAGAAAATGCCATCAAGAACGTTCAGCAAGCCa                                                                                                                                                                           >  1:139352/1‑251 (MQ=255)                                                                              cATTTGATGCTAAAGTGCCCTACTTGGGTATGCTAATTACGTAAGAGCGGTAATTATATAAAACGAATTTAAATAAACAGAAAATTATATGAAATGGGGAAGTATCCCACTCGATTATTTATCGGCGACATAATGCAATGATTTATCAATGAATAAAATTGATATTGTTTATTCTATACCGTCATTTTTACCGTCTTTGGAAGAAAATGCCATCAAGAACGTTCAGCAAGCCAAGTATCAatat                                                                                                                                                                <  1:663575/244‑1 (MQ=255)                                                                                    aTGCTAAAGTGCCCTACTTGGGTATGCTAATTACGTAAGAGCGGTAATTATATAAAACGAATTTAAATAAACAGAAAATTATATGAAATGGGGAAGTATCCCACTCGATTATTTATCGGCGACATAATGCAATGATTTATCAATGAATAAAATTGATATTGTTTATTCTATACCGTCATTTTTACCGTCTTTGGAAGAAAATGCCATCAAGAACGTTCAGCAAGCCAAGTATCAATATCAGATTTTCTCCa                                                                                                                                                   >  1:51365/1‑251 (MQ=255)                                                                                       cTAAAGTGCCCTACTTGGGTATGCTAATTACGTAAGAGCGGTAATTATATAAAACGAATTTAAATAAACAGAAAATTATATGAAATGGGGAAGTATCCCACTCGATTATTTATCGGCGACATAATGCAATGATTTATCAATGAATAAAATTGATATTGTTTATTCTATACCGTCATTTTTACCGTCTTTGGAAGAAAATGCCATCAAGAACGTTCAGCAAGCCAAGTATCAATATCAGATTTTCTCCACg                                                                                                                                                 >  1:209946/1‑250 (MQ=255)                                                                                               cccTACTTGGGTATGCTAATTACGTAAGAGCGGTAATTATATAAAACGAATTTAAATAAACAGAAAATTATATGAAATGGGGAAGTATCCCACTCGATTATTTATCGGCGACATAATGCAATGATTTATCAATGAATAAAATTGATATTGTTTATTCTATACCGTCATTTTTACCGTCTTTGGAAGAAAATGCCATCAAGAACGTTCAGCAAGCCAAGTATCAATATCAGATTTTCTCCACGCTACACgcc                                                                                                                                        <  1:189809/251‑3 (MQ=255)                                                                                                 cTACTTGGGTATGCTAATTACGTAAGAGCGGTAATTATATAAAACGAATTTAAATAAACAGAAAATTATATGAAATGGGGAAGTATCCCACTCGATTATTTATCGGCGACATAATGCAATGATTTATCAATGAATAAAATTGATATTGTTTATtc                                                                                                                                                                                                                                      <  1:737034/155‑1 (MQ=255)                                                                                                 cTACTTGGGTATGCTAATTACGTAAGAGCGGTAATTATATAAAACGAATTTAAATAAACAGAAAATTATATGAAATGGGGAAGTATCCCACTCGATTATTTATCGGCGACATAATGCAATGATTTATCAATGAATAAAATTGATATTGTTTATTCTATACCGTCATTTTTACCGTCTTTGGAAGAAAATGCCATCAAGAACGTTCAGCAAGCCAAGTATCAATATCAGATTTTCTCCACGCTACACgccc                                                                                                                                       >  1:469383/1‑247 (MQ=255)                                                                                                         gTATGCTAATTACGTAAGAGCGGTAATTATATAAAACGAATTTAAATAAACAGAAAATTATATGAAATGGGGAAGCATCCCACTCGATTATTTATCGGCGACATAATGCAATGATTTATCAATGAATAAAATTGATATTGTTTATTCTATAc                                                                                                                                                                                                                                 >  1:490196/1‑152 (MQ=255)                                                                                                             gCTAATTACGTAAGAGCGGTAATTATATAAAACGAATTTAAATAAACAGAAAATTATATGAAATGGGGAAGTATCCCACTCGATTATTTATCGGCGACATAATGCAATGATTTATCAATGAATAAAATTGATATTGTTTATTCTATACCGTCATTTTTACCGTCTTTGGAAGAAAATGCCATCAAGAACGTTCAGCAAGCCAAGTATCAATATCAGATTTTCTCCACGCTACACGCCCCATGGACGCTTTa                                                                                                                          >  1:267202/1‑251 (MQ=255)                                                                                                                                                   taaataaaCAGAAAATTATATGAAATGGGGAAGTATCCCACTCGATTATTTATCGGCGACATAATGCAATGATTTATCAATGAATAAAATTGATATTGTTTATTCTATACCGTCATTTTTACCGTctttggaagaaaat                                                                                                                                                                                                    >  1:521727/1‑126 (MQ=255)                                                                                                                                                                aaTTATATGAAATGGGGAAGTATCCCACTCGATTATTTATCGGCGACATAATGCAATGATTTATCAATGAATAAAATTGATATTGTTTATTCTATACCGTCATTTTTACCGTCTTTGGAAGAAAATGCCATCAAGAACGTTCAGCAAGCCAAGTATCAATATCAGATTTTCTCCACGCTACACGCCCCATGGACGCTTTAACACTTTTGGGGAATTTACCTTGTCTCATCATACGCCAAATTGTTGTGTCg                                                                       >  1:353600/1‑251 (MQ=255)                                                                                                                                                                 aTTATATGAAATGGGGAAGTATCCCACTCGATTATTTATCGGCGACATAATGCAATGATTTATCAATGAATAAAATTGATATTGTTTATTCTATACCGTCATTTTTACCGTCTTTGGAAGAAAATGCCATCAAGAACGTTCAGCAAGCCAAGTATCAATATCAGATTTtctc                                                                                                                                                     <  1:151176/172‑1 (MQ=255)                                                                                                                                                                         aaaTGGGGAAGTATCCCACTCGATTATTTATCGGCGACATAATGCAATGATTTATCAATGAATAAAATTGATATTGTTTATTCTATACCGTCATTTTTACCGTCTTTGGAAGAAAATGCCATCAAGAACGTTCAGCAAGCCAAGTAt                                                                                                                                                                      >  1:94481/1‑147 (MQ=255)                                                                                                                                                                                gAAGTATCCCACTCGATTATTTATCGGCGACATAATGCAATGATTTATCAATGAATAAAATTGATATTGTTTATTCTATACCGTCATTTTTACCGTCTTTGGAAGAAAATGCCATCAAGAACGTTCAGCAAGCCAAGTATCAATATCAGATTTTCTCCACGCTACACGCCCCATGGACGCTTTAACACTTTTGGGGAATTTACCTTGTCTCATCATACGCCAAATTGTTGTGTCGCTCAAGCTGGTGACtt                                                       >  1:534059/1‑251 (MQ=255)                                                                                                                                                                                 aaGTATCCCACTCGATTATTTATCGGCGACATAATGCAATGATTTATCAATGAATAAAATTGATATTGTTTATTCTATACCGTCATTTTTACCGTCTTTGGAAGAAAATGCCATCAAGAACGTTCAGCAAGCCAAGTATCAATATCAGATTTTCTCCACGCTACACGCCCCATGGACGCTTTAACACTTTTGGGGAATTTACCTTGTCTcatcat                                                                                          >  1:573166/1‑215 (MQ=255)                                                                                                                                                                                             cGATTATTTATCGGCGACATAATGCAATGATTTATCAATGAATAAAATTGATATTGTTTATTCTATACCGTCATTTTTACCGTCTTTGGAAGAAAATGCCATCAAGAACGTTCAGCAAGCCAAGTATCAATATCAGATTTTCTCCACGCTACACGCCCCATGGACGCTTTAACACTTTTGGGGAATTTACCTTGTCTCATCATACGCCAAATTGTTGTGTCGCTCAAGCTGGTGACTTCAATCACCTTTg                                           >  1:811968/1‑250 (MQ=255)                                                                                                                                                                                                  atttatCGGCGACATAATGCAATGATTTATCAATGAATAAAATTGATATTGTTTATTCTATACCGTCATTTTTACCGTCTTTGGAAGAAAATGCCATCAAGAACGTTCAGCAAGCCAAGTATCAATATCAGATTTTCTCCACGCTACACGCCCCATGGACGCTTTAACACTTTTGGGGAATTTACCTTGTCTCATCATACGCCAAATTGTTGTGTCGCTCAAg                                                                 <  1:186895/223‑1 (MQ=255)                                                                                                                                                                                                  atttatCGGCGACATAATGCAATGATTTATCAATGAATAAAATTGATATTGTTTATTCTATACCGTCATTTTTACCGTCTTTGGAAGAAAATGCCATCAAGAACGTTCAGCAAGCCAAGTATCAATATCAGATTTTCTCCACGCTACACGCCCCATGGACGCTTTAACACTTTTGGGGAATTTACCTTGTCTCATCATACGCCAAATTGTTGTGTCGCTCAAGCTGGTGACTTCAATCACCTTTGCTTTa                                      >  1:113098/1‑250 (MQ=255)                                                                                                                                                                                                        cGGCGACATAATGCAATGATTTATCAATGAATAAAATTGATATTGTTTATTCTATACCGTCATTTTTACCGTCTTTGGAAGAAAATGCCATCAAGAACGTTCAGCAAGCCAAGTATCAATATCAGATTTTCTCCACGCTACACGCCCCATGGACGCTTTAACACTTTTGGGGAATTTACCTTGTCTCATCATACGCCAAATTGTTGTGTCGCTCAAGCTGGTGACTTCAATCACCTTTGCTTTACTtagta                               >  1:787709/1‑251 (MQ=255)                                                                                                                                                                                                                                    gAATAAAATTGATATTGTTTATTCTATACCGTCATTTTTACCGTCTTTGGAAGAAAATGCCATCAAGAACGTTCAGCAAGCCAAGTATCAATATCAGATTTTCTCCACGCTACACGCCCCATGGACGCTTTAACACTTTTGGGGAATTTACCTTGtctcatca                                                                                           <  1:558209/163‑8 (MQ=255)                                                                                                                                                                                                                                     aaTAAAATTGATATTGTTTATTCTATACCGTCATTTTTACCGTCTTTGGAAGAAAATGCCATCAAGAACGTTCAGCAAGCCAAGTATCAATATCAGATTTTCTCCACGCTACACGCCCCATGGACGCTTTAACACTTTTGGGGAATTTACCTTGTCTCATCATACGCCAAATTGTTGTGTCGCTCAAGCTGGTGACTTCAATCACCTTTGCTTTACTTAGTAGCGCTGGTTGTTGTTCCATTATTCTtctc  <  1:791741/251‑1 (MQ=255)                                                                                                                                                                                                                                            |                                                                                                                                                                                                                                                    TTAGGTGTTAGAACCTTAGATATCTGATGGCACTGGAGTAAAAACCTGATATAGATCGTGTGTTGAACGATGTTGTACATTTGATGCTAAAGTGCCCTACTTGGGTATGCTAATTACGTAAGAGCGGTAATTATATAAAACGAATTTAAATAAACAGAAAATTATATGAAATGGGGAAGTATCCCACTCGATTATTTATCGGCGATATAATGCAATGATTTATAAATGAATAAAAATGATATTGTTTATTCTCTACCGTCATTTTTACCGTCACCAGAAGAAAATGCCATCAAGAACGTTCAGCAAGCCAAGTATCAATATCAGATTTTCTCCACGCTACACGGCCCATAGACGCTTTAACACTTTTGGGGAATTTACCTTGTTTAATCATACGCCAAATTGTTGTGTCGCTCAAGCTGGTGACTTCAATCACCTTTGCTTTACTTAGTAGCGCTGGTTGTTGTTCCATTATTCTTCTC  >  1/1285060‑1285538

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