Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	NC_000913	2,032,423	G→A	V219V (GTC→GTT)	yedJ ←	putative HD superfamily phosphohydrolase YedJ

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	2,032,423	0	G	A	98.1%	167.0 / NA	53	V219V (GTC→GTT)	yedJ	putative HD superfamily phosphohydrolase YedJ
Reads supporting (aligned to +/- strand):  ref base G (0/0);  major base A (33/19);  minor base C (1/0);  total (34/19)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.74e-01

GCTCCAGTGATTCTCACCTACACCATTAAGCTGCGCCACCAACGTTTTAACATCATAAATTTGCAGCAGTTTCTCCAGCATTGCCTGTGCGGCATTCCCGGTGCTGTATGAATCGGTTACTGATATATTTTCCTGCATAGATTTCACCGGCCATCTGAGATGGCCGACAGGTTAACATAATTACAGCCATTTAGGGCCACGCGCTCAGCCAGCGGATGAAAACGCATCTATCACCTTGTGATCGACACC‑TTCATTCT‑CCCCCGCCAGTTCGGCACTGAGCTTCGCCATAAACTCGACTAAAAAATGCGCGTTGTGCTGCGCCAGCTGCTTGCCTCTTGCGGTTTGCATGGTTTGCGGTAGCTTCAGCAACTTAGTCTGGAAATGATCCAGCGCGTAGCGTTTATCATCAAGCGGGCGATGCTGTGCGAACGGATCTTCACCATCAAACAGCGCCACGCCCAATGCCCCTGAAACGGCAAACACGC  >  NC_000913/2032180‑2032664                                                                                                                                                                                                                                                    |                                                                                                                                                                                                                                                    gcTCCAGTGATTCTCACCTACACCATTAAGCTGCGCCACCAACGTTTTAACATCATAAATTTGCAGCAGTTTCTCCAGCATTGCCTGTGCGGCATTCCCGGTGCTGTATGAATCGGTTACTGATATATTTTCATGCATAGATTTCACCGGCCATCTGAGATGGCCGACAGGTTAACATAATTACAGCCATTTAGGGCAACGCGATCAGCCAGCGGATGAAAACGCATCTATCACCTTGTGATCAAAACC‑tt                                                                                                                                                                                                                                             >  1:318328/1‑251 (MQ=255)    ccAGTGATTCTCACCTACACCATTAAGCTGCGCCACCAACGTTTTAACATCATAAATTTGCAGCAGTTTCTCCAGCATTGCCTGTGCGGCATTCCCGGTGCTGTATGAATCGGTTACTGATATATTTTCCTGCATAGATTTCACCGGCCATCTGAGATGGCCGACAGGTTAACATAATTACAGCCATTTCGGGCCACGCGCTCAGCCAGCGGATGAAAACGCATCTATCCCCTTGTGATCAACACC‑TTCAt                                                                                                                                                                                                                                          >  1:136749/1‑251 (MQ=255)                    acacCATTAAGCTGCGCCACCAACGTTTTAACATCATAAATTTGCAGCAGTTTCTCCAGCATTGCCTGTGCGGCATTCCCGGTGCTGTATGAATCGGTTACTGATATATTTTCCTGCATAGATTTCACCGGCCATCTGAGATGGCCGACAGGTTAACATAATTACAGCCATTTAGGGCCACGCGCTCAGCCAGCGGATGAAAACGCATCTATCACCTTGTGATCAACACC‑TTCATTCT‑CCCCCGCCAGTTc                                                                                                                                                                                                                         <  1:208470/251‑1 (MQ=255)                    acacCATTAAGCTGCGCCACCAACGTTTTAACATCATAAATTTGCAGCAGTTTCTCCAGCATTGCCTGTGCGGCATTCCCGGTGCTGTATGAATCGGTTACTGATATATTTTCCTGCATAGATTTCACCGGCCATCTGAGATGGCCGACAGGTTAACATAATTACAGCCATTTAGGGCCACGCGCTCAGCCAGCGGATGAAAACGCATCTATCACCTTGTGATCAACACC‑TTCATTCT‑CCCCCGCCAGTTc                                                                                                                                                                                                                         >  2:59098/1‑251 (MQ=255)                     cacCATTAAGCTGCGCCACCAACGTTTTAACATCATAAATTTGCAGCAGTTTCTCCAGCATTGCCTGTGCGGTATTCCCGGTGCTGTATGAATCGGTTACTGATATATTTTCCTGCATAGATTTCACCGGCCATCTGAGATGGCCGACAGGTTAACATAATTACAGCCATTTAGGGCCACGCGCTCAGCCAGCGGATGAAAACGCATCTATCACCTTGTGATCAACACC‑TTCATTCT‑CCCCCGCCAGTTCg                                                                                                                                                                                                                        <  2:356460/251‑1 (MQ=255)                      acCATTAAGCTGCGCCACCAACGTTTTAACATCATAAATTTGCAGCAGTTTCTCCAGCATTGCCTGTGCGGAATTCCAGGTGCTGTATGAATCGGTTACTGATATATTTTCCTGCATAGATTTCACCGGCCATCTGAGATGGCCGACAGGTTAACATAATTACAGCCATTTAGGGCCACGCGCTCAGCCAGCGGATGAAAACGCATCTATCACCTTGTGATCAACACA‑TTCATTCT‑CCCCCGCCAGTTCgg                                                                                                                                                                                                                       >  1:636064/1‑251 (MQ=255)                            aaGCTGCGCCACCAACGTTTTAACATCATAAATTTGCAGCAGTTTCTCCAGCATTGCCTGTGCGGCATTCCCGGTGCTGTATGAATCGGTTACTGATATATTTTCCTGCATAGATTTCACCGGCCATCTGAGATGGCCGACAGGTTAACATAATTACAGCCATTTAGGGCCACGCGCTCAGCCAGCGGATGAAAACGCATCTATCACCTTGTGATCAACACC‑TTCATTCT‑CCCCCGCCAGTTCGGCACTg                                                                                                                                                                                                                  <  2:616737/250‑1 (MQ=255)                                 gcgcacccaACGTTTTAACATCATAAATTTGCAGCAGTGTCTCCATCATTGCCTGTGCGGCATCCCCGGTGCTGTATGAATCGGTTACTGATATATTGTCCTGCATAGATTTCACCGGCCATCTGAGATGGCCGACAGGTTAAAATAATTACAGCCATTTAGGGCCACGCGCTCAGCCAGCGGATGAAAACGCATCTATCACCTTGTGATCAACACC‑TTCATTCT‑CCCCCGCCAGTTCGGCACTGAGCTTc                                                                                                                                                                                                            <  2:560368/245‑1 (MQ=255)                                           cgatttAACATCATAAATTTGCAGCAGTTTCTCAAGCATTGCATGTGCGGCATTACCGGAGCTGTATGAATCGGTTACTGATATATTTTCCTGCATAGATTTCACCGGCCATCTGAGATGGCCGACAGGTGAACATACTTACAGCCATTTAGGGCCACGCGCTCAGCCTGCGGCTGAAAACGCCTCTATCACCTTGTGATCAACACC‑TTCATTCT‑CACACACCAGTTCGTCACTGAGCTTCGCCATAAACt                                                                                                                                                                                                  >  2:664902/4‑251 (MQ=255)                                                     tcatAAATTTGCAGCAGTTTCTCCAGCATTGCCTGTGCGGCATTCCCGGTGCTGTATGAATCGGTTACTGATATATTTTCCTGCATAGATTTCACCGGCCATCTGAGATGGCCGACAGGTTAACATAATTACAGCCATTTAGGGCCACGCGCTCAGCCAGCGGATGAAAACGCATCTATCACCTTGTGATCAACACC‑TTCATTCT‑CCCCCGCCAGTTCGGCACTGAGCTTCGCCATAAACTCGACTaaaa                                                                                                                                                                                         >  1:671336/1‑250 (MQ=255)                                                     tcatAAATTTGCAGCAGTTTCTCCAGCATTGCCTGTGCGGCATTCACGGTGCTGTATGAATCGGTTACTGATATATTTTCCTGCATAGATTTCACCGGCCATCTGAGATGGCCGACAGGTTAACCTAATTACAGCCATTTAGGGCCACGCGCTCAGCCAGCGGATGAAAACGCATCTATCACCTTGTGATCAACACC‑TTCATTCT‑CCCCCGCCAGTTCGGCACTGAGCTTCGCCATAAACTCGAATaaaa                                                                                                                                                                                         >  2:310019/1‑250 (MQ=255)                                                            tttGCAGCAGTTTCTCCAGCATTGCCTGTGCGGCATTCCCGGTGCTGTATGAATCGGTTACTGATATATTTTCCTGCATAGATTTCACCGGCCATCTGAGATGGCCGACAGGTTAACATAATTACAGCCATTTAGGGCCACGCGCTCAGCCAGCGGATGAAAACGCATCTATCACCTTGTGATCAACACC‑TTCATTCT‑CCCCCGCCAGTTCGGCACTGAGCTTCGCCATAAACTCGACTAAAAAATGcgcg                                                                                                                                                                                 >  2:102547/1‑251 (MQ=255)                                                                     gTTTCTCCAGCATTGCCTGTGCGGCATTCCCGGTGCTGTATGAATCGGTTACTGATATATTTTCCTGCATAGATTTCCCCGGCCATCTGAGATGGCCGACAGGTTAACATAATTACAGCCATTTAGGGCCACGCGCTCAGCCAGCGGATGAAAACGCATCTATCACCTTGTGATCAACACC‑TTCATTCT‑CCCCCGCCAGTTCGGCACTGAGCTTCGCCATAAACTCGACTAAAAAATGCGCGTTGtgctgc                                                                                                                                                                        <  1:610130/251‑1 (MQ=255)                                                                        tctcCAGCATTGCCTGTGCGGCATTCCCGGTGCTGTATGAATCGGTTACTGATATATTTTCCTGCATAGATTTCACCGGCCATCTGAGATGGCCGACAGGTTAACATAATTACAGCCATTTAGGGCCACGCGCTCAGCCAGCGGATGAAAACGCATCTATCACCTTGTGATCAACACC‑TTCATTCT‑CCCCCGCCAGTTCGGCACTGAGCTTCGCCATAAACTCGACTAAAAAATGCGCGTTGTGCTGCGcc                                                                                                                                                                     >  1:361524/1‑251 (MQ=255)                                                                           ccAGCATTGCCTGTGCTGCATTCCCGGTGCTGTATGAATCGGTTACTGATATATTTTCCTGCATAGATTTCACCGGCCATCTGAGATGGCCGACAGGTTAACAAAATTACAGCCATTTAGGGCCACGCGCTCAGCCAGCGGATGAAAACGCATCTATCCCCTTGTGATCAACACC‑TGCATTCT‑CCCCCGCCAGTTCGGCACTGAGCTTCGCCATAAACTCGACTAAAAAATGCGCGTTGGGCTGCGCCa                                                                                                                                                                    <  1:602132/249‑1 (MQ=255)                                                                                       gtgCGGCATTCCCGGTGCTGTATGAATCGGTTACTGATATATTTTCCTGCATAGATTTCACCGGCCATCTGAGATGGCCGACAGGTTAACATAATTACAGCCATTTAGGGCCACGCGCTCAGCCAGCGGATGAAAACGCATCTATCACCTTGTGATCAACACC‑TTCATTCT‑CCCCCGCCAGTTCGGCACTGAGCTTCGCCATAAACTCGACTAAAAAATGCGCGTTGTGCTGCGCCAGCTGCTTGCCTCtt                                                                                                                                                      >  1:155557/1‑251 (MQ=255)                                                                                       gtgCGGCATTCACGGTGCTGTATGAATCGGTTACTGATATATTTTCCTGCATAGATTTCACCGGCCATCTGCGATGGCCGACAGGTTAACATAATTACAGCCATTTAGGGCCAAGCGCTCAGCCAGCGGATGAAAACGCATCTATCACCTTGTGATCAACACC‑TTCATTCT‑CCCCCGCCAGTTCGGCACTGAGCTTCGCCATAAACTCGACTAAAAAATGCGCGTTGTGCTGCGCCAGCTGCTTGCCTCtt                                                                                                                                                      >  2:271537/1‑251 (MQ=255)                                                                                          cGGCATTCCCGGTGCTGTATGAATCGGTTACTGATATATTTTCCTGCATAGATTTCACCGGCCATCTGAGATGGCCGACAGGTTAACATAATTACAGCCATTTAGGGCCACGCGCTCAGCCAGCGGATGAAAACGCATCTATCCCCTTGTGATCAACCCC‑TTCATTCT‑CCCCCGCCAGTTCGGCACTGAGCTTCGCCATAAACTCGACTAAAAAATGCGCGTTGTGCTGCGCCAGCTGCTTGCCTCTTGCg                                                                                                                                                   <  1:59098/251‑1 (MQ=255)                                                                                           ggCATTCCCGGTGCTGTATGAATCGGTTACTGATATATTTTCCTGCATAGATTTCACCGGCCATCTGAGATGGCCGACAGGTTAACATAATTACAGCCATTTAGGGCCACGCGCTCAGCCAGCGGATGAAAACGCATCTATCACCTTGAGATCAACACC‑TTCATTCT‑CCCCCGCCAGTTCGGCACTGAGCTTCGCCATAAACTCGACTAAAAAATGCGCGTTGTGCTGCGCCAGCTGCTTGCCTCTTGCgg                                                                                                                                                  >  1:287813/1‑251 (MQ=255)                                                                                                tCCCGGTGCTGTATGAATCGGTTACTGATATATTTTCCTGCATAGATTTCACCGGCCATCTGAGATGGCCGACAGGTTAACATAATTACAGCCATTTAGGGCCACGCGCTCAGCCAGCGGATGAAAACGCATCTATCACCTTGTGATCAACACC‑TTCATTCT‑CCCCCGCCAGTTCGGCACTGAGCTTCGCCATAAACTCGACTAAAAAATGCGCGTTGTGCTGCGCCAGCTGCTTGCCTCTTGCGGTTTGc                                                                                                                                             >  1:691936/1‑251 (MQ=255)                                                                                                 cccGGTGCTGTATGAATCGGTTACTGATATATTTTCCTGCATAGCTTTCACCGGCCATCTGAGATGGCCGACCGGTAAACATAATTACAGCCATTTAGGGCCCCGCGCTCAGCCAGCGGATGAAAACGCATCTATCACCTTGTGCTCAACACC‑TTCATTCT‑CCCCCGCCAGTTCAGCACTGAGCTTCGCCATAAACTCGACTAACAAATGCGCGTTGTGCTGCGCCAGCTGCTTGCCTCTTGCGGTTTGc                                                                                                                                             >  1:675240/1‑250 (MQ=255)                                                                                                 cccGGTGCTGTATGAATCGGTTACTGATATATTTTCCTGCATAGATTTCACCGGCCATCTGAGATGGCCGACAGGTTAACATAATTACAGCCATTTAGGGCCACGCGCTCAGCCAGCGGGTGAAAACGCATCTATCACCTTGTGATCAACCCC‑TTCATCCT‑CCCCCGCCAGTTCGGCACTGAGCTTCGCCATAAACTCGAGTAAAAAATGCGCGTTGTGCTGCGCCAGCTGCTTGCCTCTTGCGGTTTGCa                                                                                                                                            <  1:270742/251‑1 (MQ=255)                                                                                                    ggTGCTGTATGAAGCGGTTTCTGAAATATTTTCTTGCATAGATTTCACCGGCCATCGGAGATGGCCGAAAGGTTAACATAATTACAGCCATTTAGGGCCACGCGCTCAGCCAGCGGATGAAAACGCTTCTATCACCTTGTTATCAACCCC‑TTCATTCC‑CCCCCGCCAGTTCGGCACTGAGCTTCGCCATAAACTCGACTAAAAAATGCGCGTTGTGCTGCGCCAGCTGCTTGCCTCTTGCGGTTTGCATgg                                                                                                                                         <  2:505767/251‑1 (MQ=255)                                                                                                      tGCTGTATGAATCGGTTACTGATATCTTTTCCTGAATAGCTTTCACCGGCAATCTGAGATGGCCGAACGTTTTACATAATTACATCCATTTAGGTCCACACGCCCAGCAGGCGGATGAAAACGCCTATATCCCCTTCTGATCAAAACC‑TTCATTCT‑ACCCCGCCAGTTCGGCACTTGGCTTCGCCATACACTCGACTACAAAATGCGCGTTGTGCTGCGCCAGCTGCTTGCCTCATTCGGGTTGCATGGtt                                                                                                                                       >  2:270647/1‑251 (MQ=255)                                                                                                                          gATATATTTTCCTGCATAGATTTCACCGGCCATCTGAGATGGCCGACAGGTTAACATAATTACAGCCATTTAGGGCCACGCGCTCAGCCAGCGGATGAAACCGCATCTATCACCTTGTGATCAACACC‑TGCATTCT‑CCCCCGCCAGTTCGGCACTGAGCTTCGCCATAAACTCGACTAAAAAATGCGCGTTGTGCTGCGCCAGCTGCTTGCCTCTTGCGGTTTGCATGGTTTGCGGTAGCTTCAGCAACt                                                                                                                    <  1:98230/250‑1 (MQ=255)                                                                                                                            tatatTTTCCTGCATAGATTTCACCGGCCATCTGAGATGGCCGACAGGTTAACATAATTACAGCCATTTAGGGCAACGCGCTCAGCCAGCGGATGAAAACGCATCTATCCCCTTGTGATCAACACC‑TTCATTCT‑CCCCCGCCAGTTCGGCACTGAGCTTCGCCATAAACTCGACTAAAAAATGCGCGTTGTGCTGCGCCAGCTGCTTGCCTCTTGCGGTTTGCATGGTTTGCGGTAGCTTCAGCAACtt                                                                                                                   >  1:667088/1‑249 (MQ=255)                                                                                                                                 tttCCTGCAGAGAGTACCCCGGCCATTTGAGATGGCCGACAGGTTAACAAAATTACACCCATTTAGGGCCCCGCGCTCAGCCAGCGGATGAAAACTCATCTATCCCCTTGTGATCAACACC‑TTCATTCC‑CCCCCGCCAGTTCGGCACTGAGCTTCGCCATAAACTCGACTAAAAAATGCGCGTTGTGCTGCGCCAGCGGCTTGCCTCTTGCGGTTTGCATGGTTTGCGGTAGCCTCAGCAACGTAGTCTg                                                                                                             <  2:55737/250‑1 (MQ=255)                                                                                                                                                   ccggacaatgtagaTGGCCGCCAGGTTAACATAATACCAGCCATTTAGGGCCACGCGCTCAGCCAGCGGATGAAACCGCATCTTTCACCTTGTGATCAACACC‑TTCAT‑CTACCCCCGCCAGTTCGTCTCTGAGCTTCGCCATAAACTCGACTAAAAAATGCGCGTTGTGCTGCGCCAGCTGCTTGCCTCTTGCGGTTGGCATGGTTTGCGGTAGCTTCAGCAACTTAGTCTTGAAATGATCCAGCGCGTAg                                                                                          <  2:287813/240‑1 (MQ=255)                                                                                                                                                     ggCCATCTGAGATGGCCGACAGGTTAACATAATTACAGCCATTTAGGGCCACGCGCTCAGCCAGCGGATGAAAACGCATCTATCACCTTGTGATCAACACC‑TTCATTCT‑CCCCCGCCAGTTCGGCACTGAGCTTCGCCATAAACTCGACTAAAAAATGCGCGTTGTGCTGCGCCAGCTGCTTGCCTCTTGCGGTTTGCATGGTTTGCGGTAGCTTCAGCAACTTAGTCTGGAAATGATCCAGCGCGTAGCg                                                                                        >  1:482227/1‑251 (MQ=255)                                                                                                                                                      gCCATCAGAGATGGCCGACAGGTTAACATAATTACAGCCATTTAGGGCCACGCGCTCAGCCAGCGGATGAAAACGCATCTATCACCTTGTGATCAACACC‑TTCATTCT‑CCCCCGCCAGTTCGGCACTGAGCTTCGCCATAAACTCGACTAAAAAATGCGCGTTGTGCTGCGCCAGCTGCTTGCCTCTTGCGGTTTGCATGGTTTGCGGTAGCTTCAGCAACTTAGTCTGGAAATGATCCAGCGCGTAGCGt                                                                                       >  2:519506/1‑251 (MQ=255)                                                                                                                                                             gagaTGGCCGACAGGTTAACATAATTACAGCCATTTAGGGCCACGCGCTCAGCCAGCGGATGAAAACGCATCTATCACCTTGTGATCAACACC‑TTCATTCT‑CCCCCGCCAGTTCGGCACTGAGCTTCGCCATAAACTCGACTAACAAATGCGCGTTGTGCTGCGCCAGCTGCTTGCCTCTTGCGGTTTGCATGGTTTGCGGTAGCTTCAGCAACTTAGTCTGGAAATGATCCAGCGCGTAGCGTTTatcat                                                                                >  1:356005/1‑251 (MQ=255)                                                                                                                                                                    ccGACAGGTTAACATAATTACAGCCATTTAGGGCCACGCGCTCAGCCAGCGGATGAAAACGCATCTATCACCTTGTGATCAACACC‑TTCATTCT‑CCCCCGCCAGTTCGGCACTGAGCTTCGCCATAAACTCGACTAAAAAATGCGCGTTGTGCTGCGCCAGCTGCTTGCCTCTTGCGGTTTGCATGGTTTGCGGTAGCTTCAGCAACTTAGTCTGGAAATGATCCAGCGCGTAGCGTTTATCATCAAGCgg                                                                         >  1:78642/1‑251 (MQ=255)                                                                                                                                                                      gACAGGTTAACATAATTACAGCCATTGAGGGCCACGCGCCCAGCCAGCGGATGAAAACGCAGCTATCACCTTTGGACCAACACC‑TTCATTCC‑CCCCCGCCAGTTCGGCACTGAGCTTCGTCATAAACTCGACTAAAAAATGCGCGTTGTGCTGCGCCAGCTGCTTGCCTCTTGCGGTTTGCATGGTTTGCGGTAGCTTCAGCAACTTAGTCTGGAAATGATCCAGCGCGTAGCGTTTATCATCAAGCGGGc                                                                       <  1:704505/251‑1 (MQ=255)                                                                                                                                                                          ggTTAACATAATTACAGCCATTTAGGGCCACGCGCTCAGCCAGCGGATGAAAACGCATCTATCACCTTGTGATCAACACC‑TTCATTCT‑CCCCCGCCAGTTCGGCACTGAGCTTCGCCATAAACTCGACTAAAAAATGCGCGTTGTGCTGCGCCAGCTGCTTGCCTCTTGCGGTTTGCATGGTTTGCGGTAGCTTCAGCAACTTAGTCTGGAAATGATCCAGCGCGTAGCGTTTATCATCAAGCGGGCGATg                                                                   <  1:102547/251‑1 (MQ=255)                                                                                                                                                                               aCATAATTACAGCCATTTAGGGCCACGCGCTCAGCCAGCGGATGAAAACGCATCTATCACCTTGTGATCAACACC‑TTCATTCT‑CCCCCGCCAGTTCGGCACTGAGCTTCGCCATAAACTCGACTAAAAAATGCGCGTTGTGCTGCGCCAGCTGCTTGCCTCTTGCGGTTTGCATGGTTTGCGGTAGCTTCAGCAACTTAGTCTGGAAATGATCCAGCGCGTAGCGTTTATCATCAAGCGGGCGATGCtgtg                                                              >  2:274360/1‑251 (MQ=255)                                                                                                                                                                                    aTTACAGCCTTTTAGGGCCACGCGCTCAGCCAGCGGATGAAAACGCATCTATCACCTTGGGATCAACACC‑TTTATGCC‑CCCCCGCCAGTTCGGCACTGAGCTTCGCCATAAACTCTACTAAAAAATGCGCGTTGTGCTGCGCCAGCTGCTTGCCTCTTGCGGTTTGCATGGTTTGCGGTAGCTTCAGCAACTTAGTCTGGAAATGATCCAGCGCGTAGCGTTTATCATCAAGCGGGCGATGCTTTGCGAAc                                                         <  1:463151/251‑1 (MQ=255)                                                                                                                                                                                    aTTACAGCCCTTTAGGGCCACGCGCTCAGCCAGCGGATGACAACGCATCTATCACCTTGTGATCAACACC‑TTCATTCT‑CCCCCGCCAGTTCGGCACTGAGCTTCGCCATAAACTCGACTAAAACATGCGAGTTGTGCTGCGCCAGCTGCTTGCCTCTTGCGGTTTGCATGGTTTGCGGTAGCTTCAGCAACTTAGTCTGGAAATGATCCAGCGCGTAGCGTTTATCATCAAGCGGGCGAAGCTGTGCGAAc                                                         >  2:352056/1‑251 (MQ=255)                                                                                                                                                                                    aTTACAGCCATTTCGAGCCACGCGCTCAGCCAGCGGATGCAAACGCATCTATCACCTTGTGATCAACACC‑TTCATTCT‑CCCCCGCCAGTTCGGCACTGAGCTTCGCCATAAACTCGACTAAAAAATGCGCGTTGTGCTGCGCCAGCTGCTTGCCTCTTGCGGTTTGCATGGTTTGCGGTAGCTTCAGCAACTTAGTCTGGAAATGATCCAGCGCGTAGCGTTTATCATCAAGCGGGCGATGCTGTGCGAAc                                                         >  2:365591/1‑251 (MQ=255)                                                                                                                                                                                             aTTTAGGACCACGCGCTCAGCCAGCGGCCGAAAACGCATCTATCACCTTGTGATCAACACC‑TTCATTCT‑CCCCCGCCAGTTCGGCACTGAGCTTCGCCATAAACTCGACTAAAAAATGCGCGTTGTGCTGCGCCAGCTGCTTGCCTCTTGCGGTTTGCATGGTTTGCGGAAGCTTCAGCAACTTAGTCTGGAAATTATCCAGCGCGTAGCGTTTATCATCCAGCGGGCGATGCTGTGCGCACGGATCTTc                                                 >  2:700302/1‑250 (MQ=255)                                                                                                                                                                                                 aGGGCCACGCGCTCAGCCAGCGGATGAAAACGCATCTATCACCTTGTGATCCACACC‑TTCATTCT‑CCCCCGCCAGTTCGGCACTGAGCTTCGCCATAAACTCGACTAAAAAATGCGCGTTGTGCTGCGCCAGCTGCTTGCCTCTTGCGGTTTGCATGGTTTGCGGTAGCTTCAGCAACTTAGTCTGGAAATg                                                                                                       >  2:453353/1‑192 (MQ=255)                                                                                                                                                                                                 aGGGCCACGCGCTCAGCCAGCGGATGAAAACGCATCTATCACCTTGTGATCAACCCC‑TTCATTCT‑CCCCCGCCAGTTCGGCACTGAGCTTCGCCATAAACTCGACTAAAAAATGCGCGTTGTGCTGCGCCAGCTGCTTGCCTCTTGCGGTTTGCATGGTTTGCGGTAGCTTCAGCAACTTAGTCTGGAAATg                                                                                                       <  1:453353/192‑1 (MQ=255)                                                                                                                                                                                                 aGGGCCACGCGCGCAGCCAGCGGATGAAACCGCATCTATCACCTGGTGAACAACACC‑TTCATTCG‑CCCCCGCCAGTTCGGCACTGAGCTTCGCCATAAACTCGACTAAAAAATGCGCGTTGTGCTGCGCCAGCTGCTTGCCTCTTGCGGTTTGCATGGTTTGCGGTAGCTTCAGCAACTTAGTCTGGAAATGATCCAGCGCGTAGCGTTTATCATCAAGCGGGCGATGCTGTGCGAACGGATCTTCTCCa                                             <  2:616325/250‑1 (MQ=255)                                                                                                                                                                                                    gCCACGCGCTCAGCCAGCGGATGAAAACGCATCTATCACCTTGTGATCAACACC‑TTCATTCT‑CCCCCGCCAGTTCGGCACTGAGCTTCGCCATAAACTCGACTAAAAAATGCGCGTTGTGCTGCGCCAGCTGCTTGCCTCTTGCGGTTTGCATGGTTTGCGGTAGCTTCAGCAACTTAGTCTGGAAATGATCCAGCGCGTAGCGTTTATCATCAAGCGGGCGATGCTGTGCGAACGGATCTTCACCATCa                                          >  2:180356/1‑250 (MQ=255)                                                                                                                                                                                                      cACGCGCTCAGCCAGCGGATGAAAACGCATCTATCACCTTGTGATCAACACC‑TTCATTCT‑CCCCCGCCAGTTCGGCACTGAGCTTCGCCATAAACTCGACTAAAAAATGCGCGTTGTGCTGCGCCAGCTGCTTGCCTCTTGCGGTTTGCATGGTTTGCGGTAGCTTCAGCAACTTAGTCTGGAAATGATCCAGCGCGTAGCGTTTATCATCAAGCGGGCGATGCTGTGCGGACGGATCTTCACCATCAAAc                                       >  2:663193/1‑251 (MQ=255)                                                                                                                                                                                                          cgcTCAGCCAGCGGATGAAAACGCATCTATCCCCTTGTGATCAACACC‑TTCATTCT‑CCCCCGCCAGTTCGGCACTGAGCTTCGCCATAAACTCGACTAAAAAATGCGCGTTGTGCTGCGCCAGCTGCTTGCCTCTTGCGGTTTGCATGGTTTGCGGTAGCTTCAGCAACTTAGTCTGGAAATGATCCAGCGCGTAGCGTTTATCATCAAGCGGGCGATGCTGTGCGAACGGATCTTCACCATCAAACAgcg                                   <  1:597268/251‑1 (MQ=255)                                                                                                                                                                                                              cagccagcGGATGAAAACGCATCTATCACCTTGTGATCAACACC‑TTCATTCT‑CCCCCGCCAGTTCGGCACTGAGCTTCGCCATAAACTCGACTAAAAAATGcgcg                                                                                                                                                                                 >  1:463352/1‑105 (MQ=255)                                                                                                                                                                                                               agccagcGGATGAAAACTCATCTATCGACT‑GTGGACAACCCCCTTCATGCC‑CCCCCGCCAGTTCGGCACTGAGCTTCGCCATAAACTCGACTAAAAAATGCGCGTTGTGCTGCGCCAGCTGCTTGCCTCTTGCGGTTTGCATGGTTTGCGGTAGCTTCAGCAACTTAGTCTGGAAATGATCCAGCGCGTAGCGTTTATCATCAAGCGGGCGATGCTGTGCGAACGGATCTTCACCATCAAACAGCGCCACg                              <  2:180384/251‑1 (MQ=255)                                                                                                                                                                                                                          gAAAACGCATCTATCACCTTGTGATCAACACC‑TTCATTCT‑CCCCCGCCAGTTCGGCACTGAGCTTCGCCATAAACTCGACTAAAAAATGCGCGTTGTGCTGCGCCAGCTGCTTGCCTCTTGCGGTTTGCATGGTTTGCGGTAGCTTCAGCAACTTAGTCTGGAAATGATCCAGCGCGTAGCGTTTATCATCAAGCGGGCGATGCTGTGCGAACGGATCTTCACCATCAAACAGCGCCACGCCCAATGcccc                   >  1:350620/1‑251 (MQ=255)                                                                                                                                                                                                                          gAAAACGCATCTATCACCTTGTGATCAACACC‑TTCATTCT‑ACCCCGCCAGTTCGGCACTGAGCTTCGCCATACACTCGACTAAAAACTGCGCGTTGTGCTGCGCCCGCTGCTTGCCTCTTGCGGTTTGCATGGTTTGCGGTAGCTTCAGCAACTTAGTCTGGAAATGATCCAGCGCGTAGCGTTTATCATCAAGCGGGCGATGCTGTGCGCACGGATCTTCACCATCAAACAGCGCCACGCCCAATGcccc                   >  2:132561/1‑251 (MQ=255)                                                                                                                                                                                                                           aaaaCGCATCTATCACCTTGTGATCAACACC‑TTCATTCT‑CCCCCGCCAGTTCGGCACTGAGCTTCGCCATAAACTCGACTAAAAAATGCGCGTTGTGCTGCGCCAGCTGCTTGCCTCTTGCGGTTTGCATGGTTTGCGGTAGCTTCAGCAACTTAGTCTGGAAATGATCCAGCGCGTAGCGTTTATCATCAAGCGGGCGATGCTGTGCGAACGGATCTTCACCATCAAACAGCGCCACGCCCACTGCCCCt                  >  1:332435/1‑251 (MQ=255)                                                                                                                                                                                                                            aaaCGCATCTATCACCTTGTGATCAACACC‑TTCATTCT‑CCCCCGCCAGTTCGGCACTGAGCTTCGCCATAAACTCGACTAAAAAATGCGCGTTGTGCTGCGCCAGCTGCTTGCCTCTTGCGGTTTGCATGGTTTGCGGTAGCTTCAGCAACTTAGTCTGGAAATGATCCAGCGCGTAGCGTTTATCATCAAGCGGGCGATGCTGTGCGAACGGATCTTCACCATCAAACAGCGCCACGCCCAATGCCCCTg                 >  1:556192/1‑251 (MQ=255)                                                                                                                                                                                                                                     tATCACCTTGTGATCAACACC‑TTCATTCT‑CCCCCGCCAGTTCGGCACTGAGCTTCGCCATAAACTCGACTAAAAAATGCGCGTTGTGCTGCGCCAGCTGCTTGCCTCTTGCGGTTTGCATGGTTTGCGGTAGCTTCAGCAACTTAGTCTGGAAATGATCCAGCGCGTAGCGTTTATCATCAAGCGGGCGATGCTGTGCGAACGGATCTTCACCATCAAACAGCGCCACGCCCAATGACCCTGAAACGGCa         >  2:231464/1‑250 (MQ=255)                                                                                                                                                                                                                                           cTTGTGATCAACACC‑TTCATTCT‑CCCCCGCCAGTTCGGCACTGAGCTTCGCCATAAACTCGACTAAAAAATGCGCGTTGTGCTGCGCCAGCTGCTTGCCTCTTGCGGTTTGCATGGTTTGCGGTAGCTTCAGCAACTTAGTCTGGAAATGATCCAGCGCGTAGCGTTTATCATCAAGCGGGCGATGCTGTGCGAACGGATCTTCACCATCAAACAGCGCCACGCCCAATGCCCCTGAACCGGCAAACAcgc  >  2:295430/1‑251 (MQ=255)                                                                                                                                                                                                                                                    |                                                                                                                                                                                                                                                    GCTCCAGTGATTCTCACCTACACCATTAAGCTGCGCCACCAACGTTTTAACATCATAAATTTGCAGCAGTTTCTCCAGCATTGCCTGTGCGGCATTCCCGGTGCTGTATGAATCGGTTACTGATATATTTTCCTGCATAGATTTCACCGGCCATCTGAGATGGCCGACAGGTTAACATAATTACAGCCATTTAGGGCCACGCGCTCAGCCAGCGGATGAAAACGCATCTATCACCTTGTGATCGACACC‑TTCATTCT‑CCCCCGCCAGTTCGGCACTGAGCTTCGCCATAAACTCGACTAAAAAATGCGCGTTGTGCTGCGCCAGCTGCTTGCCTCTTGCGGTTTGCATGGTTTGCGGTAGCTTCAGCAACTTAGTCTGGAAATGATCCAGCGCGTAGCGTTTATCATCAAGCGGGCGATGCTGTGCGAACGGATCTTCACCATCAAACAGCGCCACGCCCAATGCCCCTGAAACGGCAAACACGC  >  NC_000913/2032180‑2032664

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