Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	NC_000913	3,021,885	C→T	C190C (TGC→TGT)	xdhD →	fused putative xanthine/hypoxanthine oxidase: molybdopterin‑binding subunit and Fe‑S binding subunit

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	3,021,885	0	C	T	98.4%	204.5 / NA	62	C190C (TGC→TGT)	xdhD	fused putative xanthine/hypoxanthine oxidase: molybdopterin‑binding subunit and Fe‑S binding subunit
Reads supporting (aligned to +/- strand):  ref base C (0/0);  major base T (25/36);  minor base G (1/0);  total (26/36)
Fisher's exact test for biased strand distribution p-value = 4.19e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

CACCGATCTTCTCGATCGCATCGCCGCACCTACCCGCGAAGAGATCGACGACGCGCTTTCTGGTTTGTTCAGCCGCGATGCTGGCTGGCAGCAATACTATCAGGTCATTGAACTGGCGGTTGCACGTAAAAATAATCCGCAGGCCACCATTGATATCGCTCCGACTTT‑CCGTGACGACCTAGAAGTCATTGGTAAGCATTATCCTAAAACTGATGCCGCGAAAATGGTGCAGGCGAAACCCTGCTATGTTGAAGACCGCGTAACGGCTGACGCCTGCGTCATTAAAATGTTACGTAGCCCACACGCTCACGCACTGATTACTCATCTGGATGTCAGCAAAGCTGAAGCCTTACCGGGCGTCGTTCACGTTATTACTCACCTGAATTGCCCGGATATCTACTATACCCCGGGTGGTCAGAGCGCACCGGAACCGTCACCGCTTGACCGCCGTATGTTCGGCAAGAAAATGCGTCACGTCGGCGATC  >  NC_000913/3021642‑3022126                                                                                                                                                                                                                                                     |                                                                                                                                                                                                                                                  cACCGATCTTCTCGATCGCATCGCCGCACCTACCCGCGAAGAGATCGACGACGCGCTTTCTGGTTTGTTCAGCCGCGATGCTGGCTGGCAGCAATACTATCAGGTCATTGAACTGGCGGTTGCACGTAAAAATAATCCGCAGGCCACCATTGATATCGCTCCGACTTT‑CCGTGACGACCTAGAAGTCATTGGTACGCATTATCCTAAAACTGATGCCGCGAAAATGGTGCAGGCGAAACCCTGTTATGtt                                                                                                                                                                                                                                             <  1:324534/250‑1 (MQ=255)          tctcGATCGCATCGCCGCACCTACCCGCGAAGAGATCGACGACGCGCTTTCTGGTTTGTTCAGCCGCGATGCTGGCTGGCAGCAATACTATCAGGTCATTGAACTGGCGGTTGCACGTAAAAATAATCCGCAGGCCACCATTGATATCGCTCCGACTTT‑CCGTGACGACCTAGAAGTCATTGGTAAGCATTATCCTAAAACTGATGCCGCGAAAATGGTGCAGGCGAAACCCTGTTATGTTGAAGACcgc                                                                                                                                                                                                                                    <  1:389786/250‑1 (MQ=255)              gATCGCATCGCCGCACCTACCCGCGAAGAGATCGACGACGCGCTTTCTGGTTTGTTCAGCCGCGATGCTGGCTGGCAGCAATACTATCAGGTCATTGAACTGGCGGTTGGCCGTAAAAATAATCCGCCGGCCCCCCTTTATATTCCTCCGACTTT‑TCGTGACGACCCTGAAGGCATTGGTAAGCATTTTCCTAAAACTGATGTCCCGGCAATGGTGCAGGCGAAACCCTGTTTTGTTGAAGACCGCGTAAc                                                                                                                                                                                                                               >  1:428678/1‑251 (MQ=255)                tcgcatcgcCGCACCTACCCGCGAAGAGATCGACGAAGCGCTTTCTGGTTTGTTCAGCCGCGATGCTGGCTGGCAGCAATACTATCAGGTCATTGAACTGGCGGTTGCCCGGAAAAATAAACCGCAAGGCACCCTTGATATCGCCCCGACTTTTCCCGGGCG‑CCTAGGAGTCATTGGAAAACATTATCCTAAAACCGATAGCGCGGAAATGGAGCAGGGGAAAACCTGTTTTGTTGAAGAccccggaaagg                                                                                                                                                                                                                             >  1:149121/1‑242 (MQ=255)                      cgccgcACCTACCCGCGAAGAGATCGACGACGCGCTTTCTGGTTTGTTCAGCCGCGATGCTGGCTGGCAGCAATACTATCAGGTCATTGAACTGGCGGTTGCACGTAAAAATAATCCGCAGGCCACCATTGATATCGCTCCGACTTT‑CCGTGACGACCTAGAAGTCATTGGTAAGCATTATCCTAAAACTGATGCCGCGAAAATGGTGCAGGCGAAACCCTGTTATGTTGAAGACCGCGTAACGGCTGAc                                                                                                                                                                                                                        <  1:551392/250‑1 (MQ=255)                           cACCTACCCGCGAAGAGATCGACGACGCGCTTTCTGGTTTGTTCAGCCGCGATGCTGGCTGGCAGCAATACTATCAGGTCATTGAACTGGCGGTTGCACGTAAAAATAATCCGCAGGCCACCATTGATATCGCTCCGACTTT‑CCGTGACGACCTAGAAGTCATTGGTAAGCATTATCCTAAAACTGATGCCGCGAAAATGGTGCAGGCGAAACCCTGTTATGTTGAAGACCGCGTAACGGCTGACGCCTGc                                                                                                                                                                                                                  <  2:484761/251‑1 (MQ=255)                           cACCTACCCGCGAAGAGATCGACGACGCGCTTTCTGGTTTGTTCAGCCGCGATGCTGGCTGGCAGCAATACTATCAGGTCATTGAACTGGCGGTTGCACGTAAAAATAATCCGCAGGCCACCATTGATATCGCTCCGACTTT‑CCGTGACGACCTAGAAGTCATTGGTAAGCATTATCCTAAAACTGATGCCGCGAAAATGGTGCAGGCGAAACCCTGTTATGTTGAAGACCGCGTAACGGCTGACGCCTGc                                                                                                                                                                                                                  <  2:177627/251‑1 (MQ=255)                                 cccGCGAAGAGATCGACGACGCGCTTTCTGGTTTGTTCAGCCGCGATGCTGGCTGGCAGCAATACTATCAGGTCATTGAACTGGCGGTTGCACGTAAAAATAATCCGCAGGCCACCATTGATATCGCTCCGACTTT‑CCGTGACGACCTAGAAGTCATTGGTAAGCATTATCCTAAAACTGATGCCGCGAAAATGGTGCAGGCGAAACCCTGTTATGTTGAAGACCGCGTAACGGCTGACGCCTGCGTCAtt                                                                                                                                                                                                            <  1:489572/251‑1 (MQ=255)                                    gcgAAGAGATCGACGACGCGCTTTCTGGTTTGTTCAGCCGCGATGCTGGCTGGCAGCAATACTATCAGGTCATTGAACTGGCGGTTGCACGTAAAAAAAAACCGCAGGCCCCCATTGATATCGCTCCGACTTT‑CCGTGACGACCTCGAAGTCATTGGTAAGCATTATCCTAAAACTGATGGCGCGAAAAAGGTGCAGGCGAAACCCTGTTATGTTGAAGACCGCGGAACGGCTGGAGCCTGCGTCATTaa                                                                                                                                                                                                          >  1:159567/1‑250 (MQ=255)                                         gagaTCGACGACGCGCTTTCTGGTTTGTTCAGCCGCGATGCTGGCTGGCAGCAATACTATCAGGTCATTGAACTGGCGGTTTCACGGAAAAATAATCCGCAGGCCCCCCTTTATATAGCCCCCCCTTT‑CCGTGACGCCCCACAAATCCTTTGTAAGCATTTTCCTCAAACTGATGCCCCGAAAATAGTGCAGGCGCAACCCTGTTATGTTGAAGACCGCGTAAACGCTGGCGCCCGCGGCATTAAAATGtt                                                                                                                                                                                                    >  1:290198/1‑251 (MQ=255)                                          agaTCGACGACGCGCGTTCTGGTTTGTTCAGCCGCGATGCTGGCTGGCAGCAATACTATCAGGTCATTGAACTGGCGGTTGCACGTAAAAATAATCCGCAGGCCACCATTGATATCGCTCCGACTTT‑CCGTGACGACCTAGAAGTCATTGGTAAGCATTATCCTAAAACTGATGCCGCGAAAATGGTGCAGGCGAAACCCTGTTATGTTGAAGACCGCGTAACGGCTGACGCCTGCGTCATTAAAATGTTa                                                                                                                                                                                                   <  2:428678/251‑1 (MQ=255)                                            aTCGACGACGCGCTTTCTGGTTTGTTCAGCCGCGATGCTGGCTGGCAGCAATACTATCAGGTCATTGAACTGGCGGTTGCACGTAAAAATAATCCGCAGGCCACCATTGATATCGCTCCGACTTT‑CCGTGACGACCTAGAAGTCATTGGTAAGCATTATCCTAAAACTGATGCCGCGAAAATGGTGCAGGCGAAACCCTGTTATGTTGAAGACCGCGTAACGGCTGACGCCTGCGTCATTAAAATGTTACg                                                                                                                                                                                                 <  1:23534/251‑1 (MQ=255)                                                    cgcgCTTTCTGGTTTGTTCAGCCGCGATGCTGGCTGGCAGCAATACTATCAGGTCATTGAACTGGCGGTTGCACGTAAAAATAATCCGCAGGCCACCATTGATATCGCTCCGACTTT‑CCGTGACGACCTAGAAGTCATTGGTAAGCATTATCCTAAAACTGATGCCGCGAAAATGGTGCAGGCGAAACCCTGTTATGTTGAAGACCGCGTAACGGCTGACGCCTGCGTCATTAAAATGTTACGTAGCCcac                                                                                                                                                                                         <  1:85779/251‑1 (MQ=255)                                                          ttCTGGTTTGTTCAGCCGCGATGCTGGCTGGCAGCAATACTATCAGGTCATTGAACTGGCGGTTGCACGTAAAAATAATCCGCAGGCCACCATTGATATCGCTCCGACTTT‑CCGTGACGACCTAGAAGTCATTGGTAAGCATTATCCTAAAACTGATGCCGCGAAAATGGTGCAGGCGAAACCCTGTTATGTTGAAGACCGCGTAACGGCTGACGCCTGCGTCATTAAAATGTTACGTAGCCCACACGc                                                                                                                                                                                     <  2:616905/249‑1 (MQ=255)                                                                 ttGTTCAGCCGCGATGCTGGCTGGCAGCAATACTATCAGGTCATTGAACTGGCGGTTGCACGTAAAAATAATCCGCAGGCCACCATTGATATCGCTCCGACTTT‑CCGTGACGACCTAGAAGTCATTGGTAAGCATTATCCTAAAACTGATGCCGCGAAAATGGTGCAGGCGAAACCCTGTTATGTTGAAGACCGCGTAACGGCTGACGCCTGCGTCATTAAAAt                                                                                                                                                                                                       <  2:43994/224‑1 (MQ=255)                                                                      cggacgcgATGCTGGCTGGCAGCAATACTATCACGTCATTGAACTGGCGGTTGCACGTAAAAATAATCCGCAGGCCACCATTGATATCGCTCCGGCTTT‑CCGTGACGACCTAGAAGTCATTGGTAAGCCTTATCCTAAAACTGATGCCGCGAAAATGGTGCAGGCGAAACCCTGTTTTGTTGAAGACCCCGTAACGGCTGACGCCTTCGTCATTAAAATGTTACCTAGCCCCCACGCTCACTCACTGAtt                                                                                                                                                                        <  2:454379/246‑1 (MQ=255)                                                                        gCCGCGATGCTGGCTGGCAGCAATACTATCAGGTCATTGAACTGGCGGTTGCACGTAAAAATAATCCGCAGGCCACCATTGATATCGCTCCGACTTT‑CCGTGACGACCTAGAAGTCATTGGTAAGCATTATCCTAAAACTGATGCCGCGAAAATGGTGCAGGCGAAACCCTGTTATGTTGAAGACCGCGTAACGGCTGACGCCTGCGTCATTAAAATGTTACGTAGCCCACACGCTCACGCACTGATTACt                                                                                                                                                                     <  2:211200/251‑1 (MQ=255)                                                                             gATGCTGGCTGGCAGCAATACTATCAGGTCATTGAACTGGCGGTTGCACGTAAAAATAATCCGCAGGCCACCATTGATATCGCTCCGACTTT‑CCGTGACGACCTAGAAGTCATTGGTAAGCATTATCCTAAAACTGATGCCGCGAAAATGGTGCAGGCGAAACCCTGTTATGTTGAAGACCGCGTAACGGCTGACGCCTGCGTCATTAAAATGTTACGTAGCCCACACGCTCACGCACTGATTACTCATCt                                                                                                                                                                <  2:186395/251‑1 (MQ=255)                                                                                  tggctggcAGCAATACTATCAGGTCATTGAACTGGCGGTTTCACGGAAAAAAAATCCCCAGGGCCCCATTTGTATTGCTCCCACTTTTCCG‑GACGAACTAGAAATCAATGGGAAGCATTATCCTAAAACTTATGCCGCGAAAATGGGGGAGGCGAAACCCCGGTATTTTGAAAACCGCGGAAAGGGTGACACCTGCGTCATTAAAAAGTTTCGTAGCCCCCCCGCTCCCGCACTGATTACTCATCTGGATg                                                                                                                                                           >  1:136041/1‑251 (MQ=255)                                                                                            caATACTATCAGGTCATTGAACTGGCGGTTGCACGTAAAAATAATCCGCAGGCCACCATTGATATCGCTCCGACTTT‑CCGTGACGACCTAGAAGTCAGTGGTAAGCATTATCCTAAAACTGATGCCGCGAAAATGGTGCAGGCGAAACCCTGTTATGTTGAAGACCGCGTAACGGCTGACGCCTGCGTCATTAAAATGTTACGTAGCCCACACGCTCACGCACTGATTACTCATCTGGATGTCAGCAAAGc                                                                                                                                                 <  2:23023/251‑1 (MQ=255)                                                                                                aCTATCAGGTCATTGAACTGGCGGTTGCACGTAAAAATAATCCGCAGGCCACCATTGATATCGCTCCGACTTT‑CCGTGACGACCTAGAAGTCATTGGTAAGCATTATCCTAAAACTGATGCCGCGAAAATGGTGCAGGCGAAACCCTGTTATGTTGAAGACCGCGTAACGGCTGACGCCTGCGTCATTAAAATGTTACGTAGCCCACACGCTCACGCACTGATTACTCATCTGGATGTCAGCAAAGCTGaa                                                                                                                                             <  2:105048/251‑1 (MQ=255)                                                                                                   aTCAGGTCATTGAACTGGCGGTTGCACGTAAAAATAATCCGCAGGCCACCATTGATATCGCTCCGACTTT‑CCGTGACGACCTAGAAGTCATTGGTAAGCATTATCCTAAAACTGATGCCGCGAAAATGGTGCAGGCGAAACCCTGTTATGTTGAAGACCGCGTAACGGCTGACGCCTGCGTCATTAAAATGTTACGTAGCCCACACGCTCACGCACTGATTACTCATCTGGATGTCAGCAAAGCTGAAg                                                                                                                                            <  1:299671/249‑1 (MQ=255)                                                                                                   aTCAGGTCATTGAACTGGCGGTTGCACGTAAAAATAATCCGCAGGCCACCATTGATAACGCTCCGACTTT‑CCGTGACGACCTAGAAGTCATTGGTAAGCATTATCCTAAAACTGATGCCGCGAAAATGGTGCAGGCGAAACCCTGTTATGTTGAAGACCGCGTCACGGCTGGCGCCTGCGTCATTAAAATGTTACGTAGCCCACACGCTCACGCACTGATTACTCATCTGGATGTCAGcgacgctgtagc                                                                                                                                           <  2:386539/250‑12 (MQ=255)                                                                                                          cANTGAACTGGCGGTTGCACGTAAAAATAATCCGCAGGCCACCATTGATATCGCTCCGACTTT‑CCGTGACGACCTAGAAGTCATTGGTAAGCATTATCCTAAAACTGATGCCGCGAAAATGGTGCAGTCGAAACCCTGTTATGTTGAAGACCGCGTAACGGCTGACGCCTGCGTCATTAAAATGTTACGTAGCCCACACGCTCACGCACTGATTACTCATCTGTATGTCAGCAAAGCTGAAGCCTTACCg                                                                                                                                    <  2:603414/250‑1 (MQ=255)                                                                                                                  tGGCGGTTGCACGCAAAAATAATCCGCAGGCCACCATTGATATCGCTCCGACTTT‑CCGTGACGACCTAGAAGTCATTGGTAAGCATTATCCTAAAACTGATGCCGCGAAAAAGGTGCAGGCGAAACCCTGTTATGTTGAAGACCGCGTAACGGCTGACGCCTGCGTCATTAAAATGTTACGTAGCCCACACGCTCACGCACTGATTACTCATCTGGATGTCAGCAAAGCTGAAGCCTTACCGGGCGTCGtt                                                                                                                           >  2:651411/1‑251 (MQ=255)                                                                                                                     cGGTTGCACGTAAAAATAATCCGCAGGCCACCATTGATATCGCTCCGACTTT‑CCGTGACGACCTAGAAGTCATTGGTAAGCATTATCCTAAAACTGATGCCGCGAAAATGGTGCAGGCGAAACCCTGTTATGTTGAAGACCGCGTAACGGCTGACGCCTGCGTCATTAAAATGTTACGTAGCCCACACGCTCACGCACTGATTACTCATCTGGATGTCAGCAAAGCTGAAGCCTTACCGGGCGTCGTTCAc                                                                                                                        <  1:351696/251‑1 (MQ=255)                                                                                                                         tGCACGTAAAAATAATCCGCAGGCAACCATTGATATCGCTCCGACTTT‑CCGTGACGACCTAGAAGTCATTGGTAAGCATTATCCTAAACCTGATGCCGCGAAAATGGTGCAGGCGAAACCCTGTTATGTTGAAGACCGCGTAACGGCTGACGCCTGCGTCATTAAAATGTTACGTAGCCCACACGCTCACGCACTGATTACTCATCTGGATGTCAGCAAAGCTGAAGCCTTACCGGGCGTCGTTCACGtt                                                                                                                     >  2:59523/1‑250 (MQ=255)                                                                                                                         tGCACGTAAAAAGAATCCGCAGGCCACCATTGATATCGCTCCGACTTT‑CCGTGACGACCTAGAAGTCATTGGTAAGCATTATCCTAAAACTGATGCCGCGAAAATGGTGCAGGCGAAACCCTGTTATGTTGAAGACCGCGTAACGGCTGACGCCTGCGTCATTAAAATGTTACGTAGCCCACACGCTCACGCACTGATTACTCATCTGGATGTCAGCAAAGCTGAAGCCTTACCGGGCGTCGTTCACGtt                                                                                                                     <  1:108485/250‑1 (MQ=255)                                                                                                                           cACGTAAAAATAATCCGCAGGCCACCATTGATATCGCTCCGACATT‑CCGTGACGACATAGAAGTCATTGGTAAGCATTATCCTAAAACTGATGCCGAGAAAATGGTGCAGGCGAAACCCTGTTATGTTGAAGACCGCGTAACGGCTGACGCCTGCGTCATTAAAATGTTACGTAGCACACACGCTCACGCACTGATTAATCATCTGGATGTCAGCAAAGCTGAAGGCTTACCGGGCGTCGTTCGCGttatt                                                                                                                  <  1:451161/251‑1 (MQ=255)                                                                                                                             cGTAAAAATAATCCGCAGGCCACCATTGATATCGCTCCGACTTT‑CCGTGACGACCTAGAAGTCATTGGTAAGCATTATCCTAAAACTGATGCCGCGAAAATGGTGCAGGCGAAACCCTGTTATGTTGAAGACCGCGTAACGGCTGACGCCTGCGTCATTAAAATGTTACGTAGCCCACACGCTCACGCACTGATTACTCATCTGGATGTCAGCAAAGCTGAAGCCTTACCGGGCGTCGTTCACGTTATTAc                                                                                                                >  2:371104/1‑251 (MQ=255)                                                                                                                                aaaaaTAATCCGCAGGCCACCATTGATATCGCTCCGACTTT‑CCGTGACGACCTAGAAGTCATTGGTAAGCATTATCCTAAAACTGATGCCGCGAAAATGGTGCAGGCGAAACCCTGTTATGTTGAAGACCGCGTAACGGCTGACGCCTGCGTCATTAAAATGTTACGTAGCCCACACGCTCACGCACTGATTACTCATCTGGATGTCAGCAAAGCTGAAGCCTTACCGGGCGTCGTTCACGTTATTACTc                                                                                                              >  2:115265/1‑250 (MQ=255)                                                                                                                                  aaaTAATCCGCAGGCCACCATTGATATCGCTCCGACTTT‑CCGTGACGACCTAGAAGTCATTGGTAAGCATTATCCTAAAACTGATGCCGCGAAAATGGTGCAGGCGAAACCCTGTTATGTTGAAGACCGCGTAACGGCTGACGCCTGCGTCATTAAAATGTTACGTAGCCCACACGCTCACGCACTGATTACTCATCTGGATTTCAGCAAAGCTGAAGCCTTACCGGGCGTCGTTCACGTTATTACTCAcc                                                                                                           <  1:372541/251‑1 (MQ=255)                                                                                                                                  aaaTAATCCGCAGGCCACCATTGATATCGCTCCGACTTT‑CCGTGACGACCTAGAAGTCATTGGTAAGCATTATCCTAAAACTGATGCCGCGAAAATGGTGCAGGCGAAACCCTGTTATGTTGAAGACCGCGTAACGGCTGACGCCTGCGTCATTAAAATGTTACGTAGCCCACACGCTCACGCACTGATTACTCATCTGGATGTCAGCAAAGCTGAAGCCTTACCGGGCGTCGTTCACGTTATTACTCAcc                                                                                                           <  1:649489/251‑1 (MQ=255)                                                                                                                                  aaaTAATCCGCAGGCCACCATTGATATCGCTACGACTTT‑CCGTGACGACCTAGAAGTCATTGGTAAGCATTATCCTAAAACTGATGCCGCGAAAATGGTGCAGGCGAAACCCTGTTATGTTGAAGACCGCGTAACGGCTGACGCCTGCGTCATTAAAATGTTACGTAGCCCACACGCTCACGCACTGATTACTCATCTGGATGTCAGCAAAGCTGAAGCCTTACCGGGCGTCGTTCACGTTATTACTCAcc                                                                                                           >  1:251749/1‑251 (MQ=255)                                                                                                                                         ccGCAGGCCACCATTGATATCGCTCCGACTTT‑CCGTGACGACCTAGAAGTCATTGGTAAGCATTATCCTAAAACTGATGCCGCGAAAATGGTGCAGTCGAAACCCTGTTATGTTGAAGACCGCGTAACGGCTGACGCCTGCGTCATTAAAATGTTACGTAGCCCACACGCTCACGCACTGATTACTCATCTGGATGTCAGCAAAGCTGAAGCCTTACCGTGCGTCGTTCACGTTATTACTCACCTGAATTg                                                                                                    <  1:352676/251‑1 (MQ=255)                                                                                                                                         ccGCAGGCCACCATTGATATCGCTCCGACTTT‑CCGTGACGACCTAGAAGTCATTGGTAAGCATTATCCTAAAACTGATGCCGCGAAAATGGTGCAGGCGAAACCCTGTTATGTTGAAGACCGCGTAACGGCTGACGCCTGCGTCATTAAAATGTTACGTAGCCCACACGCTCACGCACTGATTACTCATCTGGATGTCAGCAAAGCTGAAGCCTTACCGGGCGTCGTTCACGTTATTACTCACCTGAATTg                                                                                                    <  2:156233/251‑1 (MQ=255)                                                                                                                                                ccaccaTTGATATCGCTCCGACTTT‑CCGTGACGACCTAGAAGTCATTGGTAAGCATTATCCTAAAACTGATGCCGCGAAAATGGTGCAGGCGAAACCCTGTTATGTTGAAGACCGCGTAACGGCTGACGCCTGCGTCATTAAAATGTTACGTAGCCCACACGCTCACGCACTGATTACTCATCTGGATGTCAGCAAAGCTGAAGCCTTACCGGGCGTCGTTCACGTTATTACTCACCTGAATTGCCCGGat                                                                                             >  1:185888/1‑251 (MQ=255)                                                                                                                                                  accaTTGATATCGCTCCGACTTT‑CCGTGACGACCTAGAAGTCATTGGTAAGCATTATCCTAAAACTGATGCCGCGAAAATGGTGCAGGCGAAACCCTGTTATGTTGAAGACCGCGTAACGGCTGACGCCTGCGTCATTAAAATGTTACGTAGCCCACACGCTCACGCACTGATTACTCATCTGGATGTCAGCAAAGCTGAAGCCTTACCGGGCGTCGTTCACGTTATTACTCACCTGAATTGCCCGGatat                                                                                           >  1:295010/1‑251 (MQ=255)                                                                                                                                                           atCGCTCCGACTTT‑CCGTGACGACCTAGAAGTCATTGGTAAGCATTATCCTAAAACTGATGCCGCGAAAATGGTGCAGGCGAAACCCTGTTATGTTGAAGACCGCGTAACGGCTGACGCCTGCGTCATTAAAATGTTACGTAGCCCACACGCTCACGCACTGATTACTCATCTGGATGTCAGCAAAGCTGAAGCCTTACCGGGCGTCGTTCACGTTATTACTCACCTGAATTGCCCGGATATCTACTATAc                                                                                  >  1:659878/1‑251 (MQ=255)                                                                                                                                                             cGCTCCGACTTT‑CCGTGACGACCTAGAAGTCATTGGTAAGCATTATCCTAAAACTGATGCCGCGAAAATGGTGCAGGCGAAACCCTGTTATGTTGAAGACCGCGTAACGGCTGACGCCTGCGTCATTAAAATGTTACGTAGCCCACACGCTCACGCACTGATTACTCATCTGGATGTCAGCAAAGCTGAAGCCTTACCGGGCGTCGTTCACGTTATTACTCACCTGAATTGCCCGGATATCTACTATAccc                                                                                <  2:185888/251‑1 (MQ=255)                                                                                                                                                              gCTCCGACTTT‑CCGTGACGACCTAGAAGTCATTGGTAAGCATTATCCTAAAACTGATGCCGCGAAAATGGTGCAGGCGAAACCCTGTTATGTTGAAGACCGCGTAACGGATGACGCCTGCGTCATTAAAATGTTACGTAGCCCACACGCTCACGCACTGATTACTCCTCTGGATGTCAGCAAAGCTGAAGCCTTACCGGGCGTCGTTCACGTTATTACTCACCTGAATTGCCCGGATATCTACTATAcccc                                                                               >  1:619297/1‑251 (MQ=255)                                                                                                                                                                 ccGACTTT‑CCGTGACGACCTAGAAGTCATTGGTAAGCATTATCCTAAAACTGATGCCGCGAAAATGGTGCAGGCGAAACCCTGTTATGTTGAAGACCGCGTAACGGCTGACGCCTGCGTCATTAAAATGTTACGTAGCCCACACGCTCACGCACTGATTACTCATCTGTATGTCAGCAAAGCTGAAGCCTTACCGGGCGTCGTGCACGTTATTACTCACCTGAATTGCCCGTATATCTACTATACCCCg                                                                              <  2:644550/249‑1 (MQ=255)                                                                                                                                                                            gTGACGTCCTATAAGTCATTGGTAAGCATTATCCTAAAACTGATGCCGCGAAAATGGTGCAGGCGAAACCCTGTTATGTTGAAGACCGCGTAACGGCTGACGCCTGCGTCATTAAAATGTTACGTAGCCCACACGCTCACGCACTGATTACGCATCTGGATGTCAGCAAAGCTGAAGCCTTACCGGGCGTCGTTCACGTTATTACTCACCTGAATTGCCCGGATATCTACTATACCCCGGGTGGTCAGAgc                                                                  <  2:251749/251‑1 (MQ=255)                                                                                                                                                                              gacgacCTAGAAGTCATTGGTAAGCATTATCCTAAAACTGATGCCGCGAAAATGGTGCAGGCGAAACCCTGTTATGTTGAAGACCGCGTAACGGCTGACGCCTGCGTCATTAAAATGTTACGTAGCCCACACGCTCACGCACTGATTACTCATCTGGATGTCAGCAAAGCTGAAGCCTTACCGGGCGTCGTTCACGTTATTACTCTCCTGAATTGCCCGGATATCTACTATACCCCGGGTGGTCAGAgcgc                                                                <  1:257969/251‑1 (MQ=255)                                                                                                                                                                                cgacCTAGAAGTCATTGGTAAGCATTATCCTAAAACTGATGCCGCGAAAATGGTGCAGGCGAAACCCTGTTATGTTGAAGACCGCGTAACGGCTGACGCCTGCGTCATTAAAATGTTACGTAGCCCACACGCTCACGCACTGATTACTCATCTGGATGTCAGCAAAGCTGAAGCCTTACCGGGCGTCGTTCACGTTATTACTCACCTGAATTGCCCGGATATCTACTATACCCCGGGTGGTCAGAGCGCAc                                                              >  1:482265/1‑251 (MQ=255)                                                                                                                                                                                 gacCTAGAAGTCATTGGTAAGCATTATCCTAAAACTGATGCCGCGAAAATGGTGCAGGCGAAACCCTGTTATGTTGAAGACCGCGTAACGGCTGACGCCTGCGTCATTAAAATGTTACGTAGCCCACACGCTCACGCACTGATTACTCATCTGGATGTCAGCAAAGCTGAAGCCTTACCGGGCGTCGTTCACGTTATTACTCACCTGAATTGCCCGGATATCTACTATACCCCGGGTGGTCAGAGCGCAcc                                                             >  2:200422/1‑251 (MQ=255)                                                                                                                                                                                                          aTCCTAAAACTGATGCCGCGAAAATGGTGCAGGCGAAACCCTGTTATGTTGAAGACCGCGTAACGGCTGACGCCTGCGTCATTAAAATGTTACGTAGCCCACACGCTCACGCACTGATTACTCATCTGGATGTCAGCAAAGCTGAAGCCTTACCGGGCGTCGTTCACGTTATTACTCACCTGAATTGCCCGGATATCTACTATACCCCGGGTGGTCAGAGCGCACCGGAACCGTCACCGCTTGACCGCCGt                                    >  1:639954/1‑251 (MQ=255)                                                                                                                                                                                                          aTCCTAAAACTGATGCCGCGAAAATGGTGCAGGCGAAACCCTGTTATGTTGAAGACCGCGTAACGGCTGACGCCTGCGTCATTAAAATGTTACGTAGCCCACACGCTCACGCACTGATTACTCATCTGGATGTCAGCAAAGCTGAAGCCTTACCGGGCGTCGTTCACGTTATTACTCACCTGAATTGCCCGGATATCTACTATACCCCGGGTGGTCAGAGCGCACCGGAACCGTCACCGCTTGACCGCCGt                                    >  2:160306/1‑251 (MQ=255)                                                                                                                                                                                                              tAAAACTGATGCCGCGAAAATGGAGCAGGAGAAACCCTGTTATGTTGAAGACCGCGTAACGGCTGACGCCTGCGTCATTAAAATGTTACGTAGCCCACACGCTCACGCACTGATTACTCATCTGGATGTCAGCAAAGCTGAAGCCTTACCGGGCGTCGTTCACGTTATTACTCACCTGAATTGCCCGGATATCTACTATACCCCGGGTGGTCAGAGCGCACCGGAACCGTCACCGCTTGACCGCCGTATGt                                >  1:573368/1‑251 (MQ=255)                                                                                                                                                                                                                   cTGATGCCGCGAAAATGGTGCAGGCGAAACCCTGTTATGTTGAAGACCGCGTAACGGCTGACGCCTGCGTCATTAAAATGTTACGTAGCCCACACGCTCACGCACTGATTACTCATCTGGAt                                                                                                                                                            >  1:292337/1‑122 (MQ=255)                                                                                                                                                                                                                   cTGATGCCGCGAAAATGGTGCAGGCGAAACCCTGTTATGTTGAAGACCGCGTAACGGCTGACGCCTGCGTCATTAAAATGTTACGTAGCCCACACGCTCACGCACTGATTACTCATCTGGAt                                                                                                                                                            <  2:292337/122‑1 (MQ=255)                                                                                                                                                                                                                     gATGCCGCGAAAATGGTGCAGGCGAAACCCTGTTATGTTGAAGACCGCGTAACGGCTGACGCCTGCGTCATTAAAATGTTACGTAGCCCACACGCTCACGCACTGATTACTCATCTGGATGTCAGCAAAGCTGAAGCCTTACCGGGCGTCGTTCACGTTATTACTCACCTGAATTGCCCGGATATCTACTATCCCCCGGGTGGTCAGAGCGCACCGGAACCGTCCCCGCTTGACCGCCGTATGTTCGGCaa                         <  2:430890/251‑1 (MQ=255)                                                                                                                                                                                                                          cgcgAAAATGGTGCAGGCGAAACCCTGTTATGTTGAAGACCGCGTAACGGCTGACGCTTGCGTCATTAAAATGTTACGTAGCCCACACGCTCACGCACTGATTACTCATCTGGATGTCAGCAAAGCTGAAGCCTTACCGGGCGTCGTTCACGTTATTACTCACCTGAATTGCCCGGATATCTACTATCCCCCGGGTGGTCAGAGCGCACCGGAACCGTCACCGCTTGACCGCCGTATGTTCGGCAAGaaaa                    <  2:131965/251‑1 (MQ=255)                                                                                                                                                                                                                               aaaTGGAGCAGGCGAAACACTGTTATGTTGAAGACCGCGTAACGGCTGACGCCTGCGTCATTAAAATGTTACGTAGCCCACACGCTCACGCACTGATTACTCATCTGGATGTCAGCAAAGCTGAAGCCTTACCGGGCGTCGTTCACGTTATTACTCACCTGAATTGCCCGGATATCTACTATACCCCGGGTGGTCAGAGCGCACCGGAACCGTCACCGCTTGACCGCCGTATGTTCGGCAAGAAAATGCGt               >  1:521772/1‑251 (MQ=255)                                                                                                                                                                                                                                 aTGGTGCAGGCGAAACCATGTTATGTTGAAGACCGCGTAACGGCTGACGCCTGCGTCATTAAAATGTTACGTAGCCCACACGCTCACGCACTGATTACTCATCTGGATGTCAGCAAAGCTGAAGCCTTACCGGGCGTCGTTCACGTTATTACTCACCTGAATTGCCCGGATATCTACTATACCCCGGGTGGTCAGAGCGCACCGGAACCGTCACCGCTTGACCGCCGTATGTTCGGCAAGAAAATGCGTc              >  1:558090/1‑250 (MQ=255)                                                                                                                                                                                                                                     tGCAGGCGAAACCCTGTTATGTTGAAGACCGCGTAACGGCTGACGCCTGCGTCATTAAAATGTTACGTAGCCCACACGCTCACGCACTGATTACTCATCTGGATGTCAGCAAAGCTGAAGCCTTACCGGGCGTCGTTCACGTTATTACTCACCTGAATTGCCCGGATATCTACTATACCCCGGGTGGTCAGAGCGCACCGGAACCGTCACCGCTTGACCGCCGTATGTTCGGCAAGAAAATGCGTCACGTc         >  1:167675/1‑251 (MQ=255)                                                                                                                                                                                                                                       cAGGCGAAACACTGTTATGTTGATGACCGCGTAACGGCTGACGCCTGCGTCATTAAAATGTGACGTATCCCACACGCTCACGCCCTGAGTCCTCATCGGTACGTCAGAAAAGCTGACGCCTTACCGGGCGTCGTCCTCTTTATAACTCACTGGAATGGCCCGGCTATCTCCTATCCCCCGGGCGGACCGATCGCACCGGACCCGTCCTCGCTCGACCGCCGTATGTTCGGCAAGAAAATGCGTCACGTCgg       <  2:353969/251‑1 (MQ=255)                                                                                                                                                                                                                                         ggCGAAACCCTGTTATGTTGAAGACCGCGTAACGGCTGACGCCTGCGTCATTAAAATGTTACGTAGCCCACACGCTCACGCACTGATTACTCATCTGGATGTCAGCAAAGCTGAAGCCTTACCGGGCGTCGTTCACGTTATTACTCACCTGAATTGCCCGGATATCTACTATACCCCGGGTGGTCAGAGCGCACCGGAACCGTCACCGCTTGAcc                                         <  2:620667/215‑1 (MQ=255)                                                                                                                                                                                                                                         ggCGAAACCCTGTTATGTTGAAGACCGCGTAACGGCTGACGCCTGCGTCATTAAAATGTTACGTAGCCCACACGCTCACGCACTGATTACTCATCTGGATGTCAGCAAAGCTGAAGCCTTACCGGGCGTCGTTCACGTTATTACTCACCTGAATTGCCCGGATATCTACTATACCCCGGGTGGTCAGAGCGCACCGGAACCGTCACCGCTTGAcc                                         >  1:620667/1‑215 (MQ=255)                                                                                                                                                                                                                                         ggCGAAACCCTGTTATGTTGAAGACCGCGTAACGGCTGACGCCTGCGTCATTAAAATGTTACGTAGCCCACACGCTCACGCACTGATTACTCATCTGGATGTCAGCAAAGCTGAAGCCTTACCGGGCGTCGTTCACGTTATTACTCACCTGAATTGCCCGGATATCTACTATACCCCGGGTGGTCAGAGCGCACCGGAACCGTCACCGCTTGACCGCCGTAt                                  >  2:115787/1‑222 (MQ=255)                                                                                                                                                                                                                                         ggCGAAACCCTGTTATGTTGAAGACCGCGTAACGGCTGACGCCTGCGTCATTAAAATGTTACGTAGCCCACACGCTCACGCACTGATTACTCATCTGGATGTCAGCAAAGCTGAAGCCTTACCGGGCGTCGTTCACGTTATTACTCACCTGAATTGCCCGGATATCTACTATACCCCGGGTGGTCAGAGCGCACCGGAACCGTCACCGCTTGACCGCCGTAt                                  <  1:115787/222‑1 (MQ=255)                                                                                                                                                                                                                                            gAAACCCTGTTATGTTGAAGACCGCGTAACGGCCGACGCCTGCGTCATTAAAATGTTACGTAGCCCACACGCTCACGCACTGATTACTCATCTGGATGTCAGCAAAGCTGAAGCCTTACCGGGCGTCGTTCACGTTATTACTCACCTGAATTGCCCGGATATCTACTATACCCCGGGTGGTCAGAGCGCACCGGAACCGTCACCGCTTGACCGCCGTATGTTCGGCAAGAAAATGCGTCACGTCGGCGATc  >  1:14371/1‑251 (MQ=255)                                                                                                                                                                                                                                                     |                                                                                                                                                                                                                                                  CACCGATCTTCTCGATCGCATCGCCGCACCTACCCGCGAAGAGATCGACGACGCGCTTTCTGGTTTGTTCAGCCGCGATGCTGGCTGGCAGCAATACTATCAGGTCATTGAACTGGCGGTTGCACGTAAAAATAATCCGCAGGCCACCATTGATATCGCTCCGACTTT‑CCGTGACGACCTAGAAGTCATTGGTAAGCATTATCCTAAAACTGATGCCGCGAAAATGGTGCAGGCGAAACCCTGCTATGTTGAAGACCGCGTAACGGCTGACGCCTGCGTCATTAAAATGTTACGTAGCCCACACGCTCACGCACTGATTACTCATCTGGATGTCAGCAAAGCTGAAGCCTTACCGGGCGTCGTTCACGTTATTACTCACCTGAATTGCCCGGATATCTACTATACCCCGGGTGGTCAGAGCGCACCGGAACCGTCACCGCTTGACCGCCGTATGTTCGGCAAGAAAATGCGTCACGTCGGCGATC  >  NC_000913/3021642‑3022126

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: ‑