Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	NC_000913	4,531,415	C→T	E79E (GAG→GAA)	sgcX ←	putative endoglucanase with Zn‑dependent exopeptidase domain

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	4,531,415	0	C	T	100.0%	141.7 / NA	42	E79E (GAG→GAA)	sgcX	putative endoglucanase with Zn‑dependent exopeptidase domain
Reads supporting (aligned to +/- strand):  ref base C (0/0);  new base T (23/19);  total (23/19)

GACAGGCGTACCGACCTGAATGCCCATCCGTATAGCGTCGTCTTTGTCTTTGGCACCAATATCAATCCACAGTTTGTCGACAGAAGGTGACTGCGTGCGCTCGTCGCCTTTGGCAAAGTGGTAGGACTTAATGCCGATACACCCCATGACTGGCCCTTTGTCCCCGGTGAGGGTGACGATGGAACCAGCCATAGTGACCTGCGCAGGACCGCCTACGCGTTCAAAGCGTAAAAATCCAGACGGCTCTATTTTGCGCACCATAAAACCGACTTCATCCATATGCGCAAAAATCATCAGTCGCAGCGCATCGGGTTTATCACTACCGTAGCGCGCGACGACATTCCCAAGCCTGTCTCGCCAGACCTCTTTTGCCTGACGCCTGAATTCGCACAGCATGACGTCAGCGACAGCGTTTTCGTGTCCTGAAATCGCATTGTGCTGCAAAAGCGAGAACAGGGTTTCCTGCACAGAAAATG  >  NC_000913/4531172‑4531647                                                                                                                                                                                                                                                    |                                                                                                                                                                                                                                         gACAGGCGTACCGACCTGAATGCCCATCCGTATAGCGTCGTCTTTGTCTTTGGCACCAATATCAATCCACAGTTTGTCGACAGAAGGTGACTGCGTGCGCTCGTCGCCTTTGGCAAAGTGGTAGGACTTAATGCCGATACACCCCATGACTGGCCCTTTGTCCCCGGTGAGGGTGACGATGGAACCAGCCATAGTGACCTGCGCAGGACCGCCTACGCGTTCAAAGCGTAAAAATCCAGACGGTTCTAttt                                                                                                                                                                                                                                   >  1:54907/1‑251 (MQ=255)            cGACCTGAATGCCCATCCGTATAGCGTCGTCTTTGTCTTTGGCACCAATATCAATCCACAGTTTGTCGACAGAAGGTGACTGCGTGCGCTCGTCGCCTTTGGCAAAGTGGTAGGACTTAATGCCGATACACCCCATGACTGGCCCTTTGTCCCCGGTGAGGGTGACGATGGAACCAGCCATAGTGACCTGCGCAGGACCGCCTACGCGTTCAAAGCGTAAAAATCCAGACGGTTCTATTTTGCGCACCATa                                                                                                                                                                                                                        <  1:201805/251‑1 (MQ=255)               ccTGAATGCCCATCCGTATAGCGTCGTCTTTGTCTTTGGCACCAATATCAATCCACAGTTTGTCGACAGAAGGTGACTGCGTGCGCTCGTCGCCTTTGGCAAAGTGGTAGGACTTAATGCCGATACACCCCATGACTGGCCCTTTGTCCCCGGTGAGGGTGACGATGGAACCAGCCATAGTGACCTGCGCAGGACCGCCTACGCGTTCAAAGCGTAAAAATCCAGACGGTTCTATTTTGCGCACCTTaaa                                                                                                                                                                                                                      >  2:77545/1‑250 (MQ=255)                   aaTGCCCATCCGTATAGCGTCGTCTTTGTCTTTGGCACCAATATCAATCCACAGTTTGTCGACAGAAGGTGACTGCGTGCGCTCGTCGCCTTTGGCAAAGTGGTAGGACTTAATGCCGATACACCCCATGACTGGCCCTTTGTCCCCGGTGAGGGTGACGATGGAACCAGCCATAGTGACCTGCGCAGGACCGCCTACGCGTTCAAAGCGTAAAAATCCAGACGGTTCTATTTTGCGCACCATAAAACCGa                                                                                                                                                                                                                 <  1:173485/251‑1 (MQ=255)                                ataGCGTCGTCTTTGTCTTTGGCACCAATATCAATCCACAGTTTGTCGACAGAAGGTGACTGCGTGCGCTCGTCGCCTTTGGCAAAGTGGTAGGACTTAATGCCGATACACCCCATGACTGGCCCTTTGTCCCCGGTGAGGGTGACGATGGAACCAGCCATAGTGACCTGCGCAGGACCGCCTACGCGTTCAAAGCGTAAAAATCCAGACGGTTCTATTTTGCGCACCATAAAACCGACTTCATCCATATg                                                                                                                                                                                                    >  2:154712/1‑251 (MQ=255)                                        gtctttgtctttGGCACCAATATCAATCCACAGTTTGTCGACAGAAGGTGACTGCGTGCGCTCGTCGCCTTTGGCAAAGTGGTAGGACTTAATGCCGATACACCCCATGACTGGCCCTTTGTCCCCGGTGAGGGTGACGATGGAACCAGCCATAGTGACCTGCGCAGGACCGCCTACGCGTTCAAAGCGTAAAAATCCAGACGGTTCTATTTTGCGCACCAt                                                                                                                                                                                                                         >  2:210479/1‑222 (MQ=255)                                        gtctttgtctttGGCACCAATATCAATCCACAGTTTGTCGACAGAAGGTGACTGCGTGCGCTCGTCGCCTTTGGCAAAGTGGTAGGACTTAATGCCGATACACCCCATGACTGGCCCTTTGTCCCCGGTGAGGGTGACGATGGAACCAGCCATAGTGACCTGCGCAGGACCGCCTACGCGTTCAAAGCGTAAAAATCCAGACGGTTCTATTTTGCGCACCAt                                                                                                                                                                                                                         <  1:210479/222‑1 (MQ=255)                                            ttgtctttgGCACCAATATCAATCCACAGTTTGTCGACAGAAGGTGACTGCGTGCGCTCGTCGCCTTTGGCAAAGTGGTAGGACTTAATGCCGATACACCCCATGACTGGCCCTTTGTCCCCGGTGAGGGTGACGATGGAACCAGCCATAGTGACCTGCGCAGGACCGCCTACGCGTTCAAAGCGTAAAAATCCAGACGGTTCTATTTTGCGCACCATAAAACCGACTTCATCCATATGCGCAAAAat                                                                                                                                                                                           <  1:251857/248‑1 (MQ=255)                                                                          ttGTCGACAGTAGGTGACTGCGCGCGCTCGTCGCCTTTGGCAAAGTGGTAGGACTTAATGCCGATACACCCCATGACTGGCCCTTTGTCCCCGGTGAGGGTGACGATGGAACCAGCCATAGTGACCTGCGCAGGACCGCCTACGCGTTCAAAGCGTAAAAATCCAGACGGTTCTATTTTGCGCACCATAAAACCGACTTCCTCCATATGCGCAAAAATCATCAGTCGCAGCGCATCGGGTTTATCACTAcc                                                                                                                                                          >  1:82342/1‑251 (MQ=255)                                                                              cGACAGAAGGTGACTGCGTGCGCTCGTCGCCTTTGGCAAAGTGGTAGGACTTAATGCCGATACACCCCATGACTGGCCCTTTGTCCCCGGTGAGGGTGACGATGGAACCAGCCATAGTGACCTGCGCAGGACCGCCTACGCGTTCAAAGCGTAAAAATCCAGACGGTTCTATTTTGCGCACCATAAAACCGACTTCATCCATATGCGCAAAAATCATCAGTCGCAGCGCATCGGGTTTATCACTACCGt                                                                                                                                                        >  1:263536/1‑249 (MQ=255)                                                                                  aGAAGGTGACTGCGTGCGCTCGTCGCCTTTGGCAAAGTGGTAGGACTTAATGCCGATACACCCCATGACTGGCCCTTTGTCCCCGGTGAGGGTGACGATGGAACCAGCCATAGTGACCTGCGCAGGACCGCCTACGCGTTCAAAGCGTAAAAATCCAGACGGTTCTATTTTGCGCACCATAAAACCGACTTCATCCATATGCGCAAAAATCATCAGTCGCAGCGCATCGGGTTTATCACTACCGTAgcgcg                                                                                                                                                  <  1:260478/251‑1 (MQ=255)                                                                                      ggTGACTGCGAGCGCTCGTCGCCTTTGGCAAAGTGGTAGGACTTAACGCAGATACACCCCATGACTGGCCCTTTGTCCACGGTGAGGGTGACGATGGAACCAGCCATAATGACCTGGGCAGGACCGCCTACGCGTTCAAAGCGGAACAATCCAGCCGGTTCTATCTTGAGCACCATAAAACAGACTTCATACATATGCGCCAAAATCATCAGTCGCAGCGCATCGGGTTTATCACTAACGAAACGCGcga                                                                                                                                               >  1:124934/1‑250 (MQ=255)                                                                                           cTGCGTGCGCACGTCGCATCTGGCAAAGTGGTAGGACTTAAACCCGATACACCCCATGAATGGACCTTTGCCCACGGTGAGGGTGACGATGGAACCAGCCATAGTGACCTGCGCAGGACCGCCTACGCGTTCAAAGCGTAAAAATCCAGACGGTTCTATTTTGCGCACCATAAAACCGACTTCATCCATATGCGCAAAAATCATCAGTCGCAGCGCATCGGGTTTATCACTTCCGTACCGCGCgacgacct                                                                                                                                         >  2:175597/1‑249 (MQ=255)                                                                                            tgcgtgcgCTCGTCGACTTTGGCAAAGTGGTAGGACTTAATGCCGATACACCCCATGACTGGCCCTTTGTCCCCGGTGAGGGTGACGATGGAACCCGCCATAGTGACCTGCGCAGGACCGCCTACGCGCTCAAAGCGTAAAAATCCAGACGGTTCTATTTTGCGAACCATAAAACCGCCATCATCCATATGCGCAAAAATCATCAGTCGAAGCGCATCGGGATTATCACTGACGTAGCGAGCAACGTCAtt                                                                                                                                        >  2:195356/1‑251 (MQ=255)                                                                                             gcgtgcgCTCGTCGCCTTTGGCAAAGTGGTAGGACTTAATGCCGATACACCCCATGACTGGCCCTTTGTCCCCGGTGAGGGTGACGATGGAACCAGCCATAGTGACCTGCGCAGGACCGCCTACGCGTTCAAAGCGTAAAAATCCAGACGGTTCTATTTTGCGCACCATAAAACCGACTTCATCCATATGCGCAAAAATCATCAGTCGCAGCGCATCGGGTTTATCACTACCGTAGCGCGCGACGACAtt                                                                                                                                        <  1:77545/250‑1 (MQ=255)                                                                                                  cgcTCGTCGCCTTTGGCAAAGTGGTAGGACTTAATGCCGATACACCCCATGACTGGCCCTTTGTCCCCGGTGAGGGTGACGATGGAACCAGCCATAGTGACCTGCGCAGGACCGCCTACGCGTTCAAAGCGTAAAAATCCAGACGGTTCTATTTTGCGCACCATAAAACCGACTTCATCCATATGCGCAAAAATCATCAGTCGCAGCGCATCGGGTTTATCACTACCGTAGCGCGCGACGACATTCCCa                                                                                                                                    >  1:241787/1‑249 (MQ=255)                                                                                                                 gCAAAGTGGTAGGACTTAATGCCGATACACCCCAGGACTGGCCCTTTGTCCCCGGTGAGGGTGACGATGGAACCAGCCATAGTGACCTGCGCAGGACCGCCTACGCGTTCAAAGCGTAAAAATCCAGACGGTTCTATTTTGCGCACCATAAAACCGACTTCATCCATATGCGCAAAAATCATCAGTCGCAGCGCATCGGGTTTATCACTACCGTAGCGCGCGACGACATTCCCAAGCCTGTCTCGCCAGAc                                                                                                                   <  2:66891/251‑1 (MQ=255)                                                                                                                            ggACTTAATGCCGATACACCCCATGACTGGCCCTTTGTCCCCGGTGAGGGTGACGATGGAACCAGCCATAGTGACCTGCGCAGGACCGCCTACGCGTTCAAAGCGTAAAAATCCAGACGGTTCTATTTTGCGCACCATAAAACCGACTTCATCCATATGCGCAAAAATCATCAGTCGCAGCGCATCGGGTTTATCACTACCGTAGCGCGCGACGACATTCCCAAGCCTGTCTCGCCAGACCTCTTTTGCCt                                                                                                        >  2:184276/1‑251 (MQ=255)                                                                                                                            ggACTTAATGCCGATACACCCCATGACTGGCCCTTTGTCCCCGGTGAGGGTGACGATGGAACCAGCCATAGTGACCTGCGCAGGACCGCCTACGCGTTCAAAGCGTAAAAATCCAGACGGTTCTATTTTGCGCACCATAAAACCGACTTCATCCATATGCGCAAAAATCATCAGTCGCAGCGCATCGGGTTTATCACTACCGTAGCGCGCGACGACATTCCCAAGCCTGTCTCGCCAGACCTCTTTTGCCt                                                                                                        >  1:203716/1‑251 (MQ=255)                                                                                                                             gACTTAATGCCGATACACCCCATGACTGGCCCTTTGTCCCCGGTGAGGGTGACGATGGAACCAGCCATAGTGACCTGCGCAGGACCGCCTACGCGTTCAAAGCGTAAAAATCCAGACGGTTCTATTTTGCGCACCATAAAACCGACTTCATCCATATGCGCAAAAATCATCAGTCGCAGCGCATCGGGTTTATCACTACCGTAGCGCGCGACGACATTCCCAAGCCTGTCTCGCCAGACCTCTTTTGCCTg                                                                                                       <  1:208527/251‑1 (MQ=255)                                                                                                                                         aTACACCCCATGACTGGCCCTTTGTCCCCGGTGAGGGTGACGATGGAACCAGCCATAGTGACCTGCGCAGGACCGCCTACGCGTTCAAAGCGTAAAAATCCAGACGGTTCTATTTTGCGCACCATAAAACCGACTTCATCCATATGCGCAAAAATCATCAGTCGCAGCGCATCGGGTTTATCACTACCGTAGCGCGCGACGACATTCCCAAGCCTGTCTCGCCAGACCTCTTTTGCCTGACGCCTGAATTc                                                                                           >  1:13675/1‑251 (MQ=255)                                                                                                                                          aacacCCCATGACTGGCCCTTTGTCCCCGGTGAGGGTGACGATGGAACCAGCCATAGTGACCTGCGCAGGACCGCCTACGCGTTCAAAGCGTAAAAATCCAGACGGTTCTATTTTGCGCACCATAAAACCGACTTCATCCATATGCGCAAAAATCATCAGGCGCAGCGCATCGGGTTTATCACTACCGTAGCGCGCGCCGACATTCCCAAGCCTGTCTCGCCAGACCTCTTTTGCCTGACGCCTGAATTCg                                                                                          <  2:241787/250‑1 (MQ=255)                                                                                                                                                ccattatcgcctcTTTGTCCCCGGTGAGGGTGACGATGGAACCAGCCATAGGTACCTGCTCAGGACCGCCTACGTGTTCAAAGCGTAAAAACCCAGACGGTTCTATTTTGCTCACCATCAACCCTACTTCATCCATATGCGCAAAAATCATCAGTCGCAGCGCATCGGGTTTATCACTACCGTAGCGCGCGACGACATTCCCAAGCCTGTCTCGCCAGACCTCTTTTGCCTGACGCCTGAATTCGCACAGc                                                                                    <  1:241359/239‑1 (MQ=255)                                                                                                                                                  aTGACTGGCCCTTTGTCCCCGGTGAGGGTGACGATGGAACCAGCCATAGTGACCTGCGCAGGACCGCCTACGCGTTCAAAGCGTAAAAATCCAGACGGTTCTATTTTGCGCACCATAAAACCGACTTCATCCATATGCGCAAAAATCATCAGTCGCAGCGCATCGGGTTTATCACTACCGTAGCGCGCGACGACATTCCCAAGCCTGTCTCGCCAGACCTCTTTTGCCTGACGCCTGAATTCGCACAGCAt                                                                                  >  2:178399/1‑251 (MQ=255)                                                                                                                                                    gACTGGCCCCTTATACCCGGAGAGGGTGACGATGGAACCAGCCATAGTGCCCTGCGCAGGACCGCCTACGCGTTCAAAGCGTAAAAATCCAGCCGGTTCTATTTTGCGCACCATAAAACCGACTTCCTCCATATGCGCAAAAATCATCAGTCGCAGCGCATCGGGTTTATCACTACCGTAGCGCGCGACGACATTCACAAGCCTGACTCGCCAGACCTCTTTTGCCTGGCGCCTGAATTCGCACAGCATg                                                                                 >  2:246938/1‑250 (MQ=255)                                                                                                                                                         gCCCTTTGTCCCCGGCGAGGGAGACGATGGAACCAGCCATAGTGACCTGCGCAGGACCGCCTACGCGTTCAAAGCGTAAACATCCAGACGGTTCTATTTTGCGCACCATAAAAACGACTTCATCCATATGCGCAAAAATCATCAGTCGCAGCGCATCGGGTTTATCACTACCGTAGCGCGCGACGACATTCCCAAGCCTGTCTCGCCAGACCTCTTTTGCCTGACGCCTGAATTCGCACAGCATGACGTc                                                                            >  2:260008/1‑250 (MQ=255)                                                                                                                                                          cccTTTGTCCCCGGTGAGGGTGACGATGGAACCAGCCATAGTGACCTGCGCAGGACCGCCTACGCGTTCAAAGCGTAAAAATCCAGACGGTTCTATTTTGCGCACCATAAAACCGACTTCATCCATATGCGCAAAAATCATCAGTCGCAGCGCATCGGGTTTATCACTACCGTAGCGCGCGACGACATTCCCAAGCCTGTCTCGCCAGACCTCTTTTGCCTGACGCCTGAATTCGCACAGCATGACGTc                                                                            >  2:24413/1‑249 (MQ=255)                                                                                                                                                                       gTGAGGGTGACGATGGAACCAGCCATAGTGACCTGCGCAGGACCGCCTACGCGTTCAAAGCGTAAAAATCCAGACGGTTCTATTTTGCGCACCATAAAACCGACTTCATCCATATGCGCAAAAATCATCAGTCGCAGCGCATCGGGTTTATCACTACCGTAGCGCGCGACGACATTCCCAAGCCTGTCTCGCCAGACCTCTTTTGCCTGACGCCTGAATTCGCACAGCATGACGTCAGCGACAGCGTTTTc                                                             <  1:1600/251‑1 (MQ=255)                                                                                                                                                                                   aTGGAACCAGCCATAGAGACCTGCGCAGGACCGCCTACGCGTTCAAAGCGTAAAAATCCAGACGGTTCTATTTTGCGCACCATAACACCGACTTCATCCATATGCGCAAAAATCATCAGTCGCAGCGCATCGGGTTTATCACTACCGTAGCGCGCGACGACATTCCCAAGCCTGTCTCGCCAGACCTCTTTTGCCTGACGCCTGAATTCGCACAGCATAACGTCAGCGACAGCGTTTTCGTGTCCTGAAAt                                                 >  1:210381/1‑251 (MQ=255)                                                                                                                                                                                      gAACCAGCCATAGTGACCTGCGCAGGACCGCCTACGCGTTCAAAGCGCAAAAATCCAGACGGTTCTATTTTGCGCACCATAAAACCGACTTCATCCATATGCGCAAAAATCATCAGTCGCAGCGCATCGGGTTTATCACTACCGTAGCGCGCGACGACATTCCCAAGCCTGTCTCGCCAGACCTCTTTTGCCTGACGCCTGAATTCGCACAGCATGACGTCAGCGACAGCGTTTTCGTGTCCTGAAATCGc                                              >  1:219556/1‑251 (MQ=255)                                                                                                                                                                                                  gTGACCTGCGCAGGACCGCCTACGCGTTCAAAGCGTAAAAATCCAGACGGTTCTATTTTGCGCACCATAAAACCGACTTCATCCATATGCGCAAAAATCATCAGTCGCAGCGCATCGGGTTTATCACTACCGTAGCGCTCGACGACATTCCCAAGCCTGTCTCGCCAGACCTCTTTTGCCTTACGCCTGAATTCGCACAGCATGACGTCAGCGACAGCGTTTTCGTGTCCTGAAATCGCATTTTGCTGCaa                                  <  2:126733/251‑1 (MQ=255)                                                                                                                                                                                                      ccTGCGCAGGACCGCCTACGCGTTCAAAGCGTAAAAATCCAGACGGTTCTATTTTGCGCACCATAAAACCGACTTCATCCATATGCGCAAAAATCATCAGTCGCAGCGCATCGGGTTTATCACTACCGTAGCGCGCGACGACATTCCCAAGCCTGTCTCGCCAGACCTCTTTTGCCTGACGCCTGAATTCGCACAGCATGACGTCAGCGACAGCGTTTTCGTGTCCTGAAATCGCATTGTGCTGCAAAAGc                              >  1:182764/1‑251 (MQ=255)                                                                                                                                                                                                         gcacAGGACCGCCTAAGCGTTAAAAGCTTAAAAATCCAGACGGTTCTATTTTGCGCACCATAAAACCGACTTCATCCATATGCGCAAAAATCATCAGTCGCAGCTCATCTGGTTTATCACTACCGTAGCGCGCGACGACATTCCCAAGCCTGTCTCGCCAGAACTCTTTGGCCTGCCGCCTGAATTCGCACAGCATGACGTCAGCGACAGCGTTTTCGTGTCCTGAAATCGCATTGTGCTGCAAAATCgag                           <  2:199904/248‑1 (MQ=255)                                                                                                                                                                                                           gcAGGACCGCCTACGCGTTCAAAGCGTAAAAATCCAGACGGTTCTATTTTGCGCACCATAAAACCGACTTCATCCATATGCGCAAAAATCATCAGTCGCAGCGCATCGGGTTTATCACTACCGTAGCGCGCGACGACATTCCCAAGCCTGTCTCGCCAGACCTCTTTTGCCTGACGCCTGAATTCGCACAGCATGACGTCAGCGACAGCGTTTTCGTGTCCTGAAATCGCATTGTGCTGCAAAAGCGAGaa                         <  1:190144/251‑1 (MQ=255)                                                                                                                                                                                                                aCCGCCTACGCGTTCAAAGCGTAAAAATCCAGACGGTTCTATTTTGCGCACCATAAAACCGACTTCATCCATATGCGCAAAAATCATCAGTCGCAGCGCATCGGGTTTATCACTACCGTAGCGCGCGACGACATTCCCAAGCCTGTCTCGCCAGACCTCTTTTGCCTGACGCCTGAATTCGCACAGCATGACGTCAGCGACAGCGTTTTCGTGTCCTGAAATCGCATTGTGCTGCAAAAGCGAGAACAggg                    >  1:177101/1‑251 (MQ=255)                                                                                                                                                                                                                   gCCTACGCGTTCAAAGCGTAAAAATCCAGACGGTTCTATTTTGCGCACCATAAAACCGACTTCATCCATATGCGCAAAAATCATCAGTCGCAGCGCATCGGGTTTATCACTACCGTAGCGCGCGACGTCATTCCCAAGCCTGTCTCGCCAGACCTCTTTTGCCTGACGCCTGAATTCGCACAGCATGACGTCAGCGACAGCGTTTTCGTGTCCTGAAATCGCATTGTGCTGCAAAAGCGAGAACAGGGttt                 <  1:63714/251‑1 (MQ=255)                                                                                                                                                                                                                    ccTACGCGTTCAAAGCGTAAAAATCCAGTCGGTTCTATTTTGCGCACCATAAAACCGACTTCATCCATATGCGCAAAAATCATCAGTCGCAGCGCATCGGGTTTATCACTACCGTAGCGCGCGACGACATTCCCAAGCCTGTCTCGCCAGACCTCTTTTGCCTGACGCCTGAATTCGCACAGCATGACGTCAGCGACAGCGTTTTCGTGTCCTGAAATCGCATTGTGCTGCAAAAGCGAGAACAGGGTTTc                <  2:82342/251‑1 (MQ=255)                                                                                                                                                                                                                          cgTTCAAAGCGTAAAAATCCAGACGGTTCTATTTTGCGCACCATAAAACCGACTTCATCCATATGCGCAAAAATCATCAGTCGCAGCGCATCGGGTTTATCACTACCGTAGCGCGCGACGACATTCCCAAGCCTGTCTCGCCAGACCTCTTTTGCCTGACGCCTGAATTCGCACAGCATGACGTCAGCGACAGCGTTTTCGTGTCCTGAAATCGCATTGTGCTGCAAAAGCGAGAACAGGGTTTCCTGcac          >  2:89338/1‑251 (MQ=255)                                                                                                                                                                                                                          cgTTCAAAGCGTAAAAATCCAGACGGTTCTATTTTGCGCACCATAAAACCGACTTCATCCATATGCGCAAAAATCATCAGTCGCAGCGCATCGGGTTTATCACTACCGTAGCGCGCGACGACATTCCCAAGCCTGTCTCGCCAGACCTCTTTTGCCTGACGCCTGAATTCGCACAGCATGACGTCAGCGACAGCGTTTTCGTGTCCTGAAATCGCATTGTGCTGCAAAAGCGAGAACAGGGTTTCCTGcac          >  2:55818/1‑251 (MQ=255)                                                                                                                                                                                                                             tCAAAGCGTAAAAATCCAGACGGTTCTATTTTGCGCACCATAAAACCGACTTCATCCATATGCGCAAAAATCATCAGTCGCAGCGCATCGGGTTTATCACTACCGTAGCGCGCGACGACATTCCCAAGCCTGTCTCGCCAGACCTCTTTTGCCTGACGCCTGAATTCGCACAGCATGACGTCAGCGACAGCGTTTTCGTGTCCTGAAATCGCATTGTGCTGCAAAAGCGAGAACAGGGTTTCCTGCACAg        <  1:14199/250‑1 (MQ=255)                                                                                                                                                                                                                              cAAAGCTTAAAAATCCAGACGGTTCTATTTGGCGCACCATAAAACCGACTTCATCCATATGCGCAAAAATCATCAGTCCCAGCGCATCGGGTTTATCACTACCGTAGCGCGCGACGACATTCCCAAGCCTTGCTCGCCAGACCTCTTTTGCCTGACGCCTGAATTCGCACAGCATGACGTCAGCGACAGCGTTTTCGTGTCCTGAAATCGCATTGTGCTGCAAAAGCGAGAACAGGGTTTCCTGCACAGa       <  2:161586/250‑1 (MQ=255)                                                                                                                                                                                                                                  gCGTAAAAATCCAGTCGGTTCTATTTTGCGCACCATAAAACCGACTTCATCCATATGCGCAAAAATCATCAGTCGCAGCGCATCGGGTTTATCACTACCGTATCGCGCGACGACATTCCCAAGCCTGTCTCGCCAGACCTCTTTTGCCTGACGCCTGAATTCGCACAGCATGACGTCAGCGACAGCGTTTTCGTGTCCTGAAATCGCATTGTGCTGCAAAAGCGAGAACAGGGTTTCCTGCACAGAAAATg  <  1:130633/251‑1 (MQ=255)                                                                                                                                                                                                                                                    |                                                                                                                                                                                                                                         GACAGGCGTACCGACCTGAATGCCCATCCGTATAGCGTCGTCTTTGTCTTTGGCACCAATATCAATCCACAGTTTGTCGACAGAAGGTGACTGCGTGCGCTCGTCGCCTTTGGCAAAGTGGTAGGACTTAATGCCGATACACCCCATGACTGGCCCTTTGTCCCCGGTGAGGGTGACGATGGAACCAGCCATAGTGACCTGCGCAGGACCGCCTACGCGTTCAAAGCGTAAAAATCCAGACGGCTCTATTTTGCGCACCATAAAACCGACTTCATCCATATGCGCAAAAATCATCAGTCGCAGCGCATCGGGTTTATCACTACCGTAGCGCGCGACGACATTCCCAAGCCTGTCTCGCCAGACCTCTTTTGCCTGACGCCTGAATTCGCACAGCATGACGTCAGCGACAGCGTTTTCGTGTCCTGAAATCGCATTGTGCTGCAAAAGCGAGAACAGGGTTTCCTGCACAGAAAATG  >  NC_000913/4531172‑4531647

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