Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	NC_000913	3,725,176	T→G	E48A (GAG→GCG)	glyQ ←	glycine‑‑tRNA ligase subunit alpha

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	3,725,176	0	T	G	100.0%	209.8 / NA	65	E48A (GAG→GCG)	glyQ	glycine‑‑tRNA ligase subunit alpha
Reads supporting (aligned to +/- strand):  ref base T (0/0);  new base G (27/38);  total (27/38)

GTTTTCCCAGTTATCTTCCACGAAACGGATGTCGTGAATAGTCGGGTCCATGCCCAGCTCTTTCAGAGAACCGAGGTACAGCTCCTGAATATTGTCCGGCGATGGCTTAATGACCACCTGGAACTGATAGTAGTGCTGTAAACGGTTGGGGTTTTCGCCGTAGCGACCATCGGTCGGGCGACGAGAAGGCTGAACATAAGCAGCCGCCATCGGTTCTGGCCCCAGCTCGCGCAGACAGGTCATTGGGTGAGAGGTTCCCGCGCCGACTTCCATGTCCAATGGTTGAACAATGGTGCAGCCCTGGCGAGCCCAGTAATCCTGTAAGGTCAGGATCAAGCCCTGGAAGGTCCTGGTATCAAACTTTTGCATATTATTTCGTGCTGGATACGTGTGGATTTAAAGGAAGGGATCAGTATACCCGCTGGATGGAAGATATACAGTACGAAACGGGAAAAAGCAGGGCTTAACGCATG  >  NC_000913/3724950‑3725422                                                                                                                                                                                                                                   |                                                                                                                                                                                                                                                       gTTTTACCAGTTATCTTCCACGAAACGGATGTCGTGAAAAGACGGGTCCATGCCCAGCTCTTTCAGAGAACCGAGGTACAGATCCTGAATATTGTCCGGCGCTGGCTTAATGACAACCTGGAACTGATAGTAGTGCTGTAAACGGTTGGGGTTTTCGCCGTAGCGACCATCGGTCGGTCGACGAGAAGGCTGAACATAACCAGCAGCCATCGGTTCTGGCCCCAGCGCGCGCAGCCAGGtacttggggg                                                                                                                                                                                                                                  >  2:496816/1‑240 (MQ=255)       ccAGTTATCTTCCACGAAACGGATGTCGTGAATAGTCGGGTCCATGCCCAGCTCTTTCAGAGAACCGAGGTACAGCTCCTGAATATTGTCCGGCGATGGCTTAATGACCACCTGGAACTGATAGTAGTGCTGTAAACGGTTGGGGTTTTCGCCGTAGCGACCATCGGTCGGGCGACGAGAAGGCTGAACATAAGCAGCCGCCATCGGTTCTGGCCCCAGCGCGCGCAGACAGGTCATTGGGTGAGAGGTTc                                                                                                                                                                                                                          <  1:121631/251‑1 (MQ=255)       ccAGTTATCTTCCACGAAACGGATGTCGAGAATAATCGGGTCCATGCCCAGCTCTTTCAGAGAACCGAGGTACAGCTCCTGAATATTGTCCGGCGATGGCTTAATGACCACCTGGAACTGATAGTAGTGCTGTAAACGGTTGGGGTTTTCGCCGTAGCGACCTTCGGTCAGGCGACGAGAAGGCTGAACATAAGCAGCTGCCATCGGTTCTGGCCCCAGCGGGCGCAGACAGGTACTTGGGGGAGCGGTTc                                                                                                                                                                                                                          >  1:134172/1‑251 (MQ=255)           ttATCTTACACGAAACGGCTGGCGTGAATAGTCGGGTCCATGCCCAGCTCTTTCAGAGAACCGAGGTACAGCTCCTGACTATTGTCCGGCGCTGGCTTAATGACCACCTGGAACTGATAGTAGTGCTGTAAACGGTTGGGGTTTTCGCCGTAGCGACCATCGGTCGGGCGACGAGCAGGCTGACCATAAGCAGCCGCCATCGGTTCTGGCCCCAGCGCTGGCAGACAGGTCATTGGGTGCGAGGGTCCcgc                                                                                                                                                                                                                      >  2:75824/1‑251 (MQ=255)             aTCTTCCACGAAACGGATGTCGTGAATAGTCGGGTCCATGCCCAGCTCTTTCAGAGAACCGAGGTACAGCTCCTGAATATTGTCCGGCGATGGCTTAATGACCACCTGGAACTGATAGTAGTGCTGTAAACGGTTGGGGTTTTCGCCGTAGCGACCATCGGTCGGGCGACGAGAAGGCTGAACATAAGCAGCCGCCATCGGTTCTGGCCCCAGCGCGCGCAGACAGGTCATTGGGGGAGAGGTTCCCgcgc                                                                                                                                                                                                                    >  1:334592/1‑251 (MQ=255)               catCCACAAAACGGATGTCGTGAATTGTCGGGTCCATGCACAGCTCTTTCAGAGAACCGAGGGACAGCGCCTGAATATTGTCCGGCGAGGGCTTAATTACCACCTGGATCTGATAGTAGTGCTGTAAACGGTTGGGGTTTTCGCCGACGCGACCATCGGTCGGGCGACGAGAAGGCTGAACATAAGCAGCCGCCATCGGTTCTGGCCCCAGCGCGCGCAGACTTGTCATTGGGTGAGAGGTTCCCGCGCCg                                                                                                                                                                                                                  <  2:314440/249‑1 (MQ=255)                 tCCACGAAACGGATGTCGTGAATAGTCGGGTCCATGCCCAGCTCTTTCAGAGAACCGAGGTACAGCTCCTGAATATTGTCCGGCGATGGCTTAATGACCACCTGGAACTGATAGTAGTGCTGTAAACGGTTGGGGTTTTCGCCGTAGCGACCATCGGTCGGGCGACGAGAAGGCTGAACATAAGCAGCCGCCATCGGTTCTGGCCCCAGCGCGCGCAGACAGGTCATTGGGTGAGAGGTTCCCGCGcc                                                                                                                                                                                                                   <  1:118764/248‑1 (MQ=255)                      gAAACGGATGTCGTGAATAGTCGGGTCCATGCCCAGCTCTTTCAGAGAACCGAGGTACAGCTCCTGAATATTGTCCGGCGATGGCTTAATGACCACCTGGAACTGATAGAAGTGCTGTAAACGGTTGGGGGTTTCGCCGTAGCGACCATCGGTCGGGCGACGAGAAGGCTGAACATAAGCAGCCGCCAACGGGTCTGGCCCCAGCGCGCGCAGACAGGTCATTGGGAGAGAGGGTCCCGCGCCGACTTccc                                                                                                                                                                                                           >  2:579046/1‑250 (MQ=255)                       aatcGGATGTCGTGAATAGGCGGGTCCATTCCCAGCGCTTGCATGGAACCGAGGTCCAGCTCCTGATGAGTGGTCGGGGATGGTTTAATGGCCACCTGGATCTGATGGTAGTGCTGGAAACGGTTGGGGCTTTCGCCGTGGCGACCATCGGTCGGGCGACGAGAAGGCTGAACATAAGCGGCCGCCATCGGTTCTTGCCCCAGCGCGCGCAGACAGGTCATTGAGTGAGAGGTTCTCGGGCCGACTTCCAt                                                                                                                                                                                                          <  2:33862/248‑1 (MQ=255)                       aaaCGGATGTCGTGAATAGTCGGGTCCATGCCCAGCTCTTTCAGAGAACCGAGGTACAGCTCCTGAATATTGTCCGGCGATGGCTTAATGACCACCTGGAACTGATAGTAGTGCTGTAAACGGTTGGGGTTTTCGCCGTAGCGACCATCGGTCGGGCGACGAGAAGGCTGAACATAAGCAGCCGCCATCGGTTCTGGCCCCAGCGCGCGCAGACAGGTCATTGGGTGAGAGGTTCCCGCGCCGACTTCCAt                                                                                                                                                                                                          >  1:176781/1‑251 (MQ=255)                           ggATGTCGTGAATAGTCGGGTCCATGCCCAGCTCTTTCAGAGAACCGAGGTACAGCTCCTGAATATTGTCCGGCGATGGCTTAATGACCACCTGGAACTGATAGTAGTGCTGTAAACGGTTGGGGTTTTCGCCGTAGCGACCATCGGTCGGGCGACGAGAAGGCTGAACATAAGCAGCCGCCATCGGTTCTGGCCCCAGCGCGCGCAGACAGGTCATTGGGTGAGAGGTTCCCGCGCCGACTTCCATGTcc                                                                                                                                                                                                      >  1:374021/1‑251 (MQ=255)                               gTCGTGAAGAGTTGTGTCCATGCCCAGCTCTTGCAGAGAACCTAGGGACAGCCCCGGAATATTGTCCGGCGATGGCTTAATTGCCTCCTGCAACTGATAGTAGTGCTGTAAACGGTTGGGGTTTTCGCCGTAGCGACCATCTGTCGGGCTACTAGAAGGCTGACCATCAGCAGCCGCCATCGGTTCTGGCCCCAACGCGCGCAGACAGGTCATTGGGTTAGATGTTCCCGCGCCTACTTCCATGTCCAAt                                                                                                                                                                                                   <  2:520365/250‑1 (MQ=255)                                 cGTGAATAGCCGGGTCAATGCCCAGCTCTTTCAGAGAACCGAGGTACAGCTCCTTAATATTGTCCGGCGATGGCTTAATGACCACCGGGAACTGATAGTAGTGCTGGAATCGGTTGGGGTTTTCGCCGTAGCGACCATCGGTCGGGCGACGAGAAGGCTGAACATAAGCAGCCGCCATCGGTTCTGGCCCCAGCGCGCGCAGACAGGTCATTGGGTGAGAGGTTCCCGCGCCGACTTCCATGTCCAATgg                                                                                                                                                                                                 <  2:7347/250‑1 (MQ=255)                                   tGAATAGTCGGGTCCATGCCCAGCTCTTTCAGAGAACCGAGGTACAGCTCCTGAATATTGTCCGGCGATGGCTTAATGACCACCTGGAACTGATAGTAGTGCTGTAAACGGTTGGGGTTTTCGCCGTAGCGACCATCGGTCGGGCGACGAGAAGGCTGAACATAAGCAGCAGCCATCGGTTCTGGCCCCAGCGCGCGCAGACAGGTCATTGGGGGAGAGGTTCCCGCGCCGACTTCCATGTCCACTGGTTg                                                                                                                                                                                              >  1:221886/1‑251 (MQ=255)                                        aGTCGGGTCCATGCCCAGCGCTTTCAGAGAACCGAGGTACAGCTCCTGAATATTGTCCGGCGATGGCTTAATGACCACCTGGAACTGATAGTAGTGCTGTAAACGGTTGGGGTTTTCGCCGTAGCGACCATCGGTCGGGCGACGAGAAGGCTGAACATAAGCAGCCGCCATCGGTTCTGGCCCCAGCGCGCGCAGACAGGTCATTGGGTGAGAGGTTCCCGCGCCGACTTCCATGTCCAATGGTTGAACaa                                                                                                                                                                                         <  1:28398/251‑1 (MQ=255)                                                    gCCCAGATCTTTCAGAGAACCGAGGAACAGCTCCTTAATAGTGGCCGGCGATGGCGTAATGACCACCTGGAACTGATAGGAGTGCTGTAAACGGTTGGGGTTTTCGCCGTAGCGACCAGCGGTCGGGCGACGAGAAGGCTGAACATAAGCAGCCGCCATCGGTTCTGGCCCCAGCGCGCGCAGACAGGTCATTGGGTGAGAGGTTCCCGCGCCGACTTCCATGTCCAATGGTTGAACAATGGTGCAGCCCt                                                                                                                                                                             <  2:147575/251‑1 (MQ=255)                                                        aGCTCTTTCAGAGAACCGAGGTACAGCTCCTGAATATTGTCCGGCGATGGCTTAATGACCACCTGGAACTGATAGTAGTGCTGTAAACGGTTGGGGTTTTCGCCGTAGCGACCATCGGTCGGGCGACGAGAAGGCTGAACATAAGCAGCCGCCATCGGTTCTGGCCCCAGCGCGCGCAGACAGGTCATTGGGTGAGAGGTTCCCGCGCCGACTTCCATGTCCAATGGTTGAACAATGGTGCAGCCCTGGCg                                                                                                                                                                         >  2:490103/1‑251 (MQ=255)                                                        aGCTCTTTCAGAGAACCGAGGTACAGCTCCGGAATATTGTCCGGCGATGGCTTAATGACCACCTGGAACTGATAGTAGTGCTGTAAACGGTTGGGGTTTTCGCCGTAGCGACCATCGGTCGGGCGACGAGAAGGCTGAACATAAGCAGCCGCCATCGGTTCTGGCCCCAGCGCGCGCAGACAGGTCATTGGGTGAGAGGTTCCCGCGCCGACTTCCATGTCCAATGGTTGAACAATGGTGCAGCCCTGGCg                                                                                                                                                                         <  2:169182/251‑1 (MQ=255)                                                           tctTTCAGAGAACCGAGGTACAGCTCCTGAATATTGTCAGGCGATGGCTTAATGACCACCTGGAACTGATAGTAGTGCTGTAAACGGTTGGGGTTTTCGCCGTAGCGACCATCGGTCGGGCGACGAGAAGGCTGAACATAAGCAGCCGCCATCGGTTCTGGCCCCCGCGCGCGCAGACAGGTCATTGGGTGAGAGGTTCCCGCGCCGACTTCCATGTCCAATGGTTGAACAATGGGGCAGCCCTGGCGAGc                                                                                                                                                                      >  2:203764/1‑251 (MQ=255)                                                            ctTTCAGAGAACCGAGGTACAGCTCCTGAATATTGTCCGGCGATGGCTTAATGACCACCTGGAACTGATAGTAGTGCTGTAAACGGTTGGGGTTTTCGCCGTAGCGACCATCGGTCGGGCGACGAGAAGGCTGAACATAAGCAGCCGCCATCGGTTCTGGCCCCAGCGCGCGCAGACAGGTCATTGGGTGAGAGGTTCCCGCGCCGACTTCCATGTCCAATGGTTGAACAATGGTGCAGCCCTGGCGAg                                                                                                                                                                       <  1:9439/249‑1 (MQ=255)                                                               tCAGAGAACCGAGGTACAGCTCCTGAATATTGTCCGGCGATGGCTTAATGACCACCTGGAACTGATAGTAGTGCTGTAAACGGTTGGGGTTTTCGCCGTAGCGACCATCGGTCGGGCGACGAGAAGGCTGAACATAAGCAGCCGCCATCGGTTCTGGCCCCAGCGCGCGCAGACAGGTCATTGGGTGAGAGGTTCCCGCGCCGACTTCCATGTCCAATGGTTGAACAATGGTGCAGCCCTGGCGAGCCCa                                                                                                                                                                   <  1:167048/250‑1 (MQ=255)                                                                 agagAACCGAGGTACAGCTCCTGAATATTGTCCGGCGATGGCTTAATGACCACCTGGAACTGATAGTAGTGCTGTAAACGGTTGGGGTTTTCGCCGTAGCGACCATCGGTCGGGCGACGAGAAGGCTGAACATAAGCAGCCGCCATCGGTTCTGGCCCCAGCGCGCGCAGACAGGTCATTGGGTGAGAGGTTCCCGCGCCGACTTCCATGTCCAATGGTTGAACAATGGTGCAGCCCTGGCGAGCCCAGt                                                                                                                                                                 >  1:170207/1‑250 (MQ=255)                                                                     acccGAGGTACAGCTCCTGAATATTGCCCGGCGATGGCGTAATTACCACCTGGAACTGATAGTAGTGCTGTAAACGGTTGGGGTTTTCTCCGTAGCGACCATCGGGCGGGCGACGAGAAGGCTGAACATAAGCAGCCGCCATCGGTTCTGGCCCCAGCGCGCGCAGTCAGGTCATTTGGTGAGAGGTTCCCGCGCCGACTTCCATGTCCAATGGTTGAACAATGGTGCAGCCCTGGCGAGCCCAGTAATcc                                                                                                                                                            <  2:24034/249‑1 (MQ=255)                                                                              aCAGCTCCTGAATATTTTCCGGCGATGACTTAATGACCCCCTTGAACGGATAGTAGTGCTGTAAACGGTTGGGGTTTTCGCCGTAGCGACCATCGGTCGGGCGACGAGAAGGCTGAACATAAGCAGCCGCCATCGGTTCTGGCCCCAGCGCGCGCAGACAGGTCATTGGGTGAGAGGTTCCCGCGCCGACTTCCATGTCCAATGGTTGAACAATGGTGCAGCCCTGGCGAGCCCAGTAATCCTGTAAGGTc                                                                                                                                                   <  2:395102/251‑1 (MQ=255)                                                                                 gCTCCTGAATATTGTCCGGCGATGGCTTAATGACCACCTGGAACTGATAGTAGTGCGGTAAACGGTTGGGGTTTTCGCCGTAGCGACCATCGGTCGGGCGACGAGAAGGCTGAACATAAGCAGCCGCCATCGGTTCTGGCCCCAGCGCGCGCAGACAGGTCATTGGGTGAGAGGTTCCCGCGCCGACTTCCATGTCCAATGGTTGAACAATGGTGCAGCCCTGGCGAGCCCAGTAATCCTGTAAGGTCAgg                                                                                                                                                <  1:518876/251‑1 (MQ=255)                                                                                    ccTGAATATTGTCCGGGGATGGCTTAATGACCACCTGGAACTGATAGTAGTGATGTAAACGGTTGGGGTTTTCGCCGTAGCGACCATCGGTCGGGCGACGAGAAGGCTGAACATAAGCAGCCGCCATCGGTTCTGGCCCCAGCGCGCGCAGACAGGTCATTGGGTGAGAGGTTCCCGCGCCGACTTCCATGTCCAATGGTTGAACAATGGTGCAGCCCTGGCGAGCCCAGTAATCCTGTAAGTTCAGGATc                                                                                                                                             <  2:290626/251‑1 (MQ=255)                                                                                        aaTATTGTCCGGCGATGGCTTAATGCCCACCTGGAACTGATAGTAGTGCTGTAAACGGTTGGGGTTTTCGCCGTAGCGACCATCGGTCGGGCGACGAGAAGGCTGCACATAAGCAGCCGCCATCGGTTCTGGCCCCAGCGCGCGCAGACAGGTAATTGGGTGAGAGGTTCCCGCGCCGACTTCCATGTCCAATGGTTGACCAATGGTGCAGCCCTGGCGAGCCAAGTAATCCTGTAAGGTCAGGATCAAGc                                                                                                                                         >  2:363797/1‑251 (MQ=255)                                                                                         atatTGTCCGGCGATGGCTTAATGACCACCTGGAACTGCTAGTAGTGCTGTAAACGGTTGGGGTTTTCGCCGTAGCGACCATCGGTCGGGCGACGAGAAGGCTGAACATAAGCAGCCGCCATCGGTTCTGGCCCCAGCGCGCGCAGACAGGTCCTTGGGTGAGAGGTTCCCGCGCCGACTTCCATGTCCAATGGTTGAACAATGGTGCAGCCCTGGCGAGCCCAGTCATCCTGTAAGGTCAGGATCATGcc                                                                                                                                        >  2:188020/1‑251 (MQ=255)                                                                                           atTGTCCGGCGATGGCTTAATGACCACCTGGAACTGATAGTAGTGCTGTAAACGGTTGGGGTTTTCGCCGTAGCGACCATCGGTCGGGCGACGAGAAGGCTGAACATAAGCAGCCGCCATCGGTTCTGGCCCCAGCGCGCGCAGACAGGTCATTGGGTGAGAGGTTCCCGCGCCGACTTCCATGTCCAATGGTTGAACAATGGTGCAGCCCTGGCGAGCCCAGTAATCCTGTAAGGTCAGGATCAAGCCCt                                                                                                                                      >  1:404339/1‑251 (MQ=255)                                                                                             tGTCCGGCGATGGCTTAATGACCACCTGGAACTGATAGTAGTGCTGTAAACGGTTGGGGTTTTCGCCGTAGCGACCATCGGTCGGGCGACGAGAAGGCTGAACATAAGCAGCCGCCATCGGTTCTGGCCCCAGCGCGCGCAGACAGGTCATTGGGTGAGAGGTTCCCGCGCCGACTTCCATGTCCAATGGTTGAACAATGGTGCAGCCCTGGCGAGCCCAGTAATCCTGTAAGGTCAGGATCAAGCCCTgg                                                                                                                                    >  1:553911/1‑251 (MQ=255)                                                                                              gTCCGGCGATGGCTTAATGACCACCTGGAACTGATAGTAGTGCTGTAAACGGTTGGGGTTTTCGCCGTAGCGACCATCGGTCGGGCGACGAGAAGGCTGAACATAAGCAGCCGCCATCGGTTCTGGCCCCAGCGCGCGCAGACAGGTCATTGGGTGAGAGGTTCCCGCGCCGACTTCCATGTCCAATGGTTGAACAATGGTGCAGCCCTGGCGAGCCCAGTAATCCTGTAAGGTCAGGATCAAGCCCTGGa                                                                                                                                   <  1:252889/251‑1 (MQ=255)                                                                                                ccGGCGATGGCTTAATGACCACCTGGAACTGATAGTAGTGCTGTAAACGGTTGGGGTTTTCGCCGTAGCGACCATCGGTCGGGCGACGAGAAGGCTGAACATAAGCAGCCGCCATCGGTTCTGGCCCCAGCGCGCGCAGACAGGTCATTGGGTGAGAGGTTCCCGCGCCGACTTCCATGTCCAATGGTTGAACAATGGTGCAGCCCTGGCGAGCCCAGTAATCCTGTAAGGTCAGGATCAAGCCCTGGAAg                                                                                                                                 <  1:490103/251‑1 (MQ=255)                                                                                                    cGATGGCTTAATGACCCCCTGGAACTGATAGTAGTGCTGTAAACGGTTGGGGTTTTCGCCGTAGCGACCATCGGTCGGGCGACGAGAAGGCTGAACATAAGCAGCCGCCATCGGTTCTGGCCCCAGCGCGCGCAGACAGGTCATTGGGTGAGAGGTTCCCGCGCCGACTTCCATGTCCAATGGTTGAACAATGGTGCAGCCCTGGCGAGCCCAGTAATCCTGTAAGGTCAGGATCAAGCCCTGGAAGGTcc                                                                                                                             <  2:352797/251‑1 (MQ=255)                                                                                                            taatggccacccgacacggATAGTGGTGCTGTAAACGGTTGGGGTTTTCGCCGTAGCGACCAGCGGGCGGGCGACGAGAAGGCTGTACAGAAGCAGCCGCCATCGGTTCTGGCCCCAGCGCTCGCAGACAGGTCATTGGGTGAGAGGGTCCCGCGCCGACTTCCATGTCCAACGGTTTAACAATGGTGCAGCCCTGGCGAGCCCAGTAATCCTATAATGTCAGGATCAAGTCCGGTAAGGTCTTGGTAt                                                                                                                       <  2:244139/231‑1 (MQ=255)                                                                                                            tAATGACCACCTGGAACTGATAGTAGTGCTGTAAACGGTTGGGGTTTTCGCCGTAGCGACCATCGGTCGGGCGACGAGAAGGCTGAACATAAGCAGCCGCCATCGGTTCTGGCCCCAGCGCGCGCAGACAGGTCATTGGGTGAGAGGTTCCCGCGCCGACTTCCATGTCCAATGGTTGAACAATGGTGCAGCCCTGGCGAGCCCAGTAATCCTGTAAGGTCAGGATCAAGCCCTGGAAGGTCCTGGTATc                                                                                                                      >  1:275325/1‑250 (MQ=255)                                                                                                             aaTGACCACCTGGAACTGATAGTAGTGCTGTAAACGGTTGGGGTTTTCGCCGTAGCGACCATCGGTCGGGCGACGAGAAGGCTGAACATAAGCAGCCGCCATCGGTTCTGGCCCCAGCGCGCGCAGACAGGTCATTGGGTGAGAGGTTCCCGCGCCGACTTCCATGTCCAATGGTTGAACAATGGTGCAGCCCTGGCGAGCCCAGTAATCCTGTAAGGTCAGGATCAAGCCCTGGAAGGTCCTGGTATCa                                                                                                                     >  1:296478/1‑250 (MQ=255)                                                                                                                gACCACCTGGAACTGATAGTAGTGCTGTAAACGGTTGGGGTTTTCGCCGTAGCGACCATCGGTCGGGCGACGAGAAGGCTGAACATAAGCAGCCGCCATCGGTTCTGGCCCCAGCGCGCGCAGACAGGTCATTGGGTGAGAGGTTCCCGCGCCGACTTCCATGTCCAATGGTTGAACAATGGTGCAGCCCTGGCGAGCCCAGTAATCCTGTAAGGTCAGGATCAAGCCCTGGAAGGTCCTGGTATCAAACt                                                                                                                 <  1:276320/251‑1 (MQ=255)                                                                                                                 accaccTGGAACTGATAGTAGTGCTGTAAACGGGTGGGGTTTTCGCCGGAGCGACCAGCGGTCGGGCGACGAGAAGGCTGAACATAAGCAGCCGCCATCGGTTCTGGCCCCAGCGCGCGCAGACAGGTCATTGGGTGAGAGGTGCCCGCGCCGACTTCCATGTCCAATGGTTGAACAATGGTGCAGCCCTGGCGAGCCCAGTAATCCTGTAAGGTCAGGATCAAGCCCTGGAAGGTCCTGGTATCAAACtt                                                                                                                <  2:260429/251‑1 (MQ=255)                                                                                                                        ggAACTGATAGTAGTGCTGTAAACGGTTGGGGTTTTCGCCGTAGCGACCATCGGTCGGGCGACGAGAAGGCTGAACATAAGCAGCCGCCATCGGTTCTGGCCCCAGCGCGCGCAGACAGGTCATTGGGTGAGAGGTTCCCGCGCCGACTTCCATGTCCAATGGTTGAACAATGGTGCAGCCCTGGCGAGCCCAGTAATCCTGTAAGGTCAGGATCAAGCCCTGGAAGGTCCTGGTATCAAACTTTTGCata                                                                                                         <  1:579046/251‑1 (MQ=255)                                                                                                                         gAACTGATAGTAGTGCTGTAAACGGTTGGGGTTTTCGCCGTAGCGACCATCGGTCGGGCGACGAGAAGGCTGAACATAAGCAGCCGCCATCGGTTCTGGCCCCAGCGCGCGCAGACAGGTCATTGGGTGAGAGGTTCCCGCGCCGACTTCCATGTCCAATGGTTGAACAATGGTGCAGCCCTGGCGAGCCCAGTAATCCTGTAAGGTCAGGATCAAGCCCTGGAAGGTCCTGGCATCAACCTTTTGCatat                                                                                                        >  2:246761/1‑251 (MQ=255)                                                                                                                              gATAGTAGTGCTGTAAACGGTTGGGGTTTTCGCCGTAGCGACCATCGGTCGGGCGACGAGAAGGCTGAACATAAGCAGCCGCCATCGGTTCTGGCCCCAGCGCGCGCAGACAGGTCATTGGGTGAGAGGTTCCCGCGCCGACTTCCATGTCCAATGGTTGAACAATGGTGCAGCCCTGGCGAGCCCAGTAATCCTGTAAGGTCAGGATCAAGCCCTGGAAGGTCCTGGTATCAAACTTTTGCAtattat                                                                                                     <  1:554124/249‑1 (MQ=255)                                                                                                                                aagtagtGCTGTAAACGGTTGGGGTTTTCGCCGTAGCGACCATCGGTCGGGCGACGAGAAGGCTGAACATAAGCAGCCGCCATCGGTTCTGGCCCCAGCGCGCGCAGACAGGTCATTGGGTGAGAGGTTCCCGCGCCGACTTCCATGTCCAATGGTTGAACAATGGTGCAGCCCTGGCGAGCCCAGTAATCCTGTAAGGTCAGGATCAAGCCCTGGAAGGTCCTGGTATCAAACTTTTGCATATTATTTCg                                                                                                 <  1:188020/250‑1 (MQ=255)                                                                                                                                 agtatggCTGTAAACGGTTGGGGTTTGCGCCGTAGCGACCATCGGTCGGGCGCCGAGAAGGCTGAACATAAGCATCCGCCATCGGTTCTGGCCCCAGCGCGCGCAGACAGGTCATTGGGTGAGAGGTTCCCGCGCCGGCTTCCATGTCCAATGGTTGAACAATGGTGCAGCCCTGGCGAGCCCAGTAATCCTTTAAGGTCAGGATCAAGCCCTGGAAGGTTCTGGTATCAAACTTTTGCATATTATTTCg                                                                                                 <  1:327663/244‑1 (MQ=255)                                                                                                                                               cgtttGGGGTTTTCGCCGTAGCGACCAGCGGTCGGGCGACGAGAAGGCTTAACATAAGCAGCCGTCATCGGTTCTGGCCCCAGCGCGCGCAGACAGGTCATTGGGTGAGAGGTTCCCGCGCCGACTTCCATGTCCAATGGTTGAACAATGGTGCAGCCCTGGCGAGCCCAGTAATCCTGTAAGGTCAGGATCAAGCCCTGGAAGGTCCTGGTATCAAACTTTTGCATATTATTTCGTGCTGGATACGtgtg                                                                                  <  2:172456/248‑1 (MQ=255)                                                                                                                                                   tGGGGTTTTCGCCGTAGCGACCATCGGTCGGGCGACGAGAAGGCTGAACATAAGCAGCCGCCATCGGTTCTGGCCCCAGCGCGCGCAGACAGGTCATTGGGTGAGAGGTTCCCGCGCCGACTTCCATGTCCAATGGTTGAACAATGGTGCAGCCCTGGCGAGCCCAGTAATCCTGTAAGGTCAGGATCAAGCCCTGGAAGGTCCTGGTATCAAACTTTTGCATATTATTTCGTGCTGGATACGTGTGGAt                                                                               <  1:93257/250‑1 (MQ=255)                                                                                                                                                    ggggTTTTCGCCGTAGCGACCATCGGTCGGGCGACGAGAAGGCTGAACATAAGCAGCCGCCATCGGTTCTGGCCCCAGCGCGCGCAGACAGGTCATTGGGTGAGAGGTTCCCGCGCCGACTTCCATGTCCAATGGTTGAACAATGGTGCAGCCCTGGCGAGCCCAGTAATCCTGTAAGGTCAGGATCAAGCCCTGGAAGGTCCTGGTATCAAACTTTTGCATATTATTTCGTGCTGGATACGTGTGGAttt                                                                             >  1:80883/1‑251 (MQ=255)                                                                                                                                                          ttCGCCGTAGCGACCATCGGTCGGGCGACGAGAAGGCTGAACATAAGCAGCCGCCATCGGTTCTGGCCCCAGCGCGCGCAGACAGGTCATTGGGTGAGAGGTTCCCGCGCCGACTTCCATGTCCAATGGTTGAACAATGGTGCAGCCCTGGCGAGCCCAGTAATCCTGTAAGGTCAGGATCAAGCCCTGGAAGGTCCTGGTATCAAACTTTTGCATATTATTTCGTGCTGGATACGTGTGGATTTAaagg                                                                        >  2:31429/1‑250 (MQ=255)                                                                                                                                                                        cagcggtcggGGGACGAGAAGGCTGAACATAAGCAGCCGCCATTGGTGCTGGCCCCAGCGCGCGCAGACAGGTCATTGTGTGAGAGGTGCCCGCGCCGACTTCCATGTCCAATGGTTGAACAATGGTGCAGCCCTGGCGAGCCCAGTAATCCTGTAAGGTCAGGATCAAGCCCTGGAAGGTCCTGGTATCAAACTTTTGCATATTATTTCGTGCTGGATACGTGTGGATTTAAAGGAAGGGATCAGTATAc                                                         <  2:462481/248‑1 (MQ=255)                                                                                                                                                                          tcggtcggGCGACGAGAAGGCTGAACATAAGCAGCCGCCATCGGTTCTGGCCCCAGCGCGCGCAGACAGGTCATTGGGTGAGAGGTTCCCGCGCCGACTTCCATGTCCAATGGTTGAACAATGGTGCAGCCCTGGCGAGCCCAGTAATCCTGTAAGGTCAGGATCAAGCCCTGGAAGGTCCTGGTATCAAACTTTTGCATATTATTTCGTGCTGGATACGTGTGGATTTAAAGGAAGGGATCAGTATAccc                                                       >  2:467687/1‑251 (MQ=255)                                                                                                                                                                             gtcggGCGACGAGTAGGCTGAACACAAGCAGCCGCCAGCGGTTCTGGCCCCAGCGCGCTCAGACAGGTCATTGGGTGAGAGGTTCCCGCGCCGACTGCCATGTCCAATGGGTGAACAATGGTGCAGCCCTGGCGAGCCCAGTAATCCTGTAAGGTCAGGATCAAGCCCTGGAAGGTCCTTGTATCAAACTTTTGCATATTATTTCGTGCTGGATACGTGTGGATTTAAAGGAAGGGATCAGTATTccc                                                       <  2:393480/248‑1 (MQ=255)                                                                                                                                                                                gggCGACGAGAAGGCTGAACATAAGCAGCCGCCATCGGTTCTGGCCCCAGCGCGCGCAGACAGGTCATTGGGTGAGAGGTTCCCGCGCCGACTTCCATGTCCAATGGTTGAACAATGGTGCAGCCCTGGCGAGCCCAGTAATCCTGTAAGGTCAGGATCAAGCCCTGGAAGGTCCTGGTATCAAACTTTTGCATATTATTTCGTGCTGGATACGTGTGGATTTAAAGGAAGGGATCAGTATACCCGCtgg                                                  >  1:193805/1‑250 (MQ=255)                                                                                                                                                                                 ggCGACGAGAAGGCTGAACATAAGCAGCCGCCATCGGTTCTGGCCCCAGCGCGCGCAGACAGGTCATTGGGTGAGAGGTTCCCGCGCCGACTTCCATGTCCACTGGTTGAACAATGGTGCAGCCCTGGCGAGCCCAGTAATCCTGTAAGGTCAGGATCAAGCCCTGGAAGGTCCTGGTATCAAACTTTTGCATATTATTTCGTGCTGGATACGTGTGGATTTAAAGGAAGGGATCAGTATACCCGCtggat                                                >  1:559202/1‑251 (MQ=255)                                                                                                                                                                                    gacgaGAAGGCTGAACATAAGCAGCCGCCATCGGTTCTGGCCCCAGCGCGCGCAGACAGGTCATTGGGTGAGAGGTTCCCGCGCCGACTTCCATGTCCAATGGTTGAACAATGGTGCAGCCCTGGCGAGCCCAGTAATCCTGTAAGGTCAGGATCAAGCCCTGGAAGGTCCTGGTATCAAACTTTTGCATattatt                                                                                                    >  1:231590/1‑196 (MQ=255)                                                                                                                                                                                    gacgaGAAGGCTGAACATAAGCAGCCGCCATCGGTTCTGGCCCCAGCGCGCGCAGACAGGTCATTGGGTGAGAGGTTCCCGCGCCGACTTCCATGTCCAATGGTTGAACAATGGTGCAGCCCTGGCGAGCCCAGTAATCCTGTAAGGTCAGGATCAAGCCCTGGAAGGTCCTGGTATCAAACTTTTGCATATTATTTCGTGCTGGATACGTGTGGATTTAAAGGAAGGGATCAGTATACCCGCtggatgga                                             <  1:84142/251‑1 (MQ=255)                                                                                                                                                                                    gacgaGAAGGCTGAACATAAGCAGCCGCCATCGGTTCTGGCCCCAGCGCGCGAAGACAGGACATTGGGTGAGAGGTTCCCGCGCCGACTTCCATGTCCAATGGTTGAACAATGGTGCAGCCCTGGCGAGCCCAGTAATCCTGTAAGGTCAGGATCAAGCCCTGGAAGGTCCTGGTATCAAACTTTTGCATATTATTTCGTGCTGGATACGTGTGGATTTAAAGGAAGGGATCAGTATACCCGCtggatgg                                              >  1:551494/1‑250 (MQ=255)                                                                                                                                                                                    gacgaGAAGGCGGAACATAAGCAGCCGCCATCGGTTCTGGCCCCAGCGCGCGCAGACAGGTCATTGGGTGAGAGGTTCCCGCGCCGACTTCCATGTCCAATGGTTGAACAATGGTGCAGCCCTGGCGAGCCCAGTAATCCTGTAATGTCAGGATCAAGCCCTGGAAGGTCCTGGTATCAAACTTTTGCATattatt                                                                                                    <  2:231590/196‑1 (MQ=255)                                                                                                                                                                                         gaAGGCTGAACATAAGCAGCCGCCATCGGTTCTGGCCCCAGCGCGCGCAGACAGGTCATTGGGTGAGAGGTTCCCGCGCCGACTTCCATGTCCAATGGTTGAACAATGGTGCAGCCCTGGCGAGCCCAGTAATCCTGTAAGGTCAGGATCAAGCCCTGGAAGGTCCTGGTATCAAACTTTTGCATATTATTTCGTGCTGGATACGTGTGGATTTAAAGGAAGGGATCAGTATACCCTCTGGATGGAAGata                                        <  1:467687/251‑1 (MQ=255)                                                                                                                                                                                                gAACATAAGCAGCCGCCATCGGTTCTGGCCCCAGCGCGCGCAGACAGGTCATTGGGTGAGAGGTTCCCGCGCCGACTTCCATGTCCAATGGTTGAACAATGGTGCAGCCCTGGCGAGCCCAGTAATCCTGTAAGGTCAGGATCAAGCCCTGGAAGGTCCTGGTATCAAACTTTTGCATATTATTTCGTGCTGGATACGTGTGGATTTAAAGGAAGGGATCAGTATACCCGCTGGATGGAAGATATACAGTa                                 <  2:374021/251‑1 (MQ=255)                                                                                                                                                                                                 aaCATAAGCAGCCGCCATCGGTTCTGGCCCCAGCGCGCGCAGACAGGTCATTGGGTGAGAGGTTCCCGCGCCGACTTCCATGTCCAATGGTTGAACAATGGTGCAGCCCTGGCGAGCCCAGTAATCCTGTAAGGTCAGGATCAAGCCCTGGAAGGTCCTGGTATCAAACTTTTGCATATTATTTCGTGCTGGATACGTGTGGATTTAAAGGAAGGGATCAGTATACCCGCTGGCTGGAAGATATACAGTAc                                >  1:516923/1‑251 (MQ=255)                                                                                                                                                                                                  aCATAAGCAGCCGCCATCGGTTCTGGCCCCAGCGCGCGCAGACAGGTCATTGGGTGAGAGGTTCCCGCGCCGACTTCCATGTCCAATGGTTGAACAATGGTGCAGCCCTGGCGAGCCCAGTAATCCTGTAAGGTCAGGATCAAGCCCTGGAAGGTCCTGGTATCAAACTTTTGCATATTATTTCGTGCTGGATACGTGTGGATTTAAAGGAAGGGATCAGTATACCCGCTGGATGGAAGATATACAGTACg                               <  1:295702/251‑1 (MQ=255)                                                                                                                                                                                                        gcagcCGCCATCGGTTCTGGCCCCAGCGCGCGCAGACAGGTCATTGGGTGAGAGGTTCCCGCGCCGACTTCCATGTCCAATGGTTGAACAATGGTGCAGCCCTGGCGAGCCCAGTAATCCTGTAAGGTCAGGATCAAGCCCTGGAAGGTCCTGGTATCAAACTTTTGCATATTATTTCGTGCTGGATACGTGTGGCTTTAAAGGAAGGGATCAGTATCCCCGCTGGATGGAGGTTATACAGTACGAAACgg                         >  2:314395/1‑251 (MQ=255)                                                                                                                                                                                                                   cGGTTCTGGCCCCAGCGCGCGCAGACAGGTCATTGGGTGAGAGGTTCCCGCGCCGACTTCCATGTCCAATGGTTGAACAATGGTGCAGCCCTGGCGAGCCCAGTAATCCTGTAAGGTCAGGATCAAGCCCTGGAAGGTCCTGGTATCAAACTTTTGCATATTATTTCGTGCTGGATACGTGTGGATTTAAAGGAAGGGATCAGTATACCCGCTGGATGGAAGATATACAGTACGAAACGGGAAAAAGCAgg              <  1:31429/251‑1 (MQ=255)                                                                                                                                                                                                                          ggCCCCAGCGCGCGCAGACAGGTCATTGGGTGAGAGGTTCCCGCGCCGACTTCCATGTCCAATGGTTGAACAATGGTGCAGCCCTGGCGAGCCCAGTAATCCTGTAAGGTCAGGATCAAGCCCTGGAAGGTCCTGGTATCAAACTTTTGCATATTATTTCGTGCTGGATACGTGTGGATTTAAAGGAAGGGATCAGTATACCCGCTGGATGGAAGATATACAGTACGAAACGGGAAAAAGCAGGGCtt          <  2:246972/248‑1 (MQ=255)                                                                                                                                                                                                                             cccAGCGCGCGCAGACAGGTCATTGGGTGAGAGGTTCCCGCGCCGACTTCCATGTCCAATGGTTGAACAATGGTGCAGCCCTGGCGAGCCCAGTAATCCTGTAAGGTCAGGATCAAGCCCTGGAAGGTCCTGGTATCAAACTTTTGCATATTATTTCGTGCTGGATACGTGTGGATTTAAAGGAAGGGATCAGTATACCCGCTGGATGGAAGATATACAGTACGAAACGGGATAAAGCAGGGCTTAACGCa    <  2:404339/251‑1 (MQ=255)                                                                                                                                                                                                                               cAGCGCGCGCAGACAGGTCATTGGGTGAGAGGTTCCCGCGCCGACTTCCATGTCCAATGGTTGAACAATGGTGCAGCCCTGGCGAGCCCAGTAATCCTGTAAGGTCAGGATCAAGCCCTGGAAGGTCCTGGTATCAAACTTTTGCATATTATTTCGTGCTGGATACGTGTGGATTTAAAGGAAGGGATCAGTATACCCGCTGGATGGAAGATATACAGTACGAAACGGGAAAAAGCAGGGCTTAACGCATg  <  1:36237/251‑1 (MQ=255)                                                                                                                                                                                                                                   |                                                                                                                                                                                                                                                       GTTTTCCCAGTTATCTTCCACGAAACGGATGTCGTGAATAGTCGGGTCCATGCCCAGCTCTTTCAGAGAACCGAGGTACAGCTCCTGAATATTGTCCGGCGATGGCTTAATGACCACCTGGAACTGATAGTAGTGCTGTAAACGGTTGGGGTTTTCGCCGTAGCGACCATCGGTCGGGCGACGAGAAGGCTGAACATAAGCAGCCGCCATCGGTTCTGGCCCCAGCTCGCGCAGACAGGTCATTGGGTGAGAGGTTCCCGCGCCGACTTCCATGTCCAATGGTTGAACAATGGTGCAGCCCTGGCGAGCCCAGTAATCCTGTAAGGTCAGGATCAAGCCCTGGAAGGTCCTGGTATCAAACTTTTGCATATTATTTCGTGCTGGATACGTGTGGATTTAAAGGAAGGGATCAGTATACCCGCTGGATGGAAGATATACAGTACGAAACGGGAAAAAGCAGGGCTTAACGCATG  >  NC_000913/3724950‑3725422

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