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.3%	98.4 / NA	37	A87S (GCG→TCG)	metC	cystathionine beta‑lyase/L‑cysteine desulfhydrase/alanine racemase
Reads supporting (aligned to +/- strand):  ref base G (0/0);  major base T (23/13);  minor base C (1/0);  total (24/13)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

AGCTTGATACTCAACTGGTGAATGCAGGACGCAGCAAAAAATACACTCTCGGCGCGGTAAATAGCGTGATTCAGCGCGCTTCTTCGCTGGTCTTTGACAGTGTAGAAGCCAAAAAACACGCGACACGTAATCGCGCCAATGGAGAGTTGTTCTATGGACGGCGCGGAACGTTAACCCATTTCTCCTTACAACAAGCGATGTGTGAACTGGAAGGTGGCGCAGGCTGCGTGCT‑ATTTCCCTGCGGGGCGGCAGCGGTTG‑CTAATTCCATTCTTGCTTTTATCGAACAGGGCGATCATGTGTTGATGACCAACACCGCCTATGAACCGAGTCAGGATTTCTGTAGCAAAATCCTCAGCAAACTGGGCGTAACGACATCATGGTTTGATCCGCTGATTGGTGCCGATATCGTTAAGCATCTGCAGCCAAACACTAAAATCGTGTTTCTGGAATCGCCAGGCTCCATCACCATGGAAGTCCACGACGTTCCGGCG  >  NC_000913/3152249‑3152739                                                                                                                                                                                                                                                       |                                                                                                                                                                                                                                                       aGCTTGATACTCAACTGGTGAATGCAGGACGCAGCAAAAAATACCCTCTCGGCGCGGTAAATAGCGTGATTCAGCGCGCTTCTTCGCTGGTCTTTGACAGTGTAGAAGCCAAAAAACACGCGACACGTAATCGCGCCAATGGAGAGTTGTTCTATGGACGGCGCGGAACGTTAACCCATTTCTCCTTACAACAAGCGATGTGTGAACTGGAAGGTGGCGCAGGCTGCGTGCT‑ATTTCCCTGCGGGTCGGCa                                                                                                                                                                                                                                                   <  1:19997/251‑1 (MQ=255)       aTACTCAACTGGTGAATGCAGGACGCAGCAAAAAATACACTCTCGGCGCGGTAAATAGCGTGATTCAGCGCGCTTCTTCGCTGGTCTTTGACAGTGTAGAAGCCAAAAAACACGCGACACGTAATCGCGCCAATGGAGAGTTGTTCTATGGACGGCGCGGAACGTTAACCCATTTCTCCTTACAACAAGCGATGTGTGAACTGGAAGGTGGCGCAGGCTGCGTGCT‑ATTTCCCTGCGGGTCGGCAGCGGtt                                                                                                                                                                                                                                             <  2:7548/251‑1 (MQ=255)                 ggTGAATTCAGTACGCAGCAAAAAATACCCTCTCGGCGCGGGAAAGAGCGTTATTCAGCTCGCTTCTGCGCTGGTCTTTGACAGTGTAGAAGCCAAAAAACACGCGGCACGTAATCGCGCCAATGGTGAGGTGTTCTATGGACGGCGCGGAACGTTAACCCATTTCGCCTTACAACAATCGATGTGTGAACGGGAAGTTGGCGCAGGCTGCTTGCT‑ATATCCCTGCGGGTCGGCGTCGGGTG‑CTAATTCCa                                                                                                                                                                                                                                  <  2:40257/251‑1 (MQ=255)                  gTGAATGCAGGACGCAGCAAAAAATACACTCTCGGCGCGGTAAATAGCGTGATTCAGCGCGCTTCTTCGCTGGTCTTTGACAGTGTAGAAGCCAAAAAACACGCGACACGTAATCGCGCCAATGGAGAGTTGTTCTATGGACGGCGCGGAACGTTAACCCATTTCTCCTTACAACAAGCGATGTGTGAACTGGAAGGTGGCGCAGGCTGCGTGCT‑ATTTCCCTGCGGGTCGGCAGCGGTTG‑CTAATTCCAt                                                                                                                                                                                                                                 >  1:287650/1‑251 (MQ=255)                     aaTGCAGGACGCAGTAAAAAATACACACGCGGCGCGGTAAATAGCGTGATTCAGCGCGCTTCTTCGCTGGTCTGTGACAGGGTAGAAGCCAAAAAACACGCGACACTTAATCGCGCCAATGGAGAGTTGTTCTATGTACGGCGCGGAACGGTAACCCATTGCTCCTTACAACAAGCGATGTGTGAACTGGAAGGTGTCTCAGGCGGCGTGCT‑ATTTCCCTGCGGGTCGGCAGCATTTG‑CTAATTCCATTCt                                                                                                                                                                                                                              <  2:487209/251‑1 (MQ=255)                         cAGGACGCAGCAAAAAATACACTCTCGGCGCTGTAAATAGCGTGATTCAGCGCGCTTCTTCGCTGGTCTTTGACAGTGTAGAAGCCAAAAAACACGCGACACGTAATCGCTCCAATGGAGAGTTGTTCTATGGACGGCGCGGAACGTTAACCCATTTCTCCTTACAACAAGCGATGTGTGAACTGGAAGGTGGCGCAGGCTGCGTGCT‑ATTTCCCTGCGGGTCGTCAGCGGTTG‑CTAATTCCATTCTTGCt                                                                                                                                                                                                                          <  1:256548/251‑1 (MQ=255)                             aCGCAGCAAAAAATACACTCTCGGCGCGGTAAATAGCGTGATTCAGCGCGCTTCTTCGCTGGTCTTTGACAGTGTAGAAGCCAAAAAACACGCGACACGTAATCGCGCCAATGGAGAGTTGTTCTATGGACGGCGCGGAACGTTAACCCATTTCTCCTTACAACAAGCGATGTGTGAACTGGAAGGTGGCGCAGGCTGCGTGCT‑ATTTCCCTGCGGGTCGGCAGCGGTTG‑CTAATTCCATTCTTGCTTTTa                                                                                                                                                                                                                      <  1:360279/251‑1 (MQ=255)                               gcagcaAAAAATACACTCTCGGCGCGGTAAATAGCGTGATTCAGCGCGCTTCTTCGCTGGTCTTTGACAGTGTAGAAGCCAAAAAACACGCGACACGTAATCGCGCCAATGTAGAGTTGTTCTATGGACGGCGCGGAACGTTAACCCATTTCTCCTTACAACAAGCGATGTGTGAACTGGAAGGTGGCGCAGGCTGCGTGCT‑ATTTCCCTGCGGGTCGGCAGCGGTTG‑CTAATTCCATTCTTGCTTTTATc                                                                                                                                                                                                                    >  1:389722/1‑251 (MQ=255)                                             actctagacgcgGTAAATAGCGAGATTCAGCGCGCTTCTTCGCTGGTCTTTGAAAGTGTAGAAGCCAAAAAACACGCGACAAGTAATCGCGACAATGGAGAGTTGTTCTATGGACGGCGCGGAACGTTAACCCATTTCTCCTTACAACAAGCGATGTGTGAACTGGACGGTGGCGCAGGCAGCGTGCT‑ATTTCCCTGCGGGCCGGCAGCGGTTG‑CTAATTCCATTCTTGCTTTTATCGAACAGGGCGATc                                                                                                                                                                                                       >  2:339834/9‑250 (MQ=255)                                                   ggCGCGGTAAATAGCGTGATTCAGCGCGCTTCTTCGCTGGTCTTTGACAGTGTAGAAGCCAAAAAACACGCGACACGTAATCGCGCCAATGGAGAGTTGTTCTATGGACGGCGCGGAACGTTAACCCATTTCTCCTTACAACAAGCGATGTGTGAACTGGAAGGTGGCGCAGGCTGCGTGCT‑ATTTCCCTGCGGGTCGGCAGCGGTTG‑CTAATTCCATTCTTGCTTTTATCGAACAGGGCGATCATGTGtt                                                                                                                                                                                                >  2:483673/1‑251 (MQ=255)                                                            aaTAGCGTGATTCAGCGCGCTTCTTCGCTGGTCTTTGACAGTGTAGAAGCCAAAAAACACGCGACACGTAATCGCGCCAATGGAGAGTTGTTCTATGGACGGCGCGGAACGTTAACCCATTTCTCCTTACAACAAGCGATGTGTGAACTGGAAGGTGGCGCAGGCTGCGTGCT‑ATTTCCCTGCGGGTCGGCAGCGGTTG‑CTAATTCCATTCTTGCTTTTATCGAACAGGGCGATCATGTGTTGATGACCa                                                                                                                                                                                        >  1:371293/1‑250 (MQ=255)                                                               aGCGTGATTCAGCGCGCTTCTTCGCTGGTCTTTGACAGTGTAGAAGCCAAAAAACACGCGACACGTAATCGCGCCAATGGAGAGTTGTTCTATGGACGGCGCGGAACGTTAACCCATTTCTCCTTACAACAAGCGATGTGTGAACTGGAAGGTGGCGCAGGCTGCGTGCT‑ATTTCCCTGCGGGTCGGCAGCGGTTG‑CTAATTCCATTCTTGCTTTTATCGAACAGGGCGATCATGTGTTGATGACCAacac                                                                                                                                                                                    >  1:137598/1‑251 (MQ=255)                                                                gCGTGATTCAGCGCGCTTCTTCGCTGGTCTTTGACAGTGTAGAAGCCAAAAAACACGCGACACGTAATCGCGCCAATGGAGAGTTGTTCTATGGACGGCGCGGAACGTTAACCCATTTCTCCTTACAACAAGCGATGTGTGAACTGGAAGGTGGCGCAGGCTGCGTGCT‑ATTTCCCTGCGGGTCGGCAGCGGTTG‑CTAATTCCATTCTTGCTTTTATCGAACAGGGCGATCATGTGTTGATGACCAACAcc                                                                                                                                                                                   >  1:436022/1‑251 (MQ=255)                                                                 cGTGATTCAGAGCGTTTCTTCGCTGGTCTTTGACAGTGTAGAAGCCAAAAAACACGCGACACGTAATCGCGCCAATGGAGAGTTGTTCTATGGACGGCGCGGAACGTTAACCCATTTCTCCTTACAACAAGCGATGTGTGAACTGGAAGGTGGCGCAGGCTGCGTGCT‑ATATCCGTGCGGGTCGGCAGCGGCTG‑CTAATTCCATTATTGCTTTTATCGAACAGGACGATCATGTGTTGATGACCaaaactg                                                                                                                                                                                  >  1:379297/1‑246 (MQ=255)                                                                    gATTCAGCGCGCTTCTTCGCTGGTCTTTGACAGGGTAGAAGCCAAAAAACACGCGACACGTAATCGCGCCAATGGAGAGTTGTTCTATGGACGGCGCGGAACGTTAACCCATTTCTCCTTACAACAAGCGATGTGTGAACTGGAAGGTGGCGCAGGCTGCGTGCT‑ATTTCCCTGCGGGTCGGCAGCGGTTG‑CTAATTCCATTCTTGCTTTTATCGAACAGGGCGATCATGTGTTGATGACCAACACCGCCt                                                                                                                                                                               <  2:136106/251‑1 (MQ=255)                                                                      ttCAGCGCGCTTCTTCGCTGGTCTTTGACAGTGTAGAAGCCAAAAAACACGCGACACGTAATCGCGCCAATGGAGAGTTGTTCTATGGACGGCGCGGAACGTTAACCCATTTCTCCTTACAACAAGCGATGTGTGAACTGGAAGGTGGCGCAGGCTGCGTGCT‑ATTTCCCTGCGGGTCGGCAGCGGTTG‑CTAATTCCATTCTTGCTTTTATCGAACAGGGCGATCATGTGTTGATGACCAACACCGCCTAt                                                                                                                                                                             >  1:330527/1‑251 (MQ=255)                                                                                 tcttcGCTGGTCTTTGACAGTGTAGAATCCAAAAAACACGCGACACGTAATCGCGCCAATGGAGAGTTGTTCTATGGACGGCGCGGAACGTTAACCCATTTCTCCTTACAACAAGCGATGTGTGAACTGGAAGGTGGCGCAGGCTGCGTGCT‑ATTTCCCTGCGGGTCGGCAGCGGTTG‑CTAATTCCATTCTTGCTTTTATCGAACAGGGCGATCATGTGTTGATGACCAACACCGCCTATGAACCGAGTc                                                                                                                                                                   >  1:328194/1‑250 (MQ=255)                                                                                     cGCTGGTCTTTGACAGTGTAGAAGCCAAAAAACACGCGACACGTAATCGCGCCAGTGGAGAGATGTTCTATGGACGGCGCGGAACGTTAACCCATTTCTCATTACAACAAGCGATGTGTGAACTGGAAGGTGGCGCAGGCTGCGTGCT‑ATTTCCCTGCGGGTCGGCAGCGGTTG‑CTAATTCCATTCTTGCTTTTATCGAACAGGGCGATCATGTGTTGATGACCAACACCGCCTATGAACCGAGTCAGGAt                                                                                                                                                              >  2:415674/1‑251 (MQ=255)                                                                                                                   aaCACGCGACACGTAATCGCGCCAATGGAGAGTTGTTCTATGGACGGAGCGGAACGTTAACCCATTTCTCCTTACAACAAGCGATGTGTGACCTGGAAGGTGGCGCAGGCTGCGTGCT‑ATTTCAATGCGGGTCGGCAGCGGTTG‑CTAATTCCATTCTTGCTTTTATCGAACAGGGCGATCATGTGTTGATGACCAACACCGCCTATGAACCGAGTCAGGATTTCTGTAGCAAAATCCTCAGCAAACTGGGc                                                                                                                                >  1:317566/1‑251 (MQ=255)                                                                                                                          gACACGTAATCGCGCCAATGGAGAGTTGTTCTATGGACGGCGCGGAACGTTAACCCATTTCTCCTTACAACAAGCGATGTGAGAACTGGAAGGTGGCGCAGGCTGCGTGCT‑ATTTCCCTGCGGGTCGGCAGCGGTTG‑CTAATTCCATTCTTGCTTTTATCGACCAGGGCGATCATGTGTTGATGACCACCACCGCCTCTGAACCGAGTCAGGATTTCTGTAGCAAAATCCTCAGCAAACTGGGCGTAACg                                                                                                                          >  1:485187/1‑250 (MQ=255)                                                                                                                             acGTAATCGCGCCAATGGTGAGCTGTTCTATGGACGGCGCGGAAAGTTAACCCATTTCTCCTTACAACAAGCGATGTGTGAAATGGAAGGTGGCGCAGGCTGCGTGCT‑ATTTCCCTGCGGGTCGGCAGCGGTTG‑CTAATTCCATTCTTGCTTTTATCGAACAGGGCGATCATGTGTTGATGACCAACACCGCCTATGAACCGAGTCAGGATTTCTGAAGCAACATCCTCAGCAAACTGGGCGTAACg                                                                                                                          >  2:63663/1‑247 (MQ=255)                                                                                                                                 aaTCGCGCCAATGGAGAGTTGTTCTATGGACGGCGCGGAACGTTAACCCATTTCTCCTTACAACAAGCGATGTGTGAACTGGAAGGTGGCGCAGGCTGCGTGCT‑ATTTACCTGCGGGTCGGCAGCGGTTG‑CTAATTCCATTCTTGCTTTTATCGAACAGGGCGATCATGTGTTGATGACCCACACCGCCTATGAACCGAGTCAGGATTTCTGTAGCAAAATCCTCAGCAAACTGGGCGTAACGACATCATg                                                                                                                  >  2:357913/1‑251 (MQ=255)                                                                                                                                  aTCGCGCCAATGGAGAGTTGTTCTAGGGACGGCGCGGAACGTTAACCCATTTCTCCTTACAACAAGCGATGGGTGAACGGGAAGGGGGCGCAGGCTGCGTGCT‑ATTTCCCGGCGGGTCGGCAGCGGTTG‑CTAATTCCATTCTTGCTTTTATCGAACAGGGCGAACATGTGTTGATGACCAACACCGCCTATGAACCGAGTCAGGATTTCTGTAGCAAAATCCTCAGCAAACTGGGCGTAACGACATCATgg                                                                                                                 <  1:81086/251‑1 (MQ=255)                                                                                                                                        aaaaTGGAGATTTGTTCAAGGGACGGCTCGGAACGTTAACCCTTTTCCCCTTAAAACAAGCGTTGTGTAAGCGGGAAGGGGGCGGAGGCTGGGTGTT‑TTTCCCCTGCGGGTCGGCAG‑GGTTGGTTAATTCCATTCTTGCTTTTATCGAACGGGGCGATTTTGTGTTGATGACCAACACCGCCTATGAACCGAGTCAGGATTTCTGGAGCAAAATCCTCAGCAAACTGGGCGTAACGACATCATGGTTTGAt                                                                                                           <  2:147402/249‑1 (MQ=255)                                                                                                                                         cAATGGAGAGTTGTTCTATGGACGGCGCGGAACGTTAACCCAGTTCTCCTTACAACAAGCGATGTGTGAACTGGAAGGTGGCGCAGGCTGCGTGCT‑ATTTCCCTGCGGGTCGGCAGCGGTTG‑CTAATTCCATTCTTGCTTTTATCGAACAGGGCGATCATGTGTTGATGACCAACACCGCCTATTAACCGAGTCAGGATTTCTGTAGCAAAATCCTCAGCAAACTGGGCGTAACGACATCATGGTTTGATc                                                                                                          <  2:436022/251‑1 (MQ=255)                                                                                                                                                    tgttCTATGGACGGCGCGGAACGTTAACCCATTTCTCCTTACAACAAGCGATGTGTGAACTGGAAGGTGGCGCAGGCTGCGTGCT‑ATTTCCCTGCGGGTAGGCAGCGGTTG‑CTAATTCCATTCTTGCTTTTATCGAACAGGGCGATCATGTGTTGATGACCAACACCGCCTATGCACCGAGTCAGGCTTTCTGTAGCAAAATCCTCAGCAAACTGGGCGTAACGACATCATGGTTTGATCCGCTGATTGGt                                                                                               >  1:71101/1‑251 (MQ=255)                                                                                                                                                      ttCTATGGACGGCGCGGAACGTTAACCCATTTCTCCTTACAACAAGCGATGTGTGAACTGGAAGGTGGCGCAGGCTGCGTGCT‑ATTTCCCTGCGGGTCGGCAGCGGTTG‑CTAATTCCATTCTTGCTTTTATCGAACAGGGCGATCATGTGTTGATGACCAACACCGCCTATGAACCGAGTCAGGATTTCTGTAGCAAAATCCTCAGCAAACTGGGCGTAACGACATCATGGTTTGATCCGCTGATTGGTGc                                                                                             >  1:426229/1‑251 (MQ=255)                                                                                                                                                                         cGTTAACCCATTTCTCCTTACAACAAGCGATGTGTGAACTGGAAGGTGGCGCAGGCTGCGTGCT‑ATTTCCCTGCGGGTCGGCAGCGGTTG‑CTAATTCCATTCTTGCTTTTATCGAACAGGGCGATCATGTGTTGATGACCAACACCGCCTATGAACCGAGTCAGGATTTCTGTAGCAAAATCCTCAGCAAACTGGGCGTAACGACATCATGGTTTGATCCGCTGATTGGTGCCGATATCGTTAAGCATCTg                                                                          >  1:504986/1‑251 (MQ=255)                                                                                                                                                                            taccgcattcttcCTTACAACAGGGTATGTGTGAACTGGAAGGGGGCGCAGGCTGCGTGCA‑ATTTCCTTGCGGGTCGGCAGCGGTTG‑CTAATTCCTTTCTGGCTTTTAGCGATCCGGGCCATCATGGGTTGATGCCCAACACCGCCTATGAACCGAGTCAGGATTTCTGTAGCAAAATCCTCAGCAAACTGGGCGTAACGACATCATGGTTTGATCCGCTGATTGGTGCCGATATCGTTAAGCATCTGCa                                                                        <  1:297723/239‑1 (MQ=255)                                                                                                                                                                             aaCCCATTTCTCCTTACAACCAGCGATGTGTGAACTGGAAGGTGGCGCAGGCTGCGTGCT‑ATTTCCCTGCGGGTCGGCAGCGGTTG‑CTAATTCCATTCTTGCTTTTATCGAACAGGGCGATCATGTGTTGATGACCAACACCGCCTATGAACCGAGTCAGGATTTCTGTAGCAAAATCCTCAGCAAACTGGGCGTAACGACATCATGGTTTGATCCGCTGATTGGTGCCGATATCGTTAAGCATCTGCAGc                                                                      >  1:154794/1‑251 (MQ=255)                                                                                                                                                                              aCCCATTTCTCATTACAACAAGCGATGTGTGAACTGGAAGGTGGCGCAGTCTGCGTGCT‑ATTTCACTGCGGGTCGGCAGCGGTTG‑CTAATTCCATTCTTGCTTTTATCGAACAGGGCGATCATGTGTTGATGACCAACACCGCCTATGAACCGAGTCAGGATTTCTGTAGCAAAATCCTCAGCAAACTGGGCGTAACGACATCATGGTTTGATCCGCTGATTGGTGCCGCTGTCGTTAAGCATCTGCAGcc                                                                     >  2:57815/1‑251 (MQ=255)                                                                                                                                                                                   tttctctatacaagaggcggtgtttgAACGGGAAGGTGGCGCAGG‑TGGGGGCTTATTTCCCTGCGGGTCGGGCGGGGTTG‑CTAATTCCATTCTGGTTTTTATCAAAGAGGGAGAAAATGGTTTGATGACCAACCCCGCCTATGAACCGAGTCAGGATTTCTGTAGCAAAATCCTCAGCAAACTGGGCGTAACGACATCATGGTTTGATCCGCTTATTGGTGCCGATATCGTTAAGCATCTGCAGCCAAac                                                                 <  2:55880/226‑1 (MQ=255)                                                                                                                                                                                           tACAACAAGCGATGTGTGAACTGGAAGGTGGCGCAGGCTGCGTGCT‑ATTTCCCTGCGGGTCGGCAGCGGTTG‑CTAATTCCATTCTTGCTTTTATCGAACAGGGCGATCATGTGTTGATGACCAACACCGCCTATGAACCGAGTCAGGATTTCTGTAGCAAAATCCTCAGCAAACTGGGCGTAACGACATCATGGTTTGATCCGCTGATTGGTGCCGATATCGTTAAGCATCTGCAGCCAAACACTAAAATc                                                        >  1:196603/1‑251 (MQ=255)                                                                                                                                                                                                     gATGTGGGAACTGGAAGGAGGCGCAGGCTGCGTGCT‑ATTTCACTGCGGGTCGGCAGCGGTTG‑CTAATCCCATTCTTGCTTTTATCGAACAGGTCGATCATGGGTTGAGGACCAACACCGCCTATGAACCGAGTCAGGATTTCTGTAGCAAAATCCTCAGAAAACTGGGCGTAAAGACATCATGGTTTGATCCGCTGATAGGTGCCGATCTCGTTAAGCATCTGCAGCCAAAAACTAAATTCGTGATTCTgg                                              >  2:392347/1‑251 (MQ=255)                                                                                                                                                                                                              cggggAAGGTGCCGCATGTTGCGTGCT‑ATTTCCCTGCGGGTCGGCAGCGTTTG‑CTAATTCCATGCTTGCTTTTATCGAACGGGGCGATCTTGTTTTTATGCCCAACACCGCCTATGAACCGAGTCAGGATTTCTGTAGCAAAATCCTCAGCAAACTGGGCGTAACGACATCATGGTTTGATCCGCTGATTGGTGCCGATATCGTGAAGCATCTGCAGCCAAACACTAAAATCGTGTTTCTGGAATCGCCAg                                     <  1:79566/248‑1 (MQ=255)                                                                                                                                                                                                               cTGGAAGGTGGCGCAGGCTGCGTGCT‑ATTTCCATGCGGGTCGGCAGCGGTTG‑CTAATTCCATTCTTGCTTTTATCGAACAGGGCGATCATGTGTTGATGACCAACACCGCCTATGAACCGAGTCAGGATTTCTGTAGCAAAATCCTCAGCAACCTGGGCGTAACGACATCATGGTTTGATCCGCTGATTGTTGCCGATATCGTTAAGCATCTGCAGCCAAACACTAAAATCGCGTTTCTGGAATCGCCAgg                                    >  2:101035/1‑251 (MQ=255)                                                                                                                                                                                                                                                   cGGGTCGGCAGCGGTTG‑CTAATTCCATTCTTGCTTTTATCTAACAGAGCGATCATGTGTTGATGACCAACACCGCCTATGAACCGAGTCAGGATTTCTGTAGCAAAATCCTCAGCAAACTGGGCGTAACGACATCATGGTTTGCTCCGCTGTTTGGTGCCGATATCGTTAAGCATCTGAAGCCAAACACTAAAATCGTGTTTCTGGAATCGCCAGGCTCCATCACCATGGAAGTCCACGCCGTTCCGGCg  >  1:312309/1‑250 (MQ=255)                                                                                                                                                                                                                                                       |                                                                                                                                                                                                                                                       AGCTTGATACTCAACTGGTGAATGCAGGACGCAGCAAAAAATACACTCTCGGCGCGGTAAATAGCGTGATTCAGCGCGCTTCTTCGCTGGTCTTTGACAGTGTAGAAGCCAAAAAACACGCGACACGTAATCGCGCCAATGGAGAGTTGTTCTATGGACGGCGCGGAACGTTAACCCATTTCTCCTTACAACAAGCGATGTGTGAACTGGAAGGTGGCGCAGGCTGCGTGCT‑ATTTCCCTGCGGGGCGGCAGCGGTTG‑CTAATTCCATTCTTGCTTTTATCGAACAGGGCGATCATGTGTTGATGACCAACACCGCCTATGAACCGAGTCAGGATTTCTGTAGCAAAATCCTCAGCAAACTGGGCGTAACGACATCATGGTTTGATCCGCTGATTGGTGCCGATATCGTTAAGCATCTGCAGCCAAACACTAAAATCGTGTTTCTGGAATCGCCAGGCTCCATCACCATGGAAGTCCACGACGTTCCGGCG  >  NC_000913/3152249‑3152739

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