Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	NC_000913	4,297,782	G→A	P529S (CCT→TCT)	fdhF ←	formate dehydrogenase H

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	4,297,782	0	G	A	98.0%	151.9 / ‑4.9	51	P529S (CCT→TCT)	fdhF	formate dehydrogenase H
Reads supporting (aligned to +/- strand):  ref base G (1/0);  new base A (27/23);  total (28/23)
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.38e-01

GATTTGTGCGTAGCCAGGTTCATCAGCCAGCGCCGCCAGTGCCGCACAGTTACCGGTCATCGAACGGCAAGAGTAGTGACCAACTTCACGCACCGTTGACAGTACCATCGGGTACTCGTCGGTGAGTTTGTCGATTGGCGCTACCCAGTCGCAGGTGAAGAACTGCGCCAGACCGTTCGGGGTATCAAACTTCTCTTTAAACAGATAAGAAGTCCCCTGATCGGCATCTGAAGTATCGCGGCAAGGCCACTGAATGAAGCCCAGTTCGCCCATTTTCTCGTAAGTCGCACCGTAGAAATCCGGGCACAGATGACGCAACTCATCCCAGATCTCCTGGGTGTTGTTGTAGTGCATCGGATAACCCATACGGGTGGCGATTTCACTGATGATTTGCCAGTCCGTTTTCAGATCCCATTTCGGTTCAACCGCCTTGAAGAAACGCTGGAAGCCACGGTCAGCCGCAGTAAACACGCCTTCATGCTCGCCCCACGA  >  NC_000913/4297537‑4298028                                                                                                                                                                                                                                                      |                                                                                                                                                                                                                                                       gATTTGTGCGTAGCCAGGTTCATCTGCCAGCGCCGCCAGTGCCGCACAGTTACCGGTCATCGAACGGCAAGAGTAGTGACCAACTTCACGCTCCGTTGACAGTACCATCGGGTACTCGTCGGTGAGTTTGTCGATTGGCGCTGCCCAGTCGCAGGTGAAGAACTGCGCCAGACCGTGCGGGGTATCAAACTTCTCTTTAAACAGATAAGAAGTCCCCTGATCGGCATCTGAAGTATCGCGGCAAGACCAc                                                                                                                                                                                                                                                    <  1:63725/250‑1 (MQ=255)   gttgtgcgatgcctggtctacctaccagcGCCGCCAGGGCCGCACAGTTACCGGTCATCGAACGGCAAGAGTAGTGGCCAACTTCACGCACCGTTGACAGTACCATCGGGTACTCGTCGGTGAGTTTGTCGATTGGCGCTACCCAGTCGCAGGTGAAGAACTGCGCCAGACCGTTCGGGGTATCAAACTTCTCTTTAAACAGATAAGAAGTCCCCTGATCGGCATCTGAAGTATCGTGGCAAGACCACtg                                                                                                                                                                                                                                                  <  2:341098/226‑1 (MQ=255)                     attcaGCCAGCGCCGCCAGGGCCGCACAGTTTACGGTCATCGAACGGCAAGAGTAGTGACCAACGTCACGCACCGTTGACAGTACCATCGGGTACTCGTCGGTGAGTTTGTCGATTGGCGCTTCCCAGTCGCAGGTGAAGAACTGCGCCAGCCCGTGCGGGGTATCAAACTTCTCTGTAAACAGATAAGAAGGCCCCTGATCGGCATCTGAAGTATCGCGGCAAGACCACTGAATGAAGCCCAGTTCGcc                                                                                                                                                                                                                                <  2:89063/248‑1 (MQ=255)                       tcaGCCAGCGCCGCCAGTGCCGCACAGTTACCGGTCATCGAACGGCAAGAGTAGTGACCAACTTCACGCACCGTTGACAGTACCATCGGGTACTCGTCGGTGAGTTTGTCGATTGGCGCTACCCAGTCGCAGGTGAAGAACTGCGCCAGACCGTTCGGGGTATCAAACTTCTCTTTAAACAGATAAGAAGTCCCCTGATCGGCATCTGAAGTATCGCGGCAAGACCACTGAATGAAGCCCAGTTCGCCCAt                                                                                                                                                                                                                             >  1:200354/1‑251 (MQ=255)                         agccagcGCCGCCAGTGCCGCACAGTTACCGGTCATCGAACGGCAAGAGTAGTGACCAACTTCACGCACCGTTGACAGTACCATCGGGTACTCGTCGGTGAGTTTGTCGATTGGCGCTACCCAGTCGCAGGTGAAGAACTGCGCCAGACCGTTCGGGGTATCAAACTTCTCTTTAAACAGATAAGAAGTCCCCTGATCGGCATCTGAAGTATCGCGGCAAGACCACTGAATGAAGCCCAGTTCGCCCAttt                                                                                                                                                                                                                           >  2:168737/1‑251 (MQ=255)                                             cacaGTTACCGGTCATCGAACGGCAAGAGTAGTGACCAACTTCACGCACCGTTGACAGTACCATCGGGTACTCGTCGGTGAGTTTGTCGATTGGCGCTACCCAGTCGCAGGTGAAGAACTGCGCCAGACCGTTCGGGGTATCAAACTTCTCTTTAAACAGATAAGAAGTCCCCTGATCGGCATCTGAAGTATCGCGGCAAGACCACTGAATGAAGCCCAGTTCGCCCATTTTCTCGTAAGTCGCACCGTAg                                                                                                                                                                                                       <  2:200354/251‑1 (MQ=255)                                                 gTTACCGGTCATCGAACGGCAAGAGTAGTGACCAACTTCACGCACCGTTGACAGTACCATCGGGTACTCGTCGGTGAGTTTGTCGATTGGCGCTACCCAGTCGCAGGTGAAGAACTGCGCCAGACCGTTCGGGGTATCAAACTTCTCTTTAAACAGATAAGAAGTCCCCTGATCGGCATCTGAAGTATCGCGGCAAGACCACTGAATGAAGCCCAGTTCGCCCATTTTCTCGTAAGTCGCACCGTAGAAAt                                                                                                                                                                                                   <  1:114111/251‑1 (MQ=255)                                                              gAACGGCAAGAGTAGTGACCAACTTCACGCACCGTTGACAGTACCATCGGGTACTCGTCGGTGAGTTTGTCGATTGGCGCTACCCAGTCGCAGGTGAAGAACTGCGCCAGACCGTTCGGGGTATCAAACTTCTCTTTAAACAGATAAGAAGTCCCCTGATCGGCATCTGAAGTATCGCGGCAAGACCACTGAATGAAGCCCAGTTCGCCCATTTTCTCGTAAGTCGCACCGTAGAAATCCGGGCACAGATg                                                                                                                                                                                      <  1:311988/251‑1 (MQ=255)                                                                 cGGCAAGAGTAGTGACCAACTTCACGCACCGTTGACAGTACCATCGGGTACTCGTCGGTGAGTTTGTCGATTGGCGCTACCCAGTCGCAGGTGAAGAACTGCGCCAGACCGTTCGGGGTATCAAACTTCTCTTTAAACAGATAAGAAGTCCCCTGATCGGCATCTGAAGTATCGCGGCAAGACCACTGAATGAAGCCCAGTTCGCCCATTTTCTCGTAAGTCGCACCGTAGAAATCCGGGCACAGATGACg                                                                                                                                                                                   <  1:353949/251‑1 (MQ=255)                                                                   gCAAGAGTAGTGACCAACTTCACGCGCCGTTGACAGTACCATCGGGTACTCGTCGGTGAGTTTGTCGATTGGCGCTACTCAGTCGCAGGTGAAGAACTGCGCCAGATCGTTCGGGGTATTAAACTTCTCTTTAAACAGATAAGAAGTCCCCTGATCGGCATCTGAAGTAGCGCGGCAAGACCACTGAATGAAGCCTAGTTCTCCCATTTGCTCGTAAGTCGCACCGTAGAAATCCGGGCACAGATGACGCa                                                                                                                                                                                 <  1:161139/251‑1 (MQ=255)                                                                              gACCAACTTCACGCACCGTTGACAGTACCATCGGGTACTCGTCGGTGAGTTTGTCGATTGGCGCTACCCAGTCGCAGGTGAAGAACTGCGCCAGACCGTTCGGGGTATCAAACTTCTCTTTAAACAGATAAGAAGTCCCCTGATCGGCATCTGAAGTATCGCGGCAAGACCACTGAATGAAGCCCAGTTCGCCCATTTTCTCGTAAGTCGCACCGTAGAAATCCGGGCACAGATGAc                                                                                                                                                                                    <  1:304597/237‑1 (MQ=255)                                                                              gACCAACTTCACGCACCGTTGACAGTACCATCGGGTACTCGTCGGTGAGTTTGTCGATTGGCGCTACCCAGTCGCAGGTGAAGAACTGCGCCAGACCGTTCGGGGTATCAAACTTCTCTTTAAACAGATAAGAAGTCCCCTGATCGGCATCTGAAGTATCGCGGCAAGACCACTGAATGAAGCCCAGTTCGCCCATTTTCTCGTAAGTCGCACCGTAGAAATCCGGGCACAGATGAc                                                                                                                                                                                    >  2:304597/1‑237 (MQ=255)                                                                                 cAACTTCACGCACCGTTGACAGTACCATCGGGTACTCGTCGGTGAGTTTGTCGATTGGCGCTACCCAGTCGCAGGTGAAGAACTGCGCCAGACCGTTCGGGGTATCAAACTTCTCTTTAAACAGATAAGAAGTCCCCTGATCGGCATCTGAAGTATCGCGGCAAGACCACTTAATGAAGCCCAGTTCGCCCATTTTCTCGTAAGTCGCACCGTAGAAATCCGGGCACAGATGACGCAACTCATCCCAGAtc                                                                                                                                                                   <  1:274655/251‑1 (MQ=255)                                                                                   aCTTCACGCACCGTTGACAGTACCATCGGGTACTCGTCGGTGAGTTTGTCGATTGGCGCTACCCAGTCGCAGGTGAAGAACTGCGCCAGACCGTTCGGGGTATCAAACTTCTCTTTAAACAGATAAGAAGTCCCCTGATCGGCATCTGAAGTATCGCGGCAAGACCACTGAATGAAGCCCAGTTCGCCCATTTTCTCGTAAGTCGCACCGTAGAAATCCGGGCACAGATGACGCAACTCATCCCAGAtctc                                                                                                                                                                 >  2:414717/1‑251 (MQ=255)                                                                                      tCACGCTCCGTTGACAGTACCATCGGGTACTCGGCGGTGAGTTTGTCGATTGGCGCTACCCAGTCGCAGGTGAAGAACTGCGCCAGACCGTTCGGGGTATCAAACTTCTCTGTAAACAGATAAGAAGTCCCCTGATCGGCATCTGAAGTATCGCGGCAAGACCACTGAATGAAGCCCAGTTCGCCCATTTTCTCGTAAGTCGCACCGTAGAAATCCGGGCACAGATGACGCAACTCATCCCAGATCTCCt                                                                                                                                                               <  2:148685/250‑1 (MQ=255)                                                                                      tCACGCACCGTTGACAGTACCATCGGGTACTCGTCGGTGAGTTTGTCGATTGGCGCTACCCAGTCGCAGGTGAAGAACTGCGCCAGACCGTTCGGGGTATCAAACTTCTCTTTAAACAGATAAGAAGTCCCCTGATCGGCATCTGAAGTATCGCGGCAAGACCACTGAATGAAGCCCAGTTCGCCCATTTTCTCGTAAGTCGCACCGTAGAAATCCGGGCACAGATGACGCAACTCATCCCAGATCTCCt                                                                                                                                                               <  1:279178/250‑1 (MQ=255)                                                                                          gCACCGTTGACAGTACCCTCAGGTACTAGTCGGTGAGTTTGTCGATTGGCGCTACCCAGTCGCAGGTGAAGAACTGCGACAGAACGTTCGGGGTATCAAACTTCTATTTAAACAGCTAAGAAGTCCCCTGATCGGCATCTGAAGTATCGCGGCAAGACCACTGAATGAAGCCCAGCTCGCCCCTTTTCTCGTAAGACGCACCGTACAAGAACGGGCACAGATGACACCACTCATCACCGATCTCCTggg                                                                                                                                                            >  2:436225/1‑249 (MQ=255)                                                                                          gCACCGTTGACAGTACCATCGGGTACTCGTCGGTGAGTTTGTCGATTGGCGCTACCCAGTCGCAGGTGAAGAACTGCGCCAGACCGTTCGGGGTATCAAACTTCTCTTTAAACAGATAAGAAGTCCCCTGATCGGCATCTGAAGTATCGCGGCAAGACCACTGAATGAAGCCCAGTTCGCCCATTTTCTCGTAAGTCGCACCGTAGAAATCCGGGCACAGATGACGCAACTCATCCCAGATCTCCTGGgtg                                                                                                                                                          >  1:26911/1‑251 (MQ=255)                                                                                                               ggTACTCGTCGGTGAGTTTGTCGATTGGCGCTACCCAGTCGCAGGTGAAGAACTGCGCCAGACCGTTCGGGGTATCAAACTTCTCTTTAAACAGATAAGAAGTCCCCTGATCGGCATCTGAAGTATCGCGGCAAGACCACTGAATGAAGCCCAGTTCGCCCATTTTCTCGTAAGTCGCACCGTAGAAATCCGGGCACAGATGACGCAACTCATCCCAGATCTCCTGGGTGTTGTTGTAGTGCATCGGATa                                                                                                                                      >  2:90691/1‑250 (MQ=255)                                                                                                                          gTGAGTTTTTCGATTGGCGCTACCCAGTCGCAGGGGAAGAACTGCGCCAGACCGTTCGGGGTATCAAACTTCTCCTTAAACAGATAAGAAGTCCCCTGATCGGCATCTGAAGTATCGCGGCAAGACCACTGAATGAAGCCCAGGTCGCCCATTTTCTCGTAAGTCGCACCGTAGAAATCCGGGCACAGATGACGCAACTCATCCCAGATCTCCTGGGTGTTGTTGTAGTGCATCGGATAACCCATACGGGt                                                                                                                          <  2:484068/251‑1 (MQ=255)                                                                                                                                  gTCGATTGGCGCTACCCAGTCGCAGGTGAAGAACTGCGCCAGTCCGTTCGGGGTATCAAACTTCTCTTTAAACAGATAAGAAGTCCCCTGATCGGCATCTGAAGTATCGCGGCAAGACCACTGAATGAAGCCCAGTTCGCCCATTTTCTCGTAAGTCGCACCGTAGAAATCCGGGCACAGATGACGCAACTCATCCCAGATCTCCTGGGTGTTGTTGTAGTGCATCGGATAACCCATACGGGTGGCGAttt                                                                                                                  <  2:263383/251‑1 (MQ=255)                                                                                                                                    cGATTGGCGCTACCCAGTCGCAGGTGAAGAACTGCGCCAGACCGTTCGGGGTATCAAACTTCTCTTTAAACAGATAAGAAGTCCCCTGATCGGCATCTGAAGTATCGCGGCAAGACCACTGAATGAAGCCCAGTTCGCCCATTTTCTCGTAAGTCGCACCGTAGAAATCCGGGCACAGATGACGCAACTCATCCCAGCTCTCCTGGGTGTAGTTGTAGTGCATCGGATAACCCATACGGGAGGCGATTTc                                                                                                                 >  2:135704/1‑250 (MQ=255)                                                                                                                                      aTTGGCGCTACCCAGTCGCAGGTGAAGAACTGCGCCAGACCGTTCGGGGTATCAAACTTCTCTTTAAACAGATAAGAAGTCCCCTGATCGGCATCTGAAGTATCGCGGCAAGACCACTGAATGCAGCCCAGTTCGCCCATTTTCTCGTAAGTCGCACCGTAGAAATCCGGGCACAGATGACGCAACTCATCCCAGATCTCCTGGGTGTTGTTGTAGTGCATCGGCTAACCCATACGGGTGGCGTTTTCACt                                                                                                              >  2:308015/1‑251 (MQ=255)                                                                                                                                            gcTACCCAGTCGCAGGTGAAGAACTGCGCCAGACCGTTCGGGGTATCAAACTTCTCTTTAAACAGATAAGAAGTCCCCTGATCGGCATCTGAAGTATCGCGGCAAGACCACTGAATGAAGCCCAGTTCGCCCATTTTCTCGTAAGTCGCACCGTAGAAATCCGGGCACAGATGACGCAACTCATCCCAGATCTCCTgg                                                                                                                                                             <  2:97201/198‑1 (MQ=255)                                                                                                                                            gcTACCCAGTCGCAGGTGAAGAACTGCGCCAGACCGTTCGGGGTATCAAACTTCTCTTTAAACAGATAAGAAGTCCCCTGATCGGCATCTGAAGTATCGCGGCAAGACCACTGAATGAAGCCCAGTTCGCCCATTTTCTCGTAAGTCGCACCGTAGAAATCCGGGCACAGATGACGCAACTCATCCCAGATCTCCTgg                                                                                                                                                             >  1:97201/1‑198 (MQ=255)                                                                                                                                                cccAGTCGCAGGTGAAGAACTGCGCCAGACCGTTCGGGGTATCAAACTTATCTTTAAACAGATAAGAAGTCCCCTGATCGGCATCTGAAGTATCGCGGCAAGACCACTGAATGAAGCCCAGTTCGCCCATTTTCTCGTAAGTCGCACCGTAGAAATCCGGGCACAGATGACGCAACTCATCCCAGATCTCCTGGGTGTTGTTGTAGTGCATCGGATAACCCATACGGGTGGCGATTTCACTGATGATTTGc                                                                                                    >  1:383837/1‑251 (MQ=255)                                                                                                                                                       gCAGGTGAAGAACTGCGCCAGACCGTTCGGGGTATCAAACTTCTCTTTAAACAGATAAGAAGTCCCCTGATCGGCATCTGAAGTATCGCGGCAAGACCACTGAATGAAGCCCAGTTCGCCCATTTTCTCGTAAGTCGCACCGTAGAAATCCGGGCACAGATGACGCAACTCATCCCAGATCTCCTGGGTGTTGTTGTAGTGCATCGGATAACCCATACGGGTGGCGATTTCACTGATGATTTGCCAGTCCg                                                                                             >  1:151580/1‑251 (MQ=255)                                                                                                                                                           gTGAAGAACTGCGCCAGACCGTTCGGGGTATCAAACTTCTCTTTAAACAGATAAGAAGTCCCCTGATCGGCATCTGAAGTATCGCGGCAAGACCACTGAATGAAGCCCAGTTCGCCCATTTTCTCGTAAGTCGCACCGTAGAAATCCGGGCACAGATGACGCAACTCATCCCAGATCTCCTGGGTGTTGTTGTAGTGCATCGGATAACCCATACGGGTGGCGATTTCACTGATGATTTGCCAGTCCGtttt                                                                                         <  2:62430/251‑1 (MQ=255)                                                                                                                                                             gaagaaCTGCGCCAGACCGTTCGGGGTATCAAACTTCTCTTTAAACAGATAAGAAGTCCCCTGATCGGCATCTGAAGTATCGCGGCAAGACCACTGAATGAAGCCCAGTTCGCCCATTTTCTCGTAAGTCGCACCGTAGAAATCCGGGCACAGATGACGCAACTCATCCCAGATCTCCTGGGTGTTGTTGTAGTGCATCGGATAACCCATACGGGTGGCGATTTCACTGATGATTTGCCAGTCCGTTTTc                                                                                        <  1:110878/250‑1 (MQ=255)                                                                                                                                                                         cAGACCGTTCGGGGTATCAAACTTCTCTTTAAACAGATAAGAAGTCCCCTGATCGGCATCTGAAGTATCGCGGCAAGACCACTGAATGAAGCCCAGTTCGCCCATTTTCTCGTAAGTCGCACCGTAGAAATCCGGGCACAGATGACGCAACTCATCCCAGATCTCCTGGGTGTTGTTGTAGTGCATCGGATAACCCATACGGGTGGCGCTTTCACTGATGATTTGCCAGTCCGTTTTCAGATCCCATTTCg                                                                           >  2:378107/1‑251 (MQ=255)                                                                                                                                                                          aGACCGTTCGGGGTATCAAACTTCTCTTTAAACAGATAAGAAGTCCCCTGATCGGCATCTGAAGTATCGCGGCAAGACCACTGAATGAAGCCCAGTTCGCCCATTTTCTCGTAAGTCGCACCGTAGAAATCCGGGCACAGATGACGCAACTCATCCCAGATCTCCTGGGTGTTGTTGTAGTGCATCGGATAACCCATACGGGTGGCGATTTCACTGATGATTTGCCAGTCCGTTTTCAGATCCCATTTCgg                                                                          <  1:29144/251‑1 (MQ=255)                                                                                                                                                                               gTTCGGGGTATCAAACTTCTCTTTAAACAGATAAGAAGTCCCCTGATCGGCATCTGAAGTATCGCGGCAAGGCCACTGAATGAAGCCCAGTTCGCCCATTTTCTCGTAAGTCGCACCGTAGAAATCCGGGCACAGATGACGCAACTCATCCCAGATCTCCTGGGTGTTGTTGTAGTGCCTCGGATAACCCATACGGGTGGCGATTTCACTGATGCTTTGCCCGTCCGGTTTCAGATCCCATTTCGGttaac                                                                     >  2:118344/1‑248 (MQ=255)                                                                                                                                                                                    gggTATCAAACTTCTCTTTAAACAGATAAGAAGTCCCCTGATCGGCATCTGAAGTATCGCGGCAAGACCACTGAATGAAGCCCAGTTCGCCCGTTTTCTCGTAAGTCGCACCGTAGAAATCCGGGCACAGATGACGCAACTCATCCCAGATCTCCTGGGTGTTGTTGTAGTGCATCGGATAACCCATACGGGTGGCGCTTTCACTGATGATTTGCCAGTCCGTTTTCAGATCCCATTTCGGTCCAACCGcc                                                                >  2:425507/1‑251 (MQ=255)                                                                                                                                                                                     ggTATCAAACTTCTCTTTAAACAGATAAGAAGTCCCCTGATCGGCATCTGAAGTATCGCGGCAAGACCACTGAATGAAGCCCAGTTCGCCCATTTTCTCGTAAGTCGCACCGTAGAAATCCGGGCACAGATGACGCAACTCATCCCAGATCTCCTGGGTGTTGTTGTAGTGCATCGGATAACCCATACGGGTGGCGATTTCACTGATGATTTGCCAGTCCGTTTTCAGATCCCATTTc                                                                            <  2:30365/238‑1 (MQ=255)                                                                                                                                                                                     ggTATCAAACTTCTCTTTAAACAGATAAGAAGTCCCCTGATCGGCATCTGAAGTATCGCGGCAAGACCACTGAATGAAGCCCAGTTCGCCCATTTTCTCGTAAGTCGCACCGTAGAAATCCGGGCACAGATGACGCAACTCATCCCAGATCTCCTGGGTGTTGTTGTAGTGCATCGGATAACCCATACGGGTGGCGATTTCACTGATGATTTGCCAGTCCGTTTTCAGATCCCATTTc                                                                            >  1:30365/1‑238 (MQ=255)                                                                                                                                                                                      gTATCAAACTTCTCTTTAAACAGATAAGAAGTCCCCTGATCGGCATCTGAAGTATCGCGGCAAGACCACTGAATGAAGCCCAGTTCGCCCATTTTCTCGTAAGTCGCACCGTAGAAATCCGGGCACAGATGACGCAACTCATCCCAGATCTCCTGGGTGTTGTTGTAGTGCATCGGATAACCCATACGGGTGGCGATTTCACTGATGATTTGCCAGTCCGTTTTCAGATCCCATTTAGGTTCAACCGCCtt                                                              >  2:398903/1‑251 (MQ=255)                                                                                                                                                                                      gTATCAAACTTATCTTTAAAAAGATAAGAAGTCCCCTGATCGGCATCTGAAGTATCGCGGCAAGACAAATGAATGAAGCCCAGTGCGCCCATTTTATCGTAAGTCACACCGTAGAAATCCGGGCACAGATGACGCAAATCATCCCAGATCTCCTGGGGGTTGTGGTATGGCATCGGATACCCCATACGGGTGGCGATTTCAATGATGATTTTCAAGTCCGGTTTCAGAACCCCTTTCGGTTCACCCACCtt                                                              >  2:437831/1‑251 (MQ=255)                                                                                                                                                                                               ttCTCTTTAAACAGATAAGAAGTCCCCTGATCGGCATCTGAAGTATCGCGGCAAGACCACTGAATGAAGCCCAGTTCGCCCATTTTCTCGTAAGTCGCACCGTAGAAATCCGGGCACAGATGACGCAACTCATCCCAGATCTCCTGGGTGTTGTTGTAGTGCATCGGATAACCCATACGGGTGGCGCTTTCACTGTTGATTTGCCAGTCCGTTTTCAGTTCCCATTTCGGTTCAACCGCCTTGAAGAAACg                                                     >  2:374351/1‑251 (MQ=255)                                                                                                                                                                                                        aaCAGATAAGAAGTCCCCTGATCGGCATCTGAAGTATCGCGGCAAGACCACTGAATGAAGCCCAGTTCGCCCATTTTCTCGTAAGTCGCACCGTAGAAATCCGGGCACAGATGACGCAACTCATCCCAGATCTCCTGGGTGTTGTTGTAGTGCATCGGATAACCCATACGGGTGGCGATTTCACTGATGATTTGCCAGTCCGTTTTCAGATCCCATTTCGGTTCAACCGCCTTGAAGAAACGCTGGAAGcc                                            <  1:308015/251‑1 (MQ=255)                                                                                                                                                                                                         aCAGATAAGAAGTCCCCTGATCGGCATCTGAAGTATCGCGGCAAGACCACTGAATGAAGCCCAGTTCGCCCATTTTCTCGTAAGTCGCACCGTAGAAATCCGGGCACAGATGACGCAACTCATCCCAGATCTCCTGGGTGTTGTTGTAGTGCATCGGATAACCCATACGGGTGGCGATTTCACTGATGATTTGCCAGTCCGTTTTCAGATCCCATTTCGGTTCAACCGCCTTGAAGAAACGCTGGAAGccc                                           >  1:173680/1‑250 (MQ=255)                                                                                                                                                                                                          cAGATACGAAGTCCCCTGATCGGCATCTGAAGTATCGCGGCAAGACCACTGAATGAAGCCCAGTTCGCCCATTTTCTCGTAAGTCGCACCGTAGAAATCCGGGCACAGATGACGCAACTCATCCCAGATCTCCTGGGTGTTGTTGTAGTGCATCGGATAACCCAGACGGGTGGCGATTTCCCTGATGATTTGCCAGTCCGTTTTCAGATCCCATTTCGGTTCAACCGCCTTGAAGAAACGCTGGAAGCCAc                                          >  2:139164/1‑251 (MQ=255)                                                                                                                                                                                                           aGATAAGAAGTCCCCTGATCGGCATCTGAAGTATCGCGGCAAGACCACTGAATGAAGCCCAGTTCGCCCATTTTCTCGTAAGTCGCACCGTAGAAATCCGGGCACAGATGACGCAACTCATCCCAGATCTCCTGGGTGTTGTTGTAGTGCATCGGATAACCCATACGGGTGGCGATTTCACTGATGATTTGCCAGTCCGTTTTCAGATCCCATTTCGGTTCAACCGCCTTGAAGAAACGCTGGAAGCCACg                                         >  1:425780/1‑251 (MQ=255)                                                                                                                                                                                                                   agTCCCCTGATCGGCATGTGAAGTATCGCGGCAAGACCACTGAATGAAGCCCAGTTCGCCCATTTTCTCGAAAGTAGCACCGTAGAAATACGGGCACAGATGCAGCAACTCATCCCAGATCTCCTGGGTGTTGTTGTAGTGCATCGGATAACCCATACGGGAGGCGATTTCACTGATGATTTGACAGTCCGTTTTCAGATCCCATTTCGGTTCTACCGACTTGAAGAAACGCTGGAAGCCACGGTCAGCCg                                 >  1:71661/1‑251 (MQ=255)                                                                                                                                                                                                                        ccTGATCGGCATCTGAAGTATCGCGGCAAGACCACTGAATGAAGCCCAGTTCGCCCATTTTCTCGTAAGTCGCACCGTAGAAATCCGGGCACAGATGACGCAACTCATCCCAGATCTCCTGGGTGTTGTTGTAGTGCATCGGATAACCCATACGGGTGGCGATTTCACTGATGATTTGCCAGTCCGTTTTCAGATCCCATTTCGGTTCAACCGCCTTGAAGAAACGCTGGAAGCCACGGTCAGCCGCAGt                             >  1:485833/1‑250 (MQ=255)                                                                                                                                                                                                                           gATCGGCATCTGAAGTATCGCGGCAAGACCCCTGAATGAAGCCCAGTTCGCCCATTTTCTCGTAAGTCGCACCGTAGAAATCCGGGCACAGATGACGCAACTCATCCCAGATCTCCTGGGTGTTGTTGTAGTGCATCGGATAACCCATACGGGTGGCGATTTCACTGATGATTTGCCAGTCCGTTTTCAGATCCCATTTCGGTTCAACCGCCTTGAAGAAACGCTGGAAGCCACGGTCAGCCGCAGTaaa                          <  1:398903/250‑1 (MQ=255)                                                                                                                                                                                                                            aTCGGCATCTGAAGTATCGCGGCAAGACCACTGAATGAAGCCCAGTTCGCCCATTTTCTCGTAAGTCGCACCGTAGAAATCCGGGCACAGATGACGCAACTCATCCCAGATCTCCTGGGTGTTGTTGTAGTGCATCGGATAACCCATACGGGTGGCGATTTCACTGATGATTTGCCAGTCCGTTTTCAGATCCCATTTCGGTTCAACCGCCTTGAAGAAAAGCTGGAAGCCACGGTCCGCCGCAGTAAac                         >  2:789/1‑250 (MQ=255)                                                                                                                                                                                                                                      gAAGTATCGCGGCAAGACCACGGAATGAAGCCCAGTTCGCCCATTTTCTCGTAAGTCGCACCGTAGAAATCCGGGCACAGATGACGCAACTCATCCCAGATCTCCTGGGTGTTGTTGTAGTGCATCGGATAACCCATACGGGTGGCGATTTCACTGATGATTTGCCAGTCCGTTTTCAGATCCCATTTCGGTTCAACCGCCTTGAAGAAACGCTGGAAGCCACGGTCAGCCGCAGTAAACACGCCTTCAt               <  1:168737/250‑1 (MQ=255)                                                                                                                                                                                                                                        aGTATCGCGGCAAGACCACTGAATGAAGCCCAGTTCGCCCATTTTCTCGTAAGTCGCACCGTAGAAATCCGGGCACAGATGACGCAACTCATCCCAGATCTCCTGGGTGTTGTTGTAGTGCATCGGATAACCCATACGGGTGGCGATTTCACTGATGATTTGCCAGTCCGTTTTCAGATCCCATTTCGGTTCAACCGCCTTGAAGAAACGCTGGAAGCCACGGTCAGcc                                  >  2:379463/1‑229 (MQ=255)                                                                                                                                                                                                                                        aGTATCGCGGCAAGACCACTGAATGAAGCCCAGTTCGCCCATTTTCTCGTAAGTCGCACCGTAGAAATCCGGGCACAGATGACGCAACTCATCCCAGATCTCCTGGGTGTTGTTGTAGTGCATCGGATAACCCATACGGGTGGCGATTTCACTGATGATTTGCCAGTCCGTTTTCAGATCCCATTTCGGTTCAACCGCCTTGAAGAAACGCTGGAAGCCACGGTCAGcc                                  <  1:379463/229‑1 (MQ=255)                                                                                                                                                                                                                                                  cAAGACCACTGAATGAAGCCCAGTTCGCCCATTTTCTCGTAAGTCGCACCGTAGAAATCCGGGCACAGATGACGCAACTCATCCCAGATCTCCTGGGTGTTGTTGTAGTGCATCGGATAACCCATACGGGTGGCGATTTCACTGATGATTTGCCAGTCCGTTTTCAGATCCCATTTCGGTTCAACCGCCTTGAAGAAACGCTGGAAGCCACGGTCAGCCGCAGTAAACACGCCTTCATGCTCGCCCCacgc  >  1:362836/1‑250 (MQ=255)                                                                                                                                                                                                                                                  cAAGACCACTGAATGAAGCCCAGTTCGCCCATTTTCTCGTAAGTCGCACCGTAGAAACCCGGGCACAGATGACGCAACTCATCCCAGATCTCCTGGGTGTTGTTGTAGTGCATCGGATAACCCATACGGGTGGCGATTTCTATGATGTTTTGCCAGTCCGTTTTCAGATCCCATTTCGGTTCAACCGCCTTGAAGAAACGCTGGAAGCCACGGTCCGCCGCAGTAACCACGCCTACAAGCTAGCCCCacg   >  2:1596/1‑250 (MQ=255)                                                                                                                                                                                                                                                      |                                                                                                                                                                                                                                                       GATTTGTGCGTAGCCAGGTTCATCAGCCAGCGCCGCCAGTGCCGCACAGTTACCGGTCATCGAACGGCAAGAGTAGTGACCAACTTCACGCACCGTTGACAGTACCATCGGGTACTCGTCGGTGAGTTTGTCGATTGGCGCTACCCAGTCGCAGGTGAAGAACTGCGCCAGACCGTTCGGGGTATCAAACTTCTCTTTAAACAGATAAGAAGTCCCCTGATCGGCATCTGAAGTATCGCGGCAAGGCCACTGAATGAAGCCCAGTTCGCCCATTTTCTCGTAAGTCGCACCGTAGAAATCCGGGCACAGATGACGCAACTCATCCCAGATCTCCTGGGTGTTGTTGTAGTGCATCGGATAACCCATACGGGTGGCGATTTCACTGATGATTTGCCAGTCCGTTTTCAGATCCCATTTCGGTTCAACCGCCTTGAAGAAACGCTGGAAGCCACGGTCAGCCGCAGTAAACACGCCTTCATGCTCGCCCCACGA  >  NC_000913/4297537‑4298028

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