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%	258.2 / NA	70	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 (19/51);  total (19/51)

GTTCGCCAACCTGACGCACGAAACCCTGCTGTTGCATCGTGGTTAGCAGGCGGTGGGTCGTGGAATTGGGTAACCCGGCTTGTTGCGCCAGTTCCGTGAGTGCCACACTGCCATTGGATTCGGCAATCCACTCCAGTAATTTCAGGCCACGCGTTAAAGACTGAACCTGTCCAGTCGCTGGTGCGGTGGCAACGGCGGGTTTTCTGCCGCGTTTCGCGGGAATGGGTGCGACCATGACAGTCTCTTTTTTCTGTATCGTGGAAATCATTTTCATTTTTATTGTTAGCTAAT‑GCAATAGTTACTGAACTGATCCGATGAGTTAATGTTG‑AACAAATCTCATGTTGCGTGGTGGTCGCTTTTACCACAGATGCGTTTATGCCAGTATGGTTTGTTGAATTTTTATTAAATCTGGGTTGAGCGTGTCGGGAGCAAGTGTGAGCAGCAAAGTGGAACAACTGCGTGCGCAGTTAAATGAACGTATTCTGGTGCT  >  NC_000913/4223394‑4223883                                                                                                                                                                                                                                                     |                                                                                                                                                                                                                                                        gTTCGCCAGCCTGACGCACGAAACCCTGCTGTTGCATCGTGGTTAGCAGGCGGTGGGTCGTGGAATTGGGTAACCCGGCTTGTTGCGCCAGTTCCGTGAGTGCCACACTGCCATTGGATTCGGCAATCCACTCCAGTAATTTCAGGCCACGCATTAAAGACTGAACCTGTCCAGTGGCTGGTGCGGTGGCAACGGCGGGTTTTCTGCCGCGTTTCTCGTGAATGGGTGCGACCATGACAGTCTCCTTTTTc                                                                                                                                                                                                                                                   <  1:337861/251‑1 (MQ=255) gTTCGCCAACCTGACGCACGAAACCCTGCTGTTGCATCGTGGTTAGCAGGCGGTGGGTCGTGGAATTGGGTAACCCGGCTTGTTGCGCCAGTTCCGTGAGTGCCACACTGCCATTGGATTCGGCAATCCACTCCAGTAATTTCAGGCCACGCGTTAAAGACTGAACCTGTCCAGTCGCTGGTGCGGTGGCAACGGCGGGTTTTCTGCCGCGTTTCGCGGGAATGGGTGCGACCATGACAGTCTCCTTTTTc                                                                                                                                                                                                                                                   <  1:444550/251‑1 (MQ=255)    cGCCAACCTGACGCACGAAACCCTGCTGTTGCATCGTGGGTAGCAGGCGGTGGGTCGTGGAATTGGGTAACCCGGCTTGTTGCGCCAGTTCCGTGAGTGCCACACTGCCATTGGATTCGGCAATCCACTCCAGTAATTTCAGGCCACGCGTTAAAGACTGAACCTGTCCAGTCGCTGGTGCGGTGGCAACGGCGGGTTTTCTGCCGCGTTTCGCGGGAATGGGTGCGACCATGACAGTCTCCTTTTTCTg                                                                                                                                                                                                                                                 <  2:435758/250‑1 (MQ=255)        aaCCTGACGCACGAAACCCTGCTGTTGCATCGTGGTTAGCAGGCGGTGGGTCGTGGAATTGGGTAACCCGGCTTGTTGCGCCAGTTCCGTGAGTGCCACACTGCCATTGGATTCGGCAATCCACTCCAGTAATTTCAGGCCACGCGTTAAAGACTGAACCTGTCCAGTCGCTGGTGCGGTGGCAACGGCGGGGTTTTTGCCGCGGTTTGCGGGGAAGGGGGGGGCCATGACAGGCTCCTTTTTTTGTTTCg                                                                                                                                                                                                                                            >  1:460789/1‑251 (MQ=255)         aCCTGACGCACGAAACCCTGCTGTTGCATCGTGGTTAGCAGGCGGTGGGTCGTGGAATTGGGTAACCCGGCTTGTTGCGCCAGTTCCGTGAGTGCCACACTGCCATTGGATTCGGCAATCCACTCCAGTAATTTCAGGCCACGCGTTAAAGACTGAACCTGTCCAGTCGCTGGTGCGGTGGCAACGGCGGGTTTTCTGCCGCGTTTCGCGGGAATGGGTGCGACCATGACAGTCTCCTTTTTCTGTATCGt                                                                                                                                                                                                                                           <  2:424592/251‑1 (MQ=255)          ccTGACGCACGAAACCCTGCTGTTGCATCGTGGTTAGCAGGCGGTGGGTCGTGGAATTGGGTAACCCGGCTTGTTGCGCCAGTTCCGTGAGTGCCACACTGCCATTGGATTCGGCAATCCACTCCAGTAATTTCAGGCCACGCGTTAAAGACTGAACCTGTCCAGTCGCTGGTGCGGTGGCAACGGCGGGTTTTCTGCCGCGTTTCGCGGGAATGGGTGCGACCATGACAGTCTCCTTTTTCTGTATCGTg                                                                                                                                                                                                                                          <  1:34366/251‑1 (MQ=255)          ccTGACGCACGAAACCCTGCTGTTGCATCGTGGTTAGCAGGCGGTGGGTCGTGGAATTGGGTAACCCGGCTTGTTGCGCCAGTTCCGTGAGTGCCACACTGCCATTGGATTCGGCAATCCACTCCAGTAATTTCAGGCCACGCGTTAAAGACTGAACCTGTCCAGTCGCTGGTGCGGTGGCAACGGCGGGTTTTCTGCCGCGTTTCGCGGGAATGGGTGCGACCATGACAGTCTCCTTTTTCTGTATCGTg                                                                                                                                                                                                                                          <  1:330283/251‑1 (MQ=255)             gACGCACGAAACCCTGCTGTTGCATCGTGGTTAGCAGGCGGTGGGTCGTGGAATTGGGTAACCCGGCTTGTTGCGCCAGTTCCGTGAGTGCCACACTGCCATTGGATTCGGCAATCCACTCCAGTAATTTCAGGCCACGCGTTAAAGACTGAACCTGTCCAGTCGCTGGTGCGGTGGCAACGGCGGGTTTTCTGCCGCGTTTCGCGGGAATGGGTGCGACCATGACAGTCTCCTTTTTCTGTATCGTgg                                                                                                                                                                                                                                         <  1:12731/249‑1 (MQ=255)             gACGCACGAAACCCTGCTGTTGCATCGTGGTTAGCAGGCGGTGGGTCGTGGAATTGGGTAACCCGGCTTGTTGCGCCAGTTCCGTGAGTGCCACACTGCCATTGGATTCGGCAATCCACTCCAGTAATTTCAGGCCACGCGTTAAAGACTGAACCTGTCCAGTCGCTGGTGCGGTGGCAACGGCGGGTTTTCTGCCGCGTTTCGCGGGAATGGGTGCGACCATGACAGTCTCCTTTTTCTGTATCGTgg                                                                                                                                                                                                                                         <  2:295264/249‑1 (MQ=255)               cGCACGAAACCCTGCTGTTGCATCGTGGTTAGCAGGCGGTGGGTCGTGGAATTGGGTAACCCGGCTTGTTGCGCCAGTTCCGTGAGTGCCACACTGCCATTGGATTCGGCAATCCACTCCAGTAATTTCAGGCCACGCGTTAAAGACTGAACCTGTCCAGTCGCTGGTGCGGTGGCAACGGCGGGTTTTCTGCCGCGTTTCGCGGGAATGGGTGCGACCATGACAGTCTCCTTTTTCTGTATCGTGGAAAt                                                                                                                                                                                                                                     <  1:53397/251‑1 (MQ=255)                    gAAACCCTGCTGTTGCATCGTGGTTAGCAGGCGGTGGGTCGTGGAATTGGGTAACCCGGCTTGTTGCGCCAGTTCCGTGAGTGCCACACTGCCATTGGATTCGGCAATCCACTCCAGTAATTTCAGGCCACGCGTTAAAGACTGAACCTGTCCAGTCGCTGGTGCGGTGGCAACGGCGGGTTTTTTGCCGCGGTTTGCGGGAATGGGTGGGGCCATGACAGTCTCCTTTTTTTGTATCGTGGAAAtcattt                                                                                                                                                                                                                                >  1:15751/1‑251 (MQ=255)                         ccTGCTGTTGCATCGTGGTTAGCAGGCGGTGGGTCGTGGAATTGGGTAACCCGGCTTGTTGCGCCAGTTCCGTGAGTGCCACACTGCCATTGGATTCGGCAATCCACTCCAGTAATTTCAGGCCACGCGTTAAAGACTGAACCTGTCCAGTCGCTGGTGCGGTGGCAACGGCGGGTTTTCTGCCGCGTTTCGCGGGAATGGGTGCGACCATGACAGTCTCCTTTTTCTGTATCGTGGAAAtcattttcatt                                                                                                                                                                                                                           <  1:514481/251‑1 (MQ=255)                                 tGCATCGTGGTTAGCAGGCGGTGGGTCGTGGAATTGGGTAACCCGGCTTGTTGCGCCAGTTCCGTGAGTGCCACACTGCCATTGGATTCGGCAATCCACTCCAGTAATTTCAGGCCACGCGTTAAAGACTGAACCTGTCCAGTCGCTGGTGCGGTGGCAACGGCGGGTTTTCTGCCGCGTTTCGCGGGAATGGGTGCGACCATGACAGTCTCCTTTTTCTGTATCGTGGAAATCATTTTCATTTTTATTg                                                                                                                                                                                                                    <  1:531988/250‑1 (MQ=255)                                  gcatcggtgTTAGCAGGCGGTGGGTCGTGGAATTGGGTAACCCGGCTTGTTGCGCCAGTTCCGTGAGTGCCACACTGCCATTGGATTCGGCAATCCACTCCAGTAATTTCAGGCCACGCGTTAAAGACTGAACCTGTCCAGTCGCTGGTGCGGTGGCAACGGCGGGTTTTCTGCCGCGTTTCGCGGGAATGGGTGCGACCATGACAGTCTCCTTTTTCTGTATCGTGGAAATCATTTTCATTTTTATTGt                                                                                                                                                                                                                   <  2:423543/242‑1 (MQ=255)                                    aTCGTGGTTAGCAGGCGGTGGGTCGTGGAATTGGGTAACCCGGCTTGTTGCGCCAGTTCCGTGAGTGCCACACTGCCATTGGATTCGGCAATCCACTCCAGTAATTTCAGGCCACGCGTTAAAGACTGAACCTGTCCAGTCGCTGGTGCGGTGGCAACGGCGGGTTTTCTGCCGCGTTTCGCGGGAATGGGTGCGACCATGACAGTCTCCTTTTTCTGTATCGTGGAAATCATTTTCATTTTTATTGTTAg                                                                                                                                                                                                                <  1:566809/251‑1 (MQ=255)                                       gTGGTTAGCAGGCGGTGGGTCGTGGAATTGGGTAACCCGGCTTGTTGCGCCAGTTCCGTGAGTGCCACACTGCCATTGGATTCGGCAATCCACTCCAGTAATTTCAGGCCACGCGTTAAAGACTGAACCTGTCCAGTCGCTGGTGCGGTGGCAACGGCGGGTTTTCTGCCGCGTTTCGCGGGAATGGGTGCGACCATGACAGTCTCCTTTTTCTGTATCGTGGAAATCATTTTCATTTTTATTGTTAGCTa                                                                                                                                                                                                             <  1:511671/251‑1 (MQ=255)                                          gTTAGCAGGCGGTGGGTCGTGGAATTGGGTAACCCGGCTTGTTGCGCCAGTTCCGTGAGTGCCACACTGCCATTGGATTCGGCAATCCACTCCAGTAATTTCAGGCCACGCGTTAAAGACTGAACCTGTCCAGTCGCTGGTGCGGTGGCAACGGCGGGTTTTCTGCCGCGTTTCGCGGGAATGGGTGCGACCATGACAGTCTCCTTTTTCTGTATCGTGGAAATCATTTTCATTTTTATTGTTAGCTAAT‑g                                                                                                                                                                                                         <  1:230865/251‑1 (MQ=255)                                              gCAGGCGGTGGGTCGTGGAATTGGGTAACCCGGCTTGTTGCGCCAGTTCCGTGAGTGCCACACTGCCATTGGATTCGGCAATCCACTCCAGTAATTTCAGGCCACGCGTTAAAGACTGAACCTGTCCAGTCGCTGGTGCGGTGGCAACGGCGGGTTTTCTGCCGCGTTTCGCGGGAATGGGTGCGACCATGACAGTCTCCTTTTTCTGTATCGTGGAAATCATTTTCATTTTTATTGTTAGCTAAT‑GCAAt                                                                                                                                                                                                     <  2:392626/251‑1 (MQ=255)                                                     gTGGGTCGTGGAATTGGGTAACCCGGCTTGTTGCGCCAGTTCCGTGAGTGCCACACTGCCATTGGATTCGGCAATCCACTCCAGTAATTTCAGGCCACGCGTTAAAGACTGAACCGGTCCAGTCGCTGGTGCGGTGGCAACGGCGGGTTTTCTGCCGCGTTTCGCGGGAATGGGTGCGACCATGACAGTCTCCTTTTTCTGTATCGTGGAAATCATTTTCATTTTTATTGTTAGCTAAT‑GCAATAGTTact                                                                                                                                                                                              <  1:29898/251‑1 (MQ=255)                                                         ggcGTGGAATTGGGTAACCCGGCTTGTTGCGCCAGGTCCGTGAGTGCCACACTGCCATTGGATTCGGCAATCCAGTCCAGTAATTTCAGGCCACGCGTTAAAGACTGAACCTGTCCAGTCGCTGGTGCGGTGGTAACGGCGGGTTTTCTGCCGCGTTTCGCGGGAATGGGTGCGACCATGACAGTCTCCTGTTTCTGTATCGTTGAAATCATTTTCATTTTTAGTGTTAGCTAAA‑CCGATAGTTactgaac                                                                                                                                                                                          <  2:295169/249‑1 (MQ=255)                                                         gTCGTGGAATTGGGTAACCCGGCTTGTTGCGCCAGTTCCGTGAGTGCCACACTGCCATTGGATTCGGCAATCCACTCCAGTAATTTCAGGCCACGCGTTAAAGACTGAACCTGTCCAGTCGCTGGTGCGGTGGCAACGGCGGGTTTTCTGCCGCGTTTCGCGGGAATGGGTGCGACCATGACAGTCTCCTTTTTCTGTATCGTGGAAATCATTTTCATTTTTATTGTTAGCTAAT‑GCAATAGTTactgaac                                                                                                                                                                                          <  1:582301/251‑1 (MQ=255)                                                           cGTGGAATTGGGTAACCCGGCTTGTTGCGCCAGTTCCGTGAGTGCCACACTGCCATTGGATTCGGCAATCCACTCCAGTAATTTCAGGCCACGCGTTAAAGACTGAACCTGTCCAGTCGCTGGTGCGGTGGCAACGGCGGGTTTTCTGCCGCGTTTCGCGGGAATGGGTGCGACCATGACAGTCTCCTTTTTCTGTATCGTGGAAATCATTTTCATTTTTATTGTTAGCTAAT‑GCATTAGTTactgaactg                                                                                                                                                                                        <  1:46520/251‑1 (MQ=255)                                                           cGTGGAATTGGGTAACCCGGCTTGTTGCGCCAGTTCCGTGAGTGCCACACTGCCATTGGATTCGGCAATCCACTCCAGTAATTTCAGGCCACGCGTTAAAGACTGAACCTGTCCAGTCGCTGGTGCGGTGGCAACGGCGGGTTTTCTGCCGCGTTTCGCGGGAATGGGTGCGACCATGACAGTCTCCTTTTTCTGTATCGTGGAAATCATTTTCATTTTTATTGTTAGCTAAT‑GCAATAGTTactgaactg                                                                                                                                                                                        <  2:82209/251‑1 (MQ=255)                                                           cGTGGAATTGGGTAACCCGGCTTGTTGCGCCAGTTCCGTGAGTGCCACACTGCCATTGGATTCGGCAATCCACTCCAGTAATTTCAGGCCACGCGTTAAAGACTGAACCTGTCCAGTCGCTGGTGCGGTGGCAACGGCGGGTTTTCTGCCGCGTTTCGCGGGAATGGGTGCGACCATGACAGTCTCCTTTTTCTGTATCGTGGAAATCATTTTCATTTTTATTGTTAGCTAAT‑GCAATAGTTactgaact                                                                                                                                                                                         <  1:95443/250‑1 (MQ=255)                                                            gTGGAATTGGGTAACCCGGCTTGTTGCGCCAGTTCCGTGAGTGCCACACTGCCATTGGATTCGGCAATCCACTCCAGTAATTTCAGGCCACGCGTTAAAGACTGAACCTGTCCAGTCGCTGGTGCGGTGGCAACGGCGGGGTTTTTGCCGCGGTTCGCGGGGATGGGGGCGGCCATGACAGTCTCCTTTTTCTTTTTTGTGGGAATCATTTTTCTTTTTTTTTTTTGCTAAT‑GCAATAGTTactgaactg                                                                                                                                                                                        >  1:414210/1‑250 (MQ=255)                                                                      gTAACCCGGCTTGTTGCGCCAGTTCCGTGAGTGCCACACTGCCATTGGATTCGGCAATCCACTCCAGTAATTTCAGGCCACGCGTTAAAGACTGAACCTGTCCAGTCGCTGGTGCGGTGGCAACGGCGGGGTTTTTTCCGCGTTTTGCGGGGATGGGGGCGGCCCTGACAGTCTCCTTTTTTTGTTTCGTGGGAATCATTTTCATTTTTTTTTTTTGCTAAT‑GCAATAGTTACTTAACTGATCCGATgtgt                                                                                                                                                                             >  1:67317/1‑248 (MQ=255)                                                                         aCCCGGCTTGTTGCGCCAGTTCCGTGAGTGCCACACTGCCATTGGATTCGGCAATCCACTCCAGTAATTTCAGGCCACGCGTTAAAGACTGAACCTGTCCAGTCGCTGGTGCGGTGGCAACGGCGGGTTTTCTGCCGCGTTTCGCGGGAATGGGTGCGACCATGACAGTCTCCTTTTTCTGTATCGTGGAAATCATTTTCATTTTTATTGTTAGCTAAT‑GCAATAGTTACTGAACTGATCCGATGAGTTaa                                                                                                                                                                          <  2:364541/251‑1 (MQ=255)                                                                                        ccAGTTCCGTGATTGCCACACTGCCATTGGATTCGGCAATCCACTCCAGTAATTTCAGGCCACGCGTTAAAGACTGAACCTGTCCAGTCGCTGGTGCGGTGGCAACGGCGGGTTTTCTGCCTCGTTTCGCGGGAATGGGTGCGACCATGACAGTCTCCTTTTTCTGTATCGTGGAAATCATTTTCATTTTTATTGTTAGCTAAT‑GCAATAGTTACTGAACTGATCCGATGAGTTAATGTTG‑AACAAAtctc                                                                                                                                                          <  1:473376/251‑1 (MQ=255)                                                                                          aGTTCCGTGAGTGCCACACTGCCATTGGATTCGGCAATCCACTCCAGTAATTTCAGGCCACGCGTTAAAGACTGAACCTGTCCAGTCGCTGGTGCGGTGGCAACGGCGGGTTTTCTGCCGCGTTTCGCGGGAATGGGTGCGACCATGACAGTCTCCTTTTTCTGTATCGTGGAAATCATTTTCATTTTTATTGTTAGCTAAT‑GCAATAGTTACTGAACTGATCCGATGAGTTAATGTTG‑AACAAATCTCAt                                                                                                                                                        <  2:527717/251‑1 (MQ=255)                                                                                              ccGTGAGTGCCACACTGCCATTGGATTCGGCAATCCACTCCAGTAATTTCAGGCCACGCGTTAAAGACTGAACCTGTCCAGTCGCTGGTGCGGTGGCAACGGCGGGTTTTCTGCCGCGTTTCGCGGAAATGGGTGCGACCATGACAGTCTCCTTTTTCTGTAGCGTGGAAATCATTTTCATTTTTATTGTTAGCTAAT‑GCAATAGTTACTGAACTGATCCGATGAGTTAATGTTG‑AACAAATCTCATGTTg                                                                                                                                                    <  2:395921/251‑1 (MQ=255)                                                                                                       ccACACTGCCATTGGATTCGGCAATCCACTCCAGTAATTTCAGGCCACGCGTTAAAGACTGAACCTGTCCAGTCGCTGGTGCGGTGGCAACGGCGGGTTTTTTGCCGCGTTTCGCGGGAATGGGGGCGACCCTGACAGTCTCCTTTTTTTGTTTCGTGGGAATCATTTTTATTTTTTTTTTTAGCTAAT‑GCAATAATTACTGAACTGATCCGATGAGTTTATGTTT‑AACAAATTTCCTGTTGCGTggtggt                                                                                                                                           >  1:521679/1‑251 (MQ=255)                                                                                                       ccACACTGCCATTGGATTCGGCAATCCACTCCAGTAATTTCAGGCCACGCGTTAAAGACTGAACCTGTCCAGTCGCTGGTGCGGTGGCAACGGCGGGTTTTCTGCCGCGTTTCGCGGGAATGGGTGCGACCATGACAGTCTCCTTTTTCTGTATCGTGGAAATCATTTTCATTTTTATTGTTAGCTAAT‑GCAATAGTTACTGAACTGATCCGATGAGTTAATGTTG‑AACAAATCTCATGTTGCGTggtggt                                                                                                                                           <  2:218331/251‑1 (MQ=255)                                                                                                                cATTGGAATCGGCAATCCACTCCATTAATTTCAGGCCACGCGTGAAAGACTGAACCTGTCCAGTCGCTGGTGCGGTGGCAACGTCCGGTTTTCTGCCGCGTTTCGCAGGAATTGGTTCGACCATGACAGTCACCTTTTTCTGTATCTTGGAAATCATTTTCATTTTTATTGTTAGATAAT‑GCAATAGTTACTGAACTGATCCGATGAGTTAACGTGG‑AACAAATCTCATGTTGCGTGGTGGTCGCTTTTAc                                                                                                                                  <  2:81035/251‑1 (MQ=255)                                                                                                                  ttGGATTCGGCAATCCACTCCAGTAATTTCAGGCCACGCGTTAAAGACTGAACCTGTCCAGTCGCTGGTGCGGTGGCAACGGCGGGGTTTTTTCCGGGGTTCGCGGGGATGGGGGCGGGCAAGAGAGACTCCCTTTTTTTTTTCGTGGGAATTATTTTTTTTTTTTTTTTTTGCTAATTGCA‑TAGTTAATGAAATGATTCGATGGATTAATTTTG‑AAAAAATATTATTTTGCgtggggggtgcttttttcc                                                                                                                                >  2:546865/1‑236 (MQ=255)                                                                                                                   tGGATTCGGCAATCCACTCCAGTAATTTCAGGCCACGCGTTAAAGACTGAACCTGTCCAGTCGCTGGTGCGGTGGCAACGGCGGGTTTTCTGCCGCGTTTCGCGGGAATGGGTGCGACCATGACAGTCTCCTTTTTCTGTATCGTGGAAATCATTTTCATTTTTATTGTTAGCTAAT‑GCAATAGTTACTGAACTGATCCGATGAGTTAATGTTG‑AACAAATCTCATGTTGCGTGGTGGTCGCTTTTACca                                                                                                                                <  1:214994/250‑1 (MQ=255)                                                                                                                    ggATTCGGCAATCCACTCCAGTAATTTCAGGCCACGCGTTAAAGACTGAACCTGTCCAGTCGCTGGTGCGGTGGCAACGGCGGGTTTTCTGCCGCGTTTCGCGGGAATGGGTGCGACCATGACAGTCTCCTTTTTCTGTATCGTGGAAATCATTTTCATTTTTATTGTTAGCTAAT‑GCAATAGTTACTGAACTGATCCGATGAGTTAATGTTG‑AACAAATCTCATGTTGCGTGGTGGTCGCTTTTACcaca                                                                                                                              <  2:460789/251‑1 (MQ=255)                                                                                                                    ggATTCGGCAATCCACTCCAGTAATTTCAGGCCACGCGTTAAAGACTGAACCTGTCCAGTCGCTGGTGCGGTGGCAACGGCGGGGTTTTTGCCGCGGTTTGCGGGGATGGGTGCGGCCATGAAAGTCTCCTTTTTCTTTTTCGTGGGAAACATTTTCATTTTTATTTTTAGGTTAA‑GCAAAAGTTTCTGAAATGATCCGAAGAGTTTATGTTG‑GAAAAAACCCATGTTGCggggggggggcttttttcacc                                                                                                                              >  2:522894/1‑231 (MQ=255)                                                                                                                                 cACTCCAGTAATTTC‑GGCCACGCGTTAAAGACTGAACCTGTCCAGTCGCTGGTGCGGTGGCAACGGCGGGTTTTCTGCCGCGTTTCGCGGGAATGGGTGCGACCATGACAGTCTCCTTTTTCTGTATCGTGGAAATCATTTTCAGTTTTATTGTTAGCTAAT‑GCAATAGTTACTGAACTGATCCGATGAGTTAATGTTG‑AACAAATCTCATGTTGCGTGGTGGTCGCTTTTACCACAGATGCGTTTATGcc                                                                                                                <  1:12235/251‑1 (MQ=255)                                                                                                                                    tCCAGTAATTTCAGGCCACGCGTTAAAGACTGAACCTGTCCAGTCGCTGGTGCGGTGGCAACGGCGGGGTTTTTGCCGCGTTTCGCGGGGATTGGTGGGGCCCTGACAGTCTCCTTTTTTTGTATCGTGGGAATCATTTTTATTTTTTTTGTTAGCTAAT‑GCAAAAGTTAATGAACTAATCCGATGAGTTAATGTTG‑AACAAAAATCATGTTGCGTGGGGGTCGCTTTTTCCCCAGAAGCGTTTAttcc                                                                                                                >  2:269513/1‑246 (MQ=255)                                                                                                                                       aGTAATTTCAGGCTACGCGTGAAAGACTGAAACTGTTCAGTCGCTGGTGCGGTGGCAACGGAGGGTTTTATGCCGCGTTTCGCGGGAATGGGTGCGACGATGACAGTCTCCGTGTTCTGTAGCGTGGAAATCATTTTCATTTTTATGGTTAGAGAAT‑GCAATAGTTACGGAACTGATCCGATGAGTTAATGTTG‑AACAAATCTCATGTTGCGTGGTGGTCGGTTTTACCACAGATGCGTTTAGGGCAGTAt                                                                                                           <  2:539727/251‑1 (MQ=255)                                                                                                                                          aaTTTCAGGCCACGCGTTAAAGACTGAACCTGTCCAGTCGCTGGTGCGGTGGCAACGGCGGGTTTTTTGCCGCGTTTTGCGGGGATGGGTGCGGCCATGACAGTCTCCTTTTTTTGTATCGGGGAAATCATTTTTATTTTTTTTTTTAGCTTAT‑GCAATAGTTTCTGAACTGATCCGATGAGTTTATGTTG‑TAAAAAACTCATGTTGGGTGGTGGGCGCTTTTTCCACAGCTGCGGTTAttcccggttggt                                                                                                        >  1:428774/1‑240 (MQ=255)                                                                                                                                          aaTTTCAGGCCACGCGTTAAAGACTGAACCTGTCCAGTCGCTGGTGCGGTGGCAACGGCGGGTTTTCTGCCGCGTTTCGCGGGAATGGGTGCGACCATGACAGTCTCCTTTTTCTGTATCGTGGAAATCATTTTCATTTTTATTGTTAGCTAAT‑GCAATAGTTACTGAACTGATCCGATGAGTTAATGTTG‑AACAAATCTCATGTTGCGTGGTGGTCGCTTTTACCACAGATGCGTTTATGCCAGTATgg                                                                                                         <  1:214569/250‑1 (MQ=255)                                                                                                                                          aaTTTCAGGCCACGCGTTAAAGACTGAACCTGTCCAGTCGCTGGTGCGGTGGCAACGGCGGGTTTTCTGCCGCGTTTCGCGGGAATGGGTGCGACCATGACAGTCTCCTTTTTCTGTATCGTGGAAATCATTTTCATTTTTATTGTTAGCTAAT‑GCAATAGTTACTGAACTGATCCGATGAGTTAATGTTG‑AACAAATCTCATGTTGCGTGGTGGTCGCTTTTACCACAGATGCGTTTATGCCAGTATgg                                                                                                         <  1:25433/250‑1 (MQ=255)                                                                                                                                            tttCAGGCCACGCGTTAAAGACTGAACCTGTCCAGTCGCTGGTGCGGTGGCAACGGCGGGTTTTCTGCCGCGTTTCGCGGGAATGGGGGCGACCATGACAGTCTCCTTTTTTTGTATCGTGGGAATCATTTTTATTTTTTTTGTTAGCTAAT‑GCAATAGTTACTGAACTTATCCGATGAGTTTATGTTT‑AACAAATATCATTTTGCGTGGGGGTCGCTTTTTTCTCAGATGCGTTTATGCCCGTATGGttt                                                                                                      >  1:150361/1‑251 (MQ=255)                                                                                                                                                  gCCACGCGTTAAAGACTTAACCTGTCCAGTCGCTGGTGCGGTGGCAACGGCGGGTTTTCTGCCGCGTTTCGCGGGAATGGGTGCGACCATGACAGTCTCCTTTTTCTGTATCGTGGAAATCATTTTCAGTTTTATTGTTAGCTAAT‑GCAATAGTTACTGAACTGATCCGATGAGTTAATGTTG‑AACAAAGCTCATGTTGCGTGGTGGTCGCTTTTACCACAGATGCGTTTATGCCAGTATGGTTCGTTGaa                                                                                                <  2:477089/251‑1 (MQ=255)                                                                                                                                                            aaaGACTGAACCTGTCCAGTCGCTGGTGCGGTGACATCTGCGGGTTTTCTGCCGCGTTTCGCGGGAATTGGTTCGACCATGACAGTTTCCTTTGTCTGTATCGTGGAATTCATTTTCATTTTTATTGTTAGCTAAT‑GCAATAGTTCCTGAACTGATCCGATGAGTTAATGTTG‑AACAAATCTCATGTTGCGTGGTGGGCGCTTTTCCCACAGATGCGTTTATGCCAGTATGGTTTGTTGAAGTTTTATTaa                                                                                      <  1:371453/251‑1 (MQ=255)                                                                                                                                                                    aaCCTGTCCAGTCGCTGGTGCGGTGGCAACGGCGGGTTTTCTGCCGCGTTTTGCGGGGATGGGGGCGACCATGACAGTCTCCTTTTTTTGTATCGTGGGAATCATTTTTATTTTTTTTTTTAGCTAAT‑GCAATAGTTACTGAACTGATCCGATGAGTTTATGTTG‑AAAAAAACTCATGTTGCGTGGTGGTCGCTTTTTCCACAAAAGGGTTTTTGCCAGTATGGTTTTTTGAATTTTTTTTAAATTTgggg                                                                              >  1:418131/1‑250 (MQ=255)                                                                                                                                                                      cctttCCAGTCTCTTGTGGGGTGGCAACGGCTGGTTTTCTGCCGTGTTTCGCGGGGATGGGTGAGTCCATTACAGTCTCCTTTTTCTGTTTCGTGGAAATCATTTTCATTTTTATTGGTAACTGAT‑GCAAAAGGTGCTGCACTGATACGATGAGTTAATGCTT‑AACAAAGCTCTTGTTGCGTGGTGGTCGCTTTTACCACAGATGCGTTTATGCCAGGATGTTTTGTTGAATGTTTATTAAATCTGGGtgg                                                                            <  1:467918/247‑3 (MQ=255)                                                                                                                                                                      ccTGTCCAGTCGCTGGTGCGGTGGCAACGGCGGGTTTTCTGCCGCGTTTCGCGGGAATGGGTGCGACCATGACAGTCTCCTTTTTCTGTATCGTGGAAATCATTTTCATTTTTATTGTTAGCTAAT‑GCAATAGTTACTGAACTGATCCGATGAGTTAATGTTG‑AACAAATCTCATGTTGCGTGGTGGTCGCTTTTACCACAGATGCGTTTATGCCAGTATGGTTTGTTGAATTTTTATTAAATCTGGGTTg                                                                            <  1:467916/251‑1 (MQ=255)                                                                                                                                                                            cAGTCGCTGGTGCGGTGGCAACGGCGGGTTTTCTGCCGCGTTTCGCGGGAATGGGTGCGACCATGACAGTCTCCTTTTTCTGTATCGTGGAAATCATTTTCATTTTTATTGTTAGCTAAT‑GCAATAGTTACTGAACTGATCCGATGAGTTAATGTTG‑AACAAATCTCATGTTGCGTGGTGGTCGCTTTTACCACAGATGCGTTTATGCCAGTATGGTTTGTTGAATTTTTATTAAATCTGGGTTGAGCgtg                                                                      <  1:132120/251‑1 (MQ=255)                                                                                                                                                                               tCGCTGGTGCGGTGGCAACGGCGGGTTTTTTGCCGCGCTTTGCGGGGATGGGGGGGAACATGAAAGTCTCCCTTTTTTTTTTCGTGGGAATCATTTTCATTTTTTTTGTTTGCTAAT‑GTAATAGTTACTTAACTGATTCGATGTGTTTATGTTTTAACAA‑TATCTTGTTTTGTGGGGGGCTCTTTTTCCACCGATGCGTTTATTCCAGTGTGGTTTGTTGTATTTTTGTTTAATCTGGGGTTAGCGTGTCg                                                                   >  2:538584/1‑251 (MQ=255)                                                                                                                                                                                             gCAACGGCGGGTTTTCTGCCGCGTTTCGCGGGAATGGGTGCGACCATGACAGTCTCCTTTTTCTGTATCGTGGAAATCATTTTCATTTTTATTGTTAGCTAAT‑GCAATAGTTACTGAACTGATCCGATGAGTTAATGTTG‑AACAAATCTCATGTTGCGTGGTGGTCGCTTTTACCACAGATGCGTTTATGCCAGTATGGTTTGTTGAATTTTTATTAAATCTGGGTTGAGCGTGTCGGGAGCAAGTGTGag                                                     <  2:386600/251‑1 (MQ=255)                                                                                                                                                                                                    cggGTTTTCTGCCGCGTTTCGCGGGAATGGGTGCGACCATGACAGTCTCCTTTTTCTGTATCGTGGAAATCATTTTCATTTTTATTGTTAGCTAAT‑GCAATAGTTACTGAACTGATCCGATGAGTTAATGTTG‑AACAAATCTCATGTTGCGTGGTGGTCGCTTTTACCACAGATGCGTTTATGCCAGTATGGTTTGTTGAATTTTTATTAAATCTGGGTTGAGCGTGTCGGGAGCAAGTGTGAGCAGCaa                                               >  2:127605/1‑250 (MQ=255)                                                                                                                                                                                                       gTTTTCTGCCGCGTTTCGCGGGAATGGGTGCGACCATGACAGTCTCCTTTTTCTGTATCGTGGAAATCATTTTCATTTTTATTGTTAGCTAAT‑GCAATAGTTACTGAACTGATCCGATGAGTTAATGTTG‑AACAAATCTCATGTTGCGTGGTGGTCGCTTTTACCACAGATGCGTTTATGCCAGTATGGTTTGTTGAATTTTTATTAAATCTGGGTTGAGCGTGTCGGGAGCAAGTGTGAGCAGCAAAga                                            >  2:153544/1‑249 (MQ=255)                                                                                                                                                                                                       gTTTTCTGCCGCGTTTCGCGGGAATGGGTGCGACCATGACAGTCTCCTTTTTCTGTATCGTGGAAATCATTTTCATTTTTATTGTTAGCTAAT‑GCAATAGTTACTGAACTGATCCGATGAGTTAATGTTG‑AACAAATCTCATGTTGCGTGGTGGTCGCTTTTACCACAGATGCGTTTATGCCAGTATGGTTTGTTGAATTTTTATTAAATCTGGGTTGAGCGTGTCGGGAGCAAGTGTGAGCAGCAAAGTg                                           <  2:109486/251‑1 (MQ=255)                                                                                                                                                                                                         tttCTGCCGCGTTTCGCGGGAATGGGTGCGACCATGACAGTCTCCTTTTTCTGTATCGTGGAAATCATTTTCATTTTTATTGTTAGCTAAT‑GCAATAGTTACTGAACTGATCCGATGAGTTAATGTTG‑AACAAATCTCATGTTGCGTGGTGGTCGCTTTTACCACAGATGCGTTTATGCCAGTATGGTTTGTTGAATTTTTATTAAATCTGGGTTGAGCGTGTCGGGAGCAAGTGTGAGCAGCAAAGTgg                                          <  1:305053/250‑1 (MQ=255)                                                                                                                                                                                                          ttCTGCCGCGTTTCGCGGGAATGGGTGCGACCATGACAGTCTCCTTTTTCTGTATCGTGGAAATCATTTTCATTTTTATTGTTAGCTAAT‑GCAATAGTTACTGAACTGATCCGATGAGTTAATGTTG‑AACAAATCTCATGTTGCGTGGTGGTCGCTTTTACCACAGATGCGTTTATGCCAGTATGGTTTGTTGAATTTTTATTAAATCTGGGTTGAGCGTGTCGGGAGCAAGTGTGAGCAGCAAAGTGGa                                         <  2:414210/250‑1 (MQ=255)                                                                                                                                                                                                          ttCTGCCGCGTTTCGCGGGAATGGGTGCGACCATGACAGTCTCCTTTTTCTGTATCGTGGAAATCATTTTCATTTTTATTGTTAGCTAAT‑GCAATAGTTACTGAACTGATCCGATGAGTTAATGTTG‑AACAAATCTCATGTTGCGTGGTGGTCGCTTTTACCACAGATGCGTTTATGCCAGTATGGTTTGTTGAATTTTTATTAAATCTGGGTTGAGCGTGTCGGGAGCAAGTGTGAGCAGCAAAGTGGa                                         <  2:64056/250‑1 (MQ=255)                                                                                                                                                                                                              gCCGCGTTTCGCGGGAATGGGTGCGACCATGACAGTCTCCTTTTTCTGTATCGTGGAAATCATTTTCATTTTTATTGTTAGCTAAT‑GCAATAGTTACTGAACTGATCCGATGAGTTAATGTTG‑AACAAATCTCATGTTGCGTGGTGGTCGCTTTTACCACAGATGCGTTTATGCCAGTATGGTTTGTTGAATTTTTATTAAATCTGGGTTGAGCGTGTCGGGAGCAAGTGTGAGCAGCAAAGTGGaac                                       <  1:542868/248‑1 (MQ=255)                                                                                                                                                                                                                         cgGGAATGGGTGCGACCATGACAGTCTCCTTTTTCTGTATCGTGGAAATCATTTTCATTTTTATTGTTAGCTAAT‑GCAATAGTTACTGAACTGATCCGATGAGTTAATGTTG‑AACAAATCTCATGTTGCGTGGTGGTCGCTTTTACCACAGATGCGTTTATGCCAGTATGGTTTGTTGAATTTTTATTAAATCTGGGTTGAGCGTGTCGGGAGCAAGTGTGAGCAGCACAGTGGAACAACTGCGTGCGCAg                         <  2:70211/251‑1 (MQ=255)                                                                                                                                                                                                                         cgGGAATGGGTGCGACCATGACAGTCTCCTTTTTCTGTATCGTGGAAATCATTTTCATTTTTATTGTTAGCTAAT‑GCAATAGTTACTGAACTGATCCGATGAGTTAATGTTG‑AACAAATCTCATGTTGCGTGGTGGTCGCTTTTACCACAGATGCGTTTATGCCAGTATGGTTTGTTGAATTTTTATTAAATCTGGGTTGAGCGTGTCGGGAGCAAGTGTGAGCAGCAAAGTGGAACAACTGCGTGCGCAg                         <  2:196103/251‑1 (MQ=255)                                                                                                                                                                                                                          gggAATGGGTGCGACCATGACAGTCTCCTTTTTCTGTATCGTGGAAATCATTTTCATTTTTATTGTTAGCTAAT‑GCAATAGTTACTGAACTGATCCGATGAGTTAATGTTG‑AACAAATCTCATGTTGCGTGGTGGTCGCTTTTACCACAGATGCGTTTATGCCAGTATGGTTTGTTGAATTTTTATTAAATCTGGGTTGAGCGTGTCGGGAGCAAGTGTGAGCAGCAAAGTGGAACAACTGCGTGCGCAGt                        >  1:552166/1‑251 (MQ=255)                                                                                                                                                                                                                           ggAATGGGTGCGACCATGACAGTCTCCTTTTTCTGTATCGTGGAAATCATTTTCATTTTTATTGTTAGCTAAT‑GCAATAGTTACTGAACTGATCCGATGAGTTAATGTTG‑AACAAATCTCATGTTGCGTGGTGGTCGCTTTTACCACAGATGCGTTTATGCCAGTATGGTTTGTTGAATTTTTATTAAATCTGGGTTGAGCGTGTCGGGAGCAAGTGTGAGCAGCAAAGTGGAACAACTGCGTGCGCAGtt                       <  1:206959/251‑1 (MQ=255)                                                                                                                                                                                                                                 ggTGCGACCATGACAGTCTCCTTTTTCTGTATCGTGGAAATCATTTTCATTTTTATTGTTAGCTAAT‑GCAATAGTTACTGAACTGATCCGATGAGTTAATGTTG‑AACAAATCTCATGTTGCGTGGTGGTCGCTTTTACCACAGATGCGTTTATGCCAGTATGGTTTGTTGAATTTTTATTAAATCTGGGTTGAGCGTGTCGGGAGCAAGTGTGAGCAGCAAAGTGGAACAACTGCGTGCGCAGTTAAATg                  >  1:244085/1‑250 (MQ=255)                                                                                                                                                                                                                                   tGCGACCATGACAGTCTCCTTTTTCTGTATCGTGGAAATCATTTTCATTTTTATTGTTAGCTAAT‑GCAATAGTTACTGAACTGATCCGATGAGTTAATGTTG‑AACAAATCTCATGTTGCGTGGTGGTCGCTTTTACCACAGATGCGTTTATGCCAGTATGGTTTGTTGAATTTTTATTAAATCTGGGTTGAGCGTGTCGGGAGCAAGTGTGAGCAGCAAAGTGGAACAACTGCGTGCGCAGTTAAATGAAc               >  1:442560/1‑251 (MQ=255)                                                                                                                                                                                                                                      gACCATGACAGTCTCCTTTTTCTGTATCGTGGAAATCATTTTCATTTTTATTGTTAGCTAAT‑GCAATAGTTACTGAACTGATCCGATGAGTTAATGTTG‑AACAAATCTCATGTTGCGTGGTGGTCGCTTTTACCACAGATGCGTTTATGCCAGTATGGTTTGTTGAATTTTTATTAAATCTGGGTTGAGCGTGTCGGGAGCAAGTGTGAGCAGCAAAGTGGAACAACTGCGTGCGCAGTTAAATGAACGTa            <  2:249515/251‑1 (MQ=255)                                                                                                                                                                                                                                        ccATGACAGTCTCCTTTTTCTGTATCGTGGAAATCATTTTCATTTTTATTGTTAGCTAAT‑GCAATAGTTACTGAACTGATCCGATGAGTTAATGTTG‑AACAAATCTCATGTTGCGTGGTGGTCGCTTTTACCACAGATGCGTTTATGCCAGTATGGTTTGTTGAATTTTTATTAAATCTGGGTTGAGCGTGTCGGGAGCAAGTGTGAGCAGCAAAGTGGAACAACTGCGTGCGCAGTTAAATGAACGTAtt          <  2:321527/251‑1 (MQ=255)                                                                                                                                                                                                                                           tGACAGTCTCCTTTTTCTGTATCGTGGAAATCATTTTCATTTTTATTGTTAGCTAAT‑GCAATAGTTACTGAACTGATCCGATGAGTTAATGTTG‑AACAAATCTCATGTTGCGTGGTGGTCGCTTTTACCACAGATGCGTTTATGCCAGTATGGTTTGTTGAATTTTTATTAAATCTGGGTTGAGCGTGTCGGGAGCAAGTGTGAGCAGCAAAGTGGAACAACTGCGTGCGCAGTTAAATGAACGTATTCTg       >  1:22787/1‑251 (MQ=255)                                                                                                                                                                                                                                             aCAGTCTCCTTTTTCTGTATCGTGGAAATCATTTTCATTTTTATTGTTAGCTAAT‑GCAATAGTTACTGAACTGATCCGATGAGTTAATGTTG‑AACAAATCTCATGTTGCGTGGTGGTCGCTTTTACCACAGATGCGTTTATGCCAGTATGGTTTGTTGAATTTTTATTAAATCTGGGTTGAGCGTGTCGGGAGCAAGTGTGAGCAGCAAAGTGGAACAACTGCGTGCGCAGTTAAATGAACGTATTCTGGt     <  1:227998/251‑1 (MQ=255)                                                                                                                                                                                                                                                gTCTCCTTTTTCTGTATCGTGGAAATCATTTTCATTTTTATTGTTAGCTAAT‑GCAATAGTTACTGAACTGATCCGATGAGTTAATGTTG‑AACAAATCTCATGTTGCGTGGTGGTCGCTTTTACCACAGATGCGTTTATGCCAGTATGGTTTGTTGAATTTTTATTAAATCTGGGTTGAGCGTGTCGGGAGCAAGTGTGAGCAGCACAGTGGAACAACTGCGTGCGCAGTTAAATGAACGTATTCTGGTGCt  >  2:517220/1‑251 (MQ=255)                                                                                                                                                                                                                                                     |                                                                                                                                                                                                                                                        GTTCGCCAACCTGACGCACGAAACCCTGCTGTTGCATCGTGGTTAGCAGGCGGTGGGTCGTGGAATTGGGTAACCCGGCTTGTTGCGCCAGTTCCGTGAGTGCCACACTGCCATTGGATTCGGCAATCCACTCCAGTAATTTCAGGCCACGCGTTAAAGACTGAACCTGTCCAGTCGCTGGTGCGGTGGCAACGGCGGGTTTTCTGCCGCGTTTCGCGGGAATGGGTGCGACCATGACAGTCTCTTTTTTCTGTATCGTGGAAATCATTTTCATTTTTATTGTTAGCTAAT‑GCAATAGTTACTGAACTGATCCGATGAGTTAATGTTG‑AACAAATCTCATGTTGCGTGGTGGTCGCTTTTACCACAGATGCGTTTATGCCAGTATGGTTTGTTGAATTTTTATTAAATCTGGGTTGAGCGTGTCGGGAGCAAGTGTGAGCAGCAAAGTGGAACAACTGCGTGCGCAGTTAAATGAACGTATTCTGGTGCT  >  NC_000913/4223394‑4223883

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