Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	NC_000913	4,297,469	C→T	G633D (GGC→GAC)	fdhF ←	formate dehydrogenase H

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	4,297,469	0	C	T	98.1%	160.3 / ‑6.3	54	G633D (GGC→GAC)	fdhF	formate dehydrogenase H
Reads supporting (aligned to +/- strand):  ref base C (0/1);  new base T (31/21);  total (31/23)
Fisher's exact test for biased strand distribution p-value = 4.15e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.01e-01

GCCAGTGCCGCTTCGCGCAGGCGAGTTTTCAACTTGTTGTACTCGTCAATCACGTACTGCTCGGCGGCGCGCTGATCGGCGATCGGCTCGACGCGAACGGCGCAGTATTTGTACTCCGGCGTTTTCGTAATCGGGCTTAAGTTTTCGGTAACCAGCTCGTTACAGGCACCAATCCACCACTGGTAGGTCATGTAAATCGCCCCTTTGTTCGGACGATCGCTGACCTGCGCACGGGTGATAATTTTGCCTTTACGCGAGTGCACCCAAACCAATGCCTCATCTTCAATACCCAGACGTTTGGCGTCTTCGGTATTGATTTGTGCGTAGCCAGGTTCATCAGCCAGCGCCGCCAGTGCCGCACAGTTACCGGTCATCGAACGGCAAGAGTAGTGACCAACTTCACGCACCGTTGACAGTACCATCGGGTACTCGTCGGTGAGTTTGTCGATTGGCGCTACCCAGTCGCAGGTGAAGAACTGCGCCAGACCGTTC  >  NC_000913/4297223‑4297714                                                                                                                                                                                                                                                       |                                                                                                                                                                                                                                                      gCCAGTGCCGCATCGCGAAGGCGCGTTATCAACTTGTTGTACTCGTCAATCACGTACTGCTCGGCGGCGCGATGCTAGGCGATCGGCTCGACGCGAACGGCGCCGTATTTGTACTCCGGCGTTTTCGTACCCGGGCTTAAGTTTTCGGTAACCAGCTCGTTACAGGCACCAATCCACCACTGGTAGGTCATGTAAATCGCCCCTTTGTTCGGACGATCGCTGACCTGCGCACGGGTGATAATTTTGTCttt                                                                                                                                                                                                                                                   >  2:313495/1‑251 (MQ=255)  ccAGTGCCGCTTCGCGCAGGCGAGTTTTCAACTTGTTGTACTCGTCAATCACGTACTGCTCGGCGGCGCGCTGATCGGCGATCGGCTCGACGCGAACGGCGCAGTCTTTGTACTCCGGCGTTTTCGTAATCGGGCTTAAGTTTTCGGTAACCAGCTCGTTACAGGCACCAATTCACCACTGGTAGGTCATGTAAATCGCCCCTTTGTTCGGACGATCGCTGACCTGCGCACGGGTGATAATTTTGTCttt                                                                                                                                                                                                                                                   >  2:139930/1‑250 (MQ=255)  ccAGTGCCGCTTCGCGCAGGCGAGTTTTCAACTTGTTGTACTCGTAAATCACGTACTGCTCGGCGGCGCGCTGCTCGGCGATCGGCTCGACGCGAACGGCGCAGTCTTTGTACTCCGGCGTTTTCGTAATCGGGCTTAAGTTTTCGGTAACCAGCTCGTTACAGGCACCAATCCACCACTGGTAGGTCATGTAAATCGCCCCTTTGTTCGGACGATCGCTGACCTGCGCACGGGTGATAATTTTGTCttt                                                                                                                                                                                                                                                   >  1:277213/1‑250 (MQ=255)          gCTTCGCGCAGGCGAGTTTTCAACTTGTTGTACTCGTCAATCACGTACTGCTCGGCGGCGCGCTGATCGGCGATCGGCTCGACGCGAACGGCGCAGTATTTGTACTCCGGCGTTTTCGTAATCGGGCTTAAGTTTTCGGTAACCAGCTCGTTACAGGCACCAATCCACCACTGGTAGGTCATGTAAATCGCCCCTTTGTTCGGACGATCGCTGACCTGCGCACGGGTGATAATTTTGTCTTTACGCGAGTg                                                                                                                                                                                                                                          >  2:309949/1‑251 (MQ=255)                 gcAGTCGAGTTTTCAACTTGTTGTACTCGTCAATCACGTACTGCTCGGCGGCGCGCTGATCGGCGATCGGCTCGACGCGAACGGCGCAGGATTTGTACTCCGGCGTTTTCGTAATCGGGCTTAAGTTTTCGGTAACCAGCTCGTTACAGGCACCAATCCGCCACTGGTAGGTCATGTAAATCGCCCCTTTGTTCGGACGATCGCGGACCTGCTCACGGGTGATAATTTTGTCTTTACGCGAGTGCACCCaa                                                                                                                                                                                                                                   <  2:411947/251‑1 (MQ=255)                  cAGGCGAGTTTTCAACTTGTTGTACTCGTCAATCACGTACTGCTCGGCGGCGCGCTGATCGGCGATCGGCTCGACGCGAACGGCGCAGCATTTGTACTCCGGCGTTTTCGTAATCGGGATTAAGTTTTCGGTAACCAGCTCGTTCCAGGCACCAATCCACCACTGGTAGGTCATGTAAATCGCCCCTTTGTTCGGACGATCGCTGACCTGCGCACGGGTGATAATTTTGTCTTTACGCGAGTGCACCCaa                                                                                                                                                                                                                                   >  1:209468/1‑250 (MQ=255)                         gTTTTCAACTTGTTGTACTCGTCAATCACGTACTGCTCGGCGGCGCGCTGATCGGCGATCGGCTCGACGCGAACGGCGCAGTATTTGTACTCCGGCGTTTTCGTAATCGGGCTTAAGTTTTCGGTAACCAGCTCGTTACAGGCACCAATCCACCACTGGTAGGTCATGTAAATCGCCCCTTTGTTCGGACGATCGCTGACCTGCGCACGGGTGATAATTTTGTCTTTACGCGAGTGCACCCAAACCAATGc                                                                                                                                                                                                                           >  2:331651/1‑251 (MQ=255)                         gTTTTCAACTTGTTGTACTCGTCAATCACGTACTGCTCGGCGGCGCGCTGATCGGCGATCGGCTCGACGCGAACGGCGCAGTATTTGTACTCCGGCGTTTTCGTAATCGGGCTTAAGTTTTCGGTAACCAGCTCGTTACAGGCACCAATCCACCACTGGTAGGTCATGTAAATCGCCCCTTTGTTCGGACGATCGCTGACCTGCGCACGGGTGATAATTTTGTCTTTACGCGAGTGCACCCAAACCAATGc                                                                                                                                                                                                                           <  1:436806/251‑1 (MQ=255)                            ttCAACTTGTTGTACTCGTCAATCACGTACTGCTCGGCGGCGCGCTGATCGGCGATCGGCTCGACGCGAACGGCGCAGTATTTGTACTCCGGCGTTTTCGTAATCGGGCTTAAGTTTTCGGTAACCAGCTCGTTACAGGCACCAATCCACCACTGGTAGGTCATGTAAATCGCCCCTTTGTTCGGACGATCGCTGACCTGCGCACGGGTGATAATTTTGTCTTTACGCGAGTGCACCCAAACCAATGCCTc                                                                                                                                                                                                                        >  2:110444/1‑251 (MQ=255)                               aaCTTGTTGTACTCGTCAATCACGTACTGCTCGGCGGCGCGCTGATCGGCGATCGGCTCGACGCGAACGGCGCAGTATTTGTACTCCGGCGTTTTCGTAATCGGGCTTAATTTTTCGGTAACCAGCTCGTTACAGGCACCAATCCACCACTGGTAGGTCATGTAAATCGCCCCTTTGTTCGGACGATCGCTGACCTGCGCACGGGTGATAATTTTGTCTTTACGCGAGTGCACCCAAACCAATGCCTCATc                                                                                                                                                                                                                     <  1:228560/251‑1 (MQ=255)                                   tgttgtACTCGTCAATCACGTACTGCTCGGCGGCGCGCTGATCGGCGATCGGCTCGACGCGAACGGCGCAGTATTTGTACTCCGGCGTTTTCGTAATCGGGCTTAAGTTTTCGGTAACCAGCTCGTTACAGGCACCAATCCACCACTGGTAGGTCATGTAAATCGCCCCTTTGTTCGGACGATCGCTGACCTGCGCACGGGTGATAATTTTGTCTTTACGCGAGTGCACCCAAACCAATGCCTCATCTTc                                                                                                                                                                                                                  >  2:154343/1‑250 (MQ=255)                                                    aCGTACTGCTCGGCGGCGCGCTGATCGGCGATCGGCTCGACGCGAACGGCGCAGTATTTGTACTCCGTCGTTTTCGTAATCGGGCTTAAGTTTTCGGTAACCAGCTCGTTACAGGCACCAATCCACCACTGGTAGGTCATGTAAATCGCCCCTTTGTTCGGACGATCGCTGACCTGCGCACGGGTGATAATTTTGTCTTTACGCGAGTGCACCCAAACCAATGCCTCATCTTCAATACCCAGACGTTTGGc                                                                                                                                                                                                <  2:20635/251‑1 (MQ=255)                                                                 cggcgCGCTGATCGGCGATCGGCTCGACGCGAACGGCGCAGTATTTGTACTCCGGCGTTTTCGTAATCGGGCTTAAGTTTTCGGTAACCAGCTCGTTACAGGCACCAATCCACCACTGGTAGGTCATGTAAATCGCCCCTTTGTTCGGACGATCGCTGACCTGCGCACGGGTGATAATTTTGTCTTTACGCGAGTGCACCCAAACCAATGCCTCATCTTCAATACCCAGACGTTTGGCGTCTTCGGTATTg                                                                                                                                                                                   <  1:351780/251‑1 (MQ=255)                                                                    cgcgcTGATCGGCGATCGGCTCGACGCGAACGGCGCAGTATTTGTACTCCGGCGTTTTCGTAATCGGGCTTAAGTTTTCGGTAACCAGCTCGTTACAGGCACCAATCCACCACTGGTAGGTCATGTAAATCGCCCCTTTGTTCGGACGATCGCTGACCTGCGCACGGGTGATAATTTTGTCTTTACGCGAGTGCACCCAAACCAATGCCTCATCTTCAATACCCAGACGTTTGGCGTCTTCGGTATTGAtt                                                                                                                                                                                >  2:370279/1‑251 (MQ=255)                                                                     gcgcTGATCGGCGATCGGCTCGACGCGAACCGCGCAGGATTTGTACTCCGGCGTTTTCGTAATCGGGCTTAAGTTTTCGGTAACCAGCTCGTTACAGGCACCAATCCACCACTGGTAGGTCATGTAAATCGCCCCTTTGTTCGGACGATCGCTGACCTGCGCACGGGTGATAATTTTGTCTTTACGCGAGTGCACCCAAACCAATGCCTCATCTTCAATACCCAGACGTTTGGCGTCTTCGGTATTGAttt                                                                                                                                                                               >  2:122160/1‑251 (MQ=255)                                                                         tGATCGGCGATCGGCTCGACGCGAACGGCGCAGTATTTGTACTCCGGCGTTTTCGTAATCGGGCTTAAGTTTTCGGTAACCAGCTCGTTACAGGCACCAATCCACCACTGGTAGGTCCTGTAAATCGCCCCTTTGTTCGGACGATCGCTGACCTGCGCACGGGTGATAATTTTGTCTTTACGCGCGTGCACCCAAACCAATGCCTCATCTTCAATAACCAGACGTTTGGCGTCTTCGGTATTGATTTGTGc                                                                                                                                                                           >  1:464372/1‑251 (MQ=255)                                                                                 gatcggcTCGCCGCGAACGGCGCAGTATTTGTACTCCGGCGTTTTCGTAATCGGGCTTAAGTTTTCGGTAACCAGCTCGTTACAGGCACCAATCCACCACTGGTAGGTCATGTAAATCGCCCCTTTGTTCGGACGATCGCTGACCTGCGCACGGGTGATAATTTTGTCTTTACGCGAGTGCACCCCAACCAATGCCTCATCTTCAATACCCAGACGTTTGGCGTCTTCGGTATTGATTTGTGCGTAGcc                                                                                                                                                                     >  2:374310/1‑249 (MQ=255)                                                                                 gatcggcTCGACGCGAACGGCGCAGTATTTGTAATCCGGCGTTTTCGTAATCGGGCTTAAGTTTTCGGTAACCAGCTCGTTACAGGCACCAATCCACCACTGGTAGGTCATGTAAATCGCCCCTTTGTTCGGACGATCGCTGACCTGCGCACGGGTGATAATTTTGTCTTTACGCGAGTGCACCCAAACCAATGCCTCATCTTCAATACCCAGACGTTTGGCGTCTTCGGTATTGATTTGTGCGTAGcc                                                                                                                                                                     >  1:143861/1‑249 (MQ=255)                                                                                   tcggcTCGACGCGAACGGCGCAGTATTTGTACTCCGGCGTTTTCGTAATCGGGCTTAAGTTTTCGGTAACCAGCTCGTTACAGGCACCAATCCACCACTGGTAGGTCATGTAAATCGCCCCTTTGTTCGGACGATCGCTGACCTGCGCACGGGTGATAATTTTGTCTTTACGCGAGTGCACCCAAACCAATGCCTCATCTTCAATACCCAGACGTTTGGCGTCTTCGGTATTGATTTGTGCGTAGCCAGGt                                                                                                                                                                 >  1:449279/1‑251 (MQ=255)                                                                                         cGACGCGAACGGCGCAGTATTTGTACTCCGGCGTTTTCGTAATCGGGCTTAAGTTTTCGGTAACCAGCTCGTTACAGGCACCAATCCACCACTGGTAGGTCATGTAAATCGCCCCTTTGTTCGGACGATCGCTGACCTGCGCACGGGTGATAATTTTGTCTTTACGCGAGTGCACCCAAACCAATGCCTCATCTTCAATACCCAGACGTTTGGCGTCTTCGGTATTGATTTGTGCGTAGCCAGGTtcatca                                                                                                                                                           <  1:388538/251‑1 (MQ=255)                                                                                                   ggCGCAGTATTTGTACTCCGGCGTTTTCGTAATCGGGCTTAAGTTTTCGGTAACCAGCTCGTTACAGGCACCAATCCACCACTGGTAGGTCATGTAAATCGCCCCTTTGTTCGGACGATCGCTGACCTGCGCACGGGTGATAATTTTGTCTTTACGCGAGTGCACCCAAACCAATGCCTCATCTTCAATACCCAGACGTTTGGCGTCTTCGGTATTGATTTGTGCGTAGCCAGGTTCATCCGCCAGcgccg                                                                                                                                                 >  2:6040/1‑251 (MQ=255)                                                                                                       cAGTATTTGTACTCCGGCGTTTTCGTAATCGGGCTTAAGTTTTCGGTAACCAGCTCGTTACAGGCACCAATCCACCACTGGTAGGTCATGTAAATCGCCCCTTTGTTCGGACGATCGCTGACCTGCGCACGGGTGATAATTTTGTCTTTACGCGAGTGCACCCAAACCAATGCCTCATCTTCAATACCCAGACGTTTGGCGTCTTCGGTATTGATTTGTGCGTAGCCAGGTTCATCAGCCAGCGCCGCCAg                                                                                                                                             <  1:154343/251‑1 (MQ=255)                                                                                                        aGTATTTGTACTCCGGCGTTTTCGTAATCGGGCTTAAGTTTTCGGTAACCAGCTCGTTACAGGCACCAATCCACCACTGGTAGGTCATGTAAATCGCCCCTTTGTTCGGACGATCGCTGACCTGCGCACGGGTGATAATTTTGTCTTTACGCGAGTGCACCCAAACCAATGCCTCATCTTCAATACCCAGACGTTTGGCGTCTTCGGTATTGATTTGTGCGTAGCCAGGTTCATCAGCCAGCGCCGCCAGt                                                                                                                                            >  1:129870/1‑251 (MQ=255)                                                                                                                   cccGGCGTGTTCGTAATCGGTCTTAAGTTGTCGGTAACCCGCTCGTTACAGGCACCAATCCACCACTGGGAGGTCATGTAAATCGCCCCTTTGTTCGGACGATCGCTGACCGGCGCACGGGTGATAATTTTGCCTTTACGCGAGTGCACCCAAACCAAGGCCTCATCTTCAATACCCAGACGTTGTGCATCTGCGGTTTTGATTGGTGCATAGCCAGGTGCATCAGCCAGCGCCGCCAGTGCCTCACAGtt                                                                                                                                 <  2:303076/250‑1 (MQ=255)                                                                                                                    ccGGCGTTTTCGTAATCGGGCTTAAGTTTTCGGTAACCAGCTCGTTACAGGCACCAATCCCCCACTGGTAGGTCATGTAAATCGCCCCTTTGTTCGGACGATCGCTGACCTGCGCACGGGTGATAATTTTGTCTTTACGCGAGTGCACCCAAACCAATGCCTCATCTTCAATACCCAGACGTTTGGCGTCTTCGGTATTGATTTGTGCGTAGCCAGGTTCATCAGCCAGCGCCGCCAGTGCCGCACAGtt                                                                                                                                 <  1:120480/250‑1 (MQ=255)                                                                                                                           tttCGTAATCGGGCTTAAGTTTGCGGTAACCAGCTCGTTACAGGCACCAATCCCCCACTGGTAGGTCATGTAAATCGCCCCTTTGTTCGGACGATCGCTGCCCTGCGCACGGGTGATAATTTTGTCTTTACGCGAGTGCACCCAAACCAATGCCTCATCTTCAATACCCAGACGTTTGTCGTCTTCGGTATTGATTTGGGCGAAGCCAGGTTCATCAGCCAGGGCCGCCAGTGCCGCACAGGTACCGGTc                                                                                                                          <  2:43433/250‑1 (MQ=255)                                                                                                                              cGTAATCGGGCTTAAGTTTGCGGTAACCAGCTCGTTACAGGCACCAATCCACCACTGGTAGGTCAGGGAAATCGCCCCTTGGTTCGGACGATCGCTGACCGGCGCACGGGTGATAATTTTGTCTTTACGCGAGTGCACCCAAACCAAGGCCCGATCTTCAATGCCCAGCCGTTGGGCGTCTTCGGTAGTGATTGGTGCGTAGCCAGGTTAAGCAGCCAGCGCTGCCAGTGCCTCACAGTTACCTGTCATCg                                                                                                                      <  2:171019/251‑1 (MQ=255)                                                                                                                               gtagtCGGGCTAAAGTTTTCGGTAACCAGCTCGTTACAGACACCAATCCACCACTGGTAGGTCATGTAAATCGCCCCTTTGTTCGGACGATCGCTGACCTGCGCACGGGTGATAATTTTGTCTTTACGCGCGTGCACCCAAACCAATGCCTCATCTTCAATACCCAGACGTTTGGCGTCTTCGGTATTGATTTGTGCGTAGCCAGGTTCATCAGCCAGCGCCGCCAGTGCCGCACAGTTACCGGTCCTCg                                                                                                                      >  2:40712/5‑250 (MQ=255)                                                                                                                                               tttCGGTAACCAGCTCGTTACAGCCACCAATCCACCACTGGTAGGTCATGTAAATCGCCCCTTTGTTCGGACGATCGCTGACCTGCGCACGGGTGATAATTTTGTCTTTACGCGAGTGCACCCAAACCAATGCCTCATCTTCAATACCCAGACGTTTGGCGTCTTCGGTATTGATTTGTGCGTAGCCAGGTTCATCAGCCAGCGCCGCCAGTGCCGCACATTTACCGGTCATCGAACGTCAAGAGTAGTg                                                                                                      <  2:298843/250‑1 (MQ=255)                                                                                                                                               tttCGGAAACCAGCTCGTTACAGGCACCAATCCACCACTGATAGGTCATGTAAATCGCCCCTTTGTTCGGACGATCGCTGACCTGCGCACGGGTGATAATTTTGTCTTTACGCGAGTGCACCCAAACCAATGCCTCATCTTCAATACCCAGACGTTTGGCGTCTTCGGTATTGATTTGTGCGTAGCAAGGTTCATCAGCCAGCGCCGCCAGTGCCGCACAGTTACCGGTCATCGAACGGCAAGAGTCGTg                                                                                                      >  2:108872/1‑250 (MQ=255)                                                                                                                                                    gTAACCAGCTCGTTACAGGCACCAATCCACCACTGGTAGGTCATGTAAATCGCCCCTTTGTTCGGACGATCGCTGACCTGCGCACGGGTGATAATTTTGTCTTTACGCGAGTGCACCCAAACCAATGCCTCATCTTCAATACCCAGACGTTTGGCGTCTTCGGTATTGATTTGTGCGTAGCCAGGTTCATCAGCCAGCGCCGCCAGTGCCGCACAGTTACCGGTCATCGAACGGCAAGAGTAGTGACCAAc                                                                                                >  1:84072/1‑251 (MQ=255)                                                                                                                                                      aaCCAGCTCGTTACAGGCACCAATCCACCACTGGTAGGTCATGTAAACCGCCCCTTTGTTCGGACGATCGCTGACCTGCGCCCGGGTGATAATTTTGTCTTTACGCGAGTGCACCCAAACCAAGGCCTCATCTTCAATACCCAGACGTTTGGCGTCTTCGGTAGTGATTTGTGCGTAGCCAGGTTCATCAGCCAGCGCCGCCAGTGCCGCACAGTTACCGGTCATCGAACGGCAAGAGTATGGACCAACt                                                                                               <  1:146631/250‑1 (MQ=255)                                                                                                                                                       aCCAGATCGTTCCAGACCCCAATCCACCACTGGTAGGTCATGTAAATCGCCCCTTTGTTCGGACGATCGCTGACCTGCGCACGGGTGATAATTTTGTCTTTACGCGAGTGCACCCAAACCAATGCCTCATCTTCAATACCCAGACGTTTGGCGTCTTCGGTATTGATTTTTGCGTAGCCAGGTTCATCAGCCAGCGCCGCCAGTGCCGCACAGTTACCGGTCATCGAACGGCAAGAGTAGTGACCAACTTc                                                                                             >  2:275678/1‑251 (MQ=255)                                                                                                                                                           gCTCGTTACAGGCACCAATCCACCACTGGTAGGTCATGTAAATCGCCCCTTTGTTCGGACGATCGCTGACCTGCGCACGGGTGATAATTTTGTCTTTACGCGAGTGCACCCAAACCAATGCCTCATCTTCAATACCCAGACGTTTGGCGTCGTCGGTATTGATGTGTGCGTAGCCAGGTTCATCAGCCAGCGCCGCCAGTGCCGCACAGTTACCGGTCATCGAACGGCAAGAGTAGTGACCAACTTCACGc                                                                                         <  1:309949/251‑1 (MQ=255)                                                                                                                                                                ttACAGGCCCCAATCCACCCCTGGTAGGTCATGTAAATCGCCCCTTTGTTCGGACGATCGCTGCCCTGCGCACGGGTGATAATTTTGTCTTTACGCGAGTGCACCCAAACCAATGCCTCATCTGCAATCCCCAGACGTTTGGCGTCTTCGGTATTGATTTGTGCGTAGCCAGGTGCATCAGCCAGCGCCGCCAGTGCCGCACAGTTACCGGTCATCGAACGGCAAGAGTAGTGACCAACTTCACGCAcc                                                                                      <  2:319147/249‑1 (MQ=255)                                                                                                                                                                  aCAGGCACCAATCCCCCACTGGTAGGTCATGTAAATCGCCCCTTTGTTCGGACGATCGCTGACCTGCGCACGGGTGATAATTTTGTCTTTACGCGAGTGCACCCAAACCAATGCCTCATCTTCAATACCCAGACGTTTGGCGTCTTCGGTATTGATTTGTGCGTAGCCAGGTTCATCAGCCAGCGCCGCCAGTGCCGCACAGTTACCGGTCATCGAACGGCAAGAGTAGTGACCAACTTCACGCACCggt                                                                                   <  2:156269/250‑3 (MQ=255)                                                                                                                                                                   cAGGCACCAATCCACCACTGGTAGGTCATGTACATCGCCCCTTTGTTCGGACGATCGCTGACCTGCGCACGGGCGATAATTTTGTCTTTACGCGAGTGCACCCAAACCAATGCCTCATCTTCAATACCCAGCCGTTTGGCGTCTTCGGTATTGCTTTGTGCGTAGCCAGGTTCATCAGCCAGCGCCGCCAGTGCCGCACAGTTACCGGTCATCGAACGGCAAGAGGAGTGCCCAACTTCACGCACCGtggc                                                                                 >  2:324772/1‑248 (MQ=255)                                                                                                                                                                    aGGCACCAATCCACCACTGGTAGGTCATGTAAATCGCCCCTTTGTTCGGACGATCGCTGACCTGCGCACGGGTGATAATTTTGTCTTTACGCGAGTGCACCCAAACCAATGCCTCATCTTCAATACCCAGACGTTTGGCGTCTTCGGTATTGATTTGTGCGTAGCCAGGTTCATCAGCCAGCGCCGCCAGTGCCGCACAGTTACCGGTCATCGAACGGCAAGAGTAGTGACCAACTTCACGCACCGTTGAc                                                                                >  1:214291/1‑251 (MQ=255)                                                                                                                                                                            aTCCACCACTGGTAGGTCATGTAAATCGCCCCTTTGTTCGGACGATCGCTGACCTGCGCACGGGTGATAATTTTGTCTTTACGCGAGTGCACCCAAACCAATGCCTCATCTTCAATACCCAGACGTTTGGCGTCTTCGGTATTGATTTGTGCGTAGCCAGGTTCATCAGCCAGCGCCGCCAGt                                                                                                                                            >  1:193320/1‑183 (MQ=255)                                                                                                                                                                                        taggaaaTGTAAATCGCCCCTTTGTTCGGACGATCGCTGCCCTGCGCACGGGAGATAATTTTGTCTTTACGCGAGTGCACACAAACCAATGCCTCATCTTCAATACCCAGACGTTTGGCGTCTTCGGTATTGATTTGTGCGTCGCCAGGTTCATCAGCCAGCGCCGCCAGTGCCGCACAGTTACCGGTCATCGAACGGCAAGAGTAGTGACCAACTTCACGCACCGTTGAAAGTACCATCGGGCACtcgtc                                                            >  2:52984/7‑251 (MQ=255)                                                                                                                                                                                             cATGTAAATCGCCCCTTTGTTCGGACGATCGCTGACCTGCGCACGGGTGATAATTTTGTCTTTACGCGAGTGCACCCAAACCAATGCCTCATCTTCAATACCCAGACGTTTGGCGTCTTCGGTATTGATTTGTGCGTAGCCAGGTTCATCAGCCAGCGCCGCCAGTGCCGCACAGTTACCGGTCATCGAACGGCAAGAGTAGTGACCAACTTCACGCACCGTTGACAGTACCATCGGGTACTCGTCGGTg                                                        >  1:30380/1‑250 (MQ=255)                                                                                                                                                                                                            tttGTTCGGACGATCGCTGACCTGCGCACGGGTGATAATTTTGTCTTTACGCGAGTGCACCCAAACCAATGCCTCATCTTCAATACCCAGACGTTTGGCGTCTTCGGTATTGATTTGTGCGTAGCCAGGTTCATCAGCCAGCGCCGCCAGTGCCGCACAGTTACCGGTCATCGAACGGCAAGAGTAGTGACCAACTTCACGCACCGTTGACAGTACCATCGGGTACTCGTCGGTGAGTTTGTCGATTGgcg                                        >  1:104363/1‑251 (MQ=255)                                                                                                                                                                                                               gTTCGGACGATCGCTGACCTGCGCACGGGTGATAATTTTGTCTTTACGCGAGTGCACCCAAACCAATGCCTCATCTTCAATACCCAGACGTTTGGCGTCTTCGGTATTGATTTGTGCGTAGCCAGGTTCATCAGCCAGCGCCGCCAGTGCCGCACAGTTACCGGTCATCGAACGGCAAGAGTAGTGACCAACTTCACGCACCGTTGACAGTACCATCGGGTACTCGTCGGTGAGTTTGTCGATTGGCGCt                                      >  1:267029/1‑250 (MQ=255)                                                                                                                                                                                                               gTTCGGACCATCGCTGCCCTCCGCACGGGTGCTAATTTTGTCTTTACGCGAGTGCACCCAACCCCATGCCTCATCTTCAATACACAGACGTTTGGCGTCTTCGGTATTGATTTGTGCGTAGCCAGGTTCATCAGCCAGCGCCGCCAGTGCCGCACAGTTACCGGTCATCGAACGGCACGAGTAGTGACCAACTTCACGCCCCGTTGACAGTACCATCGGGTACTCGTCGGTGAGTTTGTCGATTGGCGCt                                      >  2:266860/1‑250 (MQ=255)                                                                                                                                                                                                                ttCGGACGATCGCTGACCTGCGCACGGGTGATAATTTTGTCTTTACGCGAGTGCACCCAAACCAATGCCTCATCTTCAATACCCAGACGTTTGGCGTCTTCGGTATTGATTTGTGCGTAGCCAGGTTCATCAGCCAGCGCCGCCAGTGCCGCACAGTTACCGGTCATCGAACGGCAAGAGTAGTGACCAACTTCACGCACCGTTGACAGTACCATCGGGTACTCGTCGGTGAGTTTGTCGATTGGCGCTAc                                    >  2:93001/1‑251 (MQ=255)                                                                                                                                                                                                                 tCGGACGATCGCTGACCTGCGCACGGGTGATAATTTTGTCTTTACGCGAGTGCACCCAAACCAATGCCTCATCTTCAATACCCAGACGTTTGGCGTCTTCGGTATTGATTTGTGCGTAGCCAGGTTCATCAGCCAGCGCCGCCAGTGCCGCACAGTTACCGGTCATCGAACGGCAAGAGTAGTGACCAACTTCACGCACCGTTGACAGTACCATCGGGTACTCGTCGGTGAGTTTGTCGATTGGCGCTAcc                                   >  1:191750/1‑251 (MQ=255)                                                                                                                                                                                                                 ccGGAAGATCGCTGACCTGCGCACGGGTGATAATTTTGTCTTTACGCGAGTGCACCCAAACCAATGCCTCATCTTCAATACCCAGACGTTTGGCGTCTTCGGTATTGATTTGTGCGTAGCCAGGTTCATCAGCCAGCGCCGCCAGTGCCGCACAGTTACCGGTCATCGAACGGCAAGAGTAGTGACCAACTTCACGCACCGTTGACAGTACCATCGGGTACTCGTCGGTGAGTTTGTCGATTGGCGCTAcc                                   <  2:129870/250‑1 (MQ=255)                                                                                                                                                                                                                    gacggtggcggtccTGCTCACGGGTGATAATTTTGTCTTTACGCGAGTGCACCCAAGCCAATGCCTCATCTTCATTCCCCAGGCGTTTGGCGTCTTCGGTATTGATTTGTGCGTAGCCAGGTTCATCAGCCAGCGCCGCCAGTGCCGCACAGTTACCGGTCATCGAACGGCAAGAGTAGTGACCAACTTCACGCACCGTTGACAGTACCATCGGGTACTCGTCGGTGAGTGTGTCGATTGGCGCTACCCAg                                <  1:324772/239‑1 (MQ=255)                                                                                                                                                                                                                          cgcggaccggcgctggggTGTTAATTTTTGCTTTACGCGAGTCCACCCAATCCAATGCCTCATCTTCAATCCCCAGACGTTGGGCGTCTTCGGTATTGAGTTGTGCGGAGCCAGGTTCACCAGCCAGCGCCGCCAGTGCCGCACAGTTACCGGTCATCGAACGGCAAGAGTAGTGACCAACTTCACGCACCGTGGACAGTACCATCGGGTACTCGTCGGTGAGTTTGTCGATTGGCGCTACCCAGTCGCAg                          <  2:416339/236‑1 (MQ=255)                                                                                                                                                                                                                                  tGCGCACGGGTGATAATTTTGTCTTTACGCGAGTGCACCCAAACCAATGCCTCATCTTCAATACCCAGACGTTTGGCGTCt                                                                                                                                                                                            <  1:61123/81‑1 (MQ=255)                                                                                                                                                                                                                                  tGCGCACGGGAGATCATTTTGTCTTTACGAGAGTGCACCCAAAACAATGCCTCATCTTCAATACCCAGACGTTTGGCGTCt                                                                                                                                                                                            >  2:61123/1‑81 (MQ=255)                                                                                                                                                                                                                                       agggggtaTAATTTTGTCTTTACGCGAGGTCACCAAAACCAATGCCGCATCTGCAATACCCAGACGGTTGGCGGCTGCGGTATTGATTTGGGCGTAGCCAGGTTCATCAGCAAGCGCCGCCAGTGCCGCACAGTTACCGGTCATCGAACGGCAAGAGTAGTGACCAACTTCACGCACCGTTGACAGTACCATCGGGTACTCGTCGGTGAGTTTGTCGATTGGCGCTACCCAGTCGCAGGTGAAGAACTgcg             <  2:394210/244‑1 (MQ=255)                                                                                                                                                                                                                                                 aTTTTGTCTTTACGCGAGTGCCCCCAACCCAATGCCTCATCTGCAATACCCAGACGTTTGGCGTTTCCGGTATTGATTTGTTCGGAGCCAGGTTCATCAGCCAGCGCCGCCAGTGCCGCACAGTTACCGGTCATCGAACGGCAAGAGTAGTGACCAACTTCACGCACCGTTGACAGTACCATCGGGTACTCGTCGGTGAGTTTGTCGATTGGCGCTACCCAGTCGCAGGTGAAGAACTGCGCCAGACCGtt   <  1:108872/251‑1 (MQ=255)                                                                                                                                                                                                                                                   tttGTCTTTACGCGAGGGCCCCCAACCCAATGTCTCATCTTCAATACCCAGCCGTTGGGCGTCTGCGGTATTGATTTGTGCGTATCCAGGTTCTTCAGCCAGCGCCGCCAGTGCCGCACAGTTACCGGTCATCGAACGGCAAGAGTAGTGACCAACTTCACGCACCGTTGACAGTACCATCGGGTACTCGTCGGTGAGTTTGTCGATTGGCGCTACCCAGTCGCAGGTGAAGAACTGCGCCAGACCGTTc  <  1:150538/250‑1 (MQ=255)                                                                                                                                                                                                                                                       |                                                                                                                                                                                                                                                      GCCAGTGCCGCTTCGCGCAGGCGAGTTTTCAACTTGTTGTACTCGTCAATCACGTACTGCTCGGCGGCGCGCTGATCGGCGATCGGCTCGACGCGAACGGCGCAGTATTTGTACTCCGGCGTTTTCGTAATCGGGCTTAAGTTTTCGGTAACCAGCTCGTTACAGGCACCAATCCACCACTGGTAGGTCATGTAAATCGCCCCTTTGTTCGGACGATCGCTGACCTGCGCACGGGTGATAATTTTGCCTTTACGCGAGTGCACCCAAACCAATGCCTCATCTTCAATACCCAGACGTTTGGCGTCTTCGGTATTGATTTGTGCGTAGCCAGGTTCATCAGCCAGCGCCGCCAGTGCCGCACAGTTACCGGTCATCGAACGGCAAGAGTAGTGACCAACTTCACGCACCGTTGACAGTACCATCGGGTACTCGTCGGTGAGTTTGTCGATTGGCGCTACCCAGTCGCAGGTGAAGAACTGCGCCAGACCGTTC  >  NC_000913/4297223‑4297714

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