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	97.7%	149.9 / NA	44	E48A (GAG→GCG)	glyQ	glycine‑‑tRNA ligase subunit alpha
Reads supporting (aligned to +/- strand):  ref base T (0/1);  new base G (18/25);  total (18/26)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.99e-01

AGGCACCCAGCGTCGGGTTTTCCCAGTTATCTTCCACGAAACGGATGTCGTGAATAGTCGGGTCCATGCCCAGCTCTTTCAGAGAACCGAGGTACAGCTCCTGAATATTGTCCGGCGATGGCTTAATGACCACCTGGAACTGATAGTAGTGCTGTAAACGGTTGGGGTTTTCGCCGTAGCGACCATCGGTCGGGCGACGAGAAGGCTGAACATAAGCAGCCGCCATCGGTTCTGGCCCCAGCTCGCGCAGACAGGTCATTGGGTGAGAGGTTCCCGCGCCGACTTCCATGTCCAATGGTTGAACAATGGTGCAGCCCTGGCGAGCCCAGTAATCCTGTAAGGTCAGGATCAAGCCCTGGAAGGTCCTGGTATCAAACTTTTGCATATTATTTCGTGCTGGATACGTGTGGATTTAAAGGAAGGGATCAGTATACCCGCTGGATGGAAGATATACAGTACGAAACGGGAAAAAGCAGGGCTTAACG  >  NC_000913/3724934‑3725418                                                                                                                                                                                                                                                   |                                                                                                                                                                                                                                                   aGGCACCCAGCGTCGGGTTTTCCCAGTTATCTTCCACGAAACGGATGTCGTGAATAGTCGGGTCCATGCCCAGCTCTTTCAGAGAACCGAGGTACAGCTCCTGAATATTGTCCGGCGATGGCTTAAAGACCACCTGGAACTGATAGTAGTGCTGTAAACGGTTGGGGTTTTCGCCGTAGCGACCATCGGTCGGGCGACGAGAAGGCTGAACATAAGCAGCCGCCATCGGTTCTGGCCCCAGCGCGCGCAGa                                                                                                                                                                                                                                            <  2:245230/251‑1 (MQ=255)         aGCGTCGGGTTTTCCCAGTTATCTTCCACGAAACGGATGTCGTGAATAGTCGGGTCCATGCCCAGCTCTTTCAGAGAACCGAGGTACAGCTCCTGAATATTGTCCGGCGATGGCTTAATGACCACCTGGAACTGATAGTAGTGCTGTAAACGGTTGGGGTTTTCGCCGTAGCGACCATCGGTCGGGCGACGAGAAGGCGGAACATAAGCCGCCGCCATCGGTTCTGGCCCCAGCGCGCGCAGaacgggcat                                                                                                                                                                                                                                    >  2:140882/1‑243 (MQ=255)              cGGGTTTTCCCAGTTATCTTCCACGAAACGGTTGTCGTGAATAGTCGGGTCCATGCCCAGCTCTTTCAGAGAACCGAGGTACAGCTCCTGCATATTGTCCGGCGATGGCTTAATGACCCCCTGGAACTGATAGTAGTGCTGTAAACGGTTGGGGTTTTCGCCGTAGCGACCATCGGTCGGGCGACGAGAAGGCTGAACATAAGCAGCCGCCATCGGTTCTGGCCCCAGCGCGCGCAGACAGGTCATTGGGt                                                                                                                                                                                                                               >  2:266921/1‑251 (MQ=255)               gggTTTTCCCAGTTATCTTCCACGAAACGGATGTCGTGAATAGTCGGGTCCATGCCCAGCTCTTTCAGAGAACCGAGGTACAGCTCCTGAATATTGTCCGGCGATGGCTTAATGACCACCTGGAACTGATAGTAGTGCTGTAAACGGTTGGGGTTTTCGCCGTAGCGACCATCGGTCGGGCGACGAGAAGGCTGAACATAAGCAGCCGCCATCGGTTCTGGCCCCAGCGCGCGCAGACAGGTCATTGGGTg                                                                                                                                                                                                                              >  1:62476/1‑251 (MQ=255)                 gTTTTCCCAGTTATCTTCCACGAAACGGATGTCGTGAATAGTCGGGTCCATGCCCAGCTCTTTCAGAGAACCGAGGTACAGCTCCTGAATATTGTCCGGCGATGGCTTAATGACCACCTGGAACTGATAGTAGTGCTGTAAACGGTTGGGGTTTTCGCCGTAGCGACCATCGGTCGGGCGACGAGAAGGCTGAACATAAGCAGCCGCCATCGGTTCTGGCCCCAGCGCGCGCAGACAGGTCATTGGGTgag                                                                                                                                                                                                                            <  1:140882/251‑1 (MQ=255)                      cccAGTTATCTTCCACGAAACGGATGTCGTGAATAGTCGGGTCCATGCCCAGCTCTTTCAGAGAACCGAGGTACAGCTCCTGAATATTGTCCGGCGATGGCTTAATGACCACCTGGAACTGATAGTAGTGCTGTAAACGGTTGGGGTTTTCGCCGTAGCGACCATCGGTCGGGCGACGAGAAGGCTGAACATAAGCAGCCGCCATCGGTTCTGGCCCCAGCGCGCGCAGACAGGTCATTGGGTGAGAGGtt                                                                                                                                                                                                                       >  1:10798/1‑251 (MQ=255)                          gTTATCTTCCACGAAACGGATGTCGTGAATAGTCGGGTCCATGCCCAGCTCTTTCAGAGAACCGAGGTACAGCTCCTGAATATTGTCCGGCGATGGCTTAATGACCACCTGGAACTGATAGTAGTGCTGTAAACGGTTGGGGTTTTCGCCGTAGCGACCATCGGTCGGGCGACGAGAAGGCTGAACATAAGCAGCCGCCATCGGTTCTGGCCCCAGCGCGCGCAGACAGGTCATTGGGTGAGAGGTTCCcg                                                                                                                                                                                                                   >  1:82438/1‑251 (MQ=255)                              tattCCACGAAACGGATGTCGTGAATAGTCGGGTCCATGCCCAGCTCTTTCAGAGAACCGAGGTACAGCTCCTGAATATTGTCCGGCGATGGCTTAATGACCACCTGGAACTGATAGTAGTGCTGTAAACGGTTGGGGTTTTCGCCGTAGCGACCATCGGTCGGGCGACGAGAAGGCTGAACATAAGCAGCCGCCATCGGTTCTGGCCCCAGCGCGCGCAGACAGGTCATTGGGTGAGAGGTTCCCGCGcc                                                                                                                                                                                                               >  1:84496/3‑251 (MQ=255)                                 tCCACGAAACGGATGTCGTGAATAGTCGGGTCCATGCCCAGCTCTTTCAGAGAACCGAGGTACAGCTCCTGAATATTGTCCGGCGATGGCTTAATGACCACCTGGAACTGATAGTAGTGCTGTAAACGGTTGGGGTTTTCGCCGTAGCGACCATCGGTCGGGCGACGAGAAGGCTGAACATAAGCAGCCGCCATCGGTTCTGGCCCCAGCGCGCGCAGACAGGTCATTGGGTGAGAGGTTCCCGCGcc                                                                                                                                                                                                               <  1:12835/248‑1 (MQ=255)                                        aaCGGATGTCGTGAATAGTCGGGTCCATGCCCAGCTCTTTCAGAGAACCGAGGTACAGCTCCTGAATATTGTCCGGCGATGGCTTAATGACCACCTGGAACTGATAGTAGTGCTGTAAACGGTTGGGGTTTTCGCCGTAGCGACCATCGGTCGGGCGACGAGAAGGCTGAACATAAGCAGCCGCCATCGGTTCTGGCCCCAGCGCGCGCAGACAGGTCATTGGGTGAGAGGTTCCCGCGCCGACTTCCATg                                                                                                                                                                                                     >  2:317617/1‑251 (MQ=255)                                           ggATGTCGTGAATAGTCGGGTCCATGCCCAGCTCTTTCAGAGAACCGAGGTACAGCTCCTGAATATTGTCCGGCGATGGCTTAATGACCACCTGGAACTGATAGTAGTGCTGTAAACGGTTGGGGTTTTCGCCGTAGCGACCATCGGTCGGGCGACGAGAAGGCTGAACATAAGCAGCCGCCATCGGTTCTGGCCCCAGCGCGCGCAGACAGGTCATTGGGTGAGAGGTTCCCGCGCCGACTTCCATGTcc                                                                                                                                                                                                  <  1:45184/251‑1 (MQ=255)                                               gTCGTGAATAGTCGGGTCCATGCCCAGCTCTTTCAGAGAACCGAGGTACAGCTCCTGAATATTGTCCGGCGATGGCTTAATGACCACCTGGAACTGATAGTAGTGCTGTAAACGGTTGGGGTTTTCGCCGTAGCGACCATCGGTCGGGCGACGAGAAGGCTGAACATAAGCAGCCGCCATCGGTTCTGGCCCCAGCGCGCGCAGACAGGTCATTGGGTGAGAGGTTCCCGCGCCGACTTCCATGTCCAATg                                                                                                                                                                                              >  2:261053/1‑251 (MQ=255)                                                    gAATAGTCGGGTCCATGCCCAGCTCTTTCAGAGAACCGAGGTACAGCTCCTGAATATTGTCCGGCGATGGCTTAATGACCACCTGGAACTGATAGTAGTGCTGTAAACGGTTGGGGTTTTCGCCGTAGCGACCATCGGTCGGGCGACGAGAAGGCTGAACATAAGCAGCCGCCATCGGTTCTGGCCCCAGCGCGCGCAGACAGGTCATTGGGTGAGAGGTTCCCGCGCCGACTTCCATGTCCAATGGTTg                                                                                                                                                                                          <  1:271052/250‑1 (MQ=255)                                                    gAATAGTCGGGTCCATGCCCAGCTCTTTCAGAGAACCGAGGTACAGCTCCTGAATATTGTCCGGCGATGGCTTAATGACCACCTGGAACTGATAGTAGTGCTGTAAACGGTTGGGGTTTTCGCCGTAGCGACCATCGGTCGGGCGACGAGAAGGCTGAACATAAGCAGCCGCCATCGGTTCTGGCCCCAGCGCGCGCAGACAGGTCATTGGGTGAGAGGTTCCCGCGCCGACTTCCATGTCCAATGGTTg                                                                                                                                                                                          >  1:48931/1‑250 (MQ=255)                                                      aTAGTCGGGTCCATGCCCAGCTCTTTCAGAGAACCGAGGTACAGCTCCTGAATATTGTCCGGCGATGGCTTAATGACCACCTGGAACTGATAGTAGTGCTGTAAACGGTTGGGGTTTTCGCCGTAGCGACCATCGGTCGGGCGACGAGAAGGCTGAACATAAGCAGCCGCCATCGGTTCTGGCCCCAGCGCGCGCAGACAGGTCATTGGGTGAGAGGTTCCCGCGCCGACTTCCATGTCCAATGGTTGAAc                                                                                                                                                                                       <  1:215280/251‑1 (MQ=255)                                                         gTCGGGTCCATGCCCAGCTCTTTCAGAGAACCGAGGTACAGCTCCTGAATATTGTCCGGCGATGGCTTAATGACCTCCTGGAACTGATAGTAGTGCTGTAAACGGTTGGGGTTTTCGCCGTAGCGACCATCGGTCGGGCGACGAGAAGGCTGAACATAAGCAGCCGCCATCGGTTCTGGCCCCAGCGCGCGCAGACAGGTCATTGGGTGAGAGGTTCCCGCGCCGACTTCCATGTCCAATGGTTGAACAAt                                                                                                                                                                                    <  1:171204/251‑1 (MQ=255)                                                                       cAGCTCTTTCAGAGAACCGAGGTACAGCTCCTGAATATTGTCCGGCGATGGCTTAATGACCACCTGGAACTGATAGTAGTGCTGTAAACGGTTGGGGTTTTCGCCGTAGCGACCATCGGTCGGGCGACGAGAAGGCTGAACATAAGCAGCCGCCATCGGTTCTGGCCCCAGCGCGCGCAGACAGGTCATTGGGTGAGAGGTTCCCGCGCCGACTTCCATGTCCAATGGTTGAACAATGGTGCAGCCCTGGc                                                                                                                                                                      <  1:127887/251‑1 (MQ=255)                                                                               tCAGAGAACCGAGGTACAGCTCCTGAATATTGTCCGGCGATGGCTTAATGACCACCTGGAACTGATAGTAGTGCTGTAAACGGTTGGGGTTTTCGCCGTAGCGACCATCGGTCGGGCGACGAGAAGGCTGAACATAAGCAGCCGCCATCGGTTCTGGCCCCAGCGCGCGCAGACAGGTCATTGGGTGAGAGGTTCCCGCGCCGACTTCCATGTCCAATGGTTGAACAATGGTGCAGCCCTGGCGAGCCCa                                                                                                                                                               <  1:261053/250‑1 (MQ=255)                                                                                              aCAGCTCCTGAATATTGTCCGGCGATGGCTTAATGACCACCTGGAACTGATAGTAGTGCTGTAAACGGTTGGGGTTTTCGCCGTAGCGACCATCGGTCGGGCGACGAGAAGGCTGAACATAAGCAGCCGCCATCGGTTCTGGCCCCAGCGCGCGCAGACAGGTCATTGGGTGAGAGGTTCCCGCGCCGACTTCCATGTCCAACGGTTGAACAATGGTGCAGCCCTGGCGAGCCCAGTAATCCTGTAAGGTc                                                                                                                                               >  2:64802/1‑251 (MQ=255)                                                                                              aCAGCTCCTGAATATTGTCCGGCGATGGCTTAATGACCACCTGGAACTGATAGTAGTGCTGTAAACGGTTGGGGTTTTCGCCGTAGCGACCATCGGTCGGGCGACGAGAAGGCTGAACATAAGCAGCCGCCAGCGGTTCTGGCCCCAGCGCGCGCAGACAGGTCATTGGGTGAGAGGTGCCCGCGCCGACTTCCATGTCCAATGGTTGAACAATGGTGCAGCCCTGGCGAGCCCAGTAATCCTGTAAGGTc                                                                                                                                               <  1:199387/251‑1 (MQ=255)                                                                                                       gAATATTGTCCGGCGATGGCTTAATGACCACCTGGAACTGATAGTAGTGCTGTAAACGGTTGGGGTTTTCGCCGTAGCGACCATCGGTCGGGCGACGAGAAGGCTGAACATAAGCAGCCGCCATCGGTTCTGGCCCCAGCGCGCGCAGACAGGTCATTGGGTGAGAGGTTCCCGCGCCGACTTCCATGTCCAATGGTTGAACAATGGTGCAGCCCTGGCGAGCCCAGTAATCCTGTAAGGTCAGGATCAAg                                                                                                                                      <  1:6307/251‑1 (MQ=255)                                                                                                              gTCCGGCGAGGGCTTAATGACCACCTGGAACTGATAGTAGTGCTGTAAACGGTTGGGGTTTTCGCCGTAGCGACCATCGGTCGGGCGACGAGAAGGCTTAACATAAGCAGCCGCCATCGGTTCTGGCCCCAGCTCGCGCAGACATGTCATTGGGTGAGAGGTTCCCGCGCCGACTTCCATGTCCAATGGTTGAACAATGGTGCAGCCCTGGCGAGCCCAGTAATCCTGTAAGGTCAGGATCAAGCCCTgg                                                                                                                                <  1:335689/250‑1 (MQ=255)                                                                                                                 cGGCGATGGCTTAATGACCACCTGGAACTGATAGTAGTGCTGTAAACGGTTGGGGTTTTCGCCGTAGCGACCATCGGTCGGGCGACGAGAAGGCTGAACATAAGCAGCCGCCATCGGTTCTGGCCCCAGCGCGCGCAGACAGGTCATTGGGTGAGAGGTTCCCGCGCCGACTTCCATGTCCAATGGTTGAACAATGGTGCAGCCCTGGCGAGCCCAGTAATCCTGTAAGGTCAGGATCAAGCCCTGGAAgg                                                                                                                            >  1:155392/1‑251 (MQ=255)                                                                                                                      aTGGCTTAATGACCGCCTGGAACTGATAGTAGTGCTGTAAACGGTTGGGGTTTTCGCCGTAGCGACCAGCGGTCGGGCGACGAGAAGGCTGAACATAAGCAGCCGCCATCGGTTCTGGCCCCCGCGCGCGCAGACAGGTCATTGGGTGAGAGGTTCCCGCGCCGACTTCCATGTCCAATGGTTGAACAATGGTGCAGCCCTGGCGAGCCCAGTAATCCTGTAAGGTCTGGATCAAGCCCTGTAAGGTCCTg                                                                                                                       <  2:212632/251‑1 (MQ=255)                                                                                                                                      cTGGAACTGATAGTAGTGCTGTAAACGGTTGGGGTTTTCGCCGTAGCGACCATCGGTCGGGCGACGAGAAGGCTGAACATAAGCAGCCGCCATCGGTTCTGGCCCCAGCGCGCGCAGACAGGTCATTGGGTGAGAGGTTCCCGCGCCGACTTCCATGTCCAATGGTTGAACAATGGTGCAGCCCTGGCGAGCCCAGTAATCCTGTAAGGTCAGGATCAAGCCCTGGAAGGTCCTGGTATCAAACTTTTGc                                                                                                        >  1:88682/1‑250 (MQ=255)                                                                                                                                          aaCTGATAGTAGTGCTGTAAACGGTTGGGGTTTTCGCCGTAGCGACCATCGGTCGGGCGACGAGAAGGCTGAACATAAGCAGCCGCCATCGGTTCTGGCCCCAGCGCGCGCAGACAGGTCATTGGGTGAGAGGTTCCCGCGCCGACTTCCATGTCCAATGGTTGAACAATGGTGTAGCCCTGGCGATCCCAGTAATCCTGTAAGGTCAGGATCAAGCCATGGAAGGTCCTGGTATCAAACGTTTGCAtatt                                                                                                   <  2:185275/251‑1 (MQ=255)                                                                                                                                              gATAGTAGTGCTGTAAACGGTTGGGGTTTTCGCCGTAGCGACCATCGGTCGGGCGACGAGAAGGCTGAACATAAGCAGCCGCCATCGGTTCTGGCCCCAGCGCGCGCAGACAGGTCATTGGGTGAGAGGTTCCCGCGCCGACTTCCATGTCCAATGGTTGAACAATGGTGCAGCCCTGGCGAGCCCAGTAATCCTGTAAGGTCAGGATCAAGCCCTGGAAGGTCCTGGTATCAAACTTTTGCAtattat                                                                                                 <  2:48931/249‑1 (MQ=255)                                                                                                                                              gATAGTAGTGCTGTAAACGGTTGGGGTTTTCGCCGTAGCGACCATCGGTCGGGCGACGAGAAGGCTGAACATAAGCAGCCGCCATCGGTTCTGGCCCCAGCGCGCGCAGACAGGTCATTGGGTGAGAGGTTCCCGCGCCGACTTCCATGTCCAATGGTTGAACAATGGTGCAGCCCTGGCGAGCCCAGTAATCCTGTAAGGTCAGGATCAAGCCCTGGAAGGTCCTGGTATCAAACTTTTGCAtattat                                                                                                 <  1:322220/249‑1 (MQ=255)                                                                                                                                              gATAGTAGTGCTGTAAACGGTTGGGGTTTTCGCCGTAGCGACCATCGGTCGGGCGACGAGAAGGCTGAACATAAGCAGCCGCCATCGGTTCTGGCCCCAGCGCGCGCAGACAGGTCATTGGGTGAGAGGTTCCCGCGCCGACTTCCATGTCCAATGGTTGAACAATGGTGCAGCCCTGGCGAGCCCAGTAATCCTGTAAGGTCAGGATCAAGCCCTGGAAGGTCCTGGTATCAAACTTTTGCATattatt                                                                                                <  1:269798/250‑1 (MQ=255)                                                                                                                                              gATAGTAGTGCGGTAAACGGTTGGGGTTTTCGCCGTAGCGACCATCGGTCGGGCGACGAGAAGGCTGAACATAAGCAGCCGCCATCGGTTCTGGCCCCAGCGCGCGCAGACAGGTCATTGGGTGAGAGGTTCCCGCGCCGACTTCCATGTCCAATGGTTGAACAATGGTGCAGCCCTGGCGAGCCCAGTAATCCTGTAAGGTCAGGATCAAGCCCTGGAAGGTCCTGGTATCAAACTTTTGCATATTAttt                                                                                               <  2:10798/251‑1 (MQ=255)                                                                                                                                                 agtagtGCTGTAAACGGTTGGGGTTTTCGCCGTAGCGACCATCGGTCGGGCGACGAGAAGGCTGAACATAAGCAGCCGCCATCGGTTCTGGCCCCAGCGCGCGCAGACAGGTCATTGGGTGAGAGGTTCCCGCGCCGACTTCCATGTCCAATGGTTGAACAATGGTGCAGCCCTGGCGAGCCCAGTAATCCTGTAAGGTCAGGATCAAg                                                                                                                                      >  1:64485/1‑209 (MQ=255)                                                                                                                                                 agtagtGCTGTAAACGGTTGGGGTTTTCGCCGTAGCGACCATCGGTCGGGCGACGAGAAGGCTGAACATAAGCAGCCGCCATCGGTTCTGGCCCCAGCGCGCGCAGACAGGTCATTGGGTGAGAGGTTCCCGCGCCGACTTCCATGTCCAATGGTTGAACAATGGTGCAGCCCTGGCGAGCCCAGTAATCCTGTAAGGTCAGGATCAAg                                                                                                                                      <  2:64485/209‑1 (MQ=255)                                                                                                                                                  gtagtGCTGTAAACGGTTGGGGTTTTCGCCGTAGCGACCATCGGTCGGGCGACGAGAAGGCTGAACATAAGCAGCCGCCATCGGTTCTGGCCCCAGCGCGCGCAGACAGGTCATTGGGTGAGAGGTTCCCGCGCCGACTTCCATGTCCAATGGTTGAACAATGGTGCAGCCCTGGCGAGCCCAGTAATCCTGTAAGGTCAGGATCAAGCCCTGGAAGGTCCTGGTATCAAACTTTTGCATATTATTTCGTg                                                                                           <  1:54757/251‑1 (MQ=255)                                                                                                                                                     gtGCTGTAAACGGTTGGGGTTTTCGCCGTAGCGACCATCGGTCGGGCGACGAGAAGGCTGAACATAAGCAGCCGCCATCGGTTCTGGCCCCAGCGCGCGCAGACAGGTCATTGGGTGAGAGGTTCCCGCGCCGACTTCCATGTCCAATGGTTGAACAATGGTGCAGCCCTGGCGAGCCCAGTAATCCTGTAAGGTCAGGATCAAGCCCTGGAAGGTCCTGGTATCAAACTTTTGCATATTATTTCGTGCTg                                                                                        >  1:45155/1‑251 (MQ=255)                                                                                                                                                        cTGTAAACGGTTGGGGTTTTCGCCGTAGCGACCATCGGTCGGGCGACGAGAAGGCTGAACATAAGCAGCCGCCATCGGTTCTGGCCCCAGCGCGCGCAGACAGGTCATTGGGTGAGAGGTTCCCGCGCCGACTTCCATGTCCAATGGTTGAACAATGGTGCAGCCCTGGCGAGCCCAGTAATCCTGTAAGGTCAGGATCAAGCCCTGGAAGGTCCTGGTATCAAACTTTTGCATATTATTTCGTGCTGGAt                                                                                     >  2:157976/1‑251 (MQ=255)                                                                                                                                                                            cgaggaagggactaTCGGTCGGGCGACGAGAAGGCTGAACATAAGGAGCCGCCATCGGTTCCGGCACCAGCGCGCGGAGACAGGTCATTGGGTGAGAGGTTCCCGCGCCGACTTTGATGTCCAATGGTTGAACAATGGTGCAGCCATGGCGAGCCCAGTAATCATGTAAGGTCAGGATCAAGCCCTGGAAGGTACTGGTATCAAACTTTTGCATATTATTTCGTGCTGGATACGGGTGGATTTAaagga                                                                   >  2:18498/14‑249 (MQ=255)                                                                                                                                                                             gCCGTAGCGACCATCGGTCGGGCGACGAGAAGGCTGAACATAAGCAGCCGCCATCGGTTCTGGCCCCAGCGCGCGCAGACAGGTCATTGGGTGAGAGGTTCCCGCGCCGACTTCCATGTCCAATGGTTGAACAATGGTGCAGCCCTGGCGAGCCCAGTAATCCTGTAAGGTCAGGATCAAGCCCTGGAAGGTCCTGGTATCAAACTTTTGCATATTATTTCGTGCTGGATACGTGTGGATTTAaaggaagg                                                                >  2:145995/1‑251 (MQ=255)                                                                                                                                                                                     gACCATCGGTCGGGCGACGAGAAGGCTGAACATAAGCAGCCGCCATCGGTTCTGGCCCCAGCGCGCGCAGACAGGTCATTGGGTGAGAGGTTCCCGCGCCGACTTCCATGTCCAATGGTTGAACAATGGTGCAGCCCTGGCGAGCCCAGTAATCCTGTAAGGTCAGGATCAAGCCCTGGAAGGTCCTGGTATCAAACTTTTGCATATTATTTCGTGCTGGATACGTGTGGATTTAAAGGAAGGGATCAGta                                                        <  1:301841/251‑1 (MQ=255)                                                                                                                                                                                           cggtcggGCGACGAGAAGGCTGAACATAAGCAGCCGCCATCGGTTCTGGCCCCAGCGCGCGCAGACAGGTCATTGGGTGAGAGGTTCCCGCGCCGACTTCCATGTCCAATGGTTGAACAATGGTGCAGCCCTGGCGAGCCCAGTAATCCTGTAAGGTCAGGATCAAGCCCTGGAAGTTCCTGTTATCAAACTTTTGCATATTATTTCGTGCTGGATACGTGTGGATTTAAAGGAAGTGATCAGTATACCCg                                                  <  1:221076/251‑1 (MQ=255)                                                                                                                                                                                                         gaAGGCTGAACATAAGCAGCCGCCATCGGTTCTGGCCCCAGCGCGCGCAGACAGGTCATTGGGTGAGAGGTTCCCGCGCCGACTTCCATGTCCAATGGTTGAACAATGGTGCAGCCCTGGCGAGCCCAGTAATCCTGTAAGGTCAGGATCAAGCCCTGGAAGGTCCTGGTATCAAACTTTTGCATATTATTTCGTGCTGGATACGTGTGGATTTAAAGGAAGGGATCAGTATACCCGCTGGATGGAAGata                                    <  1:157976/251‑1 (MQ=255)                                                                                                                                                                                                                       agcagcCGCCATCGGTTCTGGCCCCAGCGCGCGCAGACAGGTCATTGGGTGAGAGGTTCCCGCGCCGACTTCCATGTCCAATGGTTGAACAATGGTGCAGCCCTGGCGAGCCCAGTAATCCTGTAAGGTCAGGATCAAGCCCTGGAAGGTCCTGGTATCAAACTTTTGCATATTATTTCGTGCTGGATACGTGTGGATTTAAAGGAAGGGATCAGTATACCCGCTGGATGGAAGATATACAGTACGAAACg                      >  2:316901/1‑251 (MQ=255)                                                                                                                                                                                                                            ccgccATCGGTTCTGGCCCCAGCGCGCGCAGACAGGTCATTGGGTGAGAGGTTCCCGCGCCGACTTCCATGTCCAATGGTTGAACAATGGTGCAGCCCTGGCGAGCCCAGTAATCCTGTAAGGTCAGGATCAAGCCCTGGAAGGTCCTGGTATCAAACTTTTGCATATTATTTCGTGCTGGATACGTGTGGATTTAAAGGAAGGGATCAGTATACCCGCTGGATGGAAGATATACAGTACGAAACGGGaaa                 <  2:53655/251‑1 (MQ=255)                                                                                                                                                                                                                                          ggCCCCAGCGCGCGCAGACAGGTCATTGGGTGAGAGGTTCCCGCGCCGACTTCCATGTCCAATGGTTGAACAATGGTGCAGCCCTGGCGAGCCCAGTAATCCTGTAAGGTCAGGAGCAAGCCCTGGAAGGTCCTGGTATCAAACTTTTGCATATTATTTCGTGCTGGATACGTGTGGATTTAAAGGAAGGGATCAGTATACCCGCTGGATGGAAGATATACAGTACGAAACGGGAAAAAGCAGGGCTTAAc   <  1:41437/251‑1 (MQ=255)                                                                                                                                                                                                                                           gCCCCAGCGCGCGCAGACAGGTCATTGGGTGAGAGGTTCCCGCGCCGACTTCCATGTCCAATGGTTGAACAATGGTGCAGCCCTGGCGAGCCCAGTAATCCTGTAAGGTCAGGAGCAAGCCCTGGAAGGTCCTGGTATCAAACTTTTGCATATTATTGCGTGCTGGATACGTGTGGATTTAAAGGAAGGGATCAGTATACCCGCTGGATGGAAGATATACAGTACGAAACGGGAAAAAGCAGGGCTTAACg  <  1:266921/251‑1 (MQ=255)                                                                                                                                                                                                                                                   |                                                                                                                                                                                                                                                   AGGCACCCAGCGTCGGGTTTTCCCAGTTATCTTCCACGAAACGGATGTCGTGAATAGTCGGGTCCATGCCCAGCTCTTTCAGAGAACCGAGGTACAGCTCCTGAATATTGTCCGGCGATGGCTTAATGACCACCTGGAACTGATAGTAGTGCTGTAAACGGTTGGGGTTTTCGCCGTAGCGACCATCGGTCGGGCGACGAGAAGGCTGAACATAAGCAGCCGCCATCGGTTCTGGCCCCAGCTCGCGCAGACAGGTCATTGGGTGAGAGGTTCCCGCGCCGACTTCCATGTCCAATGGTTGAACAATGGTGCAGCCCTGGCGAGCCCAGTAATCCTGTAAGGTCAGGATCAAGCCCTGGAAGGTCCTGGTATCAAACTTTTGCATATTATTTCGTGCTGGATACGTGTGGATTTAAAGGAAGGGATCAGTATACCCGCTGGATGGAAGATATACAGTACGAAACGGGAAAAAGCAGGGCTTAACG  >  NC_000913/3724934‑3725418

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