Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	NC_000913	4,223,638	T→C	intergenic (‑10/‑190)	iclR ← / → metH	DNA‑binding transcriptional repressor IclR/cobalamin‑dependent methionine synthase

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	4,223,638	0	T	C	100.0%	133.7 / NA	38	intergenic (‑10/‑190)	iclR/metH	DNA‑binding transcriptional repressor IclR/cobalamin‑dependent methionine synthase
Reads supporting (aligned to +/- strand):  ref base T (0/0);  new base C (14/24);  total (14/24)

GCCAACCTGACGCACGAAACCCTGCTGTTGCATCGTGGTTAGCAGGCGGTGGGTCGTGGAATTGGGTAACCCGGCTTGTTGCGCCAGTTCCGTGAGTGCCACACTGCCATTGGATTCGGCAATCCACTCCAGTAATTTCAGGCCACGCGTTAAAGACTGAACCTGTCCAGTCGCTGGTGCGGTGGCAACGGCGGGTTTTCTGCCGCGTTTCGCGGGAATGGGTGCGACCATG‑ACAGTCTCTTTTTTCTGTATCGTGG‑AAATCATTTTCATTTTTATTGTT‑AGCTAATGCAATAGTTACTGAACTGATCCGATGAG‑TTAATGTTG‑AACAAA‑TCTCATGTTGCGTGGTGGTCGC‑TTTTACCACAGATGCGTTTATGCCAGTATGGTTTGTTGAATTTTTATTAAATCTGGGTTGAGCGTGTCGGGAGCAAGTGTGAGCAGCAAAGTGGAACAACTGCGTGCGCAGTTAAATGAACGTATTCTGGTGC  >  NC_000913/4223398‑4223882                                                                                                                                                                                                                                                  |                                                                                                                                                                                                                                                           gCCAACCTGACGCACGAAACCCTGCTGTTGCATCGTGGTTAGCAGGCGGTGGGTCGTGGAATTGGGTAACCCGGCTTGTTGCGCCAGTTCCGTGAGTGCCACACTGCCATTGGATTCGGCAATCCACTCCAGTAATTTCAGGCCACGCGTTAAAGACTGAACCTGTCCAGTCGCTGGTGCGGTGGCAACGGCGGGTTTTCTGCCGCGTTTCGCGGGAATGGGTGCGACCATG‑ACAGTCTCCTTTTTCTGTa                                                                                                                                                                                                                                                  <  1:323276/251‑1 (MQ=255)   cAACCTGACGCACGAAACCCTGCTGTTGCATCGTGGTTAGCAGGCGGTGGGTCGTGGAATTGGGTAACCCGGCTTGTTGCGCCAGTTCCGTGAGTGCCACACTGCCATTGGATTCGGCAATCCACTCCAGTAATTTCAGGCCACGCGTTAAAGACTGAACCTGTCCAGTCGCTGGTGCGGTGGCAACGGCGGGTTTTCTGCCGCGTTTCGCGGGAATGGGTGCGACCATG‑ACAGTCTCCTTTTTCTGTATc                                                                                                                                                                                                                                                <  1:155968/251‑1 (MQ=255)   cAACCTGACGCACGAAACCCTGCTGTTGCATCGTGGTTAGCAGGCGGTGGGTCGTGGAATTGGGTAACCCGGCTTGTTGCGCCAGTTCCGTGAGTGCCACACTGCCATTGGATTCGGCAATCCACTCCAGTAATTTCAGGCCACGCGTTAAAGACTGAACCTGTCCAGTCGCTGGTGCGGTGGCAACGGCGGGTTTTCTGCCGCGTTTCGCGGGAATGGGTGCGACCATG‑ACAGTCTCCTTTTTCTGTATc                                                                                                                                                                                                                                                <  1:40517/251‑1 (MQ=255)    aaCCGGACGCACGAATCCCTGCTGTTGCATCGTGGTTAGCAGGCGGTGGGTCGTGGAATTGGGTACCCCGGCTTGTTGCGCCAGTTCCGTGAGTGCCACACTGCCATTGGATTCGGCAATCCACTCCAGTAATTTCAGGCCACGCGTTAAAGACTGAACCTGTCCAGTCGCTGGTGCGGTGGCAACGGCGGGTTTTCTGCCGCGTTTCGCGGGAATGGGTGCGACCATG‑ACAGTCTCCTTTTTCTGTATCg                                                                                                                                                                                                                                               <  1:409129/251‑1 (MQ=255)         gACGCACGAAACCCTGCTGTGGCATCGTGGTTAGCAGGCGGTGGGTCGTGGAATTGGGTAACCCGGCTTGTTGCGCCAGTTCCGTGAGTGCCACACTGCCATTGGATTCGGCAATCCACTCCAGTAATTTCAGGCCACGCGTTAAAGACTGAACCTGTCCAGTCGCTGGTGCGGTGGCAACGGCGGGTTTTCTGCCGCGTTTCGCGGGAATGGGTGCGACCATG‑ACAGTCTCCTTTTTCTGTATCGTGG‑a                                                                                                                                                                                                                                          <  2:427051/250‑1 (MQ=255)                           gTTGCATCGTGGTTAGCAGGCTGTGGGTCGTGGAATTGGGTAACCCGGCTTGTTGCGCCAGTTCCGTGAGTGCCACACTGCCATTGGATTCGGCAATCCACTCCAGTAATTTCAGGCCACGCGTTAAAGACTGAACCTGTCCAGTCGCTGGTGCGGTGGCAACGGCGGGTTTTCTGCCGCGTTTCGCGGGAATGGGTGCGACCATG‑ACAGTCTCCTTTTTCTGTATCGTGG‑AAATCATTTTCATTTTTa                                                                                                                                                                                                                         <  2:124251/249‑1 (MQ=255)                           gTTGCATCGTGGTTAGCAGGCGGTGGGTCGTGGAATTGGGTAACCCGGCTTGTTGCGCCAGTTCCGTGAGTGCCACACTGCCATTGGATTCGGCAATCCACTCCAGTAATTTCAGGCCACGCGTTAAAGACTGAACCTGTCCAGTCGCTGGTGCGGTGGCAACGGCGGGTTTTCTGCCGCGTTTCGCGGGAATGGGTGCGACCATG‑ACAGTCTCCTTTTTCTGTATCGTGG‑AAATCATTTTCATTTTTAt                                                                                                                                                                                                                        <  1:174508/250‑1 (MQ=255)                               cATCGTGGTTAGCAGGCGGTGGGTCGTGGAATTGGGTAACCCGGCTTGTTGCGCCAGTTCCGTGAGTGCCACACTGCCATTGGATTCGGCAATCCACTCCAGTAATTTCAGGCCACGCGTTAAAGACTGAACCTGTCCAGTCGCTGGTGCGGTGGCAACGGCGGGTTTTCTGCCGCGTTTCGCGGGAATGGGTGCGTCCATG‑ACAGTCTCCTTTTTCTGTATCGTGG‑AAATCAGTTTCATTTTTATTGTT‑a                                                                                                                                                                                                                  <  2:68667/251‑1 (MQ=255)                                      gTTAGCAGGCGGTGGGTCGTGGAATTGGGTAACCCGGCTTGTTGCGCCAGTTCCGTGAGTGCCACACTGCCATTGGATTCGGCAATCCACTCCAGTAATTTCAGGCCACGCTTTAAAGACTGAACCTGTCCAGTCTCTGGTGCGGTGGCAACGGCGGGTTTTCTGCCGCGTTTCGCGGGAATTGGTTTGACCATG‑ACAGTCTCCTTTTTCTGTATCGTTG‑AAATCATTTTCATTTTTATTGTT‑ATCTAATg                                                                                                                                                                                                           <  2:390715/251‑1 (MQ=255)                                        tAGCAGGCGGTGGGTCGTGGAATTGGGTAACCCGGCTTGTTGCGCCAGTTCCGTGAGTGCCACACTGCCATTGGATTCGGCAATCCACTCCAGTAATTTCAGGCCACGCGTTAAAGACTGAACCTGTCCAGTCGCTGGTGCGGTGGCAACGGCGGGTTTTCTGCCGCGTTTCGCGGGAATGGGTGCGACCATG‑ACAGTCTCCTTTTTCTGTATCGTGG‑AAATCATTTTCATTTTTATTGTG‑AGCTAATGc                                                                                                                                                                                                          <  1:326856/250‑1 (MQ=255)                                        tAGCAGGCGGTGGGTCGTGGAATTGGGTAACCCGGCTTGTTGCGCCAGTTCCGTGAGTGCCACACTGCCATTGGATTCGGCAATCCACTCCAGTAATTTCAGGCCACGCGTTAAAGACTGAACCTGTCCAGTCGCTGGTGCGGTGGCAACGGCGGGTTTTCTGCCGCGTTTCGCGGGAATGGGTGCGACCATG‑ACAGTCTCCTTTTTCTGTATCGTGG‑AAAGCATTTTCATTTTTATTGTT‑AGCTAATGc                                                                                                                                                                                                          <  2:74398/250‑1 (MQ=255)                                           cAGGCGGTGGGTCGTGGAATTGGGTAACCCGGCTTGTTGCGCCAGTTCCGTGAGTGCCACACTGCCATTGGATTCGGCAATCCACTCCAGTAATTTCAGGCCACGCGTTAAAGACTGAACCTGTCCAGTCGCTGGTGCGGTGGCAACGGCGGGTTTTCTGCCGCGTTTCGCGGGAATGGGTGCGACCATG‑ACAGTCTCCTTTTTCTGTATCGTGG‑AAATCATTTTCATTTTTATTGGT‑AGCTAATGCAATa                                                                                                                                                                                                      <  1:127054/251‑1 (MQ=255)                                               cGGTGGGTCGTGGAATTGGGTAACCCGGCTTGTTGCGCCAGTTCCGTGAGTGCCACACTGCCATTGGATTCTGCAATCCACTCCAGTAATTTCAGGCCACGCGTTAAAGACTGAACCTGTCCAGTCGCTGGTGCGGTGGCAACGGCGGGTTTTCTGCCGCGTTTCGCGGGAATGGGTGCGACCATG‑ACAGTCTCCTTTTTCTGTATCGTGG‑AAATCATTTTCATTTTTATTGTT‑AGCTAATGCAATAGTTa                                                                                                                                                                                                  <  1:162991/251‑1 (MQ=255)                                               cGGTGGGTCGTGGAATTGGGTAACCCGGCTTGTTGCGCCAGTTCCGTGAGTGCCACACTGCCATTGGATTCGGCAATCCACTCCAGTAATTTCAGGCCACGCGTTAAAGACTGAACCTGTCCAGTCGCTGGTGCGGTGGCAACGGCGGGTTTTCTGCCGCGTTTCGCGGGAATGGGTGCGACCATG‑ACAGTCTCCTTTTTCTGTATCGTGG‑AAATCATTTTCATTTTTATTGTT‑AGCTAATGCAATAGTTa                                                                                                                                                                                                  <  1:310240/251‑1 (MQ=255)                                                ggTGGGTCGTGGAATTGGGTAACCCGGCTTGTTGCGCCAGTTCCGTGAGTGCCACACTGCCATTGGATTCGGCAATCCACTCCAGTAATTTCAGGCCACGCGTTAAAGACTGAACCTGTCCAGTCGCTGGTGCGGTGGCAACGGCGGCATTTATGGCGCGGTTCTCGGGGGTGGGGGGGACCCTGGCCAG‑ATCCTTTTTCTGTATCGTGGGACAT‑ATTTTCATTTTTTTTGTT‑AGcatattcaattatta                                                                                                                                                                                                  >  1:91254/1‑236 (MQ=255)                                                                gggTAACCCGGCGTGTTGCGCCAGTTCCGTGAGTGCCACACTGCCATTGGATTCGGCAATCCACTCCAGTAATTTCAGGCCACGCGTTAAAGACGGAACCTGTCCAGTCGCTGGTGCGGTGGCAACGGCGGGTTTTCTGCCGCGTTTCGCGGGAATGGGTGCGACCATG‑ACAGTCTCCTTTTTCTGTATCGTGG‑AAATCATTTTCATTTTTATTGTT‑AGCTAATGCAATAGTTACTGAACTGATCCGATg                                                                                                                                                                                  <  2:232839/250‑1 (MQ=255)                                                                                   gcCAGTTCCGTGAGTGCCACACTGCCATTGGATTCGGCAATCCACTCCAGTAATTTCAGGCCACGCGTTAAAGACTGAACCTGTCCAGTCGCTGGTGCGGTGGCAACGGCGGGTTTTCTGCCGCGTTTCGCGGGAATGGGTGCGACCATG‑ACAGTCTCCTTTTTCTGTATCGTGG‑AAATCATTTTCATTTTTATTGTT‑AGCTAATGCAATAGTTACTGAACTGATCCGATGAG‑TTAATGTTG‑AACAAA‑tct                                                                                                                                                           <  1:112728/251‑1 (MQ=255)                                                                                         tCCGTGAGTGCCACACTGCCATTGGATTCGGCAATCCACTCCAGTAATTTCAGGCCACGCGTTAAAGACTGAACCTGTCCAGTCGCTGGTGCGGTGGCAACGGCGGGTTTTCTGCCGCGTTTCGCGGGGATGGGGGCGACCATG‑ACAGTCTCCTTTTTCTGTTTCGTGG‑AAATCCTTTTCCTTTTTTTTTTT‑AGCTTATGCCATAGTTACTGAACTGATCCGATGCG‑TTTATGTTG‑TACAAA‑TCTCATGtt                                                                                                                                                     >  1:30682/1‑251 (MQ=255)                                                                                           cGTGAGTGCCACACTGCCATTGGATTCGGCAATCCACTCCAGTAATTTCAGGCCACGCGTTAAAGACTGAACCTGTCCAGTCGCTGGTGCGGTGGCAACGGCGGGTTTTTTGCCGCGTTTTGCGGGGATGGGGGGGACCCTG‑ACAGGCTCCTTTTTCTGTATCGTGG‑AGATCATTTTCATTTTTATTGTT‑AGCTAATGTAATAGTTCCTGAACTTATCCGAAGGG‑TTAATGTTT‑AACAAA‑TCTCATGTTGc                                                                                                                                                   >  1:325659/1‑251 (MQ=255)                                                                                           cGTGAGTGCCACACTGCCATTGGATTCGGCAATCCACTCCAGTAATTTCAGGCCACGCGTTAAAGACTGAACCTGTCCAGTCGCTGGTGCGGTGGCAACGGCGGGTTTTTTGCCGCGGTTCGCGGGGATGGGGGGGAACATG‑ACAGGCTCCTTTTTTTGTATCGTGG‑AAATCATTTTCATTTTTATTGTT‑TGCTAAATCCATTGTTTCTGAACTGTTCCGATTAG‑TTAAAGTTG‑AACAAA‑tatcctttttc                                                                                                                                                   >  1:180067/1‑241 (MQ=255)                                                                                            gTGAGTGCCAAGCTGACATTGGATTCGGCAATCCACTCCAGTAATTTCAGGCCACGCGTTAAAGACTGCACCTGTCCAGTCGCTGGTGAGGTGGCAACGGCGGGGGTTCTTCCCCGCTTCTCGGGGACAGGGGCGCCCAAG‑AAAGCCGCCTTTTTCTGTCTCGGGG‑GGATCATTTTCATTTTTTTTGTTTGGCTC‑TGCAATTGGTACTGAACTGCTCCGTAGAG‑TTTATGTTGTAACaatatcctggtggg                                                                                                                                                    >  2:28734/1‑239 (MQ=255)                                                                                                    cacaCTGTCATTGGATTCGGCATTCCTCTCCAGTAATTACAGGCCACGCGTTAAAGACTGAACCTGTCCAGTCGCTGGTGCGGTGGCAACGGCGGGTTTTCTGCCGCGTTTCGCGGGAATGGGTGGGACCATG‑ACAGTCTCCTTTTTCTGTATCGTGG‑AAAGCATTTTCATTTTTATTGTT‑AGCTAATGCAATAGTTACTGAACTGATCCGATGAG‑TTAATGTTG‑AACAAA‑TCTCATGTTGCGtggtgggc                                                                                                                                          <  1:441123/251‑3 (MQ=255)                                                                                                                            ccACTCCAGTAATTTCAAGCCACGCGTTAAAGACTGAACCTGTCCAGTCGCTGGTGCGGTGGCAACGGCGGGGTTTTTGCCGCGTTTCGCGGGGAAGGGAGCGGCCATG‑ACAGTCTTCTTTTTCTTTGTCGGGG‑AAATCCTTTTTATTTTTATTTTT‑AGCTTATGCAATAGTTAATGAACTGATCCGATGAGATTAA‑GTTT‑GACAAA‑TCTTCTGTTTTGGGGTGGGCGC‑TTTTACCCCACATGCGTTTATg                                                                                                                 >  1:233749/1‑251 (MQ=255)                                                                                                                                 ccAGTAATTTCAGGCCACGCGTTAAAGACTGAACCTGTCCAGTCGCTGGTGCGGTGGCAACGGCGGGTTTTCTGCCGCGTTTCGCGGGAATGGGTGCGACCATG‑ACAGTCTCCTTTTTCTGTATCGTGG‑AAATCATTTTCATTTTTATTGTT‑AGCTAATGCAATAGTTACTGAACTGATCCGATGAG‑TTAATGTTG‑AACAAA‑GCTCATGTTGCGTGGTGGTCGC‑TTTTACCACAGATGCGTTTATGCCAGt                                                                                                            <  1:153364/251‑1 (MQ=255)                                                                                                                                          tCAGGCCACGCGTTAAAGACTGAACCTGTCCAGTCGCTGGTGCGGTGGCAACGGCGGGTTTTTTGCCGCGGTTCGCGGGAATTGGTGCCACCCTGGAC‑GTCTCCTTTTTCTTTATCGTGG‑GAATCATTTTTATTTTTTTTGTT‑TGCTAATGCAATAGTTACTGAACTGATCCGATGGG‑GTAATTTTG‑AACCAA‑TATCATGTTGCGTGGTGGTCGC‑TTTTTCCACAGAGGCGCTTTTGGCCGTATGGTttgt                                                                                                   >  1:322286/1‑251 (MQ=255)                                                                                                                                          tCAGGCCACGCGTTAAAGACTGAACCTGTCCAGTCGCTGGTGCGGTGGCAACGGCGGGTTTTCTGCCGCGTTTCGCGGGAATGGGTGCGACCATG‑ACAGTCTCCTTTTTCTGTATCGTGG‑AAATCATTTTCATTTTTATTGTT‑AGCTAATGCAATAGTTACTGAACTGATCCGATGAG‑TTAATGTTG‑AACAAA‑TCTCATGTTGCGTGGTGGTCGC‑TTTTACCACAGATGCGTTTATGCCAGTATGGTttg                                                                                                    <  2:30682/250‑1 (MQ=255)                                                                                                                                            aGGCAACGCGTTAAAGACTGAACCTGTCCAGTCGCTGGTGCGGTGGCAACGGCGGGGTTTCTGCCGAGGGTAGCGGGGAAGGGGGCGACCAAG‑ACAGGCTCCTTTTTCTGTATTGGGG‑GAATCATTTTCATTTTTATTGTT‑AGCTAATGCAAAAGTTACTTAAATGATCCGGAGGG‑TTAATTTTG‑GACAAAATCT‑ATGTGGGGGGGGGGGCGCCTTTTAC‑ACAGAAGGGTTTATGCCAGGATGGgttgttt                                                                                                 >  2:379697/1‑250 (MQ=255)                                                                                                                                                                    tGTCCAGTCGCTGGTGCGGTGGCAACGGCGGGTTTTTTGCCGCGTTTCGCGGGAATGGGGGCGACCATG‑ACAGTCTCCTTTTTCTGTATCGTGG‑GAATCATTTTTCTTTTTTTTTTT‑AGCTAATGCAATAGTTACTGAACTGATCCGATGAG‑TTAATGTTG‑AAAAAA‑TCTCCTGTTGCGTGGTGGTCGC‑TTTTACCACAGATGCGTTTATGCCCGTATGGTttg                                                                                                    >  1:267650/1‑224 (MQ=255)                                                                                                                                                                    tGTCCAGTCGCTGGTGCGGTGGCAACGGCGGGTTTTCTGCCGCGTTTCGCGGGAATGGGTGCGACCATG‑ACAGTCTCCTTTTTCTGTATCGTGG‑AAATCATTTTCATTTTTATTGTT‑AGCTAATGCAATAGTTACTGAACTGATCCGATGAG‑TTAATGTTG‑AACAAA‑TCTCATGTTGCGTGGTGGTCGC‑TTTTACCACAGATGCGTTTATGCCAGTATGGTttg                                                                                                    <  2:267650/224‑1 (MQ=255)                                                                                                                                                                                                cggGTTTTCTGCCGAGTTTCGCGGGAATGGGAGAGACCATG‑ACAGTCTCCTTTTTCTGTATCGTGG‑AAATCATTTTCATTTTTATTGTT‑AGCTAATGCAATAGTTACTGAACTGATCCGATGAG‑TTAATGTTG‑AACAAA‑TCTCATGTTGCGTGGTGGTCGC‑TTTTACCACAGATGCGTTTATGCCAGTATGGTTTGTTGAATTTTTATTAAATCTGGGTTGAGCGTGTCGGGAGCAAGTGTGAGCAGCaa                                              >  1:125148/1‑250 (MQ=255)                                                                                                                                                                                                 gggTTTTCTGCCGCGTTTCGCGGGAATGGGTGCGACCATG‑ACAGTCTCCTTTTTCTGTATCGTGG‑AAATCATTTTCATTTTTATTGTT‑AGCTAATGCAATAGTTACTGAACTGATCCGATGAG‑TTAATGTTG‑AACAAA‑TCTCATGTTGCGTGGTGGTCGC‑TTTTACCACAGATGCGTTTATGCCAGTATGGTTTGTTGAATTTTTATTAAATCTGGGTTGAGCGTGTCGGGAGCAAGTGTGAGCAGCAAAg                                            <  1:176341/251‑1 (MQ=255)                                                                                                                                                                                                  ggTTTTCTGCCGCGTTTCGCGGGAATGGGTGCGACCATG‑ACAGTCTCCTTTTTCTGTATCGTGG‑AAATCATTTTCATTTTTATTGTT‑AGCTAATGCAATAGTTACTGAACTGATCCGATGAG‑TTAATGTTG‑AACAAA‑TCTCATGTTGCGTGGTGGTCGC‑TTTTACCACAGATGCGTTTATGCCAGTATGGTTTGTTGAATTTTTATTAAATCTGGGTTGAGCGTGTCGGGAGCAAGTGGGAGCAGCAAAGt                                           >  1:283596/1‑251 (MQ=255)                                                                                                                                                                                                    ttttCTGCCGCGTTTCGCGGGAATGGGTGCGACCATG‑ACAGTCTCCTTTTTCTGTATCGTGG‑AAATCATTTTCATTTTTATTGTT‑AGCTAATGCAATAGTTACTGAACTGATCCGATGAG‑TTAATGtt                                                                                                                                                                       <  1:269555/128‑1 (MQ=255)                                                                                                                                                                                                    ttttCTGCAGCGTTTCGCGGGAATGGGTGCGACCATG‑ACAGTCTCCTTTTTCTGTATCGTGG‑AAATCATTTTCATTTTTATTGTT‑AGCTAATGCAATAGTTACTGAACTGATCCGATGAG‑TTAATGtt                                                                                                                                                                       >  2:269555/1‑128 (MQ=255)                                                                                                                                                                                                                                  gACCATG‑ACAGTCTCCTTTTTATGTATCGTGG‑AAATCCTTTTCATTTTTATTGTT‑AGCTAATGCAATAGTTACTGAACTGATCCGATGAG‑TTAATGTTG‑AACAAA‑TCTCATGTTGCGTGGTGGTCGC‑TTTTACCACAGATGCGTTTATGCCAGTATGGTTTGTTGAATTTTTATTAAATCTGGGTTGAGCGTGTCGGGAGCAAGTGTGAGCAGCAAAGTGGAACAACTGCGTGCGCAGTTAAATGAACGt            >  2:77566/1‑250 (MQ=255)                                                                                                                                                                                                                                    ccATG‑ACAGTCTCCTTTTTCTGTATCGTGG‑AAATCATTGTCATTTTTATTGTT‑AGCTAATGCAATAGTTACTGAACTGATCCGATGAG‑TTAATGTTG‑AACAAA‑TCTCATGTTGCGTGGTGGTCGC‑TTTTACCACAGATGCGTTTATGCCAGTATGGTTTGTTGAATTTTTATTAAATCTGGGTTGAGCGTGTCGGGAGCAAGTGTGAGCAGCAAAGTGGAACAACTGCGTGCGCAGTTAAATGAACGTAtt         <  2:343653/251‑1 (MQ=255)                                                                                                                                                                                                                                        g‑ACAGTCTCCTTTTTCTGTATCGTGG‑AAATCATTTTCATTTTTATTGTT‑AGCTAATGCAATAGTTACTGAACTTATCCGATGAG‑TTAATGTTG‑AACAAA‑TCTCATGTTGCGTGGTTGTCGC‑TTTTACCACAGATGCGTTTATGCCAGTATGGTTTGTGGAATTTTTATTAAATCTGGGTTGAGCGTGTCGGGAGCAAGTTTGAGCAGCAAAGTGGAACAACTGCGTGCGCAGTTAAATGAACGTATTCTgg     <  1:189291/251‑1 (MQ=255)                                                                                                                                                                                                                                            aGTCTCCTTTTTCTGTATCGTGG‑AAATCATTTTCATTTTTATTGTT‑AGCTAATGCAATAGTTACTGAACTGATCCGATGAG‑TTAATGTTG‑AACAAA‑TCTCATGTTGCGTGGTGGTCGC‑TTTTACCACAGATGCGTTTATGCCAGTATGGTTTGTTGAATTTTTATTAAATCTGGGTTGAGCGTGTCGGGAGCAAGTGTGAGCAGCAAAGTGGAACAACTGCGTGCGCAGTTAATTGAACGTATTCTGGTGc  >  1:60337/1‑251 (MQ=255)                                                                                                                                                                                                                                                  |                                                                                                                                                                                                                                                           GCCAACCTGACGCACGAAACCCTGCTGTTGCATCGTGGTTAGCAGGCGGTGGGTCGTGGAATTGGGTAACCCGGCTTGTTGCGCCAGTTCCGTGAGTGCCACACTGCCATTGGATTCGGCAATCCACTCCAGTAATTTCAGGCCACGCGTTAAAGACTGAACCTGTCCAGTCGCTGGTGCGGTGGCAACGGCGGGTTTTCTGCCGCGTTTCGCGGGAATGGGTGCGACCATG‑ACAGTCTCTTTTTTCTGTATCGTGG‑AAATCATTTTCATTTTTATTGTT‑AGCTAATGCAATAGTTACTGAACTGATCCGATGAG‑TTAATGTTG‑AACAAA‑TCTCATGTTGCGTGGTGGTCGC‑TTTTACCACAGATGCGTTTATGCCAGTATGGTTTGTTGAATTTTTATTAAATCTGGGTTGAGCGTGTCGGGAGCAAGTGTGAGCAGCAAAGTGGAACAACTGCGTGCGCAGTTAAATGAACGTATTCTGGTGC  >  NC_000913/4223398‑4223882

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: ‑