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	94.0%	199.0 / ‑5.6	69	A87S (GCG→TCG)	metC	cystathionine beta‑lyase/L‑cysteine desulfhydrase/alanine racemase
Reads supporting (aligned to +/- strand):  ref base G (0/4);  new base T (40/23);  total (40/29)
Fisher's exact test for biased strand distribution p-value = 2.29e-02
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

AAAGCTTGATACTCAACTGGTGAATGCAGGACGCAGCAAAAAATACACTCTCGGCGCGGTAAATAGCGTGATTCAGCGCGCTTCTTCGCTGGTCTTTGACAGTGTAGAAGCCAAAAAAC‑ACGCGACACGTAATCGCGCCAATGGAGAGTTGTTCTATGGACGGCGCGGAACGTTAACCCATTT‑CTCCTTACAACAAGCGATGT‑GTGAACTGGAAGGTGGCGCAGGCTG‑CGTGCTATTTCCCTGCGGGGCG‑GCAGCGGTTGCTAATTCCATTCTTGCTTTTATCGAACAGGG‑CGATCATGTGTTGATGACCAACACCGCCTATGAACCGAGTCAGGATTTCTGTAGCAAAATCCTCAGCAAACTGGGCGTAACGACATCATGGTTTGATCCGCTGATTGGTGCCGATATCGTTAAGCATCTGCAGCCAAACACTAAAATCGTGTTTCTGGAATCGCCAGGCTCCATCACCATGGAAGTCCACGACGTTCCGGCGAT  >  NC_000913/3152247‑3152741                                                                                                                                                                                                                                                            |                                                                                                                                                                                                                                                          aaaGCTTGATACTCAACTGGTGAATGCAGGACGCAGCAAAAAATACACTCTCGCCGCGGTAAATAGCGTGATTCAGCGCGCTTCTTCGCTGGTCTTTGACAGTGTAGAAGCCAAAAAAC‑ACGCGTCACGTAATCGCGCCAATGGAGAGTTGTTCTATGGCCGGCGCGGATCGTTAACCCATTT‑CTCCTTACAACAAGCGGTGT‑GTGAACTGGACGGTGGCGCGGGCTG‑CGTGCTATTTCCCTGCGGgtcg‑g                                                                                                                                                                                                                                                       >  2:272050/1‑251 (MQ=255)   aGCTTGATACTCAACTGGTGAATGCAGGACGCAGCAAAAAATACACTCTCGGCGCGGTAAATAGCGTGATTCAGCGCGCTTCTTCGCTGGTCTTTGACAGTGTAGAAGCCAAAAAAC‑ACGCGACACGTAATCGCGCCAATGGAGAGTTGTTCTATGGACGGCGCGGAACGTTAACCCATTT‑CTCCTTACAACAAGCGATGT‑GTGAACTGGAAGGTGGCGCAGGCTG‑CGTGCTATTTCCCTGCGGgtcg‑gc                                                                                                                                                                                                                                                      >  2:467221/1‑250 (MQ=255)        gATACTCAACTGGTGAATGCAGGACGCAGCAAAAAATACACTCTCGGCGCGGTAAATAGCGTGATTCAGCGCGCTTCTTCGCTGGTCTTTGACAGTGTAGAAGCCAAAAAAC‑ACGCGACACGTAATCGCGCCAATGGAGAGTTGTTCTATGGACGGCGCGGAACGTTAACCCATTT‑CTCCTTACAACAAGCGATGT‑GTGAACTGGAAGGTGGCGCAGGCTG‑CGTGCTATTTCCCTGCGGGTCG‑GCAGCGGt                                                                                                                                                                                                                                                >  1:137880/1‑251 (MQ=255)         aTACTCAACTGGTGAATGCAGGACGCAGCAAAAAATACACTCTCGGCGCGGTAAATAGCGTGATTCAGCGCGCTTCTTCGCTGGTCTTTGACAGTGTAGAAGCCAAAAAAC‑ACGCGACACGTAATCGCGCCAATGGAGAGTTGTTCTATGGACGGCGCGGAACGTTAACCCATTT‑CTCCTTACAACAAGCGATGT‑GTGAACTGGAAGGTGGCGCAGGCTG‑CGTGCTATTTCCCTGCGGGTCG‑GCAGCGGtt                                                                                                                                                                                                                                               >  1:129154/1‑251 (MQ=255)         aTACTCAACTGGTGAATGCAGGACGCAGCAAAAAATACACTCTCGGCGCGGTAAATAGCGTGATTCAGCGCGCTTCTTCGCTGGTCTTTGACAGTGTAGAAGCCAAAAAAC‑ACGCGACACGTAATCGCGCCAATGGAGAGTTGTTCTATGGACGGCGCGGAACGTTAACCCATTT‑CTCCTTACAACAAGCGATGT‑GTGAACTGGAAGGTGGCGCAGGCTG‑CGTGCTATTTCCCTGCGGGTCG‑GCAGCGGtt                                                                                                                                                                                                                                               >  1:175150/1‑251 (MQ=255)         aTACTCAACTGGTGAATGCAGGACGCAGCAAAAAATACACTCTCGGCGCGGTAAATAGCGTGATTCAGCGCGCTTCTTCGCTGGTCTTTGACAGTGTAGAAGCCAAAAAAC‑ACGCGACACGTAATCGCGCCAATGGAGAGTTGTTCTATGGACGGCGCGGAACGTTAACCCATTT‑CTCCTTACAACAAGCGATGT‑GTGAACTGGAAGGTGGCGCAGGCTG‑CGTGCTATTTCCCTGCGGGTCG‑GCAGCGGtt                                                                                                                                                                                                                                               >  2:420102/1‑251 (MQ=255)              cAACTGGTGAATGCAGGACGCAGCAAAAAATACACTCTCGGCGCGGTAAATAGCGTGATTCAGCGCTCTTCTTCGCTGGTCTTTGACAGTGTAGAAGCCAAAAAAC‑ACGCGACACGTAATCGCGCCAATGGAGAGTTGTTCTATGGACGGCGCGGAACGTTAACCCATTT‑CTCCTTACAACAAGCGATGT‑GTGAACTGGAAGGTGGCGCAGGCTG‑CGTGCTATTTCCCTGCGGGTCG‑GCAGCGGTTGCt                                                                                                                                                                                                                                            <  1:364938/249‑1 (MQ=255)                 cTGGTGAATGCAGGACGCAGCAAAAAATACACTCTCGGCGCGGTAAATAGCGTGATTCAGCGCGCTTCTTCGCTGGTCTTTGACAGTGTAGAAGCCAAAAAAC‑ACGCGACACGTAATCGCGCCAATGGAGAGTTGTTCTATGGACGGCGCGGAACGTTAACCCATTT‑CTCCTTACAACAAGCGATGT‑GTGAACTGGAGGGTGGCGCAGGCTG‑CGTGCTATTTCCCTGCGGGTCG‑GCAGCGGTTGCTAATTc                                                                                                                                                                                                                                       >  1:451677/1‑251 (MQ=255)                  tGGTGAATGCAGGACGCAGCAAAAAATACACTCTCGGCGCGGTAAATAGCGTGATTCAGCGCGCTTCTTCGCTGGTCTTTGACAGTGTAGAAGCCAAAAAAC‑ACGCGACACGTAATCGCGCCAATGGAGAGTTGTTCTATGGACGGCGCGGAACGTTAACCCATTT‑CTCCTTACAACAAGCGATGT‑GTGAACTGGAAGGTGGCGCCGGCTG‑CGTGCTATTTCCCTGCGGGTCG‑GCAGCGGGTGCTAATTcc                                                                                                                                                                                                                                      >  2:5719/1‑251 (MQ=255)                   ggTGAATGCAGGACGCAGCAAAAAATACACTCTCGGCGCGGTAAATAGCGTGATTCAGCGCGCTTCTTCGCTGGTCTTTGACAGTGTAGAAGCCAAAAAAC‑ACGCGACACGTAATCGCGCCAATGGAGAGTTGTTCTATGGACGGCGCGGAACGTTAACCCATTT‑CTCCTTACAACAAGCGATGT‑GTGAACTGGAAGGTGGCGCAGGCTG‑CGTGCTATTTCCCTGCGGGTCG‑GCAGCGGTTGCTAATTCCa                                                                                                                                                                                                                                     <  1:181459/251‑1 (MQ=255)                   ggTGAATGCAGGACGCAGCAAAAAATACACTCTCGGCGCGGTAAATAGCGTGATTCAGCGCGCTTCTTCGCTGGTCTTTGACAGGGTAGAAGCCAAAAAAC‑ACGCGACACGTAATCGCGCCAATGGAGAGTTGTTCTAGGGACGGCGCGGAACGTTAACCCATTT‑CTCCTTACAACAAGCGATGT‑GTGAACTGGAAGGTGGCGCAGGCTG‑CGTTCTATTTCCCTGCGGGTCG‑GCAGCGGTTGCTAATTc                                                                                                                                                                                                                                       <  2:309097/249‑1 (MQ=255)                   ggTGAAGGCAGGACGCAGCAAAAAATACACTCGCGGCGCGGTAAATAGCGTGATTCAGCGCGCTTCTTCGCTGGTCTTTGACAGTGGAGAAGCCAAAAAAC‑ACGCGACACGTAATCGCGCCAATGGAGAGTTGTTCTATGGACGGCGCGGAACGTTAACCCATTT‑CTCCTTACAACAAGCGATGT‑GTGAACTGGAAGGTGGCGCAGGCTG‑CGTGCTATTTCCCTGCGGGTCG‑GCAGCGGTTGCTAATTCCa                                                                                                                                                                                                                                     <  2:186060/251‑1 (MQ=255)                           cAGGACGCAGCAAAAAATACACTCTCGGCGCGGTAAATAGCGTGATTCAGCGCGCTTCTTCGCTGGTCTTTGACAGTGTAGAAGCCAAAAAAC‑ACGCGACACGTAATCGCGCCAATGGAGAGTTGTTCTATGGTCGGCGCGGAACGTTAACCCATTT‑CTCCTTACAACAAGCGATGT‑GTGAACTGGAAGGTGGCGCAGGCTG‑CGTGCTATTGCCCTGCTGGTCG‑GCAGCGGTTGCTAATTCCATTCTTGCt                                                                                                                                                                                                                             <  1:282240/251‑1 (MQ=255)                             ggACGCAGCAAAAAATACACTCTCGGCGCGGTAAATAGCGTGATTCAGCGCGCTTCTTCGCTGGTCTTTGACAGTGTAGAAGCCAAAAAAC‑ACGCGACACGTAATCGCGCCAATGGAGAGTTGTTCTATGGACGGCGCGGAACGTTAACCCATTT‑CTCCTTACAACAAGCGATGT‑GTGAACTGGAAGGTGGCGCAGGCTG‑CGTGCTATTTCCCTGCGGGTCG‑GCAGCGGTTGCTAATTCCATTCTTGCttt                                                                                                                                                                                                                           <  1:110824/251‑1 (MQ=255)                              gACGCAGCAAAAAATACACTCTCGGCGCGGTAAATAGCGTGATTCAGCGCGCTTCTTCGCTGGTCTTTGACAGTGTAGAAGCCAAAAAAC‑ACGCGACACGTAATCGCGCCAATGGAGAGTTGTTCTATGGACGGCGCGGAACGTTAACCCATTT‑CTCCTTACAACAAGCGATGT‑GTGAACTGGAAGGTGGCGCAGGCTG‑CGTGCTATTTCCCTGCGGGTCG‑GCAGCGGTTGCTAATTCCATTCTTGCtttt                                                                                                                                                                                                                          <  1:269242/251‑1 (MQ=255)                                      aaaaaaTACACTCTCGGCGCGGTAAATAGCGTGATTCAGCGCGCTTCTTCGCTGGTCTTTGACAGTGTAGAAGCCAAAAAAC‑ACGCGACACGTAATCGCGCCAATGGAGAGTTGTTCTATGGACGGCGCGGAACGTTAACCCATTT‑CTCCTTACAACAAGCGATGT‑GTGAACTGGAAGGTGGCGCAGGCTG‑CGTGCTATTTCCCTGCGGGTCG‑GCAGCGGTTGCTAATTCCATTCTTGCTTTTATCGAAc                                                                                                                                                                                                                   >  1:211776/1‑250 (MQ=255)                                       aaaaaTACACTCTCGGCGCGGTAAATAGCGTGATTCAGCGCGCTTCTTCGCTGGTCTTTGACAGTGTAGAAGCCAAAAAAC‑ACGCGACACGTAATCGCGCCAATGGAGAGTTGTTCTATGGACGGCGCGGAACGTTAACCCATTT‑CTCCTTACAACAAGCGATGT‑GTGAACTGGAAGGTGGCGCAGGCTG‑CGTGCTATTTCCCTGCGGGTCG‑GCAGCGGTTGCTAATTCCATTCTTGCTTTTATCGAAc                                                                                                                                                                                                                   >  2:12822/1‑249 (MQ=255)                                          aaTACACTCTCGGCGCGGTAAATAGCGTGATTCAGCGCGCTTCTTCGCTGGTCTTTGACAGTGTAGAAGCCAAAAAAC‑ACGCGACACGTAATCGCGCCAATGGAGAGTTGTTCTATGGACGGCGCGGAACGTTAACCCATTT‑CTCCTTACAACAAGCGATGT‑GTGAACTGGAAGGTGGCGCCGGCTG‑CGTGCTATTTCCCTGCGGGTCG‑GCAGCGGTTGCTAATTCCATTCTTGCTTTTATCGAACAGGG‑c                                                                                                                                                                                                             >  1:461557/1‑251 (MQ=255)                                           aTACAATCTCGGCGCGGTAAATAGCGTGATTCAGCGCGCTTCTTCGCTGGTCTTTGACAGTGTAGAAGCCAAAAAAC‑ACGCGACACGTAATCGCGCCAATGGAGAGTTGTTCTATGGACGGCGCGGAACGTTAACCCATTT‑CTCCTTACAACAAGCGATGT‑GTGAACTGGAAGGTGGCGCAGGCTG‑CGTGCTATTTCCCTGCGGGTCG‑GCAGCGGTTGCTAATTCCATTCTTGCTTTTATCGAACAGGG‑Cg                                                                                                                                                                                                            >  2:145409/1‑251 (MQ=255)                                              cacTCTCGGCGCGGTAAATAGCGTGATTCAGCGCGCTTCTTCGCTGGTCTTTGACAGTGTAGAAGCCAAAAAAC‑ACGCGACACGTAATCGCGCCAATGGAGAGTTGTTCTATGGACGGCGCGGAACGTTAACCCATTT‑CTCCTTACAACAAGCGATGT‑GTGAACTGGAAGGTGGCGCAGGCTG‑CGTGCTATTTCCCTGCGGGTCG‑GCAGCGGTTGCTAATTCCATTCTTGCTTTTATCGAACAGGG‑CGATc                                                                                                                                                                                                         >  1:392172/1‑251 (MQ=255)                                                    cGGCGCGGTAAATAGCGTGATTCAGCGCGCTTCTTCGCTGGTCTTTGACAGTGTAGAAGCCAAAAAAC‑ACGCGACACGTAATCGCGCCAATGGAGAGTTGTTCTATGGACGGCGCGGAACGTTAACCCATTT‑CTCCTTACAACAAGCGATGT‑GTGAACTGGAAGGTGGCGCAGGCTG‑CGTGCTATTTCCCTGCGGGTCG‑GCAGCGGTTGCTAATTCCATTCTTGCTTTTATCGAACAGGG‑CGATCAtgt                                                                                                                                                                                                     <  1:7335/249‑1 (MQ=255)                                                         cgGTAAATAGCGTGATTCAGCGCGCTTCTTCGCTGGTCTTTGACAGTGTAGAAGCCAAAAAAC‑ACGCGACACGTAATCGCGCCAATGGAGAGTTGTTCTATGGACGGCGCGGAACGTTAACCCATTT‑CTCCTTACAACAAGCGATGT‑GTGAACTGGAAGGTGGCGCAGGCTG‑CGTGCTATTTCCCTGCGGGTCG‑GCAGCGGTTGCTAATTCCATTCTTGCTTTTATCGAACAGGG‑CGATCATGTGTtgatg                                                                                                                                                                                              <  2:211776/251‑1 (MQ=255)                                                           gTAAATAGCGTGATTCAGCGCGCTTCTTCGCTGGTCTTTGACAGTGTAGAAGCCAAAAAAC‑ACGCGACACGTAATCGCGCCAATGGAGAGTTGTTCTATGGACGGCGCGGAACGTTAACCCATTT‑CTCCTTACAACAAGCGATGT‑GTGAACTGGAAGGTGGCGCAGGCTG‑CGTGCTATTTCCCTGCGGGTCG‑GCAGCGGTTGCTAATTCCATTCTTGCTTTTATCGAACAGGG‑CGATCATGTGTTGATGAc                                                                                                                                                                                            <  2:392776/251‑1 (MQ=255)                                                                  gCGTGATTCAGCGCGCTTCTTCGCTGGTCTTTGACAGTGTAGAAGCCAAAAAAC‑ACGCGACACGTAATCGCGCCAATGGAGAGTTGTTCTATGGACGGCGCGGAACGTTAACCCATTT‑CTCCTTACAACAAGCGATGT‑GTGAACTGGAAGGTGGCGCAGGCTG‑CGTGCTATTTCCCTGCGGGTCG‑GCAGCGGTTGCTAATTCCATTCTTGCTTTTATCGAACAGGG‑CGATCATGTGTTGATGACCAACAcc                                                                                                                                                                                     <  2:335112/251‑1 (MQ=255)                                                                  gCGTGATTCAGCGCGCTGCTTCGCTGGTCTTTCACAGTGTAGAAGCCAAAAAAC‑ACGCGACCCGTAATCGCGCCAATGGAGAGTTGTTCTATGGACGGCGCGGAACTTTAACCCATTT‑CTCCTTACAACAAGCGATGT‑GTTAACTGGAAGGTGGCTCAGGCTG‑CGTGCTATTTCCCTGCGGGTCG‑GCAGCGGTTGCTAATTCCATTCTTGCTTTTATCGAACAGGG‑CGATCATGTGTTGATGACCAACAcc                                                                                                                                                                                     <  1:118251/251‑1 (MQ=255)                                                                      taTTCAGCGCGCTTCTTCGCTGGTCTTTGACAGTGTAGAAGCCAAAAAAC‑ACGCGACACGTAATCGCGCCAATGGAGAGTTGTTCTATGGACGGCGCGGAACGTTAACCCATTT‑CTCCTTACAACAAGCGATGG‑GTGAACTGGGAGGTGGCGCAGGCTG‑CGTGCTATTTCCCGGCTGTTCG‑GCAGCGGTTGCTAATTCCATTCTTGCTTTTATCGATCAGGG‑CGATCATGTGTTGATGACCAACACCGcc                                                                                                                                                                                  <  1:433958/249‑1 (MQ=255)                                                                                gcTTCTTCGCTGGTCTTTGACAGTGTAGAAGCCAAAAAAC‑ACGCGACACGTAATCGCGCCAATGGAGAGATGTTCTATGGACGGCGCGGAACGTTAACCCATTT‑CTCCTTACAACAAGCGATGT‑GTGAACTGGAAGGTGGCGCAGGCTG‑CGTGCTATTTCCCTGCGGGTCG‑GCAGCGGTTGCTAATTCCATTCTTGCTTTTATCGAACAGGG‑CGATCATGTGTTGATGACCAACACCGCCTATGAACCg                                                                                                                                                                         >  1:230001/1‑249 (MQ=255)                                                                                        gCTGGTCTTTGACAGTGTAGAAGCCAAAAAAC‑ACGCGACACGTAATCGCGCCAATGGAGCGTTGTTCTATGGACGGCGCGGACCGTTAACCCAATT‑CTCATTACAACAAGCGATGT‑GTGCACTGGAAGGTGGCGCAGGCTG‑CGTGCTATTTCCCTGCGGGTCG‑GCAGCGGTTGCTAATTCCATTCTTGATTTTATCGAACAGGG‑CGATCATGTGTTGATGCCCAACACCGCCTATGAACCGAATCAggctt                                                                                                                                                               >  2:411383/1‑248 (MQ=255)                                                                                         cccgacctggACAGTGAAGAA‑CCAAAAAACCACGCGAACCTTAATCGCGCCAAGGGAGAGTTGTTCTAGGGAGGGCGCGGAACGTTA‑CCCATTTTCCCCTCACAACAAGCGATGT‑GGAAACGGGAGGGTGGCGCGGGCTG‑GGGGTTATTCCCCGGCGGGTGG‑GCAGCGGTTGCTAATTCCATTCTTGCTTTTACCGAACGGGG‑CGAGCATGTGTGGGTGCCAAACGCCGCCTATGAACCGAGTCAGGAt                                                                                                                                                                <  2:51354/240‑1 (MQ=255)                                                                                           ggTCTTTGACAGTGTAGAAGCCAAAAAAC‑ACGCGACACGTAATCGCGCCAATGGAGAGTTGTTCTATGGACGGCGCGGAACGTTAACCCATTT‑CTCCTTACAACAAGCGATGT‑GTGAACTGGAAGGTGGCGCAGGCTG‑CGTGCTATTTCCCTGCGGGTCG‑GCAGCGGTTGCTAATTCCATTCTTGCTTTTATCGAACAGGG‑CGATCATGTGTTGATGACCAACACCGCCTATGAACCGAGTCAGGATTTCt                                                                                                                                                            >  2:349755/1‑251 (MQ=255)                                                                                                ttGACTGTGTAGAAGCAAAAAAAC‑CCGCGACACGAAAGCGCGCCAATGGAGGTTTGTTCTATGGAGGGCGCGGAACGTTAACCCATTC‑CCCCTTACAACAAGCGTTGT‑GTGAACTGGAGGGGGGCGCAGGCTG‑CGTGCTTTTTCCCGGCGGGTCG‑GCAGGGGTTGTTAATTCCATTCTTGCTTTTAATGAACGGGG‑CTATCTTGTGTTGATGACCACCACCGCCTATGAACCGAGTCAGGATTTCTGt                                                                                                                                                          <  2:447235/248‑1 (MQ=255)                                                                                                 tGACAGTGTAGAAGCCAAAAAAC‑ACGCGACACGTAATCGCGCCAATGGAGAGTTGTTCTATGGACGGCGCGGAACGTTAACCCATTT‑CTCCTTACAACAAGCGATGT‑GTGAACTGGAAGGTGGCGCAGGCTG‑CGTGCTATTTCCCTGCGGGTCG‑GCAGCGGTTGCTAATTCCATTCTTGCTTTTATCGAACAGGG‑CGATCATGTGTTGATGACCAACACCGCCTATGAACCGAGTCAGGAt                                                                                                                                                                >  1:51354/1‑241 (MQ=255)                                                                                                      gtgtAGAAGCCAAAAAAC‑ACGCGACACGTAATCGCGCCAATGGAGAGTTGTTCTATGGACGGCGCGGAACGTTAACCCATTT‑CTCCTTACAACAAGCGATGT‑GTGAACTGGAAGGTGGCGCAGGCTG‑CGTGCTATTTCCCTGCGGGTCG‑GCAGCGGTTGCTAATTCCATTCTTGCTTTTATCGAACAGGG‑CGATCATGTGTTGATGACCAACACCGCCTATGAACCGAGTCAGGATTTCTGTAGCAAAATc                                                                                                                                                 >  1:421551/1‑251 (MQ=255)                                                                                                          aGAAGCCAAAAAAC‑ACGCGACACGTAATCGCGCCAATGGAGAGTTGTTCTATGGACGGCGCGGAACGTTAACCCATTT‑CTCATTACAACAAGCGATGT‑GTGAACTGGAAGGTGGCGCAGGCTG‑CGTGCTATTTCCCTGCGGGTCG‑GCAGCGGTTGCTAATTCCATTCTTGCTTTTATCGAACAGGG‑CGATCATGTGTTGCTGACCAACACCGCCTATGAACCGAGTCAGGATTTCTGTAGCAAAATCCTc                                                                                                                                              >  1:263012/1‑250 (MQ=255)                                                                                                                 aaaaaaC‑ACGCGACACGTAACCGCGCCAAGGGATAGTTGTTCTATGGACGGCGCGGAACTTTAACCCATTT‑CTCCTTACAACAAGCGATGT‑GTGACCGGGAAGGTGGCGCAGGCTG‑CGTGCTATTTCCCTGTGGGTCG‑GCAGCGGTTGCTAATTCCATTCTTGCTTTTATCGAGCAGGG‑GGATCATGTGTTGATGACCAACACCGCCTATGAACCGAGTCAGGATTTCTGTAGCAAAATCCTCAGCAAACt                                                                                                                                      <  1:461463/251‑1 (MQ=255)                                                                                                                  aaaaaC‑ACGCGACACGTAATCGCGCCAATGGAGAGTTGTTCTATGGACGGCGCGGAACGTTAACCCATTT‑CTCCTTACAACAAGCGATGT‑GTGAACTGGAAGGTGGCGCAGGCTG‑CGTGCTATTTCCCTGCGGGTCG‑GCAGCGGTTGCTAATTCCATTCTTGCTTTTATCGAACAGGG‑CGATCATGTGTTGATGACCAACACCGCCTATGAACCGAGTCAGGATTTCTGTAGCAAAATCCTCAGCAAACTg                                                                                                                                     >  1:77465/1‑251 (MQ=255)                                                                                                                    aaaC‑ACGCGACACGTAATCGCGCCAATGGAGAGTTGTTCTATGGACGGCGCGGAACGTTAACCCATTT‑CTCCTTACAACAAGCGATGT‑GTGAACTGGAAGGTGGCGCAGGCTG‑CGTGCTATTTCCCTGCGGGTCG‑GCAGCGGTTGCTAATTCCATTCTTGCTTTTATCGAACAGGG‑CGATCATGTGTTGATGACCAACACCGCCTATGAACCGAGTCAGGATTTCTGTAGCAAAATCCTCAGCAAACTggg                                                                                                                                   >  1:68488/1‑251 (MQ=255)                                                                                                                       acccgcgACCCTTAATCGCGCCAAGGGGGAGTTTTTCTATGGACGGGGCGGAACGTTAACCCATTT‑CCCCTTAAAACAGGCGTTGT‑GTGAACGGGAGGGGGGCGCGGGCTG‑TGTGCTATTCCCCGGCGGGCCG‑GCAGCGGTTGCTAATCCAATTCTTGCTTTTATCGAACGGGG‑CGATCATGTGTTGAGGCCCAACACCGCCTATGAACCGAGTCAGGATTTCTGTAGCAAAATCCTCAGCAAACTGGGCg                                                                                                                                 <  2:421551/248‑1 (MQ=255)                                                                                                                          cgcgACACGTAATCGAGACAATGGAGAGTTGGTCTCTGGACGGCGCGGAACGTTAACCCCTTT‑CTCATTAAAAAAAGCGATGT‑GTGAACAGCACGGTGGCGCAGGCTG‑CGTGCTATTTCCCTGCGGGTCG‑GCAGCGGTTGCTAATTCCATTCTTGCGTTTATCGAACAGGG‑CGATCATGTGTTGATGACCAACACTGCCTATGAACCGAGTCAGGATTTCTGTAGCAAAATCCTCATAAAAATGGGCGTa                                                                                                                               >  2:430396/1‑250 (MQ=255)                                                                                                                                 cGTAATCGCGCCAATGGAGAGTTGTTCTATGGACGGCGCGGAACGTTAACCCCTTT‑CTCCTTACAACAAGCGATGT‑GTGAACTGGAAGGTGGCGCAGGCTG‑CGTGCTATTTCCCTGCGGGTCG‑GCAGCGGTTGCTAATTCCATTCTTGCTTTTATCGAACAGGG‑CGATCATGTGTTGATGACCAACACCGCCTATGAACCGAGTCAGGATTTCTGTAGCAAAATCCTCAGCAAACTGGGCGTAACGAcatc                                                                                                                       >  1:410731/1‑251 (MQ=255)                                                                                                                                       cgcgCCAATGGAGAGTTGTTCTATGGACGGCGCGGAACGTTAACCCATTT‑CTCCTTACAACAAGCGATGT‑GTGAACTGGAAGGTGGCGCAGGCTG‑CGTGCTATTTCCCTGCGGGTCG‑GCAGCGGTTGCTAATTCCATTCTTGCTTTTATCGAACAGGG‑CGATCATGTGTTGATGACCAACACCGCCTATGAACCGAGTCAGGATTTCTGTAGCAAAATCCTCAGCAAACTGGGCGTAACGACATCATGGtt                                                                                                                 >  1:58840/1‑251 (MQ=255)                                                                                                                                          gcCAATGGAGAGTTGTTCTATGGACGGCGCGGAACGTTAACCCATTT‑CTCCTTACAACAAGCGATGT‑GTGAACTGGAAGGTGGCGCAGGCTG‑CGTGCTATTTCCCTGCGGGTCG‑GCAGCGGTTGCTAATTCCATCCTTGCTTTTATCGAACAGGG‑CGATCATGTGTTGATGACCAACACCGCCTATGAACCGAGTCAGGATTTCTGTAGCAAAATCCTCAGCAAACTGGGCGTAACGACATCATGGTTTg                                                                                                               >  2:417421/1‑250 (MQ=255)                                                                                                                                           ccAATGGAGAGTTGTTCTATGGACGGCGCGGAACGTTAACCCATTT‑CTCCTTACAACAAGCGATGT‑GTGAACTGGAAGGTGGCGCAGGCTG‑CGTGCTATTTCCCTGCGGGTCG‑GCAGCGGTTGCTAATTCCATTCTTGCTTTTATCGAACAGGG‑CGATCATGTGTTGATGACCAACACCGCCTATGAACCGCGTCAGGATTTCTGTAGCAAAATCCTCAGCAAACTGGGCGTAACGACATCATGGTTTGAt                                                                                                             >  2:16338/1‑251 (MQ=255)                                                                                                                                                ggAGAGTTGTTCTATGGACGGCGCGGAACGTTAACCCATTT‑CTCCTTACAACAAGCGATGT‑GTGAACTGGAAGGTGGCGCAGGCTG‑CGTGCTATTTCCCTGCGGGTCG‑GCAGCGGTTGCTAATTCCATTCTTGCTTTTATCGAACAGGG‑CGATCATGTGTTGATGACCAACACCGCCTATGAACCGAGTCAGGATTTCTGTAGCAAAATCCTCAGCAAACTGGGCGTAACGACATCATGGTTTGATCCGCt                                                                                                        >  2:403803/1‑251 (MQ=255)                                                                                                                                                        gttCTAGGGCCGGCGCGGAACTTTACCCCATTC‑CTCCTTACAACAAGCGAGGT‑GTGAACGGGAAGGTGGCGAAGGCTG‑CGTTCTTTTCCCCGGCGGGGCG‑GCAGCGGTTTCTAATTCCATTCTTGCTTTTATCGAACAGGG‑CGATCATGTGTGGATGACCAACACCGCCTATGAACCGAGTCAGGATTTCTGTAGCAAAATCCTCAGCAAACTGGGCGGAACGACATCATGGTTTGATCCGCTGATTGGt                                                                                                 <  1:403803/250‑1 (MQ=255)                                                                                                                                                             aTGGACGGCGCGGAACGTTAACCCATTT‑CTCCTTACAACAAGCGATGT‑GTGAACTGGAAGGTGGCGCAGGCTG‑CGTGCTATTTCCCTGCGGGTCG‑GCAGCGGTTGCTAATTCCATTCTTGCTTTTATc                                                                                                                                                                                                                       >  1:30095/1‑128 (MQ=255)                                                                                                                                                             aTGGACGGCGCGGAACGTTAACCCATTT‑CTCCTTACAACAAGCGATGT‑GTGAACTGGAAGGTGGCGCAGGCTG‑CGTGCTATTTCCCTGCGGGTCG‑GCAGCGGTTGCTAATTCCATTCTTGCTTTTATc                                                                                                                                                                                                                       <  2:30095/128‑1 (MQ=255)                                                                                                                                                               ggccggggcgGAACTTTACCCCATTT‑CTCCTTACACCAAGCG‑TGTTGTGACCGGGAGGTTGGCGCAGGCTG‑CGTGTTATTTCCCTGCGGGTCG‑GAAGCGGTTGTAAATTCCATTCTTGCTTTTATCGACCAGGG‑CGACCATGTGTTGATGACCAACACCGCCTATGAACCGAGTCAGGATTTCTGTAGCAAAATCCTCAGCAAACTGGGCGTAACGACATCATGGTTTGATCCGCTGATTGGTGCCGatat                                                                                         <  1:123213/244‑1 (MQ=255)                                                                                                                                                                  gggggcgaaactttaccccATTT‑CCCCTTACAACAAGGGGTGT‑GGGAACGGGAAGGTGGCGCAGGCTG‑CGTGCTATTTCCCTGCGGGCCG‑GCGGCGTTTGCTAATTCCATTCTTGCTTTTATCGAACAGGG‑CGATCATGTGTTGATAACCAACACCGCCTATGAACCGAGTCAGGATTTCTGTAGCAAAATCCTCAGCAAACTGGGCGTAACGACATCATGGTTTGATCCGCTGATTGGTGCCGATATCGt                                                                                      <  2:52913/235‑1 (MQ=255)                                                                                                                                                                   ggCGCGGAACGTTACCCCATTT‑CTCCTTACAACAAGCGATGT‑GTGAACTGGAAGGTGGCGCAGGATG‑CTTGCTATTTCCCTGCGGGTCG‑GCAGCGGTTGCTAATTCCATTCTTGCTTTTATCGAACAGGG‑CGATCATGTGTTGATGACCAACACCGCCTATGAACCGAGTCCGGATTTCTGTAGCAAAATCCTCAGCAAACTGGGCGTAACGACATAATGGTTTGATCAGCTGATTGGTGCCGATATCGtt                                                                                     >  2:397746/1‑251 (MQ=255)                                                                                                                                                                    gcgcGGAACGTTAACCCATTT‑CTCCTTACAACAAGCGATGT‑GTGAACTGGAAGGTGGCGCAGGCTG‑CGTGCTATTTCCCTGCGGGTCG‑GCAGCGGTTGCTAATTCCATTCTTGCTTTTATCGAACAGGG‑CGATCATGTGTTGATGACCAACACCGCCTATGAACCGAGTCAGGATTTCTGTAGCAAAATCCTCAGCAAACTGGGCGTAACGACATCATGGTTTGATCCGCTGATTGGTGCCGATATCGtt                                                                                     >  2:111935/1‑250 (MQ=255)                                                                                                                                                                        ggAACGTTACCCCATTT‑CCCCTTACAAAAAGCGATGT‑GTGAATTGGAAGGTGGCGCGGGCTG‑CGTTCTTTTCCCCTGCGGGTCG‑GCAGCGGTTGCTAATTCCATTCTTGCTTTTATCGACCAGGG‑CGATCATGTGTTGATGCCCAACACCGCCTATGAACCGAGTCAGGATTTCTGTAGCAAAATCCTCAGCAAACTGGGCGTAACGACATCATGGTTTGATCCGCTGATTGGTGCCGATATCGTTAAGCa                                                                                <  1:467221/251‑1 (MQ=255)                                                                                                                                                                         gAACGTTAACCCATTT‑CTCCTTACAACAAGCGATGT‑GTGAACTGGAAGGTGGCGCAGGCTG‑CGTGCTATTTCCCTGCGGGTCG‑GCAGCGGTTGCTAATTCCATTCTTGCTTTTATCGAACAGGG‑CGATCATGTGTTGATGACCAACACCGCCTATGAACCGAGTCAGGATTTCTGTAGCAAAATCCTCAGCAAACTGGGCGTAACGACATCATGGTTTGATCCGCTGATTGGTGCCGATATCGTTAAGCAt                                                                               >  1:158747/1‑251 (MQ=255)                                                                                                                                                                              ttAACCCATTT‑CTCCTTACAACAAGCGATGT‑GTGAACTGGAAGGTGGCGCAGGCTG‑CGTGCTATTTCCCTGCGGGTCG‑GCAGCGGTTGCTAATTCCATTCTTGCTTTTATCGAACAGGG‑CGATCATGTGTTGATGACCAACACCGCCTATGAACCGAGTCAGGATTTCTGTAGCAAAATCCTCAGCAAACTGGGCGTAACGACATCATGGTTTGATCCGCTGATTGGTGCCGATATCGTTAAGCATCTGc                                                                           >  2:8519/1‑250 (MQ=255)                                                                                                                                                                              taaaCCCATTT‑CCCCTGACACCAAGCGATGG‑GTGAACGGGAGGGTGGCGCAGGCTT‑CGTGCTATTTCCTGGCGGGGCG‑GAAGCGGTTGCTAATTCCATTCTTGCTTTTATCGAACAGGG‑CAATCATGTTTTGATGACCAACACCGCCTATGAACCGAGTCAGGATTTCTGTAGCAAAATCCTCAGCAAACTGGGCGTAACGACATCATGGTTTGATCCGCTGATTGGTGCCGATATCGTTAAGCATCTg                                                                            <  2:77465/247‑1 (MQ=255)                                                                                                                                                                                           tcCTTACAACAAGCGATGT‑GTGAACTGGAAGGTGGCGCAGGCTG‑CGTGCTATTTCCCTGCGGGTCG‑GCAGCGGTTGCTAATTCCATTCTTGCTTTTATCGAACAGGG‑CGATCATGTGTTGATGACCAACACCGCCTATGAACCGAGTCAGGATTTCTGTAGCAAAATCCTCAGCAAACTGGGCGTAACGACATCATGGTTTGATCCGCTGATTGGTGCCGATATCGTTAAGCATCTGCAGCCAAACACTa                                                               >  1:82330/1‑250 (MQ=255)                                                                                                                                                                                            ccTTACAACAAGCGATGT‑GTGAACTGAAAGGGGGCGCTGG‑TGTCTTTCTATTTCCTTGCGGGTCG‑GCAGCGTTTTCTAATTCCATTCTTGCTTTTATCGAACAGGG‑CTATCATGTGTGGATGCCCAACTCCGCCTATGAACCGAGTCAGGATTTCTGTAGCAAAATCCTCAGCAAACTGGGCGTAACGACATCATGGTTTGATCCGCTGATTGGTGCCGATATCGTTAAGCATCTGCAGCCAAACACTaaa                                                             <  1:135217/251‑1 (MQ=255)                                                                                                                                                                                                  aacaaGCGATGT‑GTGAACTGGAAGGTGGCGCAGGCTG‑CGTGCTATTTCCCTGCGGGTCG‑GCAGCGGTTGCTAATTCCATTCTTGCTTTTATCGAACAGGG‑CGATCATGTGTTGATGACCAACACCGCCTATGAACCGAGTCAGGATTTCTGTAGCAAAATCCTCAGCAAACTGGGCGTAACGACATCATGGTTTGATCCGCTGATTGGTGCCGATATCGTTAAGCATCTGCAGCCAAACACTAAAATCgtg                                                       >  2:96849/1‑251 (MQ=255)                                                                                                                                                                                                   acaaGCGATGT‑GTGAACTGGAAGGTGGCGCAGGCTG‑CGTGCTATTTCCCTGCGGGTCG‑GCAGCGGTTGCTAATTCCATTCTTGCTTTTATCGAACCGGG‑CGATCATGTGTTGATGACCAACACCGCCTATGAACCGAGTCAGGATTTCTGTAGCAAAATCCTCAGCAAACTGGGCGTAACGACATCATGGTTTGATCCGCTGATTGGTGCCGATATCGTTAAGCATCTGCAGCCAAACACTAAAATCgtgt                                                      >  2:46400/1‑251 (MQ=255)                                                                                                                                                                                                     aaGCGATGT‑GTGAACTGGAAGGTGGCGCAGGCTG‑CGTGCTATTTCCCTGCGGGTCG‑GCAGCGGTTGCTAATTCCATTCTTGCTTTTATCGAACAGGG‑CGATCATGTGTTGATGACCAACACCGCCTATGAACCGAGTCAGGATTTCTGTAGCAAAATCCTCAGCAAACTGGGCGTAACGACATCATGGTTTGATCCGCTGATTGGTGCCGATATCGTTAAGCATCTGCAGCCAAACACTAAAATCGTGttt                                                    >  1:359626/1‑251 (MQ=255)                                                                                                                                                                                                       gCGATGT‑GTGCACGGGAAGGTGGCGCTGTCTG‑CGTGTTATTTCCCTGTGGGTCG‑CCATGGGTTGCTATTTCCATTCTTGTTTTTATCGAACAGGG‑CGATAAGGTGTTGATCACCAACACCGCCTATGAACCGAGTCAGGATTTCTGTAGCAAAATCCTCAGCAAACTGGGCGTAACGACATCATGGTTTGATCCGCTGATTGGTGCCGATATCGTTAAGCATCTGCAGCCAAACACTAAAATCGTGTTTCt                                                  <  1:78783/251‑1 (MQ=255)                                                                                                                                                                                                             t‑gtgAACTGGAAGGTGGCGCAGGCTG‑CGTGCTATTTCCCTGCGGGTCG‑GCAGCGGTTGCTAATTCCATTCTTGCTTTTATCGAACAGGG‑CGATCATGTGTTGATGACCAACACCGCCTATGAACCGAGTCAGGATTTCTGTAGCAAAATCCTCAGCAAACTGGGCGTAACGACATCATGGTTTGATCCGCTGATTGGTGCCGATATCGTTAAGCATCTGCAGCCAAACACTAAAATCGTGTTTCTGGAATc                                            >  1:69981/1‑251 (MQ=255)                                                                                                                                                                                                               gtgAACTGGAAGGTGGCGCAGGCTG‑CGTGCTATTTCCCTGCGGGTCG‑GCAGCGGTTGCTAATTCCATTCTTGCTTTTATCGAACAGGG‑CGATCATGTGTTGATGACCAACACCGCCTATGAACCGAGTCAGGATTTCTGTAGCAAAATCCTCAGCAAACTGGGCGTAACGACATCATGGTTTGATCCGCTGATTGGTGCCGATATCGTTAAGCATCTGCAGCCAAACACTAAAATCGTGTTTCTGGAATCg                                           >  1:213295/1‑251 (MQ=255)                                                                                                                                                                                                                       gAAGGTGGCGCAGGCTG‑CGTGCTATTTCCCTGCGGGTCG‑GCAGCGGTTGCTAATTCCATTCTTGCTTTTATCGAACAGGG‑CGATCATGTGTTGATGACCAACACCGCCTATGAACCGAGTCAGGATTTCTGTAGCAAAATCCTCAGCAACCTGGGCGTAACGACATCATGGTTTGATCCGCTGATTGGTGCCGATATCGTTAAGCATCTGCAGCCAACCACTAAAATCGCGTTTCTGGAATCGCCAGGCTc                                   >  2:159243/1‑251 (MQ=255)                                                                                                                                                                                                                          ggggtctcAGGCTG‑CGTTTTATTTCCCGGCGGGGCG‑GCGGCTGTTGCTAATCCCATTCTTGTTTTTATCGAACGGGG‑CTATCATGTGTTAATACCCACCGCCGCCTATGAACCGAGTCAGGATGTCTGTAGCAAAATCCTCAGCAAACTGGGCGTAACGACATCATGGTTTGATCCGCTGATTGGTGCCGAGATCGTTGAGCATCTGCAGCCACACACTAAAATCGGGTTTCTGTAATCGCCAGGCTCCAt                                <  2:458466/244‑1 (MQ=255)                                                                                                                                                                                                                            tGGCGCAGGCTG‑CGTGCTATTTCCCTGCGGGTCG‑GCAGCGGTTGCTAATTCCATTCTTGCTTTTATCGAACAGGG‑CGATCATGTGTTGATGACCAACACCGCCTATGAACCGAGTCAGGATTTCTGTAGCAAAATCCTCAGCAAACTGGGCGTAACGACATCATGGTTTGATCCGCTGATTGGTGCCGATATCGTTAAGCATCTGCAGCCAAACACTAAAATCGTGTTTCTGGAATCGCCAGGCTCCATc                               >  1:373725/1‑250 (MQ=255)                                                                                                                                                                                                                                gcAGGCGG‑CGTTTTATTTCCCTGCGGGTCG‑GCAGCGGTTGCTAATTCCATTCTTGCTTTTATCGAACAGGG‑CGATCATGTGTTGATGCCCAACACCGCCTATGAACCGAGTCAGGATTTCTGTAGCAAAATCCTCAGCAAACTGGGCGTAACGACATCATGGTTTGATCCGCTGATTGGTGCCGATATCGTTAAGCATCTGCAGCCAAACACTAAAATCGTGTTTCTGGAATCGCCAGGCTCCATCACCAt                          <  1:275703/251‑1 (MQ=255)                                                                                                                                                                                                                                      tG‑CGTGCTATTTCCCTGCGGGTCG‑GCAGCGGTTGCTAATTCCATTCTTGCTTTTATCGAACAGGG‑CGATCATGTGTTGATGACCAACACCGCCTATGAACCGAGTCAGGATTTCTGTAGCAAAATCCTCAGCAAACTGGGCGTAACGACATCATGGTTTGATCCGCTGATTGGTGCCGATATCGTTAAGCATCTGCAGCCAAACACTAAAATCGTGTTTCTGGAATCGCCAGGCTCCATCACCATGGAAGt                    >  1:335751/1‑251 (MQ=255)                                                                                                                                                                                                                                        ggttGCTATTTCCCTGCGGGTCG‑GCGGCGGTTGCTAATCCCATTCTTGCTTTTATCGAACGGGG‑CGATCTTGTGTTGATGACCAACACCGCCTATGAACCGAGTCAGGATTTCTGTAGCAAAATCCTCAGCAAACTGGGCGTAACGACATCATGGTTTGATCCGCTGATTGGTGCCGATATCGTTAAGCATCTGCAGCCAAACACTAAAATCGTGTTTCTGGAATCGCCAGGCTCCATCACCATGGAAGTc                   <  1:5719/248‑1 (MQ=255)                                                                                                                                                                                                                                            tgtATTTCCCTG‑GGGGCGTGGAGCGGTTGCAAATCCCATTCTTCCTTTTAACGAGCAGGG‑CGATCGTGTGTTGATGACAAACCCCGCCTATGAACCGAGTCAGGATTTCTGTAGCAAAATCCTCAGCAAACTGGGCGTAACGACATCATGGTTTGATCCGCTGAGTGGTGCCGATATCGTGAAGCATCTGCAGCCAAACACTAAAATCGTGTTTCTGGAATCGCCAGGCTCCATCACCATGGAAGTCCac                <  1:464246/248‑1 (MQ=255)                                                                                                                                                                                                                                                 ttCCCTGCGGGTCG‑GCAGCGGTTGCTAATTCCATTCTTGCTTTTATCGAACAGGG‑CGATCATGTGTTGATGACCAACACCGCCTATGAACCGAGTCAGGATTTCTGTAGCAAAATCCTCAGCAAACTGGGCGTAACGACATCATGGTTTGATCCGCTGATTGGTGCCGATATCGTTAAGCATCTGCAGCCAAACACTAAAATCGTGTTTCTGGAATCGCCAGGCTCCATCACCATGGAAGTCCACGACGtt          >  1:210207/1‑251 (MQ=255)                                                                                                                                                                                                                                                         ggggCG‑GCAGCGGTTGTTAATTCCTTTCTTTCTTTTAT‑GTACAGGGGCGATAATGTGTGGAGGCCCAACACCGCCGATGAACCGAGTCAGGATTTCTGTAGCAAAATCCTCAGCAAACTGGGCGTAACGACATCATGGTTTGATCCGCTGATTGGTGCCGATATCGTTAAGCATCTGCAGCCAAACACTAAAATCGTGTTTCTGGAATCGCCAGGCTCCATCACCATGGAAGTCCACGACGTTCCGGCGAt  <  2:230001/251‑1 (MQ=255)                                                                                                                                                                                                                                                            |                                                                                                                                                                                                                                                          AAAGCTTGATACTCAACTGGTGAATGCAGGACGCAGCAAAAAATACACTCTCGGCGCGGTAAATAGCGTGATTCAGCGCGCTTCTTCGCTGGTCTTTGACAGTGTAGAAGCCAAAAAAC‑ACGCGACACGTAATCGCGCCAATGGAGAGTTGTTCTATGGACGGCGCGGAACGTTAACCCATTT‑CTCCTTACAACAAGCGATGT‑GTGAACTGGAAGGTGGCGCAGGCTG‑CGTGCTATTTCCCTGCGGGGCG‑GCAGCGGTTGCTAATTCCATTCTTGCTTTTATCGAACAGGG‑CGATCATGTGTTGATGACCAACACCGCCTATGAACCGAGTCAGGATTTCTGTAGCAAAATCCTCAGCAAACTGGGCGTAACGACATCATGGTTTGATCCGCTGATTGGTGCCGATATCGTTAAGCATCTGCAGCCAAACACTAAAATCGTGTTTCTGGAATCGCCAGGCTCCATCACCATGGAAGTCCACGACGTTCCGGCGAT  >  NC_000913/3152247‑3152741

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