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%	159.1 / NA	50	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 (8/42);  total (8/42)

GCCAACCTGACGCACGAAACCCTGCTGTTGCATCGTGGTTAGCAGGCGGTGGGTCGTGGAATTGGGTAACCCGGCTTGTTGCGCCAGTTCCGTGAGTGCCACACTGCCATTGGATTCGGCAATCCACTCCAGTAATTTCAGGCCACGCGTTAAAGACTGAACCTGTCCAGTCGCTGGTGCGGTGGCAACGGCGGGTTTTCTGCCGCGTTTCGCGGGAATGGGTGCGACCATGACAGTCTCTTTTTTCTGTATCGTGGAAATCATTTTCATTTTTATTGTTAGCTAATGCAATAGTTACTGAACTGATCCGATGAGTTAATGTTGAACAAATCTCATGTTGCGTGGTGGTCGCTTTTACCACAGATGCGTTTATGCCAGTATGGTTTGTTGAATTTTTATTAAATCTGGGTTGAGCGTGTCGGGAGCAAGTGTGAGCAGCAAAGTGGAACAACTGCGTGCGCAGTTAAATGAACGTAT  >  NC_000913/4223398‑4223874                                                                                                                                                                                                                                                 |                                                                                                                                                                                                                                             gCCAACCTGACGCACGAAACCCTGCTGTTGCATCGTGGTTAGCAGGCGGTGGGTCGTGGAATTGGGTAACCCGGCTTGTTGCGCCAGTTCCGTGAGTGCCACACTGCCATTGGATTCGGCAATCCACTCCAGTAATTTCAGGCCACGCGTTAAAGACTGAACCTGTCCAGTCGCTGGTGCGGTGGCAACGGCGGGTTTTCTGCCGCGTTTCGCGGGAATGGGTGCGACCATGACAGTCTCCTTTTTCTGTa                                                                                                                                                                                                                                    <  1:406155/251‑1 (MQ=255)   cAACCTGACGCACGAAACCCTTCTGTTGCATCGTGGTTAGCAGGCGGTGGGTCGTGGAATTGGGTACCCCGGCTTGTTGCGCCAGTTCCGTGAGTGCCACACTGCCATTGGATTCGGCAATCCACTCCAGTAATTTCAGGCCACGCGTTAAAGACTGAACCTGTCCAGTCGCTGGTGCGGTGGCAACGGCTGGTTTTCTGCCTCGTTTCGCGGGAATGGGTGCGACCATGACAGTCTCCTTTTTCTGTATc                                                                                                                                                                                                                                  <  1:280460/251‑1 (MQ=255)    aaCCTGACGCACGAAACCCTGCTGTTGCATCGTGGTTAGCAGGCGGTGGGTCGTGGAATTGGGTAACCCGGCTTGTTGCGCCAGTTCCGTGAGTGCCACACTGCCATTGGATTCGGCAATCCACTCCAGTAATTTCAGGCCACGCGTTAAAGACTGAACCTGTCCAGTCGCTGGTGCGGTGGCAACGGCGGGTTTTCTGCCGCGTTTCGCGGGAATGGGTGCGACCATGACAGTCTCCTTTTTCTGTATCg                                                                                                                                                                                                                                 <  1:411358/251‑1 (MQ=255)                    cccTGCTGTTGTATCGTGGTTAGCAGGGGGTGGGTCGTGGAATTTGTGAACCCGGCTTGTTGCCCCAGTTCCGTGAGGTCCACGCTGCCATTGGATTCGGCAATCCACTCCAGTAATTTCAGGCCACGCGTTAAAGACTGACCCTGTCCAGTCGCTGGTGCGGTGGCAACGGCGGGTTTTCTGCCGCGTTTCGCGGGAATGGGTGCGACCATGACAGTCTCCTTTTTCTGTATCGTGGAAAtcattttcat                                                                                                                                                                                                                 <  2:13759/251‑1 (MQ=255)                     ccggctgTTGCATCGTGGTTAGCAGGCGGTGGGGCGTGGAATTGTGTAACCCGGCTTGTTGCGCCAGTTCCGTGAGTGCCACACTGCCATTGGATTCGGCAATCCACTCCAGTAATTTCAGGCCACGCGTTAAAGACTGAACCTGTCCAGTCGCTGGTGCGGTGGCAACGGCGGGTTTTCTGCCGCGTTTCGCGGGAATGGGTGCGACCATGACAGTCTCCTTTTTCTGTATCGTGGAAAtcattttcatt                                                                                                                                                                                                                <  2:24575/248‑1 (MQ=255)                        gctgTTGCATCGTGGTTAGCAGGCGGTGGGTCGTGGAATTGGGTAACCCGGCTTGTTGCGCCAGTTCCGTGAGTGCCACACTGCCATTGGATTCGGCAATCCACTCCAGTAATTTCAGGCCACGCGTTAAAGACTGAACCTGTCCAGTCGCTGGTGCGGTGGCAACGGCGGGTTTTCTGCCGCGTTTCGCGGGAATGGGTGCGACCATGACAGTCTCCTTTTTCTGTATCGTGGAAATcattttcatttt                                                                                                                                                                                                              <  1:368937/250‑1 (MQ=255)                              gCATCGTGGTTAGCAGGCGGTGGGTCGTGGAATTGGGTAACCCGGCTTGTTGCGCCAGTTCCGTGAGTGCCACACTGCCATTGGATTCGGCAATCCACTCCAGTAATTTCAGGCCACGCGTTAAAGACTGAACCTGTCCAGTCGCTGGTGCGGTGGCAACGGCGGGTTTTCTGCCGCGTTTCGCGGGAATGGGTGCGACCATGACAGTCTCCTTTTTCTGTATCGTGGAAATCATTTTCATTTTTATTGt                                                                                                                                                                                                        <  1:360563/250‑1 (MQ=255)                               cATCGTGGTTAGCAGGCGGTGGGTCGTGGAATTGGGTAACCCGGCTTGTTGCGCCAGTTCCGTGAGTGCCACACTGCCATTGGATTCGGCAATCCACTCCAGTAATTTCAGGCCACGCGTTAAAGACTGAACCTGTCCAGTCGCTGGTGCGGTGGCAACGGCGGGTTTTCTGCCGCGTTTCGCGGGAATGGGTGCGACCATGACAGTCTCCTTTTTCTGTATCGTGGAAATCATTTTCATTTTTATTGTTa                                                                                                                                                                                                      <  1:490254/251‑1 (MQ=255)                                  cGTGGTTAGCAGGCGGTGGGTCGTGGAATTGGGTAACCCGGCTTGTTGCGCCAGTTTCGTGAGTGCCACACTGCCATTGTATTCGGCAATCCACTCCAGTAATTTCAGGCCACGCGTTAAAGACTGACCCTGTCCAGTCGCTGGTGCGGTGGCAACGGCGGGTTTTCTGCCGCGTTTCGCGGGAATGGGTGCGTCCATGACAGTCTCCTTTTTCTGTATCGTGGAAATCATTTTCATTTTTATTGTTAgtt                                                                                                                                                                                                   <  1:451043/251‑3 (MQ=255)                                         aGCAGGCGGTGGGTCGTGGAATTGGGTAACCCGGCTTGTTGCGCCAGTTCCGTGAGTGCCACACTGCCATTGGATTCGGCAATCCACTCCAGTAATTTCAGGCCACGCGTTAAAGACTGAACCTGTCCAGTCGCTGGTGCGGTGGCAACGGCGGGGTTTTTGCCGCGTTTCTCGGGGGTGGGTGCGGCCATGACAGGCTCCTTTTTCTTTATCGTGGAAATCATTTTCATTTTTTTTGTTAgattattca                                                                                                                                                                                             >  1:480216/1‑242 (MQ=255)                                            aggcgggtggctgTGGAATTGGGTAACCCGGCTTGTTGCGCCAGTTCCGGGTGTGCCACACTGCCATTGGATTCGGCAATCCACTCCAGTAATGTCAGGCCACGCGTTAAAGACTGAACCTGTCCAGTCGCTGGTGCGGGGGCAACGGCGGGTTTTCTGCCGCGTTTCGCGGGAATGGGTGCGACCATGACAGTCTCCTTTTTCTGTATCGTGGAAATCATTTTCATTTTTATTTTTAGCTAATGCAATAg                                                                                                                                                                                         <  2:287139/239‑1 (MQ=255)                                               cGGTGGGTCGTGGAATTGGGTATCCCGGTTTGTTGCGCCAGTTCCGTGAGTGCCACACTTCCATGGGATTCGGCAATCCACTCCAGTAATTTCAGGCCACGCGGTAAAGACTGAACCTGTCCAGTCGCTGGTGCGGTGGCAACGGCGGGTTTTCTGCCGCGTTGCGCGGGAATGGGTGCGACCATGACAGTCTCCTTTTTCTGTATCGTGGAAATCATTTTCATTTTTATTGTTAGCTAATGCAATAGtga                                                                                                                                                                                      <  2:362873/251‑3 (MQ=255)                                                       cgtgtgatagagtagcccGGCTTGTTGCGCCAGTACCTCGAGTGCCTCCCTTGCAGTGGATACGGCAATCCACTCCAGTAATTTCAGGCCACGCCATAAATACTGAACCCGTCCAGTCGGTGGTTCGGTGACAATGGAGGACTTTCTGCCGCGTTTCGCGGTAATGGGTGCGTATTTGACAGTCGCCTATTTCTGTATCGTGGAAACTAGTTCCATTTTTCTTGGTAGATATTGCAATACTTgctgaccc                                                                                                                                                                               <  2:382119/235‑4 (MQ=255)                                                        gTGGAATTGGGTAACCCGGCTTGTTGCGCCATCTCCGTGAGTGCCCCACTGCCATTGGATTCGGCAATCCACTCCAGTAATTTCAGGCCACGCGTTAAAGACTGAACCTGTCCAGTCGCTGGTGCGGTGGCAACGGCGGGTTTTCTGCCGCGTTTCGCGGGAATGTGTGCGACCATGACAGTCTCCTTTTTCTGTATCGTGGAAATCATTTTCATTTTTATTGTTAGCTAATGCAATAGTTcctgaactga                                                                                                                                                                             <  2:69692/251‑1 (MQ=255)                                                                              gttgCGCCAGTTCCGTTAGTGCCACACTGCCATTGGATTCGGCAATCCACTCCAGTAATTTCAGGCCACGCGTTAAAGACTGAACCTGTCCAGTCGCTGGTGCTGTGGCAACGGCGGGTTTTCTGCCGCGTTTCGCGGGAATGGGTGCGACCATGACAGTCTCCTTTTTCTGTATCGTGGAAATCATTTTCATTTTTATTGTTAGCTAATGCAATAGTTACTGAACTTATCCGATGAGTTAATGTTGAACa                                                                                                                                                       <  1:147081/251‑1 (MQ=255)                                                                              gttgCGCCAGTTCCGTGAGTGCCACACTGCCATTGGATTCGGCAATCCACTCCAGTAATTTCAGGCCACGCGTTAAAGACTGAACCTGTCCAGTCGCTGGTGCGGTGGCAACGGCGGGTTTTTTGGCGCGTTTCGCGGGAAAGGGGGGGGCCCTGACCGTCTCCCTTTTTTGTTTTCTGGGAATCATTTTCTTTTTTTTTGTTTGTTAATTCAATAGTTACTTAACTTATTCGTTTGGTTAATGTTgtac                                                                                                                                                        >  1:207895/1‑247 (MQ=255)                                                                              gttgCGCCAGTTCCGTGAGTGCCACACTGCCATTGGATTCGGCAATCCACTCCAGTAATTTCAGGCCACGCGTTAAAGACTGAACCTGTCCAGTCGCTGGTGCGGTGGCAACGGCGGGTTTTCTGCCGCGTTTCGCGGGAATGGGTGCGACCATGACAGTCTCCTTTTTCTGTATCGTGGAAATCATTTTCATTTTTATTGTTAGCTAATGCAATAGTTACTGAACTGATCCGATGAGTTAATGTTGAACa                                                                                                                                                       <  1:258410/251‑1 (MQ=255)                                                                              gttgCGCCAGTTCCGTGAGTGCCACACTGCCATTGGATTCGGCAATCCACTCCAGTAATTTCAGGCCACGCGTTAAAGACTGAACCTGTCCAGTCGCTGGTGCGGTGGCAACGGCGGGGTTTTTGCCGCGGTTCGCGGGAATGGGGGCGGCCATGACAGTCTCCTTTTTTTTTATCTTGGGAATCATTTTTATTTTTTTTTTTAGCTAATGCAATAGTTACTGAACTGATCCGTTGGGTTAATGTTGaaa                                                                                                                                                        >  1:474389/1‑249 (MQ=255)                                                                                          ccGTGAGTGCCACACTGCCATTGGATTCGGCAATCCACTCCAGTAATTTCAGGCCACGCGTTAAAGACTGAACCTGTCCAGTCGCTGGTGCGGTGGCAACGGCGGGTTTTCTGCCGCGTTTCGCGGGAATGGGTGCGACCATGACAGTCTCCTTTTTCTGTATCGTGGAAATCATTTTCATTTTTATTGTTAGCTAATGCAATAGTTACTGAACTGATCCGATGAGTTAATGTTGAACAAATCTCATGTTg                                                                                                                                           <  1:484329/251‑1 (MQ=255)                                                                                              gAGTGCCACACTGCTATTGGATTCGGCAAGCCACTCCAGTTATTTCAGGCCTCGCGTTAAAGACTGAACCTGTCCAGTCGCTGGTGCGGTGGCAACGGCGGGTTTTCTGCCGCGTTTCGCGGTAATGGGTGCGACCATGACAGTCTCCTTTTTCTGTATCGTGGAAATCATTTGCATTTTTATTGTTAGCTAATGCAATAGTTACTGAACTGATCCGATGAGTTAATGTTGAACAAATCGCATGTTGCgtg                                                                                                                                       <  2:99211/251‑1 (MQ=255)                                                                                                 ggCCACACTGCCATTGTATTCGGCAATCCACTCCATCAATTTCAGGCCACGCGTTAAAGACTGAACCTGTCCAGTCGCTTGTGCGGTGGCAACGGCGGGTTTTCTGCCGCGTTTCGCGGGAATGGGTGCGACCATGACAGTCTCCTTTTTCTGTATCGTGTAAATCATTTTCATTTTTATTGTTAGCTAATGCAATAGTTACTGAACTGATCCGATGAGTTAATGTTGAACATATCTCATGTTGCgtggtg                                                                                                                                    <  1:122504/250‑1 (MQ=255)                                                                                                          gcctttggagtctggaaTCCACTCCAGTAATTTCAGGCCACGCGGTGAAGGCTGAACCTGTCCAGACGCTGGTGGGGTGGCAACGGCTGGTTTTCTGCCGCGTTTCGCGGGAATGGGTGCGACCATGACAGTCTCCTGTTTCTGTATCGTGGAAGTCTTTTTCATTTTTATTGTTAGCTAAGGCAATAGTGACTGAACTGATCCGATGAGTTCGTGATGAACAAATTTCATGTTGCGTGGTTGTTGCtttt                                                                                                                           <  2:471868/236‑1 (MQ=255)                                                                                                          gCCATTGGATTCGGCAATCCACTCCAGTAATTTCAGGCCACGCGTTAAAGACTGAACCTGTCCAGTCGCTGGTGCGGTGGCAACGGCGGGTTTTCTGCCGCGTTTCGCGGGAATGGGTGCGACCATGACAGTCTCCTTTTTCTGTATCGTGGAAATCATTTTCATTTTTATTGTTAGCTAATGCAATAGTTACTGAACTGATCCGATGAGTTAATGTTGAACAAATCTCATGTTGCGTGGTGGTCGCttt                                                                                                                            <  2:218252/250‑1 (MQ=255)                                                                                                            cATTGGATTCGGCAATCCACTCCAGTAATTTCAGGCCTCGCGTTAA‑GACTGAACCTGTCCAGTCGCTGGTGCGGTGGCAACGGCGGGTTTTCTGCCGCGTTTCGCGGGAATGGGTGCGACCATGACAGTCTCCTTTTTCTGTATCGTGGAAATCATTTTCATTTTTATTGTTAGCTAATGCAATAGTTACTGAACTGATCCGATGAGTTAATGTTGAACAAATCTCATGTTGCGTGGTGGTCGCTTTTAcc                                                                                                                        <  1:319272/251‑1 (MQ=255)                                                                                                                ggATTCGGCAATCCACTCCAGTAATGTCAGGCCACGCGTTAAAGACTGAACCTGTCCAGTCGCTGGTGCGGTGGCAACGGCGGGTTGTCTGCCGCGTTTCGCGGGAATGGGTGCGACCATGACAGTCTCCTTTTTCTGTATCGTGGAAATCATTTTCATTTTTATTGTTAGCTAATGCAATAGTTACTGAGCTGATCCGATGAGTTAATGTTGAACAAATCTCATGTTGCGTGGTGGTCGCTTTTACcac                                                                                                                      <  1:496056/250‑1 (MQ=255)                                                                                                                 gATTCGGCAATCCACTCCAGTAATTTCAGGCCACGCGGTAAAGACTGAACCTGTCCAGTCGCTGGTGCGGTGGCAACGGCGGGTTTTCTGCCGCGTTTCGCGGGAATGGGTGCGACCATGACAGTCTCCTTTTTCTGTATCGTGGAAATCATTTTCATTTTTATTGTTAGCTAATGCAATAGTTACTGAACTGATCCGATGAGTTAATGTTGAACAAATCTCATGTTGCGTGGTGGTCGCTTTTACCACAg                                                                                                                    <  1:176356/251‑1 (MQ=255)                                                                                                                            ccACTCCAGTAATGTCAGGCCACGCGTTAAAGACTGAACCTGTCCAGTCGCTGGTGCGGTGGCAACGGCGGGTTTTCTGCCGCGTTTCGCGGGAATGGGTGCGACCATGACAGTCTCCTTTTTCTGTATCGTGGAAATCATTTTCATTTTTATTGTTAGCTAATGCAATAGTTACTGAACTGATCCGATGAGTTAATGTTGAACAAATCTCATGTTGCGTGGTGGTCGCTTTTACCACAGATGCGTTTAt                                                                                                          <  2:430398/250‑1 (MQ=255)                                                                                                                             cACTCCAGTAATTTCAGGCCACGCGTTAAAGACTGAACCTGTCCAGTCGCTGGTGCGGTGGCAACGGCTGGTTTTCTGCCGCGTTTCGCGGGAATGGGTGCGACCATGACAGTCTCCTTTTTCTGTATCGTGGAAATCATTTTCATTTTTATTGTTAGCTAATGCAATAGTTACTGAACTGATCCGATGAGTTAATg                                                                                                                                                              <  2:102390/197‑1 (MQ=255)                                                                                                                              aCTCCAGTAATTTCAGGCCACGCGTTAAAGACTGAACCTGTCCAGTCGCTGGTGCGGTGGCAACGGCGGGGTTTTTGCCGCGTTTTGCGGGGATGGGTGCGACCCTGACAGTATCCTTTTTCTATATCGTGGAAAAAATTTTCATTTTTTTTGTTAGCTTAAGCAATAATTAATGTACTGATCCGATGAGGTAATTTTGAAAAAATCTCATTTTGGGTGGGGGGGGCTTTTACCACAAATGCGttttttgc                                                                                                       >  1:34645/1‑246 (MQ=255)                                                                                                                                   aGTAATTTCAGGCCACGCGTTAAAGACTGAACCTGTCCAGTCGCTGGTGCGGTGGCAACGGCGGGTTTTCTGCCGCGTTTCGCGGGAATGGGTGCGACCATGACAGTCTCCTTTTTCTGTATCGTGGAAATCATTTTCATTTTTATTGTTAGCTAATGCAATAGTTACTGAACTGATCCGATGAGTTAATGTTGAACAAATCTCATGTTGCGTGGTGGTCGCTTTTACCACAGATGCGTTTATGCCAGTAt                                                                                                  <  1:369331/251‑1 (MQ=255)                                                                                                                                   aGTAATTTCAGGCCACGCGTTAAAGACTGAACCTGTCCAGTCGCTGGTGCGGTGGCAACGGCGGGGTTTTTGTCGCGGTTTGCGGGGATGGGGGCGACCATCACAGTATCCTTTTTCTTTATCGTGGGAATAATTTTTATTTTTATTGTTAGCTAATGCAATAGTTACTGTACTGATCCGATGAGTTAATGTTTAAAAAATCTCAAGTTGCGTGGGGGTCGCTTTTACCCCAGATGCGGTTATTCCCGTAt                                                                                                  >  1:216042/1‑251 (MQ=255)                                                                                                                                    gTAATTTCAGGCCACGCGTTAAAGACTGAACCTGTCCAGTCGCTGGTGCGGTGGCAACGGCGGGTTTTCTGCCGCGTTTCGCGGGAATGGGTGCGACCATGACAGTCTCCTTTTTCTGTATCGTGGAAATCATTTTCATTTTTATTGTTAGCTAATGCAATAGTTACTGAACTGATCCGATGAGTTAATGTTGAACAAATCTCATGTTGCGTGGTGGTCGCTTTTCCCACAGATGCGTTTATGCCAGTATg                                                                                                 <  2:299780/251‑1 (MQ=255)                                                                                                                                           cAGGCCATGCGTGAAAGACTGCACCTGTCCAGTCGCTGGTGCGGTGGCAACGGCGGGTTTTCTGCCGCGTTTCGCGGGAATGGGTGCGACCATGACAGTCTCCTTTTGCTGTATCGTGGAAATCATTTTCATTTTTATTGTTAGCTAATGCCAGAGTTACTGAACTGATCCGATGAGTTAATGTTGAACAAATCTCAGGTTGCGTGGTGGTCGCTTTTACCACAGATGCGTTGATGCCATTATGGTttggt                                                                                          <  2:75473/251‑3 (MQ=255)                                                                                                                                           cAGGCCACGCGTTAAAGACTTAACCTGTCCAGTCGCTGGTGCGGTGGCAACGGCGGGTTTTCTGCCGCGTTTCGCGGGAATGGGTGCGACCATGACAGTCTCCTTTTTCTGTATCGTGGAAATCATTTTCATTTTTATTGTTAGCTAATACAATAGTTACTGAACTGATCCGATGAGTTAATGTTGAACAAATCTCATGTTGCGTGGTGGTCGCTTTTCCCACAGATGCGTTTATGCCAGTATGGTttgtt                                                                                          <  1:275454/251‑1 (MQ=255)                                                                                                                                             ggCCACGGGTTAAAGACTGAACCTGTCCAGTCGCTGGTGCGGTGGCAACAGCGGGTTTTCTGCCGCGTTTCGCGGGAATGGGTGCGACCATGACAGTCTCCTTTTTCTGTATCGTGGAAATCATTTTCATTTTTATTGTTAGCTAATGCAATAGTTACTGAACTGATCCGATGAGTTAATGTTGAACAAATCTCATGTTGCGTGGTGGTCGCTTTTACCACAGATGCGTTTATGCCAGTATGGTTTGTTGa                                                                                        <  2:411228/251‑1 (MQ=255)                                                                                                                                             ggCCACGCGTTAAAGACTGAACCTGTCCAGTCGCTGGTGCGGTGGCAACGGCGGGTTTTCTGCCGCGTTTCGCGGGAATGGGTGCGACCATGACAGTCTCCTTTTTCTGTATCGTGGAAATCATTTTCATTTTTATTGTTAGCTAATGCAATAGTTACTGAACTGATCCGATGAGTTAATGTTGAACAAATCTCATGTTGCGTGGTGGTCGCTTTTACCACAGATGCGTTTATGCCAGTATGGTTTGTTGa                                                                                        <  1:1431/251‑1 (MQ=255)                                                                                                                                             ggCCACGCGTTAAAGACTGAACCTGTCCAGTCGCTGGTGCGGTGGCAACGGCGGGTTTTCTGCCGCGTTTCGCGGGAATGGGTGCGACCATGACAGTCTCCTTTTTCTGTATCGTGGAAATCATTTTCATTTTTATTGTTAGCTAATGCAATAGTTACTGAACTGATCCGATGAGTTAATGTTGAACAAATCTCATGTTGCGTGGTGGTCGCTTTTACCACAGATGCGTTTATGCCAGTATGGTTTGTTGa                                                                                        <  1:43202/251‑1 (MQ=255)                                                                                                                                               ccACGCGTTAAAGACTGAACCTGTCCAGTCGCTGGTGCGGTGGCAACGGCGGGTTTTCTGCCGCGTTTCGCGGGAATGGGTGCGACCATGACAGTCTCCTTTTTCTGTATCGTGGAAATCATTTTCATTTTTATTGTTAGCTAATGCAATAGTTACTGAACTGATCCGATGAGTTAATGTTGAACAAATCTCATGTTGCGTGGTGGTCGCTTTTACCACAGATGCGTTTATGCCAGTATGGTTTGTTGAAt                                                                                      <  1:70674/251‑1 (MQ=255)                                                                                                                                                   gcgTTAAAGACTGAACCTGTCCAGTCGCTGGTGCGGTGGCAACGGCGGGTTTTCTGCCGCGTTTCGCGGGAATGGGTGCGACCATGACAGTCTCCTTTTTCTGTATCGTGGAAATCATTTTCAGTTTTATTGTTAGCTAATGCAATAGTTACTGAACTGATCCGATGAGGTAATGTTGAACAAATCTCATGTTGCGTGGTGGTCGCTTTTACCACAGATGCGTTTATGCCAGTATGGTTTGTTGAAttttt                                                                                  <  1:230919/251‑1 (MQ=255)                                                                                                                                                     ggtAAAGACTGAACCTGTCCAGTCGCTGGTGCGGTGGCAACGGCGGGTTTTCTGCCGCGTTTCGCGGGAATGGGTGCGACCATGACAGTCTCCTTTTTCTGTATCGTGGAAATCATTTTCATTTTTATTGTTAGCTAATGCAATAGTTACTGAACTGATCCGATGAGTTAATGTTGAACAAATCTCATGTTGCGTGGTGGTCGCTTTTACCACAGATGCGTTTATGCCAGTATGGTTTGTTGAATTTttat                                                                                <  2:228160/249‑1 (MQ=255)                                                                                                                                                       tAAAGTCTGAAACTGTCCAGTCGATGGTGCTGTGGCTACGTAGGGTTTTCTGCCGCTTTTCGCGGGAATGGTGGCGAACATGACAGGCTCCTTTTTCTTTATCGTGGAAATCATTTTCATTTTTATTGTTAGCTAATGCAATAGTTACTGAACTGATCCGATGAGTTAATGTTGAACAAATCTCAGGTTGTGTGGTGGTCGCTTTTACCACAGATGCGTTTATGCTTGTATGGTTGGTTGAAtttttctt                                                                               <  2:77411/250‑4 (MQ=255)                                                                                                                                                           gACTGAACCTGTCCAGTCGCTGGTGCGGTGGCAACGGCGGGTTTTCTGCCGCGTTTCGCGGGAATGGGTGCGACCATGACAGTCTCCTTTTTCTGTATCGTGGAAATCATTTTCATTTTTATTGTTAGCTAATGCAATAGTTACTGAACTGATCCGATGAGTTAATGTTGAACAAATCTCATGTTGCGTGGTGGTCGCTTTTACCACAGATGCGGTTATGCCAGTATGGTTTGTTGAATTTTTATTAAATc                                                                          <  2:237295/251‑1 (MQ=255)                                                                                                                                                                                         gCAACGGCGGGTTTTCTGCCGCGTTTCGCGGGAATGGGTGCGACCATGACAGTCTCCTTTTTCTGTATCGTGGAAATCATTTTCATTTTTATTGTTAGCTAATGCAATAGTTACTGAACTGATcc                                                                                                                                                                          <  1:405257/125‑1 (MQ=255)                                                                                                                                                                                                   gTTTTCTGCCGCGTTTCGCGGTAATGGGTGCGACCATGACAGTCTCCTTTTTCTGTATCGTGGAAATCATTTTCATTTTTATTGTTAGCTAATGCAATAGTTACTGAACTGATCCGATGAGTTAATGTTGAACAAATCTCATGTTGCGTGGTGGTCGCTTTTACCACAGATGCGTTTATGCCAGTATGGTTTGTTGAATTTTTATTAAATCTGGGTTGAGCGTGTCGGGAGCAAGTGTGAGCAGCAAAGTg                                  <  2:229270/251‑1 (MQ=255)                                                                                                                                                                                                       acggCCGCGTGTCGCGGTAATAGGTGCGACCATGACAGTCTCCTTTTTCTGTATCGTGGAAATTATTTTCATTTTTATTGTCAGCTAATGCAATAGTTACTGAACTGATCCGATGAGTTAATGGTTAACAAATCTCATGTTGCGTGGTTGTCTCTGTTTCCACAGATGTGTTTATGCCAGTAAGGTTTGTTAAATTTTTATTAAAACTGGGTTGAGCGTGTCCGGATCAAGTGTGAGCAGCAGAGTTGaa                               <  2:268429/247‑1 (MQ=255)                                                                                                                                                                                                                            gggTGCGACCATGACAGTCTCCTTTTTCTGTATCGTGGAAATCATTTTCATTTTTATTGTTAGCTAATGCAATAGTTACTGAACTGATCCGATGAGTTACTGTTGCACAAATCTCATGTTGCGTGGTGGTCGCTTTTACCACAGATGCGTTTATGCCAGTATGGTTTGTTGAATTTTTATTAAATCTGGGTTGAGCGTGTCGGGAGCAAGTGTGAGCAGCAAAGTGGAACAACTGCGTGCGCAGTTAAATg         >  1:320946/1‑251 (MQ=255)                                                                                                                                                                                                                             ggTGCGACCATGACAGTCTCCTTTTTCTGTATCGTGGAAATCATTTTCATTTTTATTGTTAGCTACTGCCATAGTTACTGAACTGATCCGATGAGTTAATGTTGAACAAATCTCATGTTGCGTGGTGGTCGCTTTTACCACAGATGCGTTTATGCCAGTATGGTTTGTTGAATTTTTATTAAATCTGGGTTGAGCGTGTCGGGAGCAAGTGTGAGAAGCAAAGTGGAAAAACTGCGTGCGCAGTTAAATg         >  2:470682/1‑250 (MQ=255)                                                                                                                                                                                                                                 cGACCATGACAGTCTCCTTTTTCTGTATCGTGGAAATCATTTTCATTTTTATTGTTAGCTAATGCAATAGTTACTGAACTGATCCGATGAGTTAATGTTGAACAAATCTCATGTTGCGTGGTGGTCGCTTTTACCACAGATGCGTTTATGCCAGTATGGTTTGTTGAATTTTTATTAAATCTGGGTTGAGCGTGTCGGGAGCAAGTGTGAGCAGCAAAGTGGAACAACTGCGTGCGCAGTTAAATGAACg     <  2:83048/250‑1 (MQ=255)                                                                                                                                                                                                                                 cGACCATGACAGTCTCCTTTTTCTGTATCGTGGAAATCATTTTCATTTTTATTGTTAGCTAATGCAATAGTTACTGAACTGATCCGATGAGTTAATGTTGAACAAATCTCATGTTGCGTGGTGGTCGCTTTTACCACAGATGCGTTTATGCCAGTATGGTTTGTTGAATTTTTATTAAATCTGGGTTGAGCGTGTCGGGAGCAAGTGTGAGCAGCAAAGTGGAACAACTGCGTGCGCAGTTAAATGAACGt    >  1:66188/1‑251 (MQ=255)                                                                                                                                                                                                                                   aacATGACAGTCTCCTTTTTCTGTATCGTGGAAATCATTTTCATTTTTATTGTTAGCTAAGTCAATAGTTACTGAACTGATCCGATGAGTTAATGTTGAACAAATCTCATGTTGCGTGGTGGTCGCTTTTACCACAGATGCGTTTATGCCAGTATGGTTTGTTGAATTTTTATTAAATCTGGGTTGAGCGTGTCGGGAGCAAGTGTGAGCAGCAAAGTGGAACAACTGCGTGCGCAGTTAAATGAACGTAt  <  2:230476/249‑1 (MQ=255)                                                                                                                                                                                                                                                 |                                                                                                                                                                                                                                             GCCAACCTGACGCACGAAACCCTGCTGTTGCATCGTGGTTAGCAGGCGGTGGGTCGTGGAATTGGGTAACCCGGCTTGTTGCGCCAGTTCCGTGAGTGCCACACTGCCATTGGATTCGGCAATCCACTCCAGTAATTTCAGGCCACGCGTTAAAGACTGAACCTGTCCAGTCGCTGGTGCGGTGGCAACGGCGGGTTTTCTGCCGCGTTTCGCGGGAATGGGTGCGACCATGACAGTCTCTTTTTTCTGTATCGTGGAAATCATTTTCATTTTTATTGTTAGCTAATGCAATAGTTACTGAACTGATCCGATGAGTTAATGTTGAACAAATCTCATGTTGCGTGGTGGTCGCTTTTACCACAGATGCGTTTATGCCAGTATGGTTTGTTGAATTTTTATTAAATCTGGGTTGAGCGTGTCGGGAGCAAGTGTGAGCAGCAAAGTGGAACAACTGCGTGCGCAGTTAAATGAACGTAT  >  NC_000913/4223398‑4223874

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