Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	NC_000913	4,618,489	C→T	T87T (ACC→ACT)	deoA →	thymidine phosphorylase

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	4,618,489	0	C	T	98.4%	194.1 / NA	64	T87T (ACC→ACT)	deoA	thymidine phosphorylase
Reads supporting (aligned to +/- strand):  ref base C (0/0);  major base T (31/31);  minor base G (0/1);  total (32/32)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

AGAAATTATTCGTAAAAAACGTGATGGTCATGCGCTGAGCGATGAAGAAATTCGTTTCTTTATCAACGGTATTCGCGACAACACTATCTCCGAAGGGCAGATTGCCGCCCTCGCGATGACCATTTTCTTCCACGATATGACAATGCCTGAGCGTGTCTCGCTGACCATGGCGATGCGAGATTCAGGAACCGTTCTCGACTGGAAAAGCCTGCATCTGAATGGCCCGATTGTTGATAAACACTCCACCGGTGGCGTCGGCGATGTGACTTCGCTGATGTTGGGGCCGATGGTCGCAGCCTGCGGCGGCTATATTCCGATGATCTCTGGTCGCGGCCTCGGTCATACTGGCGGTACGCTCGACAAACTGGAATCCATCCCTGGCTTCGACATTTTCCCGGATGACAACCGTTTCCGCGAAATTATTAAAGACGTCGGCGTGGCGATTATCGGTCAGACCAGTTCACTGGCTCCGGCTGATAAACGTTTCTA  >  NC_000913/4618243‑4618731                                                                                                                                                                                                                                                       |                                                                                                                                                                                                                                                   aGAAATTATTCGTAAAAAACGTGATGGTCATGCACTGAGCGATGAAGAAATTCGTTTCTTTATCAACGGCATTCGCGACAACACTATCTCCGAAGGGCAGATTGCCGCCCTCGCGATGACCATTTTCTTCCACGATATGACAATGCCTGAGCGTGTCTCGCTGACCATGGCGATGCGAGATTCAGGAACCGTTCTCGACTGGAAAAGCCTGCATCTGAATGGCCCGATTGTTGATAAACATTCCACTGGTg                                                                                                                                                                                                                                                <  1:181985/251‑1 (MQ=255)       tattCGTAAAAAACGTGATGGTCATGCACTGCGCGATGAAGAAATTCGTTTCTTTATCAACGGCATTCGCGACAACAATATCTCCGAAGGGCAGATTGCCGCCCTCGCGATGACCATTTTCTTCCACGATATGCCAATGCCAGCGCGTGTCTTGCTGACCATGGCGATGCGCGATTCAGGAACCGTTCTCGACTGGAAACGCCTGCATCTGAATGGCCCGATTGTTGATAAAAATTCCACTGGTGGCgccg                                                                                                                                                                                                                                          >  2:362183/1‑248 (MQ=255)       tattCGTAAAAAACGTGATGGTCATGCACTGAGCGATGAAGAAATTCGTTTCTTTATCAACGGCATTCGCGACAACACTATCTCCGAAGGGCAGATTGCCGCCCTCGCGATGACCATTTTCTTCCACGATATGACAATGCCTGAGCGTGTCTCGCTGACCATGGCGATGCGAGATTCAGGAACCGTTCTCGACTGGAAAAGCCTGCATCTGAATGGCCCGATTGTTGATAAACATTCCACTGGTGGCGTc                                                                                                                                                                                                                                           <  1:266717/250‑1 (MQ=255)       tattCGTAAAAAACGTGATGGTCATGCACTGAGCGATGAAGAAATTCGTTTCTTTATCAACGGCATTCACGACAACACTATCTCCGAAGGGCAGATTGCCGCCCTCGCGATGACCATTTTCTTCCACGATATGACAATGCCTGAGCGTGTCTCGCTGACCATGGCGATGCGAGATTCAGGAACCGTTCTCGACTGGAAAAGCCTGCATCTGAATGGCCCGATTGTTGATAAACATTCCACTGGTGGCGTCg                                                                                                                                                                                                                                          >  1:506489/1‑251 (MQ=255)                  aaCGTGATGGTCATGCCCTGAGCGATGACGAAATTCGTTTCTTTATCAACGGCATTCGCGACAACACTATCTCCGAAGGGCAGATTGCCGCCCTCGCGATGACCATTTTCTTCCACGATATGACAATGCCTGAGCGTGTCTCGCTGACCATGGCGATGCGAGATTCAGGAACCGTTCTCGACTGGAAAAGCCTGCATCTGAATGGCCCGATTGTTGATAAACATTCCACTGGTGGCGTCGGCGATGTGACt                                                                                                                                                                                                                               >  1:423830/1‑251 (MQ=255)                       gATGGTCATGCACTGAGCGATGAAGAAATTCGTTTCTTTATCAACGGCATTCGCGACAACACTATCTCCGAAGGGCAGATTGCCGCCCTCGCGATGCCCATTTTCTTCCACGATATGACAATGCCTGAGCGTGTCTCGCTGACCATGGCGATGCGAGATTCAGGAACCGTTCTCGACTGGAAAAGCCTGCATCTGAATGGCCCGATTGTTGATAAACATTCCACTGGTGGCGTCGGCGATGTGACTTCGCt                                                                                                                                                                                                                          >  1:254832/1‑251 (MQ=255)                          ggTCATGCACTGAGCGATGAAGAAATTCGTTTCTTTATCAACGGCATTCGCGACAACACTATCTCCGAAGGGCAGATTGCCGCCCTCGCGATGACCATTTTCTTCCACGATATGACAATGCCTGAGCGTGTCTCGCTGACCATGGCGATGCGAGATTCAGGAACCGTTCTCGACTGGAAAAGCCTGCATCTGAATGGCCCGATTGTTGATAAACATTCCACTGGTGGCGTCGGCGATGTGACTTCGCTGAt                                                                                                                                                                                                                       >  1:399190/1‑251 (MQ=255)                             cATGCACTGAGCGATGAAGAAATTCGTTTCTTTATCAACGGCATTCGCGACAACACTATCTCCGAAGGGCAGATTGCCGCCCTCGCGATGACCATTTTCTTCCACGATATGACAATGCCTGAGCGTGTCTCGCTGACCATGGCGATGCGAGATTCAGGAACCGTTCTCGACTGGAAAAGCCTGCATCTGAATGGCCCGATTGTTGATAAACATTCCACTGGTGGCGTCGGCGATGTGACTTCGCTGATGtt                                                                                                                                                                                                                    <  2:178637/251‑1 (MQ=255)                              aTGCACTGAGCGATGAAGAAATTCGTTTCTTTATCAACGGCATTCGCGACAACACTATCTCCGAAGGGCAGATTGCCGCCCTCGCGATGACCATTTTCTTCCACGATATGACAATGCCTGAGCGTGTCTCGCTGACCATGGCGATGCGAGATTCAGGAACCGTTCTCGACTGGAAAAGCCTGCATCTGAATGGCCCGATTGTTGATAAACATTCCACTGGTGGCGTCGGCGATGTGACTTCGCTGATGTTg                                                                                                                                                                                                                   <  1:303401/251‑1 (MQ=255)                                gcacTGAACGATGAAGAAATTCGTTTCTTTATCAACGGCATTCGCGACAACACTATCTCCGAAGGGCAGATTGCCGCCCTCGCGATGACCATTTTCTTCCACGATATGACAATGCCTGAGCGTGTCTCGCTGACCATGGCGATGCGAGATTCAGGAACCGTTCTCGACTGGAAAAGCCTGCATCTGAATGGCCCGATTGGTGATAAACATTCCACTGGTGGCGTCGGCGATGTGACTACGCTGATGTTgg                                                                                                                                                                                                                  <  1:383812/247‑1 (MQ=255)                                  acTGAGCGATGAAGAAATTCGTTTCTTTATCAACGGCATTCGCGACAACACTATCTCCGAAGGGCAGATTGCCGCCCTCGCGATGACCATTTTCTTCCACGATATGACAATGCCTGAGCGTGTCTCGCTGACCATGGCGATGCGAGATTCAGGAACCGTTCTCGACTGGAAAAGCCTGCATCTGAATGGCCCGATTGTTGATAAACATTCCACTGGTGGCGTCGGCGATGTGACTTCGCTGATGTTGGGGc                                                                                                                                                                                                               >  2:119120/2‑251 (MQ=255)                                   cTGAGCGATGAAGAAATTCGTTTCTTTATCAACGGCATTCGCGACAACACTATCTCCGAAGGGCAGATTGCCGCCCTCGCGATGACCATTTTCTTCCACGATATGACAATGCCTGAGCGTGTCTCGCTGACCATGGCGATGCGAGATTCAGGAACCGTTCTCGACTGGAAAAGCCTGCATCTGAATGGCCCGTTTGTTGATAAACATTCCACTGGTGGCGTCGGCGATGTGCGTTCGCTGATGTTGGGGcc                                                                                                                                                                                                              >  1:346452/1‑251 (MQ=255)                                     gAGCGATGAAGAAATTCGTTTCTTTATCAACGGCATTCGCGACAACACTATCTCCGAAGGGCAGATTGCCGCCCTCGCGATGACCATTTTCTTCCACGATATGACAATGCCTGAGCGTGTCTCGCTGACCATGGCGATGCGAGATTCAGGAACCGTTCTCGACTGGAAAAGCCTGCATCTGAATGGCCCGATTGTTGATAAACATTCCACTGGTGGCGTCGGCGATGTGACTTCGCTGATGTTGGGGCCGa                                                                                                                                                                                                            <  1:371758/251‑1 (MQ=255)                                        cggtGAAGGAATTCGTTTCTTTATCAACGGCATTCGCGACAACACTATCGCCGAAGGGCAGATTGCCGCCCTCGCGATGGCCATTTTCTTCCACCATATGACAATGCCTGAGCGTGTCTCGCTGACCTTGGCTGTGCGAGATTCAGGAACCGTTCTCGTCTGGACACGCCTGCATCTGAATGGCCCGCTTGTTGATAAACATTCCCCAAGTGGCGACGGCGCTGTGACTTCGCTGGTGTTGGGGCCGATgg                                                                                                                                                                                                         >  2:144991/4‑251 (MQ=255)                                             aagaaATTCGTTTCTTTATCAACGGCATTCGCGACAACACTATCTCCGAAGGGCAGATTGCCGACCTCGCGATGACCATTTTCTTCCACGATATGACAATGCCTGAGCGTGTCTCGCTGACCATGGCGATGCGAGATTCAGGAACCGTTCTCGACTGGAAAAGCCTGCATCTGAATGGCCCGATTGTTGATAAACATTCCACTGGTGGCGTCGGCGATGTGACTTCGCTGATGTGGGGGCCGATGGTctcc                                                                                                                                                                                                    >  2:232072/1‑248 (MQ=255)                                                  aTTCGTTTCTTTATCAACGGCATTCGCGACAACACTATCTCCGAAGGGCAGATTGCCGCCCTCGCGATGACCATTTTCTTCCACGATATGACAATGCCTGAGCGTGTCTCGCTGACCATGGCGATGCGAGATTCATGAACCGTTCTCGACTGGAAAAGCCTGCATCTGAATGGCCCGATTGTTGCTAAACATTCCACTGGTGGCGTCGTCGCTGTGACTTCGCTGATGTTGGGGCCGATGGTCGCTGCCTg                                                                                                                                                                                               >  1:367932/1‑251 (MQ=255)                                                   ttCGTTTCTTTATCAACGGCATTCGCGACAACACTATCTCCGAAGGGCAGATTGCCGCCCTCGCGATGACCATTTTCTTCCACGATATGACAATGCCTGAGCGTGTCTCGCTGACCATGGCGATGCGAGATTCAGGAACCGTTCTCGACTGGAAAAGCCTGCATCTGAATGGCCCGATTGTTGATAAACATTCCACTGGTGGCGTCGGCGATGTGACTTCGCTGATGTTGGGGCCGATGGTCGCAGCCTgc                                                                                                                                                                                              >  1:323926/1‑251 (MQ=255)                                                   ttCGTTTCTTTATCAACGGCATTCGCGACAACACTATCTCCGAAGGGCAGATTGCCGCCATCGCGATGACCATTTTCTTCCACGATATGACAATGCCTGAGCGTGTCACGCTGACCATGGCGATGCGAGATTCAGGAACCGTTCTCGACTGGAAAAGCCTGCATCTGAATGGCCCGATTGTTGATAAACATTCCACTGGTGGCGTCGGCGATGTGACTTCGCTGATGTTGGGGCCGATGGACGCAGCCTgc                                                                                                                                                                                              >  1:146358/1‑251 (MQ=255)                                                     cGTTTCTTTATCAACGGCATTCGCGACAACACTATCTCCGAAGGGCAGATTGCCGCCCTCGCGATGACCATTTTCTTCCACGATATGACAATGCCTGAGCGTGTATCGCTGACCATGGCGATGCGAGATTCAGGAACCGTTCTCGACTGGAAAAGCCTGCATCTGAATGGCCCGATTGTTGATAAACATTCCACTGGTGGCGTCGGCGATGTGACTTCGCTGATGTTGGGGCCGATGGTCGCAGCcggcgg                                                                                                                                                                                            >  2:169944/1‑251 (MQ=255)                                                       tttCTTTATCAACGGCATTCGCGACAACACTATCTCCGAAGGGCAGATTGCCGCCCTCGCGATGACCATTTTCTTCCACGATATGACAATGCCTGAGCGTGTCTCGCTGACCATGGCGATGCGGGGTTCAGGAACCGTTCTCGACTGGAAAAACCTGCATCTGAATGGCCCGATTGTTGATACACATTACACTGGTGGCGTCGGCGATGTGACTTCGCTGATGTTGGGGCCGATGGTCGCAGCCTgcggcg                                                                                                                                                                                          >  1:89387/1‑251 (MQ=255)                                                             tATCAACGGCATTAGCGACAACACTATCTCCGAAGGGCAGATTGCCGCCCTCGCGCTGACCATTTTCTTCCACGATATGACAATGCCTGAGCGTGTCTCGCTGACCATGGCGCTGCGAGATTCAGGTACCGTTCTCGACTGGACAAGCTAGAATATGAAGGGCTCGATTGGTGATAAACATTCCACTGGTGACGTCGGCGTGGTGACTTCGCTGCTGTTGGGGCCGTTGGTCGCAGCGCGCggcggccatc                                                                                                                                                                                    >  2:188586/1‑247 (MQ=255)                                                                 aacggctttCGCGAAAACACTATCTCCGATGGGTATATTGCCGCCCTTGTGATGACTATTTTCTTCCAAGAGATGACAATGCCGGTGCGTGGCGCGCTGTCCTGGGCGATGCTAGATTCAGGAATCGTTCTCTACTGGAAAAGCCTGCATCTGAATGGCCCGATTGTTGATAAACATTCCACTGGTGGCGTCGGCGATGGGACTTCGCTGAGGTTGGGGCCGATGGTCGCAGCCTGCGGCGGCTATATTcc                                                                                                                                                                                <  2:270101/244‑1 (MQ=255)                                                                 aaCGGCATTCGCGACAACACTATCTCCGAAGGGCAGATTGCCGCCCTCGCGATGACCATTTTCTTCCACGATATGACAATGCCTGAGCGTGTCTCGCTGACCATGGCGATGCGAGATTCAGGAACCGTTCTCGACTGGAAAAGCCTGCATCTGAATGGCCCGATTGTTGATAAACATTCCACTGGTGGCGTCGGCGATGTGACTTCGCTGATGTTGGGGCCGGAGGTCGCAGCCTGCGGCGGCTATAtttc                                                                                                                                                                                >  2:336836/1‑249 (MQ=255)                                                                      caTTCGCGACAACACTATCTCCGAAGGGCAGATTGCCGCCCTCGCGATGACCATTTTCTTCCACGATATGACAATGCCTGAGCGTGTCTCGCTGACCATGGCGATGCGAGATTCAGGAACCGTTCTCGACTGGAAAAGCCTGCATCTGAATGGCCCGATTGTTGATAAACATTCCACTGGTGGCGTCGGCGATGTGACTTCGCTGATGTTGGGGCCGATGGTCGCAGCCTGCGGCGGCTATATTCCgatgt                                                                                                                                                                           >  1:70026/2‑250 (MQ=255)                                                                             gACAACACTATCTCCGAAGGGCAGATTGCCGCCCTCGCGATGACCATTTTCTTCCACGATATGACAATGCCTGAGCGTGTCTCGCTGACCATGGCGATGCGAGATTCAGGAACCGTTCTCGACTGGAAAAGCCTGCATCTGAATGGCCCGATTGTTGATAAACATTCCACTGGTGGCGTCGGCGATGTGACTTCGCTGATGTTGGGGCCGATGGTCGCAGCCTGCGGCGGCTATATTCCGATGATCTCTgg                                                                                                                                                                    <  1:161995/251‑1 (MQ=255)                                                                               caacaCTATCTCCGAAGGGCAGATTGCCGCCCTCGCGATGACCATTTTCTTCCACGATATGACAATGCCTGAGCGTGTCTCGCTGACCATGGCGATGCGAGATTCAGGAACCGTTCTCGACTGGAAAAGCCTGCATCTGAATGGCCCGATTGTTGATAAACATTCCACTGGTGGCGTCGGCGATGTGACTTCGCTGTTGTTGGGGCCGATGGTCGCAGCCTGCGGCGGCTATATTCCGTTGATCTCTgggc                                                                                                                                                                  >  2:273074/1‑249 (MQ=255)                                                                               aaacaCTATCTCCGAAGGGCAGATTGCCGCCCTCGCTATGACCATTTTCTTCCACGATATGACAAGGCCTGTGCGTGTCTCGCTGACCATGGCGATGCGAGATTCAGGAACCGTTCTCGACTGGAAAAGCCTGCATCTGAATGGCCCGATTGTTGATAAACATTCCACTGGTGGCGTCGGCGATGTGACTTCGCTGATGTTGGGGCCGATGGTCGCAGCCTGCGGCGGCTATATTCCGATGATCTCTGGTc                                                                                                                                                                  <  2:274853/250‑1 (MQ=255)                                                                                          ccGAAGGGCAGATTGCCGCCCTCGCGATGACCATTTTCTTCCACGATATGACAATGCCTGAGCGTGTCTCGCTGACCATGGCGATGCGAGATTCAGGAACCGTTCTCGACTGGAAAAGCCTGCATCTGAATGGCCCGATTGTTGATAAACATTCCACTGGTGGCGTCGGCGATGTGACTTCGCTGATGTTGGGGCCGATGGTCGCAGCCTGCGGCGGCTATATTCCGATGATCTCTGGTCGCGGCCTc                                                                                                                                                          <  1:119120/248‑1 (MQ=255)                                                                                          ccGAAGGGCAGATTGCCGCCCTCGCGATGACCATTTTCTTCCACGATATGACAATGCCTGAGCGTGTCTCGCTGACCATGGCGATGCGAGATTCAGGAACCGTTCTCGACTGGAAAAGCCTGCATCTGAATGGCCCGATTGTTGATAAACATTCCACTGGTGGCGTCGGCGATGTGACTTCGCTGATGTTGGGGCCGATGGTCGCAGCCTGCGGCGGCTATATTCCGATGATCTCTGGTCGCGGCCTCg                                                                                                                                                         <  1:336836/249‑1 (MQ=255)                                                                                            gAAGGGCAGATTGCCGCCCTCGCGATGACCATTTTCTTCCACGATATGACAATGCCTGAGCGTGTCTCGCTGACCATGGCGATGCGAGATTCAGGAACCGTTCTCGACTGGAAAAGCCTGCATCTGAATGGCCCGATTGTTGATAAACATTCCACTGGTGGCGTCGGCGATGTGACTTCGCTGATGTTGGGGCCGATGGTCGCAGCCTGCGGCGGCTATATTCCGCTGTTCTCTGGTCCCGGCCCCGGTcc                                                                                                                                                     >  2:142713/1‑250 (MQ=255)                                                                                             aaGGGCAGATAGCCGCCCTCGCGATGAACATTTTCTTCCACGATATGACAATGCCTGAGCGTGTCTCGCTGACCATGGCGATGCGAGATTCAGGAACCGTTCTCGACTGGAAAAGCCTGCATCTGAATGGCCCGATTGTTGATAAACATTCCACTGGTGGCGTCGGCGATGTGACTTCGCTGATGTTGGGGCCGATGGTCGCAGCCTGCGGCGGCTATATTCCGATGATCTCTGGTCGCGGCCTCGGTCAt                                                                                                                                                    >  1:429651/1‑251 (MQ=255)                                                                                                   agtttggcgctcgcgcgGTGACCAGTTTCTCCCACGATATGACAATGCCAGAGGGTTGCGCGCTGCACATGGCGAGGCGAGAATCAGTAAACGTTCTCGGTTGGAAAAGCCTTCACCTGAAGGGCCCGAGTGTTGATAAACATTCCACTGGTGGCGTCGGCGATGTGACTTCGCTGATGTTGGGGCCGATGGTAGCAGCCTGCGGCGGCTATATTCCGATGATCTCTGGTCGCGGCCTCGGTCAGACTGGc                                                                                                                                              <  2:262173/238‑1 (MQ=255)                                                                                                   aGATTGCCGCCCTCGCGATGACCATTTTCTTCCACGATATGACAATGCCTGAGCGTGTCTCGCTGACCATGGCGATGCGAGATTCAGGAACCGTTCTCGACTGGAAAAGCCTGCATCTGAATGGCCCGATTGTTGATAAACATTCCCCTGGTGGCGTCGGCGATGTGACTTCGCTGATGTTGGGGCCGATg                                                                                                                                                                                                          >  2:308954/1‑191 (MQ=255)                                                                                                   aGATTGCCGCCCTCGCGATGACCATTTTCTTCCACGATATGACAATGCCTGAGCGTGTCTCGCTGACCATGGCGATGCGAGATTCAGGAACCGTTCTCGACTGGAAAAGCCTGCATCTGAATGGCCCGATTGTTGATAAACATTCCACTGGTGGCGTCGGCGATGTGACTTCGCTGATGTTGGGGCCGATg                                                                                                                                                                                                          <  1:308954/191‑1 (MQ=255)                                                                                                       tGCCGCCCTAGCGATGACCTTTTTCGTCCCCGATATGACAATGGTTGATCGTGCCTCGCTGACCATGGCGATGCGAGATTCAGGATCCGTTCTCTACTGGAAAAGCCTGCATCTGAATGGCCCGATTGGTGATAAACATTCCGCTGGTGGCGCCGGCGATGTGTCTTCGCTGATGTTGTGGCCGATGGTCTCAGCCTGCGGCGGCGAGATTCCGATGAGCTCTGGTCGCGGCCCCGGTCTTACTGGCGGt                                                                                                                                           <  2:59473/250‑1 (MQ=255)                                                                                                             ccTCGCGATGACCATTTTCTTCCACGATATGACAATGCCTGAGCGTGTCTCGCTGACCATGGCGATGCGAGATTCAGGAACCGTTCTCGACTGGAAAAGCCTGCATCTGAATGGCCCGATTGTTGATAAACATTCCACTGGTGGCGTCGGCGATGTGACTTCGCTGATGTTGGGGCCGATGGTCGCAGCCTGCGGCGGCTATATTCCGATGTTCTCTGGTCGCGGCCTCGGTCATACTGGCGGTAccccgg                                                                                                                                    >  2:252259/1‑246 (MQ=255)                                                                                                                   gaggcccATTTTTTTCCACTATTTGACAATGCCTGAGCGTGTCTCGCTGACCATGGCGATGCGAGATTCAGGAACCGTTCTCGACTGGAAAAGCCTGCATCTGAATGGCCCGATTGTTGATAAACATTCCACTGGTGGCGTCGGCGATGTGACTTCGCTGATGTTGGGGCCGATGGTCGCTGCCTGCGGCGGCTATATTCCGATGATCTCTGGTCGCGGCCTCGGTCATACTGGCGGTACGCTCGACAAAc                                                                                                                              <  1:222496/246‑1 (MQ=255)                                                                                                                    aTGACCATTTTCTTCCACGATATGACAATGCCTGAGCGTGTCTCGCTGACCATGGCGATGCGAGATTCAGGAACCGTTCTCGACTGGAAAAGCCTGCATCTGAATGGCCCGATTGTTGATAAACATTCCACTGGTGGCGTCGGCGATGTGACTTCGCTGATGTTGGGGCCGATGGTCGCAGCCTGCGGCGGCTATATTCCGACGATCTCTGGTCGCGGCCTCGGTCATACTGGCGGGACGCCCGGCAAAcc                                                                                                                             >  2:389468/1‑250 (MQ=255)                                                                                                                    aTGACCATTTTCTTCCACGATATGACAATGCCTGAGCGTGTCTCGCTGACCATGGCGATGCGAGATTCAGGAACCGTTCTCGACTGGAAAAGCCTGCATCTGAATGGCCCGATTGTTGATAAACATTCCACTGGTGGCGTCGGCGATGTGACTTCGCTAATATTGGGGCCGGTGGTCGCAGCCTGCGGCGGCTATATTCCGAAGGTCTCTGGCCGAGGCCTCGGTCATAcatgcggttggctccacacaa                                                                                                                              >  1:275358/1‑230 (MQ=255)                                                                                                                     tGACCATTTTCTTCCACGATATGACAATGCCTGTGCGTGTCTCGCTGACCATGGCGATGCGAGATTCAGGAACCGTTCTCGACTGGAAAAGCCTGCATCTGAATGGCCCGATTGTTGATAAACATTCCACTGGTGGCGTCGGCGATGTGACTTCGCTGATGTTGGGGCCGATGGTCGCAGCCTGCGGCGGCTATATTCCGATGATCTCTGGGCGCGGCCTCGGTCATACTTGCgggtcgccccgcaaactc                                                                                                                            >  2:78211/1‑235 (MQ=255)                                                                                                                         cATTTTCTTCCACGATATGACAATGCCTGAGCGTGTCTCGCTGACCATGGCGATGCGAGATTCAGGAACCGTTCTCGACTGGAAAAGCCTGCATCTGAATGGCCCGATTGTTGATAAACATTCCACTGGTGGCGTCGGCGATGTGACTTCGCTGATGTTGGGGCCGATGGTCGCAGCCTGCGGCGGCTATATTCCGATGATCTCTGGTCGCGGCCTCGGTCATACTGGCGGTACGCTCGACAAACTGGAat                                                                                                                        <  1:248971/251‑1 (MQ=255)                                                                                                                                      gATATGACAATGCCTGAGCGTGTCTCGCTGACCATGGCGATGCGAGATTCAGGAACCGTTCTCGACTGGAAAAGCCTGCATCTGAATGGCCCGATTGTTGATAAACATTCCACTGGTGGCGTCGGCGATGTGACTTCGCTGATGTTGGGGCCGATGGTCGCAGCCTGCGGCGGCTATATTCCGATGATCTCTGGTCGCGGCCTCGGTCATACTGGCGGTTCGCTCGACAAACTGGAATCCATCCCTGGCtt                                                                                                           >  1:99981/1‑251 (MQ=255)                                                                                                                                           gACAATGCCTGAGCGTGTCTCGCTGACCATGGCGATGCGAGATTCAGGAACCGTTCTCGACTGGAAAAGCCTGCATCTGAATGGCCCGATTGTTGATAAACATTCCACTGGTGGCGTCGGCGATGTGACTTCGCTGATGTTGGGGCCGATGGTCGCAGCCTGCGGCGGCTATATTCCGATGATCTCTGGTCGCGGCCTCGGTCATACTGGCGGGACGCTCGGCAAACTGGAATCCATCCCTGGCTTCGAcc                                                                                                      >  1:407109/1‑250 (MQ=255)                                                                                                                                                tGCCTGAGCGTGTCTCGCTGACCATGGCGATGCGAGATTCAGGAACCGTTCTCGACTGGAAAAGCCTGCATCTGAATGGCCCGATTGTTGATAAACATTCCACTGGTGGCGTCGGCGATGTGACTTCGCTGATGTTGGGGCCGATGGTCGCAGCCTGCGGCGGCTATATTCCGATGATCTCTGGTCGCGGCCTCGGTCATACTGGCGGTACGCTCGACCAACTGGAATCCCTCCCTGGCTTCGACAttttt                                                                                                 >  1:218149/1‑250 (MQ=255)                                                                                                                                                   cagagcgggactcGCAGACCATGTCGAGGCGAATTTAAGGAACCGGTCGCGACTGGAAAAGCCTGCATCTGAATGTCCCGATTGTTGATAAACATTCCACTGGTGGCGTCGGCGATGTGACTTCGCTGAGGTTGGGGCCGATGGTCGCAGCCTGCGGCGGCTATAGTCCGATGATCTCTGGTCGCGGCATCGGGCATACTGGCGGTACGCTCGCCAAACTGGAATCCATCCCTGGCTTCTACATTTTCCCg                                                                                              <  2:146358/241‑1 (MQ=255)                                                                                                                                                         gtgtCTCGCTGCCCATGGCGATGCGAGATTCAGGAACCGTTCTCGACTGGAAAAGCCTGCATCTGAATGGCCCGATTGTTGATAAACATTCCACTGGTGGCGTCGGCGATGTGACTTCGCTGATGTTGGGGCCGATGGTCGCAGCCTGCGGCGGCTATATTCCGATGATCTCTGGTCGCGGCCTCGGTCATACTGGCGGTACGCTCGACAAACTGGAATCCATCCCTGTCTTCGACATTTTCCCGGATGAc                                                                                        <  1:389468/251‑1 (MQ=255)                                                                                                                                                         gtgtCTCGCTGACCATGGCGATGCGAGATTCAGGAACCGTTCTCGACTGGAAAAGCCTGCATCTGAATGGCCCGATTGTTGATAAACATTCCACTGGTGGCGTCGGCGATGTGACTTCGCTGATGTTGGGGCCGATGGTCGCAGCCTgcgg                                                                                                                                                                                            <  1:436891/151‑1 (MQ=255)                                                                                                                                                         gtgtCTCGCTGACCATGGCGATGCGAGATTCAGGAACCGTTCTCGACTGGAAAAGCCTGCATCTGAATGGCCCGATTGTTGATAAACATTCCACTGGTGGCGTCGGCGATGTGACTTCGCTGATGTTGGGGCCGATGGTCGCAGCCTgcgg                                                                                                                                                                                            >  2:436891/1‑151 (MQ=255)                                                                                                                                                          tgtCTCGCTGACCATTGTGATGCGAGATTCAGGAACCGTTCCCGACTGGAAAAGCCTGAAGCTGAATGGCCCGATTGGTGATAAACATTCCACTGGTGGCGTCGGCGATGTGACTTCGCTGTTGTTGGGGCCGATGGTCGCCGCCTGCGGCGGCTATATTCCGATGATCTCTGGTCGCGGCCTCGGTCATACTGGCGGTACGCTCGACAAACTGGAATCCATCCCTGGCTTCGACATTTTCCCGGATGAc                                                                                        <  2:399190/250‑1 (MQ=255)                                                                                                                                                             ctcGCTGACCATGGCGATGCGGGATTCAGGAACCGTTCTCAACTGGAAAAGCCTGCATCTGAATGGCCCGATTGTTGATAAACATTCCACTGGTGGCGTCGGCGATGTGACTTCGCTGATGTTGGGGCCGATGGTCGAAGCCTGCGGCGGCTATATTCCGATGATCTCTGGTCGCGGCCTCGGTACTACTGGCGGTAATCTCGTCCAACTAGAATTCAACCCCGGCGTCGGAATTTTCCCGGagtcaaacc                                                                                    >  1:441375/1‑243 (MQ=255)                                                                                                                                                             atcgccacccatggcggagccagcgtCAGGATCCGTTCCCGATTGTAAAAGCCTGCATTGGAATGGTCCGATTGTTGATAAACATCCCACTGGTGGCGGCGGGCATGTGACTTCGCTGATGTGGGGTCCGATGGTCGCAGCCTGCGGCGGCGATATGCCGATGATCTCTGGGCGCGTCCGCGGTCATACTGGCGGTACGCTCGACAAACTGGAATCCATCCCTGGCTTCGACATTGTACCGTATGACAAc                                                                                     <  2:383391/225‑1 (MQ=255)                                                                                                                                                                   gACCATGGCGATGCGAGATTCAGGAACCGTTCTCGACTGGAAAAGCCTGCATCTGAATGGCCCGATTGTTTATAAACATGCCACTGGTGGCGGCGGCGATGTGACTTCGCTGATGTTGGGGCCGATGGTCGCAGCCTGCGGCGGCTATATTCCGATGATCTCTGGTCGCGGCCTCGGTCATACTGGCGGTACGCTCGACAAACTGGAATCCATCCCTGGCTGCGGCATTTTCCCGGATGACAACCGTTTcc                                                                              <  2:64222/251‑1 (MQ=255)                                                                                                                                                                        tGGCGATGCGAGATTCAGGAACCGTTCTCGACTGGAAAAGCCTGCATCTGAATGGCCCGATTGTTGATAAACATTCCACTGGTGGCGTCGGCGATGTGACTTCGCTGATGTTGGGGCCGATGGTCGCAGCCTGCGGCGGCTATATTCCGATGATCTCTGGTCGCGGCCTCGGTCATACTGGCGGTACGCCCGACACACTGGAATCCCTCCCTGGGTTCGGCATTTTCCCGGGTGACAACCGGTTCcgcgg                                                                          >  1:185852/1‑249 (MQ=255)                                                                                                                                                                         ggcagagCGAGATTCAGGAACCCTTTTCGACTGGAAAATCCTTCATCGGACTGGCCCGATTGTTGATAAACATGCCACTGGTGGCGTCGGCGATGTTACTCCGCTTATGTTGGGGCCGATGGTCGCTGCCTGCGGCGGCTATATTCCGATGATCTCTGGTCGCGGCCTCGGTCATACTGGCGGTACGCGCGACAAACTTGAAGCCAACACTGGCTTCGTCATTTTCCCGGATGATAACCGTTGCCGCGaaa                                                                        <  2:55006/245‑1 (MQ=255)                                                                                                                                                                                  agaTTCAGGAACCGTTTCCGACTGGAAAAGCCTGCATCTGAATGGCCCGATTGTTGATAAACATTCCTCTGGTGGCGTCGGCGATGTGACTTCGCTGATGTTGGGGCCGATGGTCGCAGCCTGCGGCGGCTATATTCCGATGATCTCTGGTCGCGGCCTCGGTCATACTGGCGGTACGCTCGACAAACTGGAATCCATCCCTGGCTTCGACATTTTCCCGGATGACAACCGTTTCCGCGAAATTATTAAAg                                                               <  2:423830/251‑1 (MQ=255)                                                                                                                                                                                   agatccggAACCCTTCTCCACTGGAAAAGCCTGCCTGTGAATGGCCCGGTTTTTTAATAACCTTCCCCTGGTGGGGTCCGGGGTTTGAATTCCCTGATTGTGGGGCCGATGGTCGCAGCCTGCGGCGGCTATATTCCGATGATCTCTGGTCGCGGCCTCGGTCATACTGGCGGTACGCTCGACAAACTGGAATCCATCCCTGGCTTCGACATTTTCCCGGATGACAACCGTTTCCGCGAAATTATTAAAGa                                                              <  1:324351/245‑1 (MQ=255)                                                                                                                                                                                           aaCCGTTCTCGACTGGAAAAGCCTGCATCTGAATGGCCCGATTGTTGATAAACATTCCACTGGTGGCGTCGGCGATGTGACTTCGCTGATGTTGGGGCCGATGGTCGCAGCCTGCGGCGGCTATATTCCGATGATCTCTGGTCGCGGCCTCGGTCATACTGGCGGTACGCTCGACAAAATGGAATCCATCCCTGGCTTCGACATTTTCCCGGATGACAACCGTTTCCGCGAAATTATTAAAGACGTCGGCg                                                      >  1:216438/1‑251 (MQ=255)                                                                                                                                                                                                          gAAAAGCCTGCATTGGAATGGCCCGCTTGTTGATAAACATTCCTCGGGGGGCGTCGGCGATGGGACTGCGCTGATGTTGGTGCAGATGGGTGCAGCCTGCGGCGGCTATTTTCCGAGGATCTCTGGCCGCGGCCTCGGGCATACTGGCGGTACGCTCGACAAACTGGAATCCATCCCTGGCTTCGACATTTTCCCGGATGACAACCGTTGCCGCGAAATTATTAAAGTCGTCGGCGTGGCGATTATCGGTc                                       <  2:275358/251‑1 (MQ=255)                                                                                                                                                                                                              aGCCTGCATCTGAATGGCCCGATTGTTGATAAACATTCCACTGGTGGCGTCGGCGATGTGACTTCGCTGATGTTGGGGCCGATGGTCGCAGCCTGCGGCGGCTATATTCCGATGATCTCTGGTCGCGGCCTCGGTCATACTGGCGGTACGCTCGACAAACTGGAATCCATCCCTGGCTTCGACATTTTCCCAGATGACAACCGTTTCCGCGAAATTATTAAAGACGTCGGCGTGGCGATTATCGGTCAGAc                                   <  1:269478/251‑1 (MQ=255)                                                                                                                                                                                                                         gAATGGCCCGATTGTTTATAAACATTCCACTGGTGGCGTCGGCGATGTGACTTCGCTGATGTTGGGGCCGATGGTCGCAGCCTGCGGCGGCTATATTCCGATGATCTCTGGTCGCGGCCTCGGTCATACTGGCGGTACGCTCGACAAACTGGAATCCATCCCTGGCTTCGACATTTTCCCGGATGACAACCGTTTCCGCGAAATTATTAAAGACGTCGGCGTGGCGATTATCGGTCAGACCAGCTCACt                          <  2:323926/249‑1 (MQ=255)                                                                                                                                                                                                                         gAATGGCCCGATTGTTGATAAACATTCCACTGGTGGCGTCGGCGATGTGACTTCGCTGATGTTGGGGCCGATGGTCGCAGCCTGCGGCGGCTATATTCCGATGATCTCTGGTCGCGGCCTCGGTCATACTGGCGGTACGCTCGACAAACTGGAATCCATCCCTGGCTTCGACATTTTCCCGGATGACAACCGTTTCCGCGAAATTATTAAAGACGTCGGCGTGGCGATTATCGGTCAGACCAGCTCACTgg                        <  1:169944/251‑1 (MQ=255)                                                                                                                                                                                                                             ggCCCGATTGTCGATAAACATTCCACTGGTGGCGTCGGCGATGTGACTTCGCTGATGTTGGGGCCGATGGTCGCAGCCTGCGGCGGCTATATTCCGATGATCTCTGGTCGCGGCCTCGGTCATACTGGCGGTACGCTCGACAAACTGGAATCCATCCCTGGCTTCGACATTTTCCCGGATGACATCCGTTTCCGCGAAATTATTAAAGACGTCGGCGTGGCGATTATCGGTCAGACCAGCTCACTGGCTc                     <  1:397288/250‑1 (MQ=255)                                                                                                                                                                                                                                               cattCCACTGGTTGCTTCGGCGATGTGACTTCGCTGATGTTGGGGCCGATGGTCGCAGCCTGCGGCGGCTATATTCCGATGATCTCTGGTCGCGGCCTCGGTCATACTGGCGGTACGCTCGACAAACTGGAATCCATCCCTGGCTTCGACATTTTCCCGGATGACAACCGTTTCCGCGAAATTATTAAAGACGTCGGCGTGGCGATTATCGGTCAGACCAGCTCACTGGCTCCGGCGGATAAACGTTTCTa  <  1:9/248‑1 (MQ=255)                                                                                                                                                                                                                                               cattCCACTGGTGGCGTCGGCGATGTGACTTCGCTGATGTTGGGGCCGATGGTCGCAGCCTGCGGCGGCTATATTCCGATGATCTCTGGTCGCGGCCTCGGTCATACTGGCGGTACGCTCGACAAACTGGAATCCATCCCTGGCTTCGACATTTTCCCGGATGACAACCGTTTCCGCGAAATTATTAAAGACGTCGGCGTGGCGATTATCGGTCAGACCAGCTCACTGGCTCCGGCGGATAAACGTTTCTa  <  1:354455/248‑1 (MQ=255)                                                                                                                                                                                                                                                       |                                                                                                                                                                                                                                                   AGAAATTATTCGTAAAAAACGTGATGGTCATGCGCTGAGCGATGAAGAAATTCGTTTCTTTATCAACGGTATTCGCGACAACACTATCTCCGAAGGGCAGATTGCCGCCCTCGCGATGACCATTTTCTTCCACGATATGACAATGCCTGAGCGTGTCTCGCTGACCATGGCGATGCGAGATTCAGGAACCGTTCTCGACTGGAAAAGCCTGCATCTGAATGGCCCGATTGTTGATAAACACTCCACCGGTGGCGTCGGCGATGTGACTTCGCTGATGTTGGGGCCGATGGTCGCAGCCTGCGGCGGCTATATTCCGATGATCTCTGGTCGCGGCCTCGGTCATACTGGCGGTACGCTCGACAAACTGGAATCCATCCCTGGCTTCGACATTTTCCCGGATGACAACCGTTTCCGCGAAATTATTAAAGACGTCGGCGTGGCGATTATCGGTCAGACCAGTTCACTGGCTCCGGCTGATAAACGTTTCTA  >  NC_000913/4618243‑4618731

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