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	.	97.5%	175.6 / ‑3.9	40	coding (275/1149 nt)	galK	galactokinase
Reads supporting (aligned to +/- strand):  ref base C (1/0);  new base . (15/24);  total (16/24)
Fisher's exact test for biased strand distribution p-value = 4.00e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

GATCCCGCAGTTACAGCCTACAAACTGGTTTTCTGCTTCCTGACCGTTAAGCGCGATTTGTGCGCCGTCCAGCGGCAGATGATAAAGCTGCTGCAATACGGTTCCGACCGCGACTTCCAGTGAAGCGGAAGAACTTAACCCGGCACCCTGCGGCACATTGCCGCTGATCACCATGTCCACGCCGCCGAAGCTGTTGTTACGCAGTTGCAGATGTTTCACCACGCCACGAACGTAGTTAGCCCATTGATAGTTTTCATGTGCGACAATGGGCGCATCGAGGGAAAACTCGTCGAGCTGATTTTCATAATCGGCTGCCATCACGCGAACTTTACGGTCATCGCGTGGTGCACAACTGATCACGGTTTGATAATCAATCGCGCAGGGCAGAACGAAACCGTCGTTGTAGTCGGTGTGTTCACCAATCAAATTCACGCGGCCAGGCGCCTGAATGGTGTGAGTGGCAGGGTAGCCAAATGCGTTGGCAAACAG  >  NC_000913/789464‑789952                                                                                                                                                                                                                                                  |                                                                                                                                                                                                                                                        gATCCCGCAGTTACAGCCTACAAACTGGTTTTCTGCTTCCTGACCGTTAAGCGCGATTTGTGCGCCGTCCAGCGGCAGATGATAAAGCTGCTGCAAAACGGTACCGACCGCGACTTCCAGTGAAGCGGAAGAACTTTACCCTTCACCCTGCGGCACATTGCCGCTTATCCCCATGTCCACGCCGCCGAAGCTGTTGTTACGCAGTTGCGGATGGTTCACCACGCCACGAACGTAGTTAGCC‑ATTGATAGtt                                                                                                                                                                                                                                               <  1:438061/251‑1 (MQ=255)                    aCAAACTGGTTTTCTGCTTCCTGACCGTTAAGCGCGATTTGTGCGCCTGCCAGCGGCAGATGATAACGCGGCTGCAATACGGTTCCGCCCGCGCCTTCCAGTGAAGCGGAAGAACTTAACCCGGCACCCTGCGGCACATTGCCGCTGATCACCATGTCCACGCCGCCGAAGCTGTTGTTACGCAGTTGCAGATGTTACACCACGCCACGAACGTAGTTAGCC‑ATTGATAGTTTTCATGTGCGACAATGGgc                                                                                                                                                                                                                            <  2:72569/251‑1 (MQ=255)                        aCTGGTTTTCTGCTTCCTGACCGTTAAGCGCGATTTGTGCGCCGTCCAGCGGCAGATGATAAAGCTGCTGCAATACGGTTCCGACCGCGACTTCCAGTGAAGCGGAAGAACTTAACCCGGCACCCTGCGGCACATTGCCGCTGATCACCATGTCCACGCCGCCGAAGCTGTTGTTACGCAGTTGCAGATGTTTCACCACGCCACGAACGTAGTTAGCC‑ATTGATAGTTTTCATGTGCGACAATGGGCGCAt                                                                                                                                                                                                                        <  1:150801/251‑1 (MQ=255)                           ggTTTTCTGCTTCCTGACAGTTAAGCGCGATTTGTGCGCCGTCCAGCGGCAGATGATAAAGCTGCTGCAATACGGTTCCGACCGCGACTTCCAGTGAAGCGGAAGAACTTAACCCGGCACCCTGCGGCACATTGCCGCTGATCACCATGTCCACGCCGCCGAAGCTGTTGTTACGCAGTTGCAGATGTTTCACAACGCCACGAACGTAGTTAGCCC‑TTGATAGTTTTCATGTGCGACAATGGGCGCATCg                                                                                                                                                                                                                      >  2:119495/1‑250 (MQ=255)                                   gCTTCCTGTCCGTTAAGCGCGATTTGTGCGCCGTTCAGCGGCAGATGATAAAGCTGCTGCAATACGGTTCCGACCGCTACTTCCAGTGAAGCGGAAGAACTTAACCCGGCACCCTGCGGCACATTGCCGCTGATCACCATGTCCACGCCGCCGAAGCTGTTGTTACGCAGTTGCAGATGTTTCACCACGCCACGAACGTAGTTAGCC‑ATTGATAGTTGTCATGTGCGACAATGGGCGCATCGAGGGAAAAc                                                                                                                                                                                                             <  1:367819/251‑1 (MQ=255)                                   gCTTCCTGACCGTTAAGCGCGATTTGTGCGCCGTCCAGCGGCAGATGATAAAGCTGCTGCAATACGGTTCCGACCGCGACTTCCAGTGAAGCGGAAGAACTTAACCCGGCACCCTGCGGCACATTGCCGCTGATCACCATGTCCACGCCGCCGAAGCTGTTGTTACGCAGTTGCAGATGTTTCACCACGCCACGAACGTAGTTAGCC‑ATTGATAGTTTTCATGTGCGACAATGGGCGCATCGAGGGaaa                                                                                                                                                                                                               <  1:412570/249‑1 (MQ=255)                                              gTTAAGCGCGATTTGTGCGCCGTCCAGCGGCAGATGATAACGCTGCTGCAAAACGGTTCCGACCGCGACTTCCAGTGAAGCGGAAGAACTTAACCCAGCACCCTGCGGCACATTGCCGCTGATCACCATGTCCACGCCGCCGAAGCTGTTGTTACGCAGTTGCAGATGTTTCACCACGCCACGAACGAAGTTAGCC‑ATTGATAGTTTTCATGTGCGACAATGGGCGCATCGAGGGAAAACTCGTCGAGCTg                                                                                                                                                                                                  >  1:461076/1‑251 (MQ=255)                                                tAAGCGCGATTTGTGCGCCGTCCAGCGGCAGATGATAAAGCTGCTGCAATACGGTTCCGACCGCGACTTCCAGTGAAGCGGAAGAACTTAACCCGGCACCCTGCGGCACATTGCCGCTGATCACCATGTCCACGCCGCCGAAGCTGTTGTTACGCAGTTGCAGATGTTTCACCACGCCACGAACGTAGTTAGCC‑ATTGATAGTTTTCATGTGCGACAATGGGCGCATCGAGGGAAAACTCGTCGAGCTGa                                                                                                                                                                                                 <  1:113073/250‑1 (MQ=255)                                                     gcgATTTGTGCGCCGTCCAGCGGCAGATGATAAAGCTGCTGCAATACGGTTCCGACCGCGACTTCCAGTGAAGCGGAAGAACTTAACCCGGCACCCTGCGGCACATTGCCGCTGATCACCATGTCCACGCCGCCGAAGCTGTTGTTACGCAGTTGCAGATGTTTCACCACGCCACGAACGTAGTTAGCC‑ATTGATAGTTTTCATGTGCGACAATGGGCGCATCGAGGGAAAACTCGTCGAGCTGATTTTc                                                                                                                                                                                            >  2:110591/1‑250 (MQ=255)                                                            gtgCGCCGTCCAGCGGCAGATGATAAAGCTGCTGCAATACGGTTCCGACCGCGACTTCCAGTGAAGCGGAAGAACTTAACCCGGCACCCTGCGGCACATTGCCGCTGATCACCATGTCCACGCCGCCGAAGCTGTTGTTACGCAGTTGCAGATGTTTCACCACGCCACGAACGTAGTTAGCC‑ATTGATAGTTTTCATGTGCGACAATGGGCGCATCGAGGGAAAACTCGTCGAGCTGATTTTCATAATCgg                                                                                                                                                                                    <  1:79776/251‑1 (MQ=255)                                                                                 gatAAAGCTGCTGCAATACGGTTCCGACCGCGACTTCCAGTGAAGCGGAAGAACTTAACCCGGCACCCTGCGGCACATTGCCGCTGATCACCATGTCCACGCCGCCGAAGCTGTTGTTACGCAGTTGCAGATGTTTCACCACGCCACGAACGTAGTTAGCC‑ATTGATAGTTTTCATGTGCGACAATGGGCGCATCGAGGGAAAACTCGTCGAGCTGATTTTCATAATCGGCTGCCATCACGCGAACTTTa                                                                                                                                                                <  1:138861/250‑1 (MQ=255)                                                                                       gctgctGCAATACGGTTCCGACCGCGACTTCCAGTGAAGCGGAAGAACTTAACCCGGCACCCTGCGGCACATTGCCGCTGATCACCATGTCCACGCCGCCGAAGCTGTTGTTACGCAGTTGCAGATGTTTCACCACGCCACGAACGTAGTTAGCC‑ATTGATAGTTTTCATGTGCGACAATGGGCGCATCGAGGGAAAACTCGTCGAGCTGATTTTCATAATCGGCTGCCATCACGCGAACTTTACGGTCAt                                                                                                                                                         <  1:271291/251‑1 (MQ=255)                                                                                            tgcAATACGGTTCCGACCGCGACTTCCAGTGAAGCGGAAGAACTTAACACGGCACCCTGCGGCACATTGCCGCTGATCACCATGACCACGCCGCCGAAGCTGTTGTTACGCAGTTGCAGATGGTTAACCACGCCAAGAACGTAGTTAGCC‑ATTGATAGTTTTCATGTGCGACAATGGGCGCATCGAGGGAAAACTCGTCGAGCTGATTTTCATAATAGGCTGCCATCACGCGGACTTTACGGTCATCGCgt                                                                                                                                                    >  1:222052/1‑251 (MQ=255)                                                                                               aaTACGGTTCCGACCGCGACTTCCAGTGAAGCGGAAGAACTTAACCCGGCACCCTGCGGCACATTGCCGCTGATCACCATGTCCACGCCGCCGAAGCTGTTGTTACGCAGTTGCAGATGTTTCACCACGCCACGAACGTAGTTAGCC‑ATTGATAGTTTTCATGTGCGACAATGGGCGCATCGAGGGAAAACTCGTCGAGCTGATTTTCATAATCGGCTGCCATCACGCGAACTTTACGGTCATCGCgtggg                                                                                                                                                 >  1:8326/1‑250 (MQ=255)                                                                                               aaTACGGTTCCGACCGCGACTTCCAGTGAAGCGGAAGAACTTAACCCGGCAACCTGCGGCACATTGCCGCTGATCACCATGTCCACGCTGCCGAAGATGTTGTTACGCAGTTGCAGATGTTTCACCACGCCACGAACGTAGTTAGCC‑ATTGATAGTTTTCATGTGCGACAATGGGCGCATCGAGGGAAAACTCGTCGAGCTGATTTTCATAATCGGCTGCCATCACGCGAACTTTACGGTCATCGCgtggt                                                                                                                                                 >  1:462332/1‑251 (MQ=255)                                                                                                 tACGGTTCCGACCGCGACTTCCAGTGAAGCGGAAGAACTTAACCCGGCACCCTGCGGCACATTGCCGCTGATCACCATGTCCACGCCGCCGAAGCTGTTGTTACGCAGTTGCAGATGTTTCACCACGCCACGAACGTAGTTAGCC‑ATTGATAGTTTTCATGTGCGACAATGGGCGCATCGAGGGAAAACTCGTCGAGCTGATTTTCATAATCGGCTGCCATCACGCGAACTTTACGGTCATCGCGTGGTGc                                                                                                                                               >  1:210967/1‑251 (MQ=255)                                                                                                     gTTCCGACCGCGACTTCCAGTGAAGCGGAAGAACTTAACCCGGCACCCTGCGGCACATTGCCGCTGATCACCATGTCCACGCCGCCGAAGCTGTTGTTACGCAGTTGCAGATGTTTCACCACGCCACGAACGTAGTTAGCC‑ATTGATAGTTTTCATGTGCGACAATGGGCGCATCGAGGGAAAACTCGTCGAGCTGATTTTCATAATCGGCTGCCATCACGCGAACTTTACGGTCATCGCGTGGTGCACaa                                                                                                                                           <  2:6034/251‑1 (MQ=255)                                                                                                         cGACCGCGACTTCCAGTGAAGCGGAAGAACTTAACCCGGCACCCTGCGGCACATTGCCGCTGATCACCATGTCCACGCCGCCGAAGCTGTTGTTACGCAGTTGCAGATGTTTCACCACGCCACGAACGTAGTTAGCC‑ATTGATAGTTTTCATGTGCGACAATGGGCGCATCGAGGGAAAACTCGTCGAGCTGATTTTCATAATCGGCTGCCATCACGCGAACTTTACGGTCATCGCGTGGTGCACAACTGa                                                                                                                                       <  1:110591/251‑1 (MQ=255)                                                                                                            ccGCGACTTCCAGTGAAGCGGAAGAACTTAACCCGGCACCCTGCGGCACATTGCCGCTGATCACCATGTCCACGCCGCCGAAGCTGTTGTTACGCAGTTGCAGATGTTTCACCACGCCACGAACGTAGTTAGCC‑ATTGATAGTTTTCATGTGCGACAATGGGCGCATCGAGGGAAAACTCGTCGAGCTGATTTTCATAATCGGCTGCCATCACGCGAACTTTACGGTCATCGCGTGGTGCACAACTGATc                                                                                                                                     >  1:453581/1‑250 (MQ=255)                                                                                                                gACTTCCAGTGAAGCGGAAGAACTTAACCCGGCACCCTGCGGCACATTGCCGCTGATCACCATGTCCACGCCGCCGAAGCTGTTGTTACGCAGTTGCAGATGTTTCACCACGCCACGAACGTAGTTAGCC‑ATTGATAGTTTTCATGTGCGACAATGGGCGCATCGAGGGAAAACTCGTCGAGCTGATTTTCATAATCGGCTGCCATCACGCGAACTTTACGGTCATCGCGTGGTGCACAACTGATCACGGt                                                                                                                                >  1:428204/1‑251 (MQ=255)                                                                                                                    tCCAGTGAAGCGGAAGAACTTAACCCGGCACCCTGCGGCACATTGCCGCTGATCACCATGTCCACGCCGCCGAAGCTGTTGTTACGCAGTTGCAGATGTTTCACCACGCCACGAACGTAGTTAGCC‑ATTGATAGTTTTCATGTGCGACAATGGGCGCATCGAgg                                                                                                                                                                                                                   <  1:62452/164‑1 (MQ=255)                                                                                                                    tCCAGTGAAGCGGAAGAACTTAACCCGGCACCCTGCGGCACATTGCCGCTGATCACCATGTCCACGCCGCCGAAGCTGTTGTTACGCAGTTGCAGATGTTTCACCACGCCACGAACGTAGTTAGCC‑ATTGATAGTTTTCATGTGCGACAATGGGCGCATCGAgg                                                                                                                                                                                                                   >  2:62452/1‑164 (MQ=255)                                                                                                                     ccAGTGAAGCGGAAGAACTTAACCCGGCCCCCTGCGGCACAGTGCCGCTGATCACCATGTCCACGCCGCCGAAGCTGTTGTTACGCAGTTGCAGATGTTTCACCACGCCACGAACGTAGTTAGCC‑ATTGATAGTTTTCATGTGCGACAATGGGCGCATCGAGGGAAAACTCGTCGAGCTGATTTTCATAAGCGGCTGCCATCACGCGAACGTTACGGTCATCGCGTGGTGCACAACTGATCACGGTTTGa                                                                                                                            <  2:462332/250‑1 (MQ=255)                                                                                                                       aGTGAAGCGGAAGAACTTACCCCGGCACCCTGCGGCACATTGCCGCTGATCACCATGTCCACGCCGCCGAAGCTGTTGTTACGCAGTTGCAGATGTTTCACCACGCCACGAACGTAGTTAGCC‑ATTGATAGTTTTCATGTGCGACAATGGGCGCATCGAGGGAAAACTCGTCGAGCTGATTTTCATAATCGGCTGCCATCACGCGAACTTTACGGTCATCGCGTGGTGCACAACTGATCACGGTTTGATaa                                                                                                                         <  2:458791/251‑1 (MQ=255)                                                                                                                          gAAGCGGAAGAACTTAACCCGGCACCCTGCGGCACATTGCCGCTGATCACCATGTCCACGCCGCCGAAGCTGTTGTTACGCAGTTGCAGATGTTTCACCACGCCACGAACGTAGTTAGCC‑ATTGATAGTTTTCATGTGCGACAATGGGCGCATCGAGGGAAAACTCGTCGAGCTGATTTTCATAATCGGCTGCCATCACGCGAACTTTACGGTCATCGCGTGGTGCACAACTGATCACGGTTTGATaatc                                                                                                                       >  1:298531/1‑250 (MQ=255)                                                                                                                           aaGCGGAAGAACTTAACCCGGCCCCCTGCGGCACATTGCCGCTGATCACCATGTCCACGCCGCCGAAGCTGTTGTTACGCAGTTGCAGATGTTTCACCACGCCACGAACGTAGTTAGCC‑ATTGATAGTTTTCATGTGCGACAATGGGCGCATCGAGGGAAAACTCGTCGAGCTGATTTTCATAATCGGCTGCCATCACGCGAACTTTACGGTCATCGCGTGGTGCACAACTGATCACGGTTTGATaatcaa                                                                                                                     <  2:102207/251‑1 (MQ=255)                                                                                                                              cGGAAGAACTTAACCCGGCACCCTGCGGCACATTTCCGCTGATCACCATGTCCACGCCGCCGAAGCTGTTGTTACGCAGTTGCAGATGTTTCACCACGCCACGAACGTAGTTAGCC‑ATTGATAGTTTTCATGTGCGACAATGGGCGCATCGAGGGAAAACTCGTCGAGCTGATTTTCATAATCGGCTGCCATCACGCGAACTTTACGGTCATCGCGTGGTGCACAACTGATCACGGTTTGATaatcaatct                                                                                                                  <  2:394132/251‑2 (MQ=255)                                                                                                                                gaagaaCTGACCCCGGCACCCTGCGGCACATTGCCGCTGATCACCATGTCCACGCCGCCGAAGCTGTTGTTACGCAGTTGCAGATGTTTCACCACGCCACGAACGTAGTTAGCC‑ATTGATAGTTTTCATGTGCGACAATGGGCGCATCGAGGGAAAACTCGTCGAGCTGATTTTCATAATCGGCTGCCATCACGCGAACTTTACGGTCATCGCGTGGTGCACAACTGATCACGGTTTGATTATCAATcgcg                                                                                                                <  2:298531/251‑1 (MQ=255)                                                                                                                                                                                        gccgAAGCTGTTGTTACTCAGTTGCAGATGGTTCCCCACGCCACGAACGTAGTTAGCC‑ATTGATAGTTTTCATGGGCGACAATGGGCGCATCGAGGGAAAACTCGGCGAGCTGATTTTCATAATCGGCTGCGATCACGCGAACTTTACGGTCATCGCGTGGTGCACAACTGATCACGGTTTGATAATCAATCGTGCAGGGCAGAACGAAACCGTCGTTGGAGTGGGTGTGTGCACCAATCAAAATCAcgcg                                                        <  2:461076/251‑1 (MQ=255)                                                                                                                                                                                         ccgAAGCTGTTGTTACGAAGTTGCAGATGTTTCACCAAGCCACGAACGTAGTTAGCC‑ATTGATAGTTTTCATGTGCGACAATCGGCGCAACGAGGGAAAACTCGTCGAGCTGATTTTCATAATCGGCTGCCATCACGCGAACTTTACGGTCATCGCGTGGTGAACAACTGATCACGGTTTGATAATCAATCGCGAAGGGCAGAACGAAACCGTCGTTGTAGTCGGTGTGTTCACCACTCAAATTCACGCgg                                                       >  2:375656/1‑251 (MQ=255)                                                                                                                                                                                            aaGCTGTTGTTACGCAGTTGCAGATGTTTCCCCACGCCACGAACTTAGTTAGCC‑ATTGATAGTTTTCATGTTCGACAATGGGCGCATCGAGGGAAAACTCGTCGAGCTGATTTTCATAATCGGCTGCCATCACGCGAACTTTACGGTCATCGCGTGGTGCACAACTGATCACGGTTTGATAATCAATCGCGCCGGGCAGAACGAAACCGTCGTTGTAGTCGGTGTGTTCACCAATCAAATTCACTCGGCCa                                                    <  2:124435/251‑1 (MQ=255)                                                                                                                                                                                                       aCGCAGATGCAGATGTTTCACCACGCCACGAACGTAGTTAGCC‑ATTGATAGTTTTCATGTGCGACAATGGGCGCATCGAGGGAAAACTCGTCGAGCTGATTTTCATAATCGGCTGCCATCACGCGAACTTTACGGTCATCGCGTGGTGCACAACTGATCACGGTTTGATAATCAATCGCGCAGGGCAGAACGAAACCGTCGTTGTAGTCGGTGTGTTCACCAATCAAATTCACGCGGCCAGGCGCCTGAAt                                         >  2:116208/1‑251 (MQ=255)                                                                                                                                                                                                         gCAGTTGCAGATGTTTCACCACGCCACGAACGTAGTTAGCC‑ATTGATAGTTTTCATGTGCGACAATGGGCGCATCGAGGGAAAACTCGTCGAGCTGATTTTCATAATCGGCTGCCATCACGCGAACTTTACGGTCAGCGCGTGGTGCACAACTGATCACGGTTTGATAATCAATCGCGCAGGGCAGAACGAAACCGTCGTTGTAGTCGGTGTGTTCACCAATCAAATTCACGCGGCCAGGCGCCTGAATgg                                       <  1:372336/251‑1 (MQ=255)                                                                                                                                                                                                          cAGTTGCAGATGTTTCACCACGCCACGACCGTAGTTAGCC‑ATTGATAGTTTTCATGTGCGACAATGGGCGCAGCGAGGGAAAACGCGTCGAGCTGATTTTCATAATCGGCTGCCATCACGCGAACTTTACGGTCATCGCGTGGTGCACAACTGATCACGGTTGGATAATCAATCGCGCAGGGCAGAACGAAACCGTCGTTGTAGTCGGTGTGTTCACCAATCAAATTCACGCGGCCAGGCGCCTGTATGgt                                      <  1:145960/251‑1 (MQ=255)                                                                                                                                                                                                                          cccacgccacgAACGTAGTTAGCC‑ATTGATAGTTTTCATGTGCGACAATGGGCGCATCGAGGGAAAACTCGTCGAGCTGATTTTCATAATCGGCTGCCATCACGCGAACTTTACGGTCATCGCGTGGTGCACAACTGATCACGGTTTGATAATCAATCGCGCAGGGCAGAACGAAACCGTCGTTGTAGTCGGTGTGTTCACCAATCAAATTCACGCGGCCAGGCGCCTGAATGGTGTGAGTGGCAGGGTAg                      <  2:8326/250‑1 (MQ=255)                                                                                                                                                                                                                             acgccacgAACGTAGTTAGCC‑ATTGATAGTTTTCATGTGCGACAATGGGCGCATCGAGGGAAAACTCGTCGAGCTGATTTTCATAATCGGCTGCCATCACGCGAACTTTACGGTCATCGCGTGGTGCACAACTGATCACGGTTTGATAATCAATCGCGCAGGGCAGAACGAAACCGTCGTTGTAGTCGGTGTGTTCACCAATCAAATTCACGCGGCAAGGCGCCTGAATGGTGTGAGTGGAAGGGTAGcc                    >  2:102500/1‑250 (MQ=255)                                                                                                                                                                                                                             acgacacgAACGTAGTTAGCC‑ATTGATAGTTTTCATGGGCGACCATGGGCGCCTCGCGGGAAAACTCGTCGAGCTGATTTTCATAATCGGCTGCCATCACGCGAACTTTACGGTCATCGCGTGGTGCACAACTGATCACGGTTTGATAATCAATCGCGCAGGGCAGAACGAAACCGTCGTTGTAGTCGGTGTGTTCACCAATCAAATTCACGAGGCCAGGCGCCTGAATGGGGTGAGTGGCAGGGTAGcac                   >  2:257125/5‑249 (MQ=255)                                                                                                                                                                                                                                      aCGTAGTTAGCC‑ATTGATAGTTTTCATGTGCGACAATGGGCGCATCGAGGGAAAACTCGTCGAGCTGATTTTCATAATCGGCTGCCATCACGCGAACTTTACTGTCTTCGCGTGGTGCACAACTGATCACGGTTTGATAATCAATCGCGCAGGGCAGAACGAAACCGTCGTTGTAGTCGGTGTGTTCACCAATCAAATTCACGCGGCCAGGCGCCTGAATGGTGTGAGTGGCAGGGTAGCCAAATGCGttt          <  2:73591/251‑2 (MQ=255)                                                                                                                                                                                                                                        gTAGTTAGCC‑ATTGATAGTTTTCATGTGCGACAATGGGCGCATCGAGGGAAAACTCGTCGAGCTGATTTTCATAATCGGCTGCCATCACGCGAACTTTACGGTCATCGCGTGGTGCACAACTGATCACGGTTTGATAATCAATCGCGCAGGGCAGAACGAAACCGTCGTTGTAGTCGGTGTGTTCACCAATCAAATTCACGCGGCCAGGCGCCTGAATGGTGTGAGTGGCAGGGTAGCCAAATGCGTTGGc        >  1:189275/1‑251 (MQ=255)                                                                                                                                                                                                                                          aGTTAGCC‑ATTGATAGTTTTCATGTGCGACAATGGGCGCATCGAGGGAAAACTCGTCGAGCTGATTTTCATAATCGGCTGCCATCACGCGAACTTTACGGTCATCGCGTGGTGCACAACTGATCACGGTTTGATAATCAATCGCGCAGGGCAGAACGAAACCGTCGTTGTAGTCGGTGTGTTCACCAATCAAATTCACGCGGCCAGGCGCCTGAATGGTGTGAGTGGCAGGGTAGCCAAAGGCGTTGGCaa      <  2:453581/251‑1 (MQ=255)                                                                                                                                                                                                                                                  cATTGATAGTTTTCATGTGCGACAATGGGCGCATCGAGGGAAAACTCGTCGAGCTGATTTTCATAATCGGCTGCCATCACGCGAACTTTACGGTCATCGCGTGGTGCACAACTGATCACGGTTTGATAATCAATCGCGCAGGGCAGAACGAAACCGTCGTTGTAGTCGGTGTGTTCACCAATCAAATTCACGCGGCCAGGCGCCTGAATGGTGTGAGTGGCAGGGTAGCCAAATGCGTTGGCAAACag  >  1:334057/1‑248 (MQ=255)                                                                                                                                                                                                                                                  |                                                                                                                                                                                                                                                        GATCCCGCAGTTACAGCCTACAAACTGGTTTTCTGCTTCCTGACCGTTAAGCGCGATTTGTGCGCCGTCCAGCGGCAGATGATAAAGCTGCTGCAATACGGTTCCGACCGCGACTTCCAGTGAAGCGGAAGAACTTAACCCGGCACCCTGCGGCACATTGCCGCTGATCACCATGTCCACGCCGCCGAAGCTGTTGTTACGCAGTTGCAGATGTTTCACCACGCCACGAACGTAGTTAGCCCATTGATAGTTTTCATGTGCGACAATGGGCGCATCGAGGGAAAACTCGTCGAGCTGATTTTCATAATCGGCTGCCATCACGCGAACTTTACGGTCATCGCGTGGTGCACAACTGATCACGGTTTGATAATCAATCGCGCAGGGCAGAACGAAACCGTCGTTGTAGTCGGTGTGTTCACCAATCAAATTCACGCGGCCAGGCGCCTGAATGGTGTGAGTGGCAGGGTAGCCAAATGCGTTGGCAAACAG  >  NC_000913/789464‑789952

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