Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	NC_000913	3,152,494	G→T	A87S (GCG→TCG)	metC →	cystathionine beta‑lyase/L‑cysteine desulfhydrase/alanine racemase

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	3,152,494	0	G	T	97.4%	123.0 / NA	39	A87S (GCG→TCG)	metC	cystathionine beta‑lyase/L‑cysteine desulfhydrase/alanine racemase
Reads supporting (aligned to +/- strand):  ref base G (0/1);  new base T (28/9);  total (28/11)
Fisher's exact test for biased strand distribution p-value = 2.63e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

TTGATACTCAACTGGTGAATGCAGGACGCAGCAAAAAATACACTCTCGGCGCGGTAAATAGCGTGATTCAGCGCGCTTCTTCGCTGGTCTTTGACAGTGTAGAAGCCAAAAAACACGCGACACGTAATCGCGCCAATGGAGAGTTGTTCTATGGACGGCGCGGAACGTTAACCCATTTCTCCTTACAACAAGCGATGTGTGAACTGGAAGGTGGCGCAGGCT‑GCGTGCTATTTCCCTGCGGGGCGGCAGCGGTTGCTAATTCC‑ATTCTTGCTTTTATCGAACAGGGCGATCATGTGTTGATGACCAACACCGCCTATGAACCGAGTCAGGATTTCTGTAGCAAAATCCTCAGCAAACTGGGCGTAACGACATCATGGTTTGATCCGCTGATTGGTGCCGATATCGTTAAGCATCTGCAGCCAAACACTAAAATCGTGTTTCTGGAATCGCCAGGCTCCATCACCATGGAAGTCCACGACG  >  NC_000913/3152252‑3152731                                                                                                                                                                                                                                                    |                                                                                                                                                                                                                                               ttGATACTCAACTGGTGAATGCAGGACGCAGCAAAAAATACACTCTCGGCGCGGTAAATAGCGTGATTCAGCGCGCTTCTTCGCTGGTCTTTGACAGTGTAGAAGCCAAAAAACACGCGACACGTAATCGCGCCAATGGAGAGTTGTTCTATGGACGGCGCGGAACGTTAACCCATTTCTCCTTACAACAAGCGATGTGTGAACTGGAAGGTGGCGCAGGCT‑GCGTGCTATTTCCCTGCGTGTCGGCAGCg                                                                                                                                                                                                                                        >  1:276030/1‑251 (MQ=255)  tGATACTCAACTGGTGAATGCAGGCCGCAGCAAAAAATACACTCTCGGCGCGGTAAATAGCGTGATTCAGCGCGCTTCTTCGCTGGTCTTTGACAGTGTAGAAGCCAAAAAACACGCGACACGTAATCGCGCCAATGGAGAGTTGTTCTATGGACGGCGCGGAACGTTAACCCATTTCTCCTTACAACAAGCGATGTGTGAACTGGAAGGTGGCGCAGGCT‑GCGTGCTATTTCCCTGCGGGTCGGCAGCgg                                                                                                                                                                                                                                       >  2:280172/1‑251 (MQ=255)             tGGTGAATGCAGGACGCAGCAAAAAATACACTCTCGGCGCGGTAAATAGCGTGATTCAGCGCGCTTCTTCGCAGGTCTTTGACAGTGTAGAAGCCAAAAAACACGCGACACGTAATCGCGCCAATGGAGAGTTGTTCTATGGACGGCGCGGAACGTTAACCCATTTCTCCTTACAACAAGCGATGAGTGAACTGGAAGGAGGCGCAGGCT‑GCGTGCTATTTCCCTGCGGGTCGGCAGCGGTTGCTAATTcc                                                                                                                                                                                                                            >  2:19136/1‑251 (MQ=255)                        ggACGCAGCAAAAAATACACTCTCGGCGCGGTAAATAGCGTGAATCAGCGCGCTTCTTCGCTGGTCTTTGACAGTGTAGAAGCCAAAAAACACGCGACACGTAATCGCGCCAATGGAGAGTTGTTCTATGGACGGCGCGGAACGTTAACCCATTTCTCCTTACAACAAGCGATGTGTGAACTGGAAGGTGGCGCATGCT‑GCGTGCTATTTCCCTGCGGGTCGGCAGCGGTTGCTAATTCC‑ATTCTTGCttt                                                                                                                                                                                                                <  2:276030/251‑1 (MQ=255)                         gACGCAGCAAAAAATACACTCTCGGCGCGGTAAATAGCGTGATTCAGCGCGCTTCTTCGCTGGTCTTTGACAGTGTAGAAGCCAAAAAACACGCGACACGTAATCGCGCCAATGGAGAGTTGTTCTATGGACGGCGCGGAACGTTAACCCATTTCTCCTTACAACAAGCGATGTGTGAACTGGAAGGTGGCGCAGGCT‑GCGTGCTATTTCCCTGCGGGTCGGCAGCGGTTGCTAATTCC‑ATTCTTGCtttt                                                                                                                                                                                                               >  2:270994/1‑251 (MQ=255)                         gACGCAGCAAAAAATACACTCTCGGCGCGGTAAATAGCGTGATTCAGCGCGCTTCTTCGCTGGTCTTTGACAGTGTAGAAGCCAAAAAACACGCGACACGTAATCGCGCCAATGGAGAGTTGTTCTATGGACGGCGCGGAACGTTAACCCATTTCTCCTTACAACAAGCGATGTGTGAACTGGAAGGTGGCGCAGGCT‑GCGTGCTATTTCCCTGCGGGTCGGCAGCGGTTGCTAATTCC‑ATTCTTGCtttt                                                                                                                                                                                                               >  2:250430/1‑251 (MQ=255)                             cagcaAAAAATACACTCTCGGCGCGGTAAATAGCGTGATTCAGCGCGCTTCTTCGCTGGTCTTTGACAGTGTAGAAGCCAAAAAACACGCGACACGTAATCGCGCCAATGGAGAGTTGTTCTATGGACGGCGCGGAACGTTAACCCATTTCTCCTTACAACAAGCGATGTGTGAACTGGAAGGTGGCGCAGGCT‑GCGTGCTATTTCCCTGCGGGTCGGCAGCGGTTGCTAATTCC‑ATTCTTGCTTTTATCg                                                                                                                                                                                                           <  1:142917/251‑1 (MQ=255)                              agcaATAAATACACACTCGGCGCGGTAAGTAGCGTGATTCGGCGCGCTTCTTCGCAGGTCTTTGACATTGTAGAAACCAAAAAACACGCGCCACGTAATCGAGCCAATGGAGAGTTGTTCTATGGACGGCGCGGAACGTTAACCCATTTCTCCTTTCAACAAGCGATGTTGGCACTAGAAGGTGGCGCAGGCT‑GAGTGCTATTTCCCTGCGGTTCGGCAGCGGTTGCGAAGTCA‑ATTTTTGCTTGTATCg                                                                                                                                                                                                           >  2:236108/1‑250 (MQ=255)                                   aaaaTACACTCTCGGCGCGGTAAATAGCGTGATTCAGCGCGCTGCTTCGCTGGTCTTTGACAGTGTAGAAGCCAAAAAACACGCGACACGTAATCGCGCCAATGGAGAGTTGTTCTATGGACGGCGCGGAACGTTAACCCATTTCTCCTTACAACAAGCGATGTGTGAACTGGAAGGTGGCGCAGGCT‑GCGTGCTATTTCCCTGCGGGTCGGCAGCGGTTGCTGATTCC‑ATTCTTGCTTTTATCGAACAgg                                                                                                                                                                                                     <  1:128728/251‑1 (MQ=255)                                    caaTACACTCTCGGCGCGGTAAATAGCGTGATTCAGCGCGCTACTTCGCTGGTCTTTGACCGTGTAGAAGCCAAAAAACACGCGACACGTAATCGCGCCAATGGAGAGTTGTTCTATGGACGGCGCGGAACGTTAACCCATTTCTCCTTACCACCAGCGATGTGTGAACTGGAAGGTGGAGCAGACT‑GCGTGCTATTTCACTGCGGGTCGGCAGAGGCTGCTAATTAC‑ATTCTTGCTTTTATCGatccggg                                                                                                                                                                                                    >  2:225162/2‑245 (MQ=255)                                                     ggTAAATAGCGTGATTCAGCGCGCTTCTTCGCTGGTCTTTGACAGTGTAGAAGCCAAAAAACACGCGACACGTAATCGCGCCAATGGAGAGTTGTTCTATGGACGGCGCGGAACGTTAACCCATTTCTCCTTACAACAAGCGATGTGTGAACTGGAAGGTGGCGCAGGCT‑GCGTGCTATTTCCCTGCGGGTCGGCAGCGGTTGCTAATTCC‑ATTCTTGCTGTTATCGAACAGGGCGATCATGTGTtgatga                                                                                                                                                                                   <  1:250430/251‑1 (MQ=255)                                                             gCGTGATTCAGCGCGCTTCTTCGCTGGTCTTTGACAGTGTAGAAGCCAAAAAACACGCGACACGTAATCGCGCCAATGGAGAGTTGTTCTATGGACGGCGCGGAACGTTAACCCATTTCTCCTTACAACAAGCGATGTGTGAACTGGAAGGTGGCGCAGGCT‑GCGTGCTATTTCCCTGCGGGTCGGCAGCGGTTGCTAATTCC‑ATTCTTGCTTTTATCGAACAGGGCGATCATGTGTTGATGACCAACAcc                                                                                                                                                                           <  1:19136/251‑1 (MQ=255)                                                               gTGATTCAGCGCGCTTCTTCGCTGGTCTTTGACAGTGTAGAAGCCAAAAAACACGCGACACGTAATCGCGCCAATGGAGAGTTGTTCTATGGACGGCGCGGAACGTTAACCCATTTCTCCTTACAACAAGCGATGTGTGAACTGGAAGGTGGCGCAGGCT‑GCGTGCTATTTCCCTGCGGGTCGGCAGCGGTTGCTAATTCC‑ATTCTTGCTTTTATCGAACAGGGCGATCATGTGTTGATGACCAACACCGc                                                                                                                                                                         >  1:72597/1‑251 (MQ=255)                                                                tGATTCAGCGCGCTTCTTCGCTGGTCTTTGACAGTGTAGAAGCCAAAAAACACGCGACACGTAATCGCGCCAATGGAGAGTTGTTCTATGGACGGCGCGGAACGTTAACCCATTTCTCCTTACAACAAGCGATGTGTGAACTGGAAGGTGGCGCAGGCT‑GCGTGCTATTTCCCTGCGGGTCGGCAGCGGTTGCTAATTCC‑ATTCTTGCTTTTATCGAACAGGGCGATCATGTGTTGATGACCAACACCGcc                                                                                                                                                                        >  1:146795/1‑251 (MQ=255)                                                                     cAGCGCGCTTCTTCGCTGGTCTTTGACAGTGTAGAAGCCAAAAAACACGCGACACGTAATCGCGCCAATGGAGAGTTGTTCTATGGACGGCGCGGAACGTTAACCCATTTCTCCTTACAACAAGCGATGTGTGAACTGGAAGGTGGCGCAGGCT‑GCGTGCTATTTCCCTGCGGGTCGGCAGCGGTTGCTAATTCC‑ATTCTTGCTTTTATCGAACAGGGCGATCATGTGTTGATGACCAACACCGCCTATg                                                                                                                                                                    >  1:85789/1‑250 (MQ=255)                                                                                   gCTGGTCTTTGACAGTGTAGAAGCCAAAAAACACGCGACACGTAATCGCGCCAATGGAGAGTTGTTCTATGGACGGCGCGGAACGTTAACCCATTTCTCCTTACAACAAGCGATGTGTGAACTGGAAGGTGGCGCAGGCT‑GCGTGCTATTTCCCTGCGGGTCGGCAGCGGTTGCTAATTCC‑ATTCTTGCTTTTATCGAACAGGGCGATCATGTGTTGATGACCAACACCGCCTATGAACCGAGTCAGGAtt                                                                                                                                                     >  1:87238/1‑251 (MQ=255)                                                                                              aCAGTGTAGAAGCCAAAAAACACGCGACACGTAATCGCGCCAATGGAGAGTTGTTCTATGGACGGCGCGGAACGTTAACCCATTTCTCCTTACAACAAGCGATGTGTGAACTGGAAGGTGGCGCAGGCT‑GCGTGCTATTTCCCTGCGGGTCGGCAGCGGTTGCTACTTCC‑ATTCTTGCTTTTATCGAACAGGGCGATCATGTGTTGATGACCAACACCGCCTATGAACCGAGTCAGGATTTCTGTAGCaa                                                                                                                                           >  2:11365/1‑250 (MQ=255)                                                                                                               aaaCTCGCGACACGTAATCGCGCCAATGGAGAGTTGTTCTATGGTGGGCGCGGACCGTTAACCCATTTCTCCTTCCAACAAGCGATGTGTAAACGGGAAGGTGGCGCAGGCG‑GCGTGCTATTCCCCTGGGGGTCGGCAGCGGTTGCTAATTCC‑ATTCTTGCTTTTATCGACCAGGGCGATCATTTGTGGATAACCAACACCGCCTATGAACCGAGTCAGGATTTCTGTAGCAAAATCCTCAGCAAACTggg                                                                                                                         <  2:238259/251‑1 (MQ=255)                                                                                                                     gcagcacGTAATCGCGTCAAGGGAGAGGTGTTCTTTGGCTGGGGCCGACCGTTCCCCCATTTCTCCTAACAACAGGCGATGTGTGTACTGGAATGTGCCG‑TGGCTTGTGTGATATTTCCTTGGTGGTCGTTAGCGGTTCCTAATTCC‑TTTCTTGCTTTTATCCACAAGGGCTATCATGTGTTTATGCCCAACACCGCCTATGAACCGATTCAGGATTTCTGTAGCAAAATCCTCAGCAAACTGGGCATAAc                                                                                                                   <  2:4897/247‑1 (MQ=255)                                                                                                                                   cccaaattgagtgttgttCTATGGACGGCGGGGATCGTTCCCCCCTTTCCCCTTACACCACGCTTTTTGTTCACTGTAAGGTGCTCCTGGCT‑GCGTGCTTTTCCCCTGCGGTTCGGCAGCT‑TTTCTTATTCCCTTTCTTCTTTTTATCGAACGGGGCGACCATGTTTTTATCACCACCACCGCCTATGAACCGAGTCAGGATTTCTGTAGCAAAATCCTCAGCAAACTGGGCGTAACGACATCATGGTTTg                                                                                                     <  1:247318/239‑1 (MQ=255)                                                                                                                                        aTGGAGAGTTGTTCTATGGACGGCGCGGAACGTTAACCCATTTCTCCTTACAACAAGCGATGTGTGAACTGGAAGGTGGCGCAGGCT‑GCGTGCTATTTCCCTGCGGGTCGGCAGCGGTTGCTAATTCC‑ATTCTTGCTTTTATCGAACAGGGCGATCATGTGTTGATGACCAACACCGCCTATGAACCGAGTCAGGATTTCTGTAGCAAAATCCTCAGCAAACTGGGCGTAACGACATCATGGTTTGATCCg                                                                                                >  1:106957/1‑251 (MQ=255)                                                                                                                                             gagTTGTTCTATGGACGGCGCGGAACGTTAACCCATTTCTCCTTACAACAAGCGATGTGTGAACTGGAAGGTGGCGCAGGCT‑GCGTGCTATTTCCCTGCGGGTCGGCAGCGGTTGCTAATTCC‑ATTCTTGCTTTTATCGAACAGGGCGATCATGTGTTGATGACCAACACCGCCTATGAACCGAGTCAGGATTTCTGTAGCAAAATCCTCAGCAAACTGGGCGTAACGACATCATGGTTTGATCCGCTGAt                                                                                           >  1:147170/1‑251 (MQ=255)                                                                                                                                                    tCTATGGACGGCGCGGAACGTTAACCCATTTCTCCTTACAACAAGCGATGTGTGAACTGGAAGGTGGCGCAGGCT‑GCGTGCTATTTCCCTGCGGGTCGGCAGCGGTTGCTAATTCC‑ATTCTTGCTTTTATCGAACAGGGCGATCATGTGTTGATGACCAACACCGCCTATGAACCGAGTCAGGATTTCTGTAGCAAAATCCTCAGCAAACTGGGCGTAACGACATCATGGTTTGATCCGCTGATTGGTGcc                                                                                    >  1:121019/1‑251 (MQ=255)                                                                                                                                                    tCTAGGGACGGGGCGAAACGTTACCCCATTCCCCCTTCCAACAAGCGTGGTGGGAACGGGAAGGTGGCGCAGGCG‑GCTTTTTTTTTCCCGGGGGGCCGGCGGGGTTTGTTAATTCC‑ATTCTTGCTTTTATGAAACAGGGCGAGCTTGTGTGGATGACCAACACCGCCTATGAACCGAGTCAGGATTTCTGTAGCAAAATCCTCAGCAAACTGGGCGTAACGACATCATGGTTTGATCCGCGGATTGGTGc                                                                                     <  2:328224/250‑1 (MQ=255)                                                                                                                                                     cTATGGACGGCGCAGAACGTTAACCCCTTTCTCCTTAACACAAGCGACGTGTGAACTGGAAGGTGGCGCCGGCT‑GCGTGCTCTTTCCCTGCGGGTCGGCAGCGGTTGCTAATTCC‑ATTCTTGCTTTTCTCGAACAGGCCGATCATGTGTTGATGACCACCACCGCCTATACACCGAGTCAGGCTTTCAGTAGCAAAATCCTCAGCAAACTGGGCGTAACGACATCATGGTTTGATCCGCTGATTGGTGCCg                                                                                   >  2:139915/1‑251 (MQ=255)                                                                                                                                                              acgcgGAACGTTAACCAATTTATCCATCCAACAAGCGATGAGTGAACTGGAAGGTGGCGCAGGCT‑GCGTGCTATTTCCCTGCGGGTCGGCAGCGGTTGCTACTTCC‑ATTCTTGATTTTATCGAACAGGGCGATCATGAGTTGATGACCAACACCGCCTAAGAACCGAGTAAGGATTTCTGTAGCAAAATCCTCAGCAAACTGGGCGTAACGACCTCATGGTTTGATCCGCTGATTGGTGCCGTTATCGtt                                                                           >  2:201118/2‑250 (MQ=255)                                                                                                                                                                 cgGAACGTTAACCCATTTATCCTTACAACAAGCGATGTGTGAACTGGAAGGAGGCGCAGGCT‑GCGTGCTATTTCCCTCCGGGTCGGCAGCGGTTCATAATTAC‑CTTATTGCTTTTATCGCACAGGGCGATCATGTGTTGATGACCAACACCGCCTATGAACCGAGTCAGGATTTCTGTAGCAAAATCCTCAGCAACCTGGGCGTAACGACATCATGGTTTGATCCGCTGATTGGTGCCGATATCGTTAAGc                                                                       >  1:296149/1‑251 (MQ=255)                                                                                                                                                                   gAACGTTAACCCATTTCTCCTTACAACAAGCGATGTGTGAACTGGAAGGTGGCGCAGGCT‑GCGTGCTATTTCCCTGCGGGTCGGCAGCGGTTGCTAATTCC‑ATTCTTGCTTTTATCGAACAGGGCGATCATGTGTTGATGACCAACACCGCCTATGAACCGAGTCAGGATTTCTGTAGCAAAATCCTCAGCAAACTGGGCGTAACGACATCATGGTTTGATCCGCTGATTGGTGCCGATATCGTTAAGCAt                                                                     >  1:64904/1‑251 (MQ=255)                                                                                                                                                                   gAACGTTAACCCATTTCTCCTTACAACAAGCGATGTGTGAACTGGAAGGTGGCGCAGGCT‑GCGTGCTATTTCCCTGCGGGTCGGCAGCGGTTGCTAATTCC‑ATTCTTGCTTTTATCGAACAGGGCGATCATGTGTTGATGACCAACACCGCCTATGAACCGAGTCAGGATTTCTGTAGCAAAATCCTCAGCAAACTGGGCGTAACGACATCATGGTTTGATCCGCTGATTGGTGCCGATATCGTTAAGCAt                                                                     >  1:220648/1‑251 (MQ=255)                                                                                                                                                                                            aaaaGCGATGTGTGATCTGGAAGGTGGCGCAGGCT‑GCGTGCTATTCCCCTGCGGGGCGGCAGCGGTTGCTAATTCC‑ATTCTTGCTTTTATCGAACAGGGCGATCATGTTTTGATGACCAACACCGCCTATGAACCGAGTCAGGATTTCTGTAGCAAAATCCTCAGCAAACTGGGCGTAACGACATCATGGTTTGATCCGCTGATTGGTGCCGATATCGTTAAGCATCTGCAGCCAAACACTAAAATCgtgt                                            <  2:72597/249‑1 (MQ=255)                                                                                                                                                                                              aaGCGATGTGTGAACTGGAAGGTGGCGCAGGCT‑GCGTGCTATTTCCCTGCGGGTCGGCAGCGGTTGCTAATTCC‑ATTCTTGCTTTTATCGAACAGGGCGATCATGTGTTGATGACCAACACCGCCTATGAACCGAGTCAGGATTTCTGTAGCAAAATCCTCAGCAAACTGGGCGTAACGACATCATGGTTTGATCCGCTGATTGGTGCCGCTATCGTTAAGCATCTGCAGCCAAACACTAAAATCGTGttt                                          >  2:206545/1‑251 (MQ=255)                                                                                                                                                                                                        tgAACTGGAAGGTGGCGCAGGCT‑GCGTGCTATTTCCCTGCGGGTCGGCAGCGGTTGCTAATTCC‑ATTCTTGCTTTTATCGAACAGGGCGATCATGTGTTGATGACCAACTCCGCCTATGAACCGAGTCAGGATTTCTGTAGCAAAATCCTCAGCAAACTGGGCGTAACGACATCATGGTTTGATCCGCTGATTGGTGCCGTTATCGCTAAGCATCTGCAGCCAAACACTAAAATCGTGTTTCTGGAATCGc                                >  2:134242/1‑251 (MQ=255)                                                                                                                                                                                                        tgAACTGGAAGGTGGCGCAGGCT‑GCGTGCTATTTCCCTGCGGGTCGGCAGCGGTTGCTAATTCC‑ATTCTTGCTTTTATCGAACAGGGCGATCATGTGTTGATGACCAACACCGCCTATGAACCGAGTCAGGATTTCTGTAGCAAAATCCTCAGCAAACTGGGCGTAACGACATCATGGTTTGATCCGCTGATTGGTGCCGATATCGTTAAGCATCTGCAGCCAAACACTAAAATCGTGTTTCTGGAATCGc                                >  2:26054/1‑251 (MQ=255)                                                                                                                                                                                                              ggAAGGTGGCGCAGGCT‑GCGTGCTATTTCCCTGCGGGTCGGCAGCGGTTGCTAATTCC‑ATTCTTGCTTTTATCGAACAGGGCGATCATGTGTTGATGACCAACACCGCCTATGAACCGAGTCAGGATTTCTGTAGCAAAATCCTCAGCAAACTGGGCGTAACGACATCATGGTTTGATCCGCTGATTGGTGCCGATATCGTTAAGCATCTGCAGCCAAACACTAAAATCGTGTTTCTGGAATCGCCAGGCt                          >  2:207240/1‑251 (MQ=255)                                                                                                                                                                                                                  ggTGGCGCAGGCT‑GCGTGCTATTTCCCTGCGGGTCGGCAGCGGTTGCTAATTCC‑ATTCTTGCTTTTATCGAACAGGGCGATCATGTGTTGATGACCAACACCGCCTATGAACCGAGTCAGGATTTCTGTAGCAAAATCCTCAGCAAACTGGGCGTAACGACATCATGGTTTGATCCGCTGATTGGTGCCGATATCGTTAAGCATCTGCAGCCAAACACTAAAATCGTGTTTCTGGAATCGCCAGGCTCCAt                      >  2:135761/1‑251 (MQ=255)                                                                                                                                                                                                                             cT‑GCGTGCTATTTCCCTGCGGGTCGGCAGCGGGTGCTAATTCC‑ATTCTTGCTTTTATCGAACAGGGCGATCATGTGTTGATGACCAACACCACCTATGAACCGCGTCCGGATTTCAGGAGCAACATCCACAACAAACTGGGCGTAACGACATCATGGTTTGATCCGCTGATTGGTGCCGATATCGTTAAGCATCTGCAGCAAAACACTAAAATCGGGTTTCTGGAATCGCCAGGCTCCATCACCATgaac            >  1:325832/1‑247 (MQ=255)                                                                                                                                                                                                                              t‑GCGTGCTATTTCCCTGCGGGTCGGCAGCGGTTGCTAATTCC‑ATTCTTGCTTTTATCGAACAGGGCGATCATGTGTTGATGACCAACACCGCCTATGAACCGAGTCAGGATTTCTGTAGCAAAATCCTCAGCAAACTGGGCGTAACGACATCATGGTTTGATCCGCTGATTGGTGCCGATATCGTTAAGCATCTGCAGCCAAACACTAAAATCGTGTTTCTGGAATCGCCAGGCTCCATCACCATGGAAGt          >  1:212757/1‑251 (MQ=255)                                                                                                                                                                                                                               ggcttcttttttCCCTGCGGGTCGGCAGCGGTTGCTAATTCC‑TTTCTTGTTTTTATAGACCAGGGCGTTCTGGTGTTTATGACCAACACCGCCTATGAACCGAGTCAGGATTTCTGTAGCAAAATCCTCAGCAAACTGGGCGTAACGACATCATGGTTTGATCCGCTGATTGGTGCCGATATCGTTAAGCATCTGCAGCCAAACACTAAAATCGTGTTTCTGGAATCGCCAGGCTCCATCACCATGGAAGt          <  2:64904/242‑1 (MQ=255)                                                                                                                                                                                                                                       aTTTCCCTGCGGGTCGGCAGCGGTTGCTAATTCC‑ATTCTTGCTTTTATCGAACAGGGCGATCATGTGTTGATGACCAACACCGCCTATGAACCGAGTCAGGATTTCTGTAGCAAAATCCTCAGCAAACTGGGCGTAACGACATCATGGTTTGATCCGCTGATTGGTGCCGATATCGTTAAGCATCTGCAGCCAAACACTAAAATCGTGTTTCTGGAATCGCCAGGCTCCATCACCATGGAAGTCCacgacg  >  2:215488/1‑251 (MQ=255)                                                                                                                                                                                                                                                    |                                                                                                                                                                                                                                               TTGATACTCAACTGGTGAATGCAGGACGCAGCAAAAAATACACTCTCGGCGCGGTAAATAGCGTGATTCAGCGCGCTTCTTCGCTGGTCTTTGACAGTGTAGAAGCCAAAAAACACGCGACACGTAATCGCGCCAATGGAGAGTTGTTCTATGGACGGCGCGGAACGTTAACCCATTTCTCCTTACAACAAGCGATGTGTGAACTGGAAGGTGGCGCAGGCT‑GCGTGCTATTTCCCTGCGGGGCGGCAGCGGTTGCTAATTCC‑ATTCTTGCTTTTATCGAACAGGGCGATCATGTGTTGATGACCAACACCGCCTATGAACCGAGTCAGGATTTCTGTAGCAAAATCCTCAGCAAACTGGGCGTAACGACATCATGGTTTGATCCGCTGATTGGTGCCGATATCGTTAAGCATCTGCAGCCAAACACTAAAATCGTGTTTCTGGAATCGCCAGGCTCCATCACCATGGAAGTCCACGACG  >  NC_000913/3152252‑3152731

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