Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	NC_000913	789,580	C→A	E134* (GAA→TAA)	galK ←	galactokinase

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	789,580	0	C	A	97.8%	128.2 / NA	46	E134* (GAA→TAA)	galK	galactokinase
Reads supporting (aligned to +/- strand):  ref base C (1/0);  new base A (24/21);  total (25/21)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.85e-01

GTTTGAAGTTACTGTTGATGATGACGACAGCCACACCTTTGGGCATGGAAACTGCTTTGGTCCCCAGTGAGCGGCAATCGATCAGCAAGGCATGATCTTTCTTGC‑CGAGCGCGGAAATTAGCTGATCCATGATCCCGCAGTTACAGCCTACAAACTGGTTTTCTGCTTCCTGACCGTTAAGCGCGATTTGTGCGCCGTCCAGCGGCAGATGATAAAGCTGCTGCAATACGGTTCCGACCGCGACTTCCAGTGAAGCGGAAGAACTTAACCCGGCACCCTGCGGCACATTGCCGCTGATCACCATGTCCACGCCGCCGAAGCTGTTGTTACGCAGTTGCAGATGTTTCACCACGCC‑ACGAACGTAGTTAGCCCATTGATAGTTTTCATGTGCGACAATGGGCGCATCGAGGGAAAACTCGTCGAGCTGATTTTCATAATCGGCTGCCATCACGCGAACTTTACGGTCATCGCGTGGTGCACAACTGA  >  NC_000913/789334‑789819                                                                                                                                                                                                                                                        |                                                                                                                                                                                                                                                 ggtgtaatttACTGTTGATGTTGCCGCAAGCCCCACCTTGGGGCCTGGAAACTGCTTTGGTCCCCAGTGAGCGGCAACCGATCTGCAAGGCATGATCTGTCTTGC‑CGAGCGCTGAAATTAGTTGAACCAGGATACCGCAGTTACAGCCTACAAACTGGTTTTCTGCTTCCTGACCTTTAAGCGTGAGTTGTGCGCCGTCCAGCGGCAGATGATAAAGCTGCTGCAATACGGTTCCGACCGCGACTTACAgg                                                                                                                                                                                                                                              <  2:267614/243‑2 (MQ=255)     gaaataaCTGTTGATTATGACTACAGCCACACCTGTGGGCATGGAAACTGCTTTGGTCCCCAGTGGGCGGCAATCGATCAGCAAGGCTGGATCTTTCGTGC‑CGAGCGCGGAAATTAGCTGATCCATTATCCCGCAGTTACAGCCTACAAACTGGTTTTCTGCTTCCTGACCGTTGAGCGCGAGTTGTGCGCCGTCCAGCGGCAGATGATAAAGCTGCTGCAATACGGTTCCGACCGCGACTTACAGTGaa                                                                                                                                                                                                                                           <  2:350077/244‑1 (MQ=255)                  atgatgaCGACAGCCACACCTTTAGGCATGGAAACTGCTTTGGTCCCCAGTGAGCGGCAATCGATCAGCAAGGCATGATCTTTCTTGC‑CGAGCGCGGAAATTAGCTGATCCATGATCCCGCAGTTACAGCCTACAAACTGGTTTTCTGCTTCCTGACCGTTAAGCGCGATTTGTGCGCCGTCCAGCGGCAGATGATAAAGCTGCTGCAATACGGTTCCGACCGCGACTTACAGTGAAGCGGAAGAACTTa                                                                                                                                                                                                                              >  1:282658/1‑250 (MQ=255)                    gatgaCGACAGCCACACCTTTGGGCATGGAAACTGCTTTGGTCCCCAGTGAGCGGCAATCGATCAGCAAGGCATGATCTTTCTTGC‑CGAGCGCGGAAATTAGCTGATCCATGATCCCGCAGTTACAGCCTACAAACTGGTTTTCTGCTTCCTGACCGTTAAGCGCGATTTGTGCGCCGTCCAGCGGCAGATGATAAAGCTGCTGCAATACGGTTCCGACCGCGACTTACAGTGAAGCGGAAGAACTTAAc                                                                                                                                                                                                                            <  1:245494/250‑1 (MQ=255)                     atgaCGACAGCCACACCTTTGGGCATGGAAACTGCTTTGGTCCCCAGTGAGCGGCAATCGATCAGCAAGGCATGATCTTTCATGC‑CGAGCGCGGAAATTAGCTGATCCATGATCCCGCAGTTACAGCCTACAAACTGGTTTTCTGCTTCCTGACCGTTAAGCGCGATTTGTGCGCCGTCCAGCGGCAGATGATAAAGCTGCTGCAATACGGTTCCGACCGCGACTTACAGTGAAGCGGAAGAACTTAAccc                                                                                                                                                                                                                          >  1:384724/1‑251 (MQ=255)                      tgaCGACAGCCACACCTTTGGGCATGGAAACTGCTTTGGTCCCCAGTGAGCGGCAATCGATCAGCAAGGCATGATCTTTCTTGC‑CGAGCGCGGAAATTAGCTGATCCATGATCCCGCAGTTACAGCCTACAAACTGGTTTTCTGCTTCCTGACCGTTAAGCGCGATTTGTGCGCCGTCCAGCGGCAGATGATAAAGCTGCTGCAATGCGGTTCCGACCGCGACTTACAGTGAAGCGGAAGAACTTAACCCg                                                                                                                                                                                                                         >  1:255920/1‑251 (MQ=255)                                          ggCATGGAAACTGCTTTGGTCCCCAGTGAGCGGCAATCGATCAGCAAGGCATGATCTTTCTTGC‑CGAGCGCGGAAATTAGCTGATCCATGATCCCGCAGTTACAGCCTACAAACTGGTTTTCTGCTTCCTGACCGTTAAGCGCGATTTGTGCGCCGTCCAGCGGCAGATGATAAAGCTGCTGCAATACGGTTCCGACCGCGACTTACAGTGAAGCGGAAGAACTTAACCCGGCACCCTGCGGCACATTGcc                                                                                                                                                                                                     >  1:145621/1‑251 (MQ=255)                                                                ccAGTGAGCGGCAATCGATCAGCAAGGCATGATCTTTCTTGC‑CGAGCGCGGAAATTAGCTGATCCATGATCCCGCAGTTACAGCCTACAAACTGGTTTTCTGCTTCCTGACCGTTAAGCGCGATTTGTGCGCCGTCCAGCGGCAGATGATAAAGCTGCTGCAATACGGTTCCGACCGCGACTTACAGTGAAGCGGAAGAACTTAACCCGGCACCCTGCGGCACATTGCCGCTGATCACCATGTCCAcgccg                                                                                                                                                                               <  1:266641/251‑1 (MQ=255)                                                                        cGGCAATCGATCAGCAAGGCATGATCTTTCTTGCACG‑GCGCGGAAATTAGCTGATCCAAGATCCCGCAGTTACAGCCTACAACATGGTTTTCTGCTTCCTGACCGCTAAGCGCGCTTTGTGCGCCGTCCAGCGGCAGATGATAAAGCTGCTGCCATCCGGTTCCGACCGCACCTAACATTGAAGAGTAAGTACTTAACCCGGATCCCTGAAGCACATTGCAGCTGacccacatgtaaaacgagccgaagcg                                                                                                                                                                       >  2:169879/1‑226 (MQ=255)                                                                        cGGCAATCGATCAGCAAGGCATGATCTTTCTTGC‑CGAGCGCGGAAATTAGCTGATCCATGATCCCGCAGTTACAGCCTACAAACTGGTTTTCTGCTTCCTGACCGTTAAGCGCGATTTGTGCGCCGTCCAGCGGCAGATGATAAAGCTGCTGCAATACGGTTCCGACCGCGACTTACAGTGAAGCGGAAGAACTTAACCCGGCACCCTGCGGCACATTGCCGCTGATCACCATGTCCACGCCGCCGAAGCt                                                                                                                                                                       >  1:383326/1‑251 (MQ=255)                                                                                gctCAGCAAGGCATGATCTTTCTTGC‑CGAGCGCGGCCATTAGCTGATCCATGATCCCGCAGTTACAGCCTACAAACTGGATTTCTGCTTCCTGACCGTTAAGCGCGATTTGTGCGCCGTCCAGCGGCACATGCTAAAGCTGCTGCAATACGGTTCCGACCGCAACTTCCAGTGACGCGGAAGAACTTACACCGGTCCCCTGCGGCACATTTCCGCTGATCACCATGTCCACGCCGCCGAAGCAGTTGTTAc                                                                                                                                                               >  2:42407/3‑251 (MQ=255)                                                                                 aTCAGCAAGGCATGATCTTTCTTGC‑CGAGCGCGGAAATTAGCTGATCCATGATCCCGCAGTTACAGCCTACAAACTGGTTTTCTGCTTCCTGACCGTTAAGCGCGATTTGTGCGCCGTCCAGCGGCAGATGATAAAGCTGCTGCAATACGGTTCCGACCGCGACTTACAGTGAAGCGGAAGAACTTAACCCGGCACCCTGCGGCACATTGCCGCTGATCACCATGTCCACGCCGCCGAAGCTGTTGTTACg                                                                                                                                                              <  1:380475/251‑1 (MQ=255)                                                                                    aGCAAGGCATGATCTTTCTTGC‑CGAGCGCGGAAATTAGCTGATCCATGATCCCGCAGTTACAGCCTACAAACTGGTTTTCTGCTTCCTGACCGTTAAGCGCGCTTTGTGCGCCGTCCAGCGGCAGATGATAAAGCTGCTGCAATACGGTTCCGACCGCGACTTACAGTGAAGCGGAAGAACTTAACCCGGCACCCTGCGGCACATTGCCGCTGATCAACATGTCCACGCCGCCGAAGCTGTTGTTACGCa                                                                                                                                                            >  2:332907/1‑250 (MQ=255)                                                                                    aGCAAGGCATGATCTTTCTTGC‑CGAGCGCGGAAATTAGCTGATCCATGATCCCGCAGTTACAGCCTACAAACTGGTTTTCTGCTTCCTGACCGTTAAGCGCGATTTGTGCGCCGTCCAGCGGCAGATGATAAAGCTGCTGCAATACGGTTCCGACCGCGACTTACAGTGAAGCGGAAGAACTTAACCCGGCACCCTGCGGCACATTGCCGCTGATCACCATGTCCACGCCGCCGAAGCTGTTGTTACGCa                                                                                                                                                            <  1:332907/250‑1 (MQ=255)                                                                                                   ttcttGC‑CGAGCGCGGAAATTAGCTGATCCATGATCCCGCAGTTACAGCCTACAAACTGGTTTTCTGCTTCCTGACCGTTAAGCGCGATTTGTGCGCCGTCCAGCGGCAGATGATAAAGCTGCTGCAATACGGTTCCGACCGCGACTTACAGTGAAGCGGAAGAACTTAACCCGGCACCCTGCGGCACATTGCCGCTGATCACCATGTCCACGCCGCCGAAGCTGTTGTTACGCAGTTGCAGATGTTTcac                                                                                                                                            >  1:124248/1‑251 (MQ=255)                                                                                                           cGAGCGCGGAAATTAGCTGATCCATGATCCCGCAGTTACAGCCTACAAACTGGTTTTCTGCTTCCTGACCGTTAAGCGCGATTTGTGCGCCGTCCAGCGGCAGATGATAAAGCTGCTGCAATACGGTTCCGACCGCGACTTACAGTGAAGCGGAAGAACTTAACCCGGCACCCTGCGGCACATTGCCGCTGATCACCATGTCCACGCCGCCGAAGCTGTTGTTACGCAGTTGCAGATGTTTCAccacgcc                                                                                                                                      >  1:278063/1‑250 (MQ=255)                                                                                                                gcgGAAATTAGCTGATCCATGTTCCACCAGTTACAGCCTACAAACTGGGTTTCTGCTTCCTGACCGTTAAGCGCTCTTTGGGCGCCGTCCAGTGGCAGATGATAAAACTGCTGCATGACGGTTCCGACCGCGACTTACGGTGAAGCGGAAGAACTTATCCCGGCTCCCTGCGGCACATTGCCGCTGATCACCATGTCCACTCCGCCGAAGCTGTTGTTACGCAGTTGCAGATGTTGCACCACGCC‑ACGAAc                                                                                                                               <  2:148637/251‑1 (MQ=255)                                                                                                                   gAAATTAGCTGATCCATGATCACGCAGTTACAGCCTACAAACTGGATTTCTGCTTCCTGACAGTTAAGCGCGATTTGTGCGCCGTCCAGCGGCAGATGATAAAGCTGCTGCAATAAGGTTCCGACAGCGACTTACAGTGAAGCGGAAGAACTAAACCCGGCACCCTGCGGCACATTGCTGATGATCACCATGTCCACGCCGCCGAAGCAGTTGTTAAGCAGTTGATGATATTTCAcaacgccctgtacga                                                                                                                              >  2:310108/1‑242 (MQ=255)                                                                                                                      aTTAGCTGATCCATGATCCCGCAGTTACAGCCTACAAACTGGTTTTCTGCTTCCTGACCGTTAAGCGCGATTTGTGCGCCGTCCAGCGGCAGATGATAAAGCTGCTGCAATACGGTTCCGACCGCGACTTACAGTGAAGCGGAAGAACTTAACCCGGCACCCTGCGGCACATTGCCGCTGATCACCATGTCCACGCCGCCGAAGCTGTTGTTACGCAGTTGCAGATGTTTCACCACGCC‑ACGAACGTAGtt                                                                                                                         <  1:62771/251‑1 (MQ=255)                                                                                                                        tAGCTGATCCATGATCCCGCAGTTACAGCCTACAAACTGGTTTTCTGCTTCCTGACCGTTAAGCGCGATTTGTGCGCCGTCCAGCGGCAGATGATAAAGCTGCTGCAATACGGTTCCGACCGCGACTTACAGTGAAGCGGAAGAACTTAACCCGGCACCCTGCGGCACATTGCCGCTGATCACCATGTCCACGCCGCCGAAGCTGTTGTTACGCAGTTGCAGATGTTTCACCACGCC‑ACGAACGTAGtt                                                                                                                         >  1:238691/1‑249 (MQ=255)                                                                                                                         aGCTGATCCATGATCCCGCAGTTACAGCCTACAAACTGGTTTTCTGCTTCCTGACCGTTAAGCGCGATTTGTGCGCCGTCCAGCGGCAGATGATAAAGCTGCTGCAATACGGTTCCGACCGCGACTTACAGTGAAGCGGAAGAACTTAACCCGGCCCCCTGCGGCACATTGCCGCTGATCACCATGTCCACGCCGCCGAAGCTGTTGTTACGCAGTTGCAGATGTTTCACCACGCC‑ACGAACGACGTTAGc                                                                                                                      >  2:326731/1‑251 (MQ=255)                                                                                                                                    gATCCCGCAGTTACAGCCTACAAACTGGTTTTCTGCTTCCTGACCGTTAAGCGCGATTTGTGCGCCGTCCAGCGGCAGATGATAAAGCTGCTGCAATACGGTTCCGACCGCGACTTACAGTGAAGCGGAAGAACTTAACCCGGCACCCTGCGGCACATTGCCGCTGATCACCATGTCCACGCCGCCGAAGCTGTTGTTACGCAGTTGCAGATGTTTCACCACGCC‑ACGAACGTAGTTAGCCCATTGATAg                                                                                                            <  1:326731/250‑1 (MQ=255)                                                                                                                                       cccGCAGTTACAGCCTACAAACTGGTTTTCTGCTTCCTGACCGTTAAGCGCGATTTGTGCGCCGTCCAGCGGCAGATGATAAAGCTGCTGCAATACGGTTCCGACCGCGACTTACAGTGAAGCGGAAGAACTTAACCCGGCACCCTGCGGCACATTGCCGCTGATCACCATGTCCACGCCGCCGAAGCTGTTGTTACGCAGTTGCAGATGTTTCACCACGCC‑ACGAACGTAGTTAGCCCATTGATAGtttt                                                                                                        >  1:121613/1‑251 (MQ=255)                                                                                                                                        ccGCAGTTACAGCCTACAAACTGGTTTTCTGCTTCCTGACCGTTAAGCGCGATTTGTGCGCCGTCCAGCGGCAGATGATAAAGCTGCTGCAATACGGTTCCGACCGCGACTTACAGTGAAGCGGAAGAACTTAACCCGGCACCCTGCGGCACATTGCCGCTGATCACCATGTCCACGCCGCCGAAGCTGTTGTTACGCAGTTGCAGATGTTTCACCACGCC‑ACGAACGTAGTTAGCCCATTGATAGTTTTc                                                                                                       <  1:134550/251‑1 (MQ=255)                                                                                                                                                  aGCCTACAAACTGGTTTTCTGATTCCTGAACGTTAAGCGCGATTTGTGCGCCGTCCAGCGGCAGATGATAAAGCTGCTGCAATACGGTTCCGACCGCGACTTACAGTGAAGCGTAAGAACTTAACCCGGCACCCTGCGGCACATTGCCGCTGATCACCATGTCCACGCCGCCGAAGATGTTGTTACGCAGGTGCAAATGATAAAACACGCC‑ACGAACGTAATTAGCCCATTGATAGTCTTCATGTGCGAc                                                                                              >  2:93889/1‑250 (MQ=255)                                                                                                                                                        cAAACTGGTTTTCTGCTTCCTGACCGTTAAGCGCGATTTGTGCGCCGTCCAGCGACAGATGATAAAGCTGCTGCAATACGGTTCCGACCGCGACTTACAGTGAAGCGGAAGAACTTAACCCGGCACCCTGCGGCACATTGCCGCTGCTCACCATGTCCACGCCGCCGAAGCTGGTGTTACGCAGTTGCAGATGTTTCACCACGCC‑ACGACCGTAGTTAGCCCATTGATAGGTTTCATGTGCGACCATGGgc                                                                                       >  2:66703/1‑251 (MQ=255)                                                                                                                                                            cTGGTTTTCTGCTTCCTGACCGTTAAGCGCGATTTGTGCGCCGTCCAGAGGCAGATGATAAAGCTGCTGCAATACGGTTCCGACCGCGACTTACAGTGAAGCGGAAGAACTTAACCCGGCACCCTGCGGCACATTGCCGCTGATCACCATGTCCACGCCGCCGAAGCTGTTGTTACGCAGTTGCAGATGTTTCACCACGCC‑ACGAACGTAGTTAGCCCATTGATAGTTTTCATGTGCGACACTGGgcgaat                                                                                   >  1:221266/1‑248 (MQ=255)                                                                                                                                                             tGGTTTTCTGCTTCCTGACCGTTAAGCGCGATTTGTGCGCCGTCCAGCGGCAGATGATAAAACTGCTGCAATCCGGTTCCGACCGCGACTTACAGTGAAGCGGACGAACTTAACCCGGCACCCTGCGGCACATTGCCGCTGATCACCATGTCCACGCCGCCGAAGCTGTTGTTACGCAGTTGCAGATGTTTCACCACGCC‑ACGAACGTAGTTAGCCCATTGTTAGTTTTAATGTGCGACAATGGGCGCATc                                                                                  >  2:180928/1‑251 (MQ=255)                                                                                                                                                              ggTTTTCTGCTTCCTGCCCGTTAAGCGCGATTTGTGCGCCGTCCAGCGGCAGATGATAAAGCTGCTGCAATACGGTTCCGACCGCGACTTACAGTGAAGCGGAAGAACTTAACCCGGCACCCTGCGGCACATTGCCGCTGATCACCATGTCCACGCCGTCGAAGCTGTTGTTACTCAGTTGCAGATGTTTCACCACGCC‑ACGAACGTAGTTAGCCCATTGATAGTTTTCATGTGCGACAATGGGCGCAt                                                                                   <  1:205705/249‑1 (MQ=255)                                                                                                                                                              ggTTTTCTGCTTCCTGACCGTTAAGCGCGATTTGGGCGCCGTCCAGCGGCAGATGATAAAGCTGCTGCAATACGGTTCCGACCGCGACTTACAGTGAAGCGGAAGAACTTAACCCGGCCCCCTGCGGCACATTGCCGCTGATCACCATGTCCACGCCGCCGAAGCTGTTGTTACGCAGTTGCAGATGTTTCACCACGCC‑ACGAACGTAGTTAGCCCATTGATAGTTTTCATGTGCGACAATGGGCGCATCg                                                                                 <  2:211143/251‑1 (MQ=255)                                                                                                                                                               gTTTTCTGCTTACTGACCGTTAAGCGCGATTTGTGCGCCGTCAAGCGGCAGATGCTAAAGCTGCTGCAATACGGTTCCGACCGCGACTTACAGTGAAGCGGAAGAACTTAACCCGGCACCCTGCGGCACATTGCCGCTGATCACCATGTCCACGCCGCCGCAGCTGTTGTTACGCAGTTGCATCTGTTTCACCAAGCCCAC‑AACGTAGTTAGCCCATTGGTAGTTTTCATGTGCGCCAATGGGAGCagcgg                                                                                >  2:230710/1‑247 (MQ=255)                                                                                                                                                                 tttCTGCTTCCTGACCGTTAAGCGCGATTTGTGCGCCGTCCAGCGGCAGATGATAAAGCTGCTGCAATACGGTTCCGACCGCGACTTACAGTGAAGCGGAAGAACTTAACCCGGCACCCTGCGGCACATTGCCGCTGATCACCATGTCCACGCCGCCGAAGCTGTTGTTACGCAGTTGCAGATGTTTCACCACGCC‑ACGAACGTAGTTAGCCCATTGATAGTTTTCATGTGCGACAATGGGCGCATCGAg                                                                               <  1:161543/250‑1 (MQ=255)                                                                                                                                                                    cTGCTTCCTGACCGTTAAGCGCGATTTGTGCGCCGTCCAGCGGCAGATGATAAAGCTGCTGCAATACGGTTCCGACCGCGACTTACAGTGAAGCGGAAGAACTTAACCCGGCACCCTGCGGCACATTGCCGCTGATCACCATGTCCACGCCGCCGAAGCTGTTGTTACGCAGTTGCAGATGTTTCACCACGCC‑ACGAACGTAGTTAGCCCATTGATAGTTTTCATGTGCGACAATGGGCGCATCGAGGGa                                                                            >  1:262481/1‑250 (MQ=255)                                                                                                                                                                            tGACCGTTAAGCGCGATTTGTGCGCCGTCCAGCGGCAGATGATAAAGCTGCTGCAATACGGTTCCGACCGCGACTTACAGTGAAGCGGAAGAACTTAACCCGGCACCCTGCGGCACATTGCCGCTGATCACCATGTCCACGCCGCCGAAGCTGTTGTTACGCAGATGCAGATGTTTCACCACGCC‑ACGAACGTAGTTAGCCCATTGATAGTTTTCATGTGCGACAATGGGCGCATCTAGGGAAAACtcgtc                                                                   >  1:138794/1‑251 (MQ=255)                                                                                                                                                                                                       cccAGCGGCAGATGATAAAGCTGCTGCAATACGGTTCCGACCGCGACTTACAGTGAAGCGGAAGAACTTAACCCGGCACCCTGCGGCACATTGCCGCTGATCACCATGTCCACGCCGCCGAAGCTGTTGTTACGCAGTTGCAGATGTTTCACCACGCC‑ACGAACGTAGTTAGCCCATTGATAGTTTTCATGTGCGACAATGGGCGCATCGAGGGAAAACTCGTCGAGCTGATTTTCATAATCGGCTGCCAt                                        <  1:42407/250‑1 (MQ=255)                                                                                                                                                                                                         cAGCGGCAGATGATAAAGCTGCTGCAATACGGTTCCGACCGCGACTTACAGTGAAGCGGAAGAACTTAACCCGGCACCCTGCGGCACATTGCCGCTGATCACCATGTCCACGCCGCCGAAGCTGTTGTTACGCAGTTGCAGATGTTTCACCACGCC‑ACGAACGTAGTTAGCCCATTGATAGTTTTCATGTGCGACAATTGGCGCATCGAGGGAAAACTCGTCGAACTGATTTTCATAATCGGCCGCCATcc                                      >  2:142420/1‑250 (MQ=255)                                                                                                                                                                                                           gcggcatagtatAAAGCTGCTGCAATACGGTTCCGACCGCGACTTACAGGGAAGCGGAAGAACTTAACCCGGCACCCTGCGGCACATTGCCGCTGATCACCATGTCCACGCCGCCGATGCTGTTGTTACGCAGTTTCATATGTTTATCCACGCC‑ACGAACTTAGTTAGCCCATTGATAGTTTTCATGTGCGACAATGGGCGCATCGAGGGAAAACTCGTCGAGCTGATTTTCATCATCGGCTGCTATCAcg                                    <  2:43017/241‑1 (MQ=255)                                                                                                                                                                                                           gCGGCAGAGGATAAAGCTGCTGCAAAACGGTTCCGACCGCGACTTACATTGAAGCGGAAGAACTTACCCCTGCACCCTGCGGCACATTGCCGCTGATCACCATGTCCCCGCCGCCGAAGCTGTTGTTACGCAGTTGCAGATGTGTCACCACGCC‑ACGAACGTAGTTAGCCCATTGATAGTTTTCATGTGCGACAATGGGCGCATCGAGGGAAAACTCGTCGAGCTGATTTTCATAATCGGCTGCCATCAcg                                    <  2:191025/251‑1 (MQ=255)                                                                                                                                                                                                           ccGGCTGATGATAAAGCTCCTGCAATACGTTTCCGTCCGCGACTTACAGTGAAGCGGAATAATTTAACCCGGCACCCTGCGGCCCATTGCCGCTGATCACCATGTCCACGCCGCCGAAGCTGTTGTTACGCAGTTGCAGATGTTTCACCACGCC‑ACGAACGTAGTTAGCCCATTGATAGTTTTCATGTGCGACAATGGGCGCATCTAGGGAAAACGCTTCGAGCTGATTTTCATAATCGGCTGCCATCAcg                                    <  1:135786/250‑1 (MQ=255)                                                                                                                                                                                                             ggCAGATGATAAAGCTGCTGCAATACGGTTCCGACCGCGACTTACAGTGAAGCGGAAGAACTTAACCCGGCACCCTGCGGCACATTGCCGCTGATCACCATGTCCACGCCGCCGAAGCTGTTGTTACGCAGTTGCAGATGTTTCACCACGCC‑ACGAACGTAGTTAGCCCATTGATAGTTTTCATGTGCGACAATGGGCGCATCGAGGGAAAACTCGTCGAGCTGATTTTCATAATCGGCTGCCATCAcgcg                                  >  1:195138/1‑251 (MQ=255)                                                                                                                                                                                                                    gatAAAGCTGCTGCAATACGGTTCCGACCGCGACTTACAGTGAAGCGGAAGAACTTAACCCGGCACCCTGCGGCACATTGCCGCTGATCACCATGTCCACGCCGCCGAAGCTGTTGTTACGCAGTTGCAGATGTTTCACCACGCC‑ACGAACGTAGTTAGCCCATTGATAGTTTTCATGTGCGACAATGGGCGCATCGAGGGAAAACTCGTCGAGCTGATTTTCATAATCGGCTGCCATCACGCGAACTTTa                           <  1:222859/251‑1 (MQ=255)                                                                                                                                                                                                                      tAAAGCTGCTGCAATACGGTTCCGACCGCGACTTACAGTGAAGCGGAAGAACTTAACCCGGCACCCTGCGGCACATTGCCGCTGATCACCATGTCCACGCCGCCGAAGCTGTTGTTACGCAGTTGCAGATGTTTCACCACGCC‑ACGAACGTAGTTAGCCCATTGATAGTTTTCATGTGCGACAATGGGCGCATCGAGGGAAAACTCGTCGAGCTGATTTTCATAATCGGCTGCCATCACGCGAACTTTACg                         >  1:86075/1‑251 (MQ=255)                                                                                                                                                                                                                                gcAATACGGTTCCGACCGCGACTTACATTGAAGCGGAAGAACTTAACCCGGCACCCTGCGGCACATTGCCGCTGATCACCATGCCCACGCCGCCGAAGCTGTTGTTACGCAGTTGCAGATGTTTCACCACGCC‑ACGAACGTAGTTAGCCCATTGATAGTTTTCATGTGCGACAATGGGCGCATCGAGGGAAAACTCGTCGAGCTGATTTTCATAATCGGCTGCCATCACGCGAACTTTACGGTCATcgcg                <  1:42710/250‑1 (MQ=255)                                                                                                                                                                                                                                     tcGGTTCCGAACGCGACTTAAAGTGAAGCGGAAGAACTTAACCCGGCCCCCTGCGGCACATTGCCGCTGATCCCCATGTCCACGCCGCCGAAGCTGTTGTTACGCAGTTGCAGATGTTTCACCACGCC‑ACGAACGTAGTTAGCCCATTGATAGTTTTCATGTGCGACAATGGGCGCATCGAGGGAAAACTCGTCGAGCTGATTTTCATAATCGGCTGCCATCACGCGAACTTTACGGTCATCGCgtggtg           <  1:140296/249‑1 (MQ=255)                                                                                                                                                                                                                                       ggTTCCGACCGAGACTTACAGTGAAGAGGAAGAACTTCACACGGCACCCTGCGGCACATTGCCGATGATCACCATGTCCACGCCGCCGCAGCTGTTGTTACGCAGTTGCAGATGTTTCAACACGCC‑ACGACCGTAGTTAGCCCATTGATAGTTTTCATGTGCGACAATGGGCGCATCGAGGGAAAACTCGTCGCGCTGGTTTGCATAATCGGCTGACATCACGCGAACATTACGGTCATCGCGTGGGGca         >  2:16737/1‑250 (MQ=255)                                                                                                                                                                                                                                             caccccgACTTACAATGAAGCGTAAGAACTTAACCCGGCACCCTGCGGCACATTGCCGCTGATCACCATGTCCACGCCGCCGAAGCTGTTGTTACGCAGTTGCAGATGTTTCACCACGCC‑ACGAACGTAGTGAGCCCATTGATAGTTTTCATGTGCGACAATGGGCGCATCGAGGGAAAACTCGTCGAGCTGATTTTCATAATCGGCTTCCATCACGCGAACTTTACGGTCATCGCGTGGTGCACAACTGa  <  1:152108/246‑1 (MQ=255)                                                                                                                                                                                                                                                        |                                                                                                                                                                                                                                                 GTTTGAAGTTACTGTTGATGATGACGACAGCCACACCTTTGGGCATGGAAACTGCTTTGGTCCCCAGTGAGCGGCAATCGATCAGCAAGGCATGATCTTTCTTGC‑CGAGCGCGGAAATTAGCTGATCCATGATCCCGCAGTTACAGCCTACAAACTGGTTTTCTGCTTCCTGACCGTTAAGCGCGATTTGTGCGCCGTCCAGCGGCAGATGATAAAGCTGCTGCAATACGGTTCCGACCGCGACTTCCAGTGAAGCGGAAGAACTTAACCCGGCACCCTGCGGCACATTGCCGCTGATCACCATGTCCACGCCGCCGAAGCTGTTGTTACGCAGTTGCAGATGTTTCACCACGCC‑ACGAACGTAGTTAGCCCATTGATAGTTTTCATGTGCGACAATGGGCGCATCGAGGGAAAACTCGTCGAGCTGATTTTCATAATCGGCTGCCATCACGCGAACTTTACGGTCATCGCGTGGTGCACAACTGA  >  NC_000913/789334‑789819

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