Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	NC_000913	789,705	Δ1 bp	coding (275/1149 nt)	galK ←	galactokinase

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	789,705	0	C	.	94.5%	230.1 / 2.3	55	coding (275/1149 nt)	galK	galactokinase
Reads supporting (aligned to +/- strand):  ref base C (1/2);  new base . (21/31);  total (22/33)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

TGATCCCGCAGTTACAGCCTACAAACTGGTTTTCTGCTTCCTGACCGTTAAGCGCGATTTGTGCGCCGTCCAGCGGCAGATGATAAAGCTGCTGCAATACGGTTCCGACCGCGACTTCCAGTGAAGCGGAAGAACTTAACCCGGCACCCTGCGGCACATTGCCGCTGATCACCATGTCCACGCCGCCGAAGCTGTTGTTACGCAGTTGCAGATGTTTCACCACGCCACGAACGTAGTTAGCCCATTGATAGTTTTCATGTGCGACAATGGGCGCATCGAGGGAAAACTCGTCGAGCTGATTTTCATAATCGGCTGCCATCACGCGAACTTTACGGTCATCGCGTGGTGCACAACTGATCACGGTTTGATAATCAATCGCGCAGGGCAGAACGAAACCGTCGTTGTAGTCGGTGTGTTCACCAATCAAATTCACGCGGCCAGGCGCCTGAATGGTGTGAGTGGCAGGGTAGCCAAATGCGTTGGCAAACAGAG  >  NC_000913/789463‑789954                                                                                                                                                                                                                                                   |                                                                                                                                                                                                                                                          tGATCCCGCAGTTACAGCCTACAAACTGGTTTTCTGCTTCCTGACCGTTAAGCGCGATTTGTGCGCCGTCCAGCGGCAGATGATAAAGCTGCTGCAATAAGGTTCCGACCGCGACTTCCAGTGAAGCGGAAGAACTTAACCAGGCACCCTGCGGCACATTGCCGCTGATCAACATGTACACGCCGCCGAAGCTGTTGTTACGCAGTTGCAGTTGTTTCACCACGACACGAACGTAGTTAGCC‑ATTGATAGt                                                                                                                                                                                                                                                  >  1:78886/1‑251 (MQ=255)   aTCCCGCAGTTACAGCCTACAAACTGGTTTTCTGCTTCCTGACCGTTAAGCGCGATTTGTGCGCCGTCCAGCGGCAGATGATAAAGCTGCTGCAATACGGTTCCGACCGCGACTTCCAGTGAAGCGGAAGAACTTAACCCGGCACCCTGCGGCACATTGCCGCTGATCACCATGTCCACGCCGCCGAAGCTGTTGTTACGCAGTTGCAGATGTTTCACCACGCCACGAACGTAGTTAGCC‑ATTGATAGttt                                                                                                                                                                                                                                                <  1:496370/251‑1 (MQ=255)        gCAGTTACAGCATACAAACTGGTTTTCTGCTTCCTGACCGTTAAGCGCGATTTGTGCGCCGTCCAGCGGCAGATGATAAAGCTGCTGCAATACGGTTCCGACCGCGACTTCCAGTGAAGCGGAAGAACTTAACCCGGCACCCTGCGGCACCTTGCCGCTGATCACCATGTCCACGCCGCCGAAGCTGTTGTTACGCAGTTGCAGATGTTTCACCACGCCACGAACGTAGTTAGCCC‑TTGCTAGTTTTcttg                                                                                                                                                                                                                                           >  2:280536/1‑248 (MQ=255)         aaGTTACAGCCTACAAACTGGTTTTCTGCTTCCGGACCGTTAAGCGCGATTTGTGCGCCGTCCAGCGGCAGATGATAAAGCTGATGCAATACGGTTCCGACCGCGACTTCCAGTGAAGCGGAAGAACTTAACCCGGCTCACTGCGGCACATTGCCGCTGATCACCATGTCCACGCCGCCGAAGCTGTTGTTACGCAGTTGCAGATGTTTCACCACGCCACGAACGTAGTTAGCC‑ATTGATAGTTTTCAtgt                                                                                                                                                                                                                                          <  2:2199/250‑1 (MQ=255)           gTTACAGCCTACAAACTGGTTTTCTGCTTCCTGACCGTTAAGCGCGATTTGTGCGCCGTCCAGCGGCAGATGATAAAGCTGCTGCAATACGGTTCCGACCGCGACTTCCAGTGAAGCGGAAGAACTTAACCCGGCACCCTGCGGCACATTGCCGCTGATCACCATGTCCACGCCGCCGAAGCTGTTGTTACGCAGTTGCAGATGTTTCACCACGCCACGAACGTAGTTAGCC‑ATTGATAGTTTTCATGTGc                                                                                                                                                                                                                                        <  1:253540/251‑1 (MQ=255)                aGCCTACAAACTGGTTTTCTGCTTCCTGACCGTTAAGCGCGATTTGTGCGCCGTCCAGCGGCAGATGATAAAGCTGCTGAAATACGGTTCCGACCGCGACTTCCAGTGAAGCGGAAGAACTTAACCCGGCACCCTGCGGCACATTGCCGCTGATCACCATGTCCACGCCGCCGAAGCTGTTGTTACGCAGTTGCAGATGTTTCACCACGCCACGAACGTAGTTAGCC‑ATTGATAGTTTTCATGTGCGACa                                                                                                                                                                                                                                    >  1:361026/1‑250 (MQ=255)                     aCAAACTGGTTTTCTGCTTCCTGACCGTTAAGCGCGATTTGTGCGCCGTCCAGCGGCAGATGATAAAGCTGCTGCAATACGGTTCCGACCGCGACTTCCAGTGAAGCGGAAGAACTTAACCCGGCACCCTGCGGCACATTGCCGCTGATCACCATGTCCACGCCGCCGAAGCTGTTGTTACGCAGTTGCAGATGTTTCACCACGCCACGAACGTAGTTAGCC‑ATTGATAGTTTTCATGTGCGCCAATGGgc                                                                                                                                                                                                                              >  1:494329/1‑251 (MQ=255)                     aCAAACTGGTTTTCTGCTTCCTGACCGTTAAGCGCGATTTGTGCGCCGTCCAGCGGCAGATGATAAAGCTGCTGCAATACGGTTCCGACCGCGACTTCCAGTGAAGCGGAAGAACTTAACCCGGCACCCTGCGGCACATTGCCGCTGATCACCATGTCCACGCCGCCGAAGCTGTTGTTACGCAGTTGCAGATGTTTCACCACGCCACGAACGTAGTTAGCC‑ATTGATAGTTTTCATGTGCGACAATGGgc                                                                                                                                                                                                                              <  1:23159/251‑1 (MQ=255)                      cAAACTGGTTTTCTGCTTCCTGACCGTTAAGCGCGATTTGTGCGCCGTCCAGCGGCAGATGATAAAGCTGCTGCAATACGGTTCCGACCGCGACTTCCAGTGAAGCGGAAGAACTTAACCCGGCACCCTGCGGCACATTGCCGCTGATCACCATGTCCACGCCGCCGAAGCTGTTGTTACGCAGTTGCAGATGTTTCACCACGCCACGAACGTAGTTAGCC‑AGTGATAGTGTTCATGTGCGACAATGGgcg                                                                                                                                                                                                                             <  2:538327/251‑1 (MQ=255)                           cggTTTTCTGCTTCCTGACCGTTAAGCGCGATTTGGGCGCCGTCCAGCGGCAGATGATAAAGCTGCTGCAATACGGTTCCGACCGCGACTTCCAGTGAAGCGGAAGAACTTAACCCGGCACCCTGCGGCACACTCCCGCTGATCACCATGTCCACGCCTCCGAAGTTGTTGTTACGCAGTTGCAGATGTTTCACCACGCCACGAACGTAGTTAGCC‑ATTGATAGTTTTCATGTGCGACAATGGGCGCATCg                                                                                                                                                                                                                        <  2:361026/250‑1 (MQ=255)                                        ccggACCGTTAAGCGCGATTTGTGCGCCGGCCAGCGGCAGATGATAAATCTGCTGCAATACGGTTCCGACCGCGAGTTCCAGTGAAGCGGAAGAACTTAACCCGGCACCCTGCGGCACATTGCCGCTGAGCACCATGTCCACGCCGCCGAAGCTGTTTTTACGCAGTTGCAGATGTTTCACCACGCCTCGAACGTAGTTAGCC‑AGTGATAGTTTTCATGTGCGACAATGTGCGCATCGAGTGAAAACtcg                                                                                                                                                                                                            <  2:545470/247‑1 (MQ=255)                                              cGTTAAGCGCGATTGGGGCGCCGTCCAGCGGCAGATGATAAAGCTGCTGCAATACGGTTCCGACCGCGACGTCCAGTGAAGCGGAAGAACTTAACCCGGCACCCTGCGGCACATTGCCGCTGATCACCATGTCCACGCCGCCGAAGCTGTTGTTACGCAGTTGCAGATGTTTCACCACGCCACGAACGTAGTTAGCC‑ATTGATAGTTTTCATGTGCGACAATGGGCGCATCGAGGGAAAACTCGTCGAGCt                                                                                                                                                                                                     <  1:545936/251‑1 (MQ=255)                                                   aGCGCGGTTTGTCCGCCGACCAGAGGCAGATGAGAAAGCTGCTGCAATACGGTTCCGACCGCGACTTCCAGTGAAGCGGAAGAACTTAACCCGACACCCTGCGGCACATTGCCGCTGATCTCCATGTCCACGCCGCCGAAGCTGTTGTTCCGCAGTTGCAGATGTTTCACCACGCCACGAACGTAGTTAGCC‑ATTGGTAGTTTTCATGTGCGACAATGGGCGCATCGCGGGAAAACTCGTAGAGCTGAttt                                                                                                                                                                                                >  2:390164/1‑251 (MQ=255)                                                   aGCGCGATTTGTGCGCCGTCCAGCGGCAGATGATAAAGCTGCTGCAATACGGTTCCGACCGCGACTTCCAGTGAAGCGGAAGAACTTAACCCGGCACCCTGCGGCACATTGCCGCTGATCACCATGTCCACGCCGCCGAAGCTGTTGTTACGCAGTTGCAGATGTTTCACCACGCCACGAACGTAGTTAGCC‑ATTGATAGTTTTCATGTGCGACAATGGGCGCATCGAGGGAAAACTCGTCGAGCTGAttt                                                                                                                                                                                                >  2:153128/1‑251 (MQ=255)                                                               gcgcCGTCCAGCGGCATATGATAACTCTGCTGCAATACTGTTCCGACCGCGACTTACAGTGAAGCTGAAGAACTTAACCCGGCCCCCTGCGGCACATTGACGCTGATCACCATTTCCACGCCGCCTAAGCTGTTGTTACTCAGTTGCAGATGTTTCACCACGCCACGAACGTAGTTAGCC‑ATTGATAGTTTTCATGTGCGACAATGGGCGCATCGAGGGAAAACTCTTCGAGCTGATTTTCATAATCGGCt                                                                                                                                                                                    <  1:24481/251‑1 (MQ=255)                                                                   cGTCCAGCGGCAGATGGTACAGCTCCTGCAATACGGCTCCGACCGCGACTTCCAGTGAAACGGAAGAAATTAACCCGGCACATTGAGGCACATTGCCGCTGATCACAATGTCCACGCCGCAGCAGCTGTTGTTCCGCCGTTGCAGATGTTTCACCACGCCACGAACGTAGTTAGCC‑ATTGATAATTTTCATGTGCGACAATGGGCGCATAGAGGGAACACTCGACGCGCTGATTTTCATCATCGGCTGcc                                                                                                                                                                                 >  2:28640/1‑250 (MQ=255)                                                                         gCGGCAGATGATAAAGCTGCTGCAATACGGTTCCGACCGCGACTTCCAGTGCAGCGGAAGAACTTAACCCGGCACCCTGCGGCACATTGCCGCTGATCACCATGTCCACGCCGCCGAAGCTGTTGTTACGCAGTTGCAGATGTTTCACCACGCCACGAACGTAGTTAGCC‑ATTGATAGTTTTCATGTGCGACAATGGGCGCATCGAGGGAAAACTAGTCGAGCTGATTTTCATACTCGGCTGCCATCCcgc                                                                                                                                                                          >  2:360241/1‑251 (MQ=255)                                                                          cGGCAGATGATAAAGCTGCTGCAATACGGTTCCGACCGCGACTTCCAGTGAAGCGGAAGAACTTAACCCGGCACCCTGCGGCACATTGCCGCTGATCACCATGTCCACGCCGCCGAAGCTGTTGTTACGCAGTTGCAGATGTTTCACCACGCCACGAACGTAGTTAGCC‑ATTGATAGTTTTCATGTGCGACAATGGGCGCATCGAGGGAAAACTCGTCGAGCTGATTTTCATAATCGGCTGCCATCAcgcg                                                                                                                                                                         <  1:104791/251‑1 (MQ=255)                                                                                atgatAAAGCTGCTGCAATACGGTTCCGACCGCGACTTCCAGTGAAGCGGAAGAACTTAACCCGGCACCCTGCGGCACATTGCCGCTGATCACCATGTCCACGCCGCCGAAGCTGTTGTTACGCAGTTGCAGATGTTTCACCACGCCACGAACGTAGTTAGCC‑ATTGATAGTTTTCATGTGCGACAATGGGCGCATCGAGGGAAAACTCGTCGAGCTGATTTTCATAATCGGCTGCCATCACGCGAACttt                                                                                                                                                                   >  2:332931/1‑251 (MQ=255)                                                                                           gctgcAATACGGTTCCGACCGCGACTTCCAGTGAAGCGGAAGAACTTAACCCGGCACCCGGCGGCACATTGCCGCTGATCACCATGTCCACGCCGCCGAAGCTGTTGTTACGCAGTTGCAGATGTTTCACCACGCCACGAACGTAGTTAGCC‑ATTGATAGTTTTCATGTGCGACAATGGGCGCATCGAGGGAAAACTCGTCGAGCTGATTTTCATAATCGGCTGCCATCACGCGAACTTTACGGTCATcgc                                                                                                                                                        <  1:468742/251‑1 (MQ=255)                                                                                                aaTACGGTTCCGACCGCGACTTCCAGTGAAGCGGAAGAACTTAACCCGGCACCCTGCGGCACATTGCCGCTGATCACCATGTCCACGCCGCCGAAGCTGTTGTTACGCAGTTGCAGATGTTTCACCACGCCACGAACGTAGTTAGCC‑ATTGATAGTTTTCATGTGCGACAATGGGCGCATCGAGGGAAAACTCGTCGAGCTGATTTTCATAATCGGCTGCCATCACGCGAACTTTACGGTCATCGCgtggt                                                                                                                                                   >  1:252778/1‑251 (MQ=255)                                                                                                      gTTCCGACCGCCACTTACAGCGAAGCGGCAGAACTTAACCAGGCACCCAGCAGCACATTGCCGCTGATCACAATGAACACGCCGCCGACGCTGATGTTACGCAGTTGCAGATGTTTCACCACGCCCCGAACGTAGTTAGCC‑ATTGATAGTTTTCATGTGCGACAATGGGCGCATCGAGGGAAAACTCGCCGAGCTGATTTTCATACTCGGCTGCCATCACGCGACCTTTACGGTCATCGAGTGGTGcaca                                                                                                                                              >  2:549999/1‑250 (MQ=255)                                                                                                       ttCCGACCGCGACTTCCAGTGAAGCGGAAGAACTTAACCCGGCACCCTGCGGCACATTGCCGCTGATCACCATGTCCACGCCGCCGAAGCTGTTGTTACGCAGTTGCAGATGTTTCACCACGCCACGAACGTAGTTAGCC‑ATTGATAGTTTTCATGTGCGACAATGGGCGCATCGAGGGAAAACTCGTCGAGCTGATTTTCATAATCGGCTGCCATCACGCGAACTTTACGGTCATCGCGTGGTGCACAAc                                                                                                                                            >  1:119245/1‑251 (MQ=255)                                                                                                         ccGACCGCGACTTCCAGTGAAGCGGAAGAACTTAACCCGGCACCCTGCGGCACATTGCCGCTGATCACCATGTCCACGCCGCCGAAGCTGTTGTTACGCAGTTGCAGATGTTTCACCACGCCACGAACGTAGTTAGCC‑ATTGATAGTTTTCATGTGCGACAATGGGCGCATCGAGGGAAAACTCGTCGAGCTGATTTTCATAATCGGCTGCCATCACGCGAACTTTACGGTCATCGCGTGGTGCACAACTg                                                                                                                                          <  1:516205/251‑1 (MQ=255)                                                                                                          cGACCGCGACTTCCAGTGAAGCGGAAGAACTTAACCCGGCCCCCTGCGGCACATTGCCGCTGATCACCATGTCCACGCCGCCGAAGCTGTTGTTACGCAGTTGCAGATGTTTCACCACGCCACGAACGTAGTTAGCC‑ATTGATAGTTTTCATGTGCGACAATGGGCGCATCGAGGGAAAACTCGTCGAGCTGATTTTCATAATCGGCTGCCATCACGCGAACTTTACGGTCATCGCGTGGTGCACAACTGa                                                                                                                                         <  2:370097/251‑1 (MQ=255)                                                                                                           gACCGCGACTTCCAGTGAAGCGGAAGAACTTAACCCGGCACCCTGCGGCACATTGCCGCTGATCACCATGTCCACGCCGCCGAAGCTGTTGTTACGCAGTTGCAGATGTTTCACCACGCCACGAACGTAGTTAGCC‑ATTGATAGTTTTCATGTGCGACAATGGGCGCATCGAGGGAAAACTCGTCGAGCTGATTTTCATAATCGGCTGCCATCACGCGAACTTTACGGTCATCGCGTGGTGCACAACTGAt                                                                                                                                        <  2:284556/251‑1 (MQ=255)                                                                                                                     tCCAGTGAAGCGGAAGAACTTAACCCGGCACCCTGCGGCACATTGCCGCTGATCACCATGTCCACGCCGCCGAAGCTGTTGTTACGCAGTTGCAGATGTTTCACCACGCCACGAACGTAGTTAGCC‑ATTGATAGTTTTCATGTGCGACAATGGGCGCATCGAGGGAAAACTCGTCGAGCTGATTTTCATAATCGGCTGCCATCACGCGAACTTTACGGTCATCGCGTGGTGCACAACTGATCACGGTTTg                                                                                                                               <  2:192964/250‑1 (MQ=255)                                                                                                                         gTGAAGCGGAAGAACTTAACCCGGCACCCTGCGGCACATTGCCGCTGATCACCATGTCCACGCCGCCGAAGCTGTTGTTACGCAGTTGCAGATGTTTCACCACGCCACGAACGTAGTTAGCC‑ATTGATAGTTTTCATGTGCGACAATGGGCGCATCGAGGGAAAACTCGTCGAGCTGATTTTCATAATCGGCTGCCATCACGCGAACTTTACGGTCATCGCGTGGTGCACAACTGATCACGGTTTGATaat                                                                                                                          <  2:302257/251‑1 (MQ=255)                                                                                                                                  aagaaCTTAACCCGGCACCCTGCGGCACAGTGCCGCTGATCACCATGTCCCCGCCGCCGAAGCTGTTGTTACGCAGTTGCAGATGTTTCACCACGCCACGAACGTAGTTAGCC‑ATTGATAGTTGTCATGTGCGACAATGGGCGCATCGAGGGAAAACTCGTCGAGCTGATTTTCATAATCGGCTGCCATCACGCGAACTTTACGGTCATCGCGTGGTGCACAACTGATCACGGTTTGATAATCAATCgcgc                                                                                                                 <  2:470136/251‑1 (MQ=255)                                                                                                                                      aCTTAACCCGGCACCCGGCGGCACATTGCCGCTGATCACCATGTCCACGCCGCCGAAGCTGTTGTTACGCAGTTGCAGATGTTTCACCACGCCACGAACGTAGTTAGCC‑ATTGATAGTTTTCATGTGCGACAATGGGCGCATCGAGGGAAATCTCGTCGAGCTGATTTTCATAATCGGCTGCCATCACGCGAACTTTACGGTCATCGCGTGGTGCACAACTGATCACGGTTTGATAATCAATCGCGCAggg                                                                                                             <  1:188677/251‑1 (MQ=255)                                                                                                                                          aaCCCGGCCCCCTGCGGCACATTGCCGCTGATCACCATGTCCACGCCGCCGAAGCTGTTGTTACGCAGTTGCAGATGTTTCACCACGCCACGAACGTAGTTAGCC‑ATTGATAGTTTTCATGTGCGACAATGGGCGCATCGAGGGAAAACTCGTCGAGCTGATTTTCATAATCGGCTGCCATCACGCGAACTTTACGGTCATCGCGTGGTGCACAACTGATCACGGTTTGATAATCAATCGCGCAGGGCAg                                                                                                          <  2:198757/250‑1 (MQ=255)                                                                                                                                               ggCACCCTGCGGCACATTGCCGCTGATCACCATGTCCACGCCGCCGAAGCTGTTGTTACGCAGTTGCAGATGTTTCACCACGCCACGAACGTAGTTAGCC‑ATTGATAGTTTTCATGTGCGACAATGGGCGCATCGAGGGAAAACTCGTCGAGCTGATTTTCATAATCGGCTGCCATCACGCGAACTTTACGGTCATCGCGTGGTGCACAACTGATCACGGTTTGATAATCAATCGCGCAGGGCAGAACGaa                                                                                                    <  1:562016/251‑1 (MQ=255)                                                                                                                                                gCACCCTGCTGCACATTGCCGCTGATCACCATGTCCAAGCCGCCGAAGCTGTTGTTACGCAGTTGCAGATGTTTCACCACGCCACGAACATAGTTAGCC‑ATTGATAGTTTTCATGTGCGACAATGGACGCATCGATGAAACACTCGACGAGCTGATTTTCATAATCAGCTGCCATCACGCGCACTTTACGTTCATCGCGTGGTGCGCAACTGATCACGTTTTGGTAATCAATCGCGCAGGGCACAACGaac                                                                                                   >  1:298337/1‑250 (MQ=255)                                                                                                                                                       gCGGCACATTGCCGCGGATCACCAGGTCCCCGCCGCCGAAGCTGTTGTTACGCAGTTGCAGATGTTTCACCACGCCACGAACGTAGTTAGCC‑ATTGATAGTTTTCATGTGCGACAATGGGCGCAGCGAGTGAAAACTCGGCGAGCTGATTTTCGTAATCGGCTGCCAGCACGCGAACTTTACGGGCATCGCGTGGTGCACAACTGATCTCGGTGTTAGAATCAGTCGCGCAGGGCAGAACAAAACcgtc                                                                                              <  2:285711/249‑1 (MQ=255)                                                                                                                                                         ggCACATGGCCGCTGAGCACCATGTCCACGCCGCCGAAGCTGTTGTTACGCAGTTGCAGAGGTCTCACCACGCCACGAACGGAGTGAGCC‑ATTGATAGTTTTCAGGTGCGAGAATGTGCGCATCGAGGGAAAACTCGTCGAGCTGATTGTCATAATCGGCTGCCATCACGCTGACTTTACGGTCAGCGCGTGGTGCACAACTGAGCACGGTTTGATAATCAATCGCGCAGGGCTGAACGAAACCGTCGTTg                                                                                          <  1:79092/251‑1 (MQ=255)                                                                                                                                                                                     cgccgcCGAAGCTTTTGTTACGCAGTTGCAGATGTTTCACCACGCCACGAACGTAGTTAGCC‑ATTGATAGTTTTCATGTGCGACAATGGGCGCATCGAGGGAAAACTCGTCGAGCTGATTTTCATAATCGGCTGCCATCACGCGAACTTTACGGTCAGCGCGTGGTGCACAACTGATCACGGTTTGATAATCAATCGCGCAGGGCAGAACGAAACCGTCGTTGTAGTCGGTTTTTTCTCTAATCAAATTCa                                                              <  2:83452/251‑1 (MQ=255)                                                                                                                                                                                           cgAAGCTGTTGTTACGCAGTTGCAGATGTTTCACCACGCCACGAACGTAGTTAGCC‑ATTGATAGTTTTCATGTGCGACACTGGGCGCATCGAGGGAAAACTCGTCGAGCTGATTTTCATAATCGGCTGCCATCACGCGAACTTTACGGTCATCGCGTGGTGCACAACTGATCACGGTTTGATAATCAATCGCGCAGGGCAGAACGAAACCGTCGTCGTAGTCGGTGTGTTCACCAATCACATTCACGCGGc                                                        >  1:553993/1‑251 (MQ=255)                                                                                                                                                                                            gAAGCTGTTGTTACGCAGTTGCAGATGTTTCACCACGCCACGAACGTAGTTAGCC‑ATTGATAGTTTTCATGTGCGACAATGGGCGAATCGAGGGAAAACTCGTCGAGCTGATTTTCATAATCGGCTGCCATCACGCGAACTTTACGGTCATCGCGTGGTGCACAACTGATCACGGTTTGATAATCAATCGCGCAGGGCAGAACGAAACCGTCGTTGTAGTCGGTGTGTTCACCAATCAACTTCACGCGGcc                                                       >  1:472175/1‑251 (MQ=255)                                                                                                                                                                                             aagtagttgttACGCTGGTGCAGAGGTTTCACCACGGCACGAACTGAGTTAGCC‑ATTGATAGTTTTCATGTGCGACAATTGGCGCATCGAGGGAAAACTCGTCGAGCTGATGTTCATAATCGGCTGCCATCACGCGAACTTTACGGGCATCGCGGGGTGCACACCGGATCTCGGTTTGATAATCAATCGCGCAGGGCAGATCGAAAGCGTCGGTGTAGTCGGTGTGTTCACCAATCAAATTCACGCTGCCa                                                      <  2:467408/246‑1 (MQ=255)                                                                                                                                                                                             aaGCTGTTGTTCCGGAGTTGCAGATGTTTCACCACGCCACGACCGTAGTTAGCC‑ACTGATAGTTTTCATGTGCGACAATGGGCGCATCGAGGGAAAACTCGTCGAGCTGATTTTCATACTCGGCTGCCATCACGCGAACTTTACGGCCATCGCGTGGTGCACAACTGATCCCGGTTTGATAATCAATCGCGCAGGGCAGAACGAAACCGTCGTTGTAGTCGGTGTGTTCACCAATCAAATTCACGCGGcc                                                       >  2:351420/1‑250 (MQ=255)                                                                                                                                                                                              aGCTGTTGTTCCGCAGTTGCAGATGTTTCACAACGCCACGAACGTAGTTAGCC‑ATTGATAGTTTTCATGTGCGACAATGGGCGCATCGAGGGAAAACTCGTCGAGCTGATTTTCATAATCGGCTGCCATCACGCGAACTTTACGGTCATCGCGTGGTGCACAACTGATCACGGTTTGATAATCAATCGCGCAGGGCAGAACGAAACCGTCGTTGTAGTCGGTGTGTTCACCAATCAAATTCACGCGGcc                                                       >  2:375875/1‑249 (MQ=255)                                                                                                                                                                                              aGCTGTTGTTACGCAGTTGCAGATGTTTCACCACGCCACGAACGTAGTTAGCC‑ATTGATAGTTTTCATGTGCGACAATGGGCGCATCGAGGGAAAACTCGTCGAGCTGATTTTCATAATCGGCTGCCATCACGCGAACTTTACGGTCATCGCGTGGTGCACAACTGATCACGGTTTGATAATCAATCGCGCAGGGCAGAACGAAACCGTCGTTGTAGTCGGTGTGTTCACCAATCAAATTCACGCGGcc                                                       >  1:552010/1‑249 (MQ=255)                                                                                                                                                                                              aGCTGTTGTTACGCAGTTGCAGATGTTTCACCACGCCACGAACGTAGTTAGCC‑ATTGATAGTTTTCATGTGCGACAATGGGCGCATCGAGGGAAAACTCGTCGAGCTGATTTTCATAATCGGCTGCCATCACGCGAACTTTACGGTCATCGCGTGGTGCACAACTGATCACGGTTTGATAATCAATCGCGCAGGGCAGAACGAAACCGTCGTTGTAGTCGGTGTGTTCACCAATCAAATTCACGCGGCCAg                                                     <  2:394539/251‑1 (MQ=255)                                                                                                                                                                                                      ttACGCAGTTGCAGATGTTTCACCACGCCACGAACGTAGTTAGCC‑ATTGATAGTTTTCATGTGCGACACTGGGCGCATCGAGGGAAAACTCGTCGAGCTGCTTTTCATAATCGGCTGCCATCACGCGAACTTTACGGTCATCGCGTGGTGCACAACTGATCACGGATTGATAATCAATCGCGCAGGGCAGAACGAAACCGTCGTTGTAGTCGGTGTGTTCACCAATCAAATTCACGCGGCCAGGCGCCTg                                              >  1:324016/1‑250 (MQ=255)                                                                                                                                                                                                       tACGCAGTTGCAGATGTTTCACCACGCCACGAACGTAGTTAGCC‑ATTAATAGTTTTCATGTGCGACAATGGGCGCATCGAGGGAACACTCGTCGAGCTGATTTTCATAATCGGCTGCCATCACGCGAACTTTACGGTCATCGCGTGGTGCACAACTGATCACGGTTTGATAATCAATCGCGCAGGGCAGAACGCACCCGTCGTTGTAGTCGGTGTGTTCACCAATCAAATTCACGCGGCCAGGCGCCTGa                                             >  1:502674/1‑250 (MQ=255)                                                                                                                                                                                                        aCGCAGTTGCAGATGTTTCACCACGCCACGAACGTAGTTAGCC‑ATTGATAGTTTTCATGTGCGACAATGGGCGCATCGAGGGAAAACTCGTCGAGCTGATTTTCATAATCGGCTGCCATCACGCGAACTTTACGGTCATCGCGTGGTGCACAACTGATCACGGTTTGATAATCAATCGCGCAGGGCAGAACGAAACCGTCGTTGTAGTCGGTGTGTTCACCAATCAAATTCACGCGGCCAGGCGCCTGAAt                                           <  2:108463/251‑1 (MQ=255)                                                                                                                                                                                                             gTTGCAGATGTTTCACCACGCCACGAACGTAGTTAGCC‑ATTGATAGTTTTCATGTGCGACAATGGGCGCATCGAGGGAAAACTCGTCGAGCTGATTTTCATAATCGGCTGCCATCACGCGAACTTTACGGTCATCGCGTGGTGCACAACTGATCACGGTTTGATAATCAATCGCGCAGGGCAGAACGAAACCGTCGTTGTAGTCGGTGTGTTCACCAATCAAATTCACGCGGCCAGGCGCCTGAATGgtgt                                      <  1:432105/251‑1 (MQ=255)                                                                                                                                                                                                                 caggtGTTTCACCACGCTACGAACGTAGTTAGAC‑ATTGCTAGTTTTCATGTGCGACAATGGGCGCATCGAGGGAAAACTCGTCGAGCTGATTTTCATAATCGGCTGCCATCACGCGAACTTTACGGTCATCGCGAGGTGCACAACTGATCACGATTTGATAATCAATCGCGCAGGGCAGAACGAAACCGTCGTTGTAGTCGGTGTGTTCACCAATCAAATTCACGCGGCCAGGCGCCTGCATGGTGTGAgg                                  >  2:123884/5‑250 (MQ=255)                                                                                                                                                                                                                        ttCACCACGCCACGAACGTAGTTA‑TCCAGTGATAGTTTTCATGTGCGACAATGGGCGCATCGAGGGAAAAGTCGGCGAGCTGATTTTCATAATCGGCTGCCATCACGCGAACTTTACGGTCTTCGCGTGGTGCACAACTGATCTCGTTTGGATAATCAATCGCGCAGGTGAGAACGAAACCGTCGTTGTAGTCTGTGTGTTCACCAATCAAATGCACGCGGCCAGCCGTCTGAATGGTGTGAGTTGCAgg                            <  2:566981/250‑1 (MQ=255)                                                                                                                                                                                                                          caccacGCCACGAACGTAGTTAGCC‑ATTGATAGTTTTCATGTGCGACAATGGGCGCATCGAGGGAAAACTCGTCGAGCTGATTTTCATAATCGGCTGCCATCACGCGAACTTTACGGTCATCGCGTGGTGCACAACTGATCACGGTTTGATAATCAATCGCGCAGGGCAGAACGAAACCGTCGTTGTAGTCGGTGTGTTCACCAATCAAATTCACGCGGCCAGGCGCCTGAATGGTGTGAGTGGCAGGGTa                         <  1:332931/251‑1 (MQ=255)                                                                                                                                                                                                                            ccacgccacgAACGTATTTAGCC‑ATTGATAGTTTTCATGTGCGACAATGGGAGGATCGAGGGAAAACTCGTCGAGCTGATTTTCATAATCGGCTGGCATCACGCGAACGTTTCGGTCATCGCGTGGTGCACAACTGATCACGGTTTGATAATCAATCGCGCAGGGCAGAACGAAACCGTCGTTGTAGTCGGGGTGTTCTCCAATCAAATTCACGCGGCCAGGCGCCTGAATGGTGTGAGTGGCAGGGTAGc                       <  2:487694/251‑1 (MQ=255)                                                                                                                                                                                                                             cacgccacgAACGTAGTTAGCC‑ATTGATAGATTTCATGTGCGACAATGGGCGCATCGAGGGAAAACTCGTCGAGCTGATTTTCATAATCGGCTGCCATCACGCGAACTTTACGGTCATCGCGTGGTGCACAACTGATCACGGTTTGATAATCAATCGCGCAGGGCAGAACGAAACCGTCGTTGTAGTCGGTGTGTTCACCAATCAAATTCACGCGGCCAGGCGCCTGAATGGTGTGAGTGGCAGGGTatcc                      >  2:529400/1‑248 (MQ=255)                                                                                                                                                                                                                                gcctagAACGTAGTTAGCC‑ATTGATAGGTTTCATGTGCGACAATGGGCGCATCGAGGGAAAACTCGTCGAGCTGATTTTCATAATCGGCTGCCATCACGCGAACTTTACGGTCATCGCGTGGTGCACAACTGATCACGGTTTGATAATCAATCGCGCAGGGCAGAACGAAACCGTCGTTGTAGTCGGTGTGTTCACCAATCAAATTCACGCGGCCAGGCGCCTGAATGGTGTGAGTGGCAGGGTAGCCaaa                   <  1:336545/246‑1 (MQ=255)                                                                                                                                                                                                                                  cacgAACGTAGTTAGCC‑ATTGATAGTTTGCATGTGCGACAATGGGCGCATCGAGGGAAAACTCGTCGAGCTGATTTTCATAATCGGCTGCCATCACGCGAACTTTTCGGTCATCGCGTGGTGCACAACTGATCACGGTTTGATAATCAATCGCGCAGGGCAGAACGAAACCGTCGTTGTAGTCGGTGTGTTCACCAATCAAATTCACGCGGCCAGGCGCCGGAATGGTGTGAGTGGCAGGGTAGCCAAATg                 <  1:245886/251‑1 (MQ=255)                                                                                                                                                                                                                                       aCGTAGTTA‑TCCATTGATAGTTTTCATGTGCGACAATGGGCGCATCGAGGGAAAACTCGTCGAGCTGATTTTCATAATCGGCTGCCATCACGCGAACTTTACGGTCATCGCGTGGTGCACAACTGATCACGGTTTGATAATCAATCGCGCAGGGCAGAACGAAACCGTCGTTGTAGTCGGTGTGTTCACCAATCAAATTCACGCGGCCAGGCGCCTGAATGGTGTGAGTGGCAGGGTAGCCAAATGCGTTg            <  1:534335/251‑1 (MQ=255)                                                                                                                                                                                                                                                agccATGGATAGTTTTGACGTGCGACAATGGGCGCAACGAGGGAAAACTCGTCGAGCTGATTTTCATAATCGGCTGCCATCACGCGAACTTTACGGTCATCGCGTGGTGCACAACTGATCACGGTTTGATAATCAATCGCGCAGGGCAGAACGAAACCGTCGTTGTAGTCGGTGTATTCCCCAATCAAATTAACGCGGCCAGGCGCCTGAATGGTGTGAGTGGCAGGGCAACCAAAGGAGTTGGCAAAc      >  2:170156/3‑249 (MQ=255)                                                                                                                                                                                                                                                 gccATTGATAGTTTTCATGTGCGACAATGGGCGCATCGAGGGAAAACTCGTCGAGCTGATTTTCATAATCGGCTGCCATCACGCGAACTTTACGGTCATCGCGTGGTGCACAACTGATCACGGTTTGATAATCAATCGCGCAGGGCAGAACGAAACCGTCGTTGTAGTCGGTGTGTTCACCAATCAAATTCACGCGGCCAGGCGCCTGAATGGTGTGAGAGGCAGGGTAGCCACATGCGTTGGCAAAc      >  1:548412/2‑248 (MQ=255)                                                                                                                                                                                                                                                 gccATTGAAAGTTTTCATGTGCGACAATGGGCGCATCGAGGGAAAACTCGGCGGGCTGAGTTGCATAAGCGGCTGCCATCACGCGAACTTTACGGTCATCGCGTGGTGCACAACTGATCACGGTTTGATAATCAATGGCGCATGGCAGAACGGAACCGTCGTTGTAGGCGGTGTGTTCACCAATGAAATTCACGCGGCGAGGCGGCTGAATGGTGTGAGTGGCAGGGTAGCCAAATGCGTTGGCAAACaga   <  1:351420/250‑1 (MQ=255)                                                                                                                                                                                                                                                  ccATTGATAGTTTTCATGTGCGACAATGGGCGCATCGAGGGAAAACTCGTCGAGCTGATTTTCATAATCGGCTGCCATCACGCGAACTTTACGGTCATCGCGTGGTGCACAACTGATCACGGTTTGATAATCAATCGCGCAGGGCAGAACGAAACCGTCGTTGTAGTCGGTGTGTTCACCAATCAAATGCACGCGGCCAGGCGCCTGAATGGTGTGAGTGGCAGGGTAGCCAAATGCGTTGGCAAACagag  <  2:252778/251‑1 (MQ=255)                                                                                                                                                                                                                                                   |                                                                                                                                                                                                                                                          TGATCCCGCAGTTACAGCCTACAAACTGGTTTTCTGCTTCCTGACCGTTAAGCGCGATTTGTGCGCCGTCCAGCGGCAGATGATAAAGCTGCTGCAATACGGTTCCGACCGCGACTTCCAGTGAAGCGGAAGAACTTAACCCGGCACCCTGCGGCACATTGCCGCTGATCACCATGTCCACGCCGCCGAAGCTGTTGTTACGCAGTTGCAGATGTTTCACCACGCCACGAACGTAGTTAGCCCATTGATAGTTTTCATGTGCGACAATGGGCGCATCGAGGGAAAACTCGTCGAGCTGATTTTCATAATCGGCTGCCATCACGCGAACTTTACGGTCATCGCGTGGTGCACAACTGATCACGGTTTGATAATCAATCGCGCAGGGCAGAACGAAACCGTCGTTGTAGTCGGTGTGTTCACCAATCAAATTCACGCGGCCAGGCGCCTGAATGGTGTGAGTGGCAGGGTAGCCAAATGCGTTGGCAAACAGAG  >  NC_000913/789463‑789954

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