Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	1	1,339,007	A→T	P64P (CCT→CCA)	unknown ←	hypothetical protein

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	1	1,339,007	0	A	T	100.0%	92.8 / NA	36	P64P (CCT→CCA)	unknown	hypothetical protein
Reads supporting (aligned to +/- strand):  ref base A (0/0);  new base T (22/14);  total (22/14)

ATGATACACACTTTCTCCTTTAATGATTCACACCTTTTGGCAAAAAAGCGTAAATCATAATATAGGTTCCAACAACAAGAAACTTTACTGAATGCCCATTTTTAAGAGGGAAAATATGAAACTTCCCGTCGGTATAGGCATCAGCAATATAACTAAGTTTTTCTGATAGACACTCTATTTCCTCAATGAATTTGTCAGCTACTTCCTTATTTTGGCTTTGTGTGTAAATATATTCTGCTTTATCTAAAAGCTGAAATCTAGCCATTTCAGACCATTGAACAATGTAGTCTTGTTTCATCGCTTCAATGCCTTGCGCAAATCTTTCATTACATCTTCATGTGAGTAAGTC  >  1/1338901‑1339249                                                                                                           |                                                                                                                                                                                                                                                   aaataacatacaccccttcctAATGATTCACACCTTTTGGCAAAAAAGCGTAAATCATAATATACATTCCTACAACCAGAAACTTCACTGAATGCCCATTTTTAAGTGGAAAAATGTGAAATCTCCCGTCAGTATAGGCGTCAGCAATGTAACTAAGTTTTTCTGATAGACGCTCTATTTCCTCAATGAATTTGTCAGCCACTTCCTTACTTTGGCTTTGTGTGTAAATATATTc                                                                                                                    >  1:818686/21‑235 (MQ=255)      acacacaccccttcctAATGATTCACACCTTTTGGCAAAAAAGCGTAAATCATAATATACATTCCTACAACCAGAAACTTCACTGAATGCCCATTTTTAAGTGGAAAAATGTGAAATCTCCCGTCAGTATAGGCGTCAGCAATGTAACTAAGTTTTTCTGATAGACGCTCTAt                                                                                                                                                                             <  1:448188/158‑1 (MQ=255)                 ttcctAATGATTCACACCTTTTGGCAAAAAAGCGTAAATCATAATATACATTCCTACAACCAGAAACTTCACTGAATGCCCATTTTTAAGTGGAAAAATGTGAAATCTCCCGTCAGTATAGGCGTCAGCAATGTAACTAAGTTTTTCTGATAGACGCTCTATTTCCTCAATGAATTTGTCAGCCACTTCCTTACTTTGGCTTTGTGTGTAAata                                                                                                                         >  1:225443/5‑214 (MQ=255)                  tcctAATGATTCACACCTTTTGGCAAAAAAGCGTAAATCATAATATACATTCCTACAACCAGAAACTTCACTGAATGCCCATTTTTAAGTGGAAAAATGTGAAATCTCCCGTCAGTATAGGCGTCAGCAATGTAACTAAGTTTTTCTGATAGACGCTCTATTTCCTCAATGAATTTGTCAGCCACTTCCTTACTTTGGCTTTGTGTGTAAATATATTCCGCTTTATCTAAAAGCTGAAATCTAGCCAttt                                                                                    >  1:604529/4‑250 (MQ=255)                  tcctAATGATTCACACCTTTTGGCAAAAAAGCGTAAATCATAATATACATTCCTACAACCAGAAACTTCACTGAATGCCCATTTTTAAGTGGAAAAATGTGAAATCTCCCGTCAGTATAGGCGTCAGCAATGTAACTAAGTTTTTCTGATAGACGCTCTATTTCCTCAATGAATTTGTCAGCCACTTCCTTACTTTGGCTTTGTGTGTAAATATATTCCGCTTTATCTAAAAGCTGAAATCTAGCCAttt                                                                                    >  1:545586/4‑250 (MQ=255)                  tcctAATGATTCACACCTTTTGGCAAAAAAGCGTAAATCATAATATACATTCCTACAACCAGAAACTTCACTGAATGCCCATTTTTAAGTGGAAAAATGTGAAATCTCCCGTCAGTATAGGCGTCAGCAATGTAACTAAGATTTTCTGATAGACGCTCTATTTCCTCAATGAATTTGTCAGCCACTTCCTTACTTTGGCTTTGTGTGTaaa                                                                                                                           >  1:79663/4‑211 (MQ=255)                       aTGATTCACACCTTTTGGCAAAAAAGCGTAAATCATAATATACATTCCTACAACCAGAAACTTCACTGAATGCCCATTTTTAAGTGGAAAAATGTGAAATCTCCCGTCAGTATAGGCGTCAGCAATGTAACTAAGTTTTTCTGATAGACGCTCTATTTCCTCAATGAATTTGTCAGCCACTTCCTTACTTTGGCTTTGTGTGTAAATATATTCCGCTTTATCTAAAAGCTGAAATCTAGCCATTTCAGAcc                                                                              <  1:87745/251‑1 (MQ=255)                       aTGATTCACACCTTTTGGCAAAAAAGCGTAAATCATAATATACATTCCTACAACCAGAAACTTCACTGAATGCCCATTTTTAAGTGGAAAAATGTGAAATCTCCCGTCAGTATAGGCGTCAGCAATGTAACTAAGTTTTTCTGATAGACGCTCTATTTCCTCAATGAATTTGTCAGCCACTTCCTTACTTTGGCTTTGTGTGTAAATATATTCCGCTTTATCTAAAAGCTGAAATCTAGCCATTTCAGAcc                                                                              <  1:425820/251‑1 (MQ=255)                         gATTCACACCTTTTGGCAAAAAAGCGTAAATCATAATATACATTCCTACAACCAGAAACTTCACTGAATGCCCATTTTTAAGTGGAAAAATGTGAAATCTCCCGTCAGTATAGGCGTCAGCAATGTAACTAAGTTTTTCTGATAGACGCTCTATTTCCTCAATGAATTTGTCAGCCACTTCCTTACTTTGGCtt                                                                                                                                     >  1:491282/1‑194 (MQ=255)                         gATTCACACCTTTTGGCAAAAAAGCGTAAATCATAATATACATTCCTACAACCAGAAACTTCACTGAATGCCCATTTTTAAGTGGAAAAATGTGAAATCTCCCGTCAGTATAGGCGTCAGCAATGTAACTAAGTTTTTCTGATAGACGCTCTATTTCCTCAATGAATTTGTCAGCCACTTCCTTACTTTGGCTTTGTGTGTAAATATATTCCGCTTTATCTAAAAGCTGAAATCTAGCCATTTCAGACCAt                                                                            >  1:93247/1‑251 (MQ=255)                         gATTCACACCTTTTGGCAAAAAAGCGTAAATCATAATATACATTCCTACAACCAGAAACTTCACTGAATGCCCATTTTTAAGTGGAAAAATGTGAAATCTCCCGTCAGTATAGGCGTCAGCAATGTAACTAAGTTTTTCTGATAGACGCTCTATTTCCTCAATGAATTTGTCAGCCACTTCCTTACTTTGGCTTTGTGTGTAAATATATTCCGCTTTATCTAAAAGCTGAAATCTAGCCATTTCAGACCAt                                                                            >  1:629036/1‑251 (MQ=255)                                 ccTTTTGGCAAAAAAGCGTAAATCATAATATACATTCCTACAACCAGAAACTTCACTGAATGCCCATTTTTAAGTGGAAAAATGTGAAATCTCCCGTCAGTATAGGCGTCAGCAATGTAACTAAGTTTTTCTGATAGACGCTCTa                                                                                                                                                                              >  1:43191/1‑145 (MQ=255)                                      tGGCAAAAAAGCGTAAATCATAATATACATTCCTACAACCAGAAACTTCACTGAATGCCCATTTTTAAGTGGAAAAATGTGAAATCTCCCGTCAGTATAGGC‑GCAGCAATGTAACTAAGTTTTTCTGATAGACGCTCTATTTCCTCAATGAATTTGTCAGCCACTTCCTTACTTTGGCTTTGTGTGTAAATATATTCCGCTTTATCTAAAAGCTGAAATCTAGCCATTTCAGACCATTGAACAACGTaatc                                                              <  1:193324/251‑4 (MQ=255)                                                 cGTAAATCATAATATACATTCCTACAACCAGAAACTTCACTGAATGCCCATTTTTAAGTGGAAAAATGTGAAATCTCCCGTCAGTATAGGCGTCAGCAATGTAACTAAGTTTTTCTGATAGACGCTCTATTTCCTCAATGAATTTGTCAGCCACTTCCTTACTTTGGCTTTGTGTGTAAATATATTCCGCTTTATCTAAAAGCTGAAATCTAGCCATTTCAGACCATTGAACAACGTAATCCTGTTTCAt                                                     <  1:657619/250‑1 (MQ=255)                                                       tCATAATATACATTCCTACAACCAGAAACTTCACTGAATGCCCATTTTTAAGTGGAAAAATGTGAAATCTCCCGTCAGTATAGGCGTCAGCAATGTAACTAAGTTTTTCTGATAGACGCTCTATTTCCTCAATGAATTTGTCAGCCACTTCCTTACTTTGGCTTTGTGTGTAAATATATTCCGCTTTATCTAAAAGCTGAAATCTAGCCATTTCAGACCATTGAACAACGTAATCCTGTTTCATCGTTTCa                                              >  1:400564/1‑251 (MQ=255)                                                       tCATAATATACATTCCTACAACCAGAAACTTCACTGAATGCCCATTTTTAAGTGGAAAAATGTGAAATCTCCCGTCAGTATAGGCGTCAGCAATGTAACTAAGTTTTTCTGATAGACGCTCTATTTCCTCAATGAATTTGTCAGCCACTTCCTTACTTTGGCTTTGTGTGTAAATATATTCCGCTTTATCTAAAAGCTGAAATCTAGCCATTTCAGACCATTGAACAACGTAATCCTGTTTCATCGTTTCa                                              >  1:744739/1‑251 (MQ=255)                                                        cATAATATACATTCCTACAACCAGAAACTTCACTGAATGCCCATTTTTAAGTGGAAAAATGTGAAATCTCCCGTCAGTATAGGCGTCAGCAATGTAACTAAGTTTTTCTGATAGACGCTCTATTTCCTCAATGAATTTGTCAGCCACTTCCTTACTTTGGCTTTGTGTGTAAATATATTCCGCTTTATCTAAAAGCTGAAATCTAGCCATTTc                                                                                   >  1:792671/1‑213 (MQ=255)                                                         ataataTACATTCCTACAACCAGAAACTTCACTGAATGCCCATTTTTAAGTGGAAAAATGTGAAATCTCCCGTCAGTATAGGCGTCAGCAATGTAACTAAGTTTTTCTGATAGACGCTCTATTTCCTCAATGAATTTGTCAGCCACTTCCTTACTTTGGCTTTGTGTGTAAATATATTCCGCTTTATCTAAAAGCTGAAATCTAGCCATTTCAGACCATTGAACAACGTaatcc                                                             >  1:39444/1‑230 (MQ=255)                                                         ataataTACATTCCTACAACCAGAAACTTCACTGAATGCCCATTTTTAAGTGGAAAAATGTGAAATCTCCCGTCAGTATAGGCGTCAGCAATGTAACTAAGTTTTTCTGATAGACGCTCTATTTCCTCAATGAATTTGTCAGCCACTTCCTTACTTTGGCTTTGTGTGTAAATATATTCCGCTTTATCTAAAAGCTGAAATCTAGCCATTTCAGACCATTGAACAACGTAATCCTGTTTCATCGTTTCAAt                                            >  1:491142/1‑251 (MQ=255)                                                          taataTACATTCCTACAACCAGAAACTTCACTGAATGCCCATTTTTAAGTGGAAAAATGTGAAATCTCCCGTCAGTATAGGCGTCAGCAATGTAACTAAGTTTTTCTGATAGACGCTCTATTTCCTCAATGAATTTGTCAGCCACTTCCTTACTTTGGCTTTGTGTGTAAATATATTccgc                                                                                                                 <  1:269537/181‑4 (MQ=255)                                                          taataTACATTCCTACAACCAGAAACTTCACTGAATGCCCATTTTTAAGTGGAAAAATGTGAAATCTCCCGTCAGTATAGGCGTCAGCAATGTAACTAAGTTTTTCTGATAGACGCTCTATTTCCTCAATGAATTTGTCAGCCACTTCCTTACTTTGGCTTTGTGTGTAAATATATTCCGCTTTATCTAAAAGCTGAAATCTAGCCATTTCAGACCATTGAACaacgt                                                                  >  1:592328/1‑225 (MQ=255)                                                            atatacattCCTACAACCAGAAACTTCACTGAATGCCCATTTTTAAGTGGAAAAATGTGAAATCTCCCGTCAGTATAGGCGTCAGCAATGTAACTAAGTTTTTCTGATAGACGCTCTATTTCCTCAATGAATTTGTCAGCCACTTCCTTACTTTGGCTTTGTGTGTAAATATATTCCGCTTTATCTAAAAGCTGAAATCTAGCCATTTCAGACCATTGAACAACGTAATCCTGTTTCATCGTTTCAATGc                                          >  1:623153/8‑250 (MQ=255)                                                             tatacattCCTACAACCAGAAACTTCACTGAATGCCCATTTTTAAGTGGAAAAATGTGAAATCTCCCGTCAGTATAGGCGTCAGCAATGTAACTAAGTTTTTCTGATAGACGCTCTATTTCCTCAATGAATTTGTCAGCCACTTCCTTACTTTGGCTTTGTGTGTAAATATATTCCGCTTTATCTAAAAGCTGAAATCTAGCCATTTCAGACCATTGAACAACGTAATCCTGTTTCATCGTTTCAATGCCt                                        <  1:665540/245‑1 (MQ=255)                                                                     cctacaaccaGAAACTTCACTGAATGCCCATTTTTAAGTGGAAAAATGTGAAATCTCCCGTCAGTATAGGCGTCAGCAATGTAACTAAGTTTTTCTGATAGACGCTCTATTTCCTCAATGAATTTGTCAGCCACTTCCTTACTTTGGCTTtgtg                                                                                                                                 >  1:767337/4‑154 (MQ=255)                                                                     cctacaaccaGAAACTTCACTGAATGCCCATTTTTAAGTGGAAAAATGTGAAATCTCCCGTCAGTATAGGCGTCAGCAATGTAACTAAGTTTTTCTGATAGACGCTCTATTTCCTCAATGAATTTGTCAGCCACTTCCTTACTTTGGCTTTGTGTGTAAata                                                                                                                         <  1:488859/159‑1 (MQ=255)                                                                                  aCTTCACTGAATGCCCATTTTTAAGTGGAAAAATGTGAAATCTCCCGTCAGTATAGGCGTCAGCAATGTAACTAAGTTTTTCTGATAGACGCTCTATTTCCTCAATGAATTTGTCAGCCACTTCCTTACTTTGGCTTTGTGTGTAAATATATTCCGCTTTATCTAAAAGCTGAAATCTAGCCATTTCAGAc                                                                               >  1:518623/1‑191 (MQ=255)                                                                                   cttcaCTGAATGCCCATTTTTAAGTGGAAAAATGTGAAATCTCCCGTCAGTATAGGCGTCAGCAATGTAACTAAGTTTTTCTGATAGACGCTCTATTTCCTCAATGAATTTGTCAGCCACTTCCTTACTTTGGCTTTGTGTGTAAATATATTCCGCTTTATCTAAAAGCTGAAATCTAGCCATTTCAGACCATTGAACAACGTAATCCTGTTTCATCGTTTCAATGCCTTGCGCAAATCTTTCATTAcct                   >  1:857298/5‑248 (MQ=255)                                                                                         tGAATGCCCATTTTTAAGTGGAAAAATGTGAAATCTCCCGTCAGTATAGGCGTCAGCAATGTAACTAAGTTTTTCTGATAGACGCTCTATTTCCTCAATGAATTTGTCAGCCACTTCCTTACTTTGGCTTTGTGTGTAAATATATTCCGCTTTATCTAAAAGCTGAAATCTAGCCATTTCAGACCATTGAACAACGTAATCCTGTTTCATCGTTTCAATGCCTTGCGCAAATCTTTCATTACCTCTTCAtg            <  1:421895/251‑1 (MQ=255)                                                                                         tGAATGCCCATTTTTAAGTGGAAAAATGTGAAATCTCCCGTCAGTATAGGCGTCAGCAATGTAACTAAGTTTTTCTGATAGACGCTCTATTTCCTCAATGAATTTGTCAGCCACTTCCTTACTTTGGCTTTGTGTGTAAATATATTCCGCTTTATCTAAAAGCTGAAATCTAGCCATTTCAGACCATTGAACAACGTAATCCTGTTTCATCGTTTCAATGCCTTGCGCAAATCTTTCATTACCTCTTCAtg            <  1:410129/251‑1 (MQ=255)                                                                                           aaTGCCCATTTTTAAGTGGAAAAATGTGAAATCTCCCGTCAGTATAGGCGTCAGCAATGTAACTAAGTTTTTCTGATAGACGCTCTATTTCCTCAATGAATTTGTCAGCCACTTCCTTACTTTGGCTTTGTGTGTAAATATATTCCGCTTTATCTAAAAGCTGAAATCTAGCCATTTCAGACCATTGAACAACGTAATCCTGTTTCATCGTTTCAATGCCTTGCGCAAATCTTTCATTACCTCTTCAtgtg          <  1:277073/251‑1 (MQ=255)                                                                                           aaTGCCCATTTTTAAGTGGAAAAATGTGAAATCTCCCGTCAGTATAGGCGTCAGCAATGTAACTAAGTTTTTCTGATAGACGCTCTATTTCCTCAATGAATTTGTCAGCCACTTCCTTACTTTGGCTTTGTGTGTAAATATATTCCGCTTTATCTAAAAGCTGAAATCTAGCCATTTCAGACCATTGAACAACGTAATCCTGTTTCATCGTTTCAATGCCTTGCGCAAATCTTTCATTACCTCTTCAtgtg          <  1:228148/251‑1 (MQ=255)                                                                                                 cATTTTTAAGTGGAAAAATGTGAAATCTCCCGTCAGTATAGGCGTCAGCAATGTAACTAAGTTTTTCTGATAGACGCTCTATTTCCTCAATGAATTTGTCAGCCACTTCCTTACTTTGGCTTTGTGTGTAAATATATTCCGCTTTATCTAAAAGCTGAAATCTAGCCATTTCAGACCATTGAACAACGTAATCCTGTTTCATCGTTTCAATGCCTTGCGCAAATCTTTCAt                        >  1:552000/1‑231 (MQ=255)                                                                                                  aTTTTTAAGTGGAAAAATGTGAAATCTCCCGTCAGTATAGGCGTCAGCAATCTAACTAAGTTTTTCTGATAGACGCTCTATTTCCTCAATGAATTTGTCAGCCACTTCCTTACTTTGGCTTTGTGTGTAAATATATTCCGCTTTATCTAAAAGCTGAAATCTAGCCATTTCAGACCATTGAACAACGTAATCCTGTTTCATCGTTTCAATGCCTTGCGCAAATCTTTCATTACCTCTTCATGTGAGTaggt   >  1:450901/1‑248 (MQ=255)                                                                                                   tttttAAGTGGAAAAATGTGAAATCTCCCGTCAGTATAGGCGTCAGCAATGTAACTAAGTTTTTCTGATAGACGCTCTATTTCCTCAATGAATTTGTCAGCCACTTCCTTac                                                                                                                                             >  1:144322/1‑111 (MQ=255)                                                                                                   tttttAAGTGGAAAAATGTGAAATCTCCCGTCAGTATAGGCGTCAGCAATGTAACTAAGTTTTTCTGATAGACGCTCTATTTCCTCAATGAATTTGTCAGCCACTTCCTTACTTTGGCTTTGTGTGTAAATATATTCCGCTTTATCTAAAAGCTGAAATCTAGCCATTTCAGACCATTGAACAACGTAATCCTGTTTCATCGTTTCAATGCCTTGCGCAAATCTTTCATTACCTCTTCATGTGAGTAGGTc  <  1:669865/251‑1 (MQ=255)                                                                                                   tttttAAGTGGAAAAATGTGAAATCTCCCGTCAGTATAGGCGTCAGCAATGTAACTAAGTTTTTCTGATAGACGCTCTATTTCCTCAATGAATTTGTCAGCCACTTCCTTACTTTGGCTTTGTGTGTAAATATATTCCGCTTTATCTAAAAGCTGAAATCTAGCCATTTCAGACCATTGAACAACGTAATCCTGTTTCATCGTTTCAATGCCTTGCGCAAATCTTTCATTACCTCTTCATGTGAGTAGGTc  <  1:188738/251‑1 (MQ=255)                                                                                                           |                                                                                                                                                                                                                                                   ATGATACACACTTTCTCCTTTAATGATTCACACCTTTTGGCAAAAAAGCGTAAATCATAATATAGGTTCCAACAACAAGAAACTTTACTGAATGCCCATTTTTAAGAGGGAAAATATGAAACTTCCCGTCGGTATAGGCATCAGCAATATAACTAAGTTTTTCTGATAGACACTCTATTTCCTCAATGAATTTGTCAGCTACTTCCTTATTTTGGCTTTGTGTGTAAATATATTCTGCTTTATCTAAAAGCTGAAATCTAGCCATTTCAGACCATTGAACAATGTAGTCTTGTTTCATCGCTTCAATGCCTTGCGCAAATCTTTCATTACATCTTCATGTGAGTAAGTC  >  1/1338901‑1339249

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