Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	NC_000913	3,248,939	T→C	intergenic (+117/‑30)	mzrA → / → yqjC	modulator of EnvZ/OmpR regulon/DUF1090 family putative periplasmic protein

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	3,248,939	0	T	C	100.0%	272.4 / NA	75	intergenic (+117/‑30)	mzrA/yqjC	modulator of EnvZ/OmpR regulon/DUF1090 family putative periplasmic protein
Reads supporting (aligned to +/- strand):  ref base T (0/0);  new base C (34/41);  total (34/41)

TCTCATCGCTTCGGATAACTTCCTGTTATCCGAAACATATCACTCACTTTGGTGATTTCACCGTAACTGTCTATGATTAATGAGGCGGTGGTTAACACTACCGCTTCGTTGTCTGGACAATCGTTCCTTTGTAATAGGTCCGAATCACAATGGAAGGTTCAAGAATGAAATACCGCATCGCTTTAGCTGTTTCTCTCTTTGCTCTTAGTGCCGGTAGTTATGCCACTACCCTGTGTCAGGAAAAGGAGCAAAATATCCTTAA  >  NC_000913/3248805‑3249066                                                                                                                                       |                                                                                                                                tctcATCGCTTCGGATAACTTCCTGTTATCCGAAACATATCACTCACTTTGGTGATTTCACCGTAACTGTCTATGATTAATGAGGCGGTGGTTAACACTACCGCTTCGTTGTCTGGACAATCGTTCCTTTGTAACAGGt                                                                                                                             <  2:181667/139‑1 (MQ=255)         cTTCGGATAACTTCCTGTTATCCGAAACATATCACTCACTTTGGTGATTTCACCGTAACTGTCTATGATTAATGAGGCGGTGGTTAACACTACCGCTTCGTTGTCTGGACAATCGTTCCTTTGTAACAGGTCCGAATca                                                                                                                     <  2:258395/139‑1 (MQ=255)         cTTCGGATAACTTCCTGTTATCCGAAACATATCACTCACTTTGGTGATTTCACCGTAACTGTCTATGATTAATGAGGCGGTGGTTAACACTACCGCTTCGTTGTCTGGACAATCGTTCCTTTGTAACAGGTCCGAATca                                                                                                                     <  1:399981/139‑1 (MQ=255)            cGGATAACTTCCTGTTATCCGAAACATATCACTCACTTTGGTGATTTCACCGTAACTGTCTATGATTAATGAGGCGGTGGTTAACACTACCGCTTCGTTGTCTGGACAATCGTTCCTTTGTAACAGGTCCGAATCACaa                                                                                                                  >  1:246885/1‑139 (MQ=255)              gATAACTTCCTGTTATCCGAAACATATCACTCACTTTGGTGATTTCACCGTAACTGTCTATGATTAATGAGGCGGTGGTTAACACTACCGCTTCGTTGTCTGGACAATCGTTCCTTTGTAACAGGTCCGAATCACAATg                                                                                                                <  2:246885/139‑1 (MQ=255)              gATAACTTCCTGTTATCCGAAACATATCACTCACTTTGGTGATTTCACCGTAACTGTCTATGATTAATGAGGCGGTGGTTAACACTACCGCTTCGTTGTCTGGACAATCGTTCCTTTGTAACAGGTCCGAATCACAATg                                                                                                                <  1:302779/139‑1 (MQ=255)                tAACTTCCTGTTATCCGAAACATATCACTCACTTTGGTGATTTCACCGTAACTGTCTATGATTAATGAGGCGGTGGTTAACACTACCGCTTCGTTGTCTGGACAATCGTTCCTTTGTAACAGGTCCGAATCACAATGGa                                                                                                              <  1:223295/139‑1 (MQ=255)                tAACTTCCTGTTATCCGAAACATATCACTCACTTTGGTGATTTCACCGTAACTGTCTATGATTAATGAGGCGGTGGTTAACACTACCGCTTCGTTGTCTGGACAATCGTTCCTTTGTAACAGGTCCGAATCACAATGGa                                                                                                              <  2:350592/139‑1 (MQ=255)                 agcTTCCTGTTATCCGAAACATATCACTCACTTTGGTGATTTCACCGTAACTGTCTATGATTAATGAGGCGGTGGTTAACTCTACCGCTTCGTTGTCTGGACAATCGTACCTTTGTAACAGGTCCGAATCACAATGGaa                                                                                                             <  1:81026/137‑1 (MQ=255)                     tCCTGTTATCCGAAACATATCACTCACTTTGGTGATTTCACCGTAACTGTCTATGATTAATGAGGCGGTGGTTAACACTACCGCTTCGTTGTCTGGACAATCGTTCCTTTGTAACAGGTCCGAATCACAATGGAAGCtt                                                                                                         <  1:465139/139‑1 (MQ=255)                      ccTGTTATCCGAAACATATCACTCACTTTGGTGATTTCACCGTAACTGTCTATGATTAATGAGGCGGTGGTTAACACTACCGCTTCGTTGTCTGGACAATCGTTCCTTTGTAACAGGTCCGAATCACAATGGAAGGTTc                                                                                                        <  2:689354/139‑1 (MQ=255)                             tCCGAAACATATCACTCACTTTGGTGATTTCACCGTAACTGTCTATGATTAATGAGGCGGTGGTTAACACTACCGCTTCGTTGTCTGGACAATCGTTCCTTTGTAACAGGTCCGAATCACAATGGAAGGTTCAAGAATg                                                                                                 >  1:456973/1‑139 (MQ=255)                                       atCACTCACTTTGGTGATTTCACCGTAACTGTCTATGATTAATGAGGCGGTGGTTAACACTACCGCTTCGTTGTCTGGACAATCGTTCCTTTGTAACAGGTCCGAATCACAATGGAAGGTTCAAGAATGAAATACCGCa                                                                                       >  1:204391/1‑139 (MQ=255)                                       atCACTCACTTTGGTGATTTCACCGTAACTGTCTATGATTAATGAGGCGGTGGTTAACACTACCGCTTCGTTGTCTGGACAATCGTTCCTTTGTAACAGGTCCGAATCACAATGGAAGGTTCAAGAATGAAATACCGCa                                                                                       <  1:513524/139‑1 (MQ=255)                                              actTTGGTGATTTCACCGTAACTGTCTATGATTAATGAGGCGGTGGTTAACACTACCGCTTCGTTGTCTGGACAATCGTTCCTTTGTAACAGGTCCGAATCACAATGGAAGGTTCAAGAATGAAATACCGCATCGCttt                                                                                >  1:722384/1‑139 (MQ=255)                                                tttGGTGATTTCACCGTAACTGTCTATGATTAATGAGGCGGTGGTTAACACTACCGCTTCGTTGTCTGGACAATCGTTCCTTTGTAACAGGTCCGAATCACAATGGAAGGTTCAAGAATGAAATACCGCATCGCTTTAg                                                                              >  2:496283/1‑139 (MQ=255)                                                     tGATTTCACCGTAACTGTCTATGATTAATGAGTCGGTGGTTAACACTACCGCTTCGTTGTCTGGACAATCGTTCCTTTGTAACAGGTCCGAATCACAATGGAAGGTTCAAGTATGAAATACCGCATCGCTTTAGCTGtt                                                                         <  2:682919/139‑1 (MQ=255)                                                                  aCTGTCTATGATTAATGAGGCGGTGGTTAACACTACCGCTTCGTTGTCTGGACAATCGTTCCTTTGTAACAGGTCCGAATCACAATGGAAGGTTCAAGAATGAAATACCGCATCGCTTTAGCTGTTTCTCTCTTTGctc                                                            <  2:456973/139‑1 (MQ=255)                                                                          tGATTAATGAGGCGGTGGTTAACACTACCGCTTCGTTGTCTGGACAATCGTTCCTTTGTAACAGGTCCGAATCACAATGGAAGGTTCAAGAATGAAATACCGCATCGCTTTAGCt                                                                            >  2:439125/1‑115 (MQ=255)                                                                          tGATTAATGAGGCGGTGGTTAACACTACCGCTTCGTTGTCTGGACAATCGTTCCTTTGTAACAGGTCCGAATCACAATGGAAGGTTCAAGAATGAAATACCGCATCGCTTTAGCt                                                                            <  1:439125/115‑1 (MQ=255)                                                                           gATTAATGAGGCGGTGGTTAACACTACCGCTTCGTTGTCTGGACAATCGTTCCTTTGTAACAGGTCCGAATCACAATGGAAGGTTCAAGAATGAAATACCGCATCGCTTTAGCTGTTTCTCTCTTTGCTCTTAGTGCCg                                                   <  1:86392/139‑1 (MQ=255)                                                                           gATTAATGAGGCGGGGGGTAACACTACCGCTTCGTTGTCTGGACAATCGTTCCTTTGTAACCGGTGCGAATCACAATGGAAGGTTCAAGAATGAAATACCGCATCGCTTTAGCTGTTTCTCTCTTTGCTCTTAGTGccc                                                   <  2:494545/139‑2 (MQ=255)                                                                             ttAATGAGGCGGTGGTTAACACTACCGCTTCGTTGTCTGGACAATCGTTCCTTTGTAACAGGTCCGAATCACAATGGAAGGTTCAAGAATGAAATACCGCATCGCTTTAGCTGTTTCTCTCTTTGCTCTTAGTGCCGGt                                                 >  2:620005/1‑139 (MQ=255)                                                                                aTGAGGCGGTGGTTAACACTACCGCTTCGTTGTCTGGACAATCGTTCCTTTGTAACAGGt                                                                                                                             >  1:343569/1‑60 (MQ=255)                                                                                aTGAGGCGGTGGTTAACACTACCGCTTCGTTGTCTGGACAATCGTTCCTTTGTAACAGGt                                                                                                                             <  2:343569/60‑1 (MQ=255)                                                                                aTGAGGCGGTGGTTAACACTACCGCTTCGTTGTCTGGACAATCGTTCCTTTGTAACAGGTCCGAATCACAATGGAAGGTTCAAGAATGAAATACCGCATCGCTTTAGCTGTTTCTCTCTTTGCTCTTAGTGCCGGTAGt                                              >  1:346099/1‑139 (MQ=255)                                                                                aTGAGGCGGTGGTTAACACTACCGCTTCGTTGTCTGGACAATCGTTCCTTTGTAACAGGTCCGAATCACAATGGAAGGTTCAAGAATGAAATACCGCATCGCTTTAGCTGTTTCTCTCTTTGCTCTTAGTGCCGGTAGt                                              <  2:346099/139‑1 (MQ=255)                                                                                aTGAGGCGGTGGTTAACACTACCGCTTCGTTGTCTGGACAATCGTTCCTTTGTAACAGGTCCGAATCACAATGGAAGGTTCAAGAATGAAATACCGCATCGCTTTAGCTGTTTCTCTCTTTGCTCTTAGTGCCGGTAGt                                              <  2:236984/139‑1 (MQ=255)                                                                                aTGAGGCGGTGGTTAACACTACCGCTTCGTTGTCTGGACAATCGTTCCTTTGTAACAGGTCCGAATCACAATGGAAGGTTCAAGAATGAAATACCGCATCGCTTTAGCTGTTTCTCTCTTTGCTCTTAGTGCCGGTAGt                                              >  2:357190/1‑139 (MQ=255)                                                                                aTGAGGCGGTGGTTAACACTACCGCTTCGTTGTCTGGACAATCGTTCCTTTGTAACAGGTCCGAATCACAATGGAAGGTTCAAGAATGAAATACCGCATCGCTTTAGCTGTTTCTCTCTTTGCTCTTAGTGCCGGTAGt                                              >  1:236984/1‑139 (MQ=255)                                                                                aTGAGGCGGTGGTTAACACTACCGCTTCGTTGTCTGGACAATCGTTCCTTTGTAACAGGTCCGAATCACAATGGAAGGTTCAAGAATGAAATACCGCATCGCTTTAGCTGTTTCTCTCTTTGCTCTTAGTGCCGGTAGt                                              >  1:772582/1‑139 (MQ=255)                                                                                  gAGGCGGTGGTTAACACTACCGCTTCGTTGTCTGGACAATCGTTCCTTTGTAACAGGTCCGAATCACAATGGAAGGTTCAAGAATGAAATACCGCATCGCTTTAGCTGTTTCTCTCTTTGCTCTTAGTGCCGGTAGTTa                                            <  1:795046/139‑1 (MQ=255)                                                                                      cGGTGGTTAACACTACCGCTTCGTTGTCTGGACAATCGTTCCTTTGTAACAGGTCCGAATCACAATGGAAGGTTCAAGAATGAAATACCGCATCGCttt                                                                                <  1:788413/99‑1 (MQ=255)                                                                                      cGGTGGTTAACACTACCGCTTCGTTGTCTGGACAATCGTTCCTTTGTAACAGGTCCGAATCACAATGGAAGGTTCAAGAATGAAATACCGCATCGCttt                                                                                >  2:788413/1‑99 (MQ=255)                                                                                       ggtggtTAACACTACCGCTTCGTTGTCTGGACAATCGTTCCTTTGTAACAGGTCCGAATCACAATGGAAGGTTCAAGAATGAAATACCGCATCGCTTTAGCTGTTTCTCTCTTTGCTCTTAGTGCCGGTAGTTATGCCa                                       <  1:357190/139‑1 (MQ=255)                                                                                       ggtggtTAACACTACCGCTTCGTTGTCTGGACAATCGTTCCTTTGTAACAGGTCCGAATCACAATGGAAGGTTCAAGAATGAAATACCGCATCGCTTTAGCTGTTTCTCTCTTTGCTCTTAGTGCCGGTAGTTATGCCa                                       <  1:496283/139‑1 (MQ=255)                                                                                         tggtTAACACTACCGCTTCGTTGTCTGGACAATCGTTCCTTTGTAACAGGTCCGAATCACAATGGAAGGTTCaa                                                                                                      <  2:362939/74‑1 (MQ=255)                                                                                         tggtTAACACTACCGCTTCGTTGTCTGGACAATCGTTCCTTTGTAACAGGTCCGAATCACAATGGAAGGTTCaa                                                                                                      >  1:362939/1‑74 (MQ=255)                                                                                            ttAACACTACCGCTTCGTTGTCTGGACAATCGTTCCTTTGTAACAGGTCCGAATCACAATGGAAGGTTCAAGAATGAAATACCGCATCGCTTTAGCTGTTTCTCTCTTTGCTCTTAGTGTCGGTAGTTATGCCACTAcc                                  <  2:105830/139‑1 (MQ=255)                                                                                            ttAACACTACCGCTTCGTTGTCTGGACAATCGTTCCTTTGTAACAGGTCCGAATCACAATGGAAGGTTCAAGAATGAAATACCGCATCGCTTTAGCTGTTTCTCTCTTTGCTCTTAGTGCCGGTAGTTATGCCACTAcc                                  <  1:774135/139‑1 (MQ=255)                                                                                            ttAACACTACCGCTTCGTTGTCTGGACAATCGTTCCTTTGTAACAGGTCCGAATCACAATGGAAGGTTCAAGAATGAAATACCGCATCGCTTTAGCTGTTTCTCTCTTTGCTCTTAGTGCCGGTAGTTATGCCACTAcc                                  <  1:20588/139‑1 (MQ=255)                                                                                                   tACCGCTTCGTTGTCTGGACAATCGTTCCTTTGTAACAGGTCCGAATCACAATGGAAGGTTCAAGAATGAAATACCGCATCGCTTTAGCTGTTTCTCTCTTTGCTCTTAGTGCCGGTAGTTATGCCACTACCCTGTGTc                           >  1:301014/1‑139 (MQ=255)                                                                                                   tACCGCTTCGTTGTCTGGACAATCGTTCCTTTGTAACAGGTCCGAATCACAATGGAAGGTACAAGAATGAAATACCGCATCGCTTTAGCTGTTTCTCTCTTTGCTATTAGTGCCGGAAGTTATGCCACTACCCTGTGTc                           >  1:237419/1‑139 (MQ=255)                                                                                                     ccGCTTCGTTGTCTGGACAATCGTTCCTTTGTAACAGGTCCGAATCACAATGGAAGGTTCAAGAATGAAATACCGCATCGCTTTAGCTGTTTCTCTCTTTGCTCTTAGt                                                       >  1:353809/1‑109 (MQ=255)                                                                                                     ccGCTTCGTTGTCTGGACAATCGTTCCTTTGTAACAGGTCCGAATCACAATGGAAGGTTCAAGAATGAAATACCGCATCGCTTTAGCTGTTTCTCTCTTTGCTCTTAGt                                                       <  2:353809/109‑1 (MQ=255)                                                                                                     ccGCTTCGTTGTCTGGACAATCGTTCCTTTGTAACAGGTCCGAATCACAATGGAAGGTTCAAGAATGAAATACCGCATCGCTTTAGCTGTTTCTCTCTTTGCTCTTAGTGTCGGTAGTTATg                                          >  2:622699/1‑122 (MQ=255)                                                                                                     ccGCTTCGTTGTCTGGACAATCGTTCCTTTGTAACAGGTCCGAATCACAATGGAAGGTTCAAGAATGAAATACCGCATCGCTTTAGCTGTTTCTCTCTTTGCTCTTAGTGTCGGTAGTTATg                                          <  1:622699/122‑1 (MQ=255)                                                                                                     ccGCTTCGTTGTCTGGACAATCGTTCCTTTGTAACAGGTCCGAATCACAATGGAAGGTTCAAGAATGAAATACCGCATCGCTTTAGCTGTTTCTCTCTTTGCTCTTAGTGCCGGt                                                 <  1:604539/115‑1 (MQ=255)                                                                                                     ccGCTTCGTTGTCTGGACAATCGTTCCTTTGTAACAGGTCCGAATCACAATGGAAGGTTCAAGAATGAAATACCGCATCGCTTTAGCTGTTTCTCTCTTTGCTCTTAGTGCCGGt                                                 >  2:604539/1‑115 (MQ=255)                                                                                                     ccGCTTCGTTGTCTGGACAATCGTTCCTTTGTAACAGGTCCGAATCACAATGGAAGGTTCAAGAATGAAATACCGCATCGCTTTAGCTGTTTCTCTCTTTGCTCTTAGTGCCGGTAGTTATGCCa                                       >  2:697987/1‑125 (MQ=255)                                                                                                     ccGCTTCGTTGTCTGGACAATCGTTCCTTTGTAACAGGTCCGAATCACAATGGAAGGTTCAAGAATGAAATACCGCATCGCTTTAGCTGTTTCTCTCTTTGCTCTTAGTGCCGGTAGTTATGCCa                                       <  1:697987/125‑1 (MQ=255)                                                                                                     ccGCTTCGTTGTCTGGACAATCGTTCCTTTGTAACAGGTCCGAATCACAATGGAAGGTTCAAGAATGAAATACCGCATCGCTTTAGCTGTTTCTCTCTTTGCTCTTAGTGCCGGTAGTTATGCCACTACCCTGTGTCAg                         >  1:534119/1‑139 (MQ=255)                                                                                                     ccGCTTCGTTGTCTGGACAATCGTTCCTTTGTAACAGGTCCGAATCACAATGGAAGGTTCAAGAATGAAATACCGCATCGCTTTAGCTGTTTCTCTCTTTGCTCTTAGTGCCGGTAGTTATGCCACTACCCTGTGTCAg                         >  1:259283/1‑139 (MQ=255)                                                                                                     ccGCTTCGTTGTCTGGACAATCGTTCCTTTGTAACAGGTCCGAATCACAATGGAAGGTTCAAGAATGAAATACCGCATCGCTTTAGCTGTTTCTCTCTTTGCTCTTAGTGCCGGTAGTTATGCCACTACCCTGTGTCAg                         >  2:71471/1‑139 (MQ=255)                                                                                                        cTTCGTTGTCTGGACAATCGTTCCTTTGTAACAGGTCCGAATCACAATGGAAGGTTCAAGAATGAAATACCGCATCGCTTTAGCTGTTTCTCTCTTTGCTCTTAGTGCCGGt                                                 <  1:362016/112‑1 (MQ=255)                                                                                                        cTTCGTTGTCTGGACAATCGTTCCTTTGTAACAGGTCCGAATCACAATGGAAGGTTCAAGAATGAAATACCGCATCGCTTTAGCTGTTTCTCTCTTTGCTCTTAGTGCCGGt                                                 >  2:362016/1‑112 (MQ=255)                                                                                                        cTTCGTTGTCTGGACAATCGTTCCTTTGTAACAGGTCCGAATCACAATGGAAGGTTCAAGAATGAAATACCGCATCGCTTTAGCTGTTTCTCTCTTTGCTCTTAGTGCCGGTAGTTATGCCa                                       <  2:340979/122‑1 (MQ=255)                                                                                                        cTTCGTTGTCTGGACAATCGTTCCTTTGTAACAGGTCCGAATCACAATGGAAGGTTCAAGAATGAAATACCGCATCGCTTTAGCTGTTTCTCTCTTTGCTCTTAGTGCCGGTAGTTATGCCa                                       >  2:501215/1‑122 (MQ=255)                                                                                                        cTTCGTTGTCTGGACAATCGTTCCTTTGTAACAGGTCCGAATCACAATGGAAGGTTCAAGAATGAAATACCGCATCGCTTTAGCTGTTTCTCTCTTTGCTCTTAGTGCCGGTAGTTATGCCa                                       >  1:340979/1‑122 (MQ=255)                                                                                                        cTTCGTTGTCTGGACAATCGTTCCTTTGTAACAGGTCCGAATCACAATGGAAGGTTCAAGAATGAAATACCGCATCGCTTTAGCTGTTTCTCTCTTTGCTCTTAGTGCCGGTAGTTATGCCACTACCCTGTGTCAGGaa                      >  2:466978/1‑139 (MQ=255)                                                                                                        cTTCATTGTCTGGACAATCGTTCCTTTGTAACAGGTCCGAATCACAATGGAAGGTTCAAGAATGAATTACCGCATCGCTTTAGCTGTTTCTCTCTTTGCTCTTAGTGCCGGTAGTTATGccc                                       <  1:501215/122‑2 (MQ=255)                                                                                                          tCGTTGTCTGGACAATCGTTCCTTTGTAACAGGTCCGAATCACAATGGAAGGTTCAAGAATGAAATACCGCATCGCTTTAGCTGTTTCTCTCTTTGCTCTTAGTGCCGGTAGTTATGCCACTACCCTGTGTCAGGaaaa                    >  2:182892/1‑139 (MQ=255)                                                                                                          tCGTTGTCTGGACAATCGCTCCTTTGTAACAGGTCCGAATCACAATGGAAGGTTCAAGAATGAAATACCGCATCGCTTTAGCTGTTTCTCTCTTTGCTCTTAGTGCCGGTAGTTATGCCACTACCCTGTGTCAGAaaaa                    >  2:529478/1‑139 (MQ=255)                                                                                                           cGTTGTCTGGACAATCGTTCCTTTGTAACAGGTCCGAATCACAATGGAAGGTTCAAGAATGAAATACCGCATCGCTTTAGCTGTTTCTCTCTTTGCTCTTAGAGCCGGTAGTTATGCCACTACCCTGTGTCAGGAAAAg                   >  2:504001/1‑139 (MQ=255)                                                                                                             ttGTCTGGACAATCGTTCCTTTGTAACAGGTCCGAATCACAATGGAAGGTTCAAGAATGAAATACCGCATCGCTTTAGCTGTTTCTCTCTTTGCTCTTAGTGCCGGTAGTTATGCCACTACCCTGTGTCAGGAAAAGGa                 <  1:620005/139‑1 (MQ=255)                                                                                                             ttGTCTGGACAATCGTTCCTTTGTAACAGGTCCGAATCACAATGGAAGGTTCAAGAATGAAATACCGCATCGCTTTAGCTGTTTCTCTCTTTGCTCTTAGTGCCGGTAGTTATGCCACTACCCTGTGTCAGGAAAAGGa                 <  1:71471/139‑1 (MQ=255)                                                                                                              tGTCTGGACAATCGTTCCTTTGTAACAGGTCCGAATCACAATGGAAGGTTCAAGAATGACATACCGCATCGCTTTAGCTGTTTCTCTCTTTGCTCTTAGTGCCGGTAGTTATGCCACTACCCTGTGTCAGGAAAAGGAg                <  1:182892/139‑1 (MQ=255)                                                                                                                   ggACAATCGTTCCTTTGTAACAGGTCCGAATCACAATGGAAGGTTCAAGAATGAAATACCGCATCGCTTTAGCTGTTTCTCTCTTTGCTCTTAGTGCCGGTAGTTATGCCACTACCCTGTGTCAGGAAAAGGAGCaaaa           <  1:730867/139‑1 (MQ=255)                                                                                                                    gACAATCGTTCCTTTGTAACAGGTCCGAATCACAATGGAAGGTTCAAGAATGAAATACCGCATCGCTTTAGCTGTTTCTCTCTTTGCTCTTAGTGCCGGTAGTTATGCCACTACCCTGTGTCAGGAAAAGGAGCaaaa           <  1:348571/138‑1 (MQ=255)                                                                                                                    gACAATCGTTCCTTTGTAACAGGTCCGAATCACAATGGAAGGTTCAAGAATGAAATACCGCATCGCTTTAGCTGTTTCTCTCTTTGCTCTTAGTGCCGGTAGTTATGCCACTACCCTGTGTCAGGAAAAGGAGCaaaa           >  2:348571/1‑138 (MQ=255)                                                                                                                         tCGTTCCTTTGTAACAGGTCCGAATCACAATGGAAGGTTCAAGAATGAAATACCGCATCGCTTTAGCTGTTTCTCTCTTTGCTCTTAGTGCCGGTAGTTATGCCACTAcc                                  <  1:296337/110‑1 (MQ=255)                                                                                                                         tCGTTCCTTTGTAACAGGTCCGAATCACAATGGAAGGTTCAAGAATGAAATACCGCATCGCTTTAGCTGTTTCTCTCTTTGCTCTTAGTGCCGGTAGTTATGCCACTAcc                                  >  2:296337/1‑110 (MQ=255)                                                                                                                            ttCCTTTGTAACAGGTCCGAATCACAATGGAAGGTTCAAGAATGAAATACCGCATCGCTTTAGCTGTTTCTCTCTTTGCTCTTAGTGCCGGTAGt                                              <  1:484913/95‑1 (MQ=255)                                                                                                                            ttCCTTTGTAACAGGTCCGAATCACAATGGAAGGTTCAAGAATGAAATACCGCATCGCTTTAGCTGTTTCTCTCTTTGCTCTTAGTGCCGGTAGTTATGCCACTACCCTGTGTCAGGAAAAGGAGCAAAATATCCTTaa  >  1:57311/1‑139 (MQ=255)                                                                                                                            ttCCTTTGTAACAGGTCCGAATCACAATGGAAGGTTCAAGAATGAAATACCGCATCGCTTTAGCTGTTTCTCTCTATGCTCTTAGTGCCGGTAGt                                              >  2:484913/1‑95 (MQ=255)                                                                                                                                       |                                                                                                                                TCTCATCGCTTCGGATAACTTCCTGTTATCCGAAACATATCACTCACTTTGGTGATTTCACCGTAACTGTCTATGATTAATGAGGCGGTGGTTAACACTACCGCTTCGTTGTCTGGACAATCGTTCCTTTGTAATAGGTCCGAATCACAATGGAAGGTTCAAGAATGAAATACCGCATCGCTTTAGCTGTTTCTCTCTTTGCTCTTAGTGCCGGTAGTTATGCCACTACCCTGTGTCAGGAAAAGGAGCAAAATATCCTTAA  >  NC_000913/3248805‑3249066

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