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%	160.1 / NA	45	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 (10/35);  total (10/35)

CAGTTCGCCAACCTGACGCACGAAACCCTGCTGTTGCATCGTGGTTAGCAGGCGGTGGGTCGTGGAATTGGGTAACCCGGCTTGTTGCGCCAGTTCCGTGAGTGCCACACTGCCATTGGATTCGGCAATCCACTCCAGTAATTTCAGGCCACGCGTTAAAGACTGAACCTGTCCAGTCGCTGGTGCGGTGGCAACGGCGGGTTTTCTGCCGCGTTTCGCGGGAATGGGTGCGACCATGACAGTCTCTTTTTTCTGTATCG‑TGGAAATCA‑TTTTCATTTTTATTGTTAGCTAATGCAATAGTTACTGAACTGATCCGATGAGTTAATGTTG‑AACAAATCTCATGTTGCGTGGTGGTCGCTTTTACCACAGATGCGTTTATGCCAGTATGGTTTGTTGAATTTTTATTAAATCTGGGTTGAGCGTGTCGGGAGCAAGTGTGAGCAGCAAAGTGGAACAACTGCGTGCGCAGTTAAATGAACGTATTCTGGT  >  NC_000913/4223392‑4223880                                                                                                                                                                                                                                                       |                                                                                                                                                                                                                                                      cAGTTCGCCAACCTGACGCACGAAACCCTGCTGTTGCATCGTGGTTAGCAGGCGGTGGGTCGTGGAATTGGGTAACCCGGCTTGTTGCGCCAGTTCCGTGAGTGCCACACTGCCATTGGATTCGGCAATCCACTCCAGTAATTTCAGGCCACGCGTTAAAGACTGAACCTGTCCAGTCGCTGGTGCGGTGGCAACGGCGGGTTTTCTGCCGCGTTTCGCGGGAATGGGTGCGACCATGACAGTCTCctttt                                                                                                                                                                                                                                                   <  2:290761/251‑1 (MQ=255)         cAACCGGACGCACGAAACCCTGCTGTTGCATCGTGGTTAGCAGGCGGTGGGTCGTGGAATTGGGTAACCCGGCTTGTTGCGCCAGTTCCGTGAGTGCCACACTGCCATTGGATTCGGCAATCCACTCCAGTAATTTCAGGCCACGCGTTAAAGACTGAACCTGTCCAGTCGCTGGTGCGGTGGCAACGGCGGGTTTTCTGCCGCGTTTCGCGGGAATGGGTGCGACCATGACAGTCTCCTTTTTCTGTATc                                                                                                                                                                                                                                           <  2:86383/251‑1 (MQ=255)          aaCCGGACGCACGAAACCCTGCTGTTGCATCGTGGTTAGCAGGCGGTGGGTCGTGGAATTGGGTAACCCGGCTTGTTGCGCCAGTTCCGTGAGTGCCACACTGCCATTGGATTCGGCAATCCACTCCAGTAATTTCAGGCCACGCGTTAAAGACTGAACCTGTCCAGTCGCTGGTGCGGTGGCAACGGCGGGTTTTCTGCCGCGTTTCGCGGGAATGGGTGCGACCATGACAGTCTCCTTTTTCTGTATCg                                                                                                                                                                                                                                          <  2:279329/251‑1 (MQ=255)           aCCTGACGCACGAAGCCCTGCTGTTGCATCGTGGTTAGCAGGCGGTGGGTCGTGGAATTGGGTAACCCGGCTGGTTGCGCCAGTTCCGTGAGTGCCACACTGCCATTGGATTCGGCAATCCACTCCAGTAATTTCAGGCCACGCGTTAAAGACTGAACCTGTCCAGTCGCTGGTGCGGTGGCAACGGCGGGTTTTCTGCCGCTGTTCGCGGGAATGGGGGCGACCATGACAGGCGCCTTTTTCTGTATCG‑t                                                                                                                                                                                                                                        <  1:405042/251‑1 (MQ=255)                           ccggctgTTGCATCGTGGTTAGCAGGCGGTGGGTCGTGGAATTGGGTAACCCGGCTTGTTGCGCCAGTTCCGTGAGTGCCACACTGCCATTGGATTCGGCAATCCACTCCAGTAATTTCAGGCCACGCGTTAAAGACTGAACCTGTCCAGTCGCTGGTGCGGTGGCAACGGCGGGTTTTCTGCCGCGTTTCGCGGGAATGGGTGCGACCATGACAGTCTCCTTTTTCTGTATCG‑TGGa                                                                                                                                                                                                                                     <  2:127727/235‑1 (MQ=255)                           ccTGCTGTTGCATCGTGGTTAGCAGGCGGTGGGTCGTGGAATTGGGTAACCCGGCTTGTTGCGCCAGTTCCGTGAGTGCCACACTGCCATTGGATTCGGCAATCCACTCCAGTAATTTCAGGCCACGCGTTAAAGACTGAACCTGTCCAGTCGCTGGTGCGGTGGCAACGGCGGGTTTTCTGCCCCGTTTCGCGGGGATGGGTGCGGCCATGGCAGTCTCCTTTTTTTGTATTG‑Tgggagaacgggagagg                                                                                                                                                                                                                        >  1:127727/1‑237 (MQ=255)                             tgctgTTGCATCGTGGTTAGCAAGCGGTGGGTCGTGGAATTGGGTAACCCGGCTTGTTGCGCCAGTTCCGTGAGTGCCACACTGCCATTGGATTCGGCAATCCACTCCAGTCATTTCAGGCCACGCGTTAAAGACTGAACCTGTCCAGTCGCTGGTGCGGTGGCAACGGCGGGGTTTTTTCCGCGTTTCGCGGGGATGGGTGCGGCCCTGAACGTCTCCTTTTTTTTTATCGGTGGAA‑tca‑tttttttttt                                                                                                                                                                                                                     >  2:48520/1‑245 (MQ=255)                              gctgTTGCATCGTGGTTAGCAGGCGGTGGGTCGTGGAATTGGGTAACCCGGCTTGTTGCGCCAGTTCCGTGAGTGCCACACTGCCATTGGATTCGGCAATCCACTCCAGTAATTTCAGGCCACGCGTTAAAGACTGAACCTGTCCAGTCGCTGGTGCGGTGGCAACGGCGGGTTTTCTGCCGCGTTTCGCGGGAATGGGTGCGACCATGACAGTCTCCTTTTTCTGTATCG‑TGGAAAtca‑ttttcattt                                                                                                                                                                                                                      <  1:265766/249‑1 (MQ=255)                                   tgcatagaggTTATCAGGCGGTGGGACGTGGTCTTGGGTAACCCGGCTTGTTGCGCCAGTTCCGTGAGTGCCACACTGCCATTGGATTCGGCAATCCACTCCAGTAATTTCAGGCCCCGCGTTAAAGACTGAACCTGTCCAGTAGCTGGTGCGGTGGCCACGGCGGGTTTTCTGCCCAGTTTTGCGGGGAGTGGTGGGACCATGACAGCCTCCCTTTTCTTTATCG‑TGGGAAtca‑tttttctttttttttt                                                                                                                                                                                                               >  2:33553/9‑239 (MQ=255)                                                    gCGGTGGGTCGTGGAATTGGGTAACCCGGCTTGTTGCGCCAGTTCCGTGAGTGCCACACTGCCATTGGATTCGGCAATCCACTCCAGTAATTTCAGGCCACGCGGTAAAGACTGAACCTGTCCAGTCGCTGGTGCGGTGGCAACGGCGGGTTTTCTGCCGCGTTTCGCGGGAATGGGTGCGACCATGACAGGCTCCTTTTGCTGTATCG‑TGGAAATCA‑TTTTCATTTTTATTGTTAGCTAATGCAATAGtt                                                                                                                                                                                              <  1:55571/251‑1 (MQ=255)                                                         gggTCGTGGAATTGGGTAACCCGGCTTGTTGCGCCAGTTCCGTGAGTGCCACACTGCCATTGGATTCGGCAATCCACTCCAGTAATTTCAGGCCACGCGTTAAAGACTGAACCTGTCCAGTCGCTGGTGCGGTGGCAACGGCGGGTTTTCTGCCGCGTTTCGCGGGAATGGGTGCGACCATGACAGTCTCCTTTTTCTGTATCG‑TGGAAATCA‑TTTTCAGTTTTATTGTTAGCTAATGCAATAGTTactga                                                                                                                                                                                         <  1:407502/251‑1 (MQ=255)                                                                ggAATTGGGTAACCCGGCTTGTTGCGCCAGTTCCGTGAGTGCCACACTGCCATTGGATTCGGCAATCCACTCCAGTAATTTCAGGCCACGCGTTAAAGACTGAACCTGTCCAGTCGCTGGTGCGGTGGCAACGGCGGGTTTTCTGCCGCGTTTCGCGGGAATGGGTGCGACCATGACAGTCTCCTTTTTCTGTATCG‑TGTAAATCA‑TTTTCATTTTTATTGTTAGCTAATGCAATAGTTACTGAACTGAt                                                                                                                                                                                   <  1:439366/250‑1 (MQ=255)                                                                 gAATTGGGTAACCCGGCTTGTTGCGCCAGTTCCGTGAGTGCCACACTGCCATTGGATTCGGCAATCCACTCCAGTAATTTCAGGCCACGCGTTAAAGACTGAACCTGTCCAGTCGCTGGTGCGGTGGCAACGGCGGGTTTTCTGCCGCGTTTCGCGGGAATGGGTGCGACCATGACAGTCTCCTTTTTCTGTATCG‑TGGAAATCA‑TTTTCATTTTTATTGTTAGCTAATGCAATAGTTACTGAACTGATcc                                                                                                                                                                                 <  1:227636/251‑1 (MQ=255)                                                                 gAATTGGGTAACCCGGCTTGTTGCGCCAGGTCCGTGAGTGCCACACTGCCATTGGATTCGGCAATCCACTCCAGTAATTTCAGGCCACGCGTTAAAGACTGAACCTGTCCAGTCGCTGGTGCGGTGGCAACGGCGGGTTTTCTGCCGCGTTTCGCGGGAATGGGTGCGACCATGACAGTCTCCTTTTTCTGTATCG‑TGGAAATCA‑TTTTCATTGTTATTGTTAGCTAATGCAATAGTTACTGAACTGATcc                                                                                                                                                                                 <  2:330135/251‑1 (MQ=255)                                                                   aTTGGGTAACCCGGCTTGTTGCGCCAGTTCCGTGAGTGCCACACTGCCATTGGATGCGGCAATCCACTCCAGTAATTTCAGGCCACGCGTTAAAGACTGAACCTGTCCAGTCGCTGGTGCGGTGGCAACGGCGTGTTGTCTGCCGCGTTTCGCGGGAATGGGTGCGACCATGACAGTCTCCTTTTGCTGTATCG‑TGGAAATCA‑TTTTCATTTTTATTGTTAGCTAATGCAATAGTTACTGAACTGATCCGa                                                                                                                                                                               <  1:185643/251‑1 (MQ=255)                                                                              cGGCTTGTTGCGCCAGTTCCGTGATTGCCACACTGCCATTGGATTCGGCAATCCACTCCAGTAATTTCAGGCCACGCGTTAAAGACTGAACCTGTCCAGTCGCTGGTGCTGTGGCAACGGCGGGTTTTCTGCCGCGTCTCGCGGGAAGGAGTGCGACCATGACAGTCTCCTTTCTCTGTATCG‑TTGAAATCA‑TTTTCATTTGTATTGTTAGCTAATGCAATAGTTACTGAACTGATCCGATGAGTTAATGt                                                                                                                                                                    <  1:424884/251‑1 (MQ=255)                                                                                   tgttgCGCCAGTTCCGTGAGTGCCACACTGCCATTGGATTCGGCAATCCACTCCAGTAATTTCAGGCCACGCGTTAAAGACTGAACCTGTCCAGTCGCTGGTGCGGTGGCAACGGCGGGTTTTCTGCCGCGTTTCGCGGGAATGGGTGCGACCATGACAGTCTCCTTTTTCTGTATCG‑TGGAAATCA‑TTTTCATTTTTATTGTTAGCTAATGCAATAGTTACTGAACTGATCCGATGAGTTAATGTTG‑aa                                                                                                                                                               <  1:222119/250‑1 (MQ=255)                                                                                                ccGTGAGTGCCACACTGCCATTGGATTCGGCAATCCACTCCAGTAATTTCAGGCCACGCGTTAAAGACTGAACCTGTCCAGTCGCTGGTGCGGTGGCAACGGCGGGTTTTCTGCCGCGTTTCGCGGGAATGGGTGCGACCATGACAGTCTCCTTTTTCTGTATCG‑TGGAAATCA‑TTTTCATTTTTATTGTTAGCTAATGCAATAGTTACTGAACTGATCCGATGAGTTAATGTTG‑AACAAATCTCATGTTg                                                                                                                                                 <  1:231315/251‑1 (MQ=255)                                                                                                  gTGAGTGCCACACTGCCATTGGATTCGGCAATCCACCCCAGTAATTTCAGGCCACGCGTTAAAGACTGAACCTGTCCAGTCGCTGGTGCGGTGGCAACGGCGTGTTTTCTGCCGCGTTTCGCGGGAATGGGTGCGACCATGACAGTCTCCTTTTTCTGTATCG‑TGGAAATCA‑TTTTCATTTTTATTGTTAGCTAAGGCAATAGTTACTGAACTGATCCGATGAGTTa                                                                                                                                                                        <  1:150715/227‑1 (MQ=255)                                                                                                       tGCCACACTGCCATTGGATTCGGCAATCCACTCCAGTAATTTCAGGCCACGCGTTAAAGACTGAACCTGTCCAGTCGCTGGTGCGGTGGCAACGGCGGGGTTTCTGCCGCGGTTCGCGGGAATGGGTGCGGCCATGACAGTCTCCCTTTTTTGTATCG‑TGGAAATCA‑TTTTTCTTTTTTTTGTTAGCTAATTCAATTGTTACTGAACTGATCCGATGAGTTAATGTTGTACCAAATCACA‑GTTTCgtggtg                                                                                                                                          >  1:248238/1‑251 (MQ=255)                                                                                                             acTGCCATTGGATTCGGCAATCCACTCCAGTAATTTCAGGCCACGCGTTAAAGACTGAACCTGTCCAGTCGCTGGTGCGGTGGCAACGGCGGGTTTTCTGCCGCGTTTCGCGGGAATGGGTGCGACCATGACAGTCTCCTTTTTCTGTATCG‑TGGAAATCA‑TTTTCATTTTTATTGTTAGCTAATGCAATAGTTACTGAACTGATCCGATGAGTTAATGTTG‑AACAAATCTCATGTTGCGTGGTGGTCg                                                                                                                                      <  2:389365/249‑1 (MQ=255)                                                                                                             acTGCCATTGGATTCGGCAATCCACTCCAGTAATTTCAGGCCACGCGTTAAAGACTGAACCTGTCCAGTCGCTGGTGCGGTGGCAACGGCGGGGTTTCTGCCGCGTTTCGCGGGGATGGGGGCGACCATGAAAGTCTCCTTTTTCTGTATTG‑TGGAAATCA‑TTTTTTTTTTTATTTTTAGCTAATTCAATAGTTACTGAACTGATCCGATGGGTTAATGTTG‑AACAAATCTCATGTTGCGTGGTGGTCggg                                                                                                                                    >  1:389365/1‑249 (MQ=255)                                                                                                                gCCATTGGATTCGGCAATCCACTCCAGTAATTTCAGGCCACGCGTTAAAGACTGAACCTGTCCAGTCGCTGGTGCGGTGGCAACGGCGGGTTTTCTGCCGCGTTTCGCGGGAATGGGTGCGACCATGACAGTCTCCTTTTTCTGTATCG‑TGGAAATCA‑TTTTCATTTTTATTGTTAGCTAATGCAATAGTTACTGAACTGATCCGATGAGTTAATGTTG‑AACAAATCTCATGTTGCGTGGTGGTCGCttt                                                                                                                                  <  1:212691/250‑1 (MQ=255)                                                                                                                 ccATTGGATTCGGCAATCCACTCCAGTAATTTCAGGCCACGCGTTAAAGACTGAACCTGTCCAGTCGCTGGTGCGGTGGCAACGGCGGGTTTTCTGCCGCGTTTCGCGGGAATGGGTGCGACCATGACAGTCTCCTTTTTCTGTATCG‑TGGAAATCA‑TTTTCATTTTTATTGTTAGCTAATGCAATAGTTACTGAACTGATCCGATGAGTTAATGTTG‑AACAAATCTCATGTTGCGTGGTGGTCGCTTTTa                                                                                                                                <  1:405506/251‑1 (MQ=255)                                                                                                                     tGGATTCGGCAATCCACTCCAGTAATTTCAGGCCACGCGTTAAAGACTGAACCTGTCCAGTCGCTGGTGCGGTGGCAACGGCGGGTTTTTTGCCGCGGTTCGCGGGAATGGGTGCGACCCTGGCAGCTTCCTTTTTCTGTATCG‑TGGGAATCA‑TTTTCATTTTTTTTGTTAGCTCATGCAATAGTTTATGAACTGATCCGATGAGTTTATGTTT‑AACAAATCTCATGTTGCGTGGTGGTCGCTTTTTCcac                                                                                                                            >  1:69079/1‑251 (MQ=255)                                                                                                                      ggATTCGGCAATCCACTCCAGTAATTTCAGGCCACGCGTTAAAGACTGAACCTGTCCAGTCGCTGGTGCGGTGGCAACGGCGGGTTTTCTGCCGCGTTTCGCGGGAATGGGTGCGACCATGACAGTCTCCTTTTTCTGTATCG‑TGGAAAGCA‑TTTTCATTTTTATTGTTAGCTAATGCAATAGTTACTGAACTGATCCGATGAGTTAATGTTG‑AACAAATCTCATGTTGCTTGGTGGTCGCTTTTACcaca                                                                                                                           <  2:442809/251‑1 (MQ=255)                                                                                                                                                     ccACGCGTTAAAGACTGAACCTGTCCAGTCGCTGGTGCGGTGGCAACGGCGGGTTTTCTGCCGCGTTTCGCGGGAATGGGTGCGACCATGACAGTCTCCTTTTTCTGTATCG‑TGGAAATCA‑TTTTCATTTTTATTGTTAGCTAATGCAATAGTTACTGAACTGATCCGATGAGTTAATGTTG‑AACAAATCTCATGTTGCGTGGTGGTCGCTTTTACCACAGATGCGTTTATGCCAGTATGGTTTGTTGAAt                                                                                            <  1:258507/251‑1 (MQ=255)                                                                                                                                                              aaaGACTGAACCTGTCCAGTCGCTGGTGCGGTGGCAACGGCGGGTTTTCTGCCGCGTTTCGCGGGAATGGGTGCGACCATGACAGTCTCCTTTTTCTGTATCG‑TGGAAATCA‑TTTTCATTTTTATTGTTAGCTAATGCAATAGTTACTGAACTGATCCGATGAGTTAATGTTG‑AACAAATCTCATGTTGCGTGGTGGTCGCTTTTACCACAGATGCGTTTATGCCAGTATGGTTTGTTGAATTTTTATTaa                                                                                   <  1:237246/251‑1 (MQ=255)                                                                                                                                                                aGACTGAACCTGTCCAGTCGCTGGTGCGGTGGCAACGGCGGGTTTTCTGCCGCGTTTCGCGGGAATGGGTGCGACCATGACAGTCTCCTTTTTCTGTATCG‑TGGAAATCA‑TTTTCATTTTTATTGTTAGCTAATGCAATAGTTACTGAACTGATCCGATGAGTTAATGTTG‑AACAAATCTCATGTTGCGTGGTGGTCGCTTTTACCACAGATGCGTTTATGCCAGTATGGTTTGTTGAATTTTTATTAAAt                                                                                 <  1:198685/251‑1 (MQ=255)                                                                                                                                                                              cAGTCGCTGGTGCGGTGGCAACGGCGGGTTTTCTGCCGCGTTTCGCGGGAATGGGTGCGACCATGACAGTCTCCTTTTTCTGTATCG‑TGGAAATCA‑TTTTCATTTTTATTGTTAGCTAATGCAATAGTTACTGAACTGATCCGATGAGTTAATGTTG‑AACAAATCTCATGTTGCGTGGTGGTCGCTTTTACCACAGATGCGTTTATGCCAGTATGGTTTGTTGAATTTTTATTAAATCTGGGTTGAGCgtg                                                                   <  1:435496/251‑1 (MQ=255)                                                                                                                                                                               gggcGCTGGTGCGGTGTCAACGGCGGTTTTTCTGCCTCGTTTCTCGGGAATGGGTTCGACCATGACATTCTCCTTTGTTTGTATCG‑TGTAAATCA‑TTTTCATTTTTATTGTTAGCTAATGCAATAGTTACTGAACTGATCCGATGTGTGAATGTTG‑AACAAAGCTCATGTTGCGTGGTTGTCGCTTTTACCACAGATGCGTTTATGCCAGTATGGTTTGTTTAATTTTTATTAAATCTGGTTTGAGCgtga                                                                  <  2:138496/248‑2 (MQ=255)                                                                                                                                                                                         gCGGTGGCAACGGCGGGTTTTCTGCCGCGTTTCGCGGGAATGGGTGCGACCATGACAGTCTCCTTTTTCTGTATCG‑TGGAAATCA‑TTTTCATTTTTATTGTTAGCTAATGCAATAGTTACTGAACTGATCCGATGAGTTAATGTTG‑AACAAATCTCATGTTGCGTGGTGGTCGCTTTTACCACAGATGCGTTTATGCCAGTATGGTGTGTTGAAGTTTTATTAAATCTGGGTTGAGCGTGTCTGGAGCaa                                                         <  1:174925/250‑1 (MQ=255)                                                                                                                                                                                           ggTGGCAACGGCGGGTTTTCTGCCGCGTTTCGCGGGAATGGGGGCGCCCACGACAGGCTCCTTTTTCTGTATCG‑TGGAAAACAATTTTC‑TTTTTATTGTTAGCTAATGCAATAGTTAATGAACTGAACCGCTGAGATAATGTTG‑AACAAATCTCATGTTTCGTGGGGGTCGCTTTTACCCCAAATGCGTTTATGCCAGTATGGTTTGTTGGATTTTTATTTAATCTGGGTTGAGCGGGTCGGGGGCAAgtg                                                      >  1:52312/1‑251 (MQ=255)                                                                                                                                                                                                     gcggGTTTTCTGCCGCGTTTCGCGGGAATGGGTGCGACCATGACAGTCTCCTTTTTCTGTATCG‑TGGAAATCA‑TTTTCATTTTTATTGTTAGCTAATGCAATAGTTACTGAACTGATCCGATGAGTTAATGTTG‑AACAAATCTCATGTTGCGTGGTGGTCGCTTTTACCACAGATGCGTTTATGCCAGTATGGTTTGTTGAATTTTTATTAAATCTGGGTTGAGCGTGTCGGGAGCAAGTGTGAgcagca                                             >  1:128264/1‑250 (MQ=255)                                                                                                                                                                                                       gggTTTTCTGCCGCGTTTCGCGGGAATGGGTGCGACCATGACAGTCTCCTTTTTCTGTATCG‑TGGAAATCA‑TTTTCATTTTTATTGTTAGCTAATGCAATAGTTACTGAACTGATCCGATGAGTTAATGTTG‑AACAAATCTCATGTTGCGTGGTGGTCGCTTTTACCACAGATGCGTTTATGCCAGTATGGTTTGTTGAATTTTTATTAAATCTGGGTTGAGCGTGTCGGGAGCAAGTGTGAGCAGCAAAg                                          <  1:59506/251‑1 (MQ=255)                                                                                                                                                                                                           tttCTGCCGCGTTTCGCGGGAATGGGTGCGACCATGACAGTCTCCTTTTTCTGTATCG‑TGGAAAGCA‑TTTGCATTTTTATTGTTAGCTAATGCAATAGTTACTGAACTGATCCGATGAGTTAATGTTG‑AACAAATCTCATGTTGCGTGGTGGTCGCTTTTACCACAGATGCGTTTATGCCAGTATGGTTTGTTGAATTTTTATTAAATCTGGGGTGAGCGTGTCGGGAGCAAGTGTGAGCAGCAAAGTgg                                       <  2:390012/250‑1 (MQ=255)                                                                                                                                                                                                                gCCGCGTTTCGCGGGAATGGGTGCGACCATGACAGTCTCCTTTTTCTGTATCG‑TGGAAATCA‑TTTTCATTTTTATTGTTAGCTAATGCAATAGTTACTGAACTGATCCGATGAGTTAATGTTG‑AACAAATCTCATGTTGCGTGGTGGTCGCTTTTACCACAGATGCGTTTATGCCAGTATGGTTTGTTGAATTTTTATTAAATCTGGGTTGAGCGTGTCGGGAGCAAGTGTGAGCAGCAAAGTGGaacaac                                 <  2:463392/251‑1 (MQ=255)                                                                                                                                                                                                                 ccGCGTTTCGCGGGATTGGGTGCGACCATGACAGTCTCCTTTTGCTGGTTCG‑TGGAAATCA‑TTTTCATTGTTATTGTGAGCTAATGCAATAGTTACTGAACTGATCCGATGAGTTAATTTTG‑AACAAATCTCAGTTTGCGTGGTGGTCGCTTTTGCCATAGTGGCGTTTATGCCAGTATGGTTTGTTGAATGTTGATTAAATTTGTGGTGATCGGGTCGTGATCAAGTGGGAGCAGCATAGTGTAACAACt                                <  2:454027/251‑1 (MQ=255)                                                                                                                                                                                                                           cgGGAATGGGTGCGAACATGACAGTCTCCTTTTTCTGTATCG‑TGGAAATCA‑TTTTCATTTTTATTGTTAGCTAATGCAATAGTTACTGAACTGATCCTATGAGTTAATGTTG‑CACAAATCTCATGTTGCGTGGTGGTCGCTTTTCCCACAGATGCGTTTATGCCAGTATGGTTTGTTGAATTTTTATTAAATCTGGGTTGAGCGTGTCGGGAGCAAGTGTGAGAAGCAAAGTGGAACAACTGCGTgcgc                        >  2:262510/1‑249 (MQ=255)                                                                                                                                                                                                                              gAATGGGTGCGACCATGATAGTCTCCTTTTTCTTTATCG‑TGGAAATCA‑TTTTCATTTTTATTGTTAGCTAATGCAATAGTTACTGAACTGATCCGATGAGTTAATGTTG‑AACAAATCTCATGTTGCGTGGTGGTCGCTTTTACCACAGATGCGTTTATGCTAGTATGGTTTGTTGAATTTTTATTAAATCTGGGTTGAGCGTGTCGGTAGCAAGTTTGAGCAGTAAAGTGGATCAACTGCGTGGGCAGTTa                   <  2:266208/251‑1 (MQ=255)                                                                                                                                                                                                                                aTGGGTGCGACCATGACAGTCTCCTTTTTCTGTATCG‑TGGAAATCA‑TTTTCATTTTTATTGTTAGCTAATGCAATAGTTACTGAACTGATCCGATGAGTTAATGTTG‑AACAAATCTCATGTTGCGTGGTGGTCGCTTTTACCACAGATGCGTTTATGCCAGTATGGTTTGTTGAATTTTTATTAAATCTGGGTTGAGCGTGTCGGGAGCAAGTGTGAGCAGCAAAGTGGAACAACTGCGTGCGCa                       >  2:97666/1‑245 (MQ=255)                                                                                                                                                                                                                                aTGGGTGCGACCATGACAGTCTCCTTTTTCTGTATCG‑TGGAAATCA‑TTTTCATTTTTATTGTTAGCTAATGCAATAGTTACTGAACTGATCCGATGAGTTAATGTTG‑AACAAATCTCATGTTGCGTGGTGGTCGCTTTTACCACAGATGCGTTTATGCCAGTATGGTTTGTTGAATTTTTATTAAATCTGGGTTGAGCGTGTCGGGAGCAAGTGTGAGCAGCAAAGTGGAACAACTGCGTGCGCa                       <  1:97666/245‑1 (MQ=255)                                                                                                                                                                                                                                  gggTGCGACTATGACAGTATCCTTTGTCTGTATCG‑TGGAAATTA‑TTGTCAGTTTTATTGTTATCTAATGCAATAGTTACTGATCTGATCCGATGAGTTAATGTTG‑AACAAATCTCATGGTGCGTGGTGGTCGCTTTTACCACAGAGGCGGTTATGCCAGTATGGGTTGTTGAATTTTTATTAAATCTGGGTTGAGCGTGTCGGGAGCAAGTGTGAGCAGCAAAGTGGAACAACTGCGTGCGCAGTTAAAt                <  1:448367/250‑1 (MQ=255)                                                                                                                                                                                                                                  gggTGCGACCATGACAGTCTCCTTTTTCTGTATCG‑TGGAAATCA‑TTTTCATTTTTATTGTTAGCTAATGCAATAGTTACTGAACTGATCCGATGAGTTAATGTTG‑AACAAATCTCATGTTGCGTGGTGGTCGCTTTTACCACAGATGCGTTTATGCCAGTATGGTTTGTTGAATTTTTATTAAATCTGGGTTGAGCGTGTCGGGAGCAAGTGTGAGCATCAAAGTGGAACAACTGCGTGCGCAGTTAAATg               <  1:389524/251‑1 (MQ=255)                                                                                                                                                                                                                                      gCGACCATGACAGTCTCCTTTTTCTGTATCG‑TGGAAATCA‑TTTTCATTTTTATTGTTAGCTAATGCAATAGTTACTGAACTGATCCGATGAGTTAATGTTG‑AACAAATCTCATGTTGCGTGGTGGTCGCTTTTACCACAGATGCGTTTATGCCAGTATGGTTTGTTGAATTTTTATTAAATCTGGGTTGAGCGTGTCGGGAGCAAGTGTGAGCAGCAAAGTGGAACAACTGCGTGCGCAGTTAAATGAACg           <  1:407147/251‑1 (MQ=255)                                                                                                                                                                                                                                               aCAGTCTCCTTTTTCTGTATCG‑TGGAAATCA‑TTTTCATTTTTATTGTTAGCTAATGCAATAGTTACTGAACTGATCCGATGAGTTAATGTTG‑AACAAATCTCATGTTGCGTGGTGGTCGCTTTTACCACAGATGCGTTTATGCCAGTATGGTTTGTTGAATTTTTATTAAATCTGGGTTGAGCGTGTCGGGAGCAAGTGTGAGCAGCAAAGTGGAACAACTGCGTGCGCAGTTAAATGAACGTATTCTGGt  <  1:48520/251‑1 (MQ=255)                                                                                                                                                                                                                                                       |                                                                                                                                                                                                                                                      CAGTTCGCCAACCTGACGCACGAAACCCTGCTGTTGCATCGTGGTTAGCAGGCGGTGGGTCGTGGAATTGGGTAACCCGGCTTGTTGCGCCAGTTCCGTGAGTGCCACACTGCCATTGGATTCGGCAATCCACTCCAGTAATTTCAGGCCACGCGTTAAAGACTGAACCTGTCCAGTCGCTGGTGCGGTGGCAACGGCGGGTTTTCTGCCGCGTTTCGCGGGAATGGGTGCGACCATGACAGTCTCTTTTTTCTGTATCG‑TGGAAATCA‑TTTTCATTTTTATTGTTAGCTAATGCAATAGTTACTGAACTGATCCGATGAGTTAATGTTG‑AACAAATCTCATGTTGCGTGGTGGTCGCTTTTACCACAGATGCGTTTATGCCAGTATGGTTTGTTGAATTTTTATTAAATCTGGGTTGAGCGTGTCGGGAGCAAGTGTGAGCAGCAAAGTGGAACAACTGCGTGCGCAGTTAAATGAACGTATTCTGGT  >  NC_000913/4223392‑4223880

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