Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	1	1,518,947	G→T	intergenic (‑226/‑30)	unknown ← / → unknown	Pantothenate kinase (EC 2.7.1.33)/tRNA‑Thr‑TGT

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	1	1,518,947	0	G	T	100.0%	90.7 / NA	35	intergenic (‑226/‑30)	unknown/unknown	Pantothenate kinase (EC 2.7.1.33)/tRNA‑Thr‑TGT
Reads supporting (aligned to +/- strand):  ref base G (0/0);  new base T (16/19);  total (16/19)

GCATAAGCTGCTCCGAAATAGAATGAGAAACCGGAAATTCCACGTGTTGAATTACTCTTTTTAGCCAAATATAGGTTGCAAAATATACCTTATAAATAAGAGATTGAGCAGTAAAAATTTAGCATAATCAGGTAAATTAGACGAAGATAAATTGTTTTTGTTTGATAATTAAACGAATAATTGAAAAAAACTTTATTTACTAAAAAAAAACTTGTCAGGCATGAAAAAAGCCCTATAATACGCGTCACTTGTTTTCGACGCCGACTTAGCTCAGTAGGTAGAGCAACTGACTTGTAATCAGTAGGTCATCAGTTCGATTCCGATAGTCGGCACCATTTCTCGCAATACAAGTTCAAACTTAGAGTGGAGGGATTCCCGAGCGGCCAAAGGGAGCAGACTGTAAATCTGCCGGCTCAGCCTTCGAAGGTTCGAATCCTTCTCCCTCCACCATTATCTAAGAGTAACCGATGGGA  >  1/1518718‑1519190                                                                                                                                                                                                                                      |                                                                                                                                                                                                                                                    gCATAAGCTGCTCCGAAATAGAATGAGAAACCGGAAATTCCACGTGTTGAATTACTCTTTTTAGCCAAATATAGGTTGCAAAATATACCTTATAAATAAGAGATTGAGCAGTAAAAATTTAGCATAATCAGGTAAATTAGACGAAGATAAATTGTTTTTGTTTGATAATTAAACGAATAATTGAAAAAAACTTTATTTACTAAAAAAAAACTTGTCAGGCATGAAAAAATCCCTATAATACGCGTCACTTg                                                                                                                                                                                                                                >  1:292041/1‑251 (MQ=255)                      aaTGAGAAACCGGAAATTCCACGTGTTGAATTACTCTTTTTAGCCAAATATAGGTTGCAAAATATACCTTATAAATAAGAGATTGAGCAGTAAAAATTTAGCATAATCAGGTAAATTAGACGAAGATAAATTGTTTTTGTTTGATAATTAAACGAATAATTGAAAAAAACTTTATTTACTAAAAAAAAACTTGTCAGGCATGAAAAAATCCCTATAATACGCGTCACTTGTTTTCGACGCCGACTTAGCTc                                                                                                                                                                                                           <  1:42221/251‑1 (MQ=255)                                            gtgtTGAATTACTCTTTTTAGCCAAATATAGGTTGCAAAATATACCTTATAAATAAGAGATTGAGCAGTAAAAATTTAGCATAATCAGGTAAATTAGACGAAGATAAATTGTTTTTGTTTGATAATTAAACGAATAATTGAAAAAAACTTTATTTACTAAAAAAAAACTTGTCAGGCATGAAAAAATCCCTATAATACGCGTCACTTg                                                                                                                                                                                                                                >  1:102398/1‑208 (MQ=255)                                                  aaTTACTCTTTTTAGCCAAATATAGGTTGCAAAATATACCTTATAAATAAGAGATTGAGCAGTAAAAATTTAGCATAATCAGGTAAATTAGACGAAGATAAATTGTTTTTGTTTGATAATTAAACGAATAATTGAAAAAAACTTTATTTACTAAAAAAAAACTTGTCAGGCATGAAAAAATCCCTATAATACGCGTCACTTGTTTTCGACGCCGACTTAGCTCAGTAGGTAGAGCAACTGACTTGTAATc                                                                                                                                                                                <  1:437396/250‑1 (MQ=255)                                                  aaTTACTCTTTTTAGCCAAATATAGGTTGCAAAATATACCTTATAAATAAGAGATTGAGCAGTAAAAATTTAGCATAATCAGGTAAATTAGACGAAGATAAATTGTTTTTGTTTGATAATTAAACGAATAATTGAAAAAAACTTTATTTACTAAAAAAAAACTTGTCAGGCATGAAAAAATCCCTATAATACGCGTCACTTGTTTTCGACGCCGACTTAGCTCAGTAGGTAGAGCAACTGACTTGTAATc                                                                                                                                                                                <  1:84399/250‑1 (MQ=255)                                                                 ccAAATATAGGTTGCAAAATATACCTTATAAATAAGAGATTGAGCAGTAAAAATTTAGCATAATCAGGTAAATTAGACGAAGATAAATTGTTTTTGTTTGATAATTAAACGAATAATTGAAAAAAACTTTATTTACTAAAAAAAAACTTGTCAGGCATGAAAAAATCCCTATAATAc                                                                                                                                                                                                                                          <  1:724437/177‑1 (MQ=255)                                                                   aaaTATAGGTTGCAAAATATACCTTATAAATAAGAGATTGAGCAGTAAAAATTTAGCATAATCAGGTAAATTAGACGAAGATAAATTGTTTTTGTTTGATAATTAAACGAATAATTGAAAAAAACTTTATTTACTAAAAAAAAACTTGTCAGGCATGAAAAAATCCCTATAATACGCGTCACTTGTTTTCGACGCCGACTTAGCTCAGTAGGTAGAGCAACTGACTTGTAATCAGTAGGTCATCAGTTCGa                                                                                                                                                              >  1:744227/1‑251 (MQ=255)                                                                                         cTTATAAATAAGAGATTGAGCAGTAAAAATTTAGCATAATCAGGTAAATTAGACGAAGATAAATTGTTTTTGTTTGATAATTAAACGAATAATTGAAAAAAACTTTATTTACTAAAAAAAAACTTGTCAGGCATGAAAAAATCCCTATAATACGCGTCACTTGTTTTCGACGCCGACTTAGCTCAgtaggtag                                                                                                                                                                                                  <  1:164265/193‑1 (MQ=255)                                                                                            ataaataaGAGATTGAGCAGTAAAAATTTAGCATAATCAGGTAAATTAGACGAAGATAAATTGTTTTTGTTTGATAATTAAACGAATAATTGAAAAAAACTTTATTTACTAAAAAAAAACTTGTCAGGCATGAAAAAATCCCTATAATAc                                                                                                                                                                                                                                          <  1:867207/150‑1 (MQ=255)                                                                                            ataaataaGAGATTGAGCAGTAAAAATTTAGCATAATCAGGTAAATTAGACGAAGATAAATTGTTTTTGTTTGATAATTAAACGAATAATTGAAAAAAACTTTATTTACTAAAAAAAAACTTGTCAGGCATGAAAAAATCCCTATAATACGCGTCACTTGTTTTCGACGCCGACTTAGCTCAGTAGGTAGAGCAACTGACTTGTAATCAGTAGGTCATCAGTTCGATTCCGATAGTCGGCACCAttt                                                                                                                                         >  1:298869/1‑247 (MQ=255)                                                                                              aaataaGAGATTGAGCAGTAAAAATTTAGCATAATCAGGTAAATTAGACGAAGATAAATTGTTTTTGTTTGATAATTAAACGAATAATTGAAAAAAACTTTATTTACTAAAAAAAAACTTGTCAGGCATGAAAAAATCCCTATAATACGCGTCACTTGTTTTCGACGCCGACTTAGCTCAGTAGGTAGAGCAACTGACTTGTAATCAGTAGGTCATCAGTTCGATTCCGATAGTCGGCACCATTTCTCGCa                                                                                                                                   >  1:163249/1‑251 (MQ=255)                                                                                                ataaGAGATTGAGCAGTAAAAATTTAGCATAATCAGGTAAATTAGACGAAGATAAATTGTTTTTGTTTGATAATTAAACGAATAATTGAAAAAAACTTTATTTACTAAAAAAAAACTTGTCAGGCATGAAAAAATCCCTATAATACGCGTCACTTg                                                                                                                                                                                                                                <  1:208797/156‑1 (MQ=255)                                                                                                ataaGAGATTGAGCAGTAAAAATTTAGCATAATCAGGTAAATTAGACGAAGATAAATTGTTTTTGTTTGATAATTAAACGAATAATTGAAAAAAACTTTATTTACTAAAAAAAAACTTGTCAGGCATGAAAAAATCCCTATAATACGCGTCACTTGTTTTCGACGCCGACTTAGCTCAGTAGGTAGAGCAACTGACTTGTAATCAGTAGGTCATCAGTTCGATTCCGATAGTCGGCACCATTTCTCGCAAt                                                                                                                                 >  1:434779/1‑251 (MQ=255)                                                                                                     agaTTGAGCAGTAAAAATTTAGCATAATCAGGTAAATTAGACGAAGATAAATTGTTTTTGTTTGATAATTAAACGAATAATTGAAAAAAACTTTATTTACTAAAAAAAAACTTGTCAGGCATGAAAAAATCCCTATAATACGCGTCACTTGTTTTCGACGCCGACTTAGCTCAGTAGGTAGAGCAACTGACTTGTAATCAGTAGGTCATCAGTTCGATTCCGATAGTCGGCACCATTTCTCGCAATACAAg                                                                                                                            >  1:652962/1‑251 (MQ=255)                                                                                                                                tCAGGTAAATTAGACGAAGATAAATTGTTTTTGTTTGATAATTAAACGAATAATTGAAAAAAACTTTATTTACTAAAAAAAAACTTGTCAGGCATGAAAAAATCCCTATAATACGCGTCACTTGTTTTCGACGCCGACTTAGCTCAGTAGGTAGAGCAACTGACTTGTAATCAGTAGGTCATCAGTTCGATTCCGATAGTCGGCACCATTTCTCGCAATACAAGTTCAAACTTAGAGTGGAGGGATTCCCg                                                                                                 <  1:741137/251‑1 (MQ=255)                                                                                                                                            gACGAAGATAAATTGTTTTTGTTTGATAATTAAACGAATAATTGAAAAAAACTTTATTTACTAAAAAAAAACTTGTCAGGCATGAAAAAATCCCTATAATACGCGTCACTTGTTTTCGACGCCGACTTAGCTCAGTAGGTAGAGCAACTGACTTGTAATCAGTAGGTCATCAGTTCGATTCCGATAGTCGGCACCATTTCTCGCAATACAAGTTCAAACTTAGAGTGGAGGGATTCCCGAGCGGCCAAAgg                                                                                     >  1:210745/1‑251 (MQ=255)                                                                                                                                             aCGAAGATAAATTGTTTTTGTTTGATAATTAAACGAATAATTGAAAAAAACTTTATTTACTAAAAAAAAACTTGTCAGGCATGAAAAAATCCCTATAATACGCGTCACTTGTTTTCGACGCCGACTTAGCTCAGTAGGTAGAGCAACTGACTTGTAATCAGTAGGTCATCAGTTCGATTCCGATAGTCGGCACCATTTCTCGCAATACAAGTTCAAACTTAGAGTGGAGGGATTCCCGAGCGGCCAAAggg                                                                                    <  1:838737/251‑1 (MQ=255)                                                                                                                                                aaGATAAATTGTTTTTGTTTGATAATTAAACGAATAATTGAAAAAAACTTTATTTACTAAAAAAAAACTTGTCAGGCATGAAAAAATCCCTATAATACGCGTCACTTGTTTTCGACGCCGACTTAGCTCAGTAGGTAGAGCAACTGACTTGTaa                                                                                                                                                                                  >  1:846021/1‑154 (MQ=255)                                                                                                                                                aaGATAAATTGTTTTTGTTTGATAATTAAACGAATAATTGAAAAAAACTTTATTTACTAAAAAAAAACTTGTCAGGCATGAAAAAATCCCTATAATACGCGTCACTTGTTTTCGACGCCGACTTAGCTCAGTAGGTAGAGCAACTGACTTGTAATCAGTAGGTCATCAGTTCGATTCCGATAGTCGGCACCATTTCTCGCAATACAAGTTCAAACTTAGAGTGGAGGGATTCCCGAGCGGCCAAAGGGAg                                                                                  >  1:758503/1‑250 (MQ=255)                                                                                                                                                            ttttgtttGATAATTAAACGAATAATTGAAAAAAACTTTATTTACTAAAAAAAAACTTGTCAGGCATGAAAAAATCCCTATAATACGCGTCACTTGTTTTCGACGCCGACTTAGCTCAGTAGGTAGAGCAACTGACTTGTAATCAGTAGGTCATCAGTTCGATTCCGATAGTCGGCACCATTTCTCGCAATACAAGTTCAAACTTAGAGTGGAGGGATTCCCGAGCGGCCAAAGGGAGCAGACTGTAAATc                                                                     <  1:252710/251‑1 (MQ=255)                                                                                                                                                                gtttgATAATTAAACGAATAATTGAAAAAAACTTTATTTACTAAAAAAAAACTTGTCAGGCATGAAAAAATCCCTATAATACGCGTCACTTGTTTTCGACGCCGACTTAGCTCAGTAGGTAGAGCAACTGACTTGTAATCAGTAGGTCATCAGTTCGATTCCGATAGTCGGCACCATTTCTCGCAATACAAGTTCAAACTTAGAGTGGAGGGATTCCCGAg                                                                                               <  1:339692/221‑1 (MQ=255)                                                                                                                                                                       aaTTAAACGAATAATTGAAAAAAACTTTATTTACTAAAAAAAAACTTGTCAGGCATGAAAAAATCCCTATAATACGCGTCACTTGTTTTCGACGCCGACTTAGCTCAGTAGGTAGAGCAACTGACTTGTAATCAGTAGGTCATCAGTTCGATTCCGATAGTCGGCACCATTTCTCGCAATACAAGTTCAAACTTAGAGTGGAGGGATTCCCGAGCGGCCAAAGGGAGCAGACTGTAAATCTGCCGGTTCAg                                                          >  1:685977/1‑251 (MQ=255)                                                                                                                                                                               gAATAATTGAAAAAAACTTTATTTACTAAAAAAAAACTTGTCAGGCATGAAAAAATCCCTATAATACGCGTCACTTGTTTTCGACGCCGACTTAGCTCAGTAGGTAGAGCAACTGACTTGTAATCAGTAGGTCATCAGTTCGATTCCGATAGTCGGCACCATTTCTCGCAATACAAGTTCAAACTTagag                                                                                                               <  1:867232/190‑1 (MQ=255)                                                                                                                                                                                       gAAAAAAACTTTATTTACTAAAAAAAAACTTGTCAGGGATGAAAAAATCCCTATAATACGCGTCACTTGTTTTCGACGCCGACTTAGCTCAGTAGGTAGAGCAACTGACTTGTAATCAGTAGGTCATCAGTTCGATTCCGATAGTCGGCACCATTTCTCGCAATACAAGTTCAAACTTAGAGTGGAGGGATTCCCGAGCGGCCAAAGGGAGCAGACTGTAAATCTGCCGGTTCAGCCTTCGAAGGTTCGaa                                          >  1:386546/1‑251 (MQ=255)                                                                                                                                                                                             aaCTTTATTTACTAAAAAAAAACTTGTCAGGCATGAAAAAATCCCTATAATACGCGTCACTTGTTTTCGACGCCGACTTAGCTCAGTAGGTAGAGCAACTGACTTGTAATCAGTAGGTCATCAGTTCGATTCCGATAGTCGGCACCATTTCTCGCAATACAAGTTCAAACTTagag                                                                                                               >  1:466266/1‑176 (MQ=255)                                                                                                                                                                                              aCTTTATTTACTAAAAAAAAACTTGTCAGGCATGAAAAAATCCCTATAATACGCGTCACTTGTTTTCGACGCCGACTTAGCTCAGTAGGTAGAGCAACTGACTTGTAATCAGTAGGTCATCAGTTCGATTCCGATAGTCGGCACCATTTCTCGCAATACAAGTTCAAACTTAGAGTGGAGGGATTCCCGAGCGGCCAAAGGGAGCAGACTGTAAATc                                                                     <  1:288055/217‑1 (MQ=255)                                                                                                                                                                                               cTTTATTTACTAAAAAAAAACTTGTCAGGCATGAAAAAATCCCTATAATACGCGTCACTTGTTTTCGACGCCGACTTAGCTCAGTAGGTAGAGCAACTGACTTGTAATCAGTAGGTCATCAGTTCGATTCCGATAGTCGGCACCATTTCTCGCAATACAAGTTCAAACTTAGAGTGGAGGGATTCCCGAGCGGCCAAAGGGAGCAGACTGTAAATCTGCCGGTTCAGCCTTCGAAGGTTCGAATCCTtctc                                  <  1:672806/251‑1 (MQ=255)                                                                                                                                                                                                                   cTTGTCAGGCATGAAAAAATCCCTATAATACGCGTCACTTGTTTTCGACGCCGACTTAGCTCAGTAGGTAGAGCAACTGACTTGTAATCAGTAGGTCATCAGTTCGATTCCGATAGTCGGCACCATTTCTCGCAATACAAGTTCAAACTTAGAGTGGAGGGATTCCCGAGCGGCCAAAGGGAGCAGACTGTAAATCTGCCGGTTCAGCCTTCGAAGGTTCGAATCCTTCTCCCTCCACCATTATCTAagag              <  1:662172/251‑1 (MQ=255)                                                                                                                                                                                                                   cTTGTCAGGCATGAAAAAATCCCTATAATACGCGTCACTTGTTTTCGACGCCGACTTAGCTCAGTAGGTAGAGCAACTGACTTGTAATCAGTAGGTCATCAGTTCGATTCCGATAGTCGGCACCATTTCTCGCAATACAAGTTCAAACTTAGAGTGGAGGGATTCCCGAGCGGCCAAAGGGAGCAGACTGTAAATCTGCCGGTTCAGCCTTCGAAGGTTCGAATCCTTCTCCCTCCACCATTATCTAagag              <  1:787097/251‑1 (MQ=255)                                                                                                                                                                                                                     tGTCAGGCATGAAAAAATCCCTATAATACGCGTCACTTGTTTTCGACGCCGACTTAGCTc                                                                                                                                                                                                           <  1:521223/60‑1 (MQ=255)                                                                                                                                                                                                                     tGTCAGGCATGAAAAAATCCCTATAATACGCGTCACTTGTTTTCGACGCCGACTTAGCTCAGTAGGTAGAGCAACTGACTTGTAATCAGTAGGTCATCAGTTCGATTCCGATAGTCGGCACCATTTCTCGCAATACAAGTTCAAACTTAGAGTGGAGGGATTCCCGAGCGGCCAAAGGGAGCAGACTGTAAATCTGCCGGTTCAGCCTTCGAAGGTTCGAATCCTTCTCCCTccacc                          <  1:882355/237‑1 (MQ=255)                                                                                                                                                                                                                          ggCATGAAAAAATCCCTATAATACGCGTCACTTGTTTTCGACGCCGACTTAGCTc                                                                                                                                                                                                           <  1:351909/55‑1 (MQ=255)                                                                                                                                                                                                                          ggCATGAAAAAATCCCTATAATACGCGTCACTTGTTTTCGACGCCGACTTAGCTCAGTAGGTAGAGCAACTGACTTGTAATCAGTAGGTCATCAGTTCGATTCCGATAGTCGGCACCATTTCTCGCAATACAAGTTCAAACTTAGAGTGGAGGGATTCCCGAGCGGCCACAGGGAGCAGACTGTAAATCTGCCGGTTCAGCCTTCGAAGGTTCGAATCCTTCTCCCTCCACCATTATCTAAGAGTAACCGa       >  1:434838/1‑251 (MQ=255)                                                                                                                                                                                                                           gCATGAAAAAATCCCTATAATACGCGTCACTTGTTTTCGACGCCGACTTAGCTCAGTAGGTAGAGCAACTGACTTGTAATCAGTAGGTCATCAGTTCGATTCCGATAGTCGGCACCATTTCTCGCAATACAAGTTCAAACTTAGAGTGGAGGGATTCCCGAGCGGCCAAAGGGAGCAGACTGTAAATCTGCCGGTTCAGCCTTCGAAGGTTCGAATCCTTCTCCCTCCACCATTATCTAAGAGTAACCGAt      >  1:213254/1‑251 (MQ=255)                                                                                                                                                                                                                                aaaaaaTCCCTATAATACGCGTCACTTGTTTTCGACGCCGACTTAGCTCAGTAGGTAGAGCAACTGACTTGTAATCAGTAGGTCATCAGTTCGATTCCGATAGTCGGCACCATTTCTCGCAATACAAGTTCAAACTTAGAGTGGAGGGATTCCCGAGCGGCCAAAGGGAGCAGACTGTAAATCTGCCGGTTCAGCCTTCGAAGGTTCGAATCCTTCTCCCTCCACCATTATCTAAGAGTAACCGATGGGa  >  1:568464/1‑250 (MQ=255)                                                                                                                                                                                                                                      |                                                                                                                                                                                                                                                    GCATAAGCTGCTCCGAAATAGAATGAGAAACCGGAAATTCCACGTGTTGAATTACTCTTTTTAGCCAAATATAGGTTGCAAAATATACCTTATAAATAAGAGATTGAGCAGTAAAAATTTAGCATAATCAGGTAAATTAGACGAAGATAAATTGTTTTTGTTTGATAATTAAACGAATAATTGAAAAAAACTTTATTTACTAAAAAAAAACTTGTCAGGCATGAAAAAAGCCCTATAATACGCGTCACTTGTTTTCGACGCCGACTTAGCTCAGTAGGTAGAGCAACTGACTTGTAATCAGTAGGTCATCAGTTCGATTCCGATAGTCGGCACCATTTCTCGCAATACAAGTTCAAACTTAGAGTGGAGGGATTCCCGAGCGGCCAAAGGGAGCAGACTGTAAATCTGCCGGCTCAGCCTTCGAAGGTTCGAATCCTTCTCCCTCCACCATTATCTAAGAGTAACCGATGGGA  >  1/1518718‑1519190

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