Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	1	1,568,091	C→T	Q120Q (CAG→CAA)	unknown ←	hypothetical protein

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	1	1,568,091	0	C	T	100.0%	61.1 / NA	23	Q120Q (CAG→CAA)	unknown	hypothetical protein
Reads supporting (aligned to +/- strand):  ref base C (0/0);  new base T (12/11);  total (12/11)

ACTCATTATTGTGCACAATATTCAACTGATGTTTATAGCAAATATTGGTAATTTTAGAAATGAGCCCTTTATCGTCAGGGCAATCGGTCAGTAAGATTTTATTTTGAATCATAATGCGTAACATTAAGATAAAAAGAAAACCCATAAAAATGGGTTTTCTTTTATTGGGTTTTCTTTTATTAAGATAGAAATTGAATTAAATTTTAAAAGAAGATAAAGACTTACCTTTATCTAACTGTGCTTGGATTGCTTTCGGTGTTCTGCCTTGACCTGTCCAAGTTTTCGTTACGCCCTTTTCATCAACATATTGATATTTTGCAGGACGAGTATCTCTTTTCTTACGAACCGCATCAGAACCTAAAATCTCTGCTAATTCGTCTGCAGTGATACCTTCTTGCTTTAATAATTCTTTATATTTGGCAATACGTTCTTGTCGCTCGTTTTCTTCTTGCTGAGCCTTTAAAATCTCTTCTTTCTTTT  >  1/1567856‑1568335                                                                                                                                                                                                                                            |                                                                                                                                                                                                                                                     aCTCATTATTGTGCACAATATTCAACTGATGTTTATAGCAAATATTGGTAATTTTAGAAATGAGCCCTTTATCGTCAGGGCAATCGGTCAGTAAGATTTTATTTTGAATCATAATGCGTAACATTAAGATAAAAAGAAAACCCATAAAAATGGGTTTTCTTTTATTGGGTTTTCTTTTATTAAGATAGAAATTGAATTAAATTTTAAAAGAAGATAAAGACTTACCTTTATCTAATTGTGCTTGGATTg                                                                                                                                                                                                                                         >  1:371747/1‑249 (MQ=255)    cATTATTGTGCACAATATTCAACTGATGTTTATAGCAAATATTGGTAATTTTAGAAATGAGCCCTTTATCGTCAGGGCAATCGGTCAGTAAGATTTTATTTTGAATCATAATGCGTAACATTAAGATAAAAAGAAAACCCATAAAAATGGGTTTTCTTTTATTGGGTTTTCTTTTATTAAGATAGAAATTGAATTAAATTTTAAAAGAAGATAAAGACTTACCTTTATCTAATTGTGCTTGGATTGCttt                                                                                                                                                                                                                                     >  1:234212/1‑250 (MQ=255)                              tGTTTATAGCAAATATTGGTAATTTTAGAAATGAGCCCTTTATCGTCAGGGCAATCGGTCAGTAAGATTTTATTTTGAATCATAATGCGTAACATTAAGATAAAAAGAAAACCCATAAAAATGGGTTTTCTTTTATTGGGTTTTCTTTTATTAAGATAGAAATTGAATTAAATTTTAAAAGAAGATAAAGACTTACCTTTATCTAATTGTGCTTGGATTGCTTTCGGTGTTCTGCCTTGACCTGTCCAAGt                                                                                                                                                                                                          <  1:239571/251‑1 (MQ=255)                                                             tGAGCCCTTTATCGTCAGGGCAATCGGTCAGTAAGATTTTATTTTGAATCATAATGCGTAACATTAAGATAAAAAGAAAACCCATAAAAATGGGTTTTCTTTTATTGGGTTTTCTTTTATTAAGATAGAAATTGAATTAAATTTTAAAAGAAGATAAAGACTTACCTTTATCTAATTGTGCTTGGATTGCTTTCGGTGTTCTGCCTTGACCTGTCCAAGTTTTCGTTACGCCCTTTTCATCAACATATTGa                                                                                                                                                                           <  1:483079/251‑1 (MQ=255)                                                                                   aTCGGTCAGTAAGATTTTATTTTGAATCATAATGCGTAACATTAAGATAAAAAGAAAACCCATAAAAATGGGTTTTCTTTTATTGGGTTTTCTTTTATTAAGATAGAAATTGAATTAAATTTTAAAAGAAGATAAAGACTTACCTTTATCTAATTGTGCTTGGATTGCTTTCGGTGTTCTGCCTTGACCTGTCCAAGTTTTCGTTACGCCCTTTTCATCAAc                                                                                                                                                                                  <  1:299867/222‑1 (MQ=255)                                                                                       gTCAGTAAGATTTTATTTTGAATCATAATGCGTAACATTAAGATAAAAAGAAAACCCATAAAAATGGGTTTTCTTTTATTGGGTTTTCTTTTATTAAGATAGAAATTGAATTAAATTTTAAAAGAAGATAAAGACTTACCTTTATCTAATTGTGCTTGGATTGCTTTCGGTGTTCTGCCTTGACCTGTCCAAGTTTTCGTTACGCCCTTTTCATCAACATATTGATATTTTGCAGGACGAGTATCTCtttt                                                                                                                                                 >  1:276279/1‑251 (MQ=255)                                                                                                                     cGTAACATTAAGATAAAAAGAAAACCCATAAAAATGGGTTTTCTTTTATTGGGTTTTCTTTTATTAAGATAGAAATTGAATTAAATTTTAAAAGAAGATAAAGACTTACCTTTATCTAATTGTGCTTGGATTGCTTTCGGTGTTCTGCCTTGACCTGTCCAAGTTTTCGTTACGCCCTTTTCATCAACATATTGATATTTTGCAGGACGAGTATCTCTTTTCTTACGAACCGCATCAGAACCTAAAAtct                                                                                                                    >  1:510711/1‑250 (MQ=255)                                                                                                                                         aaaaCCCATAAAAATGGGTTTTCTTTTATTGGGTTTTCTTTTATTAAGATAGAAATTGAATTAAATTTTAAAAGAAGATAAAGACTTACCTTTATCTAATTGTGCTTGGATTGCTTTCGGTGTTCTGCCTTGACCTGTCCAAGTTTTCGTTACGCCCTTTTCATCAACATATTGATATTTTGCAGGACGAGTATCTCTTTTCTTACGAACCGCATCAGAACCTAAAATCTCTGCTAATTCGTCTGCAGt                                                                                                 <  1:799879/249‑1 (MQ=255)                                                                                                                                                                  ttattgggttttcttttattAAGATAGAAATTGAATTAAATTTTAAAAGAAGATAAAGACTTACCTTTATCTAATTGTGCt                                                                                                                                                                                                                                                <  1:461863/81‑1 (MQ=255)                                                                                                                                                                      tgggttttcttttattAAGATAGAAATTGAATTAAATTTTAAAAGAAGATAAAGACTTACCTTTATCTAATTGTGCTTGGATTGCTTTCGGTGTTCTGCCTTGACCTGTCCAAGTTTTCGTTACGCCCTTTTCATCAACATATTGATATTTTGCAGGACGAGTATCTCTTTTCTTACGAACCGCATCAGAACCTAAAATCTCTGCTAATTCGTCTGCAGTGATACCTTCTTGCTTTAATAATTCTTTatat                                                                  <  1:337309/251‑1 (MQ=255)                                                                                                                                                                                                                     gaTAAAGATTTACCTTTATCTAATTGTGCTTGGATTGCTTTCGGTGTTCTGCCTTGACCTGTCCAAGTTTTCGTTACGCCCTTTTCATCAACATATTGATATTTTGCAGGACGAGTATCTCTTTTCTTACGAACCGCATCAGAACCTAAAATCTCTGCTAATTCGTCTGCAGTGATAc                                                                                            <  1:439155/178‑1 (MQ=255)                                                                                                                                                                                                                     gaTAAAGACTTACCTTTATCTAATTGTGCTTGGATTGCTTTCGGTGTTCTGCCTTGACCTGTCCAAGTTTTCGTTACGCCCTTTTCATCAACATATTGATATTTTGCAGGACGAGTATCTCTTTTCTTACGAACCGCATCAGAACCTAAAATCTCTGCTAATTCGTCTGCAGTGATACCTTCTTGCTTTAATAATTCTTTATATTTGGCAATACGTTCTTTTCGCTCGTTTTCTTCTTGCTGAGCCTTTaa                   >  1:808591/1‑251 (MQ=255)                                                                                                                                                                                                                              ttACCTTTATCTAATTGTGCTTGGATTGCTTTCGGTGTTCTGCCTTGACCTGTCCAAGTTTTCGTTACGCCCTTTTCATCAACATATTGATATTTTGCAGGACGAGTATCTCTTTTCTTACGAACCGCATCAGAACCTAAAATCTCTGCTAATTCGTCTGCAGTGATACCTTCTTGCTTTAATAATTCTTTATATTTGGCAATACGTTCTTTTCGCTCGTTTTCTTCTTGCTGAGCCTTTAAAATCtcttc          <  1:362263/251‑1 (MQ=255)                                                                                                                                                                                                                                aCCTTTATCTAATTGTGCTTGGATTGCTTTCGGTGTTCTGCCTTGACCTGTCCAAGTTTTCGTTACGCCCTTTTCATCAACATATTGATATTTTGCAGGACGAGTATCTCTTTTCTTACGAACCGCATCAGAACCTAAAATCTCTGCTAATTCGTCTGCAGTGATACCTTCTTGCTTTAATAATTCTTTATATTTGGCAATACGTTCTTTTCGCTCGTTTTCTTCTTGCTGAGCCTTTAAAATCTcttctt        >  1:526310/1‑251 (MQ=255)                                                                                                                                                                                                                                 ccTTTATCTAATTGTGCTTGGATTGCTTTCGGTGTTCTGCCTTGAcc                                                                                                                                                                                                                   <  1:790233/47‑1 (MQ=255)                                                                                                                                                                                                                                 ccTTTATCTAATTGTGCTTGGATTGCTTTCGGTGTTCTGCCTTGACCTGTCCAAGTTTTCGTTACGCCCTTTTCATCAACATATTGATATTTTGCAGGACGAGTATCTCTTTTCTTACGAACCGCATCAGAACCTAAAATCTCTGCTAATTCGTCTGCAGTGATACCTTCTTGCt                                                                                   >  1:119657/1‑175 (MQ=255)                                                                                                                                                                                                                                 ccTTTATCTAATTGTGCTTGGATTGCTTTCGGTGTTCTGCCTTGACCTGTCCAAGTTTTCGTTACGCCCTTTTCATCAACATATTGATATTTTGCAGGACGAGTATCTCTTTTCTTACGAACCGCATCAGAACCTAAAATCTCTGCTAATTCGTCTGCAGTGATACCTTCTTGCTTTAATAATTCTTTATATTTGGCAATACGTTCTTTTCGCTCGTTTTCTTCTTGCTGAGCCTTTAAAATCTCttcttt       >  1:250281/1‑251 (MQ=255)                                                                                                                                                                                                                                      aTCTAATTGTGCTTGGATTGCTTTCGGTGTTCTGCCTTGACCTGTCCAAGTTTTCGTTACGCCCTTTTCATCaa                                                                                                                                                                                   <  1:897651/74‑1 (MQ=255)                                                                                                                                                                                                                                      aTCTAATTGTGCTTGGATTGCTTTCGGTGTTCTGCCTTGACCTGTCCAAGTTTTCGTTACGCCCTTTTCATCAACATATTGATATTTTGCAGGACGAGTATCTCTTTTCTTACGAACCGCATCAGAACCTAAAAtc                                                                                                                     >  1:728790/1‑136 (MQ=255)                                                                                                                                                                                                                                      aTCTAATTGTGCTTGGATTGCTTTCGGTGTTCTGCCTTGACCTGTCCAAGTTTTCGTTACGCCCTTTTCATCAACATATTGATATTTTGCAGGACGAGTATCTCTTTTCTTACGAACCGCATCAGAACCTAAAAtc                                                                                                                     <  1:9452/136‑1 (MQ=255)                                                                                                                                                                                                                                      aTCTAATTGTGCTTGGATTGCTTTCGGTGTTCTGCCTTGACCTGTCCAAGTTTTCGTTACGCCCTTTTCATCAACATATTGATATTTTGCAGGACGAGTATCTCTTTTCTTACGAACCGCATCAGAACCTAAAATCTCTGCTAATTCGTCTGCAGTGATACCTTCTTGCTTTAATAATTCTTTATATTTGGCAATACGTTCTTTTCGCTCGTTTTCTTCTTGCTGAGCCTTTAAAATCTCTTCTTTCtttt  >  1:325581/1‑251 (MQ=255)                                                                                                                                                                                                                                      aTCTAATTGTGCTTGGATTGCTTTCGGTGTTCTGCCTTGACCTGTCCAAGTTTTCGTTACGCCCTTTTCATCAACATATTGATATTTTGCAGGACGAGTATCTCTTTTCTTACGAACCGCATCAGAACCTAAAATCTCTGCTAATTCGTCTGCAGTGATACCTTCTTGCTTTAATAATTCTTTATATTTGGCAATACGTTCTTTTCGCTCGTTTTCTTCTTGCTGAGCCTTTAAAATCTCTTCTTTCtttt  >  1:112820/1‑251 (MQ=255)                                                                                                                                                                                                                                      aTCTAATTGTGCTTGGATTGCTTTCGGTGTTCTGCCTTGACCTGTCCAAGTTTTCGTTACGCCCTTTTCATCAACATATTGATATTTTGCAGGACGAGTATCTCTTTTCTTACGAACCGCATCAGAACCTAAAATCTCTGCTAATTCGTCTGCAGTGATACCTTCTTGCTTTAATAATTCTTTATATTTGGCAATACGTTCTTTTCGCTCGTTTTCTTCTTGCTGAGCCTTTAAAATCTCTTCTTTCtttt  >  1:864991/1‑251 (MQ=255)                                                                                                                                                                                                                                            |                                                                                                                                                                                                                                                     ACTCATTATTGTGCACAATATTCAACTGATGTTTATAGCAAATATTGGTAATTTTAGAAATGAGCCCTTTATCGTCAGGGCAATCGGTCAGTAAGATTTTATTTTGAATCATAATGCGTAACATTAAGATAAAAAGAAAACCCATAAAAATGGGTTTTCTTTTATTGGGTTTTCTTTTATTAAGATAGAAATTGAATTAAATTTTAAAAGAAGATAAAGACTTACCTTTATCTAACTGTGCTTGGATTGCTTTCGGTGTTCTGCCTTGACCTGTCCAAGTTTTCGTTACGCCCTTTTCATCAACATATTGATATTTTGCAGGACGAGTATCTCTTTTCTTACGAACCGCATCAGAACCTAAAATCTCTGCTAATTCGTCTGCAGTGATACCTTCTTGCTTTAATAATTCTTTATATTTGGCAATACGTTCTTGTCGCTCGTTTTCTTCTTGCTGAGCCTTTAAAATCTCTTCTTTCTTTT  >  1/1567856‑1568335

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