Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	NC_000913	534,742	G→A	G276D (GGT→GAT)	gcl →	glyoxylate carboligase

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	534,742	0	G	A	98.3%	202.5 / NA	58	G276D (GGT→GAT)	gcl	glyoxylate carboligase
Reads supporting (aligned to +/- strand):  ref base G (0/1);  new base A (35/22);  total (35/23)
Fisher's exact test for biased strand distribution p-value = 3.97e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.90e-01

AGCTGTAGAAATGTTAATCCAGGCCGAACGTCCGGTGATTGTTGCCGGGGGCGGGGTAATTAATGCTGACGCAGCTGCACTGTTACAACAGTTTGCTGAACTGACCAGCGTTCCGGTGATCCCAACGCTAATGGGCTGGGGCTGTATCCCGGACGATCATGAACTGATGGCCGGGATGGTGGGTCTGCAAACCGCGCATCGTTACGGTAACGCAACGCTGCTGGCGTCTGACATGGTGTTTGGTATCGGTAACCGTTTTGCTAACCGTCATACCGGCTCGGTAGAGAAATACACCGAAGGGCGCAAAATCGTTCATATTGATATTGAGCCGACGCAAATTGGTCGCGTGCTGTGTCCGGATCTCGGTATTGTCTCTGATGCTAAAGCGGCGCTGACACTGCTGGTTGAAGTGGCGCAGGAGATGCAAAAAGCGGGTCGTCTGCCGTGTCGTAAAGAATGGGTCGCCGACTGCCAGCAGCGTAAACG  >  NC_000913/534500‑534985                                                                                                                                                                                                                                                   |                                                                                                                                                                                                                                                    aGCTGTAGAAATGTTAATCCAGGCCGAACGTCCGGTGATTGTTGCCGGGGGCGGGGTAATTAATGCTGACGCAGCTGCACTGTTACAACAGTTTGCTGAACTGACCAGCGTTCCGGTGATCCCATCGCTAATGGGCTGGGGCTGTATCCCGGACGATCATGAACTGATGGCCGGGATGGTGGGTCTGCAAACCGCGCATCGTTACGGTAACGCAACGCTGCTGGCGTCTGACATGGTGTTTGATATCGGt                                                                                                                                                                                                                                              >  2:404243/1‑250 (MQ=255)      tAGAAATGTTAATCCAGGCCGAACGTCCGGCGATAGTTGCCGGGGGCGGGGTACTTAATGCTGACGCACCTGCACTGTTACACCAGTTTGCTGAACAGACACGCGTCCCGGTGATCCCAACGCTAATGGGATGGGGCTGTATCCCGGACGATCAAGAACTGATGGCCGGGATGGTGAGTCTGCAAACCGCGCATCGTTACGGTAACGCACCGCTGCTGGCGTCTGACATGGTGTTTGATATCGGTAACCGt                                                                                                                                                                                                                                        >  1:544830/1‑251 (MQ=255)         aaaTGTTAATCCAGGCCGAACGTCCGGTGATTGTTGCCGGGGGCGGGGTAATTAATGCTGACGCAGCTGCACTGTTACAACAGTTTGCTGAACTGACCAGCGTTCCGGTGATCCCAACGCTAATGGGCTGGGGCTGTATCCCGGACGATCATGAACTGATGGCCGGGATGGTGGGTCTGCAAACCGCGCATCGTTACGGTAACGCAACGCTGCTGGCGTCTGACATGGTGTTTGATATCGGTAACCGtttt                                                                                                                                                                                                                                     >  1:500871/1‑251 (MQ=255)          aaTGTTAATCCAGGCCGAGCGTCCGGTGATTGTTGGCGGGGGCGGGGTAATTAATGCTGACGCAGCTGCACTGTTACAACAGTTTGCTGAACGGACCAGCGTGCCGGTGATGCCAACGCTAATGGGCTGGGGCTGTATCCCGGACGAGCATGATCTGATGGCCGGGATGGTGGGTCTGCAAACCGCGCATCGTTACGTTAACGCAACGCTGTGGGCGTCTGACATGGTGTTTGGTATCGGGCACCGtttgg                                                                                                                                                                                                                                    <  2:371430/251‑3 (MQ=255)          aaTGTTAATCCAGGCCGAACGTCCGGTGATTGTTGCCGGGGGCGGGGTAATTAATGCTGACGCAGCTGCACTGTTACAACAGTTTGCTGAACTGACCAGCGTTCCGGTGATCCCAACGCTAATGGGCTGGGGCTGTATCCCGGACGATCATGAACTGATGGCCGGGATGGTGGGTCTGCAAACCGCGCATCGTTACGGTAACGCAACGCTGCTGGCGTCTGACATGGTGTTTGATATCGGTAACCGTTTTg                                                                                                                                                                                                                                    >  1:580277/1‑251 (MQ=255)               tAATCCAGGCCGAACGTCCGGTGATTGTTGCCGGGGGCGGGGTAATTAATGCTGACGCAGCTGCACTGTTACAACAGTTTGCTGAACTGACCAGCGTTCCGGTGATCCCAACGCTAATGGGCTGGGGCTGTATCCCGGACGATCATGAACTGATGGCCGGGATGGTGGGTCTGCAAACCGCGCATCGTTACGGTAACGCAACGCTGCTGGCGTCTGACATGGTGTTTGATATCGGTAACCGTTTTGCTaa                                                                                                                                                                                                                                <  1:259429/250‑1 (MQ=255)               tAATCCAGGCCGAAAATCCGGGGATACTTGCCGGGGGCGGGGTAATTAATGCTGCCGCCCCTGCACTGATACAACAGTTTGCTGAACTGACCAGCGTTCCGGTGATACCAACGCTAATGGGCTGGGGCTGTATCCCGGACGATCATGAACTGATGGCCGGGATGGAGGGTCTGCAAACCGCGCATCGTTCCGGTCACGCAACGCTGCTGGCGTCTGACATGGTGTTTGATATCGGTAACCGATTTGCTAAc                                                                                                                                                                                                                               >  2:240989/1‑251 (MQ=255)                                ccGGTGATTGTTGCCGGGGGCGGGGTAATTAATGCTGACGCAGCTGCACTGTTACAACAGTTTGCTGAACTGACCAGCGTTCCGGTGATCCCAACGCTAATGGGCTGGGGCTGTATCCCGGACGATCATGAACTGATGGCCGGGATGGTGGGTCTGCAAACCGCGCATCGTTACGGTAACGCAACGCTGCTGGCGGCTGACATGGTGTTTGATATc                                                                                                                                                                                                                                                 <  1:413812/216‑1 (MQ=255)                                ccGGGGATTGTTGCCGGGGGCGGGGTAATTAATGCTGACGCAGCTGCACTGTTACAACAGTTTGCTGAACTGACCAGCGTTCCGGTGATCCCAACGCTAATGGGCTGGGGCTGTATCCCGGACGATCATGAACTGATGGCCGGGATGGTGGGTCTGCAAACCGCGCATCGTTACGGTAACGCAACGCTGCTGGCGTCTGACATGGTGTTTGATATc                                                                                                                                                                                                                                                 >  2:413812/1‑216 (MQ=255)                                   gTGATTGTTGCCGGGGGCGGGGAAATTAATGCTGACGCAGCTGCACTGTTACAACAGTTTGCTGAACTGACCAGCGTTCCGGTGATCCCAACGCTAATGGGCTGGGGCTGTATCCCGGACGATCATGAACTGATGGCCGGGATGGTGGGTCTGCAAACCGCGCATCGTTACGGTAACGCAACGCTGCTGGCGTCTGACATGGTGTTTGATATCGGTAACCGTTTTGCTAACCGTCATACCGGCTCGGt                                                                                                                                                                                                              >  1:507379/1‑248 (MQ=255)                                                 gggCGGGGTAATTAATGCTGACGCAGCTGCACTGTTACAACAGTTTGCTGAACTGACCAGCGTTCCGGTGATCCCAACGCTAATGGGCTGGGGCTGTATCCCGGACGATCATGAACTGATGGCCGGGATGGTGGGTCTGCAAACCGCGCATCGTTACGGTAACGCAACGCTGCTGGCGTCTGACATGGTGTTTGATATCGGTAACCGTTTTGCTAACCGTCATACCGGCTCGGTAGAGAAATACACCGa                                                                                                                                                                                               >  1:130233/1‑249 (MQ=255)                                                    cGGGGTAATTCATGCTGACGCAGCTGCACTGTTACAACAGTTTGCTGAACTGACCAGCGTTCCGGTGCTCCCAACGCTAATGGGCTGGGGCTGTATCCCGGACGATCATGAACTGATGGCCGGGATGGTGGGTCTGCAAACCGCGCATCGTTACGGTAACGCAACGCTGCTGGCGTCTGACATGGTGTTTGATATCGGTAACCGTTTTGCTAACCGTCATACCGGCTCGGTAGAGAAATACACCGAAGGgc                                                                                                                                                                                          >  2:526465/1‑251 (MQ=255)                                                         taattaatGCTGACGCAGCTGCACTGTTACAACAGTATGCTGAACTGACAAGCGTTCCGGTGATCCCAACGCTAATGGGCTGGGGCTGTATCCCGGACGATCATGAACTGATGGCCGGGATGGTGGGTCTGCAAACCGCGCATCGTTACGGTAACGCAACGCTGCTGGCGTCTGACATGGTGTTTGATATCGGTAACCGTTTTGCTAACCGTCATACCGGCTCGGTAGAGAAATACAACGAAGGGCGCaa                                                                                                                                                                                      >  2:575840/1‑250 (MQ=255)                                                              aatGCTGACGCAGCTGCACTGTTACAACAGTTTGCTGAACTGACCAGCGTTCCGGTGATCCCAACGCTAATGGGCTGGGGCTGTATCCCGGACGATCATGAACTGATGGCCGGGATGGTGGGTCTGCAAACCGCGCATCGTTACGGTAACGCAACGCTGCTGGCGTCTGACATGGTGTTTGATATCGGTAACCGTTTTGCTAACCGTCATACCGGCTCGGTAGAGAAATACACCGAAGGGCGCAAAATCGt                                                                                                                                                                                >  1:108543/1‑251 (MQ=255)                                                                  cTGACGCAGCTGCACTGTTACAACAGTTTGCTGAACTGACCAGCGTTCCGGTGATCCCAACGCTAATGGGCTGGGGCTGTATCCCGGACGATCATGAACTGATGGCCGGGATGGTGGGTCTGCAAACCGCGCATCGTTACGGTAACGCAACGCTGCTGGCGTCTGACATGGTGTTTGATATCGGTAACCGTTTTGCTAACCGTCATACCGGCTCGGTAGAGAAATACACCGAAGGGCGCAACATCGTTCat                                                                                                                                                                            >  2:286673/1‑251 (MQ=255)                                                                   tGACGCAGCTGCACTGTTACAACAGTTTGCTGAACTGACCAGCGTTCCGGTGATCCCAACGCTAATGTGCTGGGGCTGTATCCCGGACGATCATGAACTGATGGCCGGGATGGTGGGTCTGCAAACCGCGCATCGTTACGGTAACGCAACGCTGCTGGCGTCTGACATGGTGTTTGATATCGGTACCCGTTTTGCTAACCGTCATACCGGCTCGGTCGAGAAATACACCGACTGGCGCCAAATCGTTCatt                                                                                                                                                                           >  2:367833/1‑250 (MQ=255)                                                                                             ttGCTGAACTGACCAGCGTTCCGGTGATCCCAACGCTAATGGGCTGGGGCTGTATCCCGGACGATCATGAACTGATGGCCGGGATGGTGGGTCTGCAAACCGCGCATCGTTACGGTAACGCAACGCTGCTGGCGTCTGACATGGTGTTTGATATCGGTAACCGTTTTGCTAACCGTCATACCGGCTCGGTAGAGAAATACACCGAAGGGCGCAAAATCGTTCATATTGATATTGAGCCGACGCAAATTgg                                                                                                                                                  <  1:7938/250‑1 (MQ=255)                                                                                                  gAACTGACCAGCGTTCCGGTGATCCCAACGCTAATGGGCTGGGGCTGTATCCCGGACGATCATGAACTGATGGCCGGGATGGTGGGTCTGCAAACCGCGCATCGTTACGGTAACGCAACGCTGCTGGCGTCTGACATGGTGTTTGATATCGGTAACCGTTTTGCTAACCGTCATACCGGCTCGGTAGAGAAATACACCGAAGGGCGCAAAATCGTTCATATTGATATTGAGCCGACGCAAATTGGTCGCGt                                                                                                                                            >  1:254705/1‑251 (MQ=255)                                                                                                  gAACTGACCAGCGTTCCGGTGATCCCAACGCTAATGGGCTGGGGCTGTATCCCGGACGATCATGAACTGATGGCCGGGATGGTGGGTCTGCAAACCGCGCATCGTTACGGTAACGCAACGCTGCTGGCGTCTGACATGGTGTTTGATATCGGTAACCGTTTTGCTAACCGTCATACCGGCTCGGTAGAGAAATACACCGAAGGGCGCAAAATCGTTCATATTGATATTGAGCCGACGCAAAGTGGTcgc                                                                                                                                              <  1:200882/249‑1 (MQ=255)                                                                                                  gAACTGACCAGCGTTCCGGTGATCCCAACGCTAATGGGCTGGGGCTGTATCCCGGACGATCATGAACTGATGGCCGGGATGGTGGGTCTGCAAACCGCGCATCGTTACGGTAACGCAACGCTGCTGGCGTCTGACATGGTGTTTGATATCGGTAACCGTTTTGCTAACCGTCATACCGGCTCGGTAGAGAAATACAACGACGGGCGCAAAATCGGTCATCTTGATATTGAGCCGACGCAAATTGGTCGCGt                                                                                                                                            >  2:67487/1‑251 (MQ=255)                                                                                                       gACCAGCGTTCCGGTGATCCCAACGCTAATGGGCTGGGGCTGTATCCCGGACGATCATGAACTGATGGCCGGGATGGTGGGTCTGCAAACCGCGCATCGTTACGGTAACGCAACGCTGCTGGCGCCTGACATGGTGTTTGATATCGGTAACCGTTTTGCTAACCGTCATACCGGCTCGGTAGAGAAATACACCGAAGGGCGCAAAATCGTTCATATTGATATTGAGCCGACGCAAATTGGTCGCGTGCtgt                                                                                                                                       >  2:394029/1‑251 (MQ=255)                                                                                                        aCCAGCGTTCCGGTGATCCCAACGCTAATGGGCTGGGGCTGTATCCCGGACGATCATGAACTGATGGCCGGGATGGTGGGTCTGCAAACCGCGCATCGTTACGGTAACGCAACGCTGCTGGCGTCTGACATGGTGTTTGATATCGGTAACCGTTTTGCTAACCGTCATACCGGCTCGGTAGAGAAATACACCGAAGGGCGCAAAATCGTTCATATTGATATTGAGCCGACGCAAATTGGTCGCGTGCtgtg                                                                                                                                      <  1:188925/251‑1 (MQ=255)                                                                                                           aGCGTTCCGGTGATCCCAACGCTAATGGGCTGGGGCTGTATCCCGGACGATCATGAACTGATGGCCGGGATGGTGGGTCTGCAAACCGCGCATCGTTACGGTAACGCAACGCTGCTGGCGTCTGACATGGTGTTTGATATCGGTAACCGTTTTGCTAACCGTCATACCGGCTCGGTAGAGAAATACACCGAAGGGCGCAAAATCGTTCATATTGATATTGAGCCGACGCAAATTGGTCGCGTGCTGTGTcc                                                                                                                                   <  1:398181/251‑1 (MQ=255)                                                                                                                cccGGTGATCCCAACGCTAATGGGCTGGGGCTGTATCCCGGACGATCATGAACTGATGGCCGGGATGGTGGGTCTGCAAACCGCGCATCGTTACGGTAACGCAACGCTGCTGGCGTCTGACATGGTGTTTGATATCGGTAACCGTTTTGCTAACCGTCATACCGGCTCGGTAGAGAAATACACCGAAGGGCGCAAAATCGTTCATATTGATATGGAGCCGACGCAAATTGGTCGCGTGCTGTGTCCGGAtc                                                                                                                              <  1:73503/250‑1 (MQ=255)                                                                                                                 ccGGTGATCCCAACGCTAATGGGCTGGGGCTGTATCCCGGACGATCATGAACTGATGGCCGGGATGGTGGGTCTGCAAACCGCGCATCGTTACGGTAACGCAACGCTGCTGGCGTCTGACATGGTGTTTGATATCGGTAACCGTTTTGCTAACCGTCATACCGGCTCGGTAGAGAAATACACCGAAGGGCGCAAAATCGTGCATATGGATATTGAGCCGACGCAAATGGGTCGCGTGCTGTGTCCGGAtct                                                                                                                             <  1:224366/251‑1 (MQ=255)                                                                                                                 ccGGTGATCCCAACGCTAAGGGGCTGGGGCTGTATCCCGGACGATCATGAACTGATAGCGGGGATGGTGGGTCTGCAAACCGCGCATCGTTACGGTAACGCAACGCTGCTGGCGTCTGACATGGTGTTTGATATCGGTAACCGTGTTGCTAACCGTCATACCGGCTCGGTAGAGAAATACACCGAAGGGCGCAAAATCGTGCATATTGATAATGAGCCGGCGCAAATTGGTGGCGTGCTGTGTTCGGAtct                                                                                                                             <  2:75137/251‑1 (MQ=255)                                                                                                                   ggTGATCCCAACGCTAATGGGCTGGGGCTGTATCCCGGACGATCCTGAACTGCTGGCCGGGATGGTGGGTCTGCAAACTGCGCATCGTTACGGTAACGCAACGCTGCTGGCGTCTGACATGGTGTTTGATATAGGTAACCGTTTTGCTAACCGTCATACCGGCTCGGTAGAGAAATACACCGAAGGGCGCAAAATCGTTCATATTGATATTGAGACGACGCAAATTGGTCGCGTGCTGTGTCCGGATCTCg                                                                                                                           >  1:351846/1‑251 (MQ=255)                                                                                                                   ggTGATCCCAACGCTAATGGGCTGGGGCTGTATCCCGGACGATCATGAACTGATGGCCGGGATGGTGGGTCTGCAAACCGCGCATCGTTACGGTAACGCAACGCTGCTGGCGTCTGACATGGTGTTTGATATCGGTAACCGTTTTGCTAACCGTCATACCGGCTCGGTAGAGAAATACACCGAAGGGCGCAAAATCGTTCATATTGATATTGAGCCGACGCAAATTGGTCGCGTGCTGTGTCCGGATCTCg                                                                                                                           >  1:549756/1‑251 (MQ=255)                                                                                                                           cAACGCTAATGGGCTGGGGCTGTATCCCGGACGATCATGAACTGATGGCCGGGATGGTGGGTCTGCAAACCGCGCATCGTTACGGTAACGCAACGCTGCTGGCGTCTGACATGGTGTTTGATATCGGTAACCGTTTTGCTAACCGTCATACCGGCTCGGTAGAGAAATACACCGAAGGGCGCAAAATCGTTCATATTGATATTGAGCCGACGCAAATTGGTCGCGTGCTGTGTCCGGATCTCGGTATTGtc                                                                                                                   >  2:87958/1‑251 (MQ=255)                                                                                                                                                gTATCCCGGACGATCATGAACTGATGGCCGGGATGGTGGGTCTGCAAACCGCGCATCGTTACGGTAACGCAACGCTGCTGGCGTCTGACATGGTGTTTGATATCGGTAACCGTTTTGCTAACCGTCATACCGGCTCGGTAGAGAAATACACCGAAGGGCGCAAAATCGTTCATATTGATATTGAGCCGACGCAAATTGGTCGCGTGCTGTGTCCGGATCTCGGTATTGTCTCTGATGCTAAAGCGGCGCTg                                                                                              <  2:553695/251‑1 (MQ=255)                                                                                                                                                     ccGGACGATCATGAACTGATGGCCGGGATGGTGGGTCTGCAAACCGCGCATCGTTACGGTAACGCAACGCTGCTGGCGTCTGACATGGTGTTTGATATCGGTAACCGTTTTGCTAACCGTCATACCGGCTCGGTAGAGAAATACACCGAAGGGCGCAAAATCGTTCATATTGATATTGAGCCGACGCAAATTGGTCGCGTGCTGTGTCCGGATCTCGGTATTGTCTCTGATGCTAAAGCGGCGCTGACAct                                                                                         >  2:414221/1‑251 (MQ=255)                                                                                                                                                      cGGACGATCATGAACTGATGGCCGGGATGGTGGGTCTGCAAACCGCGCATCGTTACGGTAACGCAACGCTGCTGGCGTCTGACATGGTGTTTGATATCGGTAACCGTTTTGCTAACCGTCATACCGGCTCGGTAGAGAAATACACCGAAGGGCGCAAAATCGTTCATATTGATATTGAGCCGACGCAAATTGGTCGCGTGCTGTGTCCGGATCTCGGTATTGTCTCTGATGCTAAAGCGGCGCTGACActg                                                                                        <  2:159909/251‑1 (MQ=255)                                                                                                                                                       ggACGATCATGAACTGAAGGCCGGGATGGTGGGTCTGCACACCGCGCATCGTTACGGTAACGCAACGCTGCTGGCGTCTGACATGGTGTTTGATATCGGTAACCGTTTTGCTAACCGTCATACCGGCTCGGTAGAGAAATACACCGAAGGGCGCAAAATCGTTCATATTGATATTGAGCCGACGCAAATTGGTCGCGTGCTGTGTCCGGATCTCGGTATTGTCTCTGATGCTAAAGCGGCGCTGACActgc                                                                                       >  2:334492/1‑251 (MQ=255)                                                                                                                                                         aCGATCATGAACTGATGGCCGGGATGGTGGGTCTGCAAACCGCGCATCGTTACGGTAACGCAACGCTGCTGGCGTCTGACATGGTGTTTGATATCGGTAACCGTTTTGCTAACCGTCATACCGGCTCGGTAGAGAAATACACCGAAGGGCGCAAAATCGTTCATATTGATATTGAGCCGACGCAAATTGGTCGCGTGCTGTGTCCGGATCTCGGTATTGTCTCTGATGCTAAAGCGGCGCTGACActgctg                                                                                     >  1:469445/1‑251 (MQ=255)                                                                                                                                                          cGATCATGAACTGATGGCCGGGATGGTGGGTCTTCAAACCGCGCATCGTTACGGTAACGCAACGCTGCTGGCGTCTGACATGGTGTTTGATATCGGTAACCGTTTTGCTAACCGTCATTCCGGCTCGGTAGAGAAATACACCGAAGGGCGCAAAATCGTTCATATTGATAGTGAGCCGACGCAAAGTGGTCGCGTGCTGTGTCCGGATCTCGGTATTGGCTCTGATGCTAAAGCGGCGCTGACACTGCTgg                                                                                    <  2:297610/251‑1 (MQ=255)                                                                                                                                                              cATGAACTGATGGCCGGGATGGTGGGTCTGCAAACCGCGCATCGTTACGGTAACGCAACGCTTCTGGCGTCTGACATGGTGTTTGATATCGGTAACCGTTTTGCTAACCGTCATACCGTCTCGGTAGAGAAATACACCGAAGGGCGCAAAATCGTTCATATTGATATTGAGCCGACGCAAATTGGTCGCGTGCTGTGTCCGGATCTCGGTATTGTCTCTGATGCTAAAGCGGCGCTTACACTTCTGTTTg                                                                                 <  2:40614/250‑1 (MQ=255)                                                                                                                                                                  aaCTGATGGCCGGGATGGTGGGTCTGCAAACCGCGCATCGTTACGGTAACGCAACGCTGCTGGCGTCTGACATGGTGTTTGATATCGGTAACCGTTTTGCTAACCGTCATACCGGCTCGGTAGAGAAATACACCGAAGGGCGCAAAATCGTTCATATTGATATTGAGCCGACGCAAATTGGTCGCGTGCTGTGTCCGGATCTCGGTATTGTCTCTGATGCTAAAGCGGCGCTGACACTGCTGGTTGAAGTg                                                                            <  1:404243/251‑1 (MQ=255)                                                                                                                                                                             gggATGGTGGGTCTGCAAACCGCGCATCGTTACGGTAACGCAACGCTGCTGGCGTCTGACATGGTGTTTGATATCGGTAACCGTTTTGCTAACCGTCATACCGGCTCGGTAGAGAAATACACCGAAGGGCGCAAAATCGTTCATATTGATATTGAGCCGACGCAAATTGGTCGCGTGCTGTGTCCGGATCTCGGTATTGTCTCTGATGCTAAAGCGGCGCTGACACTGCTGGTTGAAGTGGCGCAgg                                                                     >  1:412862/1‑247 (MQ=255)                                                                                                                                                                                   gtggGTCTGCAAACCGCGCATCGTTACGGTAACGCAACGCTGCTGGCGTCTGACATGGTGTTTGATATCGGTAACCGTTTTGCTAACCTTCATACCGGCTCGGTAGAGAAATACACCGAAGGGCGCAAAATCGTTCATATTGATATTGAGCCTACGCAAATTGGTCGCGTGCTGTGTCCGTATCTCGGTATTGTCTCTGATGCTAAAGCGGCGCTGACACTGCTGGTTGAAGTGGCGCAGGAGATGCaaaa                                                           <  1:537323/251‑1 (MQ=255)                                                                                                                                                                                   gtggGTCTGCAAACCGCGCATCGTTACGGTAACGCAACGCTGCTGGCGTCTGACATGGTGTTTGATATCGGTAACCGTTTTGCTAACCGTCATACCGGCTCGGTAGAGAAATACACCGAAGGGCGCAAAATCGTTCATATTGATATTGAGCCGACGCAAATTGGTCGCGTGCTGTGTCCGGATCTCGGTATTGTCTCTGATGCTAAAGCGGCGCTGACACTGCTGGTTGAAGTGGCGCAGGAGATGCaaaa                                                           <  1:12223/251‑1 (MQ=255)                                                                                                                                                                                     gggTCTGGAAACCGCGCATCGTTACGGTAACGCAACGCTGCTGGCGTCTGACATGGTGTTTGATATCGGTAACCGTTTTGCTAACCGTCATACCGTCTCGGTAGAGAAATACACCGAAGGGCGCAAAATCGTTCATATTGATATTGAGCCGACGCAAATTGGTCGCGTGCTGTGTCCGGATCTCGGTATTGTCTCTGATGCTAAAGCGGCGCTGACACTGCTGGTTGAAGTGGCGCAGGATATGCaaaaa                                                          <  2:582990/250‑1 (MQ=255)                                                                                                                                                                                           gCAAACCGCGCATCGTTACGGTAACGCAACGCTGCTGGCGTCTGACATGGTGTTTGATATCGGTAACCGTTTTGCTAACCGTCATACCGGCTCGGTAGAGAAATACACCGAAGGGCGCAAAATCGTTCATATTGATATTGAGCCGACGCAAATTGGTCGCGTGCTGTGTCCGGATCTCGGTATTGTCTCTGATGCTAAAGCGGCGCTGACACTGCTGGTTGAAGTGGCGCAGGATATGCAAAAAGCgg                                                      <  1:396542/248‑1 (MQ=255)                                                                                                                                                                                            cAAACCGCGCATCGTTACGGTAACGCAACGCTGCTGGCGTCTGACATGGTGTTTGATATCGGTAACCGTTTTGCTAACCGTCATACCGGCTCGGTAGAGAAATACACCGAAGGGCGCAAAATCGTTCATATTGATATTGAGCCGACGCAAATTGGTCGCGTGCTGTGTCCGGATCTCGGTATTGTCTCTGATGCTAAAGCGGCGCTGACACTGCTGGTTGAAGTGGCGCAGGAGATGCAAAAAGCGGgtcg                                                  >  2:470559/1‑251 (MQ=255)                                                                                                                                                                                                ccGCGCATCGTTACGGTAACGCAACGCTGCTGGCGTCTGACATGGTGTTTGATATCGGTAACCGTTTTGCTAACCGTCATACCGGCTCGGTAGAGAAATACACCGAAGGGCGCAAAATCGTTCATATTGATATTGAGCCGACGCAAATTGGTCGCGTGCTGTGTCCGGATCTCGGTATTGTCTCTGATGCTAAAGCGGCGCTGACACTGCTGGTTGAAGTGGCGCAGGAGATGCAAAAAGCGGGTCGTCTg                                              >  1:479775/1‑251 (MQ=255)                                                                                                                                                                                                    gcATCGTTACGGTAACGCAACGCTGCTGGCGTCTGACATGGTGTTTGATATCGGTAACCGTTTTGCTAACCGTCATGCCGGCTCGGTAGAGAAATAAACCGAAGGGCGCAAAATCGTTCATATTGATATTGAGCCGACGCAAATTGGTCGCGTGCTGTGTCCGGATCTCGGTATTGTCTCTGATGCTAAAGCGGCGCTGACACTGCTGGTTGAAGTGGCGCAGGAGATGCAAAAAGCGGGTCGTCTGCCgt                                          >  2:293383/1‑251 (MQ=255)                                                                                                                                                                                                            aCGGTAACGCAACGCTGCTGGCGTCTGACATGGTGTTTGATATCGGTAACCGTTTTGCTAACCGTCATACCGGCTCGGTAGAGAAATACACCGAAGGGCGCAAAATAGTTCATATTGATATTGAGCCGACGCAAATTGGTCGCGTGCTGTGTCCGGATCTCGGTATTGTCTCTGATGCTAAAGCGGCGCTGACACTGCTGGTTGAAGTGGCGCAGGAGATGCAAAAAGCGGGTCGTCTGCCGTGTCGTaa                                   >  2:419839/1‑250 (MQ=255)                                                                                                                                                                                                                    gcaacgcTGCTGGCGTCTGACATGGTGTTTGATATCGGTAACCGTTTTGCTAACCGTCATACCGGCTCGGTAGAGAAATACACCGAAGGGCGCAAAATCGTTCATATTGATATTGAGCCGACGCAAATTGGTCGCGTGCTGTGTCCGGATCTCGGTATTGTCTCTGATGCTAAAGCGGCGCTGACACTGCTGGTTGAAGTGGCGCAGGAGATGCAAAAAGCGGGTCGTCTGCCGTGTCGTAAAGAATGGGt                          >  1:92291/1‑251 (MQ=255)                                                                                                                                                                                                                      aacgcTGCTGGCGTCTGACATGGTGTTTGATATCGGTAACCGTTTTGCTAACCGTCATACCGGCTCGGTAGAGAAATACACCGAAGGGCGCAAAATCGTTCATATTGATATTGAGCCGACGCAAATTGGTCGCGTGCTGTGTCCGGATCTCGGTATTGTCTCTGATGCTAAAGCGGCGCTGACACTGCTGGTTGAAGTGGCGCAGGAGATGCAAAAAGCGGGTCGTCTGCCGTGTCGTAAAGAATGGGTCg                        >  1:459409/1‑251 (MQ=255)                                                                                                                                                                                                                        cgcTGCTGGCGTCTGACATGGTGTTTGATATCGGTAACCGTTTTGCTAACCGTCATACCGGCTCGGTAGAGAAATACACCGAAGGGCGCAAAATCGTTCATATTGATATTGAGCCGACGCAAATTGGTCGCGTGCTGTGTCCGGATCTCGGTATTGTCTCTGATGCTAAAGCGGCGCTGACACTGCTGGTTGAAGTGGCGCAGGAGATGCAAAAAGCGGGTCGTCTGCCGTGTCGTAAAGAATGGGTCGcc                      >  1:398483/1‑251 (MQ=255)                                                                                                                                                                                                                         gctgctGGCGTCTGACATGGTGTTTGATATCGGTAACCGTTTTGCTAACCGTCATACCGGCTCGGTAGAGAAATACACCGAAGGGCGCAAAATCGTTCATATTGATATTGAGCCGACGCAAATTGGTCGCGTGCTGTGTCCGGATCTCGGTATTGTCTCTGATGCTAAAGCGGCGCTGACACTGCTGGTTGAAGTGGCGCAGGAGATGCAAAAAGCGGGTCGTCTGCCGTGTCGTAAAGAATGGGTCGCCg                     >  1:105425/1‑251 (MQ=255)                                                                                                                                                                                                                                gCGTCTGACATGGTGTTTGATATCGGTAACCGTTTTGCTAACCGTCATACCGGCTCGGTAGAGAAATACACCGAAGGGCGCAAAATCGTTCATATTGATATTGAGCCGACGCAAATTGGTCGCGTGCTGTGTCCGGATCTCGGTATTGTCTCTGATGCTAAAGCGGCGCTGACACTGCTGGTTGAAGTGGCGCAGGAGATGCAAAAAGCGGGTCGTCTGCCGTGTCGTAAAGAATGGGTCGCCGACTGCca              <  2:500871/251‑1 (MQ=255)                                                                                                                                                                                                                                 cGTCTGACATGGTGTTTGATATCGGTAACCGTTTTGCTAACCGTCATACCGGCTCGGTAGAGAAATACACCGAAGGGCGCAAAATCGTTCATATTGATATTGAGCCGACGCAAATTGGTCGCGTGCTGTGTCCGGATCTCGGTATTGTCTCTGATGCTAAAGCGGCGCTGACACTGCTGGTTGAAGTGGCGCAGGAGATGCAAAAAGCGGGTCGTCTGCCGTGTCGTAAAGAATGGGTCGCCGACTGCcag             <  2:108543/251‑1 (MQ=255)                                                                                                                                                                                                                                     tGACATGGTGTTTGATATCGGTAACCGTTTTGCTAACCGTCATACCGGCTCGGTAGAGAAATACACCGAAGGGCGCAAAATCGTTCCTATTGATATTGAGCCGACGCAAATTGGTCGCGTGCTGTGTCCGGATCTCGGTATTGTCTCTGATGCTAAAGCGGCGCTGACACTGCTGGTTGCAGTGGCGCAGGAGATGCAAAAAGCGGGTCGTCTGCCGTGTCGTAAAGAATGGGTCGCCGACTGCCagcagc         >  2:216036/1‑251 (MQ=255)                                                                                                                                                                                                                                     tGACATGGTGTTTGATATCGGTAACCGTTTTGCTAACCGTCATACCGGCGCGGTAGAGAAATACACCGAAGGGCGCAAAATCGTTCATATTGATATTGAGCCGACGCAAATTGGTCGCGTGCTGTGTCCGGATCTCGGTATTGTCTCTGATGCTAAAGCGGCGCTGACACTGCTGTTTGAATTGTCGCAGGAGAGGCAAAAAGCGGGTCGTCGGCCGTGTCGTAAAGAGTGGGTCGCCGACAGGcagcag          <  2:544830/250‑1 (MQ=255)                                                                                                                                                                                                                                          tGGTGTTTGATATCGGTAACCGTTTTGCTAACCGTCATACCGGCTCGGTAGAGAAATACACCGAAGGGCGCAAAATCGTTCATATTGATATTGAGCCGACGCAAATTGGTCGCGTGCTGTGTCCGGATCTCGGTATTGTCTCTGATGCTAAAGCGGCGCTGACACTGCTGGTTGAAGTGGCGCAGGAGATGCAAAAAGCGGGTCGTCTGCCGTGTCGTAAAGAATGGGTCGCCGACTGCCAGCAGCGTaa     >  1:94890/1‑250 (MQ=255)                                                                                                                                                                                                                                           ggTGTTTGATATCGGTAACCGTTTTGCTAACCGTCATACCGGCTCGGTAGAGAAATACACCGAAGGGCGCAAAATCGTTCATATTGATATTGAGCCGACGCAAATTGGTCGCGTGCTGTGTCCGGATCTCGGTATTGTCTCTGATGCTAAAGCGGCGCTGACACTGCTGGTTGAAGTGGCGCAGGAGATGCAAAAAGCGGGTCGTCTGCCGTGTCGTAAAGAATGGGTCGCCGACTTCCAGCAGCGTAAAc   <  2:275126/251‑1 (MQ=255)                                                                                                                                                                                                                                            gtgtTTGATATCGGTAACCGTTTTGCTAACCGTCATACCGGCTCGGTAGAGAAATACACCGAAGGGCGCCAAATCGTTCATATTGATATTGAGCCGACGCAAATTGGTCGCGGGCTGTGTCCGGATCTCGGTATTGTCTCTGATGCTAAAGCGGCGCTGACACTGCTGGTTGAAGTGGCGCAGGAGATGCAAAAAGCGGGTCGTCTGCCGTGTCGTAAAGAATGGGTCGCCGACTGCCAGCAGCGTAAACg  >  2:374494/1‑251 (MQ=255)                                                                                                                                                                                                                                            gtgtTTGATATCGGTAACCGTTTTGCTAACCGTCATACCGGCTCGGTAGAGAAATACACCGAAGGGCGCAAAATCGTTCATATTGATATTGAGCCGACGCAAATTGGTCGCGTGCTGTGTCCGGATCTCGGTATTGTCTCTGATGCTAAAGCGGCGCTGACACTGCTGGTTGAAGTGGCGCAGGAGATGCAAAAAGCGGGTCGTCTGCCGTGTCGTAAAGAATGGGTCGCCGACTGCCAGCAGCGTAAACg  >  1:200328/1‑251 (MQ=255)                                                                                                                                                                                                                                                   |                                                                                                                                                                                                                                                    AGCTGTAGAAATGTTAATCCAGGCCGAACGTCCGGTGATTGTTGCCGGGGGCGGGGTAATTAATGCTGACGCAGCTGCACTGTTACAACAGTTTGCTGAACTGACCAGCGTTCCGGTGATCCCAACGCTAATGGGCTGGGGCTGTATCCCGGACGATCATGAACTGATGGCCGGGATGGTGGGTCTGCAAACCGCGCATCGTTACGGTAACGCAACGCTGCTGGCGTCTGACATGGTGTTTGGTATCGGTAACCGTTTTGCTAACCGTCATACCGGCTCGGTAGAGAAATACACCGAAGGGCGCAAAATCGTTCATATTGATATTGAGCCGACGCAAATTGGTCGCGTGCTGTGTCCGGATCTCGGTATTGTCTCTGATGCTAAAGCGGCGCTGACACTGCTGGTTGAAGTGGCGCAGGAGATGCAAAAAGCGGGTCGTCTGCCGTGTCGTAAAGAATGGGTCGCCGACTGCCAGCAGCGTAAACG  >  NC_000913/534500‑534985

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