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%	249.7 / NA	72	M61T (ATG→ACG)	murC	UDP‑N‑acetylmuramate‑‑L‑alanine ligase
Reads supporting (aligned to +/- strand):  ref base T (0/0);  new base C (40/32);  total (40/32)

GCCCGGGCGTAATTGTAGCGATGCCTTTTGCATCGTATGAATTTAAGAAGTTAATGGCGTAAAGAATGAATACACAACAATTGGCAAAACTGCGTTCCATCGTGCCCGAAATGCGTCGCGTTCGGCACATACATTTTGTCGGCATTGGTGGTGCCGGTATGGGCGGTATTGCCG‑AAGTTCTGGCCAATGAAGGTTATCAGATCAGTGGTTCCGATTTAGCGCCAAATCCGGTCACGCAGCAGTTAATGAATCTGGGTGCGACGATTTATTTCAACCATCGCCCGGAAAACGTACGTGATGCCAGCGTGGTCGTTGTTTCCAGCGCGATTTCTGCCGATAACCCGGAAATTGTCGCCGCTCATGAAGCGCGTATTCCGGTGATCCGTCGTGCCGAAATGCTGGCTGAGTTAATGCGTTTTCGTCATGGCATCGCCATTGCCGGAACGCACGGCAAAACGACAACCACCGCGATGGTTTCCAGCATCTACG  >  NC_000913/100700‑101188                                                                                                                                                                                                                                                        |                                                                                                                                                                                                                                                   gCCCGGGCGTAATTGTAGCGATGCCTTTTGCATCGTATGAATTTAAGAAGTTAATGGCGTAAAGAATGAATACACAACAATTGGCAAAACTGCGTTCCATCGTGCCCGAAATGCGTCGCGTTCGGCACATACATTTTGTCGGCATTGGTGGTGCCGGTATGGTCGGTATTGCCG‑AAGTTCTGGCCAATGAAGGTTATCAGATCAGTGGTTCCGATTTAGCGCCAAATCCGGTCACTCCGCAGTTAACGaa                                                                                                                                                                                                                                                 <  2:424475/250‑1 (MQ=255)  cccGGGCGTAATTGTAGCGATGCCTTTTGCATCGTATGAATTTAAGAAGTTAATGGCGTAAAGAATGAATACACAACAATTGGCAAAACTGCGTTCCATCGTGCCCGAAATGCGTCGCGTTCGGCACATACATTTTGTCGGCATTGGTGGTGCCGGTATGGGCGGTATTGCCG‑AAGTTCTGGCCAATGAAGGTTATCAGATCAGTGGTTCCGATTTAGCGCCAAATCCGGTCACGCAGCAGTTAACGAATc                                                                                                                                                                                                                                               >  2:164562/1‑251 (MQ=255)  cccGGGCGTAATTGTAGCGATGCCTTTTGCATCGTATGAATTTAAGAAGTTAATGGCGTAAAGAATGAATACACAACAATTGGCAAAACTGCGTTCCATCGTGCCCGAAATGCGTCGCGTTCGGCACATACATTTTGTCGGCATTGGTGGTGCCGGTATGGGCGGTATTGCCG‑AAGTTCTGGCCAATGAAGGTTATCAGATCAGTGGTTCCGATTTAGCGCCAAATCCGGTCACGCAGCAGTTAACGAATc                                                                                                                                                                                                                                               >  1:263387/1‑251 (MQ=255)            aTTGTAGCGATGCCTTTTGCATCGTATGAATTTAAGAAGTTAATGGCGTAAAGAATGAATACACAACAATTGGCAAAACTGCGTTCCATCGTGCCCGAAATGCGTCGCGTTCGGCACATACATTTTGTCGGCATTGGTGGTGCCGGTATGGGCGGTATTGCCG‑AAGTTCTGGCCAATGAAGGTTATCAGATCAGTGGTTCCGATTTAGCGCCAAATCCGGTCACGCAGCAGTTAACGAATCTGGGTGcgac                                                                                                                                                                                                                                     >  1:238850/1‑251 (MQ=255)              tGTAGCGATGCCTTTTGCATCGTATGAATTTAAGAAGTTAATGGCGTAAAGAATGAATACACAACAATTGGCAAAACTGCGTTCCATCGTGCCCGAAATGCGTCGCGTTCGGCACATACATTTTGTCGGCATTGGTGGTGCCGGTATGGGCGGTATTGCCG‑AAGTTCTGGCCAATGAAGGTTATCAGATCAGTGGTTCCGATTTAGCGCCAAATCCGGTCACGCAGCAGTTAACGAATCTGGGTGcgacg                                                                                                                                                                                                                                    >  1:449072/1‑250 (MQ=255)               gTAGCGATGCCTTTTGCATCGTATGAATTTAAGAAGTTAATGGCGTAAAGAATGAATACACAACAATTGGCAAAACTGCGTTCCATCGTGCCCGAAATGCGTCGCGTTCGGCACATACATTTTGTCGGCATTGGTGGTGCCGGTATGGGCGGTATTGCCG‑AAGTTCTGGCCAATGAAGGTTATCAGATCAGTGGTTCCGATTTAGCGCCAAATCCGGTCACGCAGCAGTTAACGAATCTGGGTGCGACGat                                                                                                                                                                                                                                  <  1:141421/251‑1 (MQ=255)                         cTTTTGCATCGTATGAATTTAAGAAGTTAATGGCGTAAAGAATGAATACACAACAATTGGAAAAACTGCGTTCCATCGTGCCCGCAATGCGTCGCGTTCGGCACATACATTTTGTCGGCATTGGTGGTGCCGGTATGGGCGGTATTGCCG‑AAGTTCTGGCCCATGAAGGTTCTCAGATCAGTGGTTCCGATTTAGCGCCAAATCCGGTAAAGCAGCAGTTAACGAATATGGGTACGACGATTTATTTCAAc                                                                                                                                                                                                                        >  2:430759/1‑251 (MQ=255)                                       gAATTTAAGAAGTTAATGGCGTAAAGAATGAATACACAACAATTGGCAAAACTGCGTTCCATCGTGCCCGAAATGCGTCGCGTTCGGCACATACATTTTGTCGGCATTGGTGGTGCCGGTATGGGCGGTATTGCCG‑AAGTTCTGGCCAATGAAGGATATCAGATCAGTGGTTCCGATTTAGCGCCAAATCCGGTCACGCAGCAGTTAACGAATCTGGGTGCGACGATTTATTTCAACCATCGCCCGGaaa                                                                                                                                                                                                           >  1:205177/1‑250 (MQ=255)                                       gAATTTAAGAAGTTAATGGCGTAAAGAATGAATACACAACAATTGGCAAAACTGCGTTCCATCGTGACCGAAATGCGTCGCGTTCGGCACATACATTTTGTCGGCATTGGTGGTGACGGTCTGGGCGGTATTGCCG‑AAGTTCTGGCCAATGAAGGTTATCAGATCAGTGGTTCCGATTTAGCGCCAAATCCGGTCACGCAGCAGTTAACGAATCTGGGTGCGACGATTTATTTCAACCATCGCCCggcaa                                                                                                                                                                                                           >  2:439668/1‑247 (MQ=255)                                             aagaagTTAATGGCGTAAAGAATGAATACACAACAATTGGCAAAACTGCGTTCCATCGTGCCCGAAATGCGTCGCGTTCGGCACATACATTTTGTCGGCATTGGTGGTGCCGGTATGGGCGGTATTGCCG‑AAGTTCTGGCCAATGAAGGTTATCAGATCAGTGGTTCCGATTTAGCGCCAAATCCGGTCACGCAGCAGTTAACGAATCTGGGTGCGACGATTTATTTCAACCATCGCCCGGAAAacgtacg                                                                                                                                                                                                    >  2:111160/1‑251 (MQ=255)                                             aagaagGTATTGGCGGAAAGAAAGAATACACAACAATTGGAAAAACTGCGTTCCATCGAGCCCGAAATGCGTAGCGTTCGGCACATAAATTTTGTCGGCATTGATGGAGCCGGTATGGGCGGTATTGCCG‑AAGTTCTGGCCAATGAAGGTTATCAGATCAGTGGTTCCGATTTCGCGCCAAATCGGGTCACGCAGAAGATAACGAATCTGGGTGCGATGATTTATTTAAACCATCGCCAGGaaaaagaacg                                                                                                                                                                                                    >  2:417821/1‑245 (MQ=255)                                                aagTTAATGGCGTAAAGAATGAATACACAACAATTGGCAAAACTGCGTTCCATCGTGCCCGAAATGCGTCGCGTTCGGCACATACATTTTGTCGGCATTGGTGGTGCCGGTATGGGCGGTATTGCCG‑AAGTTCTGGCCAATGAAGGTTATCAGATCAGTGGTTCCGATTTAGCGCCAAATCCGGTCACGCAGCAGTTAACGAATCTGGGGGCGACGATTTATTTCAACCATCTCCCGGCAACCGTACGTg                                                                                                                                                                                                  >  1:441156/1‑250 (MQ=255)                                                 agTTAATGGCGTAAAGAATGAATACACAACAATTGGCAAAACTGCGTTCCATCGTGCCCGAAATGCGTCGCGTTCGGCACAGACATTTTGTCGGCATTGGTGGTGCCGGTATGGGCGGTATTGCCG‑AAGTTCTGGCCAATGAAGGTTATCAGATCAGTGGTTCCGAGTTAGCGCCAAATCCGGTCACGCAGCAGTTAACGAATCTGGGTGCGACGATTTATTTCAACCATCGCGCGGAAAACGTACGTGAt                                                                                                                                                                                                <  2:221038/251‑1 (MQ=255)                                                    tAATGGCGTAAAGAATGAATACACAACAATTGGCAAAACTGCGTTCCCGCGTGCCCGAAATGCGTAGAGTTCGGCACATACATTTTGTCGGCATTGGTGGTGCCGGTATGGGCGGTATTGCCG‑AAGTTCTGGCAACTGAAGGTTATCAGATCGGGGGTTCCGATTTAGCGCCAAATCCGGTCACGCAGCAGTTACCGAATCTGGGGGCGACGATTTATTTCAACCATCGACCGGAGAACGTACGTGCTGcc                                                                                                                                                                                             >  2:501623/1‑251 (MQ=255)                                                    tAATGGCGTAAAGAATGAATACACAACAATTGGCAAAACTGCGTTCCATCGTGCCCGAAATGCGTCGCGTTCGGCACATACATTTTGTCGGCATTGGTGGGGCCGGTATGGGCGGTATTGCCG‑AAGTTCTGGCCAATGAAGGTTATCAGATCAGTGGTTCCGATTTAGCGCCAAATCCGGTCACGCAGCAGTTAACGAATCTGGGTGCGACGATTTATTTCAACCATCGCCCGGAAAACGTACGGGATGcc                                                                                                                                                                                             >  2:84556/1‑251 (MQ=255)                                                     aaTGGCGTAAAGAATGAATACACAACAATTGGCAAAACTGCGTTCCATCGTGCCCGAAATGCGTCGCGTTCGGCACATACATTTTGTCGGCATTGGTGGTGCCGGTATGGGCGGTATTGCCG‑AAGTTCTGGCCAATGAAGGTTATCAGATCAGTGGTTCCGATTTAGCGCCAAATCCGGTCACGCAGCAGTTAACGAATCTGGGTGCGACGATTTATTTCAACCATCGCCCGGAAAACGTACGTGATGCCa                                                                                                                                                                                            <  1:437293/251‑1 (MQ=255)                                                          cGTAAAGAATGAATACACAACAATTGGCAAAACTGCGTTCCATCGTGCCCGAAATGCGTCGCGTTCGGCACATACATTTTGTCGGCATTGGTGGTGCCGGTATGGGCGGTATTGCCG‑AAGTTCTGGCCAATGAAGGTTATCAGATCAGTGGTTCCGATTTAGCGCCAAATCCGGTCACGCAGCAGTTAACGAATCTGGGTGCGACGATTTATTTCAACCATCGCCCGGAAAACGTACGTGATGCCAGCGTg                                                                                                                                                                                       >  1:236888/1‑251 (MQ=255)                                                           gTAAAGAATGAATACACAACAATTGGCAAAACTGCGTTCCATCGTGCCCGAAATGCGTCGCGTTCGGCACATACATTTTGTCGGCATTGGTGGTGCCGGTATGGGCGGTATTGCCG‑ACGTTCTGGCCAATGAAGGTTATCAGATCAGTGGTTCCGATTTAGCGCCAAATCCGGTCACGCAGCAGTTAACGAATCTGGGTGCGACGATTTATTTCAACCATCGCCCGGAAAACGTACGTGATGCCAGCGTgg                                                                                                                                                                                      >  2:446486/1‑251 (MQ=255)                                                             aaaGAATGAATACACAACAATTGGCAAAACTGCGTTCCATCGTGCCCGAAATGCGTCGCGTTCGGCACATACATTTTGTCGGCATTGGTGGTGCCGGTATGGGCGGTATTGCCG‑AAGTTCTGGCCAATGAAGGTTATCAGATCAGTGGTTCCGATTTAGCGCCAAATCCGGTCACGCAGCAGTTAACGAGTCTGGGTGCGACGATTTATTTCACCCATCGCCCGGAAAACGTACGTGATGCCATCGTGGt                                                                                                                                                                                     <  2:54235/250‑1 (MQ=255)                                                               aGAATGAATACACAACAATTGGCAAAACTGCGTTCCATCGTGCCCGAAATGCGTCGCGTTCGGCACATACATTTTGTCGGCATTGGTGGTGCCGGTATGGGCGGTATTGCCG‑AAGTTCTGGCCAATGAAGGTTATCAGATCAGTGGTTCCGATTTAGCGCCAAATCCGGGCACGCAGCAGTTAACGAATCTGGGTGCGACGATTTATTTCAACCATCGCCCGGAAAACGTACGTGATGCCAGCGTGGTCgt                                                                                                                                                                                  <  1:335030/251‑1 (MQ=255)                                                                   tgaatACACAACAATTGGCAAAACTGCGTTCCATCGTGCCCGAAATGCGTCGCGTTCGGCACATACATTTTGTCGGCATTGGTGGTGCCGGTATGGGCGGTATTGCCG‑AAGTTCTGGCCAATGAAGGTTATCAGATCAGTGGTTCCGATTTAGCGCCACATCCGGTCACGCAGCAGTTAACGAATCTGGGTGCGACGATTTATTTCAACCATCGCCCGGACAACGTACGTGATGCCAGCGTGGTCgttgtt                                                                                                                                                                              >  1:439727/1‑251 (MQ=255)                                                                        acacTACAATTGGCAAAACTGCGTTCCATCGTGCCCGAAATGCGACGCATTCGGCACATACATTTTGTCGGCATTGGTGGTGCCGGTATGGGCGGTATTGCCG‑AAGTTCTGGCCAATGAAGGTTATCAGATCAGTGGTTCCGCTTTAGCGCCAAATCCGGTCACGCAGCAGTTAACGAATCTGGGTGCGACGCTTACTTTCAACCATCGCCCGGAAAACGTACGTGAAGCCAGCGTGGGCgttgttgcc                                                                                                                                                                           >  2:75831/1‑246 (MQ=255)                                                                          acaacaATTGGCAAAACTGCGTTCCATCGTGCCCGAAATGCGTCGCGTTCGGCACATACATTTTGTCGGCATTGGTGGTGCCGGTATGGGCGGTATTGCCG‑AAGTTCTGGCCAATGAAGGTTATCAGATCAGTGGTTCCGATTTAGCGCCAAATCCGGTCACGCAGCAGTTAACGAATCTGGGTGCGACGATTTATTTCAACCATCGCCCGGAAAACGTACGTGATGCCAGCGTGGTCGTTGTTTCCAgcg                                                                                                                                                                       <  1:446486/251‑1 (MQ=255)                                                                            aacaaTTGGCAAAACTGCGTTCCATCGTGCCCGAAATGCGTCGCGTTCGGCACATACATTTTGTCGGCATTGGTGGTGCCGGTATGGGCGGTATTGCCG‑AAGTTCTGGCCAATGAAGGTTATCAGATCAGTGGTTCCGATTTAGCGCCAAATCCGGTCACGCAGCAGTTAACGAATCTGGGTGCGACGATTTATTTCAACCATCGCCCGGAAAACGTACGTGATGCCAGCGTGGTCGTTGTTTCCAGcgcg                                                                                                                                                                     >  1:138323/1‑251 (MQ=255)                                                                             acaaTTGGCAAAACTGCGTTCCATCGTGCCCGAAATGCGTCGCGTTCGGCACATACATTTTGTCGGCATTGGTGGTGCCGGTATGGGCGGTATTGCCG‑AAGTTCTGGCCAATGAAGGTTATCAGATCAGTGGTTCCGATTTAGCGCCAAATCCGGTCACGCAGCAGTTAACGAATCTGGGTGCGACGATTTATTTCAACCATCGCCCGGAAAACGTACGTGATGCCAGCGTGGTCGTTGTTTCCAGcgcg                                                                                                                                                                     >  1:281715/1‑250 (MQ=255)                                                                               aaTTGGCAAAACTGCGTTCCATCGTGCCCGAAATGCGTCGCGTTCGGCACATACATTTTGTCGGCATTGGTGGTGCCGGTATGGGCGGTATTGCCG‑AAGTTCTGGCCAATGAAGGTTATCAGATCAGTGGTTCCGATTTAGCGCCAAATCCGGTCACGCAGCAGTTAACGAATCTGGGTGCGACGATTTATTTCAACCATCGCCCGGAAAACGTACGTGATGCCAGCGTGGTCGTTGTTTCCAGCGCGAtt                                                                                                                                                                  >  1:374471/1‑251 (MQ=255)                                                                                    gCAAAACTGCGTTCCCTCGTGCCCGAAATGCGTCGCGTTCGGCCCATACATTTTGTCGGCCTTGGTGGTGCCGGTATGGGCGGTATTGCCG‑ACGTTCTGGCCAATGAAGGTTCTCAGATCAGTGGTTCCGATTTAGCGCCAAATCCGGTCACGCAGCAGTTAACGAATCTGGGTGCGACGATTTATTTCAACCATCGCCCGGAAAACGTACGTGATGCCCGCGTGGTCGTTGTTTCCAGCGCGATTTCTGc                                                                                                                                                             >  2:498345/1‑251 (MQ=255)                                                                                      aaaaCTGCGTTCCATCGTGCCCGAAATGCGTCGCGTTCGGCACATACATTTTGTCGGCATTGGTGGTGCCGGAATGGGCGGTATTGCCA‑CAGTTCTGGCCCATGAAGGTTATCAGATCAGTGGTTCCGATTTAGCGCCAAATCCGGTCACGCAGCAGTTAACGAATCTGGGGGCGAAGATTTATTTCAACCATCGCCCGGAAAACGTACGTGATGCCAGCGTGGTCGTTGTTTCCAGCGCGCTTTCTGCCg                                                                                                                                                           >  1:252483/1‑251 (MQ=255)                                                                                        aaCTGCGTTCCATCGTGCCCGAAATGCGTCGCGTTCGGCACATACATTTTGTCGGCATTGGTGGTGCCGTTATGGGCGGTATTGCCG‑AAGTTCTGGCCAATGAAGGTTATCAGATCAGTGGTTCCGATTTAGCGCCAAATCCGGTCACGCAGCAGTTAACGAATCTGGGTGcg                                                                                                                                                                                                                                       <  1:252340/173‑1 (MQ=255)                                                                                        aaCTGCGTTACATCGTGCCCGAAATGCGTCGCGTTCGGCACATACATTTTGTCGGCATTGGTGGTGCCGGTCTGGGCGGTATTGACG‑AAGTTCTGGCCAATGAAGGTTATCAGATCAGTGGTTCCGATTTAGCGCCAAATCCGGTCACGCAGCAGTTAACGAATCTGGGTGcg                                                                                                                                                                                                                                       >  2:252340/1‑173 (MQ=255)                                                                                         aCTGCGTTCCATCGTGCCCGAAATGCGTCGCGTTCGGCACATACATTTTGTCGGCATTGGTGGTGCCGGTATGGGCGGTATTGCCG‑AAGTTCTGGCCAATGAAGGTTATCAGATCAGTGGTTCCGATTTAGCGCCAAATCCGGTCACGCAGCAGTTAACGAATCTGGGTGCGACGATTTATTTCAACCATCGCCCGGAAAACGTACGTGATGCCAGCGTGGTCGTTGTTTCCAGCGCGATTTCTGCCGAt                                                                                                                                                         >  2:138780/1‑250 (MQ=255)                                                                                           tGCGTTCCATCGTGCCCGAAATGCGTCGCGTTCGGCACATACATTTTGTCGGCATTGGTGGTGCCGGTATGGGCGGTATTGCCG‑AAGTTCTGGCCAATGAAGGTTATCAGATCAGTGGTTCCGATTTAGCGCCAAATCCGGTCACGCAGCAGTTAACGAATCTGGGTGCGACGATTTATTTCAACCATCGCCCGGAAAACGTACGTGATGCCAGCGTGGTCGTTGTTTCCAGCGCGATTTCt                                                                                                                                                               <  2:404943/242‑1 (MQ=255)                                                                                           tGCGTTACATCGTGCCCGAAATGCGTCGCGTTCGGCACATACCTTTTGTCGGCATTGGTGGTGCCGGTATGGGCGGTATTGCCG‑AAGTTCTGGCCAATGAAGGTTATCAGATCAGTGGTTCCGATTTAGCGCCAAATCCGGTCACGCAGCAGTTAACGAATCTGGGTGCGACGATTTATTTCAACCATCGCCCGGAAAACGTACGTGATGCCAGCGTGGTCGTTGTTTCCAGCGCGATTTCt                                                                                                                                                               >  1:404943/1‑242 (MQ=255)                                                                                              gTTCCATCGTGCCCGAAATGCGTCGCGTTCGGCACATACATTTTGTCGGCATTGGTGGTGCCGGTATGGGCGGTATTGCCG‑AAGTTCTGGCCAATGAAGGTTATCAGATCAGTGGTTCCGATTTAGCGCCAAATCCGGTCACGCAGCAGTTAACGAATCTGGGTGCGACGATTTATTTCAACCATCGCCCGGAAAACGTACGTGATGCCAGCGTGGTCGTTGTTTCCAGCGCGATTTCTGCCGATAACCCg                                                                                                                                                   >  1:104117/1‑251 (MQ=255)                                                                                                tCCATCGTGCCCGAAATGCGTCGCGTTCGGCACATACATTTTGTCGGCATTGGTGGTGCCGGTATGGGCGGTATTGCCG‑AAGTTCTGGCCAATGAAGGTTATCAGATCAGTGGTTCCGATTTAGCGCCAAATCCGGTCACGCAGCAGTTAACGAATCTGGGTGCGACGATTTATTTCAACCATCGCCCGGAAAACGTACGTGATGCCAGCGTGGTCGTTGTTTCCAGCGCGATTTCTGCCGATAACCCgg                                                                                                                                                  <  1:11043/250‑1 (MQ=255)                                                                                                  cATCGTGCCCGAAATGCGTCGCGTTCGGCACATACATTTTGTCGGCATTGGTGGTGCCGGTATGGGCGGTATTGCCG‑AAGTTCTGGTCAATGAAGGTTATCAGATCAGTGGTTCCGATTTAGCGCCAAAGCCGGTCACGCAGCAGTTAACGAATCTGGGTGCGACGATTTATTTCAACCATCGCCCGGAAAACGTACGGGATGTCAGCGTGGTCGTTGTTTCCAGCGCGATTTCTGCCGATAACCCGGaaa                                                                                                                                               <  1:257006/251‑1 (MQ=255)                                                                                                   aTCGTGCCCGAAATGCGTCGCGTTCGGCACATACATTGTGTCGGCATTGGTGGTGCCGGTATGGGCGGTATTGCCG‑AAGTTCTGGCCAATGAAGGTTATCAGATCAGTGGTTCCGATTTAGCGCCAAATCCGGTCACGCAGCAGTTAACGAATCTGGGTGCGACGATTTATTTCAACCATCGCCCGGAAAACGTACGTGATGCCAGCGTGGTCGTTGTTTCCAGCGCGATTTCTGCCGATAACCCGGAAAt                                                                                                                                              <  2:236888/251‑1 (MQ=255)                                                                                                   aTCGTGCCCGAAATACGTCGCGTACTGCACCTACATTTTGTCGGCATTGGTGGTGCCGGTATGGGCGGTATTGCCG‑AAGTTCTGGCCAATGAAGGTTATCAGATCAGTGGTTCCGATTTAGCGCCAAATCCGGTCACGCAGCAGTTGACGAATCTGGGTGCGACGATTTATTTCACCCATCGCCCGGAAAACGTACGTGAGGCCAGCGTGGTCGTTGTTTCCAGCGCGATTTCTGCCGATACCCCGGAAAt                                                                                                                                              >  2:82973/1‑251 (MQ=255)                                                                                                     cGTGCCCGAAATGCGTCGCGTTCGGCACATACATTTTGTCGGCATTGGTGGTGCCGGTATGGGCGGTATTGCCG‑AAGTTCTGGCCAATGAAGGTTATCAGATCAGTGGTTCCGATTTAGCGCCAAATCCGGTCACGCAGCAGTTAACGAATCTGGGTGCGACGATTTATTTCAACCATCGCCCGGAAAACGTACGTGATTCCAGCGTGGTCGTTGTTGCTAGCGTGAGTTCTGCCGATAACCCGGAAATTg                                                                                                                                            <  1:316117/251‑1 (MQ=255)                                                                                                        gCCCGAAATGCGTCGCGTTCGGCACATACATTTTGTCGGCATTGGTGGTGCCGGTATGGGCGGTATTGCCG‑AAGTTCTGGCCAATGAAGGTTATCAGATCAGTGGTTCCGATTTAGCGCCAAATCCGGTCACGCAGCAGTTAACGAATCTGGGTGCGACGATTTATTTCAACCATCGCCCGGAAAACGTACGTGATGCCAGCGTGGTCGTTGTTTCCAGCGCGATTTCTGCCGATAACCCGGAAATTGTc                                                                                                                                          <  1:403314/250‑1 (MQ=255)                                                                                                          ccGAAATGCGTCGCGTTCGGCACATACATTTTGTCGGCATTGGTGGTGCCGGTATGGGCGGTATTGCCG‑AAGTTCTGGCCAATGAAGGTTATCAGATCAGTGGTTCCGATTTAGCGCCAAATCCGGTCACGCAGCAGTTAACGAATCTGGGTGCGACGATTTATTTCAACCATCGCCCGGAAAACGTACGTGATGCCAGCGTGGTCGTGGTTTCCAGCGCGATTTCTGCCGAGAACCCGGAAATGGTcgcc                                                                                                                                       <  2:209782/251‑1 (MQ=255)                                                                                                                    tCGCGTTCGGCACATGCATTTTGTCGGCATATGTGGTGCCGGTGTGGGCGGTCTTGCCGCAAGT‑CTGGCCAATGAAGCTTATTAGATCAGTGGTTCCGATTTAGCGCCAAATCCTGACACTCCGCAGTGAACGAATCAGGGTGCGACGATTGATTTCAACCATCGCCCGGAAAACGTACGTGATGCCAGCGTGGTCGTGGTTTCCAGCGCGATTTCTGCCGATAACCCGGTAATTGTCGCCGCTCATGaa                                                                                                                              <  2:206113/250‑1 (MQ=255)                                                                                                                      gcgTTCGGCACATACATTTTGTCGGCATTGGTGGTGCCGGTATGGGCGGTATTGCCG‑AAGTTCTGGCCAATGAAGGTTATCAGATCAGTGGTTCCGATTTAGCGCCAAATCCGGTCACGCAGCAGTTAACGAATCTGGGTGCGACGATTTATTTCAACCATCGCCCGGAAAACGTACGTGATGCCAGCGTGGTCGTTTTTTCCAGCGCGATTTCTGCCGATAACCCGGAAATTGTCGCCGCTCATGAAgcg                                                                                                                           <  2:388946/251‑1 (MQ=255)                                                                                                                        gTTCGGCACATACATTTTGTCGGCATTGGTGGTGCCGGTATGGGCGGTATTGCCG‑AAGTTCTGGCCAATGAAGGTTATCAGATCAGTGGTTCCGATTTAGCGCCAGATCCGGTCACGCAGCAGTTAACGAATCTGGGTGCGACGATTTATTTCAACCATCGCCCGGAAAACGTACGTGATGCCAGCGTGGTCGTTGTTTCCAGCGCGATTTCTGCCGATAACCCGGAAATTGTCGCCGCTCATGAAGcgcg                                                                                                                         <  1:138780/251‑1 (MQ=255)                                                                                                                        gTTCGGCACATACATTTTGTCGGCATTGGTGGTGCCGGTATGGGCGGAATTGCCG‑AAGTTCTGGCCAATGAAGGATAACAGACCAGTGGTTCCGATTTAGCGACAAATCCGGTCACGCAGCAGATACCGAATCTGGGAGCGACGATTTATTTCCAACATCGCCCGGCAACCGTACGTGATGCCAGCGTGGTCGTTGTTTCCAGCGCGATTTATGCCGATAACCCGGAAATTGTCGCCGCTCATGAcgcgcg                                                                                                                         >  1:207495/1‑251 (MQ=255)                                                                                                                        gTTCGGCACATACATTTTGACGGCATTGGCGGTGCCGGTATGGGCGGTATTGCCG‑AAGTTATGGCCAATGACGGTTACCCAATCCGTGGTTACGATTTAGCGCCACATCCGGTCACGCAGCAGTTAACGAATCTGGGTGCGACGCTTTATTTACACCATCGCCCGGAAAACGTACGTGATGCCAGAGTGGTCGTTGTTTCCAGCGCGCTTTCTGCAGATACCCCGGAAATTGTCGCCGCTCATGCAGcgcg                                                                                                                         >  2:490400/1‑251 (MQ=255)                                                                                                                            ggCACATACATTTTGTCGGCATTGGTGGTGCCGGTATGGGCGGTATTGCCG‑AAGTTCTGGCCAATGAAGGTTATCAGATCAGTGGTTCCGATTTAGCGCCAAATCCGGTCACGCAGCAGTTAACGAATCTGGGTGCGACGATTTATTTCAACCATCGCCCGGAAAACGTACGTGATGCCAGCGTGGTCGTTGTTTCCAGCGCGATTTCTGCCGATAACCCGGAAATTGTCGCCGCTCATGAAGCGCGTAtt                                                                                                                     <  1:387317/251‑1 (MQ=255)                                                                                                                            ggCACATACATTTTGTCGGCATTGGTGGTGCCGGTATGGGCGGTATTGCCG‑AAGTTCTGGCCAATGAAGGTTATCAGATCAGTGGTTCCGATTTAGCGCCAAATCCGGTCACGCAGCAGTTAACGAATCTGGGTGCGACGATTTATTTCAACCATCGCCCGGAAAACGTACGTGATGCCAGCGTGGTCGTTGTTTCCAGCGCGATTTCTGCCGATAACCCGGAAATTGTCGCCGCTCATGAAGCGCGTAtt                                                                                                                     <  2:217144/251‑1 (MQ=255)                                                                                                                              cacaTACATTTTGTCGGCATTGGTAGTGCCGGTATGGGAGGTATTGCCG‑AAGTTCTGGCCAATGAAGGTTATCAGATCAGTGGTTCCGATTTAGCGCCAAATCCGGTCACGCAGCAGTTAACGACTCTGGGTGCGACGATTTATTTCAACCATCGCCCGGAAAACGTACGTGATGCCAGCGTGGTCGTTGTTTCCAGCGCGATTTCGGCCGATAACCGGGAAATTGTCGCCGCTCATGAAGCGCGTATTcc                                                                                                                   >  2:422946/1‑251 (MQ=255)                                                                                                                                 atacatTTTGTCGGCATTGGTGGTGCCGGTATGGGCGGTATTGCCG‑AAGTTCTGGCCAATGAAGGTTATCAGATCAGTGGTTCCGATTTAGCGCCAAATCCGGTCACGCAGCAGTTAACGAATCTGGGTGCGACGATTTATTTCAACCATCGCCCGGAAAACGTACGTGATGCCAGCGTGGTCGTTGTTTCCAGCGCGATTTCTGCCGATAACCCGGAAATTGTCGCCGCTCATGAAGCGCGTATTCCGGt                                                                                                                <  1:452057/251‑1 (MQ=255)                                                                                                                                         tGTCGGCATTGGTTGTGTCGGTATGGGCGGTATTGCCG‑AAGTTCTGGCCAATTAAGGTTATCAGATCAGTGGTTCCGATTTAGCGCCAAATCCGGTCACGCAGCAGTTAACGAATCTGGGTGCGACGATTTATTTCAACCATCGCCCGGAAAACGTACGTGATGCCAGCGTGGTCGTTGTTTCCAGCGCGATTTCTGCCGAGAACCCGGAAATTGTCGACGCTCATGAAGCGCGTATTCCGGTGATCcgt                                                                                                         <  2:398437/250‑1 (MQ=255)                                                                                                                                                   ggtggtGCCGGTATGGGCGGTATTGCCG‑AAGTTCTGGCCAATGAAGGTTATCAGATCAGTGGTTCCGATTTAGCGCCAAATCCGGTCACGCAGCAGTTAACGAATCTGGGTGCGACGATTTATTTCAACCATCGCCCGGAAAACGTACGTGATGCCAGCGTGGTCGTTGTTTCCAGCGCGATTTCTGCCGATAACCCGGAAATTGTCGCCGCTCATGAAGCGCGTATTCCGGTGATTCGTCGTGCCGAAAt                                                                                              <  2:238850/251‑1 (MQ=255)                                                                                                                                                   ggtggtGCCGGTATGGGCGGTATTGCCG‑AAGTTCTGGCCAATGAAGGTTATCAGATCAGTGGTTCCGATTTAGCGCCAAATCCGGTCACGCAGCAGTTAACGAATCTGGGTGCGACGATTTATTTCAACCATCGCCCGGAAAACGTACGTGATGCCAGCGTGGTCGTTGTTTCCAGCGCGATTTCTGCCGATAACCCGGAAATTGTCGCCGCTCATGAAGCGCGTATTCCGGTGATCCGTCGTGCCGaaag                                                                                              <  2:263387/251‑2 (MQ=255)                                                                                                                                                    gtggtgCCGGTATGGGCGGTATTGCCG‑AAGTTCTGGCCAATGAAGGTTATCAGATCAGTGGTTCCGATTTAGCGCCAAATCCGGTCACGCAGCAGTTAACGAATCTGGGTGCGACGATTTATTTCAACCATCGCCCGGAAAACGTACGTGATGCCAGCGTGGTCGTTGTTTCCAGCGCGATTTCTGCCGATAACCCGGAAATTGTCGCCGCTCATGAAGCGCGTATTCCGGTGATCCGTCGTGCCGAAATg                                                                                             <  1:164562/251‑1 (MQ=255)                                                                                                                                                            ggTATGGGCGGTATTGCCG‑AAGTTCTGGCCAATGAAGGTTATCAGATCAGTGGTTCCGATTTAGCGCCAAATCCGGTCACGCAGCAGTTACCGAATCTGGGTGCGACGATTTATTTCAACCATCGCCCGGAAAACGTACGTGATGCCAGCGTGGTCGTTGTTTCCAGCGCGATTTCTGCCGATAACCCGGAAATTGTCGCCGCTCATGAAGCGCGTATTCCGGTGATCCGTCGTGCCGAAATgctggctg                                                                                      >  1:212252/1‑250 (MQ=255)                                                                                                                                                               aTGGGCGGTATTGCCG‑AAGTTCTGGCCAATGAAGGTTATCAGATCAGTGGTTCCGATTTAGCGCCAAATCCGGTCACGCAGCAGTTAACGAATCTGGGTGCGACGATTTATTTCAACCATCGCCCGGAAAACGTACGTGATGCCAGCGTGGTCGTTGTTTCCAGCGCGATTTCTGCCGATAACCCGGAAATTGTCGCCGCTCATGAAGCGCGTATTCCGGTGATCCGTCGTGCCGAAATGCTGGCTGAGtt                                                                                  <  1:111160/251‑1 (MQ=255)                                                                                                                                                                   gCGGTATTGCCG‑AAGTTCTGGCCAATGAAGGTTATCAGATCAGTGGTTCCGATTTAGCGCCAAATCCGGTCACGCAGCAGTTAACGAATCTGGGTGCGACGATTTATTTCAACCATCGCCCGGAAAACGTACGTGATGCCAGCGTGGTCGTTGTTTCCAGCGCGATTTCTGCCGATAACCCGGAAATTGTCGCCGCTCATGAAGCGCGTATTCCGGTGATCCGTCGTGCCGAAATGCTGGCTGAGTTAATg                                                                              >  1:84342/1‑251 (MQ=255)                                                                                                                                                                      gTATTGCCG‑AAGTTCTGGCCAATGAAGGTTATCAGTTCAGTGGTTCCGATTTAGCGCCAAATCCGGTCACGCAGCAGTTAACGAATCTGGGGGCGACGATTTATTTCAACCATCGCCCGGAAAACGTACGTGATGCCAGCGTGGTCGTTGTTTCCAGCGCGATTTCGGCCGATAACCCGGAAATTGTCGCCGCTCATGAAGCGCGTATTCCGGTGATCCGTCGTGCCGAAATGCTGGCTGAGTTAATGCGt                                                                           <  2:414259/251‑1 (MQ=255)                                                                                                                                                                                  gTTCTGGCCAATGAAGGTTATCAGATCAGTGGTTCCGATTTAGCGCCAAATCCGGTCACGCAGCAGTTAACGAATCTGGGTGCGACGATTTATTTCAACCATCGCCCGGAAAACGTACGTGATGCCAGCGTGATCGTTGTTTCAAGCGCGATTTCTGCCGATAACCCGGAAATTGACGCCGCTCATGAAGCGCGTATTCCGGTGATCCGTCGTGCCGAAATGCTGGCTGAGTTAATGCGTTTTCGTCATgg                                                                >  1:332837/1‑251 (MQ=255)                                                                                                                                                                                           aaTGAAGGTTATCAGATCAGTGGTTCCGATTTAGCGCCAAATCCGGTCACGCAGCAGTTAACGAATCTGGGTGCGACGATTTATTTCAACCATCGCCCGGAAAACGTACGTGATGCCAGCGTGGTCGTTGTTGCCAGCGCGATTTCTGCCGATAACCCGGAAATTGTCGCCGCTCATGAAGCGCGTATTCCGGTGATCCGTCGTGCCGAAATGCTGGCTGAGTTAATGCGTTTTCGTCATGGCATCGCCAt                                                       <  1:317724/251‑1 (MQ=255)                                                                                                                                                                                             tGAAGGTTATCAGATCAGTGGTTCCGATTTAGCGCCAAATCCGGTCACGCAGCAGTTAACGAATCTGGGTGCGACGATTTATTTCAACCATCGCCCGGAAAACGTACGTGATGCCAGCGTGGTCGTTGTTTCCAGCGCGATTTCTGCCGATAACCCGGAAATTGTCGCCGCTCATGAAGCGCGTATTCCGGTGATCCGTCGTGCCGAAATGCTGGCTGAGTTAATGCGTTTTCGTCATGGCATCGCCATTg                                                     >  1:118131/1‑251 (MQ=255)                                                                                                                                                                                                   ttATCAGATCAGTGGTTCCGATTTAGCGCCAAATCCGGTCACGCAGCAGTTAACGAATCTGGGTGCGACGATTTATTTCAACCATCGCCCGGAAAACGTACGTGATGCCAGCGTGGTCGTTGTTTCCAGCGCGATTTCTGCCGATAACCCGGAAATTGTCGCCGCTCATGAAGCGCGTATTCCGGTGATCCGTCGTGCCGAAATGCTGGCTGAGTTAATGCGTTTTCGTCATGGCATCGCCATTGCCGGa                                                >  1:401150/1‑250 (MQ=255)                                                                                                                                                                                                     atcagatcagTGGTTCCGCTTTAGCGCCAAATCCGGTCAAGCAGCAGTTAACGAATCTGGGTGCGACGATTTATTTCAACCATCGCCCGGAAAACGTACGTGATGCCAGCGTGGTCGTTGATTCCAGCGCGATTTCTGCCGATAACCCGGAAATTGTCGCCGCTCATGAAGCGCGTATTCCGGTGATCCGTCGTGCCGAAATGCTGGCTGAGTTACTGCGTTTTCGTCATGGCATCGCCATTGCCGGAACg                                             >  2:358165/1‑251 (MQ=255)                                                                                                                                                                                                      tcagatcagTGGTTCCGATTTAGCGCCAAATCCGGTCACGCAGCAGTTAACGAATCTGGGTGCGACGATTTATTTCAACCATCGCCCGGAAAACGTACGTGATGCCAGCGTGGTCGTTGTTTCCAGCGCGATTTCTGCCGATAACCCGGAAATTGTCGCCGCTCATGAAGCGCGTATTCCGGTGATCCGTCGTGCCGAAATGCTGGCTGAGTTAATGCGTTTTCGTCATGGCATCGCCATTGCCGGAACGc                                            >  1:234674/1‑251 (MQ=255)                                                                                                                                                                                                                ggTTCCGATGTAGCGCCAAATCCGGTCACGCAGCAGGTAACGAATCTGGGTGCGACGATTTATTTCAACCATCGCCCGGAAAACGTACGTGATGCCAGCGTGGTCGTTGTTTCCAGCGCGATTTCGGCCGATAACCCGGAAATTGTCGCCGCTCATGAAGCGCTTATTCCGGTTATCCGTCGTGCCGAAATGCTGGCTGAGTTAATGCGTTTTCGTTATGGCATCGCCATTGCCGGAACGCACGGCAAAAc                                  <  2:439727/251‑1 (MQ=255)                                                                                                                                                                                                                          tAGCGCCAAATCCGGTCACGCAGCAGTTAACGAATCTGGGTGCGACGATTTATTTCAACCATCGCCCGGAAAACGTACGTGATGCCAGCGTGGTCGTTGTTTCCAGCGCGATTTCTGCCGATAACCCGGAAATTGTCGCCGCTCATGAAGCGCGTATTCCGGTGATCCGTCGTGCCGAAATGCTGGCTGAGTTAATGCGTTTTCGTCATGGCATCGCCATTGCCGGAACGCACGGCAAAACGACAaccacc                        >  1:87651/1‑251 (MQ=255)                                                                                                                                                                                                                               ccAAATCCGGTCCCTCAGCAGTTACCGAATCTAGATGCGACGATTTCTTTCAACCATCGCCCGGAAAACGTCCGTGATGCCAGCGTGGCCGTTGTTTCCAGCGCGATTTCTGCCGATAACCCGGAAATTGTCGCCGCTCATGAAGCGCGTATACCGGTGATCCGACGTGCCGAACAGCTGGCTGAGTTAATGCGTTTTCGTCATGGACTCGCCCTAGCCGGAACGCACGGCCAAACGACAACCCACGCGAt                   >  2:59618/1‑251 (MQ=255)                                                                                                                                                                                                                                cAAATCCGGTCACGCAACAGTTAACGAATCTGGGTGCGACGATTTATTTCAACCATCGCCCGGAAAACGTACGTGATGCCAGCGTGGTCGTTGTTTCCAGCGCGATTTCTGCCGCTAACCCGGAAATTTTCGCCGCTCATGAAGCGCGTATTCCGGTGATCCGTCGTGCCGAAATGCTGGCTGAGTTAATGCGTTTTCGTCATGGCATCTCCATTGCCGGCACGCACGGCAAAACGACACCCACCGCGATg                  >  1:369555/1‑251 (MQ=255)                                                                                                                                                                                                                                            cGCAGCAGTTAACGAATCTGGGTGCGACGATTTATTTCAACCATCGCCCGGAAAACGTACGTGATGCCAGCGTGGTCGTTGTTTCCAGCGCGATTTCTGCCGATAACCCGGAAATTGTCGCCGCTCATGAAGCGCGTATTCCGGTGATCCGTCGTGCCGAAATGCTGGCTGAGTTAATGCGTTTTCGTCATGGCATCGCCATTGCCGGAACGCACGGCAAAACGACAACCACCGCGATGGTTTCCAGCATc      <  2:168710/251‑1 (MQ=255)                                                                                                                                                                                                                                             gcagcaGTTAACGAATCTGGGTGCGACGATTTATTTCAACCATCGCCCGGAAAACGTACGTGATGCCAGCGTGGTCGTTGTTTCCAGCGCGATTTCTGCCGATAACCCGGAAATTGTCGCCGCTCATGAAGCGCGTATTCCGGTGATCCGTCGTGCCGAAATGCTGGCTGAGTTAATGCGTTTTCGTCATGGCATCGCCATTGCCGGAACGCACGGCAAAACGACAACCACCGCGATGGTTTCCAGCATCt     >  1:21471/1‑251 (MQ=255)                                                                                                                                                                                                                                               agcagTTAACGAATCTGGGTGCGACGATTTATTTCAACCATCGCCCGGAAAACGTACGTGATGCCAGCGTGGTCGTTGTTTCCAGCGCGATTTCTGCCGATAACCCGGAAATTGTCGCCGCTCATGAAGCGCGTATTCCGGTGATCCGTCGTGCCGAAATGCTGGCTGAGTTAATGCGTTTTCGTCATGGCATCGCCATTGCCGGAACGCACGGCAAAACGACAACCACCGCGATGGTTTCCAGCATCTAc   <  1:84556/251‑1 (MQ=255)                                                                                                                                                                                                                                                gcagTTAACGAATCTGGGTGCGACGATTTATTTCAACCATCGCCCGGAAAACGTACGTGATGCCAGCGTGGGCGTTGTTTCCAGCGCGATTTCTGCCGATAACCCGGAAATTGTCGCCGCTCATGAAGCGCGTATTCCGGTGATCCGTCGTGCCGAAATGCTGGCTGAGTTAATGCGTTTTCGTCATGGCATCGCCATTGCCGGAACGCACGGCAAAACGACAACCACCGCGATGGTTTCCAGCATCTACg  <  2:185792/251‑1 (MQ=255)                                                                                                                                                                                                                                                        |                                                                                                                                                                                                                                                   GCCCGGGCGTAATTGTAGCGATGCCTTTTGCATCGTATGAATTTAAGAAGTTAATGGCGTAAAGAATGAATACACAACAATTGGCAAAACTGCGTTCCATCGTGCCCGAAATGCGTCGCGTTCGGCACATACATTTTGTCGGCATTGGTGGTGCCGGTATGGGCGGTATTGCCG‑AAGTTCTGGCCAATGAAGGTTATCAGATCAGTGGTTCCGATTTAGCGCCAAATCCGGTCACGCAGCAGTTAATGAATCTGGGTGCGACGATTTATTTCAACCATCGCCCGGAAAACGTACGTGATGCCAGCGTGGTCGTTGTTTCCAGCGCGATTTCTGCCGATAACCCGGAAATTGTCGCCGCTCATGAAGCGCGTATTCCGGTGATCCGTCGTGCCGAAATGCTGGCTGAGTTAATGCGTTTTCGTCATGGCATCGCCATTGCCGGAACGCACGGCAAAACGACAACCACCGCGATGGTTTCCAGCATCTACG  >  NC_000913/100700‑101188

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