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	95.1%	123.0 / ‑5.3	41	M61T (ATG→ACG)	murC	UDP‑N‑acetylmuramate‑‑L‑alanine ligase
Reads supporting (aligned to +/- strand):  ref base T (0/0);  major base C (25/14);  minor base A (1/1);  total (26/15)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.90e-01

CCGGGCGTAATTGTAGCGATGCCTTTTGCATCGTATGAATTTAAGAAGTTAATGGCGTAAAGAATGAATACACAACAATTGGCAAAACTGCGTTCCATCGTGCCCGAAATGCGTCGCGTTCGGCACATACATTTTGTCGGCATTGGTGGTGCCGGTATGGGCGGTATTGCCGAAGTTCTGGCCAATGAAGGTTATCAGATCAGTGGTTCCGATTTAGCGCCAAATCCGGTCACGCAGCAGTTAATGAATCTGGGTGCGACGATTTATTTCAACCATCGCCCGGAAAACGTACGTGATGCCAGCGTGGTCGTTGTTTCC‑AGCGCGATTTCTGCCGATAACCCGGAAATTGTCGCCGCTCATGAAGCGCGTATTCCGGTGATCCGTCGTGCCGAAATGCTGGCTGAGTTAATGCGTTTTCGTCATGGCATCGCCATTGCCGGAACGCACGGCAAAACGACAACCACCGCGATGGTTTCC  >  NC_000913/100702‑101178                                                                                                                                                                                                                                                     |                                                                                                                                                                                                                                          ccGGGCGTAATTGTAGCGATGCCTTTTGCATCGTATGAATTTAAGAAGTTACTGGCGTAAAGAATGAATACACAACAATTGGCAAAACTGCGTTCCATCGTGCCCGAAATGCGTCGCGTTCGGCACATACATTTTGTCGGCATTGGTGGTGCCGGTATGGGCGGTATTGCCGAAGTTCTGGCCAATGAAGGTTATCAGATCAGTGGTTCCGATTTAGCGCCAAATCCGGTCACGCAGCAGTTAACGAATCt                                                                                                                                                                                                                                     >  1:173783/1‑251 (MQ=255)      cGTAATTGTAGCGATGCCTTTTGCATCGTATGAATTTAAGAAGTGAATGGCGTAAAGAATGAATACACAACAATTGGCAAAACTGCGTTCCATCGTGCCCGAAATGCGGCGCGTTCGGCACATACATTTTGTCGGCATTGGTGGTGCCGGTATGGGCGGTATTGCCGAAGGTCTGGCCAATGAAGGTTATCAGATCAGTGGTTCCGATTTAGCGCCAAATCCGGTTACGCAGCAGTTAACGAATCTGGGTg                                                                                                                                                                                                                                <  2:242698/251‑1 (MQ=255)           ttGTAGCGATGCCTTTTGCATCGTATGAATTTAAGAAGTTAATGGCGTAAAGAATGAATACACAACAATTGGCAAAACTGCGTTCCATCGTGCCCGAAATGCGTCGCGATCGGCACATACATTTTGTCGGCATTGGTGGTGCCGGTATGGGCGGGATTGCCGAAGTTCTGGCCAATGAAGGTTATCAGATCAGTGGTTCCGATTTAGCGCCAAATCCGGTCACGCAGCAGTTAACGAATCTGGGTGcgacg                                                                                                                                                                                                                           >  2:157888/1‑251 (MQ=255)              aagagATGACTTTTGAATCGTAAGAAATTAAGAAGTGAATGGTGTAAAGAATGAATACACAACAATTGGCAAAACTGCGTTCCTTCGTGCCCGAAATGCGTCGCGTTCGGCACATACATTTTGTCGGCATTGGTGGTGCCGGTATGGGCGGTATTGCCGAAGTTCTGGCCAATGAAGGTTATCAGATCAGTGGTTCCGATTTAGCGCCAAATCCGGTCACACAGCAGTTAACGAATCTGGGTGCGACGatt                                                                                                                                                                                                                        <  2:110924/247‑1 (MQ=255)               aGCGATGCCTTTTGCATCGTATGAATTTAAGAAGTTAATGGCGTAAAGAATGAATACACAACAATTGGCAAAACTGCGTTCCATCGTGCCCGAAATGCGTCGCGTTCGGCACATACATTTTGTCGGCATTGGTGGTGCCGGTATGGGCGGTATTGCCGAAGTTCTGGCCAATGAAGGTTATCAGATCAGTGGTTCCGATTTAGCGCCAAATCCGGTCACGCAGCAGTTAACGAATCTGGGTGCGACGattt                                                                                                                                                                                                                       <  1:181201/251‑1 (MQ=255)                 cGATGCCTTTTGCATCGTATGAATTTAAGAAGTTAATGGCGTAAAGAATGAATACACAACAATTGGCAAAACTGCGTTCCATCGTGCCCGAAATGCGTCGCGTTCGGCACATACATTTTGTCGGCATTGGTGGTGCCGGTATGGGCGGTATTGCCGAAGTTCTGGCCAATGAAGGTTATCAGATCAGTGGTTCCGATTTAGCGCCAAATCCGGTCACGCAGCAGTTAACGAATCTGGGTGCGACGatttat                                                                                                                                                                                                                     >  1:124976/1‑251 (MQ=255)                  gATGCCTTTTGCATCGTATGAATTTAAGAAGTTAATGGCGTAAAGAATGAATACACAACAATTGGCAAAACTGCGTTCCATCGTGCCCGAAATGCGTCGCGTTCGGCACATACATTTTGTCGGCATTGGTGGTGCCGGTATGGGCGGTATTGCCGAAGTTCTGGCCAATGAAGGTTATCAGATCAGTGGTTCCGATTTAGCGCCAAATCCGGTCACGCAGCAGTTAACGAATCTGGGTGCGACGatttatt                                                                                                                                                                                                                    <  1:126932/251‑1 (MQ=255)                             catcgcatgcctttAAGAAGTTAATGGCGTAAAGAATGAATACACAACAATTGACAAAACTGCGTTCCATCGTCCCCGAAATGCGTAGCGTTCGGCACATACATTTTGTCGGCATTGGTGGTGCCGGTATGGGCGGTATTGCCGAAGTTCTGGCCAATGAAGGTTATCAGATCAGTGTTTCAGCTTTAGCGCCAAATACGGTCACGCAGCAGTTAACGAATATTGGTGCGACGATTTATTTCAACCATCGc                                                                                                                                                                                                         >  2:463208/12‑251 (MQ=255)                             cATCGTATGAATTTAAGATGTTAATGGCGTAAAGAATGAATACACAACAATTGGCAAAACTGCGTTCCATCGTGCCCGAAATGCGTCGCGTTCGGCACATACATTTTGTCGGCATTGGTGGTGCCGGTATGGGCGGTATTGCCGAAGTTCTGGCCAATGAAGGTTATCAGATCAGTGGTTCCGATTTAGCGCCAACTCCGGTCACGCAGCAGTTAACGAATCTGGGTGCGACGATTTATTTCAACCATCGc                                                                                                                                                                                                         >  1:356872/1‑251 (MQ=255)                             cATCGTATGAATTTAAGAAGTTAATGGCGTAAAGAATGAATACACAACAATTGGCAAAACTGCGTTCCATCGTGCCCGAAATGCGTCGCGTGCGGCACATACATTTTGTCGGCATTGGTGGTGCCGGTATGGGCGGTATTGCCGAAGTTCTGGCCAATGAATGTTATCAGATCAGTGGTTCCGATTGAGCGCCAAATCCGGTCACGCAGCAGTTAACGAATCTGGGTGCGACGATTTATTTCAACCATCGc                                                                                                                                                                                                         <  2:34474/251‑1 (MQ=255)                              aTCGTATGAATTTAAGAAGTTAATGGCGTAAAGAATGAATACACAACAATTGGCAAAACTGCGTTCCATCGTGCCCGAAATGCGTCGCGTTCGGCACATACATTTTGTCGGCATTGGTGGTGCCGGTATGGGCGGTATTGCCGAAGTTCTGGCCAATGAAGGTTATCAGATCAGTGGTTCCGATTTAGCGCCAAATCCGGTCACGCAGCAGTTAACGAATCTGGGTGCGACGATTTATTTCAACCATCGcc                                                                                                                                                                                                        >  1:460420/1‑251 (MQ=255)                                         ttAAGAAGTTAATGGCGTAAAGAATGAATACACAACAATTGGCAAAACTGCGTTCCATCGTGCCCGAAATGCGTCGCGTTCGGCACATACATTTTGTCGGCATTGGTGGTGCCGGTATGGGCGGTATTGCCGAAGTTCTGGCCAATGAAGGTTATCAGATCAGTGGTTCCGATTTAGCGCCAAATCCGGTCACGCAGCAGTTAACGAATCTGGGTGCGACGATTTATTTCAACCATCGCCCGGAAAacgt                                                                                                                                                                                              >  1:427180/1‑250 (MQ=255)                                         ttAAGAAGTTAATGGCGTAAAGAATGAATACACAACAATTGGCAAAACTGCGTTCCATCGTGCCCGAAATGCGTCGCGTTCGGCACATACATTTTGTCGGCATTGGTGGTGCCGGTATGGGCGGTATTGCCGAAGTTCTGGCCAATGAAGGTTATCAGATCAGTGGTTCCGATTTAGCGCCAAATCCGGTCACGCAGCAGTTAAAGAATCTGGGTGCGACGATTTATTTAAACCATAGCCCGGAAAacgtc                                                                                                                                                                                             >  2:6710/1‑250 (MQ=255)                                                  tAATGGCGTAAAGAATGAATACACAACAATTGGCAAAACTGCGTTCCATCGTGCCCGAAATGCGTCGCGTTCGGCACATACATTTTGTCGGCATTGGTGGTGCCGGTATGGGCGGTATTGCCGAAGTTCTGGCCAATGAAGGTTATCAGATCAGTGGTTCCGATTTAGCGCCAAATCCGGTCACGCAGCAGTTAACGAATCTGGGTGCGACGATTTATTTCAACCATCGCCCGGAAAACGTACGTGATGcc                                                                                                                                                                                    >  1:268790/1‑251 (MQ=255)                                                   aaTGGCGTAAAGAATGAATACACAACAATTGGCAAAACTGCGTTCCATCGTGCCCGAAATGCGTCGCGTTCGGCACATACATTTTGTCGGCATTGGTGGTGCCGGTATGGGCGGTATTGCCGAAGTTCTGGCCAATGAAGGTTATCAGATCAGTGGTTCCGATTTAGCGCCAAATCCGGTCACGCAGCAGTTAACGACTCTGGGGGCGACGATTTATTTCAACCATCGCCCGGAAAACGTACGTGCTGcc                                                                                                                                                                                    >  1:88454/1‑250 (MQ=255)                                                      ggCGTAAAGAATGAATACACAACAATTGGCAAAACTGCGTTCCATCGTGCCCGAAATGCGTCGCGTTCGGCACATACATTTTGTCGGCATTGGTGGTGCCGGTATGGGCGGTATTGCCGAAGTTCTGGCCAATGAAGGTTATCAGATCAGTGGTTCCGATTTAGCGCCAAATCCGGTCACGCAGCAGTTAACGAATCTGGGTGCGACGATTTATTTCAACCATCGCCCGGAAAACGTACGTGATGCCAGCg                                                                                                                                                                                >  1:195287/1‑251 (MQ=255)                                                         gTAAAGAATGAATACACAACAATTGGCAAAACTGCGTTCCATCGTGCCCGAAATGCGTCGCGTTCGGCACATACATTTTGTCGGCATTGGTGGTGCCGGTATGGGCGGTATTGCCGAAGTTCTGGCCAATGAAGGTTATCAGATCAGTGGTTCCGATTTAGCGCCAAATCCGGTCACGCAGCAGTTAACGAATCTGGGTGCGACGATTTATTTCAACCATCGCCCGGAAAACGTACGTGATGCCAGCGTgg                                                                                                                                                                             >  1:336903/1‑251 (MQ=255)                                                             aGAATGAATACACAACAATTGGCAAAACTGCGTTCCATCGTGCCCGAAATGCGTCGCGTTCGGCACATACATTTTGTCGGCATTGGTGGTGCCGGTATGGGCGGTATTGCCGAAGTTCTGGCCAATGAAGGTTATCAGATCAGTGGTTCCGATTTAGCGCCAAATCCGGTCACGCAGCAGTTAACGAATCTGGGTGCGACGATTTATTTCAACCATCGCCCGGAAAACGTACGTGATGCCAGCGTGGTCgt                                                                                                                                                                         <  2:125805/251‑1 (MQ=255)                                                               aatgaatACACAACAATTGGCAAAACTGCGTTCCATCGTGCCCGAAATGCGTCGCGTTCGGCACATACATTTTGTCGGCATTGGTGGTGCCGGTATGGGCGGTATTGCCGAAGTTCTGGCCAATGAAGGTTATCAGATCAGTGGTTCCGATTTAGCGCCAAATCCGGTCACGCAGCAGTTAACGAATCTGGGTGCGACGATTTATTTCAACCATCGCCCGGAAAACGTACGTGATGCCAGCGTGGTCgttg                                                                                                                                                                       >  1:318675/1‑251 (MQ=255)                                                               aatgaatACACAACAATTGGCAAAACTGCGTTCCATCGTGCCCGAAATGCGTCGCGTTCGGCACATACATTTTGTCGGCATTGGTGGTGCCGGTATGGGCGGTATTGCCGAAGTTCTGGCCAATGAAGGTTATCAGATCAGTGGTTCCGATTTAGCGCCAAATCCGGTCACGCAGCAGTTAACGAATCTGGGTGCGACGATTTATTTCAACCATCGCCCGGAAAACGTACGTGATGCCAGCGTGGTCgttg                                                                                                                                                                       >  1:224637/1‑251 (MQ=255)                                                                atgaatATACAACAATTGGCAAAACTGCGTTCCATCGTCCCCGAAATGCGTCGCGATCGGCACATACATTTTGTCGGCATTGGTGGTGCCGGTATGGGCGGTATTGCCGAAGTTCTGGCCAATGAAGGTTATCAGATCAGTGGTTCCGATTTAGCGCCAAATCCGGTCACGAAGCAGTTAACGCATCTGGGGGCGACGATTTATTTCAACCATCGCCCGGAAAAAGTCCGTGATGCCAGCGTGATCgttgt                                                                                                                                                                      >  2:197447/1‑251 (MQ=255)                                                                atgaatACACAACAATTGGCAAAACTGCGTTCCATCGTGCCCGAAATGCGTCGCGTTCGGCACATACATTTTGTCGGCATTGGTGGTGCCGGTATGGGCGGTATTGCCGAAGTTCTGGCCAATGAAGGTTATCAGATCAGTGGTTCCGATTTAGCGCCAAATCCGGTCACGCAGCAGTTAACGAATCTGGGTGCGACGATTTATTTCAACCATCGCCCGGAAAACGTACGTGATGCCAGCGTGGTCgttgt                                                                                                                                                                      <  1:19593/251‑1 (MQ=255)                                                                    atACACAACAATTGGCAAAACTGCGTTCCATCGTGCCCGAAATGCGTCGCGTTCGGCACATACATTTTGTCGGCATTGGTGGTGCCGGTATGGGCGGTATTGCCGAAGTTCTGGCCAATGAAGGTTATCAGATCAGTGGTTCCGATTTAGCGCCAAATCCGGTCACGCAGCAGTTAACGAATCTGGGTGCGACGATTTATTTCAACCATCGCCCGGAAAACGTACGTGATGCCAGCGTGGTCGTTGTTTcc                                                                                                                                                                  <  1:104002/251‑1 (MQ=255)                                                                    atACACAACAATTGGCAAAACTGCGTTCCATCGTGCCCGAAATGCGTCGCGTTCGGCACATACATTTTGTCGGCATTGGTGGTGCCGGTATGGGCGGTATTGCCGAAGTTCTGGCCAATGAAGGTTATCAGATCAGTGGTTCCGATTTAGCGCCAAATCCGGTCACGCAGCAGTTAACGAATCTGGGTGCGACGATTTATTTCAACCATCGCCCGGAAAACGTACGTGATGCCAGCGTGGTCGTTGTTTcc                                                                                                                                                                  >  1:165504/1‑251 (MQ=255)                                                                      acacAACAATTGGCAAAACTGCGTTCCATCGTGCCCGAAATGCGTCGCGTTCGGCACATACATTTTGTCGGCATTGGTGGTGCCGGTATGGGCGGTATTGCCGAAGTTCTGGCCAATGAAGGTTATCAGATCAGTGGTTCCGATTTAGCGCCAAATCCGGTCACGCAGCAGTTAACGAATCTGGGTGCGACGATTTATTTCAACCATCGCCCGGAAAACGTACGTGATGCCAGCGTGGTCGTTGTTTCC‑a                                                                                                                                                                <  1:412067/250‑1 (MQ=255)                                                                                   caaaatagccttCCATCGTACCCGAATTGCGTCGCGTTCGGCACATACATTTTGTCGGCATTGGTGGCGCCGGTATGGGCGGTATTGCCGAAGTTCTGGCCAATGAAGGTTGTCAGATCCGTGGTTCCGATTTAGCGACAAATCCGGTCACGCAGCAGTTAACGAATCTGGGGGCGACGGTTTATTTCAACCCTCGCCCGGCAAACGTGCGTGCTGCCAGCATGGTCGTTGTTTCC‑AGCGCGTTTTCTGcc                                                                                                                                                  >  2:222287/11‑251 (MQ=255)                                                                                           cGTTCCAGCGGGCCCGAAATGCGTCGCGGTCGGCGCATACATTTTGTCGGCATTGGTGTTGCCGGTATGGGCGGTATTGCCGAAGTTCTGGCCAATGAAGGTTATCAGATCAGTGGGTCCGATTTAGCGCCAAATCCGTTCACGCAGCAGTTATCGAATCTGGGGGCGACGATTTATTGCAACCTACGCTCGGTAAACGTACGTGTTGCCAGCGTGGTCGT‑GTTTCCCAGCGTGATTGCTGCCGAGAAccc                                                                                                                                          <  2:483068/251‑1 (MQ=255)                                                                                                 aTCGTACCCGAAATGCGTCGCGTTCGGCACATACATTTTGTCGGCATTGGTGGTGCCGGTATGGGCGGTATTGCCGAAGTTCTGGCCAATGAAGGTTATCAGATCAGTGGTTCCGATTTAGCGCCAAATCCGGTCACGCAGCAGTTACCGAATCTGGGTGCGACGATTTATTTCAACCATCGCCCGGAAAACGTACGTGATGCCCGCGTGGTCGTTGTTTCC‑AGCGCGATTTCTGCCGATCACCCGGAAAt                                                                                                                                    >  1:181054/1‑251 (MQ=255)                                                                                                                                  catTTTGTCGGCATTGGTGGTGCCGGTATGGGCGGTATTGCCGAAGTTCTGGCCAATGAAGGTTATCAGATCAGTGGTTCCGATTTAGCGCCAAATCCGGTCACGCAGCAGTTAACGAATCTGGGTGCGACGATTTATTTCAACCATCGCCCGGAAAACGTACGTGATGCCAGCGTGGTCGTTTTTTCC‑AGCGCGATTTCGGCCGATAACCCGGAAATTGTCGCCGCTCATGAAGCGCGGATTCCGGTGAt                                                                                                   <  2:439762/251‑1 (MQ=255)                                                                                                                                                 ggtggtGCCGGTATGGGCGGTATTGCCGAAGTTCTGGCCAATGAAGGTTATCAGATCAGTGGTTCCGATTTAGCGCCAAATCCGGTCACGCAGCAGTTAACGAATCTGGGTGCGACGATTTATTTCAACCATCGCCCGGAAAACGTACGTGATGCCAGCGTGGTCGTTGTTTCC‑AGCGCGATTTCTGCCGATAACCCGGAAATTGTCGCCGCTCATGAAGCGCGTATTCCGGTGATCCGTCGTGCCGAAAt                                                                                    <  1:287134/251‑1 (MQ=255)                                                                                                                                                                   ggTATTGCCGAAGTTCTGGCCAATGAAGGTTATCAGATCAGTGGTTCCGATTTAGCGCCAAATCCGGTCACGCAGCAGTTAACGAATCTGGGTGCGACGATTTATTTCAACCATCGCCCGGAAAACGTACGTGATGCCAGCGTGGTCGTTGTTTCC‑AGCGCGATTTCTGCCGATAACCCGGAAATTGTCGCCGCTCATGAAGCGCGTATTCCGGGGATCCGTCGTGCCGAAATGCTGGCTGAGTTAATGCg                                                                  >  2:215080/1‑251 (MQ=255)                                                                                                                                                                   ggTATTGCCGAAGTTCTGGCAAATGAAGGTTATCAGATCAGTGGTTCCGATTTAGCGCCAAATCCGGTCACGCAGCAGTTAACGAATCTGGGTGCGACGATTTATTTCAACCATCGCCCGGAAAACGTACGTGATGCCAGCGTGGTCGTTGTTTCC‑AGCGCGATTTCTGCCGATAACCCGGAAATTGTCGCCGCTCATGAAGCGCGTATTCCGGTGATCCGTCGTGCCGACATGCTGGCTGAGTTAATGCg                                                                  >  1:415298/1‑251 (MQ=255)                                                                                                                                                                                 tCTGGCCAATGAAGGTTATCAGATCAGGGGTTCCGATTTAGCGCCAAATCCGGTCACGCAGCAGTTAACGAATCTGGGTGCGACGATTTATTTCAACCATCGCCCGGAAAACGTACGTGATGCCAGCGTGGTCGTTGTTTCC‑AGCGCGATTTCTGCCGATAACCCGGAAATTGTCGCCGCTCATGAAGCGCGTATTCCGGTGATCCGTCGTGCCGAAATGCTGGCTGAGTTAATGCGTTTTCGTCATGGc                                                     <  2:176873/250‑1 (MQ=255)                                                                                                                                                                                  cTGGCCAATGAAGGTTATCAGATCAGTGGTTCCGATTTAGCGCCAAATCCGGTCACGCAGCAGTTAACGAATCTGGGTGCGACGATTTATTTCAACCATCGCCCGGAAAACGTACGTGATGCCAGCGTGGTCGTTGTTTCC‑AGCGCGATTTCTGCCGATAACCCGGAAATTGTCGCCGCTCATGAAGCGCGTATTCCGGTGATCCGTCGTGCCGAAATGCTGGCTGAGTTAATGCGTTTTCggcatggcat                                                   >  2:427085/1‑251 (MQ=255)                                                                                                                                                                                        aaTGAAGGTTAGCAGATCAGTGGTTCCGATTTAGCGCAAAAGCCGGTCACGCAGCAGTTAAAGAATCTGCGTGCGACGATTTATTTCAACCATCGCCCGGAAAACTTACGTGATGCCAGCGTGGTCGTTGTTTCC‑AGCGCGATTTCTTCCGATAACCCGGAAACTGTCGCCGCTCATGAAGCGCGGATTCCGGTGATCCGCCGTGCCGAAACGCTGGCTGAGTTAATGCGTGTTCGTCCTTGCATCGCCAt                                             <  2:391775/251‑1 (MQ=255)                                                                                                                                                                                         aTGAAGGTTATCAGATCAGTGGTTCCGATTTAGCGCCAAATCCGGTCACGCAGCAGTTAACGAATCTGGGTGCGACGTTTTATTTCAACCATCGCCCGGAAAACGTACGTGATGCCAGCGTGGTCGTTGTTTCC‑AGCGCGATTTCTGCCGATAACCCGGAAATTGTCGCCGCTCATGAAGCGCGTATTCCGGTGATCCGTCGTGCCGAAATGCTGGCTGAGTTAATGCGTTTTCGTCATGGCATCGccctt                                            >  1:333257/1‑248 (MQ=255)                                                                                                                                                                                               gTTATCAGATCAGTGGTTCCGATTTAGCGCCAAATCCGGTCACGCAGCAGTTAACGAATCTGGGTGCGACGATTTATTTCAACCATCGCCCGGAAAACGTACGTGATGCCAGCGTGGTCGTTGTTTCC‑AGCGCGATTTCTGCCGATAACCCGGAAATTGTCGCCGCTCATGAAGCGCGTATTCCGGTGATCCGTCGTGCCGAAATGCTGGCGGGGCTAATGCGTTTTCGTCATGGCATCGCCATTGGCgg                                       >  1:478506/1‑250 (MQ=255)                                                                                                                                                                                                   tcagatcagTGGTTCCGATTTAGCGCCAAATCCGGTCACGCAGCAGTTAACGAATCTGGGTGCGACGATTTATTTCAACCATCGCCCGGAAAACGTACGTGATGCCAGCGTGGTCGTTGTTTCC‑AGCGCGATTTCTGCCGATAACCCGGAAATTGTCGCCGCTCATGAAGCGCGTATTCCGGTGATCCGTCGTGCCGAAATGCTGGCTGAGTTAATGCGTTTTCGTCATGGCATCGCCATTGCCGGAACGc                                  >  2:112333/1‑251 (MQ=255)                                                                                                                                                                                                           gTGGTTCCGATTGAGCGCCAAATCCGGTCACGCATCAGTGAACGAATCTGGGTGCGACGATTTATTTCAACCATCGCCCGGAAAACGTACGGGATGCCAGCGTGGCCGTTGTTTCC‑AGCGCGATTTCTGCCGATAACCCGGAAATTGTCGCTGCTCAAGAAGCGCGTATTCCGGTGATCCGTCGTGCCGAAATGCTGGCTGAGTTAATGCGTTGTCGTCATGGCAGCGCCATTGCCGGAACGCACGGCaaa                          <  2:181054/251‑1 (MQ=255)                                                                                                                                                                                                                   gATTTAGCGCCAAATCCGGTCACGCAGCAGTTAACGAATCTGGGTGCGACGATTTATTTCAACCATCGCCCGGAAAACGTACGTGATGCCAGCGTGGTCGTTGTTTCC‑AGCGCGATTTCTGCCGATAACCCGGAAATTGTCGCCGCTCATGAAGCGCGTATTCCGGTGATCCGTCGTGCCGAAATGCTGGCTGAGTTAATGCGTTTTCGACATGGCATCGCCATTGCCGGAACGCACGGCAAAACGACAac                  >  1:270335/1‑251 (MQ=255)                                                                                                                                                                                                                                     gTCACGCAGCAGTTAACGAATCTGGGTGCGACGATTTATTTCAACCATCGCCCGGAAAACGTACGTGATGCCAGCGTGGTCGTTGTTTCC‑AGCGCGATTTCTGCCGATAACCCGGAAATTGTCGCCGCTCATGAAGCGCGTATTCCGGTGATCCGTCGTGCCGAAATGCTGGCTGAGTTAATGCGTTTTCGTCATGGCATCGCCCTTGCCGGAACGCACGGCAAACCGACAACCACCGCGATGGTTTcc  >  2:418769/1‑249 (MQ=255)                                                                                                                                                                                                                                                     |                                                                                                                                                                                                                                          CCGGGCGTAATTGTAGCGATGCCTTTTGCATCGTATGAATTTAAGAAGTTAATGGCGTAAAGAATGAATACACAACAATTGGCAAAACTGCGTTCCATCGTGCCCGAAATGCGTCGCGTTCGGCACATACATTTTGTCGGCATTGGTGGTGCCGGTATGGGCGGTATTGCCGAAGTTCTGGCCAATGAAGGTTATCAGATCAGTGGTTCCGATTTAGCGCCAAATCCGGTCACGCAGCAGTTAATGAATCTGGGTGCGACGATTTATTTCAACCATCGCCCGGAAAACGTACGTGATGCCAGCGTGGTCGTTGTTTCC‑AGCGCGATTTCTGCCGATAACCCGGAAATTGTCGCCGCTCATGAAGCGCGTATTCCGGTGATCCGTCGTGCCGAAATGCTGGCTGAGTTAATGCGTTTTCGTCATGGCATCGCCATTGCCGGAACGCACGGCAAAACGACAACCACCGCGATGGTTTCC  >  NC_000913/100702‑101178

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: ‑