Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	NC_000913	700,038	G→A	H97Y (CAT→TAT)	umpH ←	UMP phosphatase

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	700,038	0	G	A	97.4%	249.7 / ‑2.2	77	H97Y (CAT→TAT)	umpH	UMP phosphatase
Reads supporting (aligned to +/- strand):  ref base G (2/0);  new base A (42/33);  total (44/33)
Fisher's exact test for biased strand distribution p-value = 5.04e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 6.71e-01

TACCAACATAGAACGGTTTGCGCCCGGAGATTTTCTCAATCCCTGCACACAACGCGCCACAAGCGGGATAAAAACCGCGCCCGTGGGTGTCCGGATTGGTGGCGATAAAACGTGCACCGTTAGCGACGAAATAGGCTGCTTTATGCATCATGTCCCAGTTGTAGGAACGCGTTTCGCCAACAATCACGAAATCAGGGTTCACATCGGTAAT‑AGTGAAACCGGCTTTGTACAGTTCATGAATCAGTGCGCCTTCGCCCACCACATACGCTTTCTTGCCTTCCTGGCGACGCAGGAAATCGGCAGTCGCCATCGCAGAGGTATAAAACACGCTGTCAGGTACATCGACACCTGCGGTGGCAAAGCGGTTCGCCAGATCTTGCCCAGTCTGCGAAGGATAGTTGGTCAGCAACACCAGCGGCAGGCCTTTATCCATAATCCCGTGCAAAAATTCCGCTGCACCCGGTACGGCGACGTTATCGTG  >  NC_000913/699801‑700281                                                                                                                                                                                                                                               |                                                                                                                                                                                                                                                    ttacAACATAGAACGGTTTGCGCCCGGAGATTTTCTCAACCCCTGCACACAACGCGCCACAAGCGGGATAAAAACCGCGCCCGTGGGTGTCCGGATTGGTGGCGATAAAACGTTCACCGTTAGCGACGAAATAGGCTGCTTTATGCATCATGTCCCAGTTGTAGGAACGCGTTTCGCCAACAATCACGAAATCAGGGTTCACATCGGTAAT‑AGTGAAACCGGCTTTGTACAGTTCATAAATCAGTGCGcc                                                                                                                                                                                                                                         <  2:321141/247‑1 (MQ=255)  aCCAACATAGAACGGTTTGCGCCCGGAGATTTTCTCAATCCCTGCACACAACGCGCCACAAGCGGGATAAAAACCGCGCCCGTGGGCTTCCGGATTGGTGGCGATAAAACGTGCACCGTTAGCGACGAAATAGGATGCTTTATCCATCATGTCCCCGATGTAAGAACGCGTTTCGCCAACAATCACGAAATCAGGGTTCACATCGGCAAT‑AGTGAAACCGGCTTTGTACAGTATATAAATCAGTGCTCCtt                                                                                                                                                                                                                                       >  2:562505/1‑251 (MQ=255)                     cgcCCGGAGATTTTCTCAATCCCTGCACACAACGCGCCACAAGCGGGATAAAAACCGCGCCCGTGGGTGTCCGGATTGGTGGCGATAAAACATGCACCGTTAGCGACGAAATAGGCTGCTTTATGCATCATGTCCCAGTTGTAGGAACGCGTTTCGCCAACAATCACGAAATCAGGGTTCACATCGGTAAT‑AGTGAAACCGGCTTTGTACAGTTCATAAATCAGTGCGCCTTCGCCCACCACATACGCttt                                                                                                                                                                                                                    >  1:27725/1‑251 (MQ=255)                      gcCCGGAGATTTTCTCAATCCCTGCACACAACGCGCCACAAGCGGGATAAAAACCGCGCCCGTGGGTGTCCGGATTGGTGGCGATAAAACGTGCACCGTTAGCGACGAAATAGGCTGCTTTATGCATCATGTCCCAGTTGTAGGAACGCGTTTCGCCAACAATCACGAAATCAGGGTTCACATCGGCAAT‑AGTGAAACCGGCTTTGTACAGTTCATAAATCAGTGCGCCTTCGCCCACCACATACGCTTTc                                                                                                                                                                                                                   >  1:381336/1‑251 (MQ=255)                        ccGGAGATTTTCTCAATCCCTGCACACAACGCGCCACAAGCGGGATAAAAACCGCGCCCGTGGGTGTCCGGATTGGTGGCGATAAAACGTGCACCGTTAGCGACGAAATAGGCTGCTTTATGCATCATGTCCCAGTTGTAGGAACGCGTTTCGCCAACAATCACGAAATCAGGGTTCACATCGGTAAT‑AGTGAAACCGGCTTTGTACAGTTCATAAATCAGTGCGCCTTCGCCCCCCACATACGCTTTCtt                                                                                                                                                                                                                 >  2:7452/1‑251 (MQ=255)                               ttttCTCAATCCCTGCACACAACGCGCCACAAGCGGGATAAAAACCGCGCCCGTGGGGGTCCGGATTGGTGGCGATAAAACGTGCCCCGTTAGCGACGAAATAGGCTGCTTTATGCATCATGTCCCAGTTGTAGGAACGCGTTTCGCCAACAATCACGCAATCAGGGTTCACATCGGTAAT‑AGTGAAACCGGCTTTGTACAGTTCATAAATCAGAGCGCCTTCGCCCACCACATACGCTTTCTTGCCTTcc                                                                                                                                                                                                          >  2:549589/1‑251 (MQ=255)                                      aaTCCCTGCACACAACGCGCCACAAGCGGGATAAAAACCGCGCCCGTGGGTGTCTGGATTGGTGGCGATAAAACGTGCACCGTTAGCGACGAAATAGGCTGCTTTATGCATCATGTCCCAGTTGTAGGAACGCGTTTCGCCAACAATCACGAAATCAGGGTTCACATCGGTAAT‑AGTGAAACCGGCTTTGTACAGTTCATAAATCAGTGCGCCTTCGCCCACCACATACGCTTTCTTGCCTTCCTGGCGAc                                                                                                                                                                                                   >  1:569907/1‑251 (MQ=255)                                       aTCCCTGCACACAACGCGCCACAAGCGGGATAAAAACCGCGCCCGTGGGTGTCCGGATTGGTGGCGATAAAACGTGCACCGTTAGCGACGAAATAGGCTGCTTTATGCATCATGTCCCAGTTGTAGGAACGCGTTTCGCCAACAATCACGAAATCAGGGTTCACATCGGTAAT‑AGTGAAACCGGCTTTGTACAGTTCATAAATCAGTGCGCCTTCGCCCACCACATACGCTTTCTTGCCTTCCTGGCgccg                                                                                                                                                                                                  >  2:92675/1‑248 (MQ=255)                                            tGCACACAACGCGCCACAAGCGGGATAAAAACCGCGCCCGTGGGTGTCCGGATTGGTGGCGATAAAACGTGCACCGTTAGCGACGAAATAGGCTGCTTTATGCATCATGTCCCAGTTGTAGGAACGCGTTTCGCCAACAATCACGAAATCAGGGTTCACATCGGTAAT‑AGTGAAACCGGCTTTGTACAGTTCATAAATCAGTGCGCCTTCGCCCACCACATACGCTTTCTTGCCTTCCTGGCGACGCAgg                                                                                                                                                                                              >  1:76491/1‑250 (MQ=255)                                             gCACACAACGCGCCACAAGCGGGATAAAAACCGCGCCCGTGGGTGTCCGGATTGGTGGCGATAAAACGTGCACCGTTAGCGACGAAATAGGCTGCTTTATGCATCATGTCCCAGTTGTAGGAACGCGTTTCGCCAACAATCACGAAATCAGGGTTCACATCGGTAAT‑AGTGAAACCGGCTTTGTACAGTTCATAAATCAGTGCGCCTTCGCCCACCACATACGCTTTCTTGCCTTCCTGGCGACGCAgg                                                                                                                                                                                              <  1:36701/249‑1 (MQ=255)                                              cacacaACGCGCCACAAGCGGGATAAAAACCGCGCCCGTGGGTGTCCGGATTGGTGGCGATAAAACGAGCACCGTTAGCGACGAAATAGGCTGCTTTATGCATCATGTCCCAGTTGTAGGAACGCGTTTCGCCAACAATCACGAAATCAGGGTTCACATCGGAAAT‑AGTGAAACCGGCTTTGTACAGTTCATAAATCAGTGCGCCTTAGACCACCACATACGCTTTCTTGCCTTCCTGGCGACGCAGGaa                                                                                                                                                                                            >  2:578867/1‑250 (MQ=255)                                                   aaCGCGCCACAAGCGGGATAAAAACCGCGCCCGTGGGTGTCCGGATTGGTGGCGATAAAACGTGCACCGTTAGCGACGAAATAGGCTGCTTTATGCATCATGTCCCAGTTGTAGGAACGCGTTTCGCCAACAATCACGAAATCAGGGTTCACATCGGTAAT‑AGTGAAACCGGCTTTGTACAGTTCATAAATCAGTGCGCCTTCGCCCACCACATACGCTTTCTTGCCTTCCTGGCGACGCAGGAAATCGGc                                                                                                                                                                                      >  1:509486/1‑251 (MQ=255)                                                     cgcgCCACAATCGGGATAAAAACCGCGCCCGTGGGTGTCCGGATTGGTGGCGATAAAACGTGCACCGTTAGCGACGAAATAGGCTGCTTTATGCATCATGTCCCAGTTGTAGGAACGCGTTTCGCCAACAATCACGAAATCAGGGTTCACAACGGTAAT‑AGTGAAACCGGCTTTGTACAGTTCATAAATCAGTGCGCCTTCGCCCACCACATACGCTTTCTTGCCTTCCTGGCGACGCAGGAAATCGGCAg                                                                                                                                                                                    <  2:438407/251‑1 (MQ=255)                                                     cgcgCCACAAGCGGGATAAAAACCGCGCCCGTGGGTGTCCGGATTGGTGGCGATAAAACGTGCACCGTTAGCGACGAAATAGGCTGCTTTATGCATCATGTCCCAGTTGTAGGAACGCGTTTCGCCAACAATCACGAAATCAGGGTTCACATCGGTAAT‑AGTGAAACCGGCTTTGTACAGTTCATAAATCAGTGCGCCTTCGCCCACCACATACGCTTTCTTGCCTTCCTGGCGACGCAGGAAATCGGCAg                                                                                                                                                                                    <  2:261650/251‑1 (MQ=255)                                                       cgcCACAAGCGGGAAAAAAACCGCGCCCGTGGGTGTCCGGATTGGTGGCGATAAAACGTGCACCGTTAGCGACGAAATAGGCTGCTTTATGCAGCATGTCCCAGTTGTAGGAACGCGTTTCGCCAACAATCACGAAATCAGGGTTCACATCGGCAAT‑AGTGAAACCGGCTTTGTACAGTTCATAAATCAGTGCGCCTTCGCCCACCACATACGCTTTCTTGCCTTCCTGGCGACGCAGGAAATCGGCAGTc                                                                                                                                                                                  <  2:381336/251‑1 (MQ=255)                                                        gcCACAAGCGGGATAAAAACCGCGCCCGTGGGTGTCCGGATTGGTGGCGATAAAACGTGCACCGTTAGCGACGAAATAGGCTGCTTTATGCATCATGTCCCAGTTGTAGGAACGCGTTTCGCCAACAATCACGAAATCAGGGTTCACATCGGTAAT‑AGTGAAACCGGCTTTGTACAGTTCATAAATCAGTGCGCCTTCGCCCCCCACATACGCTTTCTGTCCTTCCTGGCGACGCAGGAAAGCGGCAGTcc                                                                                                                                                                                 >  2:595712/1‑250 (MQ=255)                                                             aaGCGGGATAAAAACCGCGCCCGTGGGTGTCCGGATTGGTGGCGATAAAACGTGCACCGTTAGCGACGAAATAGGCTGCTTTATGCATCATGTCCCAGTTGTAGGAACGCGTTTCGCCAACAATCACGAAATCAGGGTTCACATCGGTAAT‑AGTGAAACCGGCTTTGAACAGTTCATAAATCAGTGCGCCTTCGACCACCACATACGCTTTCTTGCCTTCTTGGCGACGCAGGCAATCGGCAgcacccagc                                                                                                                                                                            >  2:657253/1‑243 (MQ=255)                                                              aGCGGGATAAAAACCGCGCCCGTGGGTGTCCGGATTGGTGGCGATAAAACGTGCACCGTTAGCGACGAAATAGGCTGCTTTATGCATCATGTCCCAGTTGTAGGAACGCGTTTCGCCAACAATCACGAAATCAGGGTTCACATCGGTAAT‑AGTGAAACCGGCTTTGTACAGTTCATAAATCAGTGCGCCTTCGCCCACCACATACGCTTTCTTGCCTTCCTGGCGACGCAGGAAATCGGCAGTCGCCATCg                                                                                                                                                                           >  1:193735/1‑251 (MQ=255)                                                                      aaaaaCCGCGCCCGTGGGTGTCCGGATTGGTGGCGATAAAACGTGCACCGTTAGCGACGAAATAGGCTGCTTTATGCATCATGTCCCAGTTGTAGGAACGCGTTTCGCCAACAATCACGAAATCAGGGTTCACATCGGTAAT‑AGTGAAACCGGCTTTGTACAGTTCATAAATCAGTGCGCCTTCGCCCACCACATACGCTTTCTTGCCTTCCTGGCGACGCAGGAAATCGGCAGTCGCCATCGCAGAGGta                                                                                                                                                                   <  1:250606/251‑1 (MQ=255)                                                                       aaaaCCGCGCCCGTGGGTGTCCGGATTGGTGGCGATAAAACGTGCACCGTTAGCGACGAAATAGGCTGCTTTATGCATCATGTCCCAGTTGTAGGAACGCGTTTCGCCAACAATCACGAAATCAGGGTTCACATCGGTAAT‑AGTGAAACCGGCTTTGTACAGTTCATAAATCAGTGCGCCTTCGCCCACCACATACGCTTTCTTGCCTTCCTGGCGACGCAGGAAATCGGCAGTCGCCATCGCagag                                                                                                                                                                      >  1:454149/1‑247 (MQ=255)                                                                       aaaaCCGCGCCCGTGGGTGTCCGGATTGGTGGCGATAAAACGTGCACCGTTAGCGACGAAATAGGCTGCTTTATGCATCATGTCCCAGTTGTAGGAACGCGTTTCGCCAACAATCACGAAATCAGGGTTCACATCGGTAAT‑AGTGAAACCGGCTTTGTACAGTTCATAAATCAGTGCGCCTTCGCCCACCACATACGCTTTCTTGCCTTCCTGGCGACGCAGGAAATCGGCAGTCGCCATCGCAGAGGtat                                                                                                                                                                  >  1:522994/1‑251 (MQ=255)                                                                       aaaaCCGCGCCCGTGGGTGTCCGGATTGGTGGCGATAAAACGTGCACCGTTAGCGACGAAATAGGCTGCTTTATGCATCATGTCCCAGTTGTAGGAACGCGTTTCGCCAACAATCACGAAATCAGGGTTCACATCGGTAAT‑AGTGAAACCGGCTTTGTACAGTTCATAAATCAGTGCGCCTTCGCCCACCACATACGCTTTCTTGCCTTCCTGGCGACGCAGGAAATCGGCAGGCGCCATCGCagag                                                                                                                                                                      <  2:454149/247‑1 (MQ=255)                                                                           ccGCGCCCGTGGGTGGCCGGATTGGTGGCGATAAAACGTGCACCGTTAGCGACGAAATAGGCTGCTTTATGCAGCATGTCCCAGTTGTAGGAACGCGTTTCGCCAACAATCACGAAATCAGGGTTCACATCGGTAAT‑AGTGAAACCGGCTTTGTACAGTTCATAAATCAGTGCGCCTTCGCCCACCACATACGCTTTCTTGCCTTCCTGGCGACGCAGGAAATCGGCAGTCGCCATCGCAGAGGTATaaaa                                                                                                                                                              <  2:165054/251‑1 (MQ=255)                                                                            cgcgCCCGTGGGTGTCCGGATTGGTGGCGATAAAACGTGCACCGTTAGCGACGAAATAGGCTGCTTTATGCATCATGTCCCAGTTGTAGGAACGCGTTTCGCCAACAATCACGAAATCAGGGTTCACATCGGTAAT‑AGTGAAACCGGCTTTGTACAGTTCATAAATCAGTGCGCCTTCGCCCACCACATACGCTTTCTTGCCTTCCTGGCGACGCAGGAAATCGGCAGTCGCCATCGCAGAGGTATAAAac                                                                                                                                                             >  1:350648/1‑251 (MQ=255)                                                                               gcCCGTGGGAGTCCGGATTGGTGGCGCTCAAACGTGCACCGTTAGCGACGGACTAGGCTGCTTTATGCATCATGTCCCAGTTGTAGGAACGCGTTTCGCCAACAATCACGAAATCAGGGTTCACATCGGTAAT‑AGTGAAACCGGCTTTGTACAGTTCATAAATCAGTGCGCCTTCGCCCACCACATACGCTTTCTTGCCTTCCAGGCGACGCAGGCAATCGCCAGCCGCCATCGCAAAGGTATAAAACACg                                                                                                                                                          >  2:585003/1‑251 (MQ=255)                                                                                     gggTGTCCGGATTGGTGGCGATAAAACGTGCACCGTTAGCGACGAAATAGGCTGCTTTATGCATCATGTCCCAGTTGTAGGAACGCGTTTCGCCAACAATCACGAAATCAGGGTTCACATCGGTAAT‑AGTGAAACCGGCTTTGTACAGTTCATAAATCAGTGCGCCTTCGCCCACCACATACGCTTTCTTGCCTTCCTGGCGACGCAGGAAATCGGCAGTCGCCATCGCAGAGGTATAAAACACGCTGTc                                                                                                                                                     <  2:345162/250‑1 (MQ=255)                                                                                      ggTGTCCGGATTGGTGGCGATAAAACGTGCACCGTTAGCGACGAAATAGGCTGCTTTATGCATCATGTCCCAGTTGTAGGAACGCGTTTCGCCAACAATCACGAAATCAGGGTTCACATCGGTAAT‑AGTGAAACCGGCTTTGTACAGTTCATAAATCAGTGCGCCTTCGCCCACCACATACGCTTTCTTGCCTTCCTGGCGACGCAGGAAATCGGCAGTCGCCATCGCAGAGGTATAAAACACGCTGTcgg                                                                                                                                                   >  2:242566/1‑249 (MQ=255)                                                                                            cGGATTGGTGGCGATAAAACGTGCACCGTTAGCGACGAAATAGGCTGCTTTATGCATCATGTCCCAGTTGTAGGAACGCGTTTCGCCAACAATCACGAAATCAGGGTTCACATCGGTAAT‑AGTGAAACCGGCTTTGTACAGTTCATAAATCAGTGCGCCTTCGCCCACCACATACGCTTTCTTGCCTTCCTGGCGACGCAGGAAATCGGCAGTCGCCATCGCAGAGGTATAAAACACGCTGTCAGGTACAt                                                                                                                                             <  1:242309/251‑1 (MQ=255)                                                                                                         aTAAAACGTGCACCGTTAGCGACGAAATAGGCTGCTTTATGCATCATGTCCCAGTTGTAGGAACGCGTTTCGCCAACAATCACGAAATCAGGGTTCACATCGGTAAT‑AGTGAAACCGGCTTTGTACAGTTCATAAATCAGTGCGCCTTCGCCCACCACATACGCTTTCTTGCCTTCCTGGCGACGCAGGAAATCGGCAGTCGCCATCGCAGAGGTATAAAACACGCTGTCAGGTACATCGACACCTGCGGt                                                                                                                                <  1:551389/251‑1 (MQ=255)                                                                                                                gTGCACCGTTAGCGACGAAATAGGCTGCTTTATGCATCATGTCCCAGTTGTAGGAACGCGTTTCGCCAACAATCACGAAATCAGGGTTCACATCGGTAAT‑AGTGAAACCGGCTTTGTACAGTTCATAAATCAGTGCGCCTTCGCCCACCACATACGCTTTCTTGCCTTCCTGGCGACGCAGGAAATCGGCAGTCGCCATCGCAGAGGTATAAAACACGCTGTCAGGTACATCGACACCTGCGGTGGCaa                                                                                                                           >  1:74102/1‑249 (MQ=255)                                                                                                                gTGCACCGTTAGCGACGAAATAGGCTGCTTTATGCATCATGTCCCAGTTGTAGGAACGCGTTTCGCCAACAATCACGAAATCAGGGTTCACATCGGTAAT‑AGTGAAACCGGCTTTGTACAGTTCATAAATCAGTGCGCCTTCGCCCACCACATACGCTTTCTTGCCTTCCTGGCGACGCAGGAAATCGGCAGTCGCCATCGCAGAGGTATAAAACACGCTGTCAGGTACATCGACACCTGCGGTGGCaa                                                                                                                           >  1:37725/1‑249 (MQ=255)                                                                                                                     ccGTTAGCGACGAAATAGGCTGCTTTATGCATCATGTCCCAGTTGTAGGAACGGGTTTCGCCAACAATCACGAAATCAGGGTTCACATCGGTAAT‑AGTGAAACCGGCTTTGTACAGTTCATAAATCAGTGCGCCTTCGCCCACCACATACGCTTTCTTGCCTTCCTGGCGACGCAGGAAATCGGCAGTCGCCATCGCAGAGGTATAAAACACGCTGTCAGGTACATCGACACCTGCGGTGGCAAAGCGGtt                                                                                                                    >  1:470194/1‑251 (MQ=255)                                                                                                                      cGTTAGCGACGAAATAGGCTGCTTTATGCATCATGTCCCAGTTGTAGGAACGCGTTTCGCCAACAATCACGAAATCAGGGTTCACATCGGTAAT‑AGTGAAACCGGCTTTGTACAGTTCATAAATCAGTGCGCCTTCGCCCACCACATACGCTTTCTTGCCTTCCTGGCGACGCAGGAAATCGGCAGTCGCCATCGCAGAGGTATAAAACACGCTGTCAGGTACATCGACACCTGCGGTGGCAAAACGGTTc                                                                                                                   >  2:347460/1‑251 (MQ=255)                                                                                                                       gTTAGCGTCTAAATAGGCTGCTTTATGCATCATGTCCCAGTTGTAGGAACGCGTTTCGCCAACAATCACGAAATCAGGGTTCACATCGGTAAT‑AGTGAAACCGGCTTTGTACAGTTCATAAAGCAGTGCGCCTTCGCCCACCACATACGCTTTCTTGCCTTCCTGGCGACGCAGGAAATCGGCAGTCGCCATCGCAGAGGTATAAAACACGCTGTCAGGTACATCGACACCTGCGGTGGCAAAGCGGTTCg                                                                                                                  <  2:552331/251‑1 (MQ=255)                                                                                                                             gacgaAATAGGCTGCTTTATGCATCATGTCCCAGTTGTAGGAACGCGTTTCGCCAACAATCACGAAATCAGGGTTCACATCGGTAAT‑AGTGAAACCGGCTTTGTACAGTTCATAAATCAGTGCGCCTTCGCCCACCACATACGCTTTCTTGCCTTCCTGGCGACGCAGGAAATCGGCAGTCGCCATCGCAGAGGTATAAAACACGCTGTCAGGTACATCTACCTCTGCGGTGGCAAAGCGGTTCGcc                                                                                                                >  2:473508/1‑247 (MQ=255)                                                                                                                             gacgaAATAGGCTGCTTTATGCATCATGTCCCAGTTGTAGGAACGCGTTTCGCCAACAATCACGAAATCAGGGTTCACATCGGTAAT‑AGTGAAACCGGCTTTGTACAGTTCATAAATCAGTGCGCCTTCGCCCACCACATACGCTTTCTTGCCTTCCTGGCGACGCAGGAAATCGGCAGTCGCCATCGCAGAGGTATAAAACACGCTGTCAGGTACATCGACACCTGCGGTGGCAAAGCGGTTCGCCAGAt                                                                                                            <  1:116840/251‑1 (MQ=255)                                                                                                                             gacgaAATAGGCTGCTTTATGCAGCATGTCCCAGTTGTAGGAACGCGGTTCGCCAACAATCACGAAAGCAGGGTTCACAGCGGTAAT‑AGTGAAACCGGCTTTGTACAGTTCATAAATCAGTGCGCCTTCGCCCACCACATACGCTTTCTTGCCTTCCTGGCGACGCAGGAAATCGGCAGTCGCCATCGCAGAGGTATAAAACACGCTGTCAGGTACATCGACACCTGCGGTGGCAAAGCGGTTCGCCAGAt                                                                                                            <  2:88316/251‑1 (MQ=255)                                                                                                                                gaAATAGGCTGCTTTATGCATCATGTCCCAGTTGTAGGAACGCGTTTCGCCAACAATCACGAAATCAGGGTTCACATCGGTAAT‑AGTGAAACCGGCTTTGTACAGTTCATAAATCAGTGCGCCTTCGCCCACCACATACGCTTTCTTGCCTTCCTGGCGACGCAGGAAATCGGCAGTCGCCATCGCCGCGGTATAAAACACGCTGTCAGGTACATCGCCACCCGCGGGGGCAACGCGGTTAGCCAGATCtt                                                                                                         >  2:420465/1‑251 (MQ=255)                                                                                                                                        ctgctTTATGCATCATGTCCCAGTTGTAGGAACGCGTTTCGCCAACAATCACGAAATCAGGGTTCACATCGGTAAT‑AGTGAAACCGGCTTTGTACAGTTCATAAATCAGTGCGCCTTCGCCCACCACATACGCTTTCTTGCCTTCCTGGCGACGCAGGAAATCGGCAGTCGCCATCGCAGAGGTATAAAACACGCTGTCAGGTACCTCGACACCTGCGGTGGCAAAGCGGTTCGCCAGATCTTGCCCAGTc                                                                                                 >  2:512535/1‑251 (MQ=255)                                                                                                                                          gctTGATGCATCATGTCCCAGTTGTAGGAACGCGTTTCGCCAACAATCACGAAATCAGGGTTCACATCGGTAAT‑AGTGAAACCGGCTTTGTACAGTTCATAAATCAGTGCGCCTTCGCCCACCACATACGCTTTCTTGCCTTCCTGGCGACGCAGGAAATCGGCAGTCGCCATCGCAGAGGTATAAAACACGCTGTCAGGTACATCGACACCTGCGGTGGCAAAGCGGTTCGCCAGATCTTGCCCAGTCt                                                                                                <  2:21403/250‑1 (MQ=255)                                                                                                                                             ttttGCATCATGTCCCAGTTGGAGGAACGCGGTTCGCCAACAATCACGAAATCAGGGTTCACATCGGTAAT‑AGTGAAACCTGCGTTGTACAGTTCATAAAGTAGTGCGCCTTCGCCCACCACATACGCTGTCTTGCCTTCCTGGCGACGCAGGAAATCGTCAGTCGCCATCGCAGAGGTATAAAACTCGCTGTCAGGTTCATCGACACCTGCTGTGGCAAAGCGGTTCGCCAGATCTTGCCTAGTCTg                                                                                               <  2:684160/245‑1 (MQ=255)                                                                                                                                             ttATGCATCATGTCCCAGTTGTAGGAACGCGTTTCGCCAACAATCACGAAATCAGGGTTCACATCGGTAAT‑AGTGAAACCGGCTTTGTACAGTTCATAAATCAGTGCGCCTTCGCCCACCACATACGCTTTCTTGCCTTCCTGGCGACGCAGGAAATCGGCAGTCGCCATCGCAGAGGTATAAAACACGCTGTCAGGTACATCGACACCTGCGGTGGCAAAGCGGTTCGCCAGATCTTGCCCAGTCTGc                                                                                              <  1:542262/249‑1 (MQ=255)                                                                                                                                             ttATGCATCATGTCCCAGTTGTAGGAACGCGTTTCGCCAACAATCACGAAATCAGGGTTCACATCGGTAAT‑AGTGAAACCGGCTTTGTACAGTTCATAAATCAGTGCGCCTTCGCCCACCACATACGCTTTCTTGCCTTCCTGGCGACGCAGGAAATCGGCAGTCGCCATCGCAGAGGTATAAAACACGCTGTCAGGTACATCGACACCTGCGGTGGCAAAGCGGTTCGCCAGATCTTGCCCAGTCTGCg                                                                                             >  2:431804/1‑250 (MQ=255)                                                                                                                                             ttATGCATCATGTCCCAGTTGTAGGAACGCGTTTCGACAACAATCACGAAATCAGGGTTCACATCGGTAAT‑AGTGAAACCGGCTTTGTACAGTTCATAAATCAGTGCGCCTTCGCCCACCACATACGCTTTCTTGCCTTCCTGGCGACGCAGGAAATCGGCAGTCGCCATCGCAGAGGAATAAAACACGCTGTAAGGTACATCGAAACCTGCGGTGGCAAAGCGGTTCGCCAGATCTTGCCCAGTATGCgg                                                                                            >  1:444631/1‑250 (MQ=255)                                                                                                                                               aTGCATCATGTCCCAGTTGTAGGAACGCGTTTCGCCAACAATCACGAAATCAGGGTTCACATCGGTAAT‑AGTGAAACCGGCTTTGTACAGTTCATAAATCAGTGCGCCTTCGCCCACCACATACGCTTTCTTGCCTTCCTGGCGACGCAGGAAATCGGCAGTCGCCATCGCAGAGGTATAAAACACGCTGTCAGGTACATCGACACCTGCGGTGGCAAAGCGGTTCGCCAGATCTTGTCCAGTCTGCGAAg                                                                                          <  2:193735/251‑1 (MQ=255)                                                                                                                                                  catcatGTCACAGTTGAAGGTACGCGTTTCGCCAAAAATCACGAAATCAGGGTTCAATTCGGTAAT‑AGAGAAACCGGCTTTGTACAGTTCATAAATAAGTGCGCCTTCACCCACCACATACGCTTTCTTGCCTTCATGGCGACGCAGACAATCGGCAGTCGCCATCGCAGAGGTATAAAACAAGCTGTCAGGTACATCGACACCTGCGGCGGCAAAGCGGTTCGCCAGATCTTGCACAGGCAGCGAAGGa                                                                                        >  2:33925/1‑250 (MQ=255)                                                                                                                                                      atGTCCCAGTTGTAGGAACGCGTTTCGCCAACAATCACGAAATCAGGGTTCACATCGGTAAT‑AGTGAAACCGGCTTTGTACAGTTCATAAATCAGTGCGCCTTCGCCCACCACATACGCTTTCTTGCCTTCCTGGCGACGCAGGAAATCGGCAGTCGCCATCGCaga                                                                                                                                                                       <  2:643215/167‑1 (MQ=255)                                                                                                                                                      atGTCCCAGTTGTAGGAACGCGTTTCGCCAACAATCACGAAATCAGGGTTCACATCGGTAAT‑AGTGAAACCGGCTTTGTACAGTTCATAAATCAGTGCGCCTTCGCCCACCACATACGCTTTCTTGCCTTCCTGGCGACGCAGGAAATCGGCAGTCGCCATCGCaga                                                                                                                                                                       >  1:643215/1‑167 (MQ=255)                                                                                                                                                       tGTCCCAGTTGTAGGAACGCGTTTCGCCAACAATCACGAAATCAGGGTTCACATCGGTAAT‑AGTGACACCGGCTTTGTACAGTTCATAAATCAGTGCGCCTTCGCCCACCACATACGCTTTCTTGCCTTCCTGGCGACGCAGGAAATCGGCAGTCGCCATCGCAGAGGTATAAAACACGCTGTCAGGTACATAGACACCTGCGGTAGCAAAGCGGTTCGCCAGATCTTGCCCAGACTGCACAGGATAgctg                                                                                  >  2:573565/1‑248 (MQ=255)                                                                                                                                                         tCCCAGTTGTAGGAACGCGTTTCGCCAACAATCACGAAATCAGGGTTCACATCGGTAAT‑AGTGACACCGGCTTTGTACAGTTCATAAATCAGTGCGCCTTCGCCCACCACATACGCTTTCTTGCCTTCCTGGCGACGCAGGAAATCGGCAGTCGCCATCGCAGAGGTATAAAACACGCTGTCAGGTACATCGACACCTGCGGTGGCAAAGCGGTTCGCCAGATCTTGCCCAGTCTGCGAAGGATAGTTGGt                                                                                >  2:60653/1‑251 (MQ=255)                                                                                                                                                              gTTGTAGGAACGCGTTTCGCCAACAATCACGAAATCAGGGTTCACATCGGTAAT‑AGTGAAACCGGCTTTGTACAGTTCATAAATCAGTGCGCCTTCGCCCACCACATACGCTTTCTTGCCTTCCTGGCGACGCAGGAAATCGGCAGTCGCCATCGCAGAGGAATAAAACACGCTGTCAGGTACATCGACACCTGCGGTGGCAAGGCGGTTCGCCAGATCTTGCCCAGTCTGCGAAGGATAGTTGGTCAGc                                                                            >  2:137935/1‑250 (MQ=255)                                                                                                                                                                  tAGGAACGCGTTTCGCCAACAATCACGAAATCAGGGTTCACATCGGTAAT‑AGTGAAACCTGCTTTGTACAGTTCATAAATCAGTGCGCCTTCGCCCACCACATACGCTTTCTTGCCTTCCTGGCGACGCAGGAAATCGGCAGTCGCCATCGCAGAgg                                                                                                                                                                     <  1:607458/157‑1 (MQ=255)                                                                                                                                                                  tAGGAACGAGTTTCGCCAACAATAACGAAATCAAGGTTCACATCGGTAAT‑AGAGAAACCGGCTTTGTACAGTTCATAAATCAGTGCGCCTTCGCCCACCACATACGCTTTCTTGCCTTCCTGGCGACGCAGGAAATCGGCAGTCGCCATCGCAGAgg                                                                                                                                                                     >  2:607458/1‑157 (MQ=255)                                                                                                                                                                   aGGAACGCGTTTCGCCAACAATCACTAAATCAGGGTTCACATCGGTAAT‑AGTGAAACCGGCGTTGTACAGTTCATAAATCAGTGCGCCTTCGCCCACCACATACGCTTTCTTGCCTTCCTGGCGACGCAGGAAATCGGCAGGCGCCATCGCAGAGGTATAAAACACGCTGTCAGGTACATCGACACCGGCGGTGGCAAAGCGGTTCGCCAGATCTTGCCCAGTCTGCGAAGGATAGTTGGGCAGCAACAcc                                                                      <  2:657086/251‑1 (MQ=255)                                                                                                                                                                   aGGAACGCGTTTCGCCAACAATCACGAAATCAGGGTTCACATCGGTAAT‑AGTGAAACCGGCTTTGTACAGTTCATAAATCAGTGCGCCTTCGCCCACCACATACGCTTTCTTGCCTTCCTGGCGACGCAGGAAATCGGCAGTCGCCATCGCAGAGGTATAAAACACGCTGTCAGGTACATCGACACCTGCGGTGGCAAAGCGGTTCGCCAGATCTTGCCCAGTCTGCGAAGGATAGTTGGTCAGCAACAcc                                                                      <  1:512535/251‑1 (MQ=255)                                                                                                                                                                        cgcgTTTCGCCAACAATCACGAAATCAGGGTTCACATCGGTAAT‑AGTGAAACCGGATTTGTACAGTTCATAAATCAGTGCGCCTTCGCCCACCACATACGCTTTCTTGCCTTCCTGGCGACGCAGGAAATCGGCAGTCGCCATCGCaga                                                                                                                                                                       >  1:412638/1‑149 (MQ=255)                                                                                                                                                                        cgcgTTTCGCCAACAATCACGAAAGCAGGGTTCACATCGGTAAT‑AGTGAAACCGGCTTTGTACAGTTCATAAATCAGGGCGCCTGCGCCCACCTCATACGCTTTCTTGCCTTCCGGGCGACGCAGGAAATCGGCAGTCGCCATCGCaga                                                                                                                                                                       <  2:412638/149‑1 (MQ=255)                                                                                                                                                                          acgcatttccgatctCACGAAATCAGGCTCCACATCGGTAAT‑AGTGAAACCGGCTTTGTACAGTTCATAAAACAGTGCGCCTTCTCCCACCACATACGCTTTCTTGCCTTCCTGGCGACGCAGGAAATCGGCAGTCGCCATCGCAGAGGTATAAAACACGCTGTCAGGTACATCGACTCCTGCGGTGGCAAAGCGGTTCGCCAGATCTTGCCCAGTCTGCGAAGGATAGTTGGTCAGCAACACCAGCGGCa                                                               <  2:488958/237‑1 (MQ=255)                                                                                                                                                                           gTTTCGCCAACAATCACGAAATCAGGGTTCACATCGGTAAT‑AGTGAAACCGGCTTTGTACAGTTCATAAATCAGTGCGCCTTCGCCCACCACATACGCTTTCTTGCCTTCCTGGCGACGCAGGAAATCGGCAGTCGCCATCGCAGAGGTATaaaa                                                                                                                                                              >  1:13080/1‑155 (MQ=255)                                                                                                                                                                           gTTTCGCCAACAATCACGAAATCAGGGTTCACATCGGTAAT‑AGTGAAACCGGCTTTGTACAGTTCATAAATCAGTGCGCCTTCGCCCACCACATACGCTTTCTTGCCTTCCTGGCGACGCAGGAAATCGGCAGTCGCCATCGCAGAGGTATaaaa                                                                                                                                                              <  2:13080/155‑1 (MQ=255)                                                                                                                                                                              ccGCCAACAATCACGAATGCAGTGGTCACATCG‑TAATTAGTGAACCCGGCTTTGTACAGTCCATAAATCAGTGCGCCTTCGCCCACCACATACGCTTTCTTGCCTTCCTGGCGACGCAGGAAATCGGCAGTCGCCATCGCAGAGGTATAAAACACGCTGTCAGGTACATCGACACCTGCGGTGGCAAAGCGGTTCGCCAGATCTTGCCCAGTCTGCGAAGGATAGTTGGTCAGCAACACCAGCGGCAGGcc                                                           <  2:12870/250‑1 (MQ=255)                                                                                                                                                                               cGCCAACAATCACGAAATCAGGGTTCACATCGGTACT‑AGTGAAACCGGCTTTGTACAGTTCATGAATCAGTGCGCCTTCGCCCACCACATACGCTTTCTTGCCTTCCTGGCGACGCAGGAAATCGGCAGTCGCCATCGCAGAGGTATAAAACACGCTGTCAGGTACATCGACACCTGCGGTGGCAAAGCGGTTCGCCAGATCTTGCCCAGTCTGCGAAGGATAGTTGGTCAGCAACACCAGCGGCAGGCCt                                                          >  1:526043/1‑251 (MQ=255)                                                                                                                                                                               cGCCAACAATCACGAAATCAGGGTTCACATCGGTAAT‑AGTGAAACCGGCTTTGTACAGTTCATAAATCAGTGCGCCTTCGCCCACCACATACGCTTTCTTGCCTTCCTGGCGACGCAGGAAATCGGCAGTCGCCATCGCAGAGGTATAAAACACGCTGTCAGGTACATCGACACCTGCGGTGGCAAAGCGGTTCGCCAGATCTTGCCCAGTCTGCGAAGGATAGTTGGTCAGCAACACCAGCGGCAGGCCt                                                          <  1:92675/251‑1 (MQ=255)                                                                                                                                                                                cccAACAATCGTGAAATCAGGTGTCACAGCGGTAAT‑AGTGAAACCGTCTTTGTACAGTACATAAATCAGTGCGCCTTCGCCCACCACATACGATTCCTTGCCTGCCTGGCGACGCAGGAAATCTGCAGTCGCCCTCGCAGAGGTATAAAACACGCTGTCAGGTACATCGACACCTGCGGTGGCAAAGCGGTTCGCCAGATCTTGCCCGGTCTGCGAAGGATAGTTGGTCAGCAACACCAGCGGCAGGCCtt                                                         <  2:60344/250‑1 (MQ=255)                                                                                                                                                                                        tCACGAAATCAGGGTTCACATCGGTAAT‑AGTGAAACCGGCTTTGTACAGTTCATAAATCAGTGCGCCTTCGCCCACCACATACGCTTTCTTGCCTTCCTGGCGACGCAGGAAATCGGCAGTCGCCATCGCAGCGGTATAAAACACGCTGTCAGGTACATCGACACCTGCGGTGGCAAAGCGGTTCGCCAGATCTTGCCCAGTCTGCGAAGGATAGTTGGTCAGCACCCCCAGCGGCAGGCCTTTAACCAtt                                                 >  2:131561/1‑250 (MQ=255)                                                                                                                                                                                         cACGAAATCAGGGTTCACATCGGTAAT‑AGTGAAACCGGCTTTGTACAGTTCATAAATCAGTGCGCCTTCGCCCACCACATACGCTTTCTTGCCTTCCTGGCGACGCAGGAAATCGGCAGTCGCCATCGCAGAGGTATAAAACACGCTGTCAGGTACATCGACACCTGCGGTGGCAAAGCGGTTCGCCAGATCTTGCCCAGTCTGCGAAGGATAGTTGGTCAGCAACACCAGCGGCa                                                               >  1:488958/1‑236 (MQ=255)                                                                                                                                                                                                             cGGTAAT‑AGTGAAACCGGCTTTGTACAGTTCATAAATCAGTGCGCCTTCGCCCACCACATACGCTTTCTTGCCTTCCTGGCGACGCAGGAAATCGGCAGTCGCCATCGCAGAGGTATAAAACACGCTGTCAGGTACATCGACACCTGCGGTGGCAAAGCGGTTCGCCAGATCTTGCCCAGTCTGCGAAGGATAGTTGGTCAGCAACACCAGCGGCAGGCCTTTATCCATAATCCCGTGCAAAAATTCCGCt                            >  1:306715/1‑251 (MQ=255)                                                                                                                                                                                                                     agggAAACCGGCTTTGTACAGTTCATAAATCAGTGCGCCTTCGCCCACCACAGTCGCTTTCTTGCCTTCCTGGCGACGCAGGAAATCGGCAGTCGCCATCGCAGAGGTATAAAACACGCTGTCAGGTACATCGACACCTGCGGTGGCAAAGCGGTTCGCCAGATCTTGCCCAGTCTGCGAAGGATAGTTGGTCAGCAACACCAGCGGCAGGCCTTTATCCATAATCCCGTGCAAAAATTCCGCTGCACCCg                     <  2:249522/248‑1 (MQ=255)                                                                                                                                                                                                                        gAAACCGGCTTTGTACAGTTCATAAATCAGTGCGCCTTCGCCCACCACATACGCTTTCTTGCCTTCCTGGCGACGCAGGAAATCGGCAGTCGCCATCGCAGAGGTATAAAACACGCTGTCAGGTACATCGACACCTGCGGTGGCAAAGCGGTTCGCCAGATCTTGCCCAGTCTGCGAAGGATAGTTGGTCAGCAACACCAGCGGCAGGCCTTTATCCATAATCCCGTGCAAAAATTCAGCTGCACCCGGt                   >  1:464686/1‑250 (MQ=255)                                                                                                                                                                                                                          aaCCGGCTTTGTACAGTTCATAAATCAGTGCGCCTTCGCCCACCACATACGCTTTCTTGCCTTCCTGGCGACGCAGGAAATCGGCAGTCGCCATCGCAGAGGTATAAAACACGCTGTCAGGTACATCGACACCTGCGGTGGCAAAGCGGTTCGCCAGATCTTGCCCAGTCTGCGAAGGATAGTTGGTCAGCAACACCAGCGGCAGGCCTTTATCCATAATCCCGTGCAAAAATTCCGCTGCACCCGGTACg                <  1:242566/251‑1 (MQ=255)                                                                                                                                                                                                                             cGGCTTTGTACAGTTCATAAATCAGTGCGCCTTCGCCCACCACATACGCTTTCTTGCCTTCCTGGCGACGCAGGAAATCGGCAGTCGCCATCGCAGAGGTATAAAACACGCTGTCAGGTACATCGACACCTGCGGTGGCAAAGCGGTTCGCCCGATCTTGCCCAGTCTGCGAAGGATAGTTGGTCa                                                                              >  2:603764/1‑186 (MQ=255)                                                                                                                                                                                                                             cGGCTTTGTACAGTTCATAAATCAGTGCGCCTTCGCCCACCACATACGCTTTCTTGCCTTCCTGGCGACGCAGGAAATCGGCAGTCGCCATCGCAGAGGTATAAAACACGCTGTCAGGTACATCGACACCTGCGGTGGCAAAGCGGTTCGCCAGATCTTGCCCAGTCTGCGAAGGATAGTTGGTCa                                                                              <  1:603764/186‑1 (MQ=255)                                                                                                                                                                                                                             cGGCTTTGTACAGTTCATAAATCAGTGCGCCTTCGCCCACCACATACGCTTTCTTGCCTTCCTGGCGACGCAGGAAATCGGCAGTCGCCATCGCAGAGGTATAAAACACGCTGTCAGGTACATCGACACCTGCGGTGGCAAAGCGGTTCGCAAGATCTTGCCCAGTCTGCGAAGGATAGTTGGTCAGCAACACCAGCGGCAGGCCTTTATCCCTAATCCCGTGCAAAAATGCCGCCGCACCCGGTACGGCg             >  2:104421/1‑251 (MQ=255)                                                                                                                                                                                                                                 tttGTACAGTTCATGAATCAGTGCGCCTTCGCCCACCACATACGCTTTCTTGCCTTCCTGGCGACGCAGGAAATCGGCAGTCGCCATCGGAGAGGTATAAAACACGCTGTCAGGTACATCGACACCTGCGGTGGCAAAGCGGTTCGCCAGATCTTGCCCAGTCTGCGAAGGATAGTTGGTCAGCAACACCAGCGGCAGGCCTTTATCCATAATCCCGTGCAAAAATTCCGCTGCACCCGGTCCGGCGACGt         >  1:315099/1‑251 (MQ=255)                                                                                                                                                                                                                                        aGTTCATAAATCAGTGCGCCTTCGCCCACCACATACGCTTTCTTGCCTTCCTGGCGACGCAGGAAATCGGCAGTCGCCATCGCAGAGGTATAAAACACGCTGTCAGGTACATCGACACCTGCGGTGGCAAAGCGGTTCGCCAGATCTTGCCCAGTCTGCGAAGGATAGTTGGTCAGCAACACCAGCGGCAGGCCTTTATCCATAATCCCGTGCAAAAATTCCGCTGCACCCGGTACGGCGACGTTATCGTg  <  1:131561/251‑1 (MQ=255)                                                                                                                                                                                                                                           tCATAAATCAGTGCGCCTTCGCCCACCACATACGCTTTCTTGCCTTCCTGGCGACGCAGGAAATCGGCAGTCGCCATCGCAGAGGTATAAAACACGCTGTCAGGTACATCGACACCTGCGGTGGCAAAGCGGTTCGCCAGATCTTGccc                                                                                                     <  1:603836/149‑1 (MQ=255)                                                                                                                                                                                                                                           tCATAAATCAGTGCGCCTTCGCCCACCACATACGCTTTCTTGCCTTCCTGGCGACGCAGGAAATCGGCAGTCGCCATCGCAGAGGTATAAAACACGCTGTCAGGTACATCGACACCTGAGGGGGCAAAGCGGTTCGCCAGGTCTTGccc                                                                                                     >  2:603836/1‑149 (MQ=255)                                                                                                                                                                                                                                               |                                                                                                                                                                                                                                                    TACCAACATAGAACGGTTTGCGCCCGGAGATTTTCTCAATCCCTGCACACAACGCGCCACAAGCGGGATAAAAACCGCGCCCGTGGGTGTCCGGATTGGTGGCGATAAAACGTGCACCGTTAGCGACGAAATAGGCTGCTTTATGCATCATGTCCCAGTTGTAGGAACGCGTTTCGCCAACAATCACGAAATCAGGGTTCACATCGGTAAT‑AGTGAAACCGGCTTTGTACAGTTCATGAATCAGTGCGCCTTCGCCCACCACATACGCTTTCTTGCCTTCCTGGCGACGCAGGAAATCGGCAGTCGCCATCGCAGAGGTATAAAACACGCTGTCAGGTACATCGACACCTGCGGTGGCAAAGCGGTTCGCCAGATCTTGCCCAGTCTGCGAAGGATAGTTGGTCAGCAACACCAGCGGCAGGCCTTTATCCATAATCCCGTGCAAAAATTCCGCTGCACCCGGTACGGCGACGTTATCGTG  >  NC_000913/699801‑700281

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