Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	NC_000913	3,012,553	Δ1 bp	intergenic (+160/+61)	yqeA → / ← yqeB	putative amino acid kinase/XdhC‑CoxI family protein with NAD(P)‑binding Rossman fold

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	3,012,552	0	G	.	100.0%	275.1 / NA	68	intergenic (+159/+62)	yqeA/yqeB	putative amino acid kinase/XdhC‑CoxI family protein with NAD(P)‑binding Rossman fold
Reads supporting (aligned to +/- strand):  ref base G (0/0);  new base . (35/33);  total (35/33)

TCCGGCCAAAGATTGTCAAGTAAGTGCGCTTTGTTGATGCCGGATGCGGCGCGAGCGCCTTATCCGTCCTACAAAACATCGCAAATTCAATAGATTGCAATGCACTTGTAGGCCTGATAAGCGTAGCGCATCAGGCAATTTGGCGTTTGTCATCAGTCTCAGGCCGGGTGGGTAATACCATCCGGCCATTTTCATTTAAGCCACTTCCAGCACTTCTTTTGTCGCTTTCACGCCCTGATGCATCAACATCATTAACGCC  >  NC_000913/3012419‑3012677                                                                                                                                      |                                                                                                                              tCCGGCCAAAGATTGTCAAGTAAGTGCGCTTTGTTGATGCCGGATGCGGCGCGAGCGCCTTATCCGTCCTACAAAACATCGCAAATTCAATAGATTGCAATGCACTTGTAGGCCTGATAAGCGTAGCGCATCAg                                                                                                                               <  1:240422/134‑1 (MQ=255) tCCGGCCAAAGATTGTCAAGTAAGTGCGCTTTGTTGATGCCGGATGCGGCGCGAGCGCCTTATCCGTCCTACAAAACATCGCAAATTCAATAGATTGCAATGCACTTGTAGGCCTGATAAGCGTAGCGCATCAg                                                                                                                               >  2:240422/1‑134 (MQ=255) tCCGGCCAAAGATTGTCAAGTAAGTGCGCTTTGTTGATGCCGGATGCGGCGCGAGCGCCTTATCCGTCCTACAAAACATCGCAAATTCAATAGATTGCAATGCACTTGTAGGCCTGATAAGCGTAGCGCATCA‑GCAAtt                                                                                                                         >  2:566443/1‑139 (MQ=255) tCCGGCCAAAGATTGTCAAGTAAGTGCGCTTTGTTGATGCCGGATGCGGCGCGAGCGCCTTATCCGTCCTACAAAACATCGCAAATTCAATAGATTGCAATGCACTTGTAGGCCTGATAAGCGTAGCGCATCA‑GCAAtt                                                                                                                         >  1:110139/1‑139 (MQ=255) tCCGGCCAAAGATTGTCAAGTAAGTGCGCTTTGTTGATGCCGGATGCGGCGCGAGCGCCTTATCCGTCCTACAAAACATCGCAAATTCAATAGATTGCAATGCACTTGTAGGCCTGATAAGCGTAGCGCATCA‑GCAAtt                                                                                                                         >  1:687604/1‑139 (MQ=255) tCCGGCCAAAGATTGTCAAGTAAGTGCGCTTTGTTGATGCCGGATGCGGCGCGAGCGCCTTATCCGTCCTACAAAACATCGCAAATTCAATAGATTGCAATGCACTTGTAGGCCTGATAAGCGTAGCGCATCA‑GCAAtt                                                                                                                         >  1:9864/1‑139 (MQ=255)   ccgCCAAAGATTGTCAAGTAAGTGCGCTTTGTTGATGCTGGATGCGGCGCGAGCGCCTTATCCGTCCTACAAAACATCGCAAATTCAATAGATTGGAATGCACTTGTAGGACTGCTAAGCGTAGCGCATCA‑GCAATTTg                                                                                                                       >  1:724532/3‑139 (MQ=25)     gCCAAAGATTGTCAAGTAAGTGCGCTTTGTTGATGCCGGATGCGGCGCGAGCGCCTTATCCGTCCTACAAAACATCGCAAATTCAATAGATTGCAATGCACTTGTAGGCCTGATAAGCGTAGCGCATCA‑GCAATTTGGc                                                                                                                     >  1:161106/1‑139 (MQ=255)      ccAAAGATTGTCAAGTAAGTGCGCTTTGTTGATGCCGGATGCGGCGCGAGCGCCTTATCCGTCCTACAAAACATCGCAAATTCAATAGATTGCAATGCACTTGTAGGCCTGATAAGCGTAGCGCATCA‑GCAATTTGGCg                                                                                                                    >  1:610918/1‑139 (MQ=255)      ccAAAGATTGTCAAGTAAGTGCGCTTTGTTGATGCCGGATGCGGCGCGAGCGCCTTATCCGTCCTACAAAACATCGCAAATTCAATAGATTGCAATGCACTTGTAGGCCTGATAAGCGTAGCGCATCA‑GCAATTTGGCg                                                                                                                    >  2:532873/1‑139 (MQ=255)      ccAAAGATTGTCAAGTAAGTGCGCTTTGTTGATGCCCGATGCGGCGCTAGCGCCTTATCCGTCGTACAAAACATCGCAAATTCAATATATTGCAATGCACTTGTAGGCCTGATAAGCGTAGCGCATCA‑GCAATTTGGCg                                                                                                                    >  1:672618/1‑139 (MQ=21)              tGTCAAGTAAGTGCGCTTTGTTGATGCCGGATGCGGCGCGAGCGCCTTATCCGTCCTACAAAACATCGCAAATTCAATAGATTGCAATGCACTTGTAGGCCTGATAAGCGTAGCGCATCA‑GCAATTTGGCGTTTGtcat                                                                                                            <  1:532873/139‑1 (MQ=255)              tGTCAAGTAAGTGCGCTTTGTTGATGCCGGATGCGGCGCGAGCGCCTTATCCGTCCTACAAAACATCGCAAATTCAATAGATTGCAATGCACTTGTAGGCCTGATAAGCGTAGCGCATCA‑GCAATTTGGCGTTTGtcat                                                                                                            <  2:1053550/139‑1 (MQ=255)                tCAAGTAAGTGCGCTTTGTTGATGCCGGATGCGGCGCGAGCGGCTTATCCGTCCTACAAAACATCGCAAATTCAATAGATTGCAATGCACGTGTAGGCCTGATAAGCGTAGCGCATCA‑GCAATTTGGCGTTTGtcatca                                                                                                          >  2:377364/1‑139 (MQ=255)                    gtaagtGCGCTTTGTTGATGCCGGATGCGGCGCGAGCGCCTTATCCGTCCTACAAAACATCGCAAATTCAATAGATTGCAATGCACTTGTAGGCCTGATAAGCGTAGCGCATCA‑GCAATTTGGCGTTTGTCATCAGtct                                                                                                      <  2:161106/139‑1 (MQ=255)                     taagtGCGCTTTGTTGATGCCGGATGCGGCGCGAGCGCCTTATCCGTCCTACAAAACATCGCAAATTCAATAGATTGCAATGCACTTGTAGGCCTGATAAGCGTAGCGCATCA‑GCAATTTGGCGTTTGTCATCAGtctc                                                                                                     >  2:354329/1‑139 (MQ=255)                            gcTTTGTTGATGCCGGATGCGGCGCGAGCGCCTTATCCGTCCTACAAAACATCGCAAATTCAATAGATTGCAATGCACTTGTAGGCCTGATAAGCGTAGCGCATCA‑GCAATTTGGCGTTTGTCATCAGTCTCAGGCCgg                                                                                              <  2:292691/139‑1 (MQ=255)                            gcTTTGTTGATGCCGGATGCGGCGCGAGCGCCTTATCCGTCCTACAAAACATCGCAAATTCAATAGATTGCAATGCACTTGTAGGCCTGATAAGCGTAGCGCATCA‑GCAATTTGGCGTTTGTCATCAGTCTCAGGCCgg                                                                                              <  2:110139/139‑1 (MQ=255)                            gcTTTGTTGATGCCGGATGCGGCGCGAGCGCCTTATCCGTCCTACAAAACATCGCAAATTCAATAGATTGCAATGCACTTGTAGGCCTGATAAGCGTAGCGCATCA‑GCAATTTGGCGTTTGTCATCAGTCTCAGGCCgg                                                                                              <  1:3125/139‑1 (MQ=255)                            gcTTTGTTGATGCCGGATGCGGCGCGAGCGCCTTATCCGTCCTACAAAACATCGCAAATTCAATAGATTGCAATGCACTTGTAGGCCTGATAAGCGTAGCGCATCA‑GCAATTTGGCGTTTGTCATCAGTCTCAGGCCgg                                                                                              <  2:334259/139‑1 (MQ=255)                             cTTTGTTGATGCCGGATGCGGCGCGAGCGCCTTATCCGTCCTACAAAACATCGCAAATTCAATAGATTGCAATGCACTTGTAGGCCTGATAAGCGTAGCGCATCA‑GCAATTTGGCGTTTGTCATCAGTCTCAGGCCggg                                                                                             >  1:714204/1‑139 (MQ=255)                                tgttgATGCCGGATGCGGCGCGAGCGCCTTATCCGTCCTACAAAACATCGCAAATTCAATAGATTGCAATGCACTTGTAGGCCTGATAAGCGTAGCGCATCA‑GCAATTTGGCGTTTGTCATCAGTCTCAGGCCgggtgg                                                                                          <  2:341054/139‑1 (MQ=255)                                      tGCCGGATGCGGCGCGAGCGCCTTATCCGTCCTACAAAACATCGCAAATTCAATAGATTGCAATGCACTTGTAGGCCTGATAAGCGTAGCGCATCA‑GCAATTTGGCGTTTGTCATCAGTCTCAGGCCGGGTGGGTAATa                                                                                    <  1:354329/139‑1 (MQ=255)                                          ggATGCGGCGCGAGCGCCTTATCCGTCCTACAAAACATCGCAAATTCAATAGATTGCAATGCACTTGTAGGCCTGATAAGCGTAGCGCATCA‑GCAATTTGGCGTTTGTCATCAGTCTCAGGCCGGGTGGGTAATACCAt                                                                                >  1:981936/1‑139 (MQ=255)                                          ggATGCGGCGCGAGCGCCTTATCCGTCCTACAAAACATCGCAAATTCAATAGATTGCAATGCACTTGTAGGCCTGATAAGCGTAGCGCATCA‑GCAATTTGGCGTTTGTCATCAGTCTCAGGCCGGGTGGGTAATACCAt                                                                                >  2:1006065/1‑139 (MQ=255)                                            aTGCGGCGCGAGCGCCTTATCCGTCCTACAAAACATCGCAAATTCAATAGATTGCAATGCACTTGTAGGCCTGATAAGCGTAGCGCATCA‑GCAATTTGGCGTTTGTCATCAGTCTCAGGCCGGGTGGGTAATACCATcc                                                                              >  1:102093/1‑139 (MQ=255)                                                   gcgAGCGCCTTATCCGTCCTACAAAACATCGCAAATTCAATAGATTGCAATGCACTTGTAGGCCTGATAAGCGTAGCGCATCA‑GCAATTTGGCGTTTGTCATCAGTCTCAGGCCGGGTGGGTAATACCATCCGGCCAtt                                                                       >  2:476439/1‑139 (MQ=255)                                                    cgAGCGCCTTATCCGTCCTACAAAACATCGCAAATTGAATAGATTGCAATGCACTTGTAGGCCTGATAAGCGTAGCGCATCA‑GCAATTTGGCGTTTGTCATCAGTCTCAGGCCGGGTGGGTAATACCATCCGGCCAttt                                                                      <  1:667203/139‑1 (MQ=255)                                                             tATCCGTCCTACAAAACATCGCAAATTCAATAGATTGCAATGCACTTGTAGGCCTGATAAGCGTAGCGCATCA‑GCAATTTGGCGTTTGTCATCAGTCTCAGGCCGGGTGGGTAATACCATCCGGCCATTTTCATTTAAg                                                             <  1:1006065/139‑1 (MQ=255)                                                             tATCCGTCCTACAAAACATCGCAAATTCAATAGATTGCAATGCACTTGTAGGCCTGATAAGCGTAGCGCATCA‑GCAATTTGGCGTTTGTCATCAGTCTCAGGCCGGGTGGGTAATACCATCCGGCCATTTTCATTTAAg                                                             <  1:950451/139‑1 (MQ=255)                                                               tCCGTCCTACAAAACATCGCAAATTCAATAGATTGCAATGCACTTGTAGGCCTGATAAGCGTAGCGCATCA‑GCAATTTGGCGTTTGTCATCAGTCTCAGGCCGGGTGGGTAATACCATCCGGCCATTTTCATTTAAGcc                                                           <  1:886008/139‑1 (MQ=255)                                                                   tCCTACAAAACATCGCAAATTCAATAGATTGCAATGCACTTGTAGGCCTGATAAGCGTAGCGCATCA‑GCAATTTGGCGTTTGTCATCAGTCTCAGGCCGGGTGGGTAATACCATCCGGCCATTTTCATTTAAGCCACtt                                                       >  1:343555/1‑139 (MQ=255)                                                                        cAAAACATCGCAAATTCAATAGATTGCAATGCACTTGTAGGCCTGATAAGCGTAGCGCATCA‑GCAATTTGGCGTTTGTCATCAGTCTCAGGCCGGGTGGGTAATACCATCCGGCCATTTTCATTTAAGCCACTTCCAGc                                                  >  1:547537/1‑139 (MQ=255)                                                                         aaaaCATCGCACATTCAATAGATTGCAATGCACTTGTAGGCCTGATAAGCGTAGCGCATCA‑GCAATTTGGCGTTTGTCATCAGTCTCAGGCCGGGTGGGTAATACCATCCGGCCATTTTCATTTAAGCCACTTCCAGCa                                                 <  1:683972/139‑1 (MQ=255)                                                                             cATCGCAAATTCAATAGATTGCAATGCACTTGTAGGCCTGATAAGCGTAGCGCATCA‑GCAATTTGGCGTTTGTCATCAGTCTCAGGCCgg                                                                                              >  1:199280/1‑90 (MQ=255)                                                                             cATCGCAAATTCAATAGATTGCAATGCACTTGTAGGCCTGATAAGCGTAGCGCATCA‑GCAATTTGGCGTTTGTCATCAGTCTCAGGCCgg                                                                                              <  2:199280/90‑1 (MQ=255)                                                                             cATCGCAAATTCAATAGATTGCAATGCACTTGTAGGCCTGATAAGCGTAGCGCATCA‑GCAATTTGGCGTTTGTCATCAGTCTCAGGCCGGGTGTGTAATACCATCCGGCCATTTTCATTTAAGCCACTTCCAGTCcttc                                             >  2:288233/1‑139 (MQ=37)                                                                                cGCAAATTCAATAGATTGCAATGCACTTGTAGGCCTGATAAGCGTAGCGCATCA‑GCAATTTGGCGTTTGTCATCAGTCTCAGGCCgg                                                                                              <  1:704159/87‑1 (MQ=255)                                                                                cGCAAATTCAATAGATTGCAATGCACTTGTAGGCCTGATAAGCGTAGCGCATCA‑GCAATTTGGCGTTTGTCATCAGTCTCAGGCCgg                                                                                              >  2:704159/1‑87 (MQ=255)                                                                                cGCAAATTCAATAGATTGCAATGCACTTGTAGGCCTGATAAGCGTAGCGCATCA‑GCAATTTGGCGTTTGTCATCAGTCTCAGGCCGGGTGGGTAATACCATCCGGCCATTTTCATTTAAGCCACTTCCAGCa                                                 <  2:955755/132‑1 (MQ=255)                                                                                cGCAAATTCAATAGATTGCAATGCACTTGTAGGCCTGATAAGCGTAGCGCATCA‑GCAATTTGGCGTTTGTCATCAGTCTCAGGCCGGGTGGGTAATACCATCCGGCCATTTTCATTTAAGCCACTTCCAGCa                                                 >  1:955755/1‑132 (MQ=255)                                                                                  cAAATTCAATAGATTGCAATGCACTTGTAGGCCTGATAAGCGTAGCGCATCA‑GCAATTTGGCGTTTGTCATCAGTCTCAGGCCGGGTGGGTAATACCATCCGGCCATTTTCATTTAAGCCACTTCCAGCACTTCTTTTg                                        >  2:418391/1‑139 (MQ=255)                                                                                  cAAATTCAATAGATTGCAATGCACTTGTAGGCCTGATAAGCGTAGCGCATCA‑GCAATTTGGCGTTTGTCATCAGTCTCAGGCCGGGTGGGTAATACCATCCGGCCATTTTCATTTAAGCCACTTCCAGCACTTCTTTTg                                        >  2:228598/1‑139 (MQ=255)                                                                                       tCAATAGATTGCAATGCACTTGTAGGCCTGATAAGCGTAGCGCATCA‑GCAATTTGGCGTTTGTCATCAGTCTCAGGCCGGGTGGGTAATACCATCCGGCCATTTTCATTTAAGCCACTTCCAGCACTTCTTTTGTCGCt                                   <  1:237226/139‑1 (MQ=255)                                                                                           tAGATTGCAATGCACTTGTAGGCCTGATAAGCGTAGCGCATCA‑GCAATTTGGCGTTTGTCATCAGTCTCAGGCCGGGTGGGTAATACCATCCGGCCATTTTCATTTAAGCCACTTCCAGCACTTCTTTTGTCGCTTTCa                               <  2:764596/139‑1 (MQ=255)                                                                                                tGCAATGCACTTGTAGGCCTGATAAGCGTAGCGCATCA‑GCAATTTGGCGTTTGTCATCAGTCTCAGGCCgg                                                                                              >  2:8734/1‑71 (MQ=255)                                                                                                tGCAATGCACTTGTAGGCCTGATAAGCGTAGCGCATCA‑GCAATTTGGCGTTTGTCATCAGTCTCAGGCCgg                                                                                              <  1:8734/71‑1 (MQ=255)                                                                                                   aaTGCACTTGTAGGCCTGATAAGCGTAGCGCATCA‑GCAATTTGGCGTTTGTCATCAGTCTCAGGCCgg                                                                                              >  1:197202/1‑68 (MQ=255)                                                                                                   aaTGCACTTGTAGGCCTGATAAGCGTAGCGCATCA‑GCAATTTGGCGTTTGTCATCAGTCTCAGGCCgg                                                                                              <  2:197202/68‑1 (MQ=255)                                                                                                   aaTGCACTTGTAGGCCTGATAAGCGTAGCGCATCA‑GCAATTTGGCGTTTGTCATCAGTCTCAGGCCGGGTGGGTAATACCATCCGGCCATTTTCATTTAAGCCACTTCCAGCACTTCTTTTGTCGCTTTCACGCCCTGa                       <  2:343555/139‑1 (MQ=255)                                                                                                     tGCACTTGTAGGCCTGATAAGCGTAGCGCATCA‑GCAATTTGGCGTTTGTCATCAGTCTCAGGCCGGGTGGGTAATACCATCCGGCCATTTTCATTTAAGCCACTTCCAGCACTTCTTTTGTCGCTTTCACGCCCTGATg                     >  2:654116/1‑139 (MQ=255)                                                                                                        aCTTGTAGGCCTGATAAGCGTAGCGCATCA‑GCAATTTGGCGTTTGtcat                                                                                                            >  2:227027/1‑49 (MQ=17)                                                                                                        aCTTGTAGGCCTGATAAGCGTAGCGCATCA‑GCAATTTGGCGTTTGtcat                                                                                                            <  1:227027/49‑1 (MQ=17)                                                                                                         cTTGTAGGCCTGATAAGCGTAGCGCATCA‑GCAATTTGGCGTTTGTCATCAGTCTCAGGCCGGGTGGGTAATACCATCCGGCCATTTTCATTTAAGCCACTTCCAGCACTTCTTTTGTCGCTTTCACGCCCTGATGCATc                 <  1:228598/139‑1 (MQ=255)                                                                                                         cTTGTAGGCCTGATAAGCGTAGCGCATCA‑GCAATTTGGCGTTTGTCATCAGTCTCAGGCCGGGTGGGTAATACCATCCGGCCATTTTCATTTAAGCCACTTCCAGCACTTCTTTTGTCGCTTTCACGCCCTGATGCATc                 <  1:377364/139‑1 (MQ=255)                                                                                                         cTTGTAGGCCTGATAAGCGTAGCGCATCA‑GCAATTTGGCGTTTGTCATCAGTCTCAGGCCGGGTGGGTAATACCATCCGGCCATTTTCATTTAAGCCACTTCCAGCACTTCTTTTGTCGCTTTCACGCCCTGATGCATc                 >  1:364818/1‑139 (MQ=255)                                                                                                         cTTGTAGGCCTGATAAGCGTAGCGCATCA‑GCAATTTGGCGTTTGTCATCAGTCTCAGGCCGGGTGGGTAATACCATCCGGCCATTTTCATTTAAGCCACTTCCAGCACTTCTTTTGTCGCTTTCACGCCCTGATGCATc                 >  2:520069/1‑139 (MQ=255)                                                                                                            gTAGGCCTGATAAGCGTAGCGCATCA‑GCAATTTGGCGTTTGTCATCAGTCTCAGGCCGGGTGGGTAATACCATCCGGCCATTTTCATTTAAGCCACTTCCAGCACTTCTTTTGTCGCTTTCACGCCCTGATGCATCAAc              <  1:566443/139‑1 (MQ=255)                                                                                                             tAGGCCTGATAAGCGTAGCGCATCA‑GCAATTTGGCGTTTGTCATCAGTCTCAGGCCGGGTGGGTAATACCATCCGGCCATTTTCATTTAAGCCACTTCCAGCACTTCTTTTGTCGCTTTCACGCCCTGATGCATCAAca             >  1:111806/1‑139 (MQ=255)                                                                                                                  cTGATAAGCGTAGCGCATCA‑GCAATTTGGCGTTTGTCATCAGTCTCAGGCCGGGTGGGTAATACCATCCGGCCATTTTCATTTAAGCCACTTCCAGCACTTCTTTTGTCGCTTTCACGCCCTGATGCATCAACATCAtt        <  1:476439/139‑1 (MQ=255)                                                                                                                    gATAAGCGTAGCGCATCA‑GCAATTTGGCGTTTGTCATCAGTCTCAGGCCGGGTGGGTAATACCATCCGGCCATTTTCATTTAAGCCACTTCCAGCACTTCTTTTGTCGCTTTCACGCCCTGATGCATCAACATCATTaa      <  1:654116/139‑1 (MQ=255)                                                                                                                    gATAAGCGTAGCGCATCA‑GCAATTTGGCGTTTGTCATCAGTCTCAGGCCGGGTGGGTAATACCATCCGGCCATTTTCATTTAAGCCACTTCCAGCACTTCTTTTGTCGCTTTCACGCCCTGATGCATCAACATCATTaa      <  1:520069/139‑1 (MQ=255)                                                                                                                     aTAAGCGTAGCGCATCA‑GCAATTTGGCGTTTGTCATCAGTCTCAGGCCGGGTGGGTAATACCATCCGGCCATTTTCATTTAAGCCACTTCCAGCACTTCTTTTGTCGCTTTCACGCCCTGATGCATCAACATCATTAAc     >  2:657714/1‑139 (MQ=255)                                                                                                                        aGCGTAGCGCATCA‑GCAATTTGGCGTTTGTCATCAGTCTCAGGCCGGGTGGGTAATACCATCCGGCCATTTTCATTTAAGCCACTTCCAGCACTTCTTTTGTCGCTTTCACGCCCTGATGCATCAACATCATTAACGcc  <  2:610918/139‑1 (MQ=255)                                                                                                                            tAGCGCATCA‑GCAATTTGGCGTTTGTCATCAGtct                                                                                                      <  1:712834/35‑1 (MQ=17)                                                                                                                            tAGCGCATCA‑GCAATTTGGCGTTTGTCATCAGtct                                                                                                      >  2:712834/1‑35 (MQ=17)                                                                                                                               cgcATCA‑GCAATTTGGCGTTTGTCATCAGTCTCAGGCCgg                                                                                              <  1:1051917/40‑1 (MQ=18)                                                                                                                               cgcATCA‑GCAATTTGGCGTTTGTCATCAGTCTCAGGCCGGGTGGGTAATACCATCCGGCCATTTTCATTTAAGCCACTTCCAGCACTTCTTTTGTCGCTTTCACGCCCTGATGc                    >  1:128685/1‑114 (MQ=255)                                                                                                                               cgcATCA‑GCAATTTGGCGTTTGTCATCAGTCTCAGGCCGGGTGGGTAATACCATCCGGCCATTTTCATTTAAGCCACTTCCAGCACTTCTTTTGTCGCTTTCACGCCCTGATGc                    <  2:128685/114‑1 (MQ=255)                                                                                                                               cgcATCA‑GCAATTGGGCGTTTGTCATCAGTCTCAGGCCgg                                                                                              >  2:1051917/1‑40 (MQ=18)                                                                                                                                      |                                                                                                                              TCCGGCCAAAGATTGTCAAGTAAGTGCGCTTTGTTGATGCCGGATGCGGCGCGAGCGCCTTATCCGTCCTACAAAACATCGCAAATTCAATAGATTGCAATGCACTTGTAGGCCTGATAAGCGTAGCGCATCAGGCAATTTGGCGTTTGTCATCAGTCTCAGGCCGGGTGGGTAATACCATCCGGCCATTTTCATTTAAGCCACTTCCAGCACTTCTTTTGTCGCTTTCACGCCCTGATGCATCAACATCATTAACGCC  >  NC_000913/3012419‑3012677

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