Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	NC_000913	4,414,625	C→T	R117R (CGC→CGT)	aidB →	putative acyl‑CoA dehydrogenase AidB

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	4,414,625	0	C	T	97.1%	202.8 / NA	70	R117R (CGC→CGT)	aidB	putative acyl‑CoA dehydrogenase AidB
Reads supporting (aligned to +/- strand):  ref base C (0/0);  major base T (36/32);  minor base G (0/2);  total (36/34)
Fisher's exact test for biased strand distribution p-value = 2.32e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

TGGCTGGGATAGCGATTTTCTAGCCAGTATTGGTCAGCAGTTAGGAACGGCTGAATCCCTTGAACTGGGGCGGCTGGCGAATGTGAATCCGCCTGAATTATTGCGCTACGATGCGCAAGGACGCCGTCTGGACGATGTGCGTTTTCACCCCGCCTGGCACCTGCTGATGCAGGCGCTATGTACCAATCGGGTGCACAATCTTGCCTGGGAAGAAGACGCTCGCTCCGGCGCATTTGTGGCGCGCGCGGCGCGTTTTATGTTACATGCGCAGGTTGAGGCAGGGTCGTTATGTCCGATAACCATGACCTT‑TG‑CCGCCAC‑GCCATTGTTGTT‑A‑CAGATGTTACCCGCGCCGTTTCAGGACTGGACCACGCCGCTATTGAGCGATCGCTACGATTCTCACTTATTGCCAGGTGGGCAAAAACGCGGTTTGTTGATTGGCATGGGAATGACGGAAAAGCAGGGCGGTTCCGATGTTATGAGCAACACCACCCGT  >  NC_000913/4414382‑4414871                                                                                                                                                                                                                                                    |                                                                                                                                                                                                                                                            tggctggGATAGCGATTTTCTAGCCAGTATTGGTCAGCAGTTAGGAACGGCTGAATCCCTTGAACTGGGGCGGCTGGCGAATGTGAATCCGCCTGAATTATTGCGCTACGATGCGCAAGGACGCCGTCTGGACGATGTGCGTTTTCACCCCGCCTGGCACCTGCTGATGCAGGCGCTATGTACCAATCGGGTGCACAATCTTGCCTGGGAAGAAGACGCTCGCTCCGGCGCATTTGTGGCGCgtgcggcgg                                                                                                                                                                                                                                                      >  1:223776/1‑250 (MQ=255)             cGATTTTCTAGCCAGTATTGGTCAGCAGTTAGGAACGGCTGAATCACTTGAACTGGGGCGGCTGGCGAATGTGAATCCGCCTGAATTATTGCGCTACGATGCGCACGGACGCCGTCTGGACGATGTGCGTTTTCACCCCGCCTGGCACCTGCTGATGCAGGCGCTATGAACCAATCGGGTGCACAATCTTGCCTGGGACGAAGACGCTCGCTCCGGCGCATTTGTGGCCCGTGCGGCGCGTTTGATGTTAc                                                                                                                                                                                                                                          >  1:326283/1‑251 (MQ=255)                  ttCTAGCCAGTATTGGTCAGCAGTTAGGAACGGCTGAATCCCTTGAACTGGGGCGGCTGGCGAATGTGAATCCGCCTGAATTATTGCGCTACGATGCGCAAGGACGCCGTCTGGACGATGTGCGTTTTCACCCCGCCTGGCACCTGCTGATGCAGGCGCTATGTACCAATCGGGTGCACAATCTTGCCTGGGAAGAAGACGCTCGCTCCGGCGCATTTGTGGCGCGTGCGGCGCGTTTTATGTTACATgc                                                                                                                                                                                                                                      <  1:492031/250‑1 (MQ=255)                             attggtctggagtatggAACGGCTGAATTCCTTGAATGGGGGCGGCTTGCGAATGTGAAGCCGCGTTAATGATTGCGCTACGATGCGCAAGGACGCCGTCTGGACGATGTGCGTTTTCACCCCGCCTGGCACCTGCTGATTCAGGCGTTATTTTCCAATCGGGGGCACAATCTTGCCTGGTAAGAAGACTCTCGCTCCGGCGCATTTGTGGCGCGTGCGGCGGGTTTTATGTTCCAGGCGCAGGTTGAgg                                                                                                                                                                                                                           <  2:386550/235‑1 (MQ=255)                             aTTGGTCAGCAGTTAGGAACGGCTGAATCCCTTGAACTGGGGCGGCTGGCGAATGTGAATCCGCCTGAATTATTGCGCTACGATGCGCAAGGACGCCGTCTGGACGATGTGCGTTTTCACCCCGCCTGGCACCTGCTGATGCAGGCGCTATGTACCAATCGGGTGCACAATCTTGCCTGGGAAGAAGACGCTCGCTCCGGCGCATTTGTGGCGCGTGCGGCGCGTTTTATGTTACATGCGCAGGTTGAGGc                                                                                                                                                                                                                          <  1:246160/251‑1 (MQ=255)                                ggTCAGCAGTTAGGAACGGCTGAATCCCTTGAACTGGGGCGGCTGGCGAATGTGAATCCGCCTGAATTATTGCGCTACGATGCGCAAGGACGCCGTCTGGACGATGTGCGTTTTCACCCCGCCTGGCACCTGCTGATGCAGGCGCTATGTACCAATCGGGTGCACAATCTTGCCTGGGAAGAAGACGCTCGCTCCGGCGCATTTGTGGCGCGTGCGGCGGGTTTTATGTTACATGCGCAGGTTGAGGCAgg                                                                                                                                                                                                                       >  1:245866/1‑251 (MQ=255)                                ggTCAGCAGTTAGGAACGGCTGAATCCCTTGAACTGGGGCGGCTGGCGAATGTGAATCCGCCTGAATTATTGCGCTACGATGCGCAAGGACGCCGTCTGGACGATGTGCGTTTTCACCCCGCCTGGCACCTGCTGATGCAGGCGCTATGTACCAATCGGGTGCACAATCTTGCCTGGGAAGAAGACGCTCGCTCCGGCGCATTTGTGGCGCGTGCGGCGCGTTTTATGTTACATGCGCAGGTTGAGGCAg                                                                                                                                                                                                                        <  1:191138/250‑1 (MQ=255)                                  gcagcagGTAGGAACGGCTGAATCCCTTTAACGGGGGCGGCTGGCGAATGTGAATCCGCCTTAATTATTGTGCGACGATGCGCAAGGACGCCGTCTGGACGATGTGCGTTTTCCCCCCGCCTGGCACCTGCTGATGCAGGCGCTATGTACCAATCGGGTGCACAATCTGGCCTGGGAAGAAGACGCTCGCTCCGGCGCATTTGTGGCGCGTGCGGCGCGTTTTATGTTACATGCGCAGGTTGAGGCAGGGt                                                                                                                                                                                                                     <  1:516101/250‑1 (MQ=255)                                      cagTTAGGAACGGCTGAATCCCTTGAACTGGGGCGGCTGGCGAATGTGAATCCGCCTGAATTATTGCGCTACGATGCGCAAGGACGCCGTCTGGACGATGTGCGTTTTCACCCCGCCTGGCACCTGCTGATGCAGGCGCTATGTACCAATCGGGTGCACAATCTTGCCTGGGAAGAAGACGCTCGCTCCGGCGCATTTGTGGCGCGTGCGGCGCGTTTTATGTTACATg                                                                                                                                                                                                                                       >  1:267410/1‑229 (MQ=255)                                      cagTTAGGAACGGCTGAATCCCTTGAACTGGGGCGGCTGGCGAATGTGAATCCGCCTGAATTATTGCGCTACGATGCGCAAGGACGCCGTCTGGACGATGTGCGTTTTCACCCCGCCTGGCACCTGCTGATGCAGGCGCTATGTACCAATCGGGTGCACAATCTTGCCTGGGAAGAAGACGCTCGCTCCGGCGCAGTTGTGGCGCGTGCGGCGCGTTTTATGTTACATg                                                                                                                                                                                                                                       <  2:267410/229‑1 (MQ=255)                                          tAGGAACGGCTGAATCCCTTGAACTGGGGCGGCTGGCGAATGTGAATCCGCCTGAATTATTGCGCTACGATGCGCAAGGACGCCGTCTGGACGATGTGCGTTTTCACCCCGCCTGGCACCTGCTGATGCAGGCGCTATGTACCAATCGGGTGCACAATCTTGCCTGGGAAGAAGACGCTCGCTCCGGCGCATTTGTGGCGCGTGCGGCGCGTTTTATGTTACATGCGCAGGTTGAGGCAGGGTCGTTATGt                                                                                                                                                                                                             >  1:110112/1‑251 (MQ=255)                                              aaCGGCTGAATCCCTTGAACTGGGGCGGCTGGCGAATGTGAATCCGCCTGAATTATTGCGCTACGATGCGCAAGGACGCCGTCTGGACGATGTGCGTTTTCACCCCGCCTGGCACCTGCTGATGCAGGCGCTATGTACCAATCGGGTGCACAATCTTGCCTGGGAAGAAGACGCTCGCTCCGGCGCATTTGTGGCGCGTGCGGCGCGTTTTATGTTACATGCGCCGGTTGAGGCAGGGTCGTTATGTCCg                                                                                                                                                                                                          >  1:541670/1‑250 (MQ=255)                                                    tGAATCCCTTGAACTGGGGCGGCTGGCGAATGTGAATCCGCCTGAATTATTGCGCTACGATGCGCAAGGACGCCGTCTGGACGATGTGCGTTTTCACCCCGCCTGGCACCTGCTGATGCAGGCGCTATGTACCAATCGGGTGCACAATCTTGCCTGGGAAGAAGACGCTCGCTCCGGCGCATTTGTGGCGCGTGCGGCGCGTTTTATGTTACATGCGCAGGGTGAGGCAGGGTCggttttgccgtatacc                                                                                                                                                                                                    >  2:93969/1‑235 (MQ=255)                                                        tCCCTTGAACTGGGGCGGCTGGCGAATGTGAATCCGCCTGAATTATTGCGCTACGATGCGCAAGGACGCCGTCTGGACGATGTGCGTTTTCACCCCGCCTGGCACCTGCTGATGCAGGCGCTATGTACCAATCGGGTGCACAATCTTGCCTGGGAAGAAGACGCTCGCTCCGGCGCATTTGTGGCGCGTGCGGCGCGTTTTATGTTACATGCGCAGGTTGAGGCAGGGTCg                                                                                                                                                                                                                   <  2:132013/231‑1 (MQ=255)                                                        tCCCTTGAACTGGGGCGGCTGGCGAATGTGAATCCGCCTGAATTATTGCGCTACGATGCGCAAGGACGCCGTCTGGACGATGTGCGTTTTCACCCCGCCTGGCACCTGCTGATGCAGGCGCTATGTACCAATCGGGTGCACAATCTTGCCTGGGAAGAAGACGCTCGCTCCGGCGCATTTGTGGCGCGTGCGGCGCGTTTTATGTTACATGCGCAGGTTGAGGCAGGGTCg                                                                                                                                                                                                                   >  1:132013/1‑231 (MQ=255)                                                            ttgaactggtttggcgggCGAATGTGAATCCGCCTGATTCATTGCGCTACGATGCGCAAGGACGCCGTCTGGCCGATGGGCGTGTTCACCCCGTCTGGCTCCGGCTGATTCAGGCGCTATGTACCATTCGGGTGCACAATCTTGCCTGGGAAGAAGACGCTCGCTCCTGTTCATTTGTGGCGCTTGCGGCGCGTTTTATTTTTCATGCGTAGGTTGAGGCAGGGTCGTTATGTCCGATATCCATGTCCt                                                                                                                                                                                             <  2:369008/233‑1 (MQ=255)                                                             tgaccgtgggCGGCTGGCGAATGTGAATCCGCCTGAATGATTGCGCTTCGCTGCGCAAGGACTCCGTCTGCACGATGGGCGTTTTCACCCCGCCTGGCCCCTGCTGATGCAGTCGCTATGTACCAATCGGGTGCACAATCTTGCCTGGGAAGAAGACGCTCGCTCCGGCGCATTTGTGGCGCTTGCGGCACGTTTTATGTTTCATGCTCAGTTTGAGGCAGGGTCGTTATGTCCGAGATCCATGACCtt‑t                                                                                                                                                                                          <  2:436623/243‑1 (MQ=255)                                                             tgacccggggCGTCTGGCGAAGGGGAATCCGTCTGCATTATTGCGCTCCGATGCGCAAGTACGCCGTCTGGACGATGTGCGTTTTCACCCCGCCTGGCACCGGCGGATGCAGGCGCTATGTACCAAGCGGGTGCACAATCTTGCCTGGGAAGAAGACGCGCGCTCCGTCTCATGTGTGTCGCGTGCGGCGCGTTTTATGTTACATGCGCAGGTTTAGGCTGGGTCGTTATGTCCGATTACCATGACCtt                                                                                                                                                                                            <  2:68061/243‑1 (MQ=255)                                                             tGAACTGGGGCGGCTGGCGAATGTGAATCCGCCTGAATTATTGCGCTACGATGCGCAAGGACGCCGTCTGGACGATGTGCGTTTTCACCCCGCCTGGCACCTGCTGATGCAGGCGCTATGTACCAATCGGGTGCACAATCTTGCCTGGGAAGAAGACGCTCGCTCCGGCGCATTTGTGGCGCGTGCGGCGCGTTTTATGTTACATGCGCAGGTTGAGGCAGGGTCGTTATGTCCGATAACCATGACCTT‑Tg                                                                                                                                                                                         >  1:138466/1‑251 (MQ=255)                                                                  tGGGGCGGCTGGCGAATGTGAATCCGCCTGAATTATTGCGCTACGATGCGCAAGGACGCCGTCTGGACGATGTGCGTTTTCACCCCGCCTGGCACCTGCTGATGCAGGCGCTATGTACCAATCGGGTGCACAATCTTGCCTGGGAAGAAGACGCTCGCTCCGGCGCATTTGTGGCGCGTGCGGCGCGTTTTATGTTACATGCGCAGGTTGAGGCAGGGTCGTTATGTCCGATAACCATTACCTT‑Tg‑ccgcc                                                                                                                                                                                   >  1:321934/1‑251 (MQ=255)                                                                           tGGCGAATGTGCATCCGCCTGAATTATTGCGCTACGATGCGCAAGGACGCCGTCTGGACGATGTGCGTTTTCACCCCGCCTGGCACCTGCTGATGCAGGCGCTATGTCCCAATCGGGTGCAAAATCTTGCCTGGGAAGACGACGATCGCTCCGGCGCATTTGTGGCGCGTTCGGCGCGTTTTATGTTACATGCGCCGGTTGAAGATGGGTCGGTTTGTCCGATAACAATTGCCTT‑TG‑Tcgtcac‑gccactg                                                                                                                                                                         >  2:511711/1‑248 (MQ=255)                                                                            ggCGAATGTGAATCCGCCTGAATTATTGCGCTACGATGCGCAAGGACGCCGTCTGGACGATGTGCGTTTTCACCCCGCCTGGCACCTGCTGATGCAGGCGCTATGTACCAATCGGGTGCACAATCTTGCCTGGGAAGAAGACGCTCGCTCCGGCGCATTTGTGGCGCGTGCGGCGCGTTTTATGTTACATGCGCAGGTTGAGGCAGGGTCGTTATGTCCGATAACCATGACCTT‑TG‑CCGCCAC‑GCCAtt                                                                                                                                                                          <  1:310124/249‑1 (MQ=255)                                                                                     gAATCCGACTGAATTATTGCGCTACGATGCGCAAGGACGCCGTCTGGCCGATGTGCGTTTTCACCCCGCCTGGCACCTGCTGATGCAGGCGCTATGTACCAATCGGGTGCACAATCTTGCCTGGGAAGAAGACGCTCGCTCCGGCGCATTTGTGGCGCGTGCGGCGCGTTTTATGTTACATGCGCAGGTTGAGGGCGGGAGGTTAAGTCAGATAAACATTACATT‑TG‑ccgcccc‑gcccaaggttgatc                                                                                                                                                                  >  2:313035/1‑237 (MQ=255)                                                                                      aaTCCGCCTGAATTATTGCGCTCCGATGCGCAAGGACGCCGTCTGGACGATGTGCGTTTTCCCCCCGCCTGGCACCTGCTGATGCAGGCGCTATGTACCAATCGGGTGCACAATCTTGCCTGGGAAGAAGACGCTCGCTCCGGCGCATTTGTGGCGCGTGCGGCGCGTTTTATGTTACATGCGCAGGTTGAGGCAGGGTCGTTATGTCCGATAACCATGACCTT‑TG‑CCGCCAC‑GCCATTGTTGTT‑A‑CAGAt                                                                                                                                                            <  1:365142/251‑1 (MQ=255)                                                                                           gactttattagttcgcTCCGATGTGTAAGGTCGCAGTCTGGCTGATGTTCTTTTTCCCCCCGCCGGGCACCTGCTGATGTAGGCGCTATGTCCCAATCGGGTGCACAATCTTGCCTGGGAAGAAGACGCTCGCGCCGGCGCATTTGTGGCGCGTGCGGCGCGTTTTATGTTACATGCGCAGGTTGAGGCAGGGTCGTTATGTCCGATAACCATGACCTT‑TG‑CCGCCAC‑GCCATTGTTGTT‑A‑CAGATGTTAc                                                                                                                                                       <  2:19293/238‑1 (MQ=255)                                                                                           gCCTGCATTATTGCGCTACGATGCGCAAGGACGCCGTCTGGACGATGTGCGTTTTCACCCCGCCTGGCACCTGCTGATGCAGGCGCTATGTACCAATCGGGTGCACAATCTTGCCTGGGAAGAAGACGCTCGCTCCGGCGCATTTGTGGCGCGTGCGGCGCGTTTTATGTTACATGCGCAGGTTGAGGCAGGGTCGTTATGTCCGATAAACATGACCTT‑TG‑CCGCCCC‑GCCATTGTTGTT‑A‑CAGATGTTAc                                                                                                                                                       >  1:414132/1‑251 (MQ=255)                                                                                           gCCTGATTTATCGCGCTGCGATGCGCAAGGACCCCTTCTGGACTATGTCCGTTTTCGCCCCGCCTGGCTCCTGCTGATGCAGGCTCTATGTACCAATCGGTTGCTCTATCTTGCCTTGGAATAAGACGCTCGCCCCGGCGCATTTGGGGCGCGTGAGGCGCGTTTTATGTTACATGCGCAGGTTGAGGCAGTGTCGTTATGTCCGATAACCATGACCTT‑TG‑CCGCCAC‑GCCATTGTTGTT‑A‑CAGATGTTAc                                                                                                                                                       <  2:3862/251‑1 (MQ=255)                                                                                               gAATTATTGCGCTACGATGCGCATGGACGCCGTCTGGACGATGTGCGTTTTCACCCCGCCTGGCACCTGCTGATGCAGGCGCTATGTACCAATCGGGTGCACAATCTTGCCTGGGAAGAAGACGCTCGCTCCGGCGCATTTGTGGCGCGTGCGGCGCGTTTTATGTTACATGCGCAGGTTGAGGCAGGGTCGTTATGTCCGATAACCATGACCTT‑TG‑CCGCCAC‑GCCATTGTTGTT‑A‑CAGATGTTACCcg                                                                                                                                                    <  1:550854/250‑1 (MQ=255)                                                                                                aaTTATTGCGCTACGATGCGCAAGGACGCCGTCTGGACGATGTGCGTTTTCACCCCGCCTGGCACCTGCTGATGCAGGCGCTATGTACCAATCGGGTGCACAATCTTGCCTGGGAAGAAGACGCTCGCTCCGGCGCATTTGTGGCGCGTGCGGCGCGTTTTATGTTACATGCGCAGGTTGAGGCAGGGTCGTTATGTCCGATAACCATGCCCTT‑TG‑CCGCCCC‑GCCATTGTTGTT‑A‑CCGATGTTAcccccc                                                                                                                                                  >  2:201683/1‑248 (MQ=255)                                                                                                aaTTATTGCGCTACGATGCGCAAGGACGCCGTCTGGACGATGTGCGTTTTCACCCCGCCTGGCACCTGCTGATGCAGGCGCTATGTACCAATCGGGTGCACAATCTTGCCTGGGAAGAAGACGCTCGCTCCGGCGCATTTGTGGCGCGTGCGGCGCGTTTTATGTTACATGCGCAGGTTGAGGCAGGGTCGTTATGTCCGATAACCATGACCTT‑TG‑CCGCCAC‑GCCATTGTTGTT‑A‑CAGATGTTACCcgcg                                                                                                                                                  <  1:363400/251‑1 (MQ=255)                                                                                                      ggcgcTACGATGCGCAAGGACCCCGTCTGGACGATGTGCGTTTTCGCCCCGCCTGGCACCTGCTGATGCAGGCGCTATGTACCAAGCGGTTGCACAATCTTGCCTGGGAAGAAGACGCTCGCTCCGGCTCATTTGTGGCGCGGGCGGCGCGTTTTAGGTTAGATGCGCAGGTTGAGGTGGGGTCGTTATGTCCGATAGCCATGACCTT‑TG‑CCGCCAC‑GCCATTGTTGTT‑A‑CAGATGTTACCCGCGCCGttt                                                                                                                                            <  1:266863/250‑1 (MQ=255)                                                                                                       gcgcTACTATGCGCAAATACTCAGTCTGGACGATGTGCGTTTTCCCACCTCCGGGCTCCTGCTGATGCAGGCGGTATGTCCCAATGGGGTGCACAATCTTGCCTGGGAAGCAGACGCTCGCTCCGGCGCATTTGTGGCGCGTGCGGCGCGTTTTATGTTACATGCGCAGGTTGAGGCAGAGTCGTTATGTCCGATAACCATGACCTT‑GG‑CCGCCAC‑GCCATTGTTGTT‑A‑CAGATGTTACCCGCGCCGttt                                                                                                                                            <  1:378944/250‑1 (MQ=255)                                                                                                       gcgcTACGATGCGCAAGGACGCCGTCTGGACGATGTGCGTTTTCACCCCGCCTGGCACCTGCTGATGCAGGCGCTATGTACCAATCGGGTGCACAATCTTGCCTGGGAAGAAGACGCTCGCTCCGGCGCATTTGTGGCGCGTGCGGCGCGTTTTATGTTACATGCGCAGGTTGAGGCAGGGTCGTTATGTCCGATAACCATGACCTT‑TG‑CCGCCAC‑GCCATTGTTGTT‑A‑CAGATGTTACCCGCGCCGttt                                                                                                                                            <  1:286920/250‑1 (MQ=255)                                                                                                          cTACGATGCGCAAGGACGCCGTCTGGACGATGTGCGTTTTCACCCCGCCTGGCACCTGCTGATGCAGGCGCTATGTACCAATCGGGTGCACAATCTTGCCTGGGAAGAAGACGCTCGCTCCGGCGCATTTGTGGCGCGTGCGGCGCGTTTTATGTTACATGCGCAGGTTGAGGCAGGGTCGTTATGTCCGATAACCATGACCTT‑TG‑CCGCCAC‑GCCATTGTTGTT‑A‑CAGATGTTACCCGCGCCGTTTCAgg                                                                                                                                        >  1:88416/1‑251 (MQ=255)                                                                                                             cGATGCGCAAGGACGCCGTCTGGACGATGTGCGTTTTCACCCCGCCTGGCACCTGCTGATGCAGGCGCTATGTACCAATCGGGTGCACAATCTTGCCTGGGAAGAAGACGCTCGCTCCGGCGCATTTGTGGCGCGTGCGGCGCGTTTTATGTTACATGCGCAGGTTGAGGCAGGGTCGTTATGTCCGATAACCATGACCTT‑TG‑CCGCCAC‑GCCATTGTTGTT‑A‑CAGATGTTACCCGCGCCGTTTCAGGACt                                                                                                                                     <  1:529116/251‑1 (MQ=255)                                                                                                              gatgcgaaagtgcGCCTTCTGGACGAGGTGCGTTTTTCCCCCGCTGGGCGCCTTCTGATGCAGTCGCTATGTACCAAGCGGGTGCACAATTTTGCGTGGGATTAATAAGCTCGCTCCGGCGCATTTGTGGCGCGTGCGGCGCGTTTTATGTTACGTGCGAAGGTTGAGGCAGGGTCGTTATGTCCGATACCCATG‑CCTT‑TGGCAGCCAC‑GCCATTGTTGTT‑A‑CAGATGTTACCCGCGCCGTTTCAGGACTg                                                                                                                                    <  2:406133/239‑1 (MQ=255)                                                                                                                    cAAGGACGCCGTCTGGACGATGTGCGTTTTCACCCCGCCTGGCACCTGCTGATGCAGGCGCTATGTACCAATCGGGTGCACAATCTTGCCTGGGAAGAAGACGCTCGCTCCGGCGCATTTGTGGCGCGTGCGGCGCGTTTTATGTTACATGCGCAGGTTGAGGCAGGGTCGTTATGTCCGATAACCATGAc                                                                                                                                                                                               <  1:157595/191‑1 (MQ=255)                                                                                                                     atgtacgcccgctagtcGATGTGCGTTCTCACCCCGCCGGGCACCTGCTGATGCAGGCGCAGGGTGCAGAGCGTGTGCATAAATATGCCGGGGAAGAAGACGCTCGTTCTGGCGAAGTTGTGGCTCGTGCGCCGAGTTTTATGTTAGATGGGCAGGTTGAGGAAGGGTCGTTAGGACCAATAACCATGCCCGT‑TG‑CGGCCAC‑GCCATTGTTGTT‑A‑CAGATGTTACCCGCGCCGTTTCAGGACTTGTCCAcg                                                                                                                             <  2:360499/235‑1 (MQ=255)                                                                                                                       ggACGCCGTCTGGACGATGTGCGTTTTCACCCCGCCTGGCACCTGCTGATGCAGGCGCTATGTACAAATCGGGTGCACAATCTTGCCTGGGAAGAAGACGCTCGCTCCGGCGCATTTGAGGCGCGTGCAGCGCGTTTTATGTTACATGCGCAGGTTGAGGTAGGGTCGTTATGGCCGATAACCAAGCACTT‑TG‑CCGCCAC‑GCCATGGTTGTG‑A‑CAGTAGTTACCCGGGCAGGTTCCGGAGTGGAACAcgcc                                                                                                                           >  2:548417/1‑251 (MQ=255)                                                                                                                          cGCCGTCTGGACGATGTGCGTTTTCACCCCGCCTGGCACCTGCTGATGCAGGCGCTATGTACCAATCGGGTGCACAATCTTGCCTGGGAAGAAGACGCTCGCTCCGGCGCATTTGTGGCGCGTGCGGCGAGTTTTATGTTACATGCCCAGGTTGAGGCAGGGTCGTTTTGTCCGATAACCATGACCTT‑TG‑CCGCCAC‑GCCATTGTTGTT‑AACGGA‑GTTACCCGCGCCGTTTCAAGAATGGACACCGCCGCt                                                                                                                        >  2:183912/1‑251 (MQ=255)                                                                                                                              gTCTGGACGATGTGCGTTTTCACCCCGCCTGGCACCTGCTGATGCAGGCGCTATGTACCAATCGGGTGCACAATCTTGCCTGGGAAGAAGACGCTCGCTCCGGCGCATTTGTGGCGCGTGCGGCGCGTTTTATGTTACATGCGCAGGTTGAGGCAGGGTCGTTATGTCCGATAACCATGACCTT‑TG‑CCGCCAC‑GCCATTGTTGTT‑A‑CAGATGTTACCCGCGCCGTTTCAGGACTGGACCACGCCGCTATTg                                                                                                                    >  1:390229/1‑251 (MQ=255)                                                                                                                                 tGGACGATGTGCGTTTTCACCCCGCCTGGCACCTGCTGATGCAGGCGCTATGTACCAATCGGGTGCACAATCTTGCCTGGGAAGAAGACGCTCGCTCCGGCGCATTTGTGGCGCGTGCGGCGCGTTTTATGTTACATGCGCAGGTTGAGGCAGGGTCGTTATGTCCGATAACCATGACCTT‑TG‑CCGCCAC‑GCCATTGTTGTT‑A‑CAGATGTTACCCGCGCCGTTTCAGGGCTGGACCACGCCGCTATTGAGc                                                                                                                 >  1:382608/1‑251 (MQ=255)                                                                                                                                       aTGTGCGTTTTCACCCCGCCTGGCACCTGCTGATGCAGGCGCTATGTACCAATCGGGTGCACAATCTTGCCTGGGAAGAAGACGCTCGCTCCGGCGCATTTGTGGCGCGTGCGGCGCGTTTTATGTTACATGCGCAGGTTGAGGCAGGGTCGTTATGTCCGATAACCATGACCTT‑TG‑CCGCCAC‑GCCATTGTTGTT‑A‑CAGATGTTACCCGCGCCGGTTCAGGGCTGGACCACGCCGCTATTGAGCGATCGc                                                                                                           >  1:225428/1‑251 (MQ=255)                                                                                                                                         gtgCGTTTTCCCCCCGCCTGGCACCTGCTGATGCAGGCGCTATGTACCAATCGGGTGCACAATCTTGCCTGGGAAGAAGACGCTCGCTCCGGCGCATTTGTGGCGCGTGCGGCGCGTTTTATGTTACATGCGCAGGTTGAGGCAGGGTCGTTATGTCCGATAACCATGACCTT‑TG‑CCGCCAC‑GCCATTGTTGTT‑A‑CAGATGTTACCCGCGCCGTTTCAGGACTGGACCACGCCGCTATTGAGCGATCGCt                                                                                                          <  1:201683/250‑1 (MQ=255)                                                                                                                                                    ccccGCCTGGCACCTGCTGATGCAGGCGCTATGTACCAATCGGGTGCACAATCTTGCCTGGGAAGAAGACGCTCGCTCCGGCGCATTTGTGGCGCGTGCGGCGCGTTTTATGTTACATGCGCAGGGTGAGGCAGGGTCGTTATGTCCGATAACCATTACCTT‑TG‑CCGCCAC‑GCCATTGTTGTT‑A‑CAGATGTTACCCGCGCCGTTTCCGGACTGGACCACGCCGCTATTGAGCCATCGCTACGAATCTCACt                                                                                              >  1:21218/1‑251 (MQ=255)                                                                                                                                                           tGGCACCTGCTGATGCAGGCGCTATGTACCAATCGGGTGCACAATCTTGCCTGGGAAGAAGACGCTCGCTCCGGCGCATTTGTGGCGCGTGCGGCGCGTTTTATGTTACATGCGCAGGTTGAGGCAGGGTCGTTATGTCCGATAACCATGACCTT‑TG‑CCGCCAC‑GCCATTGTTGTT‑A‑CAGATGTTACCCGCGCCGGTTCAGGACTGGACCACGCCGCTATTGGGCGATCGCTACGATTCTCACTTATTGcc                                                                                       >  1:210783/1‑251 (MQ=255)                                                                                                                                                               aCCTGCTGATGCAGGCGCTATGTACCAATCGGGTGCACAATCTTGCCTGGGAAGAAGACGCTCGCTCCGGCGCATTTGTGGCGCGTGCGGCGCGTTTTATGTTACATGCGCAGGTTGAGGCAGGGTCGTTATGTCCGATAACCATGACCTT‑TG‑CCGCCAC‑GCCATTGTTGTT‑A‑CAGATGTTACCCGCGCCGTTTCAGGACTGGACCACGCCGCTATTGAGCGATCGCTACGATTCTCACTTATTGCCAggg                                                                                   >  1:114811/1‑250 (MQ=255)                                                                                                                                                                 ctgctgATGCAGGCGCTATGTACCAATCGGGTGCACAATCTTGCCTGGGAAGAAGACGCTCGCTCCGGCGCATTTGTGGCGCGTGCGGCGCGTTTTATGTTACATGCGCAGGTTGAGGCAGGGTCGTTATGTCCGATAACCATGACCTT‑TG‑CCGCCCC‑GCCATTGTTGTT‑A‑CAGATGTTACCCGCGCCGTTTCAGGACTGGACCACGCCGCTATTGAGCGATCGCTACGCTTCTCACTTATTTCCAGGTgg                                                                                 >  2:539318/1‑251 (MQ=255)                                                                                                                                                                 ctgctgATGCAGGCGCTATGTACCAATCGGGTGCACAATCTTGCCTGGGAAGAAGACGCTCGCTCCGGCGCATTTGTGGCGCGTGCGGCGCGTTTTATGTTACATGCGCAGGTTGAGGCAGGGTCGTTATGTCCGATAACCATGACCTT‑TG‑CCGCCAC‑GCCATTGTTGTT‑A‑CAGATGTTACCCGCGCCGTTTCAGGACTGGACCACGCCGCTATTGAGCGATCGCTACGATTCTCACTTATTGCCAGGTgg                                                                                 >  1:120295/1‑251 (MQ=255)                                                                                                                                                                     tgaagcaggggtattGTACGAAGCGGGTGCACAATTTTGCCGGGGGAGAAGACGCGCGCTCTGGCGCATGTGTGGCGCGTGCGGAGCGTTTTATGTTCCAGGCGCAGGTTGAGGCAGGGTCGTTATGTCCGATAACCATG‑CCTTGTG‑CCGCCAC‑GCCATTGTTGTT‑A‑CAGATGTTACCCGTGCCGTTTCAGTATGGGACCACGCCGCTATTGAGCGATCGCTACGATTCTCTCTTATTGCCAGGTGGGCa                                                                              <  2:564919/236‑1 (MQ=255)                                                                                                                                                                           aGGCGCTATGTACCACTCGGGTGCACAATCTTGCCTGGGAAGAAGACACTCGCTCCGGCGCCTTTGTGGCGCGTGCGGCGCGTTTTATGTTCCATGCGCCGGTTGAGGCAGGTGCCTTCTTTCCTCTCACCATTCCCTT‑TG‑CCGCCAC‑GCCCTTGTTGTT‑C‑CACATGTTACCCGCGCCGTTTCAGGACTTGCCCCCGCCGCATTTGGGCCATCGCTCCGATTCTACCTTCTTGCCaccatcggccgcgccc                                                                       >  2:38379/1‑236 (MQ=255)                                                                                                                                                                           aGGCGCTATGTACCAATCGGGTGCACAATCTGGCCTGGGAAGAAGACGCTCGCTCCGGCGCATTTGTGGCGCGTGCGGCGCGTTTTATGTTACATGCGCAGGTGGAGGCAGGGTCGTTATGTCCGATAACCATGACCTG‑TG‑CCGCCAC‑GCCATTGTTGTT‑A‑CAGATGTTACCCGCGCCGTTTCAGGACTGGACCACGCCGCTATTGAGCGATCGCTACGATTCTCACTTATTGCCa                                                                                      <  1:38379/236‑1 (MQ=255)                                                                                                                                                                              cgcTATGTACCAATCGGGTGCACAATCTTGCCTGGGAAGAAGACGCTCGCTCCGGCGCATTTGTGGCGCGTGCGGCGCGTTTTATGTTACATGCGCAGGTTGAGGCAGGGTCGTTATGTCCGATAACCATGACCTT‑TG‑CCGCCCC‑GCCATTGTTGTT‑A‑CAGGTGTTACCCGCGCCGTTTCAGGACTGGACCACGCCGCTATTGAGCGATCGCTACGCTTCTCACTTAATGCCAGGTGGGCaaaaaacgggg                                                                    >  2:551600/1‑245 (MQ=255)                                                                                                                                                                                cTATGTACCAATCGGGTGCACAATCTTGCCTGGGAAGAAGACGCTCGCTCCGGCGCATTTGTGGCGCGTGCGGCGCGTTTTATGTTACATGCGCAGGTTGAGGCAGGGTCGTTATGTCCGATAACCATGACCTT‑TG‑CCGCCAC‑GCCAttg                                                                                                                                                                         <  2:468117/150‑1 (MQ=255)                                                                                                                                                                                cTATGTACCAATCGGGTGCACAATCTTGCCTGGGAAGAAGACGCTCGCTCCGGCGCATTTGTGGCGCGTGCGGCGCGTTTTATGTTACATGCGCAGGTTGAGGCAGGGTCGTTATGTCCGATAACCATGACCTT‑TG‑CCGCCAC‑GCCAttg                                                                                                                                                                         >  1:468117/1‑150 (MQ=255)                                                                                                                                                                                 tATGTTCCAATCGGGTGCACAATCTTGCCTGGGAAGAAGACGCTCGCTCGGGCGCATTTGTGGCGCGTGCGGCGCGTTTTATGTTACAGGCGCAGGTTGAGGCAGGGTCGTTATGTCCGATAACCATTACCTT‑TG‑CCGCCAC‑GCCATTGTTGTT‑A‑CAGATTTTACCCGCGCCGTTTCAGGACTGGACCACGCCGCTATTGAGCGATCGCTACGATTCTCACTTATTGCCAGGTGGGCAAAAACGCGGttt                                                                  <  1:537930/250‑1 (MQ=255)                                                                                                                                                                                  atgtcccagacGGGTGCACAATTTTGCCGGGGAAGAAGCCGCTCGCCCCGGCGCGTTTGTGGCGCGTGCGGCGCGTTTTTTTTTCCATGCGGAGGGTGAGGCAGGGTCGTTATGTCCTATACCCATGACTGT‑GG‑CCGCCAC‑GCCATTGTTGTT‑A‑CAGATGTTACCCGCGCCGTTTCAGGACTGGACCACGCCGCTATTGAGCGATCGCTACGATTCTCACTTATTGCCAGGTGGGCAAAAACGCGGTttgt                                                                <  2:382608/241‑1 (MQ=255)                                                                                                                                                                                   tGTACCAATCGGGTGCACAATCTTGCCTGGGAAGAAGACGCTCGCTCCGGCGCATTTGTGGCGCGTGCGGCGCGTTTTATGTTACATGCGCAGGTTGAGGCAGGGTCGTTATGTCCGGTAACCATGACCTT‑TG‑CCGC‑ACAGCCATTGTTGTT‑A‑CAGATGTTACCCGCGCCGTTTCAGGACTGGACCACGCCGCTATTGAGCGTTCGCTACGATTCTCACTTATTGCCAGGGGGGCAAAAACGCGGTttgtt                                                               >  2:289895/1‑251 (MQ=255)                                                                                                                                                                                     tACCAATCGGGTGCACAATCTTGCCTGGGAAGAAGACGCTCGCTCCGGCGCATTTGTGGCGCGTGCGGCGCGTTTTATGTTACATGCGCAGGTTGAGGCAGGGTCGTTATGTCCGCTAACCATTACCTT‑TG‑CCGCCAC‑GCCCTTGGTGTT‑A‑CAGATGTTACCCGCGCCGTTTCAGGACTGGACCACGCCGCTATTGAGCGATCGCTACGAATCTCCCTTATTGCCAGGAGGGTAAACACGCGGTttgttgg                                                             >  2:315749/1‑250 (MQ=255)                                                                                                                                                                                        cAATCGGGTGCACAATCTTGCCTGGGAAGAAGACGCTCGCTCCGGCGCATTTGTGGCGCGTGCGGCGCGTTTTATGTTACATGCGCAGGTTGAGGCAGGGTCGTTATGTCCGATAACCATGACCTT‑TG‑CCGCCAC‑GCCATTGTTGTT‑A‑CAGATGTTACCCGCGCCGTTTCAGGACTGGACCACGCCGCTATTGAGCGATCGCTACGATTCTCACTTATTGCCAGGTGGGCAAAAACGCGGTTTGTTGATTg                                                          >  1:219381/1‑251 (MQ=255)                                                                                                                                                                                         aaTCGGGTGCACCATCTTGCCTGGGAAGAAGACGCTCGCTCCGGCGCATTTGTGGCGCGTGCGGCGCGTTTTCTGTTACCTGCGCAGGGTGAGGCAGGGTCGTTATGTCCGATAACCATGGCCTT‑TG‑CCGCCCC‑GCCCTTGTTGTTCA‑CA‑ATGTTACCCGCGCCGTTTCAGGACTCGACCAAGCCGCTATTGAGCGATCGCTAAGATTCGCACTTCATGCCAGGGGGGCCATACCGCGGgttgttcagtgg                                                         >  2:267429/1‑245 (MQ=255)                                                                                                                                                                                          aTCGGGTGCACAATCTTGCCTGGGAAGAAGACGCTCGCTCCGGCGCATTGGTGGCGCGTGCGGCGCGTTTTATGTTACATGCGCAGGTTGAGGCAGGGTCGTTATGTCCGATAACCATGACCTT‑TG‑CCGCCAC‑GCCATTGTTGTT‑A‑CAGATGTTACCCGCGCCGTTTCAGGACTGGACCACGCCGCTATTGAGCGATCGCTGCGATGCTCACTTATTGTCAGGTGGGCAAAAACGCGGTTTGTTGATTGGc                                                        <  1:483347/251‑1 (MQ=255)                                                                                                                                                                                            cGGGTGCACAATCTTGCCTGGGAAGAAGACGCTCGCTCCGGCGCATTTGTGGCGCGTGCGGCGCGTTTTATGTTACATGCGCAGGTTGAGGCAGGGTCGTTATGTCCGATAACCATGACCTT‑TG‑CCGCCAC‑GCCATTGTTGTT‑A‑CAGATGTTACCCGCGCCGTTTCAGGACTGGACCACGCCGCTATTGAGCGATCGCTACGATTCTCACTTATTGCCAGGTGGGCAAAAACGCGGTTTGTTGATTGGCAt                                                      >  1:496149/1‑251 (MQ=255)                                                                                                                                                                                             gggTGCACAAGCTTGCCGGGGATGAAGCCGCGCGCTCCGGCGCTTTGGTGGCGTGTGCGGCGCTTTTTATGTTACAGTCGCAGGGTGAGGCAGGTTCGTTATGTCCGATAACCTTGGCCTT‑TG‑CCGCCAC‑GCCATTGTTGTT‑A‑CAGATGTTACCCGCGCCGTTTCAGGACTGGACCACGCCGCTATTGAGCGATCGCTACGATTCTCACTTATTGCCAGGTGGGCAAAAACGCGGTTTTTTGATTGGCATg                                                     <  2:257270/251‑1 (MQ=255)                                                                                                                                                                                                          tGCCTGGGAAGAAGACGCTCGCTCCGGCGCATTTGTGGCGCGTGCGGCGCGTTTTATGTTACATGCGCAGGTTGAGGCAGGGTCGTTATGTCCGATAACCATGACCTT‑TG‑CCGCCAC‑GCCATTGTTGTT‑A‑CAGATGTTACCCGCGCCGTTTCAGGACTGGACCACGCCGCTATTGAGCGATCGCTACGATTCTCACTTATTGCCAGGTGGGCCAAAACGCGGTTTGTTGATTGGCATGGGAATGACGGaa                                         >  1:248887/1‑250 (MQ=255)                                                                                                                                                                                                               gggAAGAAGACGCTCGCTCCGGCGCATTTGTGGCGCGTGCGGCGCGTTTTATGTTACATGCGCAGGTTGAGGCAGGGTCGTTATGTCCGATAACCATGACCTT‑TG‑CCGCCAC‑GCCATTGTTGTT‑A‑CAGATGTTACCCGCGCCGTTTCAGGACTGGACCACGCCGCTATTGAGCGATCGCTACGATTCTCACTTATTGCCAGGTGGGCAAAAACGCGGTTTGTTGATTGGCATGGGAATGACGGAAAAGCAg                                   <  1:93969/251‑1 (MQ=255)                                                                                                                                                                                                                  aagaagcacctcgctcCGGCGACTTTGTGGCGCGTGCGGCGCGTTTTATGTTTCATGCGCAGGTTGAGGCAGGGTCTTTATGTCCGATACCCATGACCTT‑TG‑CCGCCAC‑GCCATTGTTGTT‑A‑CAGATGTTACCCGCGCCGTTTCAGGACTGGACCACGCCGCTATTGAGCGATCGCTACGATTCTCACTTATTGCCAGGTGGGCAAAAACGCGGTTTGTTGATTGGCATGGGAATGACGGAAAAGCAGGGc                                <  2:496149/242‑1 (MQ=255)                                                                                                                                                                                                                      agaCGCTCGCTCCGGCGCATTTGTGGCGCGTGCGTCGCGTTTTATGTTACATGCGCAGGTTGAGGCAGGGTCGTTATGTCCGATAACCATGACCTT‑TG‑CCGCCAC‑GCCATTGTTGTT‑A‑CAGATGTTACCCGCGCCGTTTCAGGACTGGCCCCAGCCCCTATTGAGCGATCGCTACGCTTATCACTTATTTCCAGGGGGGCCAAAACTCGGTTTGTTGTTTGGCATGGGTATGCGGGAAAAGCAGGGGGGtt                            >  2:223680/1‑251 (MQ=255)                                                                                                                                                                                                                      agaCGCTCGCTCCGGCGCATTTGTGGCGCGTGCGGCGCGTTTTATGTTACATGCGCAGGTTGAGGCAGGGTTGTTATGTCCGATAACCATGACCTT‑TG‑CCGCCAC‑GCCATTGTTGTT‑A‑CAGATGTTACCCGCGCCGTTTCAGGACTGGACCACGCCGCTATTGAGCGATCGCTACGATTCTCACTTATTGCCAGGTGGGCAAAAACGCGGTTTGTTGATTGGCATGGGAATGACGGAAAAGCAGGGCGGtt                            >  1:281023/1‑251 (MQ=255)                                                                                                                                                                                                                              gctcCGGCGCATTTGTGGCGCGTGCGGCGCGTTTTATGTTACATGCGCAGGTTGAGGCAGGGTCGTTATGTCCGATAACCATGACCTT‑TG‑CCGCCAC‑GCCATTGTTGTT‑A‑CAGATGTTACCCGCGCCg                                                                                                                                               >  1:548691/1‑128 (MQ=255)                                                                                                                                                                                                                              gctcCGGCGCATTTGTGGCGCGGGCGGCGAGTTTTATGTTACATGCGCAGGTTGAGGAAGGGTCGTTATGTCCGATAACCATGACCTT‑TG‑CCGCCAC‑GCCATTGTTGTT‑A‑CAGATGTTACCCGCGCCg                                                                                                                                               <  2:548691/128‑1 (MQ=255)                                                                                                                                                                                                                                                cgcgtgcgGCGCGTTTTATGTTACATGCGCAGGTTGAGGCAGGGTCGTTATGTCCGATAACCATGACCTT‑TG‑CCGCCAC‑GCCATTGTTGTT‑A‑CAGATGTTACCCGCGCCGTTTCAGGACTGGACCACGCCGCTATTGAGCGATCGCTACGATTCTCACTTATTGCCAGGTGGGCAAAAACGCGGTTTGTTGATTGGCATGGGAATGACGGAAAAGCAGGGCGGTTCCGATGTTATGAGCAACACCACCCgg  >  1:225945/1‑250 (MQ=255)                                                                                                                                                                                                                                                    |                                                                                                                                                                                                                                                            TGGCTGGGATAGCGATTTTCTAGCCAGTATTGGTCAGCAGTTAGGAACGGCTGAATCCCTTGAACTGGGGCGGCTGGCGAATGTGAATCCGCCTGAATTATTGCGCTACGATGCGCAAGGACGCCGTCTGGACGATGTGCGTTTTCACCCCGCCTGGCACCTGCTGATGCAGGCGCTATGTACCAATCGGGTGCACAATCTTGCCTGGGAAGAAGACGCTCGCTCCGGCGCATTTGTGGCGCGCGCGGCGCGTTTTATGTTACATGCGCAGGTTGAGGCAGGGTCGTTATGTCCGATAACCATGACCTT‑TG‑CCGCCAC‑GCCATTGTTGTT‑A‑CAGATGTTACCCGCGCCGTTTCAGGACTGGACCACGCCGCTATTGAGCGATCGCTACGATTCTCACTTATTGCCAGGTGGGCAAAAACGCGGTTTGTTGATTGGCATGGGAATGACGGAAAAGCAGGGCGGTTCCGATGTTATGAGCAACACCACCCGT  >  NC_000913/4414382‑4414871

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