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	89.5%	88.7 / 0.8	38	G276D (GGT→GAT)	gcl	glyoxylate carboligase
Reads supporting (aligned to +/- strand):  ref base G (3/1);  new base A (17/17);  total (20/18)
Fisher's exact test for biased strand distribution p-value = 6.06e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 7.06e-01

AGCTGTAGAAATGTTAATCCAGGCCGAACGTCCGGTGATTGTTGCCGGGGGCGGGGTAATTAATGCTGACGCAGCTGCACTGTTACAACAGTTTGCTGAACTGACCAGCGTTCCGGTGATCCCAACGCTAATGGGCTGGGGCTGTATCCCGGACGATCATGAACTGATGGCCGGGATGGTGGGTCTGCAAACCGCGCATCGTTACGGTAACGCAACGCTGCTGGCGTCTGACATGGTGTTTGGTAT‑CGGTAACCGTTTTGCTAACCGTCATACCGGCTCGGTAGAGAAATACACCGAAGGGCGCAAAATCGTTCATATTGATATTGAGCCGACGCAAATTGGTCGCGTGCTGTGTCCGGATCTCGGTATTGTCTCTGATGCTAAAGCGGCGCTGACACTGCTGGTTGAAGTGGCGCAGGAGATGCAAAAAGCGGGTCGTCTGCCGTGTCGTAAAGAATGGGTCGCCGACTGCCAGCAGCGTAAACGCA  >  NC_000913/534500‑534987                                                                                                                                                                                                                                                   |                                                                                                                                                                                                                                                       cgCTGAAGAAATGTTAATCCAGGCCGAACGTCCGGTGATTGTTGCCGGGGGCGGGGGACTTAATGCTGACGAAGCCGCACTGTTACAACAGTTTGCTGAACTGACCAGCGTTCCGGTGATCCCAACGCTACTGGGCTGGGGCTGTATCCCGGACGATCATGAACCGATGACCGGGATGGTGGGGCTGCAAATCGCGCATCGTTACGGTAACGCGACGCTACGGGGGTCTGAACTGGAGTTTGATAT‑CGGt                                                                                                                                                                                                                                                >  2:96791/2‑250 (MQ=255)           aTGTTAACCCCGTCCGAACGGCCGGTGATTGTTGCCGGGGGCGGGGTAATTAATGCTGACGCAGCTGCACTGTTACAACGGTTTGCTGACCTGACCAGCGTTCCGGTGATCCCAACGCTAATGGGCTGGGGCTGTATCCCGGACGATCATGAACTGATGGCCGGGATGGTGGGTCTGCAAACCGCGCATCGTTACGGTAACGCAACGCGGCTGGCGTCTGACATGGTGTTTGATAT‑CGGTATCCGTTTTGc                                                                                                                                                                                                                                     <  2:360841/251‑1 (MQ=255)                   ccAGGCCGAACGTCCGGTGATTGTTGCCGGGGGCGGGGTAATTAATGCTGACGCAGCTGCACGGGTACAACAGTTTGCTGAACTGACCAGCGTTCCGGTGAGCCCACAGCTAATGGGCTGGGGCTGTATCCCGGACGAGCATGAACTGATGGCCGGGGTGGTGGGCCTGCAAACCGCGCATCTTTACGGCAACACAACGCTGCTGGCGTCTGACATGGTGTT‑GGTTTACGGTAACCGTTTTGCTCACCGTc                                                                                                                                                                                                                             >  2:74549/1‑251 (MQ=255)                      ggCCGAACGTCCGGTGATTGTTGCCGGGGGCGGGGTAATTAATGCTGACGCAGCTGCACTGTTACAACAGTTTGCTGAACCGACCAGCGTTCCGGTGATCCCAACGCTAATGGGCTGGGGCTGTATCCCGGACGATCATGAACTGATGGACGGGATGGTGGGTCAGCAAACCGCGAATCGTTACGGAAACGAAACGCTGCTGGGGCCTGACCGGGTGTTTGATTT‑CCGTAACCGTTTTGCTACACGTCAt                                                                                                                                                                                                                           >  2:307386/1‑250 (MQ=255)                         cGAACGTCCGGTGATTGTTGCCGGGGGCGAGGAAATTAATGCTGACGCAGCTGCACTGTTACAACAGGTTGCTGAACTGACCAGCGTTCCGGTGATCCCAACGCTAATGGGCTGGGGCTGTATCCCGGACGATCATGCACTGATGGCCGGGATGGGGGGTCTGCAAACCGCGCATCGTTACGGTAACGCAACGCTGCTGGCGTCTGACATGGTGTTTGATAT‑CGGTAACCGTTTTGCTAACCGTcctcacg                                                                                                                                                                                                                       >  2:79867/1‑245 (MQ=255)                          gAACGTCCGGTGATTGTTGCCGGGGGCGGGGTAATTAATGCTGACGCAGCTGCACTGTTACAACAGTTTGCTGAACTGACCAGCGTTCCGGTGATCCCAACGCTAATGGGCTGGGGCTGTATCCCGGACGATCATGAACTGATGGCCGGGATGGTGGGTCTGCAAACCGCGCATCGTTACGGTAACGCAACGCTGCTGGCGTCTGACATGGTGTTTGATAT‑CGGTAACCGTTTTGCTAACCGTCATACCgg                                                                                                                                                                                                                      <  1:238116/251‑1 (MQ=255)                             cGTCCGGTGATTGTTGCCGGGGGCGGGGTAATTAATGCTGACGCAGCTGCACTGTTACAACAGTTTGCTGAACTGACCAGCGTTCCGGTGATCCCAACGCTAATGGGCTGGTGCTGTATCCCGGACGATCATGAACTGATGGCCGGGATGGTGGGTCTGCAAACCGCGCATCGTCACGGTAACGCAACGCTGCTGGCGTCTGACATGGTGTTTGATAT‑CGGTAACCGTTTTGCTAACCGTCATACCGGCTc                                                                                                                                                                                                                   <  1:18628/251‑1 (MQ=255)                                                                tGCTGACGCAGCTGCACTGTTACAACAGTTTGCTGAACTGACCAGCGTTCCGGTGATCCCAACGCTAATGGGCTGGGGCTGTATCCCGGACGATCATGAACTGATGGCCGGGATGGTGGGTCTGCAAACCGCGCATCGTTACGGTAACGCAACGCTGCTGGCGTCTGACATGGTGTTTGATAT‑CGGTAACCGTTTTGCTAACCGTCATACCGGCTCGGTAGAGAAATACACCGAAGGGCGCAAAATCGTTc                                                                                                                                                                                >  1:82770/1‑251 (MQ=255)                                                                                 tGTTACAACAGTTTGCTGAACTGACCAGCGTTCCGGTGATCCCAACGCTAATGGGCTGGGGCTGTATCCCGGACGATCATGAACTGATGGCCGGGATGGTGGGTCTGCAAACCGCGCATCGTTACGGTAACGCAACGCTGCTGGCGTCTGACATGGTGTTTGATAT‑CGGTAACCGTTTTGCTAACCGTCATACCGGCTCGGTAGAGAAATACACCGAAGGGCGCAAAATCGTTCATATTGATATTGAACCg                                                                                                                                                               >  1:316860/1‑251 (MQ=255)                                                                                           gTTTGCTGAACTGACCAGCGTTCCGGTGATCCCAACGCTAATGGGCTGGGGCTGTATCCCGGACGATCATGAACTGATGGCCGGGATGGTGGGTCTGCAAACCGCGCATCGTTACGGTAACGCAACGCTGCTGGCGTCTGACATGGTGTTTGATAT‑CGGTAACCGTTTTGCTAACCGTCATACCGGCTCGGTAGAGAAATACACCGAAGGGCGCAAAATCGTTCATATTGATATTGAGCCGACGCAAATTg                                                                                                                                                     >  1:320794/1‑251 (MQ=255)                                                                                             ttGCTGAACTGACCAGCGTTCCGGTGATCCCAACGCGAATGGGCTGGGGCTGTATCCCGGACGATCATGAACTGATGGCCGGGATGGTGGGTCTGCAAACCGCGCATCGTTACGGTAACGCAACGCTGCTGGCGTCTGACATGGTGTTTGGTAT‑CGGTAACCGTTTTGCTAACCGTCATACCGGCTCGGTAGAGAAATACACCGAAGGGCGCAAAAGCGTTCATATTGATATTGAGCCGACGCAAATTgg                                                                                                                                                    <  2:128846/250‑1 (MQ=255)                                                                                              ggctttactgaccagtttaccagTGATCCCAACGCGGATGGGCTGGGTCTGTATCTCGGACGATCAGGGACTGAGGGCCGTGAGGGTGGGTCTGCAAACCGCGCATCGTTACGTTAACGCAACGCTGCGGGCGTCTGACATGGTGTTTGATAT‑CGGTAACCGATTTGCTAACCGGCATACCTGCTCGGTAGAGAAATACACCGAAGGGCGCAAAATCGTTCATATTGATATTGAGACGACTCAAATTGGTc                                                                                                                                                  <  2:316860/229‑1 (MQ=255)                                                                                                      tGACCAGCGTTCCGGTGATCCCAACGCTAATGGGCTGGGGCTGTATCCCGGACGATCATGAACTGATGGCCGGGATGGTGGGTCTGCAAACCGCGCATCGTTACGGTAACGCAACGCTGCTGGCGTCTGACATGGTGTTTGATAT‑CGGTAACCGTTTTGCTAACCGTCATACCGGCTCGGTAGAGAAATACACCg                                                                                                                                                                                                  <  1:347036/195‑1 (MQ=255)                                                                                                        aCCAGCGTTCCGGTGATCCCAACGCTAATGGGCTGGGGCTGTATCCCGGACGATCATGAACTGATGGCCGGGATGGTGGGTCTGCAAACCGCGCATCGTTACGGTAACGCAACGCTGATGGCGTCTGACATGGTGTTTGATAT‑CGGTAACCGTTTTGCTAACCGTCATACCGGCTCGGTAGAGAAATACACCGAAGGGCGCAAACTCGTTCATATTGATATTGAGCCGACGCAAATTGGTCGCATGCtgtg                                                                                                                                        >  1:372867/1‑251 (MQ=255)                                                                                                           aGCGTTCCGGTGATCCCAACGCTAATGGGCTGGGGCTGTATCCCGGACGATCATGAACTGATGGCCGGGATGGTGGGTCTGCAAACCGCGCATCGTTACGGTAACGCACCGCTGCTGGCGTCTGACATGGTGTTTGATAT‑CGGTAACCGTTTTGCTAACCGTCATACCGGCTCGGTAGAGAAATACACCGAAGGGCGCAAAATCGTTCATATTGATATTGAGCCGACGCAAATTGGTCGCGTGCTGTGTcc                                                                                                                                     >  1:364837/1‑251 (MQ=255)                                                                                                                      gATCCCAACGCTAATGGGCTGGGGCTGTATCCCGGACGATCATGAACTGATGGCCGGGATGGTGGGTCTGCAAACCGCGCATCGTTACGGTAACGCAACGCTGCTGGCGTCTGACATGGTGTTTGATAT‑CGGTAACCGTTTTGCTAACCGTCATACCGGCTCGGTAGAGAAATACACCGAAGGGCGCAAAATCGTTCATATTGATATTGAGCCGACGCAAATTGGTCGCGTGCTGTGTCCGGATCTCGGTa                                                                                                                          <  1:278902/251‑1 (MQ=255)                                                                                                                       aTCCCAACGCTAATGGGCTGGGGCTGTATCCCGGACGATCATGAACTGATGGCCGGGATGGTGGGTCTGCAAACCGCGCATCGTTACGGTAACGCAACGCTGCTGGCGTCTGACATGGTGTTTGATAT‑CGGTAACCGTTTTGCTAACCGTCATACCGGCTCGGTAGAGAAATACACCGAAGGGCGCAAAATCGTTCATATTGATATTGAGCCGACGCAAATTGGTCGCGTGCtgtg                                                                                                                                        <  1:270979/236‑1 (MQ=255)                                                                                                                       aTCCCAACGCTAATGGGCTGGGGCTGTATCCCGGACGATCATGAACTGAAGGCCGGGATGGTGGGTCTGCAAACCGCGCATCGTTACGGTAACGCAACGCTGCTGGCGTCTGACATGGTGTTTGATAT‑CGGTAACCGTTTTGCTAACCGTCATACCGGCTCGGGAGAGAAATACACCCAAGGGGGCAAAATCGTTCATATTGATATTGAGCCGACGCAAACTGGTCGCTTGCtgtgagatcggaagagtgg                                                                                                                         >  2:270979/1‑236 (MQ=255)                                                                                                                         cccAACGCTAATGGGCTGTGGCTGTATCCCGGACGATCATGAACTGATGGCCGGGATGGTGGGTCTGCAAACCGCGCATCGTTACGGTAACGCAACGCTGCTGGCGTCTGACATGGTGTTTGATAT‑CGGTAACCGTTTTGCTAACCGTCATACCGGCTCGGTAGAGAAATACACCGAAGGGCGCAAAATCGTTCATATTGATATTGAGCCGACGCAAATTGGTCGCGTGCTGTGTCCGGATCTCGGTATTg                                                                                                                       <  1:82401/251‑1 (MQ=255)                                                                                                                                  aaTGGGCTGTGGCTGTATCCCGGACGATCATGAACTGATGGCCGGGATGGTGGGTCTGCAAACCGCGCATCGTTACGGTAACGCAACGCTGCTGGCGTCTGACATGGTGTTTGATAT‑CGGTAACCGTTTTGCTAACCGTCATACCGGCTCGGTAGAGAAATACACCGAAGGGCGCAAAATCGTTCATATTGATATTGAGCCGACGCAAATTGGTCGCGTGCTGTGTCCGGATCTCGGTATTGTCTCTGAt                                                                                                               <  1:56429/250‑1 (MQ=255)                                                                                                                                   aTGGGCTGGGGCTGTATCCCGGACGATCATGAACTGATGGCCGGGATGGTGGGTCTGCAAACCGCGCATCGATACGGTAACGCAACGCTGCTGGCGTCTGACATGGTGTTTGATAT‑CGGTAACCGTTTTGCTAACCGTCATACCGGCTCGGTAGAGAAATACACCGAAGGGCGCAAAATCGTTCATATTGATATTGAGCCGACGCAAATTGGTCGCGTGCTGTGTACGGATCTCGGTATTGTATCTGATGc                                                                                                             >  1:154002/1‑251 (MQ=255)                                                                                                                                              ctgTATCCCGGACGATCATGAACTGCTGGCCGGGATGGTGGGTCTGCAAACCGCGCATCGTTACGGTAACGCAACGCTGCTGGCGTCTGACATGGTGTTTGATAT‑CGGTAACCGTTTTGCTAACCGTCATACCGGCTCGGTAGAGAAATACACCGAAGGGCGCAAAATCGTTCATATTGATATTGAGCCGACGCAAATTGGTCGCGTGCTGTGTCCGGATCTCGGTATTGTCTCTGATGCTAAAGCGgcgc                                                                                                  >  1:254231/1‑251 (MQ=255)                                                                                                                                                         aCGATCATGAACTGATGGACGGGATGGTGGGTCTGCAAACCGCGCATCGTTACGGTAACGCAACGCTGCTGGCGTCTGACATGGTGTTTGATAT‑CGGTAACCGTTTTGCTAACCGTCATACCGGCTCGGTAGAGAAATACACCGAAGGGCGCAAAATCGTTCATATTGATATTGAACCGACGCAAATTGGTCGCGTGCTGTGTCCGGAACTCGGTATTGACTCTGATGCTAAAGCGGCGATGacaatcctg                                                                                       >  1:277430/1‑245 (MQ=255)                                                                                                                                                           gATCATGAACTGATGGCCGGGATGGCGGGTCTGCAAACCGCGCATCGTTAAGGTAACGCAACGCTGCTGGCGTCTGACATGGTGTTTGATAT‑CGGTAACCGTTTTGCTAACCGTCATACCGGCTCGGTAGAGAAATACACCGAAGGGCGCAAAATCGTTCATATTGATATTGAGCCGACGCAAATTGGTCGCGTGCTGTGTCCGGATCTCGGTATTGTCTCTGATGCTAAAGCGGCGCTGACACTGCTGGt                                                                                     >  1:279351/1‑251 (MQ=255)                                                                                                                                                            aTCATGAACTGATGGCCGGGATGGTGGGTCTGCAAACCGCGCATCGTTACGGTAACGCAACGCTGCTGGTGTCTGACATGGTGTTTGATAT‑CGGTAACCGTTTTGCTAACCGTCATACCGGCTCGGTAGAGAAATACACCGAAGGGCGCAAAATCGTTCATATTGATATTGAGCCGACGCAAATTGGTCGCGTGCTGTGTCCGGATCTCGGTATTGTCTCTGATGCTAAAGCGGCGCTGACACTGCTGGtt                                                                                    <  1:202991/251‑1 (MQ=255)                                                                                                                                                                   acggATGGCCGGGATGGTGGGTCTGCAAACCGCGCATCGTTACGGTAACGCAACGCTGCTGGCGTCTGACATGGTGTTTGATAT‑CGGTAACCGTTTTGCTAACCGTCATACCGGCTCGGTAGAGAAATACACCGAAGGGCGCAAAATCGTTCATATTGATATTGAGCCGACGCAAATTGGTCGCGTGCTGTGTCCGGATCTCGGTATTGTCTCTGATGCTAAAGCGGCGCTGACACTGCTGGTTGAAGTgg                                                                             <  1:56019/248‑1 (MQ=255)                                                                                                                                                                             gggATGGTGGGTCTGCAAACCGCGCATCGTTACGGTAACGCAACGCTGCTGGCGTCTGACATGGTGTTTGATAT‑CGGTAACCGTTTTGCTAACCGTCATACCGGCTCGGTAGAGAAATACACCGAAGGGCGCAAAATCGTTCATATTGATATTGAGCCGACGCAAATTGGTCGCGTGCTGTGTCCGGATCTCGGTATTGTCTCTGATGCTAAAGCGGCGCTGACACTGCTGGTTGAAGTGGCGCAgg                                                                       >  1:220678/1‑247 (MQ=255)                                                                                                                                                                               gATGGTGGGTCTGCAAACCGCGCATCGTTACGGTAACGCAACGCTGCTGGCGTCTGACATGGTGTTTGATAT‑CGGTAACCGTTTTGCTAACCGTCATACCGGCTCGGTAGAGAAATACACCGAAGGGCGCAAAATCGTTCATATTGATATTGAGCCGACGCAAATTGGTCGCGTGCTGTGTCCGGATCTCGGTATTGTCTCTGATGCTAAAGCGGCGCTGACACTGCTGGTTGAAGTGGCGCAGGAGATGc                                                                 >  1:219667/1‑251 (MQ=255)                                                                                                                                                                                  ggtggGTCTGCAAACCGCGCATCGTTACGGTAACGCAACGCTGCTGGCGTCTGACATGGTGTTTGATAT‑CGGTAACCGTTTTGCTAACCGTCATACCGGCTCGGTAGAGAAATACACCGAAGGGCGCAAAATCGTTCATATTGATATTGAGCCGACGCAAATTGGTCGCGTGCTGTGTCCGGATCTCGGTATTGTCTCTGATGCTAAAGCGGCGCTGACACTGCTGGTTGAAGTGGCGCAGGAGATGCaaa                                                              <  1:10181/251‑1 (MQ=255)                                                                                                                                                                                      ggTCTGCAAACCGCGCATCGTTACGGTAACTCAACGCTGCTTGCGTCTGACATGGTGTTTGAGAT‑CGGTAACCGTTTTGCTATCCGTCATACCGGCTCGGTAGAGAAATACACCGAAGGGCGCAAAATCGTTCATATTGATATTGAGCCTACGCAAATTGGTCGCGTGCTGTGTCCGGATCTGGGTATTGTCTCTGATGCTAAAGCGGCGCTGACACTGCTGGTTGAAGTGTCGCAGGAGATGCAAAAAGc                                                          <  2:219667/251‑1 (MQ=255)                                                                                                                                                                                             aaaCCGCGCATCGTTACGGTAACGCAACGCTGCTGGCGTCTGACATGGTGTTTGATAT‑CGGTAACCGTTTTGCTAACCGTCATACCGGCTCGGTAGAGAAATACACCGAAGGGCGCAAAATCGTTCATATTGATATTGAGCCGACGCAAATTGGTCGCGTGCTGTGTCCGGATCTCGGTATTGTCTCTGATGCTAAAGCGGCGCTGACACTGCTGGTTGAAGTGGCGCAGGAGATGCAAAAAGCGGgtcgt                                                   <  1:79867/251‑1 (MQ=255)                                                                                                                                                                                                    gcATCGTTACGGTAACGCAACGCTTCTGGCGTCTGACATGGTGTTTGATAT‑CGGTAACCGTTTTGCTAACCGTCATACCGGCTCGGTAGAGAAATACACCGAAGGGCGCAAAATCGTTCATATTGATATTGAGCCGACGCAAATTGGTCGCGTGCTGTGTCCGGATCTCGGTATTGTCTCTGATGCTAAAGCGGCGCTGACACTGCTGGTTGAAGTGGCGCAGGAGATGCAAAAAGCGGGTCGTCTGCCgt                                            <  1:258861/251‑1 (MQ=255)                                                                                                                                                                                                                               ggCGTCTGACATGGTGTTTGGTAT‑CGGTAACCGTTTTGCTAACCGTCATACCGGCTCGGTAGAGAAATACACCGAAGGGCGCAAAATCGTTCATATTGATATTGAGCCGACGCAAATTGGTCGCTTGCTGTGTCCGGATCTCGGTATTGTCTCTGATGCTAAAGCGGCGCTGACACTGCTGGTTGAAGTGGCGCAGGAGATGCAAAAAGCGGGTCGTCTGCCGTGTCGTAAAGAATGGGTCGCCGACTGcc                 >  1:177179/1‑251 (MQ=255)                                                                                                                                                                                                                                  gTCTGACATGGTGTTTGATAT‑CGGTAACCGTTTTGCTAACCGTCATACCGGCTCGGTAGAGAAATACACCGAAGGGCGCAAAATCGTTCATATTGATATTGAGCCGACGCAAATTGGTCGCGTGCTGTGTCCGGATCTCGGTATTGTCTCTGATGCTAAAGCGGCGCTGACACTGCTGGTTGAAGTGGCGCAGGAGATGCAAAAAGCGGGTCGTCTGCCGTGTCGTAAAGAATGGGTCGcc                        <  1:350972/241‑1 (MQ=255)                                                                                                                                                                                                                                  gTCTGACATGGTGTTTGATAT‑CGGTAACCGTTTTGCTAAACGTCATACCGGCTCGGTAGAGAAATACACCGCAGGGCGCAAAATCGTTCATATTGATATTGAGCCGAAGAAAATTGGTCGCGTGCTGTGTCCGGATCTCGGTATTGTCTCTGATGCTAAAGCCGCGCTGACACTGCTGGTTGAAGTGGCGCAGGAGATGCAAAAAGCGGGTCGTCTGCCGTGTCGTAAAGAATGGGTCGcc                        >  2:350972/1‑241 (MQ=255)                                                                                                                                                                                                                                    cTGACATGGTGTTTGGTAT‑CGGTAACCGTTTTGCTAACCGACATACCGGCTCGGTAGAGAAATACACCCAAGGGCGACAAATCGTTCATATTGATATTGAGCCGAAGCAAATTGGTCGCGTGCTGTGTCCGGATCTCGGTATTGTCTCTCATGCTAAAGCGGCGCTGATACTGCTGGTTGAAGTGGCGCAGGGGCTGCATAAAACGGGGCGACTGCCAGGTCGGTAAGAATggggagccgcccgacaaccg            >  2:170625/1‑234 (MQ=255)                                                                                                                                                                                                                                        cATGGTGTTTGATAT‑CGGTAACCGTTTTGCTAACCGTCATACCGGCTCGGTAGAGAAATACACCGAAGGGCGCAAAATCGTTCATATTGATATTGAGCCGACGCAAATTGGTCGCGTGCTGTGTCCGGATCTCGGTATTGTCTCTGATGCTAAAGCGGCGCTGACACTGCTGGTTGAAGTGGCGCAGGAGATGCAAAAAGCGGGTCGTCTGCCGTGTCGTAAAGAATGGGTCGCCGACTGCAAGCAGCGt         >  1:343787/1‑250 (MQ=255)                                                                                                                                                                                                                                              gtTTGATAT‑CGGTAACCGTTTTGCTAACCGTCATACCGGCTCGGTAGAGAAATACACCGAAGGGCGCAAAATCGTTCATATTGATATTGAGCCGACGCAAATTGGTCGCGTGCTGTGTCCGGATCTCGGTATTGTCTCTGATGCTAAAGCGGCGCTGACACTGCTGGTTGAAGTGGCGCAGGAGATGCAAAAAGCGGGTCGTCTGCCGTGTCGTAAAGAATGGGTCGCCGACTGCCAGCAGCGTAAACGCa  <  1:48570/251‑1 (MQ=255)                                                                                                                                                                                                                                                   |                                                                                                                                                                                                                                                       AGCTGTAGAAATGTTAATCCAGGCCGAACGTCCGGTGATTGTTGCCGGGGGCGGGGTAATTAATGCTGACGCAGCTGCACTGTTACAACAGTTTGCTGAACTGACCAGCGTTCCGGTGATCCCAACGCTAATGGGCTGGGGCTGTATCCCGGACGATCATGAACTGATGGCCGGGATGGTGGGTCTGCAAACCGCGCATCGTTACGGTAACGCAACGCTGCTGGCGTCTGACATGGTGTTTGGTAT‑CGGTAACCGTTTTGCTAACCGTCATACCGGCTCGGTAGAGAAATACACCGAAGGGCGCAAAATCGTTCATATTGATATTGAGCCGACGCAAATTGGTCGCGTGCTGTGTCCGGATCTCGGTATTGTCTCTGATGCTAAAGCGGCGCTGACACTGCTGGTTGAAGTGGCGCAGGAGATGCAAAAAGCGGGTCGTCTGCCGTGTCGTAAAGAATGGGTCGCCGACTGCCAGCAGCGTAAACGCA  >  NC_000913/534500‑534987

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