Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	NC_000913	2,038	C→T	L568L (CTG→TTG)	thrA →	fused aspartate kinase/homoserine dehydrogenase 1

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	2,038	0	C	T	89.5%	85.4 / 0.3	38	L568L (CTG→TTG)	thrA	fused aspartate kinase/homoserine dehydrogenase 1
Reads supporting (aligned to +/- strand):  ref base C (3/1);  new base T (16/18);  total (19/19)
Fisher's exact test for biased strand distribution p-value = 6.04e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 6.39e-01

GGCTGAAGAATAAACATATCGACTTACGTGTCTGCGGTGTTGCCAACTCGAAGGCTCTGCTCACCAATGTACATGGCCTTAATCTGGAAAACTGGCAGGAAGAACTGGCGCAAGCCAAAGAGCCGTTTAATCTCGGGCGCTTAATTCGCCTCGTGAAAGAATATCATCTGCTGAACCCGGTCATTGTTGACTGCACTTCCAGCCAGGCAGTGGCGGATCAATATGCCGACTTCCTGCGCGAAGGTTTCCACGTTGTCACGCCGAACAAAAAGGCCAACACCTCGTCGATGGATTACT‑ACCATCAGTTGCGTTATGCGGCGGAAAAATCGCGGCGTAAATTCCTCTATGACACCAACGTTGGGGCTGGATTACCGGTTATTGAGAACCTGCAAAATCTGCTCAATGCAGGTGATGAATTGATGAAGTTCTCCGGCATTCTTTCTGGTTCGCTTTCTTATATCTTCGGCAAGTTAGAC  >  NC_000913/1805‑2280                                                                                                                                                                                                                                          |                                                                                                                                                                                                                                                    ggctgAAGAATAAACATATCGACTTACGTGTCTGCGGTGGTGCCAACTCGAAGGCACTGCGCCCCAATGTGCATGGCCTAAATCTGGAAAACTGGCAGGAAGAACTGGCGCAAGCCAAAGAGCCGTTTAATCTCGGGCGCTTAATTCGCCTCGTGAAAGAATATCATCTGCTGAACCCGGTCATTGTTGACTGCACTTCCAGCCAGGCAGTGGCGGATCAATATGCCGACTTCTTGCGCGAAGGTTTCCAc                                                                                                                                                                                                                                    <  2:246174/251‑1 (MQ=255)           tAAACATATCGACTTACGTGTCTGCGGTGTTGCCAACTCGAAGGCACTGCTCACCAATGTGCATGGCCTAAATCTGGAAAACTGGCAGGAAGAACTGGCGCAAGCCAAAGAGCCGTTTAATCTCGGGCGCTTAATTCGCCTCGTGAAAGAATATCATCTGCTGAACCCGGTCATTGTTGACTGCACTTCCAGCCAGGCAGTGGCGGATCAATATGCCGACTTCTTGCGCGATGGTTTCCACGTTGTCACGc                                                                                                                                                                                                                          >  1:50645/1‑251 (MQ=255)               cATATCGACTTACGTGTCTGCGGTGTTGCCAACTCGAAGGCACTGCTCACCAATGTGCATGGCCTAAATCTGGAAAACTGGCAGGAAGAACTGGCGCAAGCCAAAGAGCCGTTTAATCTCGGGCGCTTAATTCGCCTCGTGAAAGAATATCATCTGCTGAACCCGGTCATTGTTGACTGCACTTCCAGCCAGGCAGTGGCGGATCAATATGCCGACTTCTTGCGCGAAGGGTTCCACGTTGTCACGCCGa                                                                                                                                                                                                                       >  1:2116/1‑250 (MQ=255)                 tatCGACTTACGTGTCTGCGGTGTTGCCAACTCGAAGGCACTGTTCACCAATGTGAATGGCCTAAATCTGGAACACTGGAAGGAAGAACTGGCGCAAGCCAAAGAGCCGTTTAATCTCGGGCGCTTAATTCGCCTCGTGAAAGAATATCATCTGCTGAACCCGGTCATTGTTGACTGCACTTCCAGCCAGGCAGTGGCGGATCAATATGCCGCCTTCTTGCGCGGAGGGTTCCCCGTTTTCACGCCGAAcc                                                                                                                                                                                                                    >  1:130967/1‑250 (MQ=255)                                    ggTGTTGCCAACTCGAAGGCACTGCTCACCAATGTGCATGGCCTAAATCTGGAAAACTGGCAGGAAGATCTGGCGCAAGCCAAAGAGCCGTTTAAGCTCGGGCGCTTAATTCGCCTCGTGAAAGAATATCATCTGCTGAACCCGGTCATTGTTGACTGCACTTCCAGCCAGGCAGTGGCGGATCAATATGCCGTCTTCTTGCGCGAAGGTTTCCACGTTGTCACGCCGAACAAAAAGGCCATCACCtcgtc                                                                                                                                                                                                 <  1:252619/251‑1 (MQ=255)                                          gCCAACTCGAAGGGACTGCTCACCAATGTGCATGGCCTAAATCTGGAAAACTGGCAGGAAGAACTGGCGCAAGCCAAAGAGCCGTTTAATCTCGGGCGCTTAATGCGCCTCGTGAAAGAATATCATCTTCTGAACCCGGTCATTGTTGACTGCACTTCCAGCCAGGCATTGGCGGATCAATATGCCTACTTCTTGCGCGAAGGTTTCCACGTTGTCACGCCGAACAAAAAGGCCATCACCTCGTCGATGGa                                                                                                                                                                                           <  1:184039/251‑1 (MQ=255)                                                     ggcaatgcgcACCAATGTGTAAGGCCTAAATCTGGAAAACTGGCTGGAAGAACGTGCGTAAGCCAAAGAGCCGTGTAATCTCGGGCGCTTAATTCGCCTCGTGAAAGAATATTATCTGCTGAACCCGGTCATTGTTGACTGCACTCCCAGCCAGGCAGTGGCGGATCAATATGCCGACTTCTAGCGCTAAGGTTTCCACGTTGTCACGCCGAACAAAAAGTCAAACGCCTCGGCGATCGATGACT‑ACCATc                                                                                                                                                                               <  2:165951/242‑1 (MQ=255)                                                         ctGCTCACCAATGTGCATGGCCTAAATCTGGAAAACTGGCAGGAAGAACTGGCGCAAGCCAAAGAGCCGTTTAATCTCGGGCGATTAATTCGCCTCGTGAAAGAATATCATCTGCTGAACCCGGTCATTGTTGACTGCACTTCCAGCCAGGCAGTGGCGGATCAATATGCCCACTTCTTGCGCGAAGGTTTCCACGTTGTCACGCCGAACAAAAAGGCCAACACCTCGTCGATGGATTACT‑ACCATctgat                                                                                                                                                                           >  1:67864/1‑247 (MQ=255)                                                           gCTCACCAATGTGCATGGCCTAAATCTGGAAAACTGGCAGGAAGAACTGGCGCAAGCCAAAGAGCCGTTTAATCTCGGGCGCTTAATTCGCCTCGTGAAAGAATATCATCTGCTGAACCCGGTCATTGTTGACTGCACTTCCAGCCAGGCAGTGGCGGATCAATATGCCGACTTCTTGCGCGAAGGTTTCCACGTTGTCACGCCGAACAAAAAGGCCAACACCTCGTCGATGGATTACT‑ACCATCTGtt                                                                                                                                                                           <  1:376086/249‑1 (MQ=255)                                                           gCTCACCAATGTGCATGGCCTAAATCTGGAAAACTGGCAGGAAGAACTGGCGCAAGCCAAAGAGCCGTTTAATCTCGGGCGCTTAATTCGCCTCGTGAAAGAATATCATCTGCTGAACCCGGTCATTGTTGACTGCACTTCCAGCCAGGCAGTGGCGGATCAATATGCCGACTTCTTGCGCGAAGGTTGCCACGTTGTCACGCCGAACAAAAAGGCCAACACCTCGTCGATGGATTACT‑ACCATCTGtt                                                                                                                                                                           <  1:171287/249‑1 (MQ=255)                                                                ccAATGTGCATGGCCTAAATCTGGAAAACTGGCAGGAAGAACTGGCGCAAGCCAAAGAGCCGTTTAATCTCGGGCGCTTAATTCGCCTCGTGAAAGAATATCATCTGCTGAACCCGGTCATTGTTGACTGCACTTCCAGCCAGGCAGTGGCGGATCAATATGCCGACTTCTTGCGCGAAGGTTTCCACGTTGTCACGCCGAACAAAAAGGCCAACACCTCGTCGATGGATTACT‑ACCATCTGTTGCGtcat                                                                                                                                                                    <  1:159096/251‑4 (MQ=255)                                                                       gcATGGCCTAAATCTGGAAAACTGGCAGGAAGAACTGGCGCAAGCCAAAGAGCCGTTTAATCTCGGGCGCTTAATTCGCCTCGTGAAAGAATATCATCTGCTGAACCCGGTCATTGTTGACTGCACTTCCAGCCAGGCAGTGGCGGATCAATATGCCGACTTCTTGCGCGAAGGTTTCCACGTTGTCACGCCGAACAAAAAGGCCAACACCTCGTCGATGGATTACT‑ACCATCTGTTGCGTCATGcggcgg                                                                                                                                                             <  1:328530/250‑1 (MQ=255)                                                                                aaaTCTGGAAAACTGGCAGGAAGAACTGGCGCAAGCCAAAGAGCCGTTTAATCTCGGGCGCTTAATTCGCCTCGTGAAAGAATATCATCTGCTGAACCCGGTCATTGTTGACTGCACTTCCAGCCAGGCAGTGGCGGATCAATATGCCGACTTCTTGCGCGAAGGTTTCCACGTTGTCACGCCGAACAAAAAAGCCAAAACCTCGTCGATGGATTACT‑ACCATCTGTTGCGTCATGCGGCGGGAAAATcgc                                                                                                                                                    >  1:3434/2‑251 (MQ=255)                                                                                       gacatcgggCAGGAAGAACGGGCGCATGCCAAAGAGCGTTTTAACCTCTGGTCATTAATTTGCCACGTTAACTAATATCATCTGCTGAACCCGGTCAGTGGTGACTGCACTTCCAGCCAGGCAGTGGCGGAGCAATATGCCGACTTCTTGCGCTAATGTTTCCACGTTGGCCCGCCGAACAAAAAGGCCAACACCGCGTCGATGTATTACT‑ACCATCTGTTGCGTCATGCGGCGGAAAAATCGCGGCGTaa                                                                                                                                             <  2:338677/244‑1 (MQ=255)                                                                                       gAAAACTGGCAGGAAGAACTGGCGCAAGCCAAAGAGCCGTTTAATCTCGGGCGCTTAATTCGCCTCGTGAAAGAATATCATCTGCTGAACCCGGTCATTGTTGACTGCACTTCCAGCCAGGCAGTGGCGGATCAATATGCCGACTTCTTGCGCGAAGGTTTCCACGTTGTCACGCCGAACAAAAAGGCCAACACCTCGTCGATGGATTACT‑ACCATCTGTTGCGTCATGCGGCGGAAAAATCGCGGCGTaa                                                                                                                                             <  1:51367/251‑1 (MQ=255)                                                                                                      gaaCTGGCGCAAGCCAAAGAGCCGTTTAATCTCGGGCGCTTAATTCGCCTCGTGAAAGAATATCATCTGCTGAACCCGGTCATTGTTGACTGCACTTCCAGCCAGGCAGTGGCGGATCAATATGCCGACTTCTTGCGCGAAGGTTTCCACGTTGTCACGCCGAACAAAAAGGCCAACACCTCGTCGATGGATTACT‑ACCATCTGTTGCGTCATGCGGCGGAAAAATCGCGGCGTAAATTCCTCTATGacac                                                                                                                              >  1:204573/1‑251 (MQ=255)                                                                                                              gcAAGCCAAAGAGCCGTTTAATCTCGGGCGCTTAATTCGCCTCGTGAAAGAATATCATCTGCTGAACCCGGTCATTGTTGACTGCACTTCCAGCCAGGCAGTGGCGGATCAATATGCCGACTTCTTGCGAGAAGGTTTCCACGGTGTTACGCCGAACAAAAAGGCCCAAACCTCGTCGATGGATTACT‑ACCATCTGTTGCGTCATGCGGCGGAAAAATCGCGGCGTAAATTCCTTTATGACACCAACGTTg                                                                                                                      >  1:333450/1‑251 (MQ=255)                                                                                                                     aaaGAGCCGTTTAATCGCGGGCGCTTAATTCGCCTCGTGAAAGAATATCATCTGCTGAACCCGGTCATTGTTGACTGCACTTCCAGCCAGGCAGTGGCGGATCAATATGCCGACTTCTTGCGCGAAGGTTTCCACGTTGTCACGCCGAACAAAAAGGCCAACACCTCGTCGATGGATTACT‑ACCATCTGTTGCGTCATGCGGCGGAAAAATCGCGGCGTAAATTCCTCTATGACACCAACGTTGGGGCTgg                                                                                                               <  2:163939/251‑1 (MQ=255)                                                                                                                       agagCCGTTTAATCTCGGGCGCTTAATTCGCCTCGTGAAAGAATATCATCTGCTGAACCCGGTCATTGTTGACTGCACTTCCAGCCAGGCAGTGGCGGCACAATAAGTCGGCTTCTTGCGCGAAGGTTTCCACGTTGTCACGACGGACAAAAAAGCCAAACCCTCGTAGCTGGGTTACTTACCA‑CTGTTGCGTCATGACGCGGAAAAATCGCGGCGTAACTTCCACTATGACACCAACGTTGGggatggat                                                                                                             >  2:349742/1‑251 (MQ=255)                                                                                                                         agCCGTTTAATCTCGGGCGCTTAATTCGCCTCGTGAAAGAATATCATCTGCTGAACCCGGTCATTGTTGACTGCACTTCCAGCCAGGCAGTGGCGGATCAATATGCCGACTTCTTGCGCGAAGGTTTCCAAGTTGTCACGCCGAACAAAAAGGCCAACACCCCGTCGCTGGCTTACT‑AACATCTGTTGCGTCATGCGGCGGAAAAATAGCGGCGTAAATTCCTCTATGACCACAACGTTGGGGCTGGAtt                                                                                                            >  1:50060/1‑250 (MQ=255)                                                                                                                                 aaTCTCGGGCGCTTAATTCGCCTCGTGAAAGAATATCATCTGCTGAACCCGGTCATTGTTGACTGCACTTCCAGCCAGGCAGTGGCGGATCAATATGCCGACTTCTTGCGCGAAGGTGTCCACGTTGTCACGCCGAACAAAAAGGCCAACTCCTCGTCGATTGATTACT‑ACCTTCTGTTGCGTCATGCGGCGGAAAAATCGCGGCGTAAAGTCCTCGATGACTCCAACGTTGGGGCTGGATTACCGGTTa                                                                                                    <  1:306139/250‑1 (MQ=255)                                                                                                                                        ggCGCTTAATTCGCCTCGTGAAAGAATATCATCTGCTGAACCCGGTCATTGTTGACTGCACTTCCAGCCAGGCAGTGGCGGATCAATATGCCGACTTCCTGCGCGAAGGTTGCCACGTTGTCACGCCGAACAAAAAGGCCAACACCTCGTCGATGGATTACT‑ACCATCAGtt                                                                                                                                                                           <  1:168848/172‑1 (MQ=255)                                                                                                                                        ggCGCTTAATTCGCCTCGTGAAAGAATATCATCTGCTGAACCCGGTCATTGTTGACTGCACTTCCAGCCAGACAGTGGCGGCTCAATATGCCGACTTCCTGCGCGAAGGTTTCCACGTTGTCACGCCGAACAAAAAGGCCAACACCCCGTCGATGGATTACT‑AACATccggt                                                                                                                                                                           >  2:168848/1‑168 (MQ=255)                                                                                                                                                              ataaTATCATCTGCTGAACCCGGTCATTGTTGACTGCACTTCCAGCCAGGCAGTGGCGGATCAATATGCCGACTTCTTGCTCGAAGGTTTCCACGTTGTCACGCCGAACAAAAAGGCCAACACCTCGTCGATGGATTACT‑ACCATCTGTTGCGTCATGCGGCGGAAAAATCGCGGCGTAAATTCCTCTATGACACCAACGTTGGGGCTGGATTACCTGTTATTGAGAACCTGCAAAATCTGCTCAATGctg                                                                      <  1:116819/249‑3 (MQ=255)                                                                                                                                                               gAATATCATCTGCTGAACCCGGTCATTGTTGACTGCACTTCCAGCCAGGCAGTGGCGGATCAATATGCCGACTTCTTGCGCGAAGGTTTCCACGTTGTCACGCCGAACAAAAAGGCCAAAACCTCGTCGGTGGGTTACT‑ACCATCTGTTTCGTCATACGGAGGAAAAATAGCGGCGTAAAATCATCTATGGCACCAAAGTTGGGGCTGGATTACCGGTTATTGAGAACCTGACAAATCGGCTCAATGctgg                                                                     >  1:179705/1‑248 (MQ=255)                                                                                                                                                               gAATATCATCTGCTGAACCCGGTCATTGTTGACTGCACTTCCAGCCAGGCAGTGGCGGATCAACATGCCGACTTCTTTCGCGAAGGTTTCCCCGTTGTCACGCCGAACAAAAAAGCCACCACCTCGTCGGTGGATTCAT‑CACATCTGTTGCGGCATGCGGCGGACAAATCCCGGCGCAACTTACTCTTTAACACCAACGTTGGGGCGGGATTTCCGGTTATTTATACACTGCAAAAACTGCCCAATGctgg                                                                     >  2:37166/1‑248 (MQ=255)                                                                                                                                                                          gctgAACCCGGTCATTGTTGACCGCACGTCCAGCCAGGCAGTGGCGGATCAATATGCCGACTTCTTGCGCGAAGGTTTCCACGTTGTCACGCCGAACAAAAAGGCCAACACCTCGTCGATGGATTACT‑ACCATCTGTTGCGTCATGCGGCGGAAAAATCGCGGCGTAAATTCCTCTATGACACCAACGTTGGGGCTGGATTACCGGTTATTGAGAACCTGCAAAATCTGCTCAATGCTGGTGATGAAtt                                                            <  1:72705/249‑1 (MQ=255)                                                                                                                                                                            tgAACCCGGTCATTGTTGACTGCACTTCCAGCCAGGCAGTGGCGGATCAATATGCCGACTTCTTGCGCGAAGGTTTCCACGTTGTCACGCCGAACCAAAAGGCCAACACCTCGTCGATGGATTAAT‑ACCATCTGTTGCGTCATGCGGCGGAAAAATCGCGGCGTAAATTCCTCTATGACACCAACGTTGGGGCTGGATTACCGGTTATTGAGAACCTGCAAAATCTGCTCAATGCTGGGGATGAATtgatg                                                        >  1:365732/1‑251 (MQ=255)                                                                                                                                                                                cccGGTCATTGTTGACTGCACTTCCAGCCAGGCAGTGGCGGATCAATATGCCGACTTCCTGCGCGAAGGTTTCCCCGTTGTCACGCCGACCAAAAAGGCCAACACCTCGTCGATGGATTACT‑CCCATCAGTTGTGTTATGCGGCGGAAAAATCGCGGCGTAAATTCCTCTCTGACACCAACGTTGGGGCTGGGTTACCGGTTATTGAGACCCTGCAAAATCTGCTCCATGTAGGTGATGCATTGATGAAGt                                                    >  1:47531/1‑251 (MQ=255)                                                                                                                                                                                   ggTCATTGTTGACTGCACTTCCAGCCAGGCAGTGGCGGATCAATATGCCGACTTCTTGCGCGAAGGTGTCCACGTTGTCACGCCGAACAAAAAGGCCAACACCTCGTCGATGGATTACT‑ACCATCTGTTGCGTCATGCGGCGGAAAAATCGCGGCGTAAATTCCTCTATGACACCAACGTTGGGGCTGGATTACCGGTTATTGAGAACCTGCAAAATCTGCTCAATGCTGGTGATGTATTGATGAAGTtct                                                 <  1:305939/251‑1 (MQ=255)                                                                                                                                                                                    ttCATTGTTTACTGCACTTCCAGCCAGGCAGTGGCGGATCAATATACCGACTTCTTGCGCGAAGGGTTCCACGTTGTCACGCCAAACAAAAAGGCCAAAACCTCGTCGTTGTATTAAT‑ACCATCTGTTAGGTCATACTGCGGAAAAATCGCGGAGTATATTACTATTAGAAAACAACGTTGGGACAGGATTAACGGTGATTGAGAAACCGAACAATCTGCACAGTGCTGGTAGTGAATTTATGAAgagctc                                                >  2:264943/2‑246 (MQ=255)                                                                                                                                                                                                  cACTTCCAGCCAGGCAGTGGCGGATCAATATGCCGACTTCTTGCGCGAAGGTTTCCACGTTGTCACGCCGAACAAAAAGGCCAACACCTCGTCGATGGATTACT‑ACCATCTGTTGCGTCATGCGGCGGAAAAATCGCGGCGTAAATTCCTCTATGACACCAACGTTGGGGCTGGATTACCGGTTATTGAGAACCTGCAAAATCTGCTCAATGCTGGTGATGAATTGATGAAGTTCTCCGGCAttctttcag                                  <  1:154121/251‑3 (MQ=255)                                                                                                                                                                                                            cagGCAGTGGCGGATCAATATGCCGACTTCCTGCGCGAAGGTTTCCACGTTGTCACGCCGAACAAAAAGGCCAACACCTCGTCGATGGATTACT‑ACCATCAGTTGCGTTATGCGGCGGAAAAATCGCGGCGTAAATTCCTCTATGACACCAACGTTGGGGCTGGATTACCGGTTATTGAGAACCTGCAAAATCTGCTCAATGCAGGTGATGAATTGATGAAGTTCTCCGGCATTCTTTCTGGTTCGCTTTc                        >  1:193580/1‑251 (MQ=255)                                                                                                                                                                                                                             atatGCCGACTTCTTGCGCGAAGGTTTCCACGTTGTCACGCCGAACAAAAAGGCCAACACCTCGTCGATGGATTACT‑ACCATCTGTTGCGTCATGCGGCGGAAAAATCGCGGCGTAAATTCCTCTATGACACCAACGTTGGGGCTGGATTACCGGTTATTGAGAACCTGCAAAATCTGCTCAATGCTGGTGATGAATTGATGAAGTTCTCCGGCATTCTTTCAGGTTCGCTTTCTTATATCTTCGGCAAGt       >  1:62486/1‑251 (MQ=255)                                                                                                                                                                                                                              aatGCCGACTTCTTGCGCGAAGGTTTCCACGTTGTCACGCCGAACAAAAAGGCCAACACCTCGTCGATTGATTACT‑ACCATCTGTTGCGTCATGCGGCGGAAAAATCGCGGCGTAAATTCCTCTATGACACCAACGTTGGGGCTGGATTACCGGTTATTGAGAACCTGCAAAATCTGCTCAATGCTGGTGATGAATTGATGAAGTTCTCCGGCATTCTTTCAGGTTCGCTTTCTTATATCTTCGGCAAGtt      <  2:62486/250‑1 (MQ=255)                                                                                                                                                                                                                               atGCCGACTTCTTGCGCGAAGGTTTCCACGTTGTCACGCCGAACAAAAAGGCCAACACCTCGTCGATGGATTACT‑ACCATCTGTTGCGTCATGCGGCGGAAAAATCGCGGCGTAAATTCCTCTATGACACCAACGTTGGGGCTGGATTACCGGTTATTGAGAACCTGCAAAATCTGCTCAATGCTGGTGATGAATTGATGAAGTTCTCCGGCATTCTTTCAGGTTCGCTTTCTTATATCTTCGGCAAGtt      >  1:100147/1‑250 (MQ=255)                                                                                                                                                                                                                               atGCCGAATTCTTGCGCGAAGGTTTCCACGTTGTCACGCCGAACAAAAAGGCCAACACCTCGTCGATGGATTACT‑ACCATCTGTTGCGTCATGCGGCGGAAAAATCGCGGCGTAAATTCCTCTATGACACCAACGTTGGGGCTGGATTACCGGTTATTGAGAACCTGCAAAATATGCTCAATGCTGGTGATGAATTGATGAAGTTCTCCGGCATTCTTTCAGGTTCGCTTTCTTATATCTTCGGCACGttt     >  2:161573/1‑250 (MQ=255)                                                                                                                                                                                                                                  ccGACTTCTTGCGCGAAGGTTTCCACGTTGTCACGCCGAACAAAAAGGCCAACACCTCGTCGATGGATTACT‑ACCATCTGTTGCGTCATGCGGCGGAAAAATCGCGGCGTAAATTCCTCTATGACACCAACGTTGGGGCTGGATTACCGGTTATTGAGAACCTGCAAAATCTGCTCAATGCTGGTGATGAATTGATGAAGTTCTCCGGCATTCTTTCAGGTTCGCTTTCTTATATCTTCGGCAAGTTAGAc  <  1:37166/251‑1 (MQ=255)                                                                                                                                                                                                                                          |                                                                                                                                                                                                                                                    GGCTGAAGAATAAACATATCGACTTACGTGTCTGCGGTGTTGCCAACTCGAAGGCTCTGCTCACCAATGTACATGGCCTTAATCTGGAAAACTGGCAGGAAGAACTGGCGCAAGCCAAAGAGCCGTTTAATCTCGGGCGCTTAATTCGCCTCGTGAAAGAATATCATCTGCTGAACCCGGTCATTGTTGACTGCACTTCCAGCCAGGCAGTGGCGGATCAATATGCCGACTTCCTGCGCGAAGGTTTCCACGTTGTCACGCCGAACAAAAAGGCCAACACCTCGTCGATGGATTACT‑ACCATCAGTTGCGTTATGCGGCGGAAAAATCGCGGCGTAAATTCCTCTATGACACCAACGTTGGGGCTGGATTACCGGTTATTGAGAACCTGCAAAATCTGCTCAATGCAGGTGATGAATTGATGAAGTTCTCCGGCATTCTTTCTGGTTCGCTTTCTTATATCTTCGGCAAGTTAGAC  >  NC_000913/1805‑2280

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