Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	NC_000913	2,355,848	C→T	P110S (CCA→TCA)	glpC →	anaerobic glycerol‑3‑phosphate dehydrogenase subunit C

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	2,355,848	0	C	T	100.0%	118.5 / NA	40	P110S (CCA→TCA)	glpC	anaerobic glycerol‑3‑phosphate dehydrogenase subunit C
Reads supporting (aligned to +/- strand):  ref base C (0/0);  new base T (12/28);  total (12/28)

ATCCAGGGCCAAAACAAGCCGGGCCGGATGGCGAGCGTCTGCGTTTGAAAGATGGCGCACTGTATGACGAGGCGCTGAAATATTGCATCAACTGCAAACGTTGTGAAGTCGCCTGCCCGTCCGATGTGAAGATTGGCGATATTATCCAGCGCGCGCGGGCGAAATATGACACCACGCGCCCGTCGCTGCGTAATTTTGTGTTGAGTCATACCGACCTGATGGGTAGCGTTTCCACGCCGTTCGCACCAATCGTCAACACCGCTACCTCGCTGAAACCGGTGCGGCAGCTGCTTGATGCGGCGTTAAAAATCGATCATCGCCGCACGCTACCGAAATACTCCTTCGGCACGTTCCGTCGCTGGTATCGCAGCGTGGCGGCTCAGCAAGCACAATATAAAGACCAGGTCGCTTTCTTTCACGGCTGCTTCGTTAACTACAACCATCCGCAGTTAGGTAAAGATTTAATTAAAGTGCTCAACGCAATGGGT  >  NC_000913/2355603‑2356090                                                                                                                                                                                                                                                      |                                                                                                                                                                                                                                                   aTCCACGGCCCAAACAAGCCGGGCCGGCTAGCGAGCGTCTGCGCTAGCAAGATGGCGCACTGCATGACGAAGCGCTGAAATATTGCATCAACTGCAAACGTTGTGAAGTCGCCTGCCCGTCCGATGTGAAGATTGGCGATATTATCCAGCGCGCGCGGGCGAAATATGACACCACGCGCCCGTCGCCACGTAATTTTGTGTTGAGTCATACCGACCTGATGGGTAGCGTTTCCACGCCGTTCGCATCAATc                                                                                                                                                                                                                                               >  2:176572/1‑251 (MQ=255)        gCCAAAACAAGCCGGGCCGGATGGCGAGCGTCTGCGTTTGAAAGATGGCGCACTGTATGACGAGGCGCTGAAATATTGCATCAACTGCAAACGTTGTGAAGTCGCCTGCCCGTCCGATGTGAAGATTGGCGATATTATCCAGCGCGCGCGGGCGAAATATGACACCACGCGCCCGTCGCTGCGTAATTTTGTGTTGAGTCATACCGACCTGATGGGTAGCGTTTCCACGCCGTTCGCATCAATCGTCAaca                                                                                                                                                                                                                                        <  2:74800/251‑1 (MQ=255)                aaGCCGGGCCGGATGGCGAGCGTCTGCGTTTGAAAGATGGCGCACTGTATGACGAGGCGCTGAAATATTGCATCAACTGCAAACGTTGTGAAGTCGCCTGCCCGTCCGATGTGAAGATTGGCGATATTATCCAGCGCGCGCGGGCGAAATATGACACCACGCGCCCGTCGCTGCGTAATTTTGTGTTGAGTCATACCGACCTGATGGGTAGCGTTTCCACGCCGTTCGCATCAATCGTCAACACCGCTAcc                                                                                                                                                                                                                                >  1:117533/1‑251 (MQ=255)                           gATGGCGAGCGTCTGCGTTTGAAAGATGGCGCACTGTATGACGAGGCGCTGAAATATTGCATCAACTGCAAACGTTGTGAAGTCGCCTGCCCGTCCGATGTGAAGATTGGCGATATTATCCAGCGCGCGCGGGCGAAATATGACACCACGCGCCCGTCGCTGCGTAATTTTGTGTTGAGTCATACCGACCTGATGGGTAGCGTTTCCACGCCGTTCGCGTCAATCGTCAACACCGCTACCTCGCTGAAAcc                                                                                                                                                                                                                     <  2:94550/251‑1 (MQ=255)                                  aGCGTCTGCGTTTGAAAGATGGCGCACTGTATGACGAGGCGCTGAAATATTGCATCAACTGCAAACGTTGTGAAGTCGCCTGCCCGTCCGATGTGAAGATTGGCGATATTATCCAGCGCGCGCGGGCGAAATATGACACCACGCGCCCGTCGCTGCGTAATTTTGTGTTGAGTCATACCGACCTGATGGGTAGCGTTTCCACGCCGTTCGCATCAATCGTCAACACCGCTACCTCGCTGAAACCGGTGCgg                                                                                                                                                                                                              <  1:171359/251‑1 (MQ=255)                                         gctgttGAAAGATGGCGCACTGTATGACGAGGCGCTGAAATATTGCATCAACTGCAAACGTTGTGAAGTCGCCTGCCCGTCCGATGTGAAGATTGGCGATATTATCCAGCGCGCGCGGGCGAAATATGACACCACGCGCCCGTCGCTGCGTAATTTTGTGTTGAGTCATACCGCCCTGATGGGTAGCGTTTCCACGCCGTTCGCATCAATCGTCAACACCGCTACCTCGCTGAAACCGGTGCGGCAgcg                                                                                                                                                                                                         <  2:350115/245‑2 (MQ=255)                                         gCGTTTGAAAGATGGCGCACTGTATGCCGAGGCGCTGAAATATTGCATCAACTGCAAACGTTGTGAAGTCGCCTGCCCGTCCGATGTGAAGATTGGCGATATTATCCAGCGCGCGCGGGCGAAATATGACACCACGCGCCCGTCGCTGCGTAATTTTGTGTTGAGTCATACCGACCTGATGGGTAGCGTTTCCACGCCGTTCGCATCAATCGTCAACACCGCTACCTCGCTGAAACCGGTGCGGCAgc                                                                                                                                                                                                          <  2:304139/248‑1 (MQ=255)                                             ttGAAAGATGGCGCACTGTATGACGAGGCGCTGAAATATTGCATCAACTGCAAACGTTGTGAAGTCGCCTGCCCGTCCGATGTGAAGATTGGCGATATTATCCAGCGCGCGCGGGCGAAATATGACACCACGCGCCCGTCGCTGCGTAATTTTGTGTTGAGTCATACCGACCTGATGGGTAGCGTTTCCACGCCGTTCGCATCAATCGTCAACACCGCTACCTCGCTGAAACCGGTGCGGCAGCTGCTTg                                                                                                                                                                                                    <  1:222314/250‑1 (MQ=255)                                              tGAAAGATGGCGCACTGTATGACGAGGCGCTGAAATATTGCATCAACTGCAAACGTTGTGAAGTCGCCTGCCCGTCCGATGTGAAGATTGGCGATATTATCCAGCGCGCGCGGGCGAAATATGCCACCACGCGCCCGTCGCTGCGTAATTTTGTGTTGAGTCATACCGACCTGATGGGTAGCGTTTCCACGCCGTTCGCATCAATCGTCAACACCGCTACCTCGCTGAAACCGGTGCGGCAGCTGCTTGa                                                                                                                                                                                                   <  1:230218/250‑1 (MQ=255)                                               gAAAGATGGCGCTCGGTATTACGAGGCGCTGAAATATTTCAGCAACTGCAAACGTTGTGAAGTCGCCTGCCCGTCCGATGTGAAGATTGGCGATATTATCCAGCGCGCGCGGGCGAAATATGACACCACGCGCCCGTCGCTGCGTAATTTTGTGTTGAGTCATACCGACCTGATGGGTAGCGTTTCCACGCCGTTCGCATCAATCGTCAACACCGCTACCTCGCTGAAACCGGTGCGGCAGCTGCTTGATg                                                                                                                                                                                                 <  2:32590/251‑1 (MQ=255)                                                                aTGACGAGGCGCTGAAATATTGCATCAACTGCAAACGTTGTGAAGTCGCCTGCCCGTCCGATGTGAAGATTGGCGATATTATCCAGCGCGCGCGGGCGAAATATGACACCACGCGCCCGTCGCTGCGTAATTTTGTGTTGAGTCATACCGACCTGATGGGTAGCGTTTCCACGCCGTTCGCATCAATCGTCAACACCGCTACCTCGCTGAAACCGGTGCGGCAGCTGCTTGATGCGGCGTTAAAAATCGat                                                                                                                                                                                <  1:347120/251‑1 (MQ=255)                                                                  gACGAGGCGTTGAAATATTGCATCAACTGCAAACGTTGTGAAGTCGCCTGCCCGCCCGATGGGAAGAGTGGCGATATTATCCAGCGCGCGCGGGCGAAATATGCCACCACGCGCCCGTCGCTGCGTAATTTTGTGTTGAGTCATACCGACCTGCTGGGTAGCGTTTCCACGCCGTGCGCATCAATCGTCAACACCGCTACCTCGCTGAAACCGGTGCGGCAGCTGCTTGATGCGGCGTTAAAAATCGatc                                                                                                                                                                               <  2:81181/250‑1 (MQ=255)                                                                                     gCATCAACTGCAAACGTTGTGAAGTCGCCTGCCCGTCCGATGTGAAGAGTGGCGATCTTATCCAGCGCGCGCGGGCGAAATATGACACCACGCGCCCGTCGCTGCGTAATTTTGTGTTGAGTCATACCGACCTGATGGGTAGCGTTTCCACGCCGTTCGCATCAATCGTCAACACCGCTACCTCGCTGAAACCGGTGCGGCAGCTGCTTGATGCGGCGTTAAAAATCGATCATCGCCGCACGCTACCGaa                                                                                                                                                            <  1:271016/250‑1 (MQ=255)                                                                                          aacggCAAACGTTTTGAAGTCGCCTGCCCGTCCGATGTGAAGATTGGCGATATTATCCAGCGCGCGCGGGCGAAATATGACACCACGCGCCCGTCGCTGCGTAATTTTGTGTTGAGTCATACCGACCTGATGGGGAGCGGTGCCACGCCGTTCGCATCAATCGTCAACACCGCGTCCTCGCTGAAGCCGGTGCGGCAGCTGCTTGATGCGGCGTTAAAAATCGATCATCGCCGCACGCTACCGAAATACTc                                                                                                                                                      <  1:405628/247‑1 (MQ=255)                                                                                             gccAAAGGTTGT‑AAGTCGCCTGCCCGTCCGATGGGAAGATTGGCGATATTATCCAGCGCGCGCGGGCGAAATATGACACCACGCGCCCGTCGCTGCGTAATTTTGTGTTGAGTCATACCGACCTGATGGGTAGCGTTTCCACGCCGTTCGCATCAATCGTCAACTACGCTACCGCGCTGAAACCGGTGCGGCAGCTGCTTGATGCGGCGTTAAAAACCGATCATCGCCGCACGTTACCGAAATACtatgtc                                                                                                                                                  <  2:330994/249‑6 (MQ=255)                                                                                              gCAAACGTTGTGAAGTCGCCTGCCCGTCCGATGTGAAGATTGGCGATATTATCCAGCGCGCGCGGGCGAAATATGACACCACGCGCCCGTCGCTGCGTAATTTTGTGTTGAGTCATACCGACCTGATGGGTAGCGTTTCCACGCCGTTCGCATCAATCGTCAACACCGCTACCTCGCTGAAACCGGTGCGGCAGCTGCTTGATGCGGCGTTAAAAATCGATCATCGCCGCACGCTACCGAAATACTCCtt                                                                                                                                                   <  2:208518/250‑1 (MQ=255)                                                                                                   cGTTGTGAAGTCGCCTGCCCGTCCGATGTGAAGATTGGCGATATTATCCAGCGCGCGCGGGCGAAATATGACACCACGCGCCCGTCGCTGCGTAATTTTGTGTTGAGTCATACCGACCTGATGGGTAGCGTTTCCACGCCGTTCGCATCAATCGTCAACACCGCTACCTCGCTGAAACCGGTGCGGCAGCTGCTTGATGCGGCGTTAAAAAt                                                                                                                                                                                    <  1:363210/212‑1 (MQ=255)                                                                                                   cGTTGTGAAGTCGCCTGCCCGTCCGATGTGAAGATTGGCGATATTATCCAGCGCGCGCGGGCGAAATATGACACCACGCGCCCGTCGCTGCGTAATTTTGTGTTGAGTCATACCGACCTGATGGGTAGCGTTTCCACGCCGTTCGCATCAATCGTCAACACCGCTACCTCGCTGAAACCGGTGCGGCAGCTGCTTGATGCGGCGTTAAAAAt                                                                                                                                                                                    >  2:363210/1‑212 (MQ=255)                                                                                                       gtgAAGTCGCCGGCCCGTCCGATGTGAAGATTGGCGATATTATCCAGCGCGCGCGGGCGAAATATGACACCACGCGCCCGTCGCTGCGTAATTTTGTGTTGAGTCATACCGACCTGATGGGTAGCGTTTCCACGCCGGTCGCATCAATCGTCAACACCGCTTCCTCGCTGAAACCGGTGCGTCAGCGGCTTGATGCGGCGTTAAAAATCGATTATCGCCGCACGCTAGCGAAATTCTCCTTCGGCTCGTTc                                                                                                                                         <  2:241468/251‑1 (MQ=255)                                                                                                                ccTGCCCGTCCGATGTGAAGATTGGCGATATTATCCAGCGCGCGCGGGCGAAATATGACACCACGCGCCCGTCGCTGCGTAATTTTGTGTTGAGTCATACCGACCTGATGGGTAGCGTTTCCACGCCGTTCGCATCAATCGTCAACACCGCTACCTCGCTGAAACCGGTGCGGCAGCTGCTTGATGCGGCGTTAAAAATCGATCATCGCCGCACGCTACCGAAATACTCATTGGGCACGTTCCGTCGCTgg                                                                                                                                >  1:395403/1‑251 (MQ=255)                                                                                                                          cGATGTGAAGATTGGCGATATTATCCAGCGCGCGCGGGCGAAATATGACACCACGCGCCCGTCGCTGCGTAATTTTGTGTTGAGTCATACCGACCTGATGGGTAGCGTTTCCACGCCGTTCGCATCAATCGTCAACACCGCTACCTCGCTGAAACCGGTGCGGCAGCTGCTTGATGCGGCGTTAAAAATCGATCATCGCCGCACGCTACCGAAATACTCCTTCGGCACGTTCCGTCGCTGGTATCGCAGCg                                                                                                                      <  1:42370/251‑1 (MQ=255)                                                                                                                                gAAGATTGGCGATATTATCCAGCGCGCGCGGGCGAAATATGACACCACGCGCCCGTCGCTGCGTAATTTTGTGTTGAGTCATACCGACCTGATGGGTAGCGTTTCCACGCCGTTCGCATCAATCGTCAACACCGCTACCTCGCTGAAACCGGTGCGGCAGCTGCTTGATGCGGCGTTAAAAATCGATCATCGCCGCACGCTACCGAAATACTCCTTCGGCACGTTCCGTCGCTGGTATCGCAGCGTggcgg                                                                                                                >  1:352330/1‑251 (MQ=255)                                                                                                                                    aTTGGCGATATTATCCAGCGCGCGCGGGCGAAATATGACACCACGCGCCCGTCGCTGCGTAATTTTGTGTTGAGTCATACCGACCTGATGGGTAGCGTTTCCACGCCGTTCGCATCAATCGTCAACACCGCTACCTCGCTGAAACCGGTGCGGCAGCTGCTTGATGCGGCGTTAAAAATCGATCATCGCCGCACGCTACCGAAATACTCCTTCGGCACGTTCCGTCGCTGGTATCGCAGCGTGGCGGCTCa                                                                                                            <  1:3329/251‑1 (MQ=255)                                                                                                                                           atatTATCCAGCGCGCGCGGGCGAAATATGGCACCACTCGCCCGTCGCTGCGTAATTTTGTGTTGAGTCATACCGACCTGATGGGTAGCGTTTCCACGCCGTTCGCATCAATCGTCAACACCGATACCTCGCTGAAACCGGTGCGGCAGCTGCTTGATGCGGCGTTAAAAATCGATCAGCGCCGCTCGCTACCGAAATACGCCTTCGGCACGTTCCGTCGCTGGTATCGCAGCGTGGCGGCTCagcaagca                                                                                                     <  2:309165/251‑1 (MQ=255)                                                                                                                                           atatTATCCAGCGCGCGCGGGCGAAATATGACACCACGCGCCCGTCGCTGCGTAATTTTGTGTTGAGTCATACCGACCTGATGGGTAGCGTTTCCACGCCGTTCGCATCAATCGTCAACACCGCTACCTCGCTGAAACCGGTGCGGCAGCTGCTTGATGCGGCGTTAAAAATCGATCATCGCCGCACGCTACCGAAATACTCCTTCGGCACGTTCCGTCGCTGGTATCGCAGCGTGGCGGCTCagcaagc                                                                                                      >  1:374287/1‑250 (MQ=255)                                                                                                                                                 tCCAGCGCGCGCGGGCGAAATATGACACCACGCGCCCGTCGCTGCGTAATTTTGTGTTGAGTCATACCGACCTGATGGGTAGCGTTTCCACGCCGTTCGCATCAATCGTCAACACCGCTACCTCGCTGAAACCGGTGCGGCAGCTGCTTGATGCGGCGTTAAAAATCGATCATCGCCGCACGCTACCGAAATACTCCTTCGGCACGTTCCGTCGCTGGTATCGCAGCGTGGCGGCTCAGCAAGCACAatat                                                                                               >  1:193096/1‑251 (MQ=255)                                                                                                                                                   catcgcgcggggcggaaatactacacCACGCGCCTGTCTCTGGGTAAGTGTGTGTTGAGTCATACCGCCGTGATGGGTAGCGTTTCCACGCCGTTCGCATCAATCGTCAACACCGCTTCCTCGCTGAAACCGGTGCGGGAGCTGCTTGAGGCGGCGTTAAAAATCGATCAGCGCCGCACGCTGTCGAAATTCTCCTGAGGCACTTTCCGTCGCTGTTATCGCAGCGTTGTGGCGCAGAAAGCACAATATaa                                                                                             <  2:405203/229‑1 (MQ=255)                                                                                                                                                          cgcgGGCGAAATATGACACCACGCGCCCGTCGCTGCGTAATTTTGTGTTGAGTCATACCGACCTGATGGGTAGCGTTTCCACGCCGTTCGCATCAATCGTCAACACCGCTACCTCGCTGAAACCGGTGCGGCAGCTGCTTGATGCGGCGTTAAAAATCGATCATCGCCGCACGCTACCGAAATACTCCTTCGGCACGTTCCGTCGCTGGTATCGCAGCGTGGCGGCTCAGCAAGCACAATATAAAGACCAg                                                                                      <  1:230721/251‑1 (MQ=255)                                                                                                                                                                              ccgcgCCCGTCGCTGCGTAATTTTGTGTTGAGTCATACCGACCTGATGGGTAGCGTTTCCACGCCGTTCGCATCAATCGTCAACACCGCTACCTCGCTGAAACCGGTGCGGCAGCTGCTTGATGCGGCGTTAAAAATCGATCATCGCCGCGCGCTACCGAAATACTCCTTCGGCACGTTCCGTCGCTGGTATCGCAGCGTGGCGGCTCAGCAAGCACAATATAAAGACCAGGTCGCTTTCTTTCACGgctg                                                                  <  2:31677/250‑1 (MQ=255)                                                                                                                                                                                           tGCGTAATTTTGTGTTGAGTCATACCGACCTGATGGGTAGCGTTTCCACGCCGTTCGCATCAATCGTCAACACCGCTACCTCGCTGAAACCGGTGCGGCAGCTGCTTGATGCGGCGTTAAAAATCGATCATCGCCGCACGCTACCGAAATACTCCTTCGGCACGTTCCGTCGCTGGTATCGCAGCGTGGCGGCTCAGCAAGCACAATATAAAGACCAGGTCGCTTTATTTCACGGCTGCTTCGTTAACTAc                                                     >  1:108868/1‑251 (MQ=255)                                                                                                                                                                                                 aTTTTGTGTTGAGTCATACCGACCTGATGGGTAGCGTTTCCACGCCGTTCGCATCAATCGTCAACACCGCTACCTCGCTGAAACCGGTGCGGCAGCTGCTTGATGCGGCGTTAAAAAGCGATCATCGCCGCACGCTACCGAAATACTCCTTCGGCACGTTCCGTCGCTGGTATCGCAGCGTGGCGGCTCAGCAAGCACAATATAAAGACCAGGTCGCTTTCTTTCACGGCTGCTTCGTTAACTACAACCAt                                               <  1:394531/251‑1 (MQ=255)                                                                                                                                                                                                   tttGTGTTGAGTCTTACCGACCTGATGGGTAGCGCTTACATGCCGTTCGCATCAATCGTCAACACCGCTACCTCGCTGAAACCGGTGCGGCTGCTGCTTGATACGGCGTTAAAAATCGATCATCGCCGCACGCTACCGAACTACTCCTTCGGCACGTTCCGGCGCTGGGATAGCAGCGTGACGGCTCAGCAAGCACAGTATAAAGACCAGGTCGCGTTCTTTCACGCCTGCTTCAGTAACTACAACCAtac                                             >  2:81967/1‑249 (MQ=255)                                                                                                                                                                                                      gtgtTGAGTCATACCGACCTGATGGGTAGCGTTGCCACGCCGTTCGCATCAATCGTCAACACCGCTACCTCGCTGAAACCGGTGCGGCAGCTGCTTGATGCGGCGTTAAAAATCGATCATCGCCGCACGCTACCGAAATACGCCTTCGGCACGTTCCGTCGCTGGTATCGCAGCGTGGCGGCTCAGCAAGCACAATATAAAGACCAGGTCGCTTTCTTTCACGGCTGCTTCGTTAACTACAACCATCCGCa                                          <  1:211800/251‑1 (MQ=255)                                                                                                                                                                                                                    cGACCTGATGGGTAGCGTTTCCACGCCGTTCGCATCAATCGTCAACACCGCTACCTCGCTGAAACCGGTGCGGCAGCTGCTTGATGCGGCGTTAAAAATCGATCATCGCCGCACGCTACCGAAATACTCCTTCGGCACGTTCCGTCGCTGGTATCGCAGCGTGGCGGCTCAGCAAGCACAATATAAAGACCAGGTCGCTTTCTTTCACGGCTGCTTCGTTAACTACAACCATCCGCAGTTAGGTAAAGAtt                            >  1:12952/1‑251 (MQ=255)                                                                                                                                                                                                                     gACCTGATGGGTAGCGTTTCCACGCCGTTCGCATCAATCGTCAACACCGCTACCTCGCTGAAACCGGTGCGGCAGCTGCTTGATGCGGCGTTAAAAATCGATCATCGCCGCACGCTACCGAAATACTCCTTCGGCACGTTCCGTCGCTGGTATCGCAGCGTGGCGGCTCAGCAAGCACAATATAAAGACCAGGTCGCTTTCTTTCACGGCTGCTTCGTTAACTACAACCATCCGCAGTTAGGTAAAGAttt                           <  1:176572/251‑1 (MQ=255)                                                                                                                                                                                                                           aggggTAGCGGTTCCACGCCGTTCGCATCAATCGTCAACACCGCTACCTCGCTGAAACCGGTGCGGCAGCTGCTTGATGCGGCGTTAAAAATCGATCATCGCCGCACGCTACCGAAATACTCCTTCGGCACGTTCCGTCGCTGGTATCGCAGCGTGGCGGCTCAGCAAGCACAATATAAAGACCAGGTCGCTTTCTTTCACGGCTGCTTCGTTAACTACAACCATCCGCAGTTAGGTAAAGAGttaattaa                     <  2:301629/249‑1 (MQ=255)                                                                                                                                                                                                                           aTGGGTAGCGTTTCTACGCCGTTCGCATCAATCGTCAACACCGCTACATCGCTGAAACCGGTGCGGCAGCTGCTTGATGCGGCGTTAAAAATCGAGCTTCGCCGCACGCTACCGAAATACTCCTTCGGCACGTTCCGTCGCTGGTATCGCAGCGTGGCGGCTCAGCAAGCACAATATAAAGACCAGGTCGCTTTCGTTAACGGCTGCTTCGTTAACTACACCCATCCGCAGTTAGGTAAAGATttaattaa                     <  2:29116/251‑1 (MQ=255)                                                                                                                                                                                                                                               gTTCGCATCAATCGTCAACACCGCTACCTCGCTGAAACCGGTGCGGCAGCTGCTTGATGCGGCGTTAAAAATCGATCATCGCCGCACGCTACCGAAATACTCCTTCGGCACGTTCCGTCGCTGGTATCGCAGCGTGGCGGCTCAGCAAGCACAATATAAAGACCAGGTCGct                                                                                >  2:299146/1‑172 (MQ=255)                                                                                                                                                                                                                                               gTTCGCATCAATCGTCAACACCGCTACCTCGCTGAAACCGGTGCGGCAGCTGCTTGATGCGGCGTTAAAAATCGATCATCGCCGCACGCTACCGAAATACTCCTTCGGCACGTTCCGTCGCTGGTATCGCAGCGTGGCGGCTCAGCAAGCACAATATAAAGACCAGGTCGCTTTCTTTCACGGCTGCTTCGTTAACTACAACCATCCGCAGTTAGGTAAAGATTTAATTAAAGTGCTCAACGCAATGGGt  >  1:374797/1‑250 (MQ=255)                                                                                                                                                                                                                                               gTTCGCATCAATCGTCAACACCGCTACCTCGCTGAAACCGGTGCGGCAGCTGCTTGATGCGGCGTTAAAAATCGATCATCGCCGCACGCTACCGAAATACTCCTTCGGCACGTTCCGTCGCTGGTATCGCAGCGTGGCGGCTCAGCAAGCACAATAGAAAGACCAGGTCGct                                                                                <  1:299146/172‑1 (MQ=255)                                                                                                                                                                                                                                                      |                                                                                                                                                                                                                                                   ATCCAGGGCCAAAACAAGCCGGGCCGGATGGCGAGCGTCTGCGTTTGAAAGATGGCGCACTGTATGACGAGGCGCTGAAATATTGCATCAACTGCAAACGTTGTGAAGTCGCCTGCCCGTCCGATGTGAAGATTGGCGATATTATCCAGCGCGCGCGGGCGAAATATGACACCACGCGCCCGTCGCTGCGTAATTTTGTGTTGAGTCATACCGACCTGATGGGTAGCGTTTCCACGCCGTTCGCACCAATCGTCAACACCGCTACCTCGCTGAAACCGGTGCGGCAGCTGCTTGATGCGGCGTTAAAAATCGATCATCGCCGCACGCTACCGAAATACTCCTTCGGCACGTTCCGTCGCTGGTATCGCAGCGTGGCGGCTCAGCAAGCACAATATAAAGACCAGGTCGCTTTCTTTCACGGCTGCTTCGTTAACTACAACCATCCGCAGTTAGGTAAAGATTTAATTAAAGTGCTCAACGCAATGGGT  >  NC_000913/2355603‑2356090

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