Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	NC_000913	1,789,888	G→A	P453S (CCT→TCT)	selO ←	protein adenylyltransferase SelO

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	1,789,888	0	G	A	97.5%	124.7 / ‑7.1	41	P453S (CCT→TCT)	selO	protein adenylyltransferase SelO
Reads supporting (aligned to +/- strand):  ref base G (1/0);  new base A (24/15);  total (26/15)
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

CCAATCACAAGACCGCTTCCCCGCAGCCTATTCGGCGAATTTCCAGCCAGACACTGTTAGGTCCGGATGGCAAACTGATTATCGATCATGACGGGCAAGAATATCTGCTCCGTAAAACTCAGGCTGGCAAGCTGCTGTTGACCAAGTAGCCTTTAACTCGAGCAGCTGACTTCCAGCCGTTTACCCCAGTCAGGTGGACGACTGACATAGTCATCATCTCTGTCGCTGAAAGGATTTCGCAACGCCTCATGCAGGCGGTGCAATTCCGTCATATC‑ACCCTTTTCTGCCGCTTCAATCGCCCGTTGCGCCAACCAGTTGCGCAATACCAGAGCGGGGTTAACGCTTTGCATCAGTTGCTGACGCTCACTATCGCTAACCTCGTCTTGTTGCAAACGCCCCCGATAACGGGCAAACCAGTCATCAAATGCCGCAC‑GATCAATAAACTCATCACGTAGCGGTGACGCCGCGCTGT  >  NC_000913/1789656‑1790127                                                                                                                                                                                                                                         |                                                                                                                                                                                                                                                  ccAATCACAAGACCGCTTCCCCGCAGCCTATTCGGCGAATTTCCAGCCAGACACTGTTAGGTCCGGATGGCAAACTGATTATCGATCATGACGGGCAAGAATATCTGCTCCGTAAAACTCAGGCTGGCAAGCTGCTGTTGACCAAGTAGCCTTTAACTCGAGCAGCTGACTTCCAGCCGTTTACCCCAGTCAGGTGGACGACTGACATAGTCATCATCTCTGTCGCTGAAAGAATTTCGCAACGCCTCATg                                                                                                                                                                                                                                 <  1:533430/251‑1 (MQ=255)    atcacaagactgcatcaaagCAGCCTATTCCGCGAATTTCCAGCCAGACACTGTTAGGACAGGATGGCAAACTGATTATCGATCATGACGGGCAAGAATATCCGCTCAGTAAAACTCAGGCAGGCAAGCTGCAGTTGACCAAGTAGTCTTTAACTCGAGCAGCTGATTTCCAGCCGTTTACCCCAGACAGTAGGAAGACTGACATAGTCATCATCTCTGTCGCTGCAAGAATTTCGCAACGCCTCATGc                                                                                                                                                                                                                                >  2:229898/20‑249 (MQ=255)              cGCTTCCCCGCAGCCTATTCGGCGAATTTCCAGCCAGACACTGTTAGGTCCGGATGGCAAACTGATTATCGATCATGACGGGCAAGAATATCTGCTCCGTAAAACTCAGGCTGGCAAGCTGCTGTTGACCAAGTAGCCTTTAACTCGAGCAGCTGACTTCCAGCCGTTTACCCCAGTCAGGTGGACGACTGACATAGTCATCATCTCTGTCGCTGAAAGAATTTCGCAACGCCTCATGCAGGCGGTGCAAt                                                                                                                                                                                                                    >  1:608502/1‑251 (MQ=255)                          gCCTATTCGGCGAATTTCCAGCCAGACACTGTTAGGTCCGGATGGCAAACTGATTATCGATCATGACGGGCAAGAATATCTGCTCCGTAAAACTCAGGCTGGCAAGCTGCTGTTGACCAAGTAGCCTTTAACTCGAGCAGCTGACTTCCAGCCGTTTACCCCAGTCAGGTGGACGACTGACATAGTCATCATCTCTGTCGCTGAAAGAATTTCGCAACGCCTCATGCAGGCGGTGCAATTCCGTCATATC‑a                                                                                                                                                                                                       <  2:1795/251‑1 (MQ=255)                              aTTCGGCGAATTTCCAGCCAGACACTGTTAGGTCCGGATGGCAAACTGATTATCGATCATGACGGGCAAGAATATCTGCTCCGTAAAACTCAGGCTGGCAAGCTGCTGTTGACCAAGTAGCCTTTAACTCGAGCAGCTGACTTCCAGCCGTTTACCCCAGTCAGGTGGACGACTGACATAGTCATCATCTCTGTCGCTGAAAGAATTTCGCAACGCCTCATGCAGGCGGTGCAATTCCGTCATATC‑ACCCt                                                                                                                                                                                                   >  1:282781/1‑251 (MQ=255)                                                cagACACTGTTAGGTCCGGATGGCAAACTGATTATCGATCATGACGGGCAAGAATATCTGCTCCGTAAAACTCAGGCTGGCAAGCTGCTGTTGACCAAGTAGCCTTTAACTCGAGCAGCTGACTTCCAGCCGTTTACCCCAGTCAGGTGGACGACTGACATAGTCATCATCTCTGTCGCTGAAAGAATTTCGCAACGCCTCATGCAGGCGGTGCAAATCCGTCATATCCACCCTTT‑CTGCCGCTTCAATCg                                                                                                                                                                                 >  1:332052/1‑251 (MQ=255)                                                 agACACTGTTCGGTCCGGATGGCAAACTGATTATCGATCATGACGGGCAAGCATATCTGCTCCGTAAAACTCAGGCTGGCAAGCTGCTGTTGACCAAGTAGCCTTTAACTCGCGCAGCTGACTTCCAGCCGGTTACCCCAGTCAGGGGGACGACTGACATAGTCATCATCTCTGTCGCTGAAAGACTTTCGCAACGCCTCATGCAGGCGGCGCAATTCCGTCCAATA‑CCCCCTTTCTGCCGCCTCAATCGc                                                                                                                                                                                >  2:241013/1‑251 (MQ=255)                                                  gACACTGTTAGGTCCGGATGGCAAACTGATTATCGATCATGACGGGCAAGAATATCTGCTCCGTAAAACTCAGGCTGGCAAGCTGCTGTTGACCAAGTAGCCTTTAACTCGAGCAGCTGACTTCCAGCCGTTTACCCCAGTCAGGTGGACGACTGACATAGTCATCATCTCTGTCGCTGCAAGAATTTCGCAACGCCTCATGCAGGCGGTGCAATTCCGTCATATC‑ACCCTTTTTTGCCGCTTCAATCGcc                                                                                                                                                                               >  1:427390/1‑251 (MQ=255)                                                   acacTGTTAGGTCCGGATGGAAAACTGATTATCGATCATGACGGGCAAGAATATCTGCTCCGTAAAACTCAGGATGGCAAGCTGCTGTTGAACAAGTAGCCTTTAACTCGAGCAGCTGACTTCCAGCCGTTTACCCCAGTCAGATGGACGACTGACATAGTCATCATCTCTGTCTCTGAACGGATTTCGTAACGCCTCATGCAGGCGGTACAATTTCGTCCTATC‑ACCCTTTTCTGGCGCTTTAATCGccc                                                                                                                                                                              >  2:351742/1‑251 (MQ=255)                                                                     ggCAAACTGATTATCGATCATGACGGGCAAGAATATCTGCTCCGTAAAACTCAGGCTGGCAAGCTGCTGTTGACCAAGTAGCCTTTAACTCGAGCAGCTGACTTCCAGCCGTTTACCCCAGTCAGGTGGACGACTGACATAGTCATCATCTCTGTCGCTGAAAGAATTTCGCAACGCCTCATGCAGGCGGTGCAATTCCGTCATATC‑ACCCTTTTCTGCCGCTTCAATCGCCCGTTGCGCCAACCAGTTgc                                                                                                                                                            <  1:251832/251‑1 (MQ=255)                                                                             gATTATCGATCATGACGGGCAAGAATATCTGCTCCGGAAAACTCAGGCTGGCAAGCTGCTGTTGACCAAGTAGCCTTTAACTCGAGCATCTGACTGCCAGCCGTTTACCCCAGTCAGGTGGACGACTGACATAGTCATCATCTCTGTCGCTGAAAGAATTTCGCAACGTCTCATGCAGGCTGTGCAATTGCTTCATATC‑TCCCCGTGCTCCCGCTTCAATCGCCCTTTGTGCCAAGGAGTTGCCCAAtcc                                                                                                                                                     <  2:297512/250‑3 (MQ=255)                                                                                 aTCGATCATGACGGGCAAGAATATCTGCTCCGTAAAACTCAGGCTGGCAAGCTGCTGTTGACCAAGTAGCCTTTAACTCGAGCAGCTGACTTCCAGCCGTTTCCCCCAGTCAGGTGGACGACTGACATAGTCATCATCTCTGTCGCTGAAAGAATTTCGCAACGCCTCATGCAGGCGGTGCAATTCCGTCATATC‑ACCCTTTTCTGCCGCTTCAATCGCCCGTTGCGCCAACCAGTTGCGCAATCCCaga                                                                                                                                                 <  2:332052/250‑1 (MQ=255)                                                                                 aTCGATCATGACGGGCAAGAATATCTGCTCCGTAAAACTCAGGCTGGCAAGCTGCTGTTGACCAAGTAGCCTTTAACTCGAGCAGCTGACTTCCAGCCGTTTACCCCAGTCAGGTGGACGACTGACATAGTCATCATCTCTGTCGCTGAAAGAATTTCGCAACGCCTCATGCAGGCGGTGCAATTCCGTCCTATC‑AACCTTTTCTTCCGCTTCAATCGCCCGTTGCGCCCACCAGTTGCGCGATACCaagg                                                                                                                                                >  1:355741/1‑248 (MQ=255)                                                                                    gATCATGACGGGCAAGAATATCTGCTCCGTAAAACTCAGGCTGGCAAGCTGCTGTTGACCAAGTAGCCTTTAACACGAGCAGCTGACTTCCAGCCGTTTACCCCAGTCAGGTGGACGACTGACATAGTCATCATCTCTGTCGCTGAAAGAATTTCGCAACGCCTCATGCAGGCGGTGCAATTCCGTCATATC‑ACCCTTTTCTGCCGCTTCAATCGCCCGTTGCGCCAACCAGTTGCGCAATACCAGAGCgg                                                                                                                                             >  1:574483/1‑251 (MQ=255)                                                                                     aTCATGACGGGCAAGAATATCTGCTCCGTAAAACTCAGGCTGGCAAGCTGCTGTTGACCAAGTAGCCTTTAACTCGAGCAGCTGACTTCCAGCCGTTTACCCCAGTCAGGTGGACGACTGACATAGTCATCATCTCTGTCGCTGAAAGAATTTCGCAACGCCTCATGCAGGCGGTGCAATTCCGTCATATC‑ACCCTTTTCTGCCGCTTCAATCGCCCGTTGCTCCAACCAGTTGCGCAATACCAGAGCggg                                                                                                                                            <  1:538498/251‑1 (MQ=255)                                                                                       cATGAAGGGCAAGCATATCTGCTCCGTAAAACTCAGGCTGGCACGCTGCTGTTGACCAAGTAGCCTTTAACTCGAGCAGCTGACTTCCAGCCGTTTACCCCAGTCAGGTGGCCGACTGACATAGTCATCATCTCTGTCGCTGAAAGAATTTCGCAACGCCTCATGCAGGCGGTGCAATTCCGTCATATC‑ACCCTTTTCTGCCGCTTCAATCGCCCGTTGCGCCAACCAGTTTCGCAATCACAGAGCGGGGt                                                                                                                                          >  2:309931/1‑251 (MQ=255)                                                                                        aTGACGGGCAAGAATATCTGCTCCGTAAAACTCAGGCTGGCAAGCTGCTGTTGACCAAGTAGCCTTTAACTCGAGCAGCTGACTTCCAGCCGTTTACCCCAGTCAGGTGGACGACTGACATAGTCATCATCTCTGTCGCTGAAAGAATTTCGCAACGCCTCATGCAGGCGGTGCAATTTCGTCATATT‑ACCCTTTTCTGCCGCTTCAATCGCCCGTTTCGCCCACCAGTTCCGCAATACCCGAGCgggggt                                                                                                                                         >  2:392421/1‑249 (MQ=255)                                                                                                   gAATATCTGCTCCGTAAAACTCAGGCTGGCAAGCTGCTGTTGACCAAGTAGCCTTTAACTCGAGCAGCTGACTTCCAGCCGTTTACCCCAGTCAGGTGGACGACTGACATAGTCATCATCTCTGTCGCTGAAAGAATTTCGCAACGCCTCATGCAGGCGGTGCAATTCCGTCATATC‑ACCCTTTTCTGCCGCTTCAATCGCCCGTTGCGCCAACCAGTTGCGCAATACCAGAGCGGGGTTAACGCTTTGc                                                                                                                               >  2:176706/1‑250 (MQ=255)                                                                                                          tGCTCCGTAAACCTCAGGCTGGCAAGCTGCTGTTGACCAAGTAGCCTTTAACTCGAGCAGCTGACTTCCAGCCGTTTACCCCAGTCAGGTGGACGACTGACATAGTCATCATCTCTGTCGCTGAAAGAATTTCGCAACGCCTCATGCAGGCGGTGCAATTCCGTCATATC‑ACCCTTTTCTGCCGCTTCAATCGCCCGTTGCGCCAACCAGTTGCGCAATACCAGAGCGGGGTTAACGCTTTGCATCAGtt                                                                                                                        <  2:494091/250‑1 (MQ=255)                                                                                                               cGTAAAACTCAGGCTGGCAAGCTGCTGTTGACCAAGTAGCCTTTAACTCGAGCAGCTGACTTCCAGCCGTTTACCCCAGTCAGGTGGACGACTGACATAGTCATCATCTCTGTCGCTGAAAGAATTTCGCAACGCCTCATGCAGGCGGTGCAATTCCGTCATATC‑ACCCTTTTCTGCCGCTTCAATCGCCCGTTGCGCCAACCAGTTGCGCAATACCAGAGCGGGGTTAACGCTTTGCATCAGTTGCTGAc                                                                                                                  >  1:399342/1‑251 (MQ=255)                                                                                                                                     tgctgTTGACCAAGTAGCCTTTAACTCGAGCAGCTGACTTCCAGCCGTTTACCCCAGTCAGGTGGACGACTGACATAGTCATCATCTCTGTCGCTGAAAGAATTTCGCAACGCCTCATGCAGGCGGTGCAATTCCGTCATATC‑ACCCTTTTCTGCCGCTTCAATCGCCCGTTGCGCCAACCAGTTTCGCAATACCAGAGCGGGGTTAACGCTTTGCATCAGTTGCTGACGCTCACTATCGCTAACCTCGTc                                                                                            >  2:244790/1‑251 (MQ=255)                                                                                                                                      gctgTTGACCAAGTAGCCTTTAACTCGAGCAGCTGACTTCCAGCCGTTTACCCCAGTCAGGTGGACGACTGACATAGTCATCATCTCTGTCGCTGAAAGAATTTCGCAACGCCTCATGCAGGCGGTGCAATTCCGTCATATC‑ACCCTTTTCTGCCGCTTCAATCGCCCGTTGCGCCAACCAGTTGCGCAAAACCCGAGCGGGGTTAACGCTTTTGATCCGTTGCTGACGCTCCCTATCGCTAACCTCggct                                                                                           >  2:489612/1‑248 (MQ=255)                                                                                                                                         gTTGACCAAGTAGCCTTTAACTCGAGCAGCTGACTTCCAGCCGTTTACCCCAGTCAGGTGGACGACTGACATAGTCATCATCTCTGTCGCTGAAAGAATTTCGCAACGCCTCATGCAGGCGGTGCAATTCCGTTATATC‑AACCTTTTCTGCCGCTTCAATCGCCCGTTGCGCCAACCAGTTGCGCAATACCAGAGCGGGGTTAACGCTTTGCATCAGTTGCTGACGCTCACTATCGATAACCTCGTCttgt                                                                                        >  1:168509/1‑251 (MQ=255)                                                                                                                                                           aaCTCGAGCAGCTGACTTCCAGCCGTTTACCCCAGTCAGGCGGACGACTGACATAGTCATCATCTCTGTCGCTGAAAGAATTTCGCAACGCCTCATGCAGGCGGTGCAATTCCGTCATCTC‑ACCCTTTTCTTCCGCTTCAATCGCCCGTTGCGCCCACCCGTTGCGCAATACCAGAGCGGGGTTAACGCTTTGCATCAGTTGCGTACGCTTCCTATCGCTAACCTCGTCTTGTTGCAAACGCCCCCGaca                                                                       >  2:54747/1‑248 (MQ=255)                                                                                                                                                           aaCTCGAGCAGCTGACTTCCAGCAGTTTACCCCAGTCAGGTGGACGACGGACATAGTCATCATCTCTGTCGCTGAAAGCATTTCGCAACGCCTCATGAAGGCGGTGCAATTCCGTCATATC‑ACCCTTTTCAGCCGCTTCAATCGCCCGTTGAGCCAACCAGTTGCGCAATACCAGAGCGGGGTTAACGCTTTGCATCAGTTGCTGACGCTCACTATCGCTAACCTCGCCTTGTTGCAAACGCACCCGATa                                                                       >  1:74303/1‑250 (MQ=255)                                                                                                                                                             cTCGAGCAGCTGACTTCCAGCCGTTTACCCCAGTCAGGTGGACGACTGACATAGTCATCATCTCTGTCGCTGAAAGAATTTCGCAACGCCTCATGCAGGCGGTGCAATTCCGTCATATC‑ACCCTTTTCTTCCGCTTCAATAGCCCGTTGCGCCAACCAGTTGCGCAATACCAGAGCGGGGTTAACGCTTTGCATCAGGTGCTGACGCTCCCTATCGCTCACCCCGACTTGTTGCAAACGCCCCCGATaa                                                                      >  2:487405/1‑249 (MQ=255)                                                                                                                                                               cGAGCAGCTGACTTCCAGCCGTTTACCCCAGTCAGGTGGACGACTGACATAGTCATCATCTCTGTCGCTGAAAGAATTTCGCAACGCCTCAGGCAGGCGGTGCAATTCCGTCATATC‑ACCCTTTTCTGCCGCTTCAATCGCCCGTTGCGCCAACCAGTTGCGCAATGCCAGAGCGGGGTTAACGCTTTGCGTCAGTTGCTGACGCTCACTAGCGCTAACCACGTCTTGTTGCAAACGCCCCCGATAACggg                                                                  <  1:241013/251‑1 (MQ=255)                                                                                                                                                                       tGACTTCCAGCCGTTTACCCCAGTCAGGTGGACGACTGACATAGTCATCATCTCTGTCGCTGAAAGAATTTCGCAACGCCTCATGCAGGCGGTGCAATTCCGTCATATC‑ACCCTTTTCTGCCGCTTCAATCGCCCGTTGCGCCAACCAGTTTCGCAATACCAGAGCGGGGTTAACGCTTTGCATCAGTTGCTGACGCTCACTATCGCCAACCTCGTCTTGTTGCAAACCCCCCCGATAACGAGAAAACCa                                                           >  2:24186/1‑250 (MQ=255)                                                                                                                                                                               aGCCGTTTACCCCAGTCAGGTGGACGACTGACATAGTCATCATCTCTGTCGCTGAAAGAATTTCGCAACGCCTCATGCAGGCGGTGCAATTCCGTCATATC‑ACCCTTTTCTGCCGCTTCAATCGCCCGTTGCGCCAACCAGTTGCGCAATACCAGAGCGGGGTTAACGCTTTGCATCAGTTGCTGACGCTCACTATCGCTAACCTCGTCTGGTTGCAAACGCCGCCGATAACGGGCAAACCAGTCATCaaa                                                  <  1:571102/251‑1 (MQ=255)                                                                                                                                                                                        ccccAGTCAGGTGGACGACTGATATAGTAATCATCTCTGTCGCTGAAAGAATTTCGCACCGCCTCATGCAGGCGGTGCAATTCCGTCATATC‑CCCCTTTTCTGCCGCTTCAATCGCCCGTTGCGCCAACCAGTTGCGCAATACCAGAGCGGGGTTAACGCTTTGCATCAGTTGCTGACGCTCACTATCGCTAACCTCGTCTTGTTGCAAACGCCCCCGATAACGGGCAAACCAGTCATCAAATGCCGCAC‑g                                        <  2:355741/251‑1 (MQ=255)                                                                                                                                                                                          ccAGTCAGGTGGACGACTGACATAGTCATCATCTCTGTCGCTGAAAGAATTTCGCAACGCCTCATGCAGGCGGTGCAATTCCGTCATATC‑ACCCTTTTCTGCCGCTTCAATCGCCCGTTGCGCCAACCAGTTGCGCAATACCAGAGCGGGGTTAACGCTTTGCATCAGTTGCTGACGCTCACTATCGCTAACCTCGTCTTGTTGCAAACGCCCCCGATAAAGGGCAAACCAGTCATCAAATGCCGCAC‑GAt                                      <  1:427230/251‑1 (MQ=255)                                                                                                                                                                                                aGGTGGACGACTGACATAGTCATCATCTCTGTCGCTGAAAGAATTTCGCAACGCCTCATGCAGGCGGTGCAATTCCGTCATATC‑ACCCTTTTCTGCCGCCTCAATCGCCCGTTGCGCCAACCAGTTGCGCAATACCAGAGCGGGGTTAACGCTTTGCATCAGTTGCTGACGCTCACTATCGCTAACCTCGTCTTGTTGCAAACGCCCCCGATAACGGGCAAACCAGTCATCAAATGCCCCAC‑GATCAATa                                 >  1:83152/1‑250 (MQ=255)                                                                                                                                                                                                       cgacTGACATAGTCATCATCTCTGTCGCTGAAAGAATTTCGCAACGCCTCATGCAGGCGGTGCCATTCCTTCAAATC‑ACCCTTTTTTGCCGCTTCAATCGCCCGTTGCGCCCACCCGTTTCGCAATACCCGAGCGGGGTTAACGCTTTGCATCAGTTGCTGACGCTCACTTACGATAACCTCGTCTTGTTTCAAACCCCCCCGATAACGAGCAAACCAGTCATCAAATGCC‑CACCGCTCAATAAACtcatc                         >  2:124478/1‑251 (MQ=255)                                                                                                                                                                                                          cTGACATAGTCATCATCTCTGTCGCTGAAAGAATTTCGCAACGCCTCATGCAGGCGGTGCAATTCCGTCATATC‑AACCTTTTCTGCCGCTTCAATCGCCCGTTGCGCCAACCAGTTGCGCAATAACAGAGCGGGGTTAACGCTTTGCATCAGTTGCTGGCGCTCACTATCGCCAACCTCGTCTTGTTGCCAACGCCCCCGATAACGGGCAAACCAGTCATCAAATGCCGCAC‑GATCAATAAACCCATCAcc                      >  2:547823/1‑250 (MQ=255)                                                                                                                                                                                                                 aGTCATCATCTCTGTCGCTGAAAGAATTTCGCAACGCCTCATGCAGGCGGTGCAATTCCGTCATATC‑ACCCTTTTCTGCCGCTTCAATCGCCCGTTGCGCCAACCAGTTGCGCAATACCAGAGCGGGGTTAACGCTTTGCATCAGTTGCTGACGCTCACTATCGCTAACCTCGTCTTGTTGCAAACGccc                                                                             <  1:642608/190‑1 (MQ=255)                                                                                                                                                                                                                 aGTCATCATCTCTGTCGCTGAAAGAATTTCGCAACGCCTCATGCAGGCGGTGCAATTCCGTCATATC‑ACCCTTTTCTGCCGCTTCAATCGCCCGTTGCGCCAACCAGTTGCGCAATACCAGAGCGGGGGTAACGCTTTGCCTCAGTTGCTGACGCTCACTATCGCTAACCTCGTCTTGTTGCAAACGccc                                                                             >  2:642608/1‑190 (MQ=255)                                                                                                                                                                                                                  gTCATCAGCTCTGTCACCGAAAGAATTTCGCAACTAATCATGCAGGCGGTGCAATTCCGTCATATC‑ACCCATTGCTGCCGCTTCAATAGCCCGTTGCCCCAACCAGTTGCACAATACCAGAGCGTGGTTAACGCTTTGCATCAGTTGCTGACTCTCACTATCGCTAACCTAGTCGTGTTGCAAACGCCCCCGATAACGGGCAAACCAGTCATCAAATGCCGCTC‑GATCAATAAACTCATCACGTAGCGGTg              <  2:612297/251‑1 (MQ=255)                                                                                                                                                                                                                    catcatCTCTGTCGCTGAAAGAATTTCGCAACGCCTCATGCAGGCGGTGCAATTCCGTCATATC‑ACCCTTTTCTGCCGCTTCAATCGCCCGTTGCGCCAACCAGTTGCGCAATACCAGAGCGGGGTTAACGCTTTGCATCAGTTGCTGACGCTCACTATCGCTAACCTCGTCTTGTTGCAAACGCCCCCGATAACGGGCAAACCAGTCATCAAATGCCGCAC‑GATCAATAAACTCATCACGTAGCGGTGa             <  1:176706/250‑1 (MQ=255)                                                                                                                                                                                                                          ctctGTCGCTGAAAGAATTTCGCAACGCCTCATGCAGGCGGTGCAATTCAGTCATATC‑ACCCTTTTCTGCCGCTTCAATCGCCCGTTGCGCCAACCAGTTGCGCAATACCAGAGCGGGGTTAACGCTTTGCATCAGTTGCTGGAGCTCACTATCGCTAACCTCGGCTTGTTGCAAACGCCCCCGATAACGGGCAAACCAGTCATCAAATGCCGCAC‑GATCAATAAACTCATCAAGTAACGGGGACGCcgcg      >  1:177358/1‑251 (MQ=255)                                                                                                                                                                                                                             tGTCGCTGAAAGAATTTCGCAACGCCTCATGCAGGCGGTGCAATTCCGTCATATC‑ACCCCTTTCTGCCGCTTCAATCGCCCGTTGCGCCAACCAGTTGCGCAATACCAACGCGGGGTTAACGCCTTGCATCAGTTGCTGCCGCTCAATATCGCTAACATCGTCTTGTTGCCAACGCCCCCGATAACGGGCACACCAGTCCTCAAATGCCGCAC‑GATCAATAAACTCATAACGTAGCGGTGACGCCGCGCtc   >  1:414585/1‑250 (MQ=255)                                                                                                                                                                                                                              gTCGCTGAAAGAATTTCGCAACGCCTCATGCAGGCGGTGCAATTCCGTCATATC‑ACCCTTTTCTGCCGCTTCAATCGCCCGTTGCGCCAACCAGTTGCGCAATACCAGAGCGGGGTTAACGCTTTGCATCAGTTGCTGACGCTCACTATCGCTAACCTCGTCTTGTTGCAAACGCCCCCGATAACGGGCAAACCAGTCATCAAATGCCGCAC‑GATCAATAAACTCATCACGTAGCGGTGACGCCGCGCtgt  <  1:378605/251‑1 (MQ=255)                                                                                                                                                                                                                                         |                                                                                                                                                                                                                                                  CCAATCACAAGACCGCTTCCCCGCAGCCTATTCGGCGAATTTCCAGCCAGACACTGTTAGGTCCGGATGGCAAACTGATTATCGATCATGACGGGCAAGAATATCTGCTCCGTAAAACTCAGGCTGGCAAGCTGCTGTTGACCAAGTAGCCTTTAACTCGAGCAGCTGACTTCCAGCCGTTTACCCCAGTCAGGTGGACGACTGACATAGTCATCATCTCTGTCGCTGAAAGGATTTCGCAACGCCTCATGCAGGCGGTGCAATTCCGTCATATC‑ACCCTTTTCTGCCGCTTCAATCGCCCGTTGCGCCAACCAGTTGCGCAATACCAGAGCGGGGTTAACGCTTTGCATCAGTTGCTGACGCTCACTATCGCTAACCTCGTCTTGTTGCAAACGCCCCCGATAACGGGCAAACCAGTCATCAAATGCCGCAC‑GATCAATAAACTCATCACGTAGCGGTGACGCCGCGCTGT  >  NC_000913/1789656‑1790127

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