Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	NC_000913	1,716,110	G→A	S371S (TCC→TCT)	tyrS ←	tyrosine‑‑tRNA ligase

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	1,716,110	0	G	A	98.2%	166.1 / NA	54	S371S (TCC→TCT)	tyrS	tyrosine‑‑tRNA ligase
Reads supporting (aligned to +/- strand):  ref base G (0/0);  major base A (31/22);  minor base C (1/0);  total (32/22)
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

GCGAGAATATTTTTCATCATGTTCCCTGTATAAAAACCAGGGGGAGTGATTTCTCACTCCCCCTTTCCACTTAATGCATTATTTCCAGCAAATCAGACAGTAATTCTTTTTACCGCGACGCAGTAAGGTAAAACGACCAAACAGACGATCTTCTTCTTTAAAGAAGTATTCAGGATCGGACTGTTTTTCACCGTTAATGGTGATGGCATTGGAGGCGATAGTTTTACGTGCCTGACCACGGGAAGGTTGCAGTTCAGAATCGACCAGTGCCTGCATCAGGTCTGCGCCCTTTTCCATCTCAACCATCGGTACGCCGTCCTGCGCCAGCTGTTCGAAGTCCGCTTCACTCAGCGCACTCAAAGAACCGCTGAACAGGCATTCGGTAATACGTTTTGCCGCCTGTAAACCTTCTTCACCGTGAACCAGACGAGTCACCTGCTCCGCCAGTACATACTGGGCGCGCGGTGCTTTACCGCTGTTTTTATCT  >  NC_000913/1715870‑1716356                                                                                                                                                                                                                                                 |                                                                                                                                                                                                                                                       gCGAGAATATTTTTCATCATGTTCCCTGTATAAAAACCAGGGGGAGTGATTTCTCACTCCCCCTTTCCACTTAATGCATTATTTCCAGCAAATCAGACAGTAATTCTTTTTACCGCGACGCAGTAAGGTAAAACGACCAAACAGACGATCTTCTTCTTTAAAGAAGTATTCAGGATCGGACTGTTTTTCACCGTTAATGGTGATGGCATTGGAGGCGATAGTTTTACGTGCCTGACCACGAGAAGGTTGc                                                                                                                                                                                                                                               >  1:166711/1‑250 (MQ=255)   gagaATATTTTTCATCATGTTCCCTGTATAAAAACCAGGGGGAGTGATTTCTCACTCCCCCTTTCCACTTAATGCATTATTTCCAGCAAATCAGACAGTAATTCTTTTTACCGCGACGCAGTAAGGTAAAACGACCAAACAGACGATCTTCTTCTTTAAAGAAGTATTCAGGATCGGACTGTTTTTCACCGTTAATGGTGATGGCATTGGAGGCGATAGTTTTACGTGCCTGACCACGAGAAGGTTGCAGt                                                                                                                                                                                                                                            >  2:178813/1‑251 (MQ=255)                atcatGTTCCCTGTATAAAAACCAGGGGGAGTGATTTCTCACTCCCCCTTTCCACTTAATGCATTATTTCCAGCAAATCAGACAGTAATTCTTTTTACCGCGACGCAGTAAGGTAAAACGACCAAACAGACGATCTTCTTCTTTAAAGAAGTATTCAGGATCGGACTGTTTTTCACCGTTAATGGTGATGGCATTGGAGGCGATAGTTTTACGTGCCTGACCACGAGAAGGTTGCAGTTCAGAATCGAcc                                                                                                                                                                                                                                >  2:121610/1‑250 (MQ=255)                atcatGTTCCCTGTATAAAAACCAGGGGGAGTGATTTCTCACTCCCCCTTTCCACTTAATGCATTATTTCCAGCAAATCAGACAGTAATTCTTTTTACCGCGACGCAGTAAGGTAAAACGACCAAACAGACGATCTTCTTCTTTAAAGAAGTATTCAGGATCGGACTGTTTTTCACCGTTAATGGTGATGGCATTGGAGGCGATAGTTTTACGTGCCTGACCACGAGAAGGTTGCAGTTCAGAATCGACCa                                                                                                                                                                                                                               <  1:228526/251‑1 (MQ=255)                 tcatGTTCCCTGTATAAAAACCAGGGGGAGTGATTTCTCACTCCCCCTTTCCACTTAATGCATTATTTCCAGCAAATCAGACAGTAATTCTTTTTACCGCGACGCAGTAAGGTAAAACGACCAAACAGACGATCTTCTTCTTTAAAGAAGTATTCAGGATCGGACTGTTTTTCACCGTTAATGGTGATGGCATTGGAGGCGATAGTTTTACGTGCCTGACCACGAGAAGGTTGCAGTTCAGAATCGAcc                                                                                                                                                                                                                                >  1:477605/1‑249 (MQ=255)                   atGTTCCCTGTATAAAAACCAGGGGGAGTGATTTCTCACTCCCCCTTTCCACTTCATGCATTATTTCCAGCAAATCAGACAGTAATTCTTTTTACCGCGACGCAGTAAGGTAAAACGACCAAACAGACGATCTTCTTCTTTAAAGAAGTATTCAGGATCGGACTGTTTTTCACCGTTAATGGTGATGGCATTGGAGGCGATAGTTTTACGTGCCTGACCACGAGAAGGTTGCAGTTCAGAATCGACCAGTg                                                                                                                                                                                                                            >  1:410108/1‑251 (MQ=255)                       tCCCTGTATAAAAACCAGGGGGAGTGATTTCTCACTCCCCCTTTCCACTTAATGCATTATTTCCAGCAAATCAGACAGTAATTCTTTTTACCGCGACGCAGTAAGGTAAAACGACCAAACAGACGATCTTCTTCTTTAAAGAAGTATTCAGGATCGGACTGTTTTTCACCGTTAATGGTGATGGCATTGGAGGCGATAGTTTTACGTGCCTGACCACGAGAAGGTTGCAGTTCAGAATCGACCAGTGCCt                                                                                                                                                                                                                         <  2:477605/250‑1 (MQ=255)                         ccTGTATAAAAACCAGGGGGAGTGATTTCTCACTCCCCCTTTCCACTTAATGCATTATTTCCAGCAAATCAGACAGTAATTCTTTTTACCGCGACGCAGTAAGGTAAAACGACCAAACAGACGATCTTCTTCTTTAAAGAAGTATTCAGGATCGGACTGTTTTTCACCGTTAATGGTGATGGCATTGGAGGCGATAGTTTTACGTGCCTGACCACGAGAAGGTTGCAGTTCAGAATCGACCAGTGCCTGCa                                                                                                                                                                                                                      <  1:209794/251‑1 (MQ=255)                            gTATAAAAACCAGGGGGAGTGATTTCTCACTCCCCCTTTCCACTTAATGCATTATTTCCAGCAAATCAGACAGTAATTCTTTTTACCGCGACGCAGTAAGGTAAAACGACCAAACAGACGATCTTCTTCTTTAAAGAAGTATTCAGGATCGGACTGTTTTTCACCGTTAATGGTGATGGCATTGGAGGCGATAGTTTTACGTGCCTGACCACGAGAAGGTTGCAGTTCAGAATCGACCAGTGCCTGCATc                                                                                                                                                                                                                    <  1:11386/250‑1 (MQ=255)                              ataAAAACCAGGGGGAGTGATTTCTCACTCCCCCTTTCCACTTAATGCATTATTTCCAGCAAATCAGACAGTAATTCTTTTTACCGCGACGCAGTAAGGTAAAACGACCAAACAGACGATCTTCTTCTTTAAAGAAGTATTCAGGATCGGACTGTTTTTCACCGTTAATGGTGATGGCATTGGAGGCGATAGTTTTACGTGCCTGACCACGAGAAGGTTGCAGTTCAGAATCGACCAGTGCCTGCATCAgg                                                                                                                                                                                                                 >  2:531128/1‑251 (MQ=255)                                aaaaaCCAGGGGGAGTGATTTCTCACTCCCCCTTTCCACTTAATGCATTATTTCCAGCAAATCAGACAGTAATTCTTTTTACCGCGACGCAGTAAGGTAAAACGACCAAACAGACGATCTTCTTCTTTAAAGAAGTATTCAGGATCGGACTGTTTTTCACCGTTAATGGTGATGGCATTGGAGGCGATAGTTTTACGTGCCTGACCACGAGAAGGTTGCAGTTCAGAATCGACCAGTGCCTGCATCAGGTc                                                                                                                                                                                                               >  2:265578/1‑251 (MQ=255)                                                       cACTCCCCCTTTCCACTTAATGCATTATTTCCAGCAAATCAGACAGTAATTCTTTTTACCGCGACGCAGTAATGTAAAACGACCAAACAGACGATCTTCTTCTTTAAAGAAGTATTCAGGATCGGACTGTTTTTCACCGTTAATGGTGATGGCATTGGAGGCGATAGTTTTACGTGCCTGACCACGAGAAGGTTGCAGTTCAGAATCGACCAGTGCCTGCATCAGGTCTGCGCCCTTTTCCATCTCAAcc                                                                                                                                                                                         <  2:177381/250‑1 (MQ=255)                                                                         aaTGCATTATTTCCAGCAAATCAGACAGTAATTCTTTTTACCGCGACGCAGTAAGGTAAAACGACCAAACAGACGCTCTTCTTCTTTAAAGAAGTATTCAGGATCGGACTGTTTTTCACCGTTAATGGTGATGGCATTGGAGGCGATAGTTTTACGTGCCTGACCACGAGAAGGTTGCAGTTCAGAATCGACCAGTGCCTGCATCAGGTCTGCGCCCTTTTCCATCTCAACCATCGGTACGCCGTCCTgcg                                                                                                                                                                      >  2:416652/1‑251 (MQ=255)                                                                                 attTCCAGCAAATCAGCCAGTAATACTTTTTACCGAGACGCAGTAACGTAAAACGACCAAACAGACGATCTTCTTCTTTAAAGCAGCATTAAGGATCGAACTGTTTTTCAACGTTAATGGTGATGGCATTGGAGGCGATAGTTTTACGTGCCTGACCACGAGAAGGTTGCAGTTCAGAATCGACCAGTGCCTGCATCAGGTCTGCGCCCTTTTCCATCTCAACCATCGGTACGCCGTCCTGCGCCAGCTGt                                                                                                                                                              >  2:400148/1‑251 (MQ=255)                                                                                             tCAGACAGTAATTCTTTTTACCGCGACGCAGTAAGGTAAAACGACCAAACAGACGATCTTCTTCTTTAAAGAAGTATTCAGGATCGGACTGTTTTTCACCGTTAATGGTGATGGCATTGGAGGCGATAGTTTTACGTGCCTGACCACGAGAAGGTTGCAGTTCAGAATCGACCAGTGCCTGCATCAGGTCTGCGCCCTTTTCCATCTCAACCATCGGTACGCCGTCCTGCGCCAGCTGTTCGAAGTCCGCt                                                                                                                                                  >  1:474917/1‑251 (MQ=255)                                                                                                     tacttCTTTTTACCGCGAAGCAGTAAGGAAAAACGACCAAACAGACGATCTTCTTCTTTAAAGAAGTATTCAGGATCGGACTGTTTTTCACCGTTAATGGTGATGGCATTGGAGGCGATAGTTTTACGTGCCTGACCACGAGAAGGTTGCAGTTCAGAATCGACCAGTGCCTGCATCAGGTCTGCGACCTTTTCCATCTCAACCATCGGTACGCCGTCCTGCGCCAGCTGTTCGAAGTCCGCTTCACTc                                                                                                                                            >  2:277005/4‑249 (MQ=255)                                                                                                          cTTTTTACCGCGACGCAGTAAGGTAAAACGACCAAACAGACGATCTTCTTCTTTAAAGAAGCATTCAGGATCGGACTGTTTTTCACCGTTAATGGTGATGGCATTGGAGGCGATAGTTTTACGTGCCTGACCACGAGAAGGTTGCAGTTCAGAATCGACCAGTGCCTGCATCAGGTCTGCGCCCTTTTCCATCTCAACCATCGGTACGCCGTCCTGCGCCAGCTGTTCGAAGTCCGCTTCACTCAGCGCAc                                                                                                                                     >  1:93591/1‑251 (MQ=255)                                                                                                              ttACCGCGACCCAGTAAGGTAAACCGACCAAACAGACGATCTTCTCCTTTCATGAAGACTGCAGGATAGGACTGTTTTTCTCAGTTCATGGTGATGGCATTGCAGGCGATAGGTTTCCGTGCGTGACAACGCGAGGATCGCAGTTCAGAATCGACCAGTGCCTGCATCATGTCTCCTCCATATTCCGTCTCAACCATAGGCACGCCGCCCTGCGCCAAATGTTCGCAGACCGCTTGACTCAGCGCAATCa                                                                                                                                  >  2:341056/1‑250 (MQ=255)                                                                                                               tACCGCGACGCAGTAAGGTAAAACGACCAAACAGACGATCTTCTTCTTTAAAGAAGTATTCAGGATCGGACTGTTTTTCACCGTTAATGGTGATGGCATTGGAGGCGATAGTTTTACGTGCCTGACCACGAGAAGGTTGCAGTTCAGAATCGACCAGTGCCTGCATCAGGTCTGCGCCCTTTTCCATCTCAACCATCGGTACGCCGTCCTGCGCCAGCTGTTCGAAGTCCGCTTCACTCAGCGCACTCaa                                                                                                                                 <  1:266176/250‑1 (MQ=255)                                                                                                                  cgcgACGCAGTAAGGTAAAACGACCAAACAGACGATCTTCTTCTTTAAAGAAGTATTCAGGATCGGACTGTTTTTCACCGTTAATGGTGATGGCATTGGAGGCGATAGTTTTACGTGCCTGACCACGAGAAGGTTGCAGTTCAGAATCGACCAGTGCCTGCATCAGGTCTGCGCCCTTTTCCATCTCAACCATCGGTACGCCGTCCTGCGCCAGCTGTTCGAAGTCCGCTTCACTCAGCGCACTCAAAGaa                                                                                                                             <  2:200139/251‑1 (MQ=255)                                                                                                                       cGCAGTAAGGTAAAACGACCAAACAGACGATCTTCGTCTTTAAAGAAGTATTCAGGATCGGACTGTTTTTCACCGTTAATGGTGATGGCATTGGAGGCGATAGTTTTACGTGCCTGACCACGAGAAGGTTGCAGTTCAGAATCGACCAGTGCCTGCATCAGGTCTGCGCCCTTTTCCATCTCAACCATCGGTACGCCGTCCTGCGCCAGCTGTTCGAAGTCCGCTTCACTCAGCGCACTCAAAGAACCGCt                                                                                                                        <  2:429609/251‑1 (MQ=255)                                                                                                                             aaGGTAAAACGACCAAACAGACGATCTTCTTCTTTAAAGAAGTATTCAGGATCGGACTGTTTTTCACCGTTAATGGTGATGGCATTGGAGGCGATAGTTTTACGTGCCTGACCACGAGAAGGTTGCAGTTCAGAATCGACCAGTGCCTGCATCAGGTCTGCGCCCTTTTCCATCTCAACCATCGGTACGCCGTCCTGCGCCAGCTGTTCGAAGTCCGCTTCACTCAGCGCACTCAAAGAACCGCTGAAc                                                                                                                    >  2:336308/1‑249 (MQ=255)                                                                                                                                  aaaaCGACCAAACAGACGATCTTCTTCTTTAAAGAAGTATTCAGGATCGGACTGTTTTTCACCGTTAATGGTGATGGCATTGGAGGCGATAGTTTTACGTGCCTGACCACGAGAAGGTTGCAGTTCAGAATCGACCAGTGCCTGCATCAGGTCTGCGCCCTTTTCCATCTCAACCATCGGTACGCCGTCCTGCGCCAGCTGTTCGAAGTCCGCTTCACTCAGCGCACTCAAAGAACCGCTGAACAGGCAtt                                                                                                             <  2:166711/251‑1 (MQ=255)                                                                                                                                    aaCGACCAAACAGACGATCTTCTTCTTTAAAGAAGTATTCAGGATCGGACTGTTTTTCACCGTTAATGGTGATGGCATTGGAGGCGATAGTTTTACGTGCCTGACCACGAGAAGGTTGCAGTTCAGACTCGACCAGTGCCTGCATCAGGTCTGCGCCCTTTTCCATCTCAACCATCGGTACGCCGTCCTGCGCCAGCTGTTCGAAGTCCGCTTCACTCAGCGCACTCAAAGAACCGCTGAACAGGCATTCg                                                                                                           >  2:252764/1‑251 (MQ=255)                                                                                                                                                              ttAAAGAAGTATTCAGGATCGGACTGTTTTTCACCGTTAATGGTGATGGCATTGGAGGCGATAGTTTTACGTGCCTGACCACGAGAAGGTTGCAGTTCAGAATCGACCAGTGCCTGCAt                                                                                                                                                                                                                     >  1:14695/1‑119 (MQ=255)                                                                                                                                                              ttAAAGAAGTATTCAGGATCGGACTGTTTTTCACCGTTAATGGTGATGGCATTGGAGGCGATAGTTTTACGTGCCTGACCACGAGAAGGTTGCAGTTCAGAATCGACCAGTGCCTGCAt                                                                                                                                                                                                                     <  2:14695/119‑1 (MQ=255)                                                                                                                                                               tAAAGAAGTATTCAGGATCGGACTGTTTTTCACCGTTAATGGTGATGGCATTGGAGGCGATAGTTTTACGTGCCTGACCACGAGAAGGTTGCAGTTCAGAATCGACCAGTGCCTGCATCAGGTCTGCGCCCTTTTCCATCTCAACCATCGGTACGCCGTCCTGCGCCAGCTGTTCGAAGTCCGCTTCACTCAGCGCACTCAAAGAACCGCTGAACAGGCATTCGGTAATACGTTTTGCCGCCTGTAAAcc                                                                                 <  1:288514/250‑1 (MQ=255)                                                                                                                                                                aaaGAAGTATTCAGGATCGGACTGTTTTTCACCGTTAATGGTGATGGCATTGGAGGCGATAGTTTTACGTGCCTGACCACGAGAAGGTTGCAGTTCAGAATCGACCAGTGCCTGCATCAGGTCTGCGCCCTTTTCCATCTCAACCATCGGTACGCCGTCCTGCGCCAGCTGTTCGa                                                                                                                                                          <  2:470848/176‑1 (MQ=255)                                                                                                                                                                aaaGAAGCATTCAGGATCGGACTGTTTTTCATCGTTAATTGAGATGGCATTGTAGGCGCTCGGTTTACGTGCCTGACCACGAGATGGTTGCAGTTCAGAATCGACCAGTGCCTGCATCAGGTCGGCGCCCTTTTCCATCCCAACCATCGGTACCCCGTCCTGCGCCAGCTGTTCGAAGCCCGCTTCACCCAGCGCACTCAAAGCACCGCTGAAAAGGCATTCGGTAATACGTTTCGCAGCCCGTAAACctt                                                                               >  2:340488/1‑251 (MQ=255)                                                                                                                                                                  agaagTATTCAGGATCGGACTGTTTTTCACCGTTAATGGTGATGGCATTGGAGGCGATAGTTTTACGTGCCTGACCACGAGAAGGTTGCAGTTCAGAATCGACCAGTGCCTGCATCAGGTCTGCGCCCTTTTCCATCTCAACCATCGGTACGCCGTCCTGCGCCAGCTGTTCGAAGTCCGCTTCACTCAGCGCACTCAAAGAACCGCTGAACAGGCATTCGGTAATACGTTTTGCCGCCTGTAAACcttct                                                                             >  2:192867/1‑251 (MQ=255)                                                                                                                                                                     agTATTCAGGATCGGACTGTTTTTCACCGTTAATGGTGATGGCATTGGAGGCGATAGTTTTACGTGCCTGACCACGAGAAGGTTGCAGTTCAGAATCGACCAGTGCCTGCATCAGGTCTGCGCCCTTTTCCATCTCAACCATCGGTACGCCGTCCTGCGCCAGCTGTTCGAAGTCCGCTTCACTCAGCGCACTCAAAGAACCGCTGAACAGGCATTCGGTAATACGTTTTGCCGCCTGTAAACCttcttc                                                                           >  2:497007/1‑250 (MQ=255)                                                                                                                                                                            aGGATCGGACTGTTTTTCACCGTTAATGGTGATGGCATTGGAGGCGATAGTTTTACGTGCCTGACCACGAGAAGGTTGCAGTTCAGAATCGACCAGTGCCTGCATCAGGTCTGCGCCCTTTTCCATCTCAACCATCGGTACGCCGTCCTGCGCCAGCTGTTCGAAGTCCGCTTCACTCAGCGCACTCAAAGAACCGCTGAACAGGCATTCGGTAATACGTTTTGCCGCCTGTAAACCTTCTTCACCGTGaa                                                                   <  1:280329/251‑1 (MQ=255)                                                                                                                                                                            aGGATCGGACTGTTTTTCACCGTTAAAGGTGATGGGATTGGCGGCGATAGTTTTACGTGCCTGACCAAGAGAAGGTTGCCGGTCAGAATCGACCACTGCCTGCATCAGGTCTGCGCCATTTTCCATCTCAACCATCGGTACGCCGTCCTGCGCCAGCTGTTCGAAGTCCGCTTCACTCAGCGCACTCAAAGAACCGCTGAACAGGCATTCGGTAATACGTTTTGCCGCCTGTAAACCTTCTTCACCGTGa                                                                    >  2:76123/1‑250 (MQ=255)                                                                                                                                                                             ggATCGGACTGTTTTTCACCGTTAATGGTGATGGCATTGGAGGCGATAGTTTTACGTGCCTGACCACGAGAAGGTTGCAGTTCAGAATCGACCAGTGCCTGCATCAGGTCTGCGCCCTTTTCCATCTCAACCATCGGTACGCCGTCCTGCGCCAGCTGTTCGAAGTCCGCTTCACTCAGCGCACTCAAAGAACCGCTGAACAGGCATTCGGTAATACGTTTTGCCGCCTGTAAACCTTCTTCACCGTGAAc                                                                  >  2:222302/1‑251 (MQ=255)                                                                                                                                                                              gATCGGACTGTTTTTCACCGTTAATGGTGATGGCATTGGAGGCGATAGTTTTACGTGCCTGACCACGAGAAGGTTGCAGTTCAGAATCGACCAGTGCCTGCATCAGGTCTGCGCCCTTTTCCATCTCAACCATCGGTACGCCGTCCTGCGCCAGCTGTTCGAAGTCCGCTTCACTCAGCGCACTCAAAGAACCGCTGCACAGGCATTCGGTAATACGTTTTGCCGCCTGTAAACCTTCTTCACCGTGAAcc                                                                 >  2:621690/1‑251 (MQ=255)                                                                                                                                                                                      tGTTTTTCACCGTTAATGGTGATGGCATTGGAGGCGATAGTTTTACGTGCCTGACCACGAGAAGGTTGCAGTTCAGAATCGACCAGTGCCTGCATCAGGTCTGCGCCCTTTTCCATCTCAACCATCGGTACGCCGTCCTGCGCCAGCTGTTCGAAGTCCGCTTCACTCAGCGCACTCAAAGAACCGCTGAACAGGCATTCGGTAATACGTTTTGCCGCCTGTAAACCTTCTTCACCGTGAACCAGACGAGt                                                         >  1:81607/1‑251 (MQ=255)                                                                                                                                                                                       gTTTTTCACCGTTAATGGTGATGGCATTGGAGGCTATAGTTTTACGTGCCTGACCACGAGAAGGTTGCAGTTCAGAATCGACCAGTGCCTGCATCAGGTCTGCGCCCTTTTCCATCTCAACCATCGGTACGCCGTCCTGCGCCAGCTGTTCGAAGTCCGCTTCACTCAGCGCACTCAAAGAACCGCTGAACAGGCATTCGGTAATACGTTTTGCCGCCTGTAAACCTTCTTCACCGTGAACCAGACGAGt                                                         <  2:382956/250‑1 (MQ=255)                                                                                                                                                                                         ttttCACCGTTAATGGTGATGGCATTGGAGGCGATAGTTTTACGTGCCTGACCACGAGAAGGTTGCAGTTCAGAATCGACCAGTGCCTGCATCAGGTCTGCGCCCTTTTCCATCTCAACCATCGGTACGCCGTCCTGCGCCAGCTGTTCGAAGTCCGCTTCACTCAGCGCACTCAAAGAACCGCTGAACAGGCATTCGTTAATACGTTTTGCCGCCTGTAAACCTTCTTCACCGTGAACCAGACGAGTCa                                                       <  2:261713/250‑1 (MQ=255)                                                                                                                                                                                         ttttCACCGTTAATGGTGATGGCATTGGAGGCGATAGTTTTACGTGCCTGACCACGAGAAGGTTGCAGTTCAGAATCGACCAGTGCCTGCATCAGGTCTGCGCCCTTTTCCATCTCAACCATCGGTACGCCGTCCTGCGCCAGCTGTTCGAAGTCCGCTTCACTCAGCGCACTCAAAGAACCGCTGAACAGGCATTCGGTAATACGTTTTGCCGCCTGTAAACCTTCTTCACCGTGAACCAGACGAGTCAc                                                      >  1:473040/1‑251 (MQ=255)                                                                                                                                                                                               ccGTTAATGGTGATGGCATTGGAGGCGATAGTTTTACGTGCCTGACCACGAGAAGGTTGCAGTTCAGAATCGACCAGTGCCTGCATCAGGTCTGCGCCCTTTTCCATCTCAACCATCGGTACGCCGTCCTGCGCCAGCTGTTCGAAGTCCGCTTCACTCAGCGCACTCAAAGAACCGCTGAACAGGCATTCGGTAATACGTTTTGCCGCCTGTAAACCTTCTTCACCGTGAACCAGACGAGTCACCTGCTc                                                >  1:463831/1‑251 (MQ=255)                                                                                                                                                                                                  ttAATGGTGATGGCATTGGAGGCGATAGTTTTACGTGCCTGACCACGAGAAGGTTGCAGTTCAGAATCGACCAGTGCCTGCATCAGGTCTGCGCCCTTTTCCATCTCAACCATCGGTACGCCGTCCTGCGCCAGCTGTTCGAAGTCCGCTTCACTCAGCGCACTCAAAGAACCGCTGAACAGGCATTCGGTAATACGTTTTGCCGCCTGTAAACCTTCTTCACCGTGAACCAGACGAGTCACCTGCTCCGc                                             >  1:9084/1‑251 (MQ=255)                                                                                                                                                                                                       ggTGATGGAATTGGAGGCGATAGTTTTACGTGCCTGCCCACGAGAAGGTTGCATTTCAAAATCGACCAGTGCCTGCATCTGGTCTGCGCCCTTTTCCATCTCAACCATCGGTACGCCGTCCTGCGCCAGCTGTTCGAAGTCCGCTTCCCTCAGCGCACTCAAAGAACCGCTGAACAGGCATTCGGTAATACGTTTTGCCGCCTGTAAACCTTCTTCACCGTGAACCAGACGAGTCACCTGCTCCGCCAGTa                                        <  2:138573/251‑1 (MQ=255)                                                                                                                                                                                                        gTGATGGCATTGGAGGCGATAGTTTTACGTGCCTGACCACGAGAAGGTTGCAGTTCAGAATCGACCAGTGCCTGCATCAGGTCTGCGCCCTTTTCCATCTCAACCATCGGTACGCCGTCCTGCGCCAGCTGTTCGAAGTCCGCTTCACTCAGCGCACTCAAAGAACCGCTGAACAGGCATTCGGTAATACGTTTTGCCGCCTGTAAACCTTCTTCACCGTGAACCAGACGAGTCACCTGCTCCGCCAGTAc                                       >  2:532058/1‑251 (MQ=255)                                                                                                                                                                                                              gCATTGGAGGCGATAGTTTTACGTGCCTGACCACGAGAAGGTTGCAGTTCAGAATCGACCAGTGCCTGCATCAGGTCTGCGCCCTTTTCCATCTCAACCATCGGTACGCCGTCCTGCGCCAGCTGTTCGAAGTCCGCTTCACTCAGCGCACTCAAAGAACCGCTGAACAGGCATTCGGTAATACGTTTTGCCGCCTGTAAACCTTCTTCACCGTGAACCAGACGAGTCACCTGCTCCGCCAGTACATACTg                                 >  1:80272/1‑251 (MQ=255)                                                                                                                                                                                                               cATTGGAGGCGATAGTTTTACGTGCCTGACCACGAGAAGGTGCCATTTCAGAATCGGCCAGTGCCTGCATCAGGTCTGCGCCCTTTTCCATGTCAACCATCGGTACGCCGTCCTGCGCCAGCTGTTCGAAGTCCGCTTCACTCAGCGCACTCAAAGAACCGCTGAACAGGCATTCGGTAATACGTTTTGCCGCCTGTAAACCTTCTTCACCGTGAACCAGACGAGTCACCTGCTCCGCCAGTACATACTgg                                <  2:269993/251‑1 (MQ=255)                                                                                                                                                                                                                  tGGAGGCGATAGTTTTACGTGCCTGACCACGAGAAGGTTGCAGTTCAGAATCGACCAGTGCCTGCATCAGGTCTGCGCCCTTTTCCATCTCAACCATCGGTACGCCGTCCTGCGCCAGCTGTTCGAAGTCCGCTTCACTCAGCGCACTCAAAGAACCGCTGAACAGGCATTCGGTAATACGTTTTGCCGCCTGTAAACCTTCTTCACCGTGAACCAGACGAGTCACCTGCTCCGCCAGTACATACTCGgcg                             >  2:173616/1‑251 (MQ=255)                                                                                                                                                                                                                               tttACGTGCCTGACCACGAGAAGGTTGCAGTTCAGAATCGACCAGTGCCTGCATAAGGTCTGCGCCCTTTTCCATCTCAACCATCGGTACGCGGCCCTGCGCCAGCTGTTCGAAGTCCGCTTCACTCAGCGCACTCAAAGAACCGCTGAACAGGCATTCGGTAATACGTTTTGCCGCCTGTAATCCTTCTTCACCGTGAACCAGACGAGTCACCTGCTCCGCCAGTACATACTGGGCGCGCGGTGCTTTa                 <  2:511497/250‑1 (MQ=255)                                                                                                                                                                                                                                   cGTGCCTGACCACGAGAAGGTTGCAGTTCAGAATCGACCAGTGCCTGCATCAGGTCTGCGCCCTTTTCCATCTCAACCATCGGTACGCCGTCCTGCGCCAGCTGTTCGAAGTCCGCTTCACTCAGCGCACTCAAAGAACCGCTGAACAGGCATTCGGTAATACGTTTTGCCGCCTGTAAACCTTCTTCACCGTGAACCAGACGAGTCACCTGCTCCGCCAGTACATACTGGGCGCGCGGTGCTTTACCGCt            <  1:265578/251‑1 (MQ=255)                                                                                                                                                                                                                                      gCCTGACCACGAGAAGGTTGCAGTTCAGAATCGACCAGTGCCTGCATCAGGTCTGCGCCCTTTTCCATCTCAACCATCGGTACGCCGTCCTGCGCCAGCTGTTCGAAGTCCGCTTCACTc                                                                                                                                            <  1:573212/120‑1 (MQ=255)                                                                                                                                                                                                                                      gCCTGACCACGAGAAGGTTGCAGTTCAGAATCGACCAGTGCCTGCATCAGGTCTGCGCCCTTTTCCATCTCAACCATCGGTACGCCGTCCTGCGCCAGCTGTTCGAAGTCCGCTTCACTc                                                                                                                                            >  2:573212/1‑120 (MQ=255)                                                                                                                                                                                                                                      gCCTGACCACGAGAAGGTTGCAGTTCAGAATCGACCAGTGCCTGCATCAGGTCTGCGCCCTTTTCCATCTCAACCATCGGTACGCCGTCCTGCGCCAGCTGTTCGAAGTCCGCTTCACTCAGCGCACTCAAAGAACCGCTGAACAGGCATTCGGTAATACGTTTTGCCGCCTGTAAACCTTCTTTACCGTGAACCAGACGAGTCACCTGCTCCGCCAGTACATACTGGGCGCGCGGTGCTTTACCGCTGtt         >  2:560932/1‑251 (MQ=255)                                                                                                                                                                                                                                      gCCTGACCACAAGAAGGTGGCAGTTCAAAACCGACCAGGGCCTGCATCAGGTCTTCGCCCTTTTCCATCCCAACCACCGGTACGCCTTCCTGCGCCAGCTGTTCTAATTCCGTTTCACTCAGCGCCCTCAAAGAACCGATGAACAGCCTTTCGGGAATACTTTTTGCCGCCTGTAAACCTTCTTCCCCGTGAGCCAGACGAGTCCCCTGCTCCGCCAGTACATACTGGGCGCGCGGTGCTTTACCGCTGt          <  2:410108/250‑1 (MQ=255)                                                                                                                                                                                                                                            ccACGAGAAGGTTGCAGTTAAGAATCGACCAGTGCATGCATCAGGCCTGCGACCTTTTCCATCTCAACCATCGGTACGCCGTCCTGCGACAGCTGTTCGAAGTCCGCTTCACTCAGCGCACTCACAGAACCGCTGAACAGGCATTCGGTAATACGTTTTGCCGCCTGTAAACCTTCTTCACCGTGCACCAGACGATTCACCTGCTCCGCCAGTACATACTGGGCGCGCGGTGCTTTACCGCTGTTTTTAtc   >  2:70826/1‑251 (MQ=255)                                                                                                                                                                                                                                             cACGAGAAGGTTGCAGTGCAGAATCGCCCAGTCCCTGCATCAGGTCTGCGCCCTTTCCCATCTCACCCATCGGTACGCCGTCCTGCGCCAGCTGTTCTAAGTCCGCTTCACTCAGCGCCCTCAAAGAACCGCTGAACAGGCATTCGTTAATACGTTTTGCCGCCTGTAAACCTTCTTCACCGTGAACCAGACGAGTCACTTGCTCCGCCAGTGCATACTGGGCGCGCTGTGCTTTTCCGCTGTTTTTAtct  <  2:388260/251‑1 (MQ=255)                                                                                                                                                                                                                                                 |                                                                                                                                                                                                                                                       GCGAGAATATTTTTCATCATGTTCCCTGTATAAAAACCAGGGGGAGTGATTTCTCACTCCCCCTTTCCACTTAATGCATTATTTCCAGCAAATCAGACAGTAATTCTTTTTACCGCGACGCAGTAAGGTAAAACGACCAAACAGACGATCTTCTTCTTTAAAGAAGTATTCAGGATCGGACTGTTTTTCACCGTTAATGGTGATGGCATTGGAGGCGATAGTTTTACGTGCCTGACCACGGGAAGGTTGCAGTTCAGAATCGACCAGTGCCTGCATCAGGTCTGCGCCCTTTTCCATCTCAACCATCGGTACGCCGTCCTGCGCCAGCTGTTCGAAGTCCGCTTCACTCAGCGCACTCAAAGAACCGCTGAACAGGCATTCGGTAATACGTTTTGCCGCCTGTAAACCTTCTTCACCGTGAACCAGACGAGTCACCTGCTCCGCCAGTACATACTGGGCGCGCGGTGCTTTACCGCTGTTTTTATCT  >  NC_000913/1715870‑1716356

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