Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	NC_000913	100,946	T→C	M61T (ATG→ACG)	murC →	UDP‑N‑acetylmuramate‑‑L‑alanine ligase

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	100,946	0	T	C	100.0%	161.0 / NA	46	M61T (ATG→ACG)	murC	UDP‑N‑acetylmuramate‑‑L‑alanine ligase
Reads supporting (aligned to +/- strand):  ref base T (0/0);  new base C (22/24);  total (22/24)

GCCCGGGCGTAATTGTAGCGATGCCTTTTGCATCGTATGAATTTAAGAAGTTAATGGCGTAAAGAATGAATACACAACAATTGGCAAAACTGCGTTCCATCGTGCCCGAAATGCGTCGCGTTCGGCACATACATTTTGTCGGCATTGGTGGTGCCGGTATGGGCGGTATTGCCGAAGTTCTGGCCAATGAAGGTTATCAGATCAGTGGTTCCGATTTAGCGCCAAATCCGGTCACGCAGCAGTTAATGAATCTGGGTGCGACGATTTATTTCAACCATCGCCCGGAAAACGTACGTGATGCCAGCGTGGTCGTTGTTTCCAGCGCGATTTCTGCCGATAACCCGGAAATTGTCGCCGCTCATGAAGCGCGTATTCCGGTGATCCGTCGTGCCGAAATGCTGGCTGAGTTAATGCGTTTTCGTCATGGCATCGCCATTGCCGGAACGCACGGCAAAACGACAACCACCGCGATGGTTTCCAGCATCTAC  >  NC_000913/100700‑101187                                                                                                                                                                                                                                                       |                                                                                                                                                                                                                                                  gCCCGGGCGTAATTGTAGCGATGCCTTTTGCATCGTATGAATTTAAGAAGTTAATGGCGTAAAGAATGAATACACAACAATTGGCAAAACTGCGTTCCATCGTGCCCGAAATGCGTCGCGTTCGGCACATACATTTTGTCGGCATTGGTGGTGCCGGTATGGGCGGTATTGCCGAAGTTCTGGCCAATGAAGGTTATCAGATCAGTGGTTCCGATTTAGCGCCAAATCCGGTCACGCAGCAGTTAACGAAt                                                                                                                                                                                                                                               >  1:125471/1‑251 (MQ=255) gCCCGGGCGTAATTGTAGCGATGCCTTTTGCATCGTATGAATTTAAGAAGTTAATGGCGTAAAGAATGAATACACAACAATTGGCAAAACTGCGTTCCATCGTGCCCGAAATGCGTCGCGTTCGGCACATACATTTTGTCGGCATTGGTGGTGCCGGTATGGGCGGTATTGCCGAAGTTCTGGCCAATGAAGGTTATCAGATCAGTGGTTCCGATTTAGCGCCAAATCCGGTCACGCAGCAGTTAACGAAt                                                                                                                                                                                                                                               >  2:262431/1‑251 (MQ=255) gCCCGGGCGTAATTGTAGCGATGCCTTTTGCATCGTATGAATTTAAGAAGTTAATGGCGTAAAGAATGAATACACAACAATTGGCAAAACTGCGTTCCATCGTGCCCGAAATGCGTCGCGTTCGGCACATACATTTTGTCGGCATTGGTGGTGCCGGTATGGGCGGTATTGCCGAAGTTCTGGCCAATGAAGGTTATCAGATCAGTGGTTCCGATTTAGCGCCAAATCCGGTCACGCAGCAGTTAACGAAt                                                                                                                                                                                                                                               <  1:235152/251‑1 (MQ=255)               gTAGCGATGCCTTTTGCATCGTATGAATTTAAGAAGTTAATGGCGTAAAGAATGAATACACAACAATTGGCAAAACTGCGTTCCATCGTGCCCGAAATGCGTCGCGTTCGGCACATACATTTTGTCGGCATTGGTGGTGCCGGTATGGGCGGTATTGCCGAAGTTCTGGCCAATGAAGGTTATCAGATCAGTGGTTCCGATTTAGCGCCAAATCCGGTAACGCAGCAGTTAACGAATCTGGGTGCGACGat                                                                                                                                                                                                                                 >  2:12358/1‑251 (MQ=255)                              gCATCGTATGAATTTAAGAAGTTAATGGCGTAAAGAATGAATACACAACAATTGGCAAAACTGCGTTCCATCGTGCCCGAAATGCGTCGCGTTCGGCACATACATTTTGTCGGCATTGGTGGTGCCGGTATGGGCGGTATTGCCGAAGTTCTGGCCAATGAAGGTTATCAGATCAGTGGTTCCGATTTAGCGCCAAATCCGGTCACGCAGCAGTTAACGAATCTGGGTGCGACGATTTATTTCAACCATCg                                                                                                                                                                                                                  <  2:76790/251‑1 (MQ=255)                                     aTGAATTTAAGAAGTTAATGGCGTAAAGAATGAATACACAACAATTGGCAAAACTGCGTTCCATCGTGCCCGAAATGCGTCGCGTTCGGCACATACATTTTGTCGGCATTGGTGGTGCCGGTATGGGCGGTATTGCCGAAGTTCTGGCCAATGAAGGTTATCAGATCAGTGGTTCCGATTTAGCGCCAAATCCGGTCACGCAGCAGTGAACGAATCTGGGTGCGACGATTTATTTCAACCATCGCCCGGaa                                                                                                                                                                                                           <  1:258841/251‑1 (MQ=255)                                        aaTTTAAGAAGTTAATGGCGTAAAGAATGAATACACAACAATTGGCAAAACTGCGTTCCATCGTGCCCGAAATGCGTCGCGTTCGGCACATACATGTTGTCGGCATTGGTGGTGCCGGTATGGGCGGTATTGCCGAAGTTCTGGCCAATGAAGGTTATCAGATCAGTGGTTCCGATTTAGCGCCAAATCCGGTCACGCAGCAGTTAACGAATCTGGGTGCGACGATTTATTTCAACCATCGCCCGGAAAac                                                                                                                                                                                                        <  2:89711/251‑1 (MQ=255)                                             aagaagTTAATGGCGTAACGAATGAATACACAACAATTGGCAAAACTGCGTTCCATCGTGCCCGAAATGCGTCGCGTTCGGCACATACATTTTGTCGGCATTGGTGGTGCCGGTATGGGCGGTATTGCCGAAGTTCTGGCCAATGAAGGTTATCAGATCAGTGGTTCCGATTTAGCGCCAAATCCGGTCACGCAGCAGTTAACGAATCTGGGTGCGACGATTTATTTCAACCATCGCCCGGAAAacgtacg                                                                                                                                                                                                   >  2:305290/1‑251 (MQ=255)                                               gaagTTAATGGCGTAAAGAATGAATACACAACAATTGGCAAAACTGCGTTCCATCGTGCCCGAAATGCGTCGCGTTCGGCACATACATTTTGTCGGCATTGGTGGTGCCGGTATGGGCGGTATTGCCGAAGTTCTGGCCAATGAAGGTTATCAGATCAGTGGTTCCGATTTAGCGCCAAATCCGGTCACGCAGCAGTTAACGAATCTGGGTGCGACGATTTATTTCAACCATCGCCCGGAAAACGTACGTg                                                                                                                                                                                                 <  1:258439/251‑1 (MQ=255)                                                 agTTAATGGCGTAAAGAATGAATACACAACAATTGGCAAAACTGCGTTCCATCGTGCCCGAAATGCGTCGCGTTCGGCACATACATTTTGTCGGCATTGGTGGTGCCGGTATGGGCGGTATTGCCGAAGTTCTGGCCAATGAAGGTTATCAGATCAGTGGTTCCGATTTAGCGCCAAATCCGGTCACGCAGCAGTTAACGAATCTGGGTGCGACGATTTATTTCAACCATCGCCCGGAAAACGTACGTGAt                                                                                                                                                                                               <  2:315937/251‑1 (MQ=255)                                                  gttaataccGTAAAGAATGAATACACAACAATTGGCAAAACTGCGTTCCATCGTGCCCGAAATGCGTCGCGTTCGGCACATACATTTTGTCGGCATTGGTGGTGCCGGTATGGGCGGTATTGCCGAAGTTCTGGCCAATGAAGGTTATCAGATCAGTGGTTCCGATTTAGCGCCAAATCCGGTCACGCAGCAGTTAACGAATCTGGGTGCGACGATTTATTTCAACCATCGCCAGGAAAacgtacgtactg                                                                                                                                                                                              >  2:272257/9‑247 (MQ=255)                                                   ttAATGGCGTCAAGAATGAATACACAACAATTGGCAAAACTGCGTTCCATCGTGCCCGAAATGCGTCGCGTTCGGCACATACATTTTGTCGGCATTGGTGGTGCCGGTATGGGCGGTATTGCCGACGTTCTGGCCAATGAAGGTTATCAGATCAGTGGTTCCGATTTAGCGCCAAATCCGGTCACGCAGCAGTTAACGAATCTGGGTGCGACGATTTATTTCAACCATCGCCCGGAAAACGTACGTGATGc                                                                                                                                                                                             >  2:309583/1‑251 (MQ=255)                                                         gCGTAAAGAATGAATACCCAACAATTGGCAAAACTGCGTTCCATCGTGCCCGAAATGCGTCGCGTTCGGCACATACATTTTGTCGGCATTGGTGGTGCCGGTATGGGCGGTATTGCCGAAGTTCTGGCCAATGAAGGTTATCAGATCAGTGGTTCCGATTTAGCGCCAAATCCGGTCACGCAGCAGTTAACGAATCTGGGTGCGACGATTTATTTCAACCATCGCCCGGAAAACGTACGTGATGCCAGCGt                                                                                                                                                                                       <  2:88552/251‑1 (MQ=255)                                                                 aatgaatACACAACAATTGGCAAAACTGCGTTCCATCGTGCCCGAAATGCGTCGCGTTCGGCCCATACATTTTGTCGGCATTGGTGGTGCCGGTATGGGCGGTATTGCCGAAGTTCTGGCCAATGAAGGTTATCAGATCAGTGGTTCCGATTTAGCGCCAAATCCGGTCACGCCGCAGTTAACGAATCTGGGTGCGACGATTTATTTCAACCATCGCCCGGACAACGTACGTGATGCCCGCGTGGTCgttg                                                                                                                                                                               >  2:136408/1‑251 (MQ=255)                                                                    gaatACACAACAATTGGCAAAACTGCGTTCCATCGTGCCCGAAATGCGTCGCGTTCCGCACATACATTTTGTGGGCATTGGAGGTGCCGGTATGGGCGGAATTGCCGAAGTTCTGGCCAATGAAGGATATCAGATCAGAGGTTCCGATTTAGCGACAAATCCGGTAACGCAGCAGTTAACGAATCTGGGTGCGACGATTTATTTCACCCATCGCCCGGAAAACGTACGTGATGCCAGCGGGATCGTTGttt                                                                                                                                                                            >  1:202271/1‑251 (MQ=255)                                                                      atACACAACAATTGGCAACACTGCGTTCCATCGTGCCCGAAATGAGTCGCGTTCGGCACATACATTTTGTCGGCATTGGTGGTGCCGGTATGGGCGGTATTGCCGAAGTTCTGGCCAATGAAGGTTATCAGATCAGTGGTTCCGATTTAGCGCCAAATCCGGTCACGCAGCAGTTAACGAATCTGGGTGCGACGATTTATTTCAACCATCGCCCGGAAAACGTACGGGATGCCAGCGTGGCCGTTGTTTcc                                                                                                                                                                          >  2:228365/1‑251 (MQ=255)                                                                              caaTTGGCAAAACTGCGGTCCATCGTGCCCGAAATGCGTCGCGTTCGGCACATACATTTTGTCGGCATTGGTGGTGCCGGTATGGGCGGTATTGCCGAAGTTCTGGCCAATGAAGGTTATCAGATCAGTGGTTCCGATTTAGCGCCAAATCCGGTCACGCAGCAGTTAACGAATCTGGGTGCGACGATTTATTTCAACCATCGCCCGGAAAACGTACGTGATGCCAGCGTGGTCGTTGTTTCCAGCGCGAt                                                                                                                                                                  <  1:7425/251‑1 (MQ=255)                                                                               aaTTGGCAAAACTGCGTTCCAGCGGGCCCGAAATGCGTCGCGTTCGGCACATACATTTTGTCGGCATTGGTGGTGCCGGTATGGGCGGTATTGCCGAAGTTCTGGCCAATGAAGGTTATCAGATCAGTGGTTCCGATTTAGCGCCAAATCCGGTCACGCAGCAGTTAACGAATCTGGGTGCGACGATTTATTTCAACCATCGCCCGGAAAACGTACGTGATGCCAGCGTGGTCGTTGTTTCCAGCGCGAtt                                                                                                                                                                 <  2:79883/251‑1 (MQ=255)                                                                                aTTGGCAAAACTGCGTTCCATCGTGCCCGAAATGCGTCGCGTTCGGCACATACATTTTGTCGGCATTGGTGGTGCCGGTATGGGCGGTATTGCCGAAGTTCTGGCCAATGAAGGTTATCAGATCAGTGGTTCCGATTTAGCGCCAAATCCGGTCACGCAGCAGTTAACGAATCTGGGTGCGACGATTTATTTCAACCATCGCCCGGAAAACGTACGTGATGCCAGCGTGGTCGTTGTTTCCAGCGCGAttt                                                                                                                                                                >  1:102831/1‑251 (MQ=255)                                                                                   ggCAAAACTGCGTTCCATCGTGCCCGAAATGCGTCGCGTTCGGCACATACATTTTGTCGGCATTGGTGGTGCCGGTATGGGCGGTATTGCCGAAGTTCTGGCCAATGAAGGTTATCAGATCAGTGGTTCCGATTTAGCGCCAAATCCGGTCACGCAGCAGTTAACGAATCTGGGTGCGACGATTTATTTCAACCATCGCCCGGAAAACGTACGTGATGCCAGCGTGGTCGTTGTTTCCAGCGCGATTTCTg                                                                                                                                                             <  1:300230/251‑1 (MQ=255)                                                                                           tGCGTTCCATCGGGACCGAAATGCGTGGCGTTCGGCACATACATTATGTCGGCATTGGTGGTGCCGGTATGGGCGGTATTGCCGCAGTTCTGGCCAATGAAGGTTATCAGATCAGTGGTTCCGATTTAGCGCCAAATCCGGACACGCAGAAGTTAACGAATATGGGTGCGACGATTTATTTCCACCATCGCCCGGAAAACGTACGTGTTGACCGCGCGGTCGTTGTTTCCAGCTCGATTACTGCCGATAAc                                                                                                                                                     >  1:8683/1‑251 (MQ=255)                                                                                            gCGTTCCATCGTGCCCGAAATGCGTCGCGTTCGGCACATACATTTTGTCGGCATTGGTGGTGCCGGTATGGGCGGTATTGCCGAAGTTCTGGCCAATGAAGGTTATCAGATCAGTGGTTCCGATTTAGCGCCAAATCCGGTCACGCAGCAGTTAACGAATCTGGGTGCGACGATTTATTTCAACCATCGCCCGGAAAACGTACGTGATGCCAGCGTGGTCGTTGTTTCCAGCGCGATTTCTGCCGATaa                                                                                                                                                      <  1:245769/249‑1 (MQ=255)                                                                                                tCCATCGAGCCCGAAATGCGTCGCGTTCGGCACATACATTTTGTCGGCATTGGTGGTGCCGGTATGGGCGGTATTGCCGAAGTTCTGGCCAATGAAGGTTATCAGATCAGTGGTTCCGATTTAGCGCCAAATCCGGTCACGCAGCAGTTAACGAATCTGGGTGCGACGATTTATTTCAACCATCGCCCGGAAAACGTACGTGATGCCAGCGTGGTCGTTGTTTCCAGCGCGATTTCTGCCGATAACCCgg                                                                                                                                                 >  2:172455/1‑250 (MQ=255)                                                                                                   aTCGTGCCCGAAATGCGTCGCGTTCGGCACATACATTTTGTCGGCATTGGTGGTGCCGGTATGGGCGGTATTGCCGAAGTTCTGGCCAATGAAGGTTATCAGATCAGTGGTTCCGATTTAGCGCCAAATCCGGTCACGCAGCAGTTAACGAATCTGGGTGCGACGATTTATTTCAACCATCGCCCGGAAAACGTACGTGATGCCAGCGTGGTCGTTGTTTCCAGCGCGATTTCTGCCGATAACCCGGAAAt                                                                                                                                             >  1:227396/1‑251 (MQ=255)                                                                                                     cGTGCCCGAAATGCGTCGCGTTCGGCACATACATTTTGTCGGCATTGGTGGTGCCGGTATGGGGGGTATTGCCGAAGTTCTGGCCAATGAAGGTTATCAGATCAGTGGTTCCGATTTAGCGCCAAATCCGGTCACGCAGCAGTTAACGCATCTGGGTGCGACGATTTATTTCAACCATCGCCCGTAAAACGTACGTGATGCCAGCGTGGTCGTCGGTTCCAGCGCGATTTCGGCCTATAACCCGGAAAtt                                                                                                                                            <  2:50495/250‑1 (MQ=255)                                                                                                     cGTGCCCGAAATGCGTCGCGTTCGGCACATACATTTTGTCGGCATTGGTGGTGCCGGTATGGGCGGTATTGCCGAAGTTCTGGCCAATGAAGGTTATCAGATCAGTGGTTCCGATTTAGCGCCAAATCCGGTCACGCAGCAGTTAACGAATCTGGGTGCGACGATTTATTTCAACCATCGCCCGGAAACCGTACGTGATGCCAGCGTGGTCGTTGTTTCCAGCGCGATTTCTGCCGATAACCCGGAAATTg                                                                                                                                           >  2:1664/1‑251 (MQ=255)                                                                                                           cGAAATGCGTCGCGTTCGGCACATACATTTTGTCGGCATTGGTGGTGCCGGTATGGGCGGTATTGCCGAAGTTCTGGCCAATGAAGGTTATCAGATCAGTGGTTCCGATTTAGCGCCAAATCCGGTCACGCAGCAGTTAACGAATCTGGGTGCGACGATTTATTTCAACCATCGCCCGGAAAACGTACGTGATGCCAGCGTGGTCGTTGTTTCCAGCGCGATTTCTGCCGATAACCCGGAAATTGTcgccg                                                                                                                                     <  1:294915/251‑1 (MQ=255)                                                                                                               aTGCGTCGCGTTCGGCACATACATTTTGTCGGCATTGGTGGTGCCGGTATGGGCGGAATTGCCGAAGTTCTGGCCAATGAAGGTTATCAGATCAGTGGTTCCGATTTAGCGCCAAATCCGGTCACGCAGCAGTTAACGAATCTGGGTGCGACGATTTATTTCAACCATCGCCCGGAAAACGTACGTGATGCCAGCGGGGTCGTTGTTTCCAGCGCGATTTCTGCCGATACCCCGGAAATTGTCGCCGCTc                                                                                                                                  >  2:78799/1‑250 (MQ=255)                                                                                                                          tCGGCACATACATTTTGTCGGCATTGGTGGTGCCGGTATGGGCGGTATTGCCGAAGTTCTGGCCAATGAAGGTTATCAGATCAGTGGTTCCGATTTAGCGCCAAATCCGGTCACGCAGCAGTTAACGAATCTGGGTGCGACGATTTATTTCAACCATCGCCCGGAAAACGTACGTGATGCCAGCGTGGTCGTTGTTTCCAGCGCGATTTCTGCCGATAACCCGGAAATTGTCGCCGCTCATGAAGCGCGt                                                                                                                       <  2:177766/250‑1 (MQ=255)                                                                                                                                   acatTTTGTCGGCATTGGTGGTGCCGGTATGGGCGGTATTTCCGAAGTTCTGGCCAAGGAAGGTTAGCAGATCAGTGGTTCCGATTTAGCGCCAAATCCGGTCACGCAGCAGTTAACGAATCGGGGTGCGACGATTTATTTCAACCATCCCCCGGAAAACGTACGTGATGCCAGCGTGTTCGTTGTTTCCAGCGCGACTTCTGCCGATACCCCGGAAATTGTCGCCGCTCATGAATCGCGTATTCCGGTGa                                                                                                             <  1:191446/251‑1 (MQ=255)                                                                                                                                              gCATTGGTGGTGCCGGTATGGGCGGTATTGCCGAAGTTCTGGCCAATGAAGGTTATCAGATCAGTTGTTCCGATTTAGCGCCAAATCCGGTCACGCAGCAGTTAACGAATCGGGGTGCGACGATTTATTTCAACCATCGCCCGGAAAACGTACGTGATGCCAGCGTGGGCGTTGTTTCCAGCGCGATTTCTGCCGATAACCCGGAAATTGTCGCCGCTCATGAAGCGCGTATTCCTGTGATCCGGCGTGcc                                                                                                  <  2:48961/251‑1 (MQ=255)                                                                                                                                                                   gCGGTATTGCCGAAGTTCTGGCCAATGAAGGTTATCAGATCAGTGGTTCCGATTTAGCGCCAAATCCGGTCACGCAGCAGTTAACGAATCTGGGTGCGACGATTTATTTCAACCATCGCCCGGAAAACGTACGTGATGCCAGCGTGGTCGTTGTTTCCAGCGCGATTTCTGCCGATAACCCGGAAATTGTCGCCGCTCATGAAGCGCGTATTCCGGTGATCCGTCGTGCCGAAATGCTGGCTGAGTTAATg                                                                             <  1:262431/251‑1 (MQ=255)                                                                                                                                                                          tGCCGAAGTTCTGGCCAATGAAGGTTATCAGATCAGTGGTTCCGATTTAGCGCCAAATCCGGTCACGCAGCAGTTAACGAATCTGGGTGCGACGATTTATTTCAACCATCGCCCGGAAAACGTACGTGATGCCAGCGTGGTCGTTGTTTCCAGCGCGATTTCTGCCGATAACCCGGAAATTGTCGCCGCTCATGAAGCGCGTATTCCGGTGATCCGTCGTGCCGAAATGCTGGCTGAGTTAATGCGtttt                                                                       <  1:124494/250‑1 (MQ=255)                                                                                                                                                                                aGTTCTGGCCAATGAAGGTTATCAGATCAGTGGTTCCGATTTAGCGCCAAATCCGGTCACGCAGCAGTTAACGAATCTGGGTGCGACGATTTATTTCAACCATCGCCCGGAAAACGTACGTGATGCCAGCGTGGTCGTTGTTTCCAGCGCGATTTCTGCCGATAACCCGGAAATTGTCGCCGCTCATGAAGCGCGTATTCCGGTGATCCGTCGTGCCGAAATGCTGGCTGAGTTAATGCGTTTTCGTCATg                                                                >  2:123286/1‑251 (MQ=255)                                                                                                                                                                                   tCTGGCCAAGGAAGGTTATCAGATCAGTGGTTCCGATTTAGCGCCAAATCCGGTCACGCAGCAGTTAACGAATCTGGGTGCGACGATTTATTTCAACCATCGCCCGGAAAACGTACGTGATGCCAGCGTGGTCGTTGTTTCCAGCGCGATTTCTGCCGATAACCCGGAAATTGTCGCCGCTCATGAAGCGCGTATTCCGGTGATCCGTCGTGCCGAAATGCTGGCTGAGTTAATGCGTTGTCGTCATGGc                                                              <  1:78799/250‑1 (MQ=255)                                                                                                                                                                                    cTGGCCAATGAAGGTTATCAGATCAGTGGTTCCGATTTAGCGCCAAATCCGGTCACGCAGCAGTTAACGAATCTGGGTGCGACGATTTATTTCAACCATCGCCCGGAAAACGTACGTGATGCCAGCGTGGTCGTTGTTTCCAGCGCGATTTCTGCCGATAACCCGGAAATTGTCGCCGCTCATGAAGCGCGTATTCCGGTGATCCGTCGTGCCGAAATGCTGGCTGAGTTAATGCGTTTTCGTCATGGCAt                                                            >  1:101498/1‑251 (MQ=255)                                                                                                                                                                                        ccAATGAAGGTTATCAGATCAGTGGTTCCGATTTAGCGCCAAATCCGGTCACGCAGCAGTTAACGAATCTGGGTGCGACGATTTATTTCAACCATCGCCCGGAAAACGTACGTGATGCCAGCGTGGTCGTTGTTTCCAGCGCGATTTCTGCCGATAACCCGGAAATTGTCGCCGCTCATGAAGCGCGTATTCCGGTGATCCGTCGTGCCGAAATGCTGGCTGAGTTAATGCGTTTTCGTCATGGCATCGcc                                                        <  1:1664/251‑1 (MQ=255)                                                                                                                                                                                                              tGGTTCAGATTTAGCGCCAAATCCGGTCACGCAGCAGTTAACGAATCTGGGTGCGACGCTTTATTTCAACCATCGCCCGGAAAACGTACGTGATGCCAGCGTGGTCGTTGTTTCCAGCGCGATTTCTGCCGATAACCCGGAAATTGTCGCCGCTCATGAAGCGCGTATTCCGGTGATCCGTCGTGCCGAAATGCTGGCTGAGTTAATGCGTTTTCGTCATGGCATCGCCATTGCCGGAACGCACGGCaaa                                   >  2:210941/1‑250 (MQ=255)                                                                                                                                                                                                                      aTTTAGCGCCAAATCCGGTCACGCAGCAGTTAACGAATCTGGGTGCGACGATTTATTTCAACCATCGCCCGGAAAACGTACGTGATGCCAGCGTGGTCGTTGTTTCCAGCGCGATTTCTGCCGATAACCCGGAAATTGTCGCCGCTCATGAAGCGCGTATTCCGGTGATCCGTCGTGCCGAAATGCTGGCTGAGTTAATGCGTTTTCGTCATGGCATCGCCATTGCCGGAACGCACGGCAAAACGACAacc                          >  1:306055/1‑251 (MQ=255)                                                                                                                                                                                                                          aGCGCCAAATCCGGTCACGCAGCAGTTAACGAATCTGGGTGCGACGATTTATTTCAACCATCGCCCGGAAAACGTACGTGATGCCAGCGTGGTCGTTGTTTCCAGCGCGATTTCTGCCGATAACCCGGAAATTGTCGCCGCTCATGAAGCGCGTATTCCGGTGATCCGTCGTGCCGAAATGCTGGCTGAGTTAATGCGTTTTCGTCATGGCATCGCCATTGCCGGAACGCACGGCAAAACGACAACCACcg                      <  1:210941/251‑1 (MQ=255)                                                                                                                                                                                                                               cAAATCCGGTCACGCAGCAGTTAACGAATCTGGGTGCGACGATTTATTTCAACCATCGCCCGGAAAACGAACGTGCTGCCAGCGTGGTCGTTGTTTCCAGCGCGATTTCTGCCGATAACCCGGAAATTGTCGCCGCTCATGAAGCGCGTATTCCGGTGATCCGTCATGCCGAAATGCTGGCTGAGTTAATGCGTTTTCGTCATGGCATCGCCATTGCCGGAACGCACGGCCAAACGCCCACCACcgcgctg                 >  2:221012/1‑248 (MQ=255)                                                                                                                                                                                                                                  aTCCGGTCACGCAGCAGTTAACGAATCTGGGTGCGACGATTTATTTCAACCATCGCCCGGAAAACGTACTTGATGCCAGCGTGGTCGTTGTTTCCAGCGCGATTTCTGCCGATAACCCGGAAATTGTCGCCGCTCATGAAGCGCGTATTCCGGTGATCCGTCGTGCCGAAATGCTGGCTGAGTTAATGCGTTTTCGTCATGGCATCGCCATTGCCGGAACGCACGGCAAAACGACAACCACCGCGATGGtt              <  2:306055/251‑1 (MQ=255)                                                                                                                                                                                                                                  aTCCGGTCACGCAGCAGTTAACGAATCTGGGTGCGACGATTTATTTCAACCATCGCCCGGAAAACGTACGTGATGCCAGCGTGGTCGTTGTTTCCAGCGCGATTTCTGCCGATAACCCGGAAATTGTCGCCGCTCATGAAGCGCGTATTCCGGTGATCCGTCGTGCCGAAATGCTGGCTGAGTTAATGCGTTTTCGTCATGGCATCGCCATTGCCGGAACGCACGGCAAAACGACAACCACCGCGATGGtt              <  1:229061/251‑1 (MQ=255)                                                                                                                                                                                                                                  aTCCGGTCACGCAGCAGTTAACGAATCTGGGTGCGACGATTTATTTCAACCATCGCCCGGAAAACGTACGTGATGCCAGCGTGGTCGTTGTTTCCAGCGCGATTTCTGCCGATAACCCGGAAATTGTCGCCGCTCATGAAGCGCGTATTCCGGTGATCCGTCGTGCCGAAATGCTGGCTGAGTTAATGCGTTTTCGTCATGGCATCGCCATTGCCGGAACGCACGGCAAAACGACAACCACCGCGATGGtt              >  1:1683/1‑251 (MQ=255)                                                                                                                                                                                                                                       gTCACGCAGCAGTTAACGAATCTGGGTGCGACGATTTATTTCAACCATCGCCCGGAAAACGTACGTGATGCCAGCGTGGTCGTTGTTTCCAGCGCGATTTCTGCCGATAACCCGGAAATTGTCGCCGCTCATGAAGCGCGTATTCCGGTGATCCGTCGTGCCGAAATGCTGGCTGAGTTAATGCGTTTTCGTCATGGCATCGCCATTGCCGGAACGCACGGCAAAACGACAACCACCGCGATGGTTTcc           >  2:70073/1‑249 (MQ=255)                                                                                                                                                                                                                                              agcagTTAACGAATCTGGGTGCGACGATTTATTTCAACCATCGCCCGGAAAACGTACGTGATGCCAGCGTGGTCGTTGTTTCCAGCGCGATTTCTGCCGATAACCCGGAAATTGTCGCCGCTCATGAAGCGCGTATTCCGGTGATCCGTCGTGCCGAAATGCTGGCTGAGTTAATGCGTTTTCGTCATGGCATCGCCATTGCCGGAACGCACGGCAAAACGACAACCACCGCGATGGTTTCCAGCATCTAc  <  1:70073/251‑1 (MQ=255)                                                                                                                                                                                                                                                       |                                                                                                                                                                                                                                                  GCCCGGGCGTAATTGTAGCGATGCCTTTTGCATCGTATGAATTTAAGAAGTTAATGGCGTAAAGAATGAATACACAACAATTGGCAAAACTGCGTTCCATCGTGCCCGAAATGCGTCGCGTTCGGCACATACATTTTGTCGGCATTGGTGGTGCCGGTATGGGCGGTATTGCCGAAGTTCTGGCCAATGAAGGTTATCAGATCAGTGGTTCCGATTTAGCGCCAAATCCGGTCACGCAGCAGTTAATGAATCTGGGTGCGACGATTTATTTCAACCATCGCCCGGAAAACGTACGTGATGCCAGCGTGGTCGTTGTTTCCAGCGCGATTTCTGCCGATAACCCGGAAATTGTCGCCGCTCATGAAGCGCGTATTCCGGTGATCCGTCGTGCCGAAATGCTGGCTGAGTTAATGCGTTTTCGTCATGGCATCGCCATTGCCGGAACGCACGGCAAAACGACAACCACCGCGATGGTTTCCAGCATCTAC  >  NC_000913/100700‑101187

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