Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	NC_000913	4,286,001	(C)5→6	coding (1371/1959 nt)	acs ←	acetyl‑CoA synthetase (AMP‑forming)

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	4,286,001	1	.	C	91.4%	110.0 / 3.8	35	coding (1371/1959 nt)	acs	acetyl‑CoA synthetase (AMP‑forming)
Reads supporting (aligned to +/- strand):  ref base . (1/2);  new base C (11/21);  total (12/23)
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.71e-01

CGGATGCGCCACCAGCGCCGACTCAATCTCTGCCGTCCCCAGACGGTGACCGGAGACGTTCAGCACGTCGTCCACACGCCCGGTTATCCAGTAATAGCCATCTTCATCGCGACGCGCGCCGTCGCCGCTGAAATACATATTTTTGAAGGTGGAGAAGTAGGTCTGTTCAAAACGTTCGTGATCGCCAAACAGCGTACGCGCCTGACCCGGCCAGGAGTCGGTGATTACCAGGCTACCTTCGGTGGCCCCC‑TCCAGCGGGTTACC‑TTCGTTATCGACCAGCGCCGGTTGCACGCCGAAGAACGGACGTGTTGCCGAACCGGCTTTCAGCTCGGTAGCGCCAGGCAGCGGGGTGATCATGAAACCGCCGGTTTCGGTCTGCCACCAGGTATCGACCACCGGACATTTCTCGTTGCCGATTTTTTTCCAGTACCACTCCCACGCTTCCGGGTTAATTGGCTCGCCCACGGAACCGAGAATGCGCAGCGACGA  >  NC_000913/4285752‑4286240                                                                                                                                                                                                                                                           |                                                                                                                                                                                                                                                 cGGATGCGCCACCAGCGCCGACTCAATCTCTGCCGTCCCCAGACGGTGACCGGAGACGTTCAGCACGTCGTCCACACGCCCGGTTATCCAGTAATAGCCATCTTCATCGCGACGCGCGCCGTCGCCGCTGAAATACATATTTTTGAAGGTGGAGAAGTAGGTCTGTTCAAAACGTTCGTGATCGCCAAACAGCGTACGCGCCTGACCCGGCCAGGAGTCGGTGATTACCAGGCTACCTTCGGTGGcccccc                                                                                                                                                                                                                                                  <  2:15608/251‑2 (MQ=255)   gATGCGCCACCAGCGCCGACTCAATCTCTGCCGTCCCCAGACGGTGACCGGAGACGTTCAGCACGTCGTCCACACGCCCGGTTATCCAGTAATAGCCATCTTCATCGCGACGCGCGCCGTCGCCGCTGAAATACATATTTTTGAAGGTGGAGAAGTAGGTCTGTTCAAAACGTTCGTGATCGCCAAACAGCGTACGCGCCTGACCCGGCCAGGAGTCGGTGATTACCAGGCTACCTTCGGTGGcccccctc                                                                                                                                                                                                                                                >  1:81987/1‑248 (MQ=255)          caccaGCGCCGACTCAATCTCTGCCGTCCCCAGACGGTGACCGGAGACGTTCAGCACGTCGTCCACACGCCCGGTTATCCAGTAATAGCCATCTTCATCGCGACGCGCGCCGTCGCCGCTGAAATACATATTTTTGAAGGTGGAGAAGTAGGTCTGTTCAAAACGTTCGTGATCGCCAAACAGCGTACGCGCCTGACCCGGCCAGGAGTCGGTGATTACCAGGCTACCGTCGGTGGCCCCCCTCCAGCggg                                                                                                                                                                                                                                         <  1:234773/251‑1 (MQ=255)           accaGCGCCGACTCAATCTCTGCCGGCCCCAGACGGTGACCGGAGACGTTCAGCACGTCGTCCACACGCCCGGTTATCCAGTAATAGCCATCTTCATCGCGACGCGCGCCGTCGCCGCTGAAATACATATTTTTGAAGGTGGAGAAGTAGGTCTGTTCAAAACGTTCGTGATCGCCAAACAGCGTACGCGCCTGACCCGGCCAGGAGTCGGTGATTACCAGGCTACCTTCGGTGGCCCCC‑TCCAGCGGGtt                                                                                                                                                                                                                                       <  1:47886/251‑1 (MQ=255)                              ctGCCGTCCCCGGACGGTGACCGGAGACGTTCAGCACGTCGTCCACACGCCCTGTTATCCGGTAATAGCCATCTTCATCGCGACGCGCGCCGTCGCGTCTGAAATGCATATGTTTTAAGTTGGAGAAGTAGGTCTGTTCAAAACGGTCGGGATCGCCAAACAGCGTACGCGCTTGACCCGTCCAGTAGGCGGTGATTACCATGCTACCTTCGGTGGCCCCCCTCCAGCGGGTTACG‑GTCGTTATCGAcc                                                                                                                                                                                                                      <  2:133476/249‑1 (MQ=255)                                                 aCCGGAGACGTTCAGCACGTCGTCCACACGCCCGGTTATCCAGTAATAGCCATCTTCATCGCGACGCGCGCCGTCGCCGCTGAAATACATATTTTTGAAGGTGGAGAAGTAGGTCTGTTCAAAACGTTCGTGATCGCCAAACAGCGTACGCGCCTGACCCGGCCAGGAGTCGGTGATTACCAGGCTACCTTCGGTGGCCCCCCTCCAGCGGGTTACC‑TTCGTTATCGACCAGCGCCGGTTGCACGCCgaag                                                                                                                                                                                                 <  1:100783/251‑1 (MQ=255)                                                 aCCGGAGACGTTCAGCACGTCGTCCACACGCCCGGTTATCCAGTAATAGCCATCTTCATCGCGACGCGCGCCGTCGCCGCTGAAATACATATTTTTGAAGGTGGAGAAGTAGGTCTGTTCAAAACGTTCGTGATCGCCAAACAGCGTACGCGCCTGACCCGGCCAGGAGTCGGTGATTACCAGGCTACCTTCGGTGGCCCCCCTCCAGCGGGTTACC‑TTCGTTATCGACCAGCGCCGGTTGCACGCCgaag                                                                                                                                                                                                 <  1:164282/251‑1 (MQ=255)                                                    ggAGACGTTCAGCACGTCGTCCACACGCCCGGTTATCCAGTAATAGCCATCTTCATCGCGACGCGCGCCGTCGCCGCTGAAATACATATTTTTGAAGGTGGAGAAGTAGGTCTGTTCAAAACGTTCGTGATCGCCAAACAGCGTACGCGCCTGACCCGGCCAGGAGTCGGTGATTACCAGGCTACCTTCGGTGGCCCCCCTCCAGCGGGTTACC‑TTCGTTATCGACCAGCGCCGGTTGCACGCCGAAGAAc                                                                                                                                                                                              <  1:52551/251‑1 (MQ=255)                                                       gaCGTTCAGCACGTCGTCCACACGCCCGGTTATCCAGTAATAGCCATCTTCATCGCGACGCGCGCCGTCGCCGCTGAAATACATATTTTTGAAGGTGGAGAAGTAGGTCTGTTCAAAACGTTCGTGATCGCCAAACAGCGTACGCGCCTGACCCGGCCAGGAGTCGGTGATTACCAGGCTACCTTCGGTGGCCCCCCTCCAGCGGGTTACC‑TTCGTTATCGACCAGCGCCGGTTGCACGCCGAAGAACgg                                                                                                                                                                                            >  1:244646/1‑250 (MQ=255)                                                            tCAGCACGTCGTCCACACGCCCGGTTATCCAGTAATAGCCATCTTCATCGCGACGCGCGCCGTCGCCGCTGAAATACATATTTTTGAAGGTGGAGAAGTAGGTCTGTTCAAAACGTTCGTGATCGCCAAACAGCGTACGCGCCTGACCCGGCCAGGAGTCGGTGATTACCAGGCTACCTTCGGTGGCCCCCCTCCAGCGGGTTACC‑TTCGTTATCGACCAGCGCCGGTTGCACGCCGAAGAACGGACgtg                                                                                                                                                                                       <  1:112787/250‑1 (MQ=255)                                                                 aCGTCGTCCACACGCCCGGTTATCCAGTAATAGCCATCTTCATCGCGACGCGCGCCGTCGCCGCTGAAATACATATTTTTGAAGGTGGAGAAGTAGGTCTGTTCAAAACGTTCGTGATCGCCAAACAGCGTACGCGCCTGACCCGGCCAGGAGTCGGTGATTACCAGGCTACCTTCGGTGGCCCCCCTCCAGCGGGTTACC‑TTCGTTATCGACCAGCGCCGGTTGCActccttattccggacgtgttgccg                                                                                                                                                                                 <  1:55261/251‑24 (MQ=255)                                                                                 cGGTTATCCAGTAATAGCCATCTTCATCGCGACGCGCGCCGTCGCCGCTGAAATACATATTTTTGAAGGTGGAGAAGTAGGTCTGTTCAAAACGTTCGTGATCGCCAAACAGCGTACGCGCCTGACCCGGCCAGGAGTCGGTGATTACCAGGCTACCTTCGGTGGCCCCCCTCCAGCGGGTTACC‑TTCGTTATCGACCAGCGCCGGTTGCACGCCGAAGAACGGACGTGTTCCCGAACCGGCTTTCAGCTc                                                                                                                                                                 >  2:141558/1‑251 (MQ=255)                                                                                           gTAATAGCCAGCTTCATCGCGACGCGCGCCGTCGCCGCTGAAATACATATTGTTGAAGGTGGAGAAGTAGGTCTGTTCAAAACGTTCGTGATCGCCAAACAGCGTACGCGCCTGACCCGGCCAGGAGTCGGTGATTACCAGGCTACCTTCGGTGGCCCCCCTCCAGCGGGTTACC‑TTCGTTATCGACCAGCGCCGGTTGCACGCCGAAGAACGGACGTGTTGCCGAACCGGCTTTCAGCTCGGTAGCGCCa                                                                                                                                                       <  2:249478/251‑1 (MQ=255)                                                                                             aaTAGCCATCTTCATCGCGACGCGCGCCGTCGCCGCTGAAATACATATTTTTGAAGGTGGAGAAGTAGGTCTGTTCAAAACGTTCGTGATCGCCAAACAGCGTACGCGCCTGACCCGGCCAGGAGTCGGTGATTACCAGGCTACCTTCGGTGGCCCCCCTCCAGCGGGTTACC‑TTCGTTATCGACCAGCGCCGGTTGCACGCCGAAGAACGGACGTGTTGCCGAACCGGCTTTCAGCTCGGTAGCGCCAg                                                                                                                                                      <  1:132805/250‑1 (MQ=255)                                                                                                      cTTCATCGCGACGCGCGCCGTCGCCGCTGCAATACATATTTTTGAAGGTGGAGAAGTAGGTCTGTTCAAAACGTTCGTGATCGCCAAACAGCGTACGCGCCTGACCCGGCCAGGAGTCGGTGATTACCAGGCTACCTTCGGTGGCCCCCCTCCAGCGGGTTACC‑TTCGTTATCGACCAGCGCCGGTTGCACGCCGAAGAACGGACGTGTTGCCGAACCGGCTTTCAGCTCGGTAGCGCCAGGCAGCggggg                                                                                                                                            >  1:211731/1‑250 (MQ=255)                                                                                                        tCATCGCGACGCGCGCCGTCGCCGCTGAAATACATATTTTTGAAGGTGGAGAAGTAGGTCTGTTCAAAACGTTCGTGATCGCCAAACAGCGTACGCGCCTGACCCGGCCAGGAGTCGGTGATTGCCAGGCTACCTTCGGTGGCCCCCCTCCAGCGGGTTACC‑TTCGTTATCGACCAGCGCCGGTTGCACGCCGAAGAACGGACGTGTTGCCGAACCGGCTTTCAGCTCGTTAGCGCCAGGCAGCggg                                                                                                                                              <  1:90369/247‑1 (MQ=255)                                                                                                        tCATCGCGACGCGCGCCGTCGCCGCTGAAATACATATTTTTGAAGGTGGAGAAGTAGGTCTGTTCAAAACGTTCGTGATCGCCAAACAGCGTACGCGCCTGACCCGGCCAGGAGTCGGTGATTACCAGGCTACCTTCGGTGGCCCCCCTCCAGCGGGTTACC‑TTCGTTATCGACCAGCGCCGGTTGCACGCCGAAGAACGGACGTGTTGCCGAACCGGCTTTCAGCTCGGTAGCGCCAGGCAGCgggg                                                                                                                                             <  1:144521/248‑1 (MQ=255)                                                                                                          aTCGCGACGCGCGCCGTCGCCGCTGAAATACATATTTTTGAAGGTGGAGAAGTAGGTCTGTTCAAAACGTTCGTGATCGCCAAACAGCGTACGCGCCTGACCCGGCCAGGAGTCGGTGATTACCAGGCTACCTTCGGTGGCCCCCCTCCAGCGGGTTACC‑TTCGTTATCGACCAGCGCCGGTTGCACGCCGAAGAACGGACGTGTTGCCGAACCGGCTTTCAGCTCGGTAGCGCCAGGCAGCGGGGTGATc                                                                                                                                        >  2:153349/1‑251 (MQ=255)                                                                                                          aTCGCGACGCGCGCCGTCGCCGCTGAAATACATATTTTTGAAGGTGGAGAAGTAGGTCTGTTCAAAACGTTCGTGATCGCCAAACAGCGTACGCGCCTGACCCGGCCAGGAGTCGGTGATTACCAGGCTACCTTCGGTGGCCCCCCTCCAGCGGGTTACC‑TTCGTTATCGACCAGCGCCGGTTGCACGCCGAAGAACGGACGTGTTGCCGAACCGGCTTTCAGCTCGGTAGCGCCAGGCAGCGGGGTGATc                                                                                                                                        <  1:195328/251‑1 (MQ=255)                                                                                                             gcgACGCGCGCCGTCGCCGCTGAAATACATATTTTTGAAGGTGGAGAAGTAGGTCTGTTCAAAACGTTCGTGATCGCCAAACAGCGTACGCGCCTGACCCGGCCAGGAGTCGGTGATTACCAGGCTACCTTCGGTGGCCCCCCTCCAGCGGGTTACC‑TTCGTTATCGACCAGCGCCGGTTGCACGCCGAAGAACGGACGTGTTGCCGAACCGGCTTTCAGCTCGGTAGCGCCAGGCAGCGGGGTGATCATg                                                                                                                                     <  2:212588/251‑1 (MQ=255)                                                                                                               gACGCGCGCCGTCGCCGCTGAAATACATATTTTTGAAGGTGGAGAAGTAGGTCTGTTCAAAACGTTCGTGATCGCCAAACAGCGTACGCGCCTGACCCGGCCAGGAGTCGGTGATTACCAGGCTACCTTCGGTGGCCCCCCTCCAGCGGGTTACC‑TTCGTTATCGACCAGCGCCGGTTGCACGCCGAAGAACGGACGTGTTGCCGAACCGGCTTTCAGCTCGGTAGCGCCAGGCAGCGGGGTGATCATGaa                                                                                                                                   <  2:81987/251‑1 (MQ=255)                                                                                                                      gcCGTCGCCGCTGAAATACATATTTTTGAAGGTGGAGAAGTAGGTCTGTTCAAACCGTTCGTGATCGCCAAACAGCGTACGCGCCTGACCCGGCCAGGAGTCGGTGATTACCAGGCTACCTTCGGTGGCCCCCCTCCAGCGGGTTACC‑TTCGTTATCGACCAGCGCCGGTTGCACGCCGAAGAACGGACGTGTTGCCGAACCGGCTTTCAGCTCGGTAGCGCCAGGCAGCGGGGTGATCATGAAAccgccg                                                                                                                            >  1:54144/1‑251 (MQ=255)                                                                                                                                cTGAAATACATATTTTTGAAGGTGGAGAAGTAGGTCTGTTCAAAACGTTCGTGATCGCCAAACAGCGTACGCGCCTGACCCGGCCAGGAGTCGGTGATTACCAGGCTACCTTCGGTGGCCCCCCTCCAGCGGGTTACC‑TTCGTTATCGACCAGCGCCGGTTGCACGCCGAAGAACGGACGTGTTGCCGAACCGGCTTTCAGCTCGGCAGCGCCAGGCAGCGGGGTGATCATGACACCGCCGGTTTCGGTCt                                                                                                                  >  2:179248/1‑251 (MQ=255)                                                                                                                                        cATATTTTTGAAGGTGGAGAAGTAGGTCTGTTCAAAACGTTCGTGCTCGCCAAACAGCGTACGCGCCTGACCCGGCCAGGAGTCGGTGATTACCAGGCTACCTTCGGTGGCCCCC‑TCCAGCGGGTTACC‑TTCGTTATCGACCAGCGCCGGTTGCACGCCGAAGAACGGACGTGTTGCCGAACCGGCTTTCAGCTCGGTAGCGCCAGGCAGCGGGGTGATCATGAAACCGCCGGTTTCGGTCTGCCACCAgg                                                                                                         >  1:37880/1‑251 (MQ=255)                                                                                                                                               ttGAAGGGGGAGAAGTAGGTCTTTTCAAAACGTTCGTAATCCCCAAACAGCGTCCGCGCCTGACCCGGCCAGGAGTCGGTGATTACCAGGCTCCTTTCGGTGGCCCCCCTCCAGCGGGTTACC‑TTCGTTATCGACCAGCGCCGGTTGCACGCCGAAGAACGGACGTGTTGCCGACCCGGCTTTCAGCTCGGTAGCGCCAGTCCGCGGGGTGATCATGAAACCGCCGGTTTCGGTCTGCCACCAGGTATCg                                                                                                    <  1:237866/250‑1 (MQ=255)                                                                                                                                                      tGGAGAAGTAGGTCTGTTCAAAACGTTCGTGATCGCCAAACAGCGTACGCGCCTGACCCGGCCAGGAGTCGGTGATTACCAGGCCACCTTCGGTGGCCCCCCTCCAGCGGGTTACC‑TTCGTTATCGACCAGCGCCGGTTGCACGCCGAAGAACGGACGTGTTGCCGAACCGGCTTTCAGCTCGGTAGCGCCAGGCAGCGGGGTGATCATGAAACAGCCGGTTTCGGTCTGCCACCAGGTATCGACCACCgg                                                                                            >  1:119255/1‑251 (MQ=255)                                                                                                                                                                     gTTCAAACCGTTCGTGATCGCCAAACAGCGTACGCGCCTGCCCCGGCCAGGAGTCGGTGATTACCAGGCTCCCTCGGGTGGCCCCCCTCCTGCGGTTTACC‑TTCGTTATCGACCGGCGCGGTTTGCACGCCGAAGAACGGACGTGTTGCCGAACGGGTTTTCAGCTCGGTAGCTCCAGGGAGCGGGGTGATCATAAAACCGCCGGTTTCGGTCTGCCACCAGGTATCGACCACCGGACATTTCTCGTTGc                                                                              <  2:100285/250‑1 (MQ=255)                                                                                                                                                                        cAAAACGTTCGTGATCGCCAAACAGCGTACGCGCCTGACCCGCCCAGGAGTCGGTGATTACCAGGCTACCTTCGGTGGCCCCCCTCCAGCGGGTTCCC‑TTCGTTATCGACCAGCGCCGTTTGCACGCCGAAGAACGGACGTGTTGCCGAACCGGCTTTCACCTCGGTAGCGCCAGGCAGCGGGGTGATCATGAAACCGCCGGTTTCGGTCTGCCACCAGGTATCGACCACCGGACAGTTCTCGTTGCCGAt                                                                          <  1:202169/251‑1 (MQ=255)                                                                                                                                                                          aaaCGTTCGTGATCGCCAAACAGCGTACGCGCCTGACCCGGCCAGGAGTCGGTGATTACCAGGCTACCTTCGGTGGCCCCCCTCCAGCGGGTTACC‑TTCGTTATCGACCAGCGCCGGTTGCACGCCGAAGAACGGACGTGTTGCCGAACCGGCTTTCAGCTCGGTAGCGCCAGGCAGCGGGGTGATCATGAAACCGCCGGTTTCGGTCTGCCACCAGGTATCGACCACCGGACATTTCTCGTTGCCGAttt                                                                        <  1:262560/251‑1 (MQ=255)                                                                                                                                                                           aaCGTTCGTGATCGCCAAACAGCGTACGCGCCTGACCCGGCCAGGAGTCGGTGATTACCAGGCTACCTTCGGTGGCCCCCCTCCAGCGGGTTACC‑TTCGTTATCGACCAGCGCCGGTTGCACGCCGAAGAACGGACGTGTTGCCGAACCGGCTTTCAGCTCGGTAGCGCCAGGCAGCGGGGTGATCATGAAACCGCCGGTTTCGGTCTGCCACCAGGTATCGACCACCGGACATTTCTCGTTGCCGAtttt                                                                       <  1:2931/251‑1 (MQ=255)                                                                                                                                                                                 cgtgagccccAAACAGCGTCCGCGCCTGACCCGGCCAGGAGTCGGTGATTACCAGGCGACCTTCGGTGGCCCCCCTCCAGCGGTTTACC‑TTCGTTATCGCCCAGCGCCGGTTGCACGCCGAAGAACGGACGTGTTGCCGAACCGGCTTTCAGCTCGGTAGCGCCAGGCAGCGGGGTGATCATGAAACCGCCGGTTTCGGTCTGCCACCAGGTATCGACCACCGGACATTTCTCGTTGCCGATTTTTTTcc                                                                  <  1:97118/242‑1 (MQ=255)                                                                                                                                                                                        gCCAAACAGCGTACGCGCCTGACCCGGCCAGGAGTCGGTGATTACCAAGCTACCTTCGGTGGCCCCCCTCCAGCGGGTTCCC‑TTCGTTATCGACCAGCGCCGGTTGCACGCCGACGAACGGAAGTGTTGCCGAACCGGCTTTAAGCTCGGTAGCGCCAGGCAGCGGGGTGATCATGAAACCGCCGGTTTCGGTCTGCCACCAGGTATCGACCACCGGACATTTCTCGTTGCCGATTTTTTTCCAGTACCAc                                                          >  1:179435/1‑251 (MQ=255)                                                                                                                                                                                            aaCAGCGTACGAGCCTGACCCGGCCAGGAGTCGGTGATTACCAGGCTACCTTCGGTGGCCCCCCTCCAGCGGGTTCCC‑TTCGTTATCGACCAGCGCCGGTTGCACGCCGAAGAACGGACGTGTTGCCGAACCGGCTTTCAGCTCGGTAGCGCCAGGCAGCGGGGTGATCATGAAACCGCCGGATTAGGTCTGCCACCAGGTATCGACCACCGGACATTTCTCGTTGCCGATTTTTTTCCAGTACCACTccc                                                      >  2:243993/1‑251 (MQ=255)                                                                                                                                                                                                         ccTGACCCGCCCAGGAGTCGGTGATTACCAGGCTACCTTCGGTGGCCCCCCTCCAGCGGGTTACC‑TTCGTTATCGACCAGCGCCGGTTGCACGCCGAAGAACGGACGTGTTGCCGAACCGGCTTTCAGCTCGGTAGCGCCAGCCAGCGGGGTGATCATGAAACCGCCGGTTTCGGTCTGCCACCAGGTATCGACCACCGGACATTTCTCGTTGCCGATTTTTTTCCAGTACCACTCCCACGCTTCCGGGtt                                         <  1:186827/251‑1 (MQ=255)                                                                                                                                                                                                                                        gcaccctgcGGTGGCCCCC‑CCCCGGGGGTTACCCTGCGTTACCACCCACCGCCGGTGGCACGATGAAGAGCGGCCTGGTTACCTAACCGGTTTTCAGCGCGGTACCGCCAGGCAGCGGGGTGATCATGAAACCGCCGGTTTCGGTCTGCCACCAGGTATCGACCACCGGACATTTCTCGTTGCCGATTTTTTTCCAGTACCACTCCCACGCTGCCGGGTGAATTGGCTCGCGCACGGAACCGAGAATgcgc          <  2:173075/243‑1 (MQ=255)                                                                                                                                                                                                                                              ttCGGTGGCCCCCCTCCAGCGGGTTACC‑TTCGTTATCGACCAGCGCCGGTTGCACGCCGAAGAACGGACGTGTTGCCGAACCGGCTTTCAGCTCGGTAGCGCCAGGCAGCGGGGTGATCATGAAACCGCCGGTTTCGGTCTGCCACCAGGTATCGACCACCGGACATTTCTCGTTGCCGATTTTTTTCCAGTACCACTCCCACGCTTCCGGGTTAATTGGCTCGCCCACGGAACCGAGAATGCGCAGcgac    >  1:26687/1‑251 (MQ=255)                                                                                                                                                                                                                                                 ggTGGCCCCCCTCCAGCGGGTTACC‑TTCGTTATCGACCAGCGCCGGTTGCACGCCGAAGAACGGACGTGTTGCCGAACCGGCTTTCAGCTCGGTAGCGCCAGGCAGCGGGGTGATCATGAAACCGCCGGTTTCGGTCTGCCACCAGGTATCGACCACCGGACATTTCTCGTTGCCGATTTTTTTCCAGTACCACTCCCACGCTTCCGGGTTAATTGGCTCGCCCACGGAACCGAGAATGCGCAGcgacga  >  1:221150/1‑250 (MQ=255)                                                                                                                                                                                                                                                           |                                                                                                                                                                                                                                                 CGGATGCGCCACCAGCGCCGACTCAATCTCTGCCGTCCCCAGACGGTGACCGGAGACGTTCAGCACGTCGTCCACACGCCCGGTTATCCAGTAATAGCCATCTTCATCGCGACGCGCGCCGTCGCCGCTGAAATACATATTTTTGAAGGTGGAGAAGTAGGTCTGTTCAAAACGTTCGTGATCGCCAAACAGCGTACGCGCCTGACCCGGCCAGGAGTCGGTGATTACCAGGCTACCTTCGGTGGCCCCC‑TCCAGCGGGTTACC‑TTCGTTATCGACCAGCGCCGGTTGCACGCCGAAGAACGGACGTGTTGCCGAACCGGCTTTCAGCTCGGTAGCGCCAGGCAGCGGGGTGATCATGAAACCGCCGGTTTCGGTCTGCCACCAGGTATCGACCACCGGACATTTCTCGTTGCCGATTTTTTTCCAGTACCACTCCCACGCTTCCGGGTTAATTGGCTCGCCCACGGAACCGAGAATGCGCAGCGACGA  >  NC_000913/4285752‑4286240

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